激光掃描和攝影測(cè)量在坡面侵蝕演變過(guò)程的適用性

羅 斌，張 勇，張志偉，倪世民，張 歆，王軍光

激光掃描和攝影測(cè)量在坡面侵蝕演變過(guò)程的適用性

羅 斌1，張 勇2，張志偉1，倪世民1，張 歆2，王軍光1※

（1. 華中農(nóng)業(yè)大學(xué)水土保持研究中心，武漢 430070；2. 長(zhǎng)江水利委員會(huì)長(zhǎng)江流域水土保持監(jiān)測(cè)中心站，武漢 430010）

為研究激光掃描和攝影測(cè)量技術(shù)在監(jiān)測(cè)坡面侵蝕演變過(guò)程中的精度及適用性，該研究利用近景攝影測(cè)量技術(shù)和三維激光掃描技術(shù)對(duì)長(zhǎng)歷時(shí)條件下坡面侵蝕演變過(guò)程進(jìn)行監(jiān)測(cè)，獲取不同時(shí)段的坡面微地形數(shù)據(jù)，基于坡面高精度數(shù)字高程模型（Digital Elevation Model，DEM）對(duì)坡面侵蝕演變過(guò)程進(jìn)行分析，探究2種非接觸式測(cè)量方法在坡面侵蝕監(jiān)測(cè)中的適用性和精確度。結(jié)果表明：1）按主導(dǎo)侵蝕方式的不同，坡面侵蝕過(guò)程可分為片蝕階段、細(xì)溝發(fā)育階段和細(xì)溝成熟階段；2）2種非接觸式測(cè)量方法均能夠精確地對(duì)坡面侵蝕產(chǎn)沙過(guò)程進(jìn)行監(jiān)測(cè)，最大相對(duì)誤差為?16.82 %，2種方法在坡面侵蝕量測(cè)量方面有很好的適用性。3）近景攝影測(cè)量技術(shù)在坡面侵蝕產(chǎn)沙監(jiān)測(cè)、細(xì)溝深度測(cè)量和坡面微地形模擬方面要優(yōu)于三維激光掃描技術(shù)。該研究可土壤坡面侵蝕監(jiān)測(cè)方法的選擇提供參考。

坡面；侵蝕；土壤；近景攝影測(cè)量技術(shù)；三維激光掃描技術(shù)；對(duì)比

0 引 言

土壤侵蝕會(huì)造成耕地破壞、農(nóng)田肥力降低，進(jìn)一步加劇現(xiàn)已十分緊張的人口、糧食和環(huán)境矛盾，嚴(yán)重影響人類的生存和發(fā)展[1]。坡面作為土壤侵蝕發(fā)生的基本單元，研究其侵蝕演變過(guò)程對(duì)深入土壤侵蝕機(jī)理研究具有重要意義[2-3]。

長(zhǎng)期以來(lái)，眾多學(xué)者就坡面侵蝕開(kāi)展了大量調(diào)查和研究。在早期，相關(guān)土壤侵蝕研究多基于野外天然降雨條件下展開(kāi)。然而，天然降雨具有很強(qiáng)的不確定性，對(duì)于進(jìn)一步深入土壤侵蝕機(jī)理研究具有較大限制。因此，目前多采用室內(nèi)人工模擬降雨進(jìn)行土壤侵蝕機(jī)理研究，在此條件下，有關(guān)學(xué)者對(duì)土壤坡面侵蝕的發(fā)生以及演化過(guò)程開(kāi)展了一系列研究，剖析了坡面侵蝕的關(guān)鍵過(guò)程和影響因子[4-5]。盡管土壤侵蝕過(guò)程和侵蝕機(jī)理等相關(guān)研究在不斷深入，然而研究中常用的傳統(tǒng)填土法[6]、直尺測(cè)量法[7]和測(cè)針板法[8]等方法，因效率低、對(duì)坡面擾動(dòng)大、易受人工讀數(shù)誤差影響等原因展現(xiàn)了較大的局限性。

近年來(lái)，隨著高新技術(shù)的發(fā)展，諸如近景攝影測(cè)量技術(shù)、三維激光掃描技術(shù)[9-10]等非接觸式測(cè)量方法因測(cè)量效率高、誤差小等特點(diǎn)被廣泛應(yīng)用于土壤侵蝕監(jiān)測(cè)和測(cè)量中。Gao等[11]基于三維激光掃描，提出基于點(diǎn)云數(shù)據(jù)和地形變化信息的侵蝕量計(jì)算方法；Jiang等[12]在室內(nèi)模擬降雨條件下，借助數(shù)字近景攝影測(cè)量技術(shù)監(jiān)測(cè)連續(xù)降雨期間坡面細(xì)溝侵蝕的動(dòng)態(tài)過(guò)程，對(duì)細(xì)溝侵蝕過(guò)程進(jìn)行研究；明旭輝等[13]利用近景攝影測(cè)量計(jì)算凍融坡面侵蝕產(chǎn)沙過(guò)程的產(chǎn)沙量，與實(shí)測(cè)值間的誤差平均精度為90.67%。李俊利等[14]利用攝影測(cè)量技術(shù)和三維激光掃描對(duì)坡面侵蝕溝進(jìn)行三維重建，通過(guò)對(duì)比2種技術(shù)估算的侵蝕量發(fā)現(xiàn)，二者相對(duì)誤差僅為8.05%。

綜上，已有研究證實(shí)了三維激光掃描和攝影測(cè)量技術(shù)在坡面侵蝕試驗(yàn)中的可行性，但對(duì)于二者在坡面侵蝕過(guò)程中精度對(duì)比及其適用性尚不明確，尤其是在細(xì)溝形態(tài)特征方面的對(duì)比研究仍有不足。

因此，本文基于室內(nèi)人工模擬降雨試驗(yàn)利用2種非接觸式測(cè)量方法對(duì)土壤坡面侵蝕過(guò)程進(jìn)行監(jiān)測(cè)，研究坡面侵蝕的動(dòng)態(tài)變化過(guò)程，分析坡面細(xì)溝形態(tài)特征的演變，同時(shí)對(duì)2種非接觸式測(cè)量方法的測(cè)量精度進(jìn)行對(duì)比，并對(duì)2種方法侵蝕監(jiān)測(cè)試驗(yàn)中的適用性進(jìn)行評(píng)價(jià)，以期為土壤坡面侵蝕監(jiān)測(cè)方法的選擇提供參考。

1 材料與方法

1.1 試驗(yàn)土壤

試驗(yàn)土壤為花崗巖母質(zhì)發(fā)育紅壤，采自福建長(zhǎng)汀縣的典型侵蝕區(qū)，根據(jù)林分調(diào)查（植被群落）、裸露程度和田間土壤剖面的診斷，在崩崗侵蝕區(qū)附近分別選擇了位于坡中位置的4個(gè)采樣點(diǎn)，采樣點(diǎn)的剖面構(gòu)型為A-B-BC-C型，由淋溶層（A，0～0.3 m）、淀積層（B，0.3～6 m）、過(guò)渡層（BC）和母質(zhì)層（C）組成。采樣點(diǎn)按發(fā)生層次自下而上分別采集母質(zhì)層、淀積層和淋溶層的土壤作為試驗(yàn)材料。采樣過(guò)程中盡可能選擇剖面結(jié)構(gòu)完整、地形、土地利用方式一致的樣地進(jìn)行采樣。淋溶層土壤采集時(shí)在采樣點(diǎn)15 m半徑范圍內(nèi)按照“S”形選擇了10個(gè)樣點(diǎn)，采集土壤后將其混合作為試驗(yàn)材料；下層土壤則在挖取的典型剖面上采集。上述土樣被帶回實(shí)驗(yàn)室后，將其進(jìn)行自然風(fēng)干，并通過(guò)10 mm篩網(wǎng)對(duì)土樣進(jìn)行篩分，去除根系等雜質(zhì)后備用，試驗(yàn)還采集了散土樣和3個(gè)重復(fù)的標(biāo)準(zhǔn)環(huán)刀樣（100 cm3）用于測(cè)量土樣基本理化性質(zhì)。試驗(yàn)土壤的理化性質(zhì)采用常規(guī)方法[15-16]測(cè)定，其基本的理化性質(zhì)如表1所示。

