利用探地雷達(dá)技術(shù)獲取與土壤鹽漬化程度相關(guān)的信息,農(nóng)業(yè)基礎(chǔ)科學(xué)論文土壤鹽漬化是干旱、半干旱地區(qū)主要的土地退化問(wèn)題之一[1-2],造成了農(nóng)作物產(chǎn)量下降、區(qū)域生態(tài)環(huán)境改變.準(zhǔn)確地獲取土壤鹽漬化鹽分信息,對(duì)于干旱區(qū)鹽漬土治理和利用,促進(jìn)干旱區(qū)農(nóng)業(yè)可持續(xù)發(fā)展具有重要意義[3-4].國(guó)內(nèi)外在鹽漬化土壤探測(cè)技術(shù)的研究中,對(duì)可見(jiàn)光、紅外和高頻微波波段的星載遙感探測(cè)技術(shù)研究較多.當(dāng)前應(yīng)用的星載雷達(dá)都是高頻雷達(dá),其對(duì)土壤的穿透深度有限,和其他非微波波段的星載遙感技術(shù)一樣,探測(cè)的主要是表層土壤的鹽漬化狀況,而由于土壤鹽漬化是一定深度土壤水鹽運(yùn)動(dòng)的結(jié)果,表層土壤鹽漬化狀況并不能完全揭示土壤鹽漬化規(guī)律,如非表聚型鹽漬化土壤等.近年,地面電磁技術(shù)等用于土壤鹽漬化程度的探測(cè)中[5],其主要途徑是建立探測(cè)信號(hào)和鹽漬化土壤實(shí)驗(yàn)室分析數(shù)據(jù)的定量關(guān)系模型,當(dāng)前用得較多的大地電導(dǎo)率儀,固然能夠探測(cè)一定深度的土壤,但一般不能直接反映土壤鹽分狀況.為了較全面、準(zhǔn)確地定量獲得土壤鹽分信息,需要發(fā)展新的土壤鹽分探測(cè)技術(shù)和方式方法.探地雷達(dá)技術(shù)是基于電磁原理的地面探測(cè)技術(shù),由于具有較高的分辨率、能夠?qū)崟r(shí)成像和探測(cè)效率較高,在土壤水分探測(cè)中得到了較多的應(yīng)用,也用于土壤鹽漬化狀況定性探測(cè)中[6].由于影響探地雷達(dá)探測(cè)結(jié)果的土壤參數(shù)較多,當(dāng)前在利用探地雷達(dá)進(jìn)行土壤參數(shù)探測(cè)多是在一定控制條件下的試驗(yàn)性研究[7].根據(jù)當(dāng)前國(guó)內(nèi)外土壤鹽漬化探測(cè)研究的現(xiàn)在狀況和存在的問(wèn)題,針對(duì)具有一定深度的鹽漬化土壤,研究了怎樣利用探地雷達(dá)技術(shù)定量獲取與土壤鹽漬化程度有關(guān)的含鹽量和電導(dǎo)率信息,以期為鹽漬化土壤定量探測(cè)和防治提供新的技術(shù)手段和方式方法.1研究區(qū)大概情況研究區(qū)位于塔里木盆地的中北部,新疆阿克蘇地區(qū)渭干河-庫(kù)車河流域的下游,中部天山南麓,地勢(shì)北高南低,屬大陸性暖溫帶極端干旱氣候,平均蒸發(fā)量2124mm/年,平均降水量52mm/年,蒸降比約40∶1.研究區(qū)以鹽巖、石膏為主體的鹽類沉積規(guī)模大,地下水位和礦化度偏高,土壤鹽漬化嚴(yán)重且普遍.土壤類型主要為亞砂土和亞黏土.土壤pH值約為8,土壤含鹽主要成份為氯化物.2研究方式方法當(dāng)前研究較多的是單層介質(zhì)的電磁波傳播速度,而準(zhǔn)確求取多層介質(zhì)中每層介質(zhì)的電磁波傳播速度較為困難.為了利用探地雷達(dá)電磁波傳播速度表示出土壤成分組成情況,結(jié)合探測(cè)的實(shí)際情況,將每層土壤的探地雷達(dá)電磁波傳播速度利用每層土壤的厚度與電磁波在該層傳播的雙程旅行時(shí)間的比值計(jì)算,即v=2h/t,式中,v為探地雷達(dá)電磁波在每層土壤中的傳播速度,h為每層土壤的厚度,t為電磁波在該層土壤的雙程旅行時(shí)間.研究中使用兩種儀器設(shè)備:拉脫維亞雷達(dá)系統(tǒng)公司zond-12e探地雷達(dá),使用100MHz屏蔽天線,用于土壤鹽漬化信息探測(cè);基于電磁感應(yīng)原理的土壤溫度、濕度和電導(dǎo)率探測(cè)儀,其利用的電磁頻率為50MHz,用于實(shí)時(shí)測(cè)量土壤溫度、濕度和電導(dǎo)率.根據(jù)屢次實(shí)地調(diào)查結(jié)果,在研究區(qū)選擇鹽漬化狀況不同的區(qū)域建立了5個(gè)模型,對(duì)每個(gè)區(qū)域設(shè)定邊長(zhǎng)30m的正方形樣方,保證每個(gè)樣方內(nèi)的土壤狀況基本一致,且和其他樣方有明顯差異.利用探地雷達(dá)進(jìn)行數(shù)據(jù)采集時(shí),首先沿著樣方的4條邊界和2條對(duì)角線,利用100MHz屏蔽天線采用剖面法進(jìn)行探測(cè),并記錄各條測(cè)線的數(shù)據(jù),然后,在樣方對(duì)角線交點(diǎn)處,即樣方中心點(diǎn)處進(jìn)行定點(diǎn)探測(cè),并把該點(diǎn)的測(cè)量數(shù)據(jù)作為代表每個(gè)樣方鹽漬化情況的樣本數(shù)據(jù);探地雷達(dá)儀器系統(tǒng)記錄完中心點(diǎn)的數(shù)據(jù)后,對(duì)探測(cè)數(shù)據(jù)進(jìn)行處理和判讀.從處理后的探地雷達(dá)剖面圖上可看出每個(gè)樣方中心點(diǎn)的探測(cè)剖面都有明顯的分層情況.圖1是樣方1的探地雷達(dá)探測(cè)成果,能夠看出,鹽漬化土壤在100MHz探地雷達(dá)探測(cè)剖面上有明顯的分層,能被清楚明晰辨識(shí)的分層情況的深度到達(dá)了近3m,淺的也到達(dá)了近0.8m;這是當(dāng)前星載遙感難以直接探測(cè)的深度,講明低頻探地雷達(dá)在鹽漬化土壤探測(cè)深度上具有較大的優(yōu)勢(shì).利用探地雷達(dá)數(shù)據(jù)處理軟件系統(tǒng),獲取每個(gè)剖面的位置和厚度情況,根據(jù)數(shù)據(jù)判讀結(jié)果開(kāi)挖土壤剖面,利用量尺在土壤剖面上定位判讀每層土壤的深度和厚度,自下而上采集每層的土壤樣本,帶回實(shí)驗(yàn)室進(jìn)行化學(xué)成分分析.在采集每個(gè)樣方的每層土壤樣本的同時(shí),利用50MHz的土壤溫度、濕度和電導(dǎo)率探測(cè)儀測(cè)定每層土壤的電導(dǎo)率,為了獲得較準(zhǔn)確的數(shù)據(jù),將每層進(jìn)行屢次測(cè)量的數(shù)據(jù)平均值作為該層真實(shí)電導(dǎo)率值.5個(gè)樣方都作類似的檢測(cè).將采集的土壤樣品帶回實(shí)驗(yàn)室自然風(fēng)干、磨碎,過(guò)0.5mm孔徑的篩子,按土水比1∶5的比例提取浸提液,分別測(cè)定土壤的含鹽量、電導(dǎo)率、礦化度、pH值、化學(xué)離子等參數(shù).3數(shù)據(jù)計(jì)算與結(jié)果分析將每個(gè)樣方每層土壤的含鹽量、電導(dǎo)率分別與之對(duì)應(yīng)的電磁波傳播速度進(jìn)行pearson相關(guān)性分析,結(jié)果見(jiàn)表1.表中土壤剖面深度指的是每個(gè)樣方能被容易辨識(shí)的分層情況的深度,也是土壤標(biāo)本的采樣深度.根據(jù)中國(guó)水利部行業(yè)標(biāo)準(zhǔn)(土壤侵蝕分類分級(jí)標(biāo)準(zhǔn)〕(SL190-96),表1中樣方1屬于鹽土,樣方2和樣方5屬于強(qiáng)度鹽漬化土壤,樣方3屬于中度鹽漬化土壤,樣方4屬于輕度鹽漬化土壤.從表1中各個(gè)樣方每層的電磁波傳播速度和土壤含鹽量或土壤電導(dǎo)率相關(guān)分析結(jié)果得出:電磁波傳播速度和樣方1、樣方2、樣方3的土壤含鹽量有顯著的相關(guān)性,與樣方3的相關(guān)系數(shù)甚至到達(dá)了0.932;樣方4和樣方5的電磁波速與土壤含鹽量的相關(guān)性不顯著;同時(shí),樣方1、樣方4和樣方5的波速與土壤電導(dǎo)率有顯著的相關(guān)性,而樣方2和樣方3的兩者相關(guān)系數(shù)較低,相關(guān)性不顯著;除此之外,樣方1的電磁波傳播速度與土壤含鹽量和電導(dǎo)率都有顯著的相關(guān)性,并且,其與土壤含鹽量的相關(guān)系數(shù)大于其與土壤電導(dǎo)率的相關(guān)系數(shù).和土壤含鹽量相關(guān)性不顯著的樣方,都和土壤電導(dǎo)率有顯著的相關(guān)性,講明探地雷達(dá)剖面圖中分層現(xiàn)象的主要影響因素是土壤鹽分或土壤電導(dǎo)率.探地雷達(dá)探測(cè)的剖面圖上顯示的是研究區(qū)土壤鹽分含量或電導(dǎo)率的分布狀況.對(duì)表1中各個(gè)樣方相關(guān)性顯著的參數(shù),以含鹽量或電導(dǎo)率為因變量,以電磁波傳播速度為自變量,進(jìn)行三次多項(xiàng)式回歸分析,結(jié)果如此圖2所示.樣方1、樣方2、樣方3的電磁波傳播速度和土壤含鹽量的回歸擬合結(jié)果的決定系數(shù)R2都在0.95以上,講明電磁波傳播速度能夠很好地反映這三個(gè)樣方的土壤含鹽量狀況.樣方1、樣方4和樣方5電磁波傳播速度和土壤電導(dǎo)率回歸結(jié)果的決定系數(shù)R2分別為0.3896、0.5958和0.3498,講明電磁波傳播速度也能反映這三個(gè)樣方的土壤電導(dǎo)率狀況.探測(cè)的5個(gè)不同程度的鹽漬化土壤樣方中,有3個(gè)土壤鹽分和探地雷達(dá)電磁波傳播速度有顯著相關(guān)性,其他2個(gè)樣方的土壤電導(dǎo)率和電磁波傳播速度有顯著相關(guān)性,講明探地雷達(dá)能夠用于土壤含鹽量定量探測(cè),但其存在一定的不確定性.4結(jié)論研究表示清楚,利用100MHz天線的探地雷達(dá)能夠較好地識(shí)別研究區(qū)鹽漬化土壤分層情況;探地雷達(dá)的電磁波傳播速度和土壤鹽分或電導(dǎo)率有顯著的相關(guān)性,利用低頻探地雷達(dá)能夠定量探測(cè)土壤鹽分或電導(dǎo)率垂直分布狀況,并且利用探地雷達(dá)電磁波傳播速度可估計(jì)土壤含鹽量或土壤電導(dǎo)率.利用探地雷達(dá)估計(jì)土壤含鹽量或土壤電導(dǎo)率具有一定的不確定性,或者講,利用低頻探地雷達(dá)進(jìn)行測(cè)鹽漬化土壤探測(cè)時(shí),影響探地雷達(dá)電磁波傳播速度的主導(dǎo)因素具有一定的不確定性.筆者只研究了100MHz天線的探地雷達(dá)在鹽漬化土壤探測(cè)中的應(yīng)用,有必要研究電磁波頻率對(duì)探測(cè)結(jié)果的影響,同時(shí),研究中只牽涉到干旱區(qū)較嚴(yán)重鹽漬化土壤的含鹽量或電導(dǎo)率參數(shù),沒(méi)有考慮土壤濕度參數(shù).因而,本研究主要是基于數(shù)學(xué)統(tǒng)計(jì)方式給出了鹽漬化土壤層探測(cè)的一個(gè)可能性的方式方法,該方式方法的適用性和深層次的機(jī)理問(wèn)題需要進(jìn)一步研究.以下為參考文獻(xiàn):[1]AmmariTG,TahhanR,AbubakerS,etal.SoilsalinitychangesintheJordanvalleypotentiallythreatensustainableirrigatedagri-culture[J].Pedosphere,2020,23(3):376-384.[2]AnilKM.Climatechangeimpact,mitigationandadaptationstrat-egiesforagriculturalandwaterresources,inGangaPlain(India)[J].MitigAdaptStrategGlobChange,2020,18:673-689.[3]YaoRJ,YangJS,ShaoHB.Accuracyanduncertaintyassess-mentongeostatisticalsimulationofsoilsalinityinacoastalfarm-landusingauxiliaryvariable[J].EnvironMonitAssess,2020,185:5151-5164.[4]AkhtarA,ShahbazK,NisarH,etal.Characterizingsoilsalinityinirrigatedagricultureusingaremotesensingapproach[J].PhysicsandChemistryoftheEarth,2020,55~57:43-52.[5]孫永猛,丁建麗,瞿娟,等.應(yīng)用電磁感應(yīng)和遙感的新疆綠洲區(qū)域尺度鹽漬土辨別[J].農(nóng)業(yè)工程學(xué)報(bào),2020,28(20):180-187.[6]薛建,曾昭發(fā),田剛,等.探地雷達(dá)在吉林西部地區(qū)探測(cè)土壤堿化層[J].物探與化探,2005,29(5):421-424.[7]LauroSE,MatteiE,BaronePM.Estimationofsubsurfacedielec-trictargetdepthforGPRplanetaryexploration:Laboratorymeas-urementsandmodeling[J].JournalofAppliedGeophysics,2020,93:93-100.[8]AndreaDM,ElenaP,EvertS.Early-timeGPRsignalattributestoestimatesoildielectricpermittivity:Atheoreticalstudy[J].IEEETransactionsonGeoscienceandRemoteSensing,2020,51:1643-1653.[9]MohammadR,MahmoudzadehA.Off-andon-groundGPRtech-niquesforfield-scalesoilmoisturemapping[J].Geo