表1 試驗(yàn)土樣的基本理化性質(zhì)

1.2 室內(nèi)試驗(yàn)設(shè)計(jì)

模擬降雨試驗(yàn)采用華中農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院水土保持研究中心降雨大廳下噴式模擬降雨器（西安清遠(yuǎn)，QYJY-503T）[17-18]。該降雨器的噴頭組有小、中、大3種規(guī)格的噴頭組成，3個(gè)噴頭在空間上重合疊加，形成一個(gè)雨強(qiáng)分布均勻的降雨區(qū)，降雨均勻度大于90%。降雨高度為10 m，降雨雨滴通過(guò)噴頭下噴后在地面可以達(dá)到雨滴終速。降雨器通過(guò)控制系統(tǒng)調(diào)節(jié)噴頭組合和上水壓力可以實(shí)現(xiàn)15～220 mm/h的降雨強(qiáng)度范圍。

本研究在模擬降雨過(guò)程中根據(jù)試驗(yàn)?zāi)康氖褂玫耐敛垡?guī)格為4.0 m（長(zhǎng)）×2.0 m（寬）×0.60 m（高），可填土深度為0.50 m，坡度可在0°～45°之間調(diào)節(jié)。槽下端安裝V型鋼槽用于收集徑流和泥沙，土槽底板均勻分布有排水孔。土槽邊框貼有鋼尺，便于快捷獲取坡面位置。土槽分為左右平行兩區(qū)，土槽兩側(cè)設(shè)置9個(gè)非對(duì)稱標(biāo)靶，用直尺量測(cè)標(biāo)靶位置，作為攝影測(cè)量建?？刂泣c(diǎn)（圖1）。試驗(yàn)土槽內(nèi)從下至上依次按環(huán)刀測(cè)得的容重填裝厚度為15 cm的母質(zhì)層、20 cm的淀積層及15 cm的淋溶層土壤。在試驗(yàn)開(kāi)始前，對(duì)土槽內(nèi)土壤進(jìn)行多次干濕循環(huán)和長(zhǎng)達(dá)2個(gè)月的靜置沉降，以使土壤層次和顆粒盡可能接近自然條件。土槽采用分層填裝的方法，邊填土邊壓實(shí)，以每層5 cm的增量進(jìn)行填裝，將土壤壓實(shí)后，使用水平儀將土壤表面調(diào)整平整。在正式試驗(yàn)開(kāi)始前，在2種土槽的土壤表面覆蓋一層紗網(wǎng)（孔徑2 mm），使用20 mm/h雨強(qiáng)對(duì)其進(jìn)行預(yù)濕潤(rùn)處理，直至坡面開(kāi)始產(chǎn)流為止。紗網(wǎng)距離坡面表面高度為10 cm，以降低雨滴打擊對(duì)表土結(jié)構(gòu)的影響。預(yù)濕潤(rùn)降雨的目的有3個(gè)：1）保持下墊面土壤初始含水量一致；2）通過(guò)降雨濕潤(rùn)作用固結(jié)表面分散孤立的土壤顆粒；3）減少下墊面表面條件的空間變異性。在預(yù)濕潤(rùn)后，用塑料布將試驗(yàn)土槽遮蓋，將其靜置24 h自由排掉重力水。當(dāng)使用時(shí)域反射儀（Time Domain Reflectometry，TDR）測(cè)定的坡面土壤含水率降至30%±0.5%時(shí)，開(kāi)始進(jìn)行試驗(yàn)。

圖1 土槽和標(biāo)靶布設(shè)位置

在試驗(yàn)開(kāi)始前，降雨強(qiáng)度設(shè)定為（90±5）mm/h，在降雨器或供水流量穩(wěn)定后率定降雨強(qiáng)度和上方匯水流量，當(dāng)降雨強(qiáng)度或流量實(shí)測(cè)值和目標(biāo)值的差值小于5%時(shí)正式試驗(yàn)，雨強(qiáng)通過(guò)放置在土槽周圍的6個(gè)雨量筒進(jìn)行監(jiān)測(cè)，匯水流量通過(guò)稱質(zhì)量法進(jìn)行測(cè)定?；趯?duì)采樣點(diǎn)集水區(qū)坡度的原位調(diào)查，降雨試驗(yàn)均在15°坡度上進(jìn)行，這個(gè)坡度也是中國(guó)低山丘陵區(qū)坡耕地的常見(jiàn)坡度[19-20]。試驗(yàn)開(kāi)始后，觀察并記錄坡面的初始產(chǎn)流時(shí)間、細(xì)溝出現(xiàn)的時(shí)間和位置等，量測(cè)細(xì)溝的溝頭位置和深度。坡面產(chǎn)流后，0～60 min內(nèi)每隔1 min，>60～180 min內(nèi)每隔3 min，>180～600 min內(nèi)每隔6 min，>600～3 000 min內(nèi)每隔12 min，>3 000～6 000 min內(nèi)每隔20 min，接取一次徑流泥沙樣，接樣時(shí)間為1 min，每組接樣裝置包括2個(gè)250 mL鋁盒和一個(gè)12 L接樣桶，其中鋁盒用于計(jì)算徑流含沙量。同時(shí)，試驗(yàn)過(guò)程通過(guò)安裝在土槽上方降雨器上的高清相機(jī)（Sony6300）每隔10 s進(jìn)行全程拍照，記錄細(xì)溝發(fā)育過(guò)程和表面形態(tài)的動(dòng)態(tài)數(shù)據(jù)。單場(chǎng)模擬降雨或徑流試驗(yàn)的持續(xù)時(shí)間為坡面產(chǎn)流后1 h，每2場(chǎng)間歇性降雨時(shí)間之間間隔24 h，降雨總計(jì)100 h。降雨后，觀測(cè)坡面表面形態(tài)，使用三維激光掃描技術(shù)和近景攝影測(cè)量技術(shù)分別記錄坡面表面微地形形態(tài)，獲取表面數(shù)字點(diǎn)云信息。同時(shí)，使用直尺法每隔10 cm測(cè)量坡面各條細(xì)溝的長(zhǎng)度、寬度和深度。

1.3 數(shù)據(jù)采集

1.3.1 數(shù)字近景攝影測(cè)量

每場(chǎng)降雨后，將土槽放平，利用數(shù)碼相機(jī)（Sony6300）從不同角度和高度對(duì)坡面進(jìn)行拍攝，采集坡面微地形信息。照片位置重復(fù)率控制在70%以上。坡面影像采集完畢后，導(dǎo)入到Agisoft Metashape Professional軟件中。通過(guò)軟件中的相機(jī)標(biāo)定工具利用試驗(yàn)相機(jī)的標(biāo)定參數(shù)完成拍攝照片的畸變標(biāo)定，借助對(duì)齊照片工具和建立密集點(diǎn)云工具完成點(diǎn)云匹配并建立試驗(yàn)所需要的高密度點(diǎn)云。

1.3.2 三維激光掃描技術(shù)

在每場(chǎng)降雨后，將土槽放平，以減少掃描死角。由于試驗(yàn)場(chǎng)地限制，三維激光掃描儀的架設(shè)位置在距土槽左右兩側(cè)最遠(yuǎn)2 m以內(nèi)。因?yàn)樵囼?yàn)土槽面積較大、高度較高等原因，對(duì)于坡面的掃描測(cè)量無(wú)法在一個(gè)測(cè)站的條件下完成。為了完整地將坡面的微地形數(shù)據(jù)記錄下來(lái)，采用分測(cè)站掃描的方法，設(shè)置5個(gè)不同位置的測(cè)站（圖2），從不同的角度對(duì)坡面進(jìn)行掃描，以獲得更加全面的坡面三維點(diǎn)云數(shù)據(jù)。為了將5個(gè)測(cè)站的數(shù)據(jù)較好地拼接到同一個(gè)坐標(biāo)系中，在每次試驗(yàn)前在土槽上放置3個(gè)標(biāo)靶，標(biāo)靶位置按照儀器標(biāo)靶放置規(guī)則固定于土槽的不同位置。根據(jù)試驗(yàn)?zāi)康暮途纫?，點(diǎn)云分辨率設(shè)置為1 mm×1 mm。在Faro Scene V4.7.3軟件中，對(duì)5站數(shù)據(jù)進(jìn)行配置，生成高精度點(diǎn)云數(shù)據(jù)。

圖2 三維激光掃描儀測(cè)站位置分布

1.3.3 點(diǎn)云數(shù)據(jù)處理

將得到的攝影測(cè)量和三維激光掃描點(diǎn)云數(shù)據(jù)輸入到Cloud Compare V2.11.3軟件中。利用軟件的降噪功能對(duì)點(diǎn)云數(shù)據(jù)進(jìn)行降噪；借助裁剪工具對(duì)坡面數(shù)據(jù)進(jìn)行裁剪，提取研究區(qū)。然后將點(diǎn)云數(shù)據(jù)輸入到ArcGIS10.2中，利用ArcGIS 10.2中的數(shù)據(jù)轉(zhuǎn)換工具將點(diǎn)云數(shù)據(jù)轉(zhuǎn)換為柵格數(shù)據(jù)。最后用地理工具對(duì)柵格數(shù)據(jù)進(jìn)行空間校準(zhǔn)，得到降雨前原始坡面和不同降雨歷時(shí)下坡面的高精度DEM數(shù)據(jù)（精度為5 mm×5 mm）。

1.4 數(shù)據(jù)處理與分析

1.4.1 土壤侵蝕計(jì)算

1）徑流和泥沙收集方法

收集徑流泥沙是土壤流失試驗(yàn)中最可靠的方法[21]。將試驗(yàn)過(guò)程中收集的鋁盒樣品進(jìn)行稱量，然后，將其靜置24 h，倒掉上層清水，放置在105 ℃烘箱烘至恒質(zhì)量，以計(jì)算產(chǎn)沙量收集了降雨試驗(yàn)的徑流泥沙樣。隨后，將計(jì)算得到的產(chǎn)沙量作為實(shí)測(cè)值計(jì)算土壤流失量，并將徑流和泥沙收集的結(jié)果與數(shù)字?jǐn)z影測(cè)量攝影測(cè)量和三位激光掃描觀測(cè)結(jié)果進(jìn)行比較。

2）基于數(shù)字圖像的計(jì)算

土壤侵蝕體積的計(jì)算原理是基于DEM的體積差。在計(jì)算中，降雨前坡面的原始DEM減去侵蝕發(fā)生后不同時(shí)間段坡面的DEM，將每個(gè)柵格作為體積計(jì)算的微分單元，用該單元的面積乘以該格網(wǎng)點(diǎn)到參考面的距離來(lái)計(jì)算土壤侵蝕量，將計(jì)算的土壤侵蝕體積與容重進(jìn)行換算得到土壤流失量[22]。其數(shù)學(xué)表達(dá)式為

式中為土壤侵蝕體積，kg；和為DEM數(shù)據(jù)的行列數(shù)；(,)為DEM格網(wǎng)上編號(hào)是(,)的點(diǎn)到參考面的距離；D為DEM單元格的面積。

1.4.2 細(xì)溝形態(tài)特征參數(shù)的計(jì)算

1）細(xì)溝寬深比

細(xì)溝寬深比是坡面某一坡位處細(xì)溝的寬度（m）和深度（m）的比值，能夠反映細(xì)溝的斷面形態(tài)在水平和豎直方向的大小變化，可以作為坡面侵蝕發(fā)育強(qiáng)度的重要參考指標(biāo)[23]。

2）細(xì)溝總表面積

細(xì)溝總表面積（m2）也可以稱為細(xì)溝實(shí)際破壞面積，是指研究區(qū)域內(nèi)所有細(xì)溝的平面面積之和，在本試驗(yàn)中基于細(xì)溝矢量進(jìn)行細(xì)溝面積計(jì)算[8]。

2 結(jié)果與分析

2.1 基于攝影測(cè)量和三維激光掃描的坡面侵蝕演變形態(tài)觀測(cè)

在降雨試驗(yàn)中，坡面經(jīng)歷了片蝕階段-細(xì)溝發(fā)育階段-細(xì)溝成熟階段3個(gè)侵蝕演化過(guò)程（圖3），攝影測(cè)量技術(shù)和三維激光掃描技術(shù)觀測(cè)的坡面侵蝕過(guò)程均出現(xiàn)相同的演變趨勢(shì)。在1～5 h期間，坡表面主要以片蝕為主，這是由于坡面徑流的能量較小，對(duì)坡面的沖刷能力較弱，只能帶走土壤表面細(xì)小顆粒和一些較為松散的物質(zhì)，而坡面中下部坡面匯水面積較大，坡面的徑流量和流速相較坡頂明顯增加，當(dāng)坡面徑流侵蝕力大于坡面抗侵蝕的能力后，坡面土壤開(kāi)始被剝離，在此階段坡中和坡下受侵蝕影響較大，坡面中部在水流沖刷作用下最先形成跌坎，然后由于坡面徑流在跌坎內(nèi)不斷聚集，跌坎進(jìn)一步發(fā)育形成斷續(xù)細(xì)溝。吳淑芳等[22]研究了坡面細(xì)溝的侵蝕演變過(guò)程中，也同樣發(fā)現(xiàn)在前期降雨試驗(yàn)后，侵蝕主要集中在坡面中下部，而頂部由于匯水面積較小，水流剪切力不足以剝離表層土壤，而隨著跌坎的出現(xiàn)也為細(xì)溝出現(xiàn)創(chuàng)造條件。從圖3中可以看出，在降雨歷時(shí)1 h后，坡面100和160 cm附近開(kāi)始出現(xiàn)跌坎。歷時(shí)5 h后，坡面100 cm處演變?yōu)槎鄺l斷續(xù)細(xì)溝。

在降雨歷時(shí)5～10 h期間，斷續(xù)細(xì)溝通過(guò)溝頭的不斷溯源連通為連續(xù)細(xì)溝，形成固定流路，細(xì)溝侵蝕開(kāi)始出現(xiàn)，水力作用從雨滴濺蝕和片狀侵蝕轉(zhuǎn)化為集中徑流侵蝕[24]，坡面通過(guò)不斷匯水產(chǎn)生集中徑流，同時(shí)細(xì)溝及細(xì)溝間的區(qū)域使得剝蝕的土壤得以快速輸移[25]。

根據(jù)細(xì)溝發(fā)育主導(dǎo)侵蝕方式的不同，可將細(xì)溝發(fā)育過(guò)程分為：細(xì)溝發(fā)育前期和細(xì)溝發(fā)育后期2個(gè)部分[26]。一般情況下，細(xì)溝集中徑流的剝蝕能力和輸移能力遠(yuǎn)大于薄層徑流[27]，本試驗(yàn)結(jié)果也與上述相一致，在細(xì)溝發(fā)育前期（5～50 h）徑流相對(duì)集中，溝內(nèi)水流流速快、流量大，對(duì)坡面的侵蝕能力和泥沙搬運(yùn)能力明顯增強(qiáng)，以溝頭溯源侵蝕和溝底下切侵蝕為主。5～10 h期間，細(xì)溝長(zhǎng)度由58.37 cm增長(zhǎng)到112.11 cm，細(xì)溝最大深度由4.35 cm增長(zhǎng)到5.45 cm，該時(shí)期是細(xì)溝發(fā)育最活躍的階段，隨后細(xì)溝前進(jìn)速度逐漸變緩，細(xì)溝下切侵蝕作用逐漸增強(qiáng)。10～50 h期間，細(xì)溝長(zhǎng)度由112.11 cm增長(zhǎng)到216.43 cm，細(xì)溝下切最大深度由5.37 cm增加為27.32 cm，細(xì)溝大部分地方已經(jīng)下切至淀積層。

圖3 連續(xù)降雨條件下的坡面侵蝕演變過(guò)程

在細(xì)溝發(fā)育后期（50～70 h），由于細(xì)溝長(zhǎng)度的不斷增加，大量坡面徑流開(kāi)始從細(xì)溝兩側(cè)流入細(xì)溝中，細(xì)溝內(nèi)徑流量明顯增加，溝壁擴(kuò)張侵蝕開(kāi)始出現(xiàn)。由于溝底已經(jīng)下切至抗蝕能力較強(qiáng)的淀積層，所以溝底下切作用逐漸減弱。溝中徑流能量被分散于溯源侵蝕和溝壁擴(kuò)張侵蝕，溝頭在坡面不斷向前推進(jìn)，溝壁不斷向兩側(cè)擴(kuò)展，細(xì)溝長(zhǎng)度由216.43 cm增長(zhǎng)到256.18 cm，最大溝深由27.32 cm增長(zhǎng)到27.86 cm，細(xì)溝最大溝寬由19.65 cm增長(zhǎng)到25.75 cm。

在細(xì)溝成熟階段，溝壁擴(kuò)張侵蝕占據(jù)主導(dǎo)地位（圖3），細(xì)溝下切侵蝕幾乎停止，溝頭溯源侵蝕發(fā)育緩慢，坡面侵蝕作用主要來(lái)自于溝壁擴(kuò)張侵蝕，溝壁土塊的崩塌成為坡面侵蝕產(chǎn)沙的主要來(lái)源，細(xì)溝發(fā)育較為穩(wěn)定。溝腳被細(xì)溝中徑流掏涮之后，當(dāng)驅(qū)動(dòng)力（土體自身重力）大于阻力（土壤顆粒間的粘結(jié)力和摩擦力等）時(shí)，溝壁內(nèi)側(cè)容易形成張力裂隙[28-29]。隨著溝腳被集中水流進(jìn)一步淘涮，懸空土體體積增大，溝壁開(kāi)始隨機(jī)性崩塌，溝寬增大[30-32]。在降雨歷時(shí)70 h后，溝頭溯源和下切基本停止。在隨后的試驗(yàn)中，坡面侵蝕主要是溝壁擴(kuò)張侵蝕引起的溝壁坍塌。在此階段，細(xì)溝長(zhǎng)度由256.18 cm演變?yōu)?87.35 cm，細(xì)溝最大寬度由25.75 cm發(fā)展為35.4 cm，細(xì)溝最大溝深由27.86 cm增長(zhǎng)到29.40 cm。

2.2 三維激光掃描和攝影測(cè)量的精度分析

2.2.1 坡面侵蝕量測(cè)量精度對(duì)比

在坡面土壤侵蝕量演變觀測(cè)試驗(yàn)中，以徑流泥沙采集法所得到的土壤侵蝕量結(jié)果作為標(biāo)準(zhǔn)，對(duì)近景攝影測(cè)量法和三維激光掃描法的測(cè)量精度進(jìn)行評(píng)價(jià)。利用每次降雨間隙采集到的DEM數(shù)據(jù)進(jìn)行解算，可得到相應(yīng)時(shí)段內(nèi)不同方法所測(cè)得的土壤侵蝕量。

表2為3種測(cè)量方法在不同降雨時(shí)間段內(nèi)測(cè)量精度對(duì)比結(jié)果。由表可知，2種非接觸式測(cè)量方法的測(cè)量值的相對(duì)誤差多為負(fù)數(shù)，測(cè)量值多小于實(shí)際采集的泥沙質(zhì)量。這說(shuō)明在花崗巖紅壤土壤侵蝕試驗(yàn)中，2種非接觸式測(cè)量方法會(huì)低估真實(shí)的土壤侵蝕量。但總的來(lái)說(shuō)2種非接觸式測(cè)量方法的測(cè)量結(jié)果和實(shí)測(cè)值比相對(duì)誤差較小，最大相對(duì)誤差為?16.82%，2種方法在坡面侵蝕量測(cè)量方面有很好的適用性[33-34]。

從試驗(yàn)過(guò)程來(lái)看，2種方法的測(cè)量誤差呈相同的變化趨勢(shì)。在試驗(yàn)前期和末期2種測(cè)量方法的測(cè)量誤差均較小，而在試驗(yàn)中期測(cè)量誤差則較大。這是因?yàn)樵谠囼?yàn)中期，坡面細(xì)溝窄而深，并且坡面中部細(xì)溝橫截面多為細(xì)頸瓶型即細(xì)溝底部較寬肩部較窄，所以會(huì)存在較多的數(shù)據(jù)采集盲區(qū)。在點(diǎn)云和DEM處理過(guò)程中，軟件在填補(bǔ)盲區(qū)時(shí)會(huì)將盲區(qū)進(jìn)行平滑化處理，所以會(huì)低估細(xì)溝的深度和細(xì)溝底部的寬度，從而影響到細(xì)溝體積的測(cè)量，造成測(cè)量誤差。在試驗(yàn)前期近景攝影測(cè)量的測(cè)量誤差要大于三維激光掃描，這是因?yàn)樵囼?yàn)前期數(shù)據(jù)采集掃描盲區(qū)較小，三維激光掃描儀采集到的點(diǎn)云數(shù)據(jù)數(shù)量更多，構(gòu)建的DEM分辨率也更高，所以測(cè)量誤差較小。然而隨著試驗(yàn)的進(jìn)行，表現(xiàn)為近景攝影測(cè)量的測(cè)量誤差小于三維激光掃描。這是因?yàn)樵谠囼?yàn)中后期細(xì)溝深且窄，三維激光掃描儀在使用中需要固定架設(shè)，所以會(huì)存在較多的數(shù)據(jù)采集盲區(qū)。而近景攝影測(cè)量所使用的數(shù)碼相機(jī)在使用時(shí)無(wú)需固定架設(shè)，可根據(jù)細(xì)溝形態(tài)選擇不同的拍攝角度，所以數(shù)據(jù)盲區(qū)相對(duì)會(huì)較小，因此測(cè)量誤差較小。其中近景攝影測(cè)量技術(shù)的平均絕對(duì)誤差為8.84%，三維激光掃描技術(shù)的平均絕對(duì)誤差為8.97%，近景攝影測(cè)量技術(shù)的測(cè)量精度要更高。并且徑流泥沙采集法、近景攝影測(cè)量法和三維激光掃描法所測(cè)得的累積土壤侵蝕量呈相同的變化趨勢(shì)，但由于誤差的存在，2種方法計(jì)算的泥沙量普遍小于徑流泥沙法所測(cè)得的泥沙量，所以在試驗(yàn)中后期2種非接觸式方法測(cè)得的累計(jì)徑流泥沙量越來(lái)越小于實(shí)際測(cè)得的泥沙量。

表2 2種非接觸式測(cè)量法的土壤侵蝕量測(cè)量誤差

2.2.2 坡面形態(tài)特征測(cè)量精度對(duì)比

1）細(xì)溝長(zhǎng)寬深測(cè)量結(jié)果對(duì)比

對(duì)于細(xì)溝長(zhǎng)度的量測(cè)，在每次試驗(yàn)間隔借助土槽兩側(cè)的直尺讀取溝頭和溝腳的位置來(lái)測(cè)定，將該實(shí)測(cè)值作為評(píng)價(jià)2種非接觸式測(cè)量方法的標(biāo)準(zhǔn)。在Cloud Compare V2.11.3軟件中讀取近景攝影測(cè)量點(diǎn)云數(shù)據(jù)和三維激光掃描點(diǎn)云數(shù)據(jù)中點(diǎn)云三維坐標(biāo)（,,），并利用測(cè)量點(diǎn)云點(diǎn)的坐標(biāo)計(jì)算細(xì)溝形態(tài)參數(shù)。

在每次試驗(yàn)間隔利用直尺每隔10 cm對(duì)細(xì)溝寬度和細(xì)溝深度進(jìn)行量測(cè)，將量測(cè)的細(xì)溝寬度和細(xì)溝深度值作為評(píng)價(jià)2種非接觸式測(cè)量方法的標(biāo)準(zhǔn)。然后在Cloud Compare V2.11.3軟件中獲取最大溝寬處點(diǎn)云點(diǎn)的三維坐標(biāo)對(duì)細(xì)溝最大寬度進(jìn)行計(jì)算，由于2種非接觸式測(cè)量方法在試驗(yàn)中期存在數(shù)據(jù)采集盲區(qū)，所以選用2種方法生成的DEM計(jì)算細(xì)溝深度?；跍y(cè)得的最大細(xì)溝寬度進(jìn)行分析，從圖4a中可以看出，3種量測(cè)方法所測(cè)得的最大溝寬并無(wú)明顯差異，但用直尺量測(cè)的最大溝寬總體上要略大于2種非接觸式量測(cè)方法的測(cè)量值。這是因?yàn)樵诩?xì)溝量測(cè)過(guò)程中，試驗(yàn)土槽高度較高且細(xì)溝兩側(cè)坡面也受到不同程度的侵蝕，所以在量測(cè)過(guò)程中會(huì)出現(xiàn)直尺兩端高低差異，進(jìn)而出現(xiàn)測(cè)量數(shù)值偏大的結(jié)果?？偟膩?lái)看，2種非接觸式測(cè)量方法在細(xì)溝寬度量測(cè)方面具有很好的適用性，并在測(cè)量效率方面要優(yōu)于直尺量測(cè)。

對(duì)于細(xì)溝長(zhǎng)度的量測(cè)，在每次試驗(yàn)間隔借助土槽兩側(cè)的直尺讀取溝頭和溝腳的位置來(lái)測(cè)定，并將該實(shí)測(cè)值作為評(píng)價(jià)兩種非接觸式測(cè)量方法的標(biāo)準(zhǔn)。讀取近景攝影測(cè)量點(diǎn)云數(shù)據(jù)和三維激光掃描點(diǎn)云數(shù)據(jù)中細(xì)溝兩端的點(diǎn)云三維坐標(biāo)（,,），將測(cè)量點(diǎn)云點(diǎn)的坐標(biāo)進(jìn)行相減就可得到對(duì)應(yīng)時(shí)間段細(xì)溝的長(zhǎng)度。從圖4b中可以看出，3種測(cè)量方法所得到的細(xì)溝長(zhǎng)度演變特征一致，細(xì)溝長(zhǎng)度差異較小，其中近景攝影測(cè)量和三維激光掃描的最大相對(duì)誤差分別為2.7%和2.3%。2種非接觸式測(cè)量方法在細(xì)溝長(zhǎng)度測(cè)量方面均具有較好的準(zhǔn)確度。

基于圖4c中最大溝深進(jìn)行分析發(fā)現(xiàn)，在試驗(yàn)前期和后期2種方法的深度測(cè)量值的最大相對(duì)誤差為18.5%，而在試驗(yàn)中期測(cè)量值的最大相對(duì)誤差為25.3%，證明在細(xì)溝深度測(cè)量方面，前后期測(cè)量結(jié)果要優(yōu)于中期。同時(shí)對(duì)比細(xì)溝寬深比后發(fā)現(xiàn)在試驗(yàn)中期最大寬深平均值為0.84，而前期和后期最大寬深比平均值分別為1.42和1.24。在試驗(yàn)中期細(xì)溝寬度較窄并且深度較深，在數(shù)據(jù)采集時(shí)會(huì)存在數(shù)據(jù)采集盲區(qū)或點(diǎn)云數(shù)據(jù)質(zhì)量較差的問(wèn)題，所以在DEM擬合過(guò)程中會(huì)出現(xiàn)較明顯的高程值計(jì)算誤差。將2種測(cè)量方法的試驗(yàn)中期對(duì)溝深的測(cè)量結(jié)果進(jìn)行對(duì)比分析發(fā)現(xiàn)，同時(shí)試驗(yàn)中期攝影測(cè)量計(jì)算的最大溝深平均相對(duì)誤差為6.8%，小于攝影測(cè)量計(jì)算的平均相對(duì)誤差（7.2%），說(shuō)明近景攝影測(cè)量技術(shù)在“窄深型”細(xì)溝的深度量測(cè)方面要優(yōu)于三維激光掃描技術(shù)。

2）細(xì)溝總表面積測(cè)量結(jié)果對(duì)比

從圖5中可以看出，2種非接觸式測(cè)量方法的測(cè)量值之間并無(wú)明顯差異，2種量測(cè)方法所得到的細(xì)溝總表面積數(shù)據(jù)呈相同的變化規(guī)律。細(xì)溝總表面積的發(fā)展演變與細(xì)溝發(fā)展密切相關(guān)。在溝壁擴(kuò)張侵蝕未出現(xiàn)之前，細(xì)溝面積僅受到溯源侵蝕的影響，細(xì)溝面積的增加主要來(lái)自于溝頭的推進(jìn)，而由于溝頭寬度較窄，所以細(xì)溝面積增長(zhǎng)較慢。溝壁擴(kuò)張侵蝕的出現(xiàn)引起細(xì)溝自上而下的寬度增加，細(xì)溝總面積增速加快。當(dāng)?shù)郊?xì)溝成熟階段，由于溝壁擴(kuò)張侵蝕占據(jù)主導(dǎo)地位，整條細(xì)溝的溝壁兩側(cè)都受到了更為強(qiáng)烈的侵蝕作用，所以細(xì)溝總表面積的增長(zhǎng)速度進(jìn)一步加快。

圖4 近景攝影測(cè)量和三維激光掃描技術(shù)測(cè)量結(jié)果對(duì)比

圖5 細(xì)溝總表面積變化

3）坡面微地形橫截面模擬對(duì)比

坡面微地形是用于描述坡面起伏變化的特征參數(shù)[33-34]。一方面地表微地形會(huì)隨著雨滴打擊和徑流沖擊引起的土壤分離而不斷發(fā)生變化，另一方面微地形的變化又會(huì)影響到土壤侵蝕產(chǎn)流產(chǎn)沙的變化。所以精確測(cè)量地表微地形對(duì)研究坡面侵蝕過(guò)程有重要意義。攝影測(cè)量和三維激光掃描在地表微地形模擬和測(cè)量中具有效率高、精度高等特點(diǎn)，近年來(lái)被廣泛應(yīng)用于坡面微地形研究中。然而2種方法在對(duì)高礫石含量的土壤進(jìn)行監(jiān)測(cè)時(shí)，模擬和測(cè)量精度是否會(huì)受礫石覆蓋層的影響，目前尚不明確。因此，在該部分本文基于2種非接觸式測(cè)量方法的高精度DEM利用ArcGIS 10.2的3D分析工具隨機(jī)提取坡面不同時(shí)間段的橫截面進(jìn)行對(duì)比分析，從而對(duì)2種測(cè)量方法的適用性進(jìn)行評(píng)價(jià)。

從圖6可以看出在坡面模擬方面，2種量測(cè)方法的模擬結(jié)果均呈現(xiàn)出前期粗糙，后期逐漸平滑的變化趨勢(shì)。這是因?yàn)樵谠囼?yàn)前期，由于降雨歷時(shí)短，坡面受侵蝕程度較小，坡面表面較少的土壤被帶走，較少的石英顆粒被滯留在坡面表面，石英顆粒覆蓋率較低，分散的石英顆粒同周圍的土壤呈現(xiàn)出明顯的高差。因此，坡面表面的橫截面抖動(dòng)較大。然而隨著降雨歷時(shí)的不斷增加，坡面土壤受侵蝕越來(lái)越嚴(yán)重，被降雨和徑流篩選出來(lái)的石英顆粒也越來(lái)越多，坡面石英覆蓋率不斷提升，石英顆粒同周圍的石英顆粒并無(wú)明顯的高差，所以坡面橫截面逐漸平滑。對(duì)比結(jié)果發(fā)現(xiàn)，近景攝影測(cè)量得到的坡面橫截面較為平滑，三維激光掃描技術(shù)所得到的坡面橫截面較為粗糙。這可能是因?yàn)槿S激光掃描儀在掃描時(shí)，儀器發(fā)射的高能激光受到了坡面石英顆粒立較強(qiáng)的反射，出現(xiàn)了較明顯的測(cè)量誤差?？偟膩?lái)看，依據(jù)近景攝影測(cè)量技術(shù)所得到的坡面模擬結(jié)果要優(yōu)于三維激光掃描技術(shù)。

圖6 不同降雨歷時(shí)坡面微地形橫截面

3 討 論

本文研究結(jié)果表明，坡面侵蝕演變呈現(xiàn)出明顯的階段特性，按主導(dǎo)侵蝕方式的不同，坡面侵蝕過(guò)程可劃分為片蝕階段、細(xì)溝發(fā)育階段和細(xì)溝成熟階段[35]。在片蝕階段（0～5 h），坡面侵蝕發(fā)育緩慢，坡面侵蝕以片蝕為主。在細(xì)溝發(fā)育階段（5～70 h），坡面侵蝕發(fā)育活躍，細(xì)溝發(fā)育受溯源侵蝕、下切侵蝕和溝壁擴(kuò)張侵蝕的共同作用影響。在細(xì)溝成熟階段（70～100 h），坡面侵蝕發(fā)育穩(wěn)定，細(xì)溝侵蝕以溝壁擴(kuò)張侵蝕為主。

細(xì)溝形態(tài)參數(shù)的演變?cè)诓煌那治g階段呈現(xiàn)出不同的特征，細(xì)溝長(zhǎng)度在細(xì)溝發(fā)育初期增長(zhǎng)速度最快，隨著降雨歷時(shí)的不斷增加，溝頭溯源侵蝕作用逐漸減弱，細(xì)溝長(zhǎng)度增長(zhǎng)速度逐漸變緩；細(xì)溝寬度在細(xì)溝發(fā)育初期增長(zhǎng)較慢，隨著溝壁擴(kuò)張侵蝕加劇，細(xì)溝寬度增長(zhǎng)不斷加快；細(xì)溝深度在細(xì)溝發(fā)育初期增長(zhǎng)速度較慢，在細(xì)溝發(fā)育中期增加最快，在細(xì)溝成熟階段由于溝底下切侵蝕逐漸消失，細(xì)溝深度增長(zhǎng)幾乎停止。不同坡位處的細(xì)溝形態(tài)呈現(xiàn)出不同的演變過(guò)程，坡上位置的細(xì)溝寬度和深度增長(zhǎng)最為迅速，坡中和坡下次之。對(duì)細(xì)溝的寬深比變化進(jìn)行分析發(fā)現(xiàn)，細(xì)溝底部（坡面50 cm）經(jīng)歷了一個(gè)“寬淺型-窄深型-寬淺型”的演變過(guò)程，細(xì)溝中部（坡面150 cm）始終為“窄深型”，細(xì)溝頭部（坡面250 cm）經(jīng)歷了一個(gè)“窄深型-寬淺型”的演變過(guò)程。由于溝壁擴(kuò)張侵蝕的出現(xiàn)和溝壁擴(kuò)張侵蝕對(duì)細(xì)溝侵蝕的影響作用不斷加強(qiáng)，細(xì)溝總表面積增長(zhǎng)加快。

通過(guò)對(duì)2種測(cè)量技術(shù)的精度對(duì)比，發(fā)現(xiàn)兩者在坡面侵蝕質(zhì)量測(cè)量方面有很好的適用性，其中坡面侵蝕量最大相對(duì)誤差為?16.82%，遠(yuǎn)小于Jiang等[12]利用近景攝影測(cè)定的坡面侵蝕量最大誤差（34.57%），同時(shí)在本研究中發(fā)現(xiàn)在細(xì)溝長(zhǎng)度量測(cè)方面，2種測(cè)量技術(shù)與直尺測(cè)量的結(jié)果差別較小，測(cè)量精度較高；在細(xì)溝寬度量測(cè)上，2種技術(shù)要優(yōu)于直尺測(cè)量；對(duì)于侵蝕量和深度量測(cè)，2種方法測(cè)量誤差均呈前期和末期小，中期大的特點(diǎn)，總體來(lái)說(shuō)近景攝影測(cè)量技術(shù)的測(cè)量精度要高于三維激光掃描技術(shù)，這是由于中后期細(xì)溝較深且窄，三維激光掃描存在較多的數(shù)據(jù)采集盲區(qū)[11]。李俊利等[14]也同樣發(fā)現(xiàn)近景攝影測(cè)量在溝壁和起伏區(qū)域重建的點(diǎn)云效果要優(yōu)于激光掃描。在細(xì)溝總表面積測(cè)量方面，2種方法的測(cè)量精度并無(wú)明顯差異。雖然2種測(cè)量技術(shù)的微地形模擬效果在試驗(yàn)前期較差，但隨著坡面侵蝕的發(fā)育，土壤本身所含的石英顆粒在坡面的不斷富集，模擬效果逐漸轉(zhuǎn)好，近景攝影測(cè)量技術(shù)的微地形模擬效果要優(yōu)于三維激光掃描技術(shù)。

綜上，近景攝影測(cè)量技術(shù)對(duì)于“窄深型”細(xì)溝的量測(cè)方面要優(yōu)于三維激光掃描，同時(shí)考慮到儀器設(shè)備價(jià)格、掃描采集靈活性等方面因素，近景攝影測(cè)量技術(shù)比三維激光掃描技術(shù)更適合于室內(nèi)人工模擬降雨條件下的坡面侵蝕演變研究。但是由于研究是基于室內(nèi)小規(guī)模坡面侵蝕試驗(yàn)展開(kāi)的。在野外中大規(guī)模坡面侵蝕測(cè)量時(shí)，2種方法的適用性、精確性是否如室內(nèi)試驗(yàn)一致，還需要進(jìn)一步驗(yàn)證。

4 結(jié) 論

本文基于激光掃描和攝影測(cè)量在坡面侵蝕演變過(guò)程的應(yīng)用研究，研究結(jié)論如下：

1）坡面侵蝕演變呈現(xiàn)出明顯的階段特性，按主導(dǎo)侵蝕方式的不同坡面侵蝕過(guò)程可分為：片蝕階段、細(xì)溝發(fā)育階段和細(xì)溝成熟階段，同時(shí)細(xì)溝形態(tài)參數(shù)的演變?cè)诓煌那治g階段呈現(xiàn)出不同的特征；

2）與實(shí)測(cè)的侵蝕量進(jìn)行對(duì)比發(fā)現(xiàn)，近景攝影測(cè)量技術(shù)的平均絕對(duì)誤差為8.84%，三維激光掃描技術(shù)的平均絕對(duì)誤差為8.97%，2種測(cè)量方法均能對(duì)坡面侵蝕產(chǎn)沙過(guò)程進(jìn)行監(jiān)測(cè)；

3）近景攝影測(cè)量技術(shù)在坡面侵蝕產(chǎn)沙監(jiān)測(cè)、細(xì)溝深度測(cè)量和坡面微地形模擬方面要優(yōu)于三維激光掃描技術(shù)，在細(xì)溝發(fā)育階段，細(xì)溝發(fā)育主要以溯源侵蝕、下切侵蝕為主，易形成“窄深型”細(xì)溝，這是由于三維激光掃描存在較多的數(shù)據(jù)采集盲區(qū)所導(dǎo)致的。

[1] Phinzi K, Ngetar N S. The assessment of water-borne erosion at catchment level using GIS-based RUSLE and remote sensing: A review[J]. International Soil and Water Conservation Research, 2019, 7(1): 27-46.

[2] Ni S M, Zhang D Q, Cai C F, et al.. Exploring rainfall kinetic energy induced erosion behavior and sediment sorting for a coarse-textured granite derived soil of south China[J]. Soil & Tillage Research, 2021, 208: 104915.

[3] Ni S M, Hui W, Wilson G V, et al. A simulated study of surface morphological evolution on coarse-textured soils under intermittent rainfall events[J]. Catena, 2022, 208: 105767.

[4] Shen H O, Zheng F L, Wen L L, et al. An experimental study of rill erosion and morphology[J]. Geomorphology, 2015, 231: 193-201.

[5] 覃超，何超，鄭粉莉，等. 黃土坡面細(xì)溝溝頭溯源侵蝕的量化研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(6)：160-167.

Qin Chao, He Chao, Zheng Fenli, et al. Quantitative research of rill head advancing process on loessial hillslope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(6): 160-167. (in Chinese with English abstract)

[6] Dong Y Q, Li F, Zhang Q W, et al. Determining ephemeral gully erosion process with the volume replacementmethod[J]. Catena, 2015, 131: 119-124.

[7] Casalí J, Loizu J, Campo M A, et al. Accuracy of methods for field assessment of rill and ephemeral gully erosion[J]. Catena, 2006, 67(2): 128-138.

[8] 徐錫蒙，鄭粉莉，武敏. 雨強(qiáng)和坡度對(duì)黃土陡坡地淺溝形態(tài)特征影響的定量研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(11)：124-132.

Xu Ximeng, Zheng Fenli, Wu Min. Quantification study of rainfall intensity and slope gradient impacts on ephemeral gully morphological characteristic on steep loessial hillslope[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(11): 124-132. (in Chinese with English abstract)

[9] Prosdocimi M, Burguet M, Di Prima S, et al. Rainfall simulation and Structure-from-Motion photogrammetry for the analysis of soil water erosion in Mediterranean vineyards[J]. Science of The Total Environment, 2017, 574: 204-215.

[10] 陶浩然，陳權(quán)，李震，等. 近景攝影測(cè)量提高裸露地表粗糙度測(cè)量精度[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(15)：162-167.

Tao Haoran, Chen Quan, Li Zhen, et al. Improvement of soil surface roughness measurement accuracy by close-range photogrammetry[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(15): 162-167. (in Chinese with English abstract)

[11] Gao C, Li P, Hu J, et al. Development of gully erosion processes: A 3D investigation based on field scouring experiments and laser scanning[J]. Remote Sensing of Environment, 2021, 265: 112683.

[12] Jiang Y, Shi H, Wen Z, et al. The dynamic process of slope rill erosion analyzed with a digital close range photogrammetry observation system under laboratory conditions[J]. Geomorphology, 2020, 350: 106893.

[13] 明旭輝，王瑄，盛思遠(yuǎn)，等. 基于近景攝影測(cè)量的徑流沖刷條件下凍融坡面侵蝕產(chǎn)沙過(guò)程[J]. 水土保持學(xué)報(bào)，2019，33(2)：39-42，48.

Ming Xuhui, Wang Xuan, Sheng Siyuan, et al. Erosion and sediment yield process on freeze-thaw slope based on close range photogrammetry under runoff scouring conditions[J]. Journal of Soil and Water Conservation, 2019, 33(2): 39-42, 48. (in Chinese with English abstract)

[14] 李俊利，李斌兵，柳方明，等. 利用照片重建技術(shù)生成坡面侵蝕溝三維模型[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2015，31(1)：125-132.

Li Junli, Li Binbing, Liu Fangming, et al. Generating 3D model of slope eroded gully based on photo reconstruction technique[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(1): 125-132. (in Chinese with English abstract)

[15] 張甘霖, 龔子同. 土壤調(diào)查實(shí)驗(yàn)室分析方法[M]. 北京：科學(xué)出版社，2012.

[16] 楊瑋，蘭紅，李民贊，等. 基于圖像處理和GBRT模型的表土層土壤容重預(yù)測(cè)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2020，51(9)：193-200.

Yang Wei, Lan Hong, Li Minzan, et al. Prediction of top soil layer bulk density based on image processing and gradient boosting regression tree mode[J]. Transactions of the Chinese Society for Agricultural Machinery, 2020, 51(9): 193-200. (in Chinese with English abstract).

[17] Guo Z L, Ma M J, Cai C F, et al. Combined effects of simulated rainfall and overland flow on sediment and solute transport in hillslope erosion[J]. Journal of Soils and Sediments, 2018, 18: 1120-1132.

[18] Ni S M, Zhang D Q, Wen H, et al.. Erosion processes and features for a coarse-textured soil with different horizons: A laboratory simulation[J]. Journal of Soils and Sediments, 2020, 20: 2997-3012.

[19] Wang L, Shi Z H, Wang J, et al. Rainfall kinetic energy controlling erosion processes and sediment sorting on steep hillslopes: a case study of clay loam soil from the Loess Plateau. China[J]. Journal of Hydrology, 2014, 512: 168-176.

[20] 方乾，趙龍山，樊春華，等. 石漠化區(qū)薄土層坡地雨水轉(zhuǎn)化及土壤侵蝕影響因素分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(8)：88-97.

Fang Qian, Zhao Longshan, Fan Chunhua, et al. Influencing factors of rainwater transformation and soil erosion in thin soil hillslope of rock desertification regions[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(8): 88-97. (in Chinese with English abstract)

[21] 郭萍，夏振堯，高峰，等. 香根草植物籬對(duì)三峽庫(kù)區(qū)坡地紫色土侵蝕的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(19)：105-112.

Guo Ping, Xia Zhenyao, Gao Feng, et al. Effects of Vetiveria zizanioides hedgerow on the erosion of purple soil of slope land in the Three Gorges Reservoir Area of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(19): 105-112. (in Chinese with English abstract)

[22] 吳淑芳，劉勃洋，雷琪，等. 基于三維重建技術(shù)的坡面細(xì)溝侵蝕演變過(guò)程研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(9)：114-120.

Wu Shufang, Liu Boyang, Lei Qi, et al. Evolution process of slope rill erosion based on 3D photo reconstruction technique[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(9): 114-120. (in Chinese with English abstract)

[23] 王睿，李鵬，韓建純，等. 溝頭高度和土壤質(zhì)地對(duì)細(xì)溝溯源侵蝕特征和形態(tài)發(fā)育的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(10)：91-99.

Wang Rui, Li Peng, Han Jianchun, et al. Effects of gully head height and soil texture on headward erosion characteristics and topography evolution[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(10): 91-99. (in Chinese with English abstract)

[24] Qin C, Zheng F L, Xu X M, et al. A laboratory study on rill network development and morphological characteristics on loessial hillslope[J]. Journal of Soils and Sediments, 2018, 18(4): 1679-1690.

[25] Bewket W, Sterk G. Assessment of soil erosion in cultivated fields using a survey methodology for rills in the Chemoga watershed, Ethiopia[J]. Agriculture, Ecosystems and Environment, 2003, 97: 81-93.

[26] 郝好鑫，郭忠錄，王先舟，等. 降雨和徑流條件下紅壤坡面細(xì)溝侵蝕過(guò)程[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(8)：134-140.

Hao Haoxin, Guo Zhonglu, Wang Xianzhou, et al. Rill erosion process on red soil slope under interaction of rainfall and scouring flow[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(8): 134-140. (in Chinese with English abstract)

[27] Jiang F S, Huang Y H, Wang M K, et al. Sediment selective behaviors of colluvial deposit materials on steep slopes and under heavy rainfall conditions in southeastern China[J]. Energy and Ecology, 2018, 73: 133-142.

[28] Istanbulluoglu E. Implications of bank failures and fluvial erosion for gully development: Field observations and modeling[J]. Journal of Geophysical Research, 2005, 110: F01014.

[29] Thomas J T, Iverson N R, Burkart M R. Bank-collapse processes in a valley-bottom gully, western Iowa[J]. Earth Surface Processes and Landforms, 2009, 34(1): 109-122.

[30] Martinez-Casasnovas J A, Ramos M C, Poesen J. Assessment of sidewall erosion in large gullies using multi-temporal DEMs and logistic regression analysis[J]. Geomorphology, 2014, 58(1/3/4): 305-321.

[31] Chen A, Zhang D, Peng H, et al. Experimental study on the development of collapse of overhanging layers of gully in Yuanmou Valley, China[J]. Catena, 2013, 109: 177-185.

[32] Qin C, Zheng F, Wells R R, et al. A laboratory study of channel sidewall expansion in upland concentrated flows[J]. Soil and Tillage Research, 2018, 178: 22-31.

[33] Guo M, Shi H, Zhao J, et al. Digital close range photogrammetry for the study of rill development at flume scale[J]. Catena, 2016, 143: 265-274.

[34] Li S, Li Q Q, Chen J, et al. Application of 3D laser image scanning technology and cellular automata model in the prediction of the dynamic process of rill erosion[J]. Remote Sensing, 2021, 13: 2586.

[35] Quan X, He J J, Cai Q G, et al. Soil erosion and deposition characteristics of slope surfaces for two loess soils using indoor simulated rainfall experiment[J]. Soil & Tillage Research, 2020, 204: 104714.

Application of laser scanning and photogrammetry in the evolution process of slope erosion

Luo Bin1, Zhang Yong2, Zhang Zhiwei1, Ni Shimin1, Zhang Xin2, Wang Junguang1※

(1.,,430070,; 2,,430010)

Soil erosion on a slope is one of the most serious soil degradation that caused by the complex dynamic process of water draining down the slope. Quantitative monitoring of the slope erosion process is of great practical significance to clarify the erosion mechanism for the better models in the recent years. A non-contact measurement has been widely used, such as the laser scanning or photogrammetry, due to the high efficiency and accuracy. However, it is still lacking on the accuracy and applicability evaluation in the process of slope erosion, especially for the rill morphology characteristics. In this study, the laser scanning and photogrammetry were selected to monitor the evolution of soil slope erosion under the indoor artificial rainfall, in order to quantitatively explore their accuracy and applicability. The red soil was collected with the granite parent material from the Changting County, Hubei Province, China. The test material was then selected as the leaching, sedimentary, and parent material layer from the bottom to the top at the sampling point. The soil trough was used with the size of 4.0 m (length) × 2.0 m (width) × 0.60 m (height), particularly with the depth of soil filling of 0.50 m. The soil was then filled with the 0.15 m thick parent material layer, 0.20 m sedimentary layer, and 0.15 cm leaching layer before the test. The soil tank was subjected to the multiple cycles of drying and wetting. The static settlement was set at least two months, in order to make the soil layers and particles as close as possible to natural conditions. The whole duration was 1 h for the single-field simulated rainfall or runoff test after the runoff on the slope surface. The total rainfall was 100 h, where the interval was 24 h between each two intermittent rainfall. The rain intensity was also set to (90±5) mm/ h. The sediment ocean was then collected during the rainfall. After that, the 3D laser scanning and close-range photogrammetry were employed to record the micro-topographic and surface morphology of the slope at each rainfall interval. As such, the digital point cloud was obtained on the slope surface after the rainfall. At the same time, the length, width and depth of each rill on the slope were measured at the intervals of 10 cm, in order to determine the evolution of the morphological characteristics of the slope rill. The characteristic parameters of rill morphology were calculated to collect the runoff sediment and measurement. A comparison was also made on the laser scanning and photogrammetry technology in the erosion monitoring test. The results show that: 1) The slope erosion process was divided into the sheet erosion, rill development and maturity stage, according to the different dominant erosion. 2) Both non-contact measurement methods performed better to accurately measure the slope erosion. Specifically, the maximum relative error was -16.82% for the process of sediment production. The better applicability was also achieved in the measurement of slope erosion quality. 3) The close-range photogrammetry technology was superior to the 3D laser scanning in the monitoring slope erosion and sediment production, rill depth measurement, and slope micro-topography simulation. Therefore, the digital expression was realized for the evolution process of slope erosion, in order to promote the high efficiency and high precision of slope erosion monitoring. The finding can also provide the strong reference in the technology selection for the soil slope erosion.

slope; erosion; soils; close-range photogrammetry technology; 3D laser scanning technology; comparison

10.11975/j.issn.1002-6819.2022.17.011

S157.1

A

1002-6819(2022)-17-0101-09

羅斌，張勇，張志偉，等. 激光掃描和攝影測(cè)量在坡面侵蝕演變過(guò)程的適用性[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(17)：101-109.doi：10.11975/j.issn.1002-6819.2022.17.011 http://www.tcsae.org

Luo Bin, Zhang Yong, Zhang Zhiwei, et al. Application of laser scanning and photogrammetry in the evolution process of slope erosion[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(17): 101-109. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2022.17.011 http://www.tcsae.org

2022-04-25

2022-07-10

國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2021YFD1500703）；國(guó)家自然科學(xué)基金（42177317）

羅斌，研究方向?yàn)橥寥狼治g與水土保持。Email：Luobin2021@163.com

王軍光，博士，副教授，研究方向?yàn)橥寥狼治g機(jī)理。Email：jgwang@mail.hzau.edu.cn