• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于地統(tǒng)計分析的河套灌區(qū)土壤水鹽空間變異規(guī)律分析

    2022-09-28 07:46:58劉燕楠孫貫芳毛威成蕭堯朱焱楊金忠
    灌溉排水學報 2022年9期
    關(guān)鍵詞:含鹽量礦化度變異

    劉燕楠,孫貫芳,毛威,成蕭堯,朱焱,楊金忠

    基于地統(tǒng)計分析的河套灌區(qū)土壤水鹽空間變異規(guī)律分析

    劉燕楠1,孫貫芳2,毛威1,成蕭堯3,朱焱1,楊金忠1*

    (1.武漢大學水資源與水電工程科學國家重點實驗室,武漢 430072;2.西北農(nóng)林科技大學 水土保持研究所,陜西 楊凌 712100;3.中國市政工程中南市政設(shè)計研究總院有限公司,武漢 430101)

    【目的】研究干旱鹽漬化灌區(qū)不同深度土壤水分和鹽分的空間變異規(guī)律。【方法】以內(nèi)蒙古河套灌區(qū)隆勝井渠結(jié)合區(qū)為研究區(qū),通過野外取樣和實驗室分析結(jié)合,應(yīng)用經(jīng)典統(tǒng)計學和地統(tǒng)計學方法研究了根系層(0~60 cm)及深層土壤(60~120、120~180 cm)土壤含水率和含鹽量的空間變異性、空間格局及二者之間的關(guān)系,并探究了地下水埋深和地下水礦化度對土壤含鹽量的影響?!窘Y(jié)果】①土壤含水率在垂向上大小關(guān)系為0~60 cm<60~120 cm<120~180 cm,土壤含鹽量為0~60 cm>60~120 cm>120~180 cm。土壤含水率和含鹽量的變異系數(shù)和塊金系數(shù)自上而下逐層減小,變程逐層增大。根系層的土壤含水率和含鹽量均具有中等程度的空間相關(guān)性,深層土壤含水率和含鹽量的空間相關(guān)性較強。②各層土壤含水率的變異系數(shù)和塊金系數(shù)小于土壤含鹽量,變程大于土壤含鹽量,二者之間的空間格局存在顯著的正相關(guān)關(guān)系。③研究區(qū)各層土壤含鹽量均與地下水埋深呈指數(shù)相關(guān)關(guān)系,與地下水礦化度呈冪函數(shù)關(guān)系。當?shù)叵滤裆畲笥?.5 m時,總體土壤含鹽量較低,且各層含鹽量變化較小。【結(jié)論】研究區(qū)土壤鹽分呈表聚型,其空間特征和分布格局與土壤含水率有著密切的關(guān)系。調(diào)控地下水埋深是灌區(qū)控制土壤鹽漬化的關(guān)鍵措施,為有效控制土壤鹽漬化,建議地下水埋深控制在2.5 m以下。

    河套灌區(qū);土壤水鹽;地統(tǒng)計學;空間變異;地下水環(huán)境

    0 引言

    【研究意義】全球有100多個國家面臨土壤鹽漬化問題,土壤鹽漬化面積占世界陸地面積的6%以上[1]。在干旱和半干旱地區(qū),近50%的灌溉土地存在一定程度的土壤鹽漬化問題[2]。土壤鹽漬化導(dǎo)致農(nóng)業(yè)生產(chǎn)力顯著下降,嚴重制約著灌區(qū)的生態(tài)安全和可持續(xù)發(fā)展。土壤含鹽量受氣候、土質(zhì)、地下水環(huán)境等自然因素及耕作、灌溉排水等人類活動的共同影響,具有高度的空間變異性,其精確預(yù)測和治理的難度較高。因此,有必要明確土壤含鹽量的空間變異規(guī)律,探明其主要影響因素,為灌區(qū)土壤鹽漬化防治提供理論依據(jù)?!狙芯窟M展】國內(nèi)外學者利用地統(tǒng)計學方法研究了不同區(qū)域、不同時間、不同深度的土壤含鹽量的空間分布特征,并取得了一系列成果。Wang等[3]利用地統(tǒng)計學和Lorenz曲線方法,確定了三工河流域表層土壤(0~20 cm)含鹽量累積的空間特征與格局。Xie等[4]運用地統(tǒng)計學和近端感應(yīng)電磁法相結(jié)合的方法分析了長江口2006—2017年表層土壤含鹽量的時空變異特征,發(fā)現(xiàn)土壤含鹽量隨著與海洋和河流距離的減小而逐漸增加。劉洪波等[5]揭示了新疆第二師31團灌區(qū)的春季不同深度土層(0~20、20~40、80~100 cm)的土壤含鹽量以及鹽基離子的空間特征,結(jié)果表明不同深度土層土壤含鹽量及其空間變異性不同。婁婷等[6]在分析河套灌區(qū)典型灌域夏灌前后0~15 cm土層的土壤含鹽量和地下水埋深空間變異特征的基礎(chǔ)上,得出土壤含鹽量隨地下水埋深增大而減小的結(jié)論。張雪等[7]認為隨土層深度的增加,孔雀河流域土壤含鹽量逐漸降低,其空間相關(guān)性逐漸增大?!厩腥朦c】已有的研究結(jié)果肯定了地統(tǒng)計學在揭示土壤含鹽量的空間變異規(guī)律方面的可行性,但是目前的研究大多只關(guān)注土壤含鹽量的空間變異,同時研究不同深度土壤含水率和含鹽量的空間變異特征及相互關(guān)系,并探究地下水埋深及礦化度對土壤含鹽量的影響的報道較少?!緮M解決的關(guān)鍵問題】以內(nèi)蒙古河套灌區(qū)隆勝試驗區(qū)為研究區(qū),通過野外監(jiān)測試驗,應(yīng)用地統(tǒng)計學方法和GIS技術(shù),分析了不同深度土壤含水率和含鹽量的空間分布規(guī)律及相關(guān)關(guān)系,進一步分析地下水埋深和地下水礦化度對土壤含鹽量的影響,為干旱區(qū)土壤鹽漬化防控提供技術(shù)支持。

    1 材料與方法

    1.1 研究區(qū)概況

    隆勝試驗區(qū)位于內(nèi)蒙古自治區(qū)巴彥淖爾市河套灌區(qū)永濟灌域內(nèi)(107°28′E—107°32′E,40°51′E—107°32′E),地理位置如圖1所示,試驗區(qū)西起永濟干渠,東至東濟支渠,北為永剛分干溝,南臨永剛分干渠,西南至東北方向長約為15.5 km,西北至東南方向?qū)捈s為8.0 km,總控制面積8 219.75 hm2,灌溉面積5 560 hm2,海拔1 029~1 043 m。

    研究區(qū)屬于溫帶大陸氣候,干燥少雨,光照時間長,蒸發(fā)強烈,多年平均日照時間3 229.9 h,主要種植作物有玉米、葵花和瓜菜等。研究區(qū)鄰近臨河氣象站,年均降雨量148.87 mm(1991—2013年),主要集中在5—9月作物生育期,1957—2013年多年平均蒸發(fā)量(20 cm蒸發(fā)皿)為2 291.1 mm,年蒸發(fā)量是年降水量的15.4倍,屬于典型的無灌溉即無農(nóng)業(yè)的地區(qū)[8]。試驗區(qū)土壤質(zhì)地主要為壤土和砂土插花分布,表層基本為壤土,土壤體積質(zhì)量1.45~1.55 g/cm3,120 cm以下以砂土為主,大部分區(qū)域的地下水礦化度<3.0 g/L,年均地下水埋深1.88 m[9]。

    1.2 采樣點布設(shè)與數(shù)據(jù)測定

    1)采樣點布設(shè)。在試驗區(qū)內(nèi)布置17口地下水觀測井和40個土壤水鹽采樣點,涵蓋荒地、番茄地、葫蘆地、葵花地和玉米地(見圖1)。取樣時間為2020年4月末,以手持全球定位系統(tǒng)(GPS)測定實地坐標,采樣深度從表層到180 cm(或至地下水位),每20 cm為1層,共9層(0~20、20~40、40~60、60~80、80~100、100~120、120~140、140~160 cm和160~180 cm)。每個采樣點重復(fù)2次,重復(fù)的均值作為該采樣點的實測結(jié)果,共采集土壤樣本數(shù)據(jù)710個。

    2)數(shù)據(jù)測定。采用烘干法測定土壤含水率[10],同步測量土壤含鹽量,研磨的風干土樣過2 mm篩后利用電導(dǎo)率儀(上海雷磁電導(dǎo)率儀DDSJ-308F)測試土水比為1∶5的土壤浸提液的電導(dǎo)率,并根據(jù)換算公式計算土壤含鹽量[9]。

    圖1 隆勝研究區(qū)和采樣點分布示意

    1.3 數(shù)據(jù)分析方法

    1)半方差函數(shù)

    本文采用地統(tǒng)計學方法研究土壤含水率和含鹽量的空間變異性。半方差函數(shù)是地統(tǒng)計學的基本工具[11-12],以采樣點之間方差的一半表征區(qū)域化變量的空間變異結(jié)構(gòu),表達式如式(1):

    式中:()為滯后距離,等于點的對數(shù);(x)為點x處變量的實測值;(x+h)為與點x偏離處變量的實測值。本文中的半方差函數(shù)采用高斯模型進行擬合,如式(2):

    式中:為滯后距離;為拱高;0為塊金值;0+為基臺值;為變程,表示空間變量的最大相關(guān)距離。通常采用塊金系數(shù)0/0+表示隨機變異(灌溉、施肥、耕作等人為活動)占總變異的大小,一般認為0/0+<0.25時,變量的空間相關(guān)性強;0/0+在0.25~0.75時,變量的空間相關(guān)性中等;0/0+>0.75時,變量的空間相關(guān)性弱;如果該比值接近1,則說明該變量在整個尺度上的變異幾乎是隨機的。

    半方差函數(shù)的擬合精度采用決定系數(shù)2以及殘差平方和RSS來表征,其中2取值在0~1之間,該值越大說明擬合效果越好;RSS表示剩余誤差,該值越小模型的模擬效果越好。

    2)普通克里金插值

    克里金插值是地統(tǒng)計中最為常用的插值法,利用原始數(shù)據(jù)和半方差函數(shù)的結(jié)構(gòu)性,實現(xiàn)區(qū)域變量的最優(yōu)無偏估計[13]。本文使用克立金家族中最常用的普通克里金來插值得到研究區(qū)土壤含水率和含鹽量的空間分布情況。假設(shè)估計點0的鄰域內(nèi)共有個測點,即1,2,…,2,普通克里金插值公式如式(3):

    式中:(x)為估計點0附近選擇進行克里金插值的觀測值;Z(0)為0處的克里金估計值;為第個觀測值對估計點0的權(quán)重,由半方差函數(shù)建立的克里金方程組求得[14]。

    3)Pearson相關(guān)分析

    Pearson相關(guān)系數(shù)表示2組數(shù)據(jù)之間的相關(guān)程度,計算如式(4):

    表1 土壤含水率和含鹽量的統(tǒng)計特征值

    2 結(jié)果與分析

    2.1 土壤含水率和含鹽量的地統(tǒng)計特征

    土壤含水率和含鹽量按照每60 cm為1層,分3層(0~60、60~120 cm和120~180 cm)進行研究,其統(tǒng)計特征值如表1所示。一般認為,變異系數(shù)小于等于0.1時為弱變異性,0.1~1為中等變異性,大于1為強變異性[15]。0~60、60~120、120~180 cm土壤含水率的均值分別為0.220、0. 253、0.259 g/g,逐層增加;土壤含鹽量均值分別為1.430、1.036、0.876 g/kg,逐層減小。根系層(0~60 cm)土壤含鹽量分別是深層土壤(60~120、120~180 cm)含鹽量的1.38倍和1.63倍,屬于典型的蒸發(fā)型土壤鹽分剖面[16]。各層土壤含水率的變異系數(shù)分別為0.209、0.142、0.133,均呈中等變異;土壤含鹽量的變異系數(shù)分別為1.217、0.745、0.516,除根系層(0~60 cm)為強變異外,深層土壤(60~120、120~180 cm)均處在中等變異程度。土壤含水率和含鹽量的變異系數(shù)逐層減小,說明土壤水鹽的空間變異性隨土壤深度的增加而減小,主要原因可能是導(dǎo)致土壤含水率和含鹽量空間變異的因素(如氣候、地形和人為活動等)隨著深度的增加其影響減小[7,17]。此外,同一深度土壤含鹽量的變異系數(shù)明顯大于土壤含水率的變異系數(shù),表明土壤含鹽量的空間變異強于土壤含水率[18]。

    表2 土壤含水率和含鹽量的理論擬合模型及特征參數(shù)

    土壤含水率的峰度系數(shù)變化范圍為0.618~1.668,偏度均小于0,數(shù)據(jù)略向左側(cè)拖尾,經(jīng)單樣本K-S正態(tài)檢驗,僅根系層(0~60 cm)土壤含水率不符合正態(tài)分布(=0.004<0.05),深層土壤含水率均呈正態(tài)分布(>0.05);土壤含鹽量的峰度系數(shù)變化范圍為2.710~21.454,偏度變化范圍為1.697~4.245,數(shù)據(jù)左側(cè)拖尾嚴重,經(jīng)單樣本K-S正態(tài)檢驗,各層土壤含鹽量均不符合正態(tài)分布(<0.05),影響半方差函數(shù)的穩(wěn)健性[19]。受人為和自然因素的影響,土壤含鹽量存在少許特異值,利用3原則將其剔除[20],隨后對呈偏態(tài)的數(shù)據(jù)進行Box-Cox變換[21],變換后的數(shù)據(jù)偏度和峰度均接近0,經(jīng)K-S檢驗均符合正態(tài)分布(>0.05),滿足半方差函數(shù)擬合的要求。

    本文利用地統(tǒng)計學軟件GS+9.0,對土壤含水率和含鹽量的半方差函數(shù)進行擬合,擬合結(jié)果如圖2、圖3和表2所示。各層土壤含水率和含鹽量均可采用高斯模型較好地擬合,其RSS在0.001~0.008,2在0.654~0.867,均已達到顯著水平。0~60、60~120、120~180 cm土壤含水率的塊金系數(shù)分別為0.303、0.138、0.023,土壤含鹽量的塊金系數(shù)分別為0.357、0.145、0.105,均逐層減小,土壤含水率的變程分別為700.0、1 282.0、1 433.0 m,土壤含鹽量的變程分別為659.0、990.3、1 333.3 m,均逐層增大。上述結(jié)果表明,土壤含水率和含鹽量的空間相關(guān)性隨土壤深度的增加而增加,空間連續(xù)性和均勻性逐漸增強。根系層(0~60 cm)的土壤含水率和含鹽量具有中等程度的空間相關(guān)性,其空間變異是結(jié)構(gòu)因素(氣候、地形、土質(zhì)等自然因素)和隨機因素(耕作、施肥、灌溉等各種人為活動)共同作用的結(jié)果,而深層土壤(60~120、120~180 cm)含水率和含鹽量空間相關(guān)性較強,其空間變異主要受自然因素的影響。不同深度的土壤含鹽量的塊金系數(shù)均大于土壤含水率,變程均小于土壤含水率,表明土壤含鹽量的空間相關(guān)性及空間均勻性較土壤含水率弱,這是因為鹽隨水動,土壤含鹽量受到土壤含水率的影響,使得土壤含鹽量的隨機變異性增強,空間相關(guān)性減弱。

    圖2 土壤含水率的半方差函數(shù)模型

    圖3 土壤含鹽量的半方差函數(shù)模型

    2.2 土壤含水率和含鹽量的空間分布

    根據(jù)前述最優(yōu)擬合模型及參數(shù),基于普通克里金插值所得結(jié)果如圖4、圖5所示。各層土壤含水率、含鹽量的空間分布具有相似性,總體上,試驗區(qū)東部和東北部的土壤含水率較高,南部的土壤含水率次之,西部和西北部的土壤含水率較低。對于土壤含鹽量而言,高值均出現(xiàn)在研究區(qū)的東部及東北部,低值分布在西南部及西部,說明土壤含水率、含鹽量的垂向分布均較為連續(xù)。此外,土壤含水率與土壤含鹽量高低值分布位置較為接近,二者關(guān)系密切,詳細分析見2.3節(jié)。

    圖4 土壤含水率空間分布

    Fig.4 Spatial distribution of soil water

    圖5 土壤含鹽量空間分布

    統(tǒng)計各土層不同含水率和含鹽量區(qū)間的面積所占比例如圖6所示,隨著深度的增加,土壤含水率逐漸增大,土壤含鹽量逐漸減小。由圖6(a)可知,僅根系層(0~60 cm)土壤含水率出現(xiàn)了0.13~0.20 g/g的低值區(qū),面積占比24%,0~60、60~120 cm和120~180 cm土壤含水率處于0.20~0.25 g/g的面積分別為66%、51%和42%,逐層減小,處于0.25~0.33 g/g的高值區(qū)域面積分別為10%、49%和58%,逐層增大。隆勝試驗區(qū)0~60、60~120 cm和120~180 cm各層土壤含鹽量≤2 g/kg的面積分別為81%、98%和100%,逐層增加,2 g/kg<土壤含鹽量≤3 g/kg的面積分別為17%、2%和0%,逐層減小,僅根系層(0~60 cm)存在3 g/kg<土壤含鹽量≤4 g/kg和土壤含鹽量>4 g/kg的區(qū)域,面積占比分別為1.5%和0.5%。

    圖6 各土層不同土壤含水率和含鹽量區(qū)間面積占比

    2.3 土壤含水率和含鹽量之間的關(guān)系

    土壤水分與鹽分之間存在著復(fù)雜的耦合關(guān)系,由圖4和圖5可知,土壤含水率與含鹽量的空間分布模式存在著較強的相似性,土壤含水率較高的區(qū)域含鹽量一般較高(如圖4、圖5的東部及東北部),土壤含水率相對較低的區(qū)域含鹽量相對較?。ㄈ鐖D4、圖5的西部)。利用ArcGIS軟件的柵格分割功能將土壤含水率和含鹽量的空間插值圖像分割成7個子區(qū),每個子區(qū)不同深度土壤含水率和含鹽量均值如圖7所示。不同深度土壤含水率和含鹽量均呈顯著差異(<0.05),各子區(qū)土壤含水率均值大小為0~60 cm<60~120 cm<120~180 cm,土壤含鹽量均值大小為0~60 cm>60~120 cm>120~180 cm。不同子區(qū)土壤含水率與含鹽量的變化趨勢一致,最大值所在子區(qū)略有差異,土壤含水率和含鹽量最高的區(qū)域編號分別為5和3,分別為研究區(qū)東北部和東部,這可能與研究區(qū)東部的地下水礦化度較高有關(guān)。不同深度土壤含水率與含鹽量的相關(guān)系數(shù)如圖8所示,各層土壤含水率呈顯著強-極強正相關(guān)(0.714<<0.938),距離越近,其相關(guān)性越強。各層土壤含鹽量呈顯著-極強正相關(guān)(0.883<<0.982),土層間距離越近,其相關(guān)性越強。各層土壤含水率與含鹽量呈顯著強相關(guān)關(guān)系(0.615<<0.793),且同一深度土壤含水率與含鹽量的相關(guān)系數(shù)最大,說明土壤含水率與含鹽量密切相關(guān)。本次試驗的取樣時間為生育期初,灌區(qū)尚未開始大規(guī)模灌溉。地下水及其攜帶的鹽分在強烈的蒸發(fā)作用下不斷向上運動,是研究區(qū)土壤含水率和含鹽量的主要補給源。地下水對土壤含水率和含鹽量的補給程度主要受地下水埋深和地下水礦化度的影響。

    圖7 各子區(qū)不同深度土壤含水率和含鹽量均值

    注 SW表示土壤含水率,SS表示土壤含鹽量,P<0.01。

    2.4 地下水埋深及礦化度對土壤含鹽量的影響

    研究區(qū)地下水埋深和礦化度的空間分布如圖9所示。研究區(qū)地下水埋深自東北向西南逐漸增加,研究區(qū)的東部和東北部地下水埋深小于2 m,為強蒸發(fā)積鹽區(qū),占研究區(qū)總面積的6.5%;西部和南部的地下水埋深大于3 m,面積占比6.5%;大多數(shù)區(qū)域的地下水埋深處于2~3 m,面積占比87%。73%的區(qū)域地下水礦化度值小于3.0 g/L,地下水礦化度較高的區(qū)域分布在研究區(qū)的東部,對應(yīng)著土壤含鹽量最高的地區(qū)。根據(jù)上文的分區(qū),統(tǒng)計不同子區(qū)的地下水埋深和礦化度均值,并與各層土壤含鹽量均值進行對比。如圖10所示,各層土壤含鹽量均隨地下水埋深的減小而增加,隨地下水礦化度的增加而增加。各層土壤含鹽量與地下水埋深以及地下水礦化度的相關(guān)系數(shù)矩陣如表3所示,各層土壤含鹽量與地下水埋深呈顯著的強-極強負相關(guān),與地下水礦化度呈極顯著且極強的正相關(guān)。各層土壤含鹽量與地下水埋深及礦化度的關(guān)系如圖11所示,各層土壤含鹽量與地下水埋深均滿足指數(shù)關(guān)系,2>0.504,與地下水礦化度均滿足冪函數(shù)關(guān)系,2>0.875。從圖11也可看出,當?shù)叵滤裆畲笥?.5 m時,總體土壤含鹽量較低,且各層含鹽量變化較小,因此,為有效控制土壤鹽漬化,應(yīng)控制地下水埋深在2.5 m以下。

    圖9 地下水埋深和礦化度空間分布

    圖10 各子區(qū)土壤含鹽量及地下水埋深和礦化度均值

    圖11 土壤含鹽量與地下水埋深及礦化度的關(guān)系

    表3 土壤含鹽量與地下水埋深及礦化度的相關(guān)矩陣

    注*表示<0.05,**表示<0.01。

    3 討論

    本文通過分析隆勝試驗區(qū)土壤含水率和含鹽量空間分布特征,發(fā)現(xiàn)隨著深度的增加,研究區(qū)土壤含水率逐層遞增,土壤含鹽量逐層遞減,屬于典型的表聚型土壤,這是因為研究區(qū)地處干旱-半干旱的河套灌區(qū),干燥少雨,蒸發(fā)強烈,土壤含水率不斷向上運動并在地表蒸發(fā),并將其攜帶的鹽分釋放、累積于表土層,研究結(jié)果與部分學者[5,15,21-22]在干旱地區(qū)的研究結(jié)論一致。除0~60 cm土層的土壤含鹽量呈強變異外,各層土壤含水率和含鹽量均為中等變異。同時,地統(tǒng)計結(jié)果也發(fā)現(xiàn),研究區(qū)0~60 cm土層的土壤含水率和含鹽量表現(xiàn)出中等程度的空間相關(guān)性,變程最小,60~120 cm和120~180 cm土層的土壤含水率和含鹽量呈現(xiàn)出強烈的空間相關(guān)性,變程逐漸增大。經(jīng)典統(tǒng)計學以及地統(tǒng)計學的結(jié)果均表明土壤含水率和含鹽量的空間相關(guān)性隨土壤深度的增加而增加,空間連續(xù)性和均勻性逐漸增強,這可能與隨機因素(耕作、施肥、灌溉等各種人為活動)對深層土壤的作用減弱有關(guān)[7,17]。

    土壤含鹽量一般伴隨著土壤水分的遷移,這使得土壤含鹽量受隨機變異性的影響較土壤水分大,表現(xiàn)為同一深度土壤含水率的變異系數(shù)和塊金系數(shù)小于土壤含鹽量,即土壤含鹽量的空間變異性強于土壤含水率,空間相關(guān)性弱于土壤含水率,這與楊勁松等[10]的結(jié)果一致。另一方面,土壤含水率和含鹽量存在顯著的正相關(guān)關(guān)系,即土壤含水率高的區(qū)域土壤含鹽量往往較大,土壤含水率低的區(qū)域土壤含鹽量往往較小,這是因為較高的土壤含水率表示地下水對土壤的補給量較多,在地下水礦化度相同的情況下,其攜帶的含鹽量越多,反之亦然。

    地下水埋深是影響地下水對土壤水補給量的重要因素,而地下水礦化度影響著隨水分運移的鹽量,因此地下水埋深和礦化度均與土壤含鹽量密切相關(guān)。本研究表明,土壤含鹽量與地下水埋深滿足指數(shù)關(guān)系,這與竇旭等[18]和常春龍等[23]等在河套灌區(qū)的結(jié)論一致,土壤含鹽量與地下水礦化度滿足冪函數(shù)關(guān)系,而王金哲等[24]發(fā)現(xiàn)環(huán)渤海平原區(qū)土壤含鹽量與地下水礦化度之間呈指數(shù)函數(shù)關(guān)系,這可能與研究區(qū)域不同有關(guān)。本研究區(qū)為井渠結(jié)合區(qū),采用地下水灌溉顯著改變了土壤-地下水含水層的補排關(guān)系,從而導(dǎo)致土壤含鹽量與地下水礦化度之間的關(guān)系發(fā)生改變。此外,當?shù)叵滤裆畲笥?.5 m時,土壤含鹽量較低,且隨著土壤深度的增大,土壤含鹽量的變化趨勢減少。因此,為有效控制土壤鹽漬化,本文建議地下水埋深控制在2.5 m以下,這與徐英等[25]和楊會峰等[26]在河套灌區(qū)的結(jié)論一致。土壤含鹽量隨地下水埋深的增加而減小,隨地下水礦化度的增加而增加,由于研究區(qū)地下水埋深在年內(nèi)的變化規(guī)律十分明顯,而地下水礦化度年內(nèi)變化規(guī)律不明顯[11],因此,合理調(diào)控地下水埋深將是灌區(qū)鹽堿化防控的關(guān)鍵。

    4 結(jié)論

    1)隨著土壤深度的增加,土壤含水率逐漸增加,其空間變異性逐漸減小,土壤含鹽量及其變異性均逐漸減小。研究區(qū)根系層土壤含水率和含鹽量的空間變化受結(jié)構(gòu)性和隨機因素的共同作用,深層土壤的含水率和含鹽量空間變異主要受結(jié)構(gòu)性因素的作用。

    2)土壤含水率與含鹽量密切相關(guān),一方面,土壤含水率的空間變異性弱于土壤含鹽量,另一方面,土壤含水率較低的區(qū)域土壤含鹽量相對較大,反之亦然,二者存在著顯著的正相關(guān)關(guān)系。

    3)研究區(qū)各層土壤含鹽量均與地下水埋深呈指數(shù)相關(guān)關(guān)系;與地下水礦化度呈冪函數(shù)關(guān)系。當?shù)叵滤裆畲笥?.5m時,土壤含鹽量低,且各層土壤含鹽量變化較小。因此,建議地下水埋深控制在2.5m以下。

    [1] LOESCGER W, CHAN Z, GRUMET R. Options for developing salt-tolerant crops[J]. Hort Science, 2011, 46(8): 1 085-1 092.

    [2] WANG Yugang, XIAO Duning, LI Yan, et al. Soil salinity evolution and its relationship with dynamics of groundwater in the oasis of inland river basins: case study from the Fubei region of Xinjiang Province, China[J]. Environmental Monitoring and Assessment, 2008, 140(1): 291-302.

    [3] WANG Yugang, DENG Caiyun, LIU Yan, et al. Identifying change in spatial accumulation of soil salinity in an inland river watershed, China[J]. Science of the Total Environment, 2018, 621: 177-185.

    [4] XIE Wenping, YANG Jingsong, YAO Rongjiang, et al. Spatial and temporal variability of soil salinity in the Yangtze River Estuary using electromagnetic induction[J]. Remote Sensing, 2021, 13(10): 1 875.

    [5] 劉洪波, 丁邦新, 白云崗, 等. 典型干旱區(qū)綠洲春季土壤鹽分空間分布特征分析[J]. 土壤通報, 2021, 52(2): 279-285.

    LIU Hongbo, DING Bangxin, BAI Yungang, et al. Spatial distribution characteristics of soil salinity in an oasis of typical arid regions in spring[J]. Chinese Journal of Soil Science, 2021, 52(2): 279-285.

    [6] 婁婷, 馮紹元, 蘇濤, 等. 河套灌區(qū)土壤鹽堿化與地下水埋深關(guān)系分析[J]. 灌溉排水學報, 2014, 33(S1): 132-134.

    LOU Ting, FENG Shaoyuan, SU Tao, et al. Relationship between soil salinization and groundwater depth in Hetao Irrigation District[J]. Journal of Irrigation and Drainage, 2014, 33(S1): 132-134.

    [7] 張雪, 陳耀軒, 張盼, 等. 孔雀河流域農(nóng)業(yè)灌區(qū)土壤鹽分空間分布及變異特征[J]. 土壤, 2016, 48(6): 1 223-1 228.

    ZHANG Xue, CHEN Yaoxuan, ZHANG Pan, et al. Soil salinity spatial distribution and variation characteristics in irrigated land of Kongque River Basin[J]. Soils, 2016, 48(6): 1 223-1 228.

    [8] 常曉敏, 王少麗, 陳皓銳, 等. 河套灌區(qū)土壤鹽分時空變化特征及影響因素[J]. 排灌機械工程學報, 2018, 36(10): 1 000-1 005.

    CHANG Xiaomin, WANG Shaoli, CHEN Haorui, et al. Spatiotemporal changes and influencing factors of soil salinity in Hetao Irrigation District[J]. Journal of Drainage and Irrigation Machinery Engineering, 2018, 36(10): 1 000-1 005.

    [9] 毛威, 楊金忠, 朱焱, 等. 河套灌區(qū)井渠結(jié)合膜下滴灌土壤鹽分演化規(guī)律[J]. 農(nóng)業(yè)工程學報, 2018, 34(1): 93-101.

    MAO Wei, YANG Jinzhong, ZHU Yan, et al. Soil salinity process of Hetao Irrigation District after application of well-canal conjunctive irrigation and mulched drip irrigation[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(1): 93-101.

    [10] 楊勁松, 姚榮江. 黃河三角洲地區(qū)土壤水鹽空間變異特征研究[J]. 地理科學, 2007, 27(3): 348-353.

    YANG Jingsong, YAO Rongjiang. Spatial variability of soil water and salt characteristics in the Yellow River Delta[J]. Scientia Geographica Sinica, 2007, 27(3): 348-353.

    [11] 姚玲, 楊洋, 孫貫芳, 等. 基于地統(tǒng)計分析的河套灌區(qū)地下水埋深與礦化度時空變異規(guī)律研究[J]. 灌溉排水學報, 2020, 39(8): 111-121.

    YAO Ling, YANG Yang, SUN Guanfang, et al. Spatiotemporal variation in depth of groundwater table and water salinity in Hetao Irrigation District estimated using geostatistics[J]. Journal of Irrigation and Drainage, 2020, 39(8): 111-121.

    [12] 杜軍, 楊培嶺, 李云開, 等. 河套灌區(qū)年內(nèi)地下水埋深與礦化度的時空變化[J]. 農(nóng)業(yè)工程學報, 2010, 26(7): 26-31.

    DU Jun, YANG Peiling, LI Yunkai, et al. Analysis of spatial and temporal variations of groundwater level and its salinity in Hetao Irrigation District[J]. Transactions of the Chinese Society of Agricultural Engineering, 2010, 26(7): 26-31.

    [13] 李亮亮, 依艷麗, 凌國鑫, 等. 地統(tǒng)計學在土壤空間變異研究中的應(yīng)用[J]. 土壤通報, 2005(2): 265-268.

    LI Liangliang, YI Yanli, LING Guoxin, et al. Utilization of geostatistics in soil spatial variability[J]. Chinese Journal of Soil Science, 2005(2): 265-268.

    [14] HU Xuhua, CHEN Mengting, LIU Dong, et al. Reference evapotranspiration change in Heilongjiang Province, China from 1951 to 2018: the role of climate change and rice area expansion[J]. Agricultural Water Management, 2021, 253: 106 912.

    [15] 化騫寂, 馮紹元, 葛洲, 等. 河套灌區(qū)典型區(qū)周年內(nèi)耕層土壤鹽分時空變異研究[J]. 灌溉排水學報, 2020, 39(8): 26-34.

    HUA Qianji, FENG Shaoyuan, GE Zhou, et al. Intra-annual spatiotemporal variation in salt content in the plough layer in Hetao Irrigation District[J]. Journal of Irrigation and Drainage, 2020, 39(8): 26-34.

    [16] 孫貫芳. 干旱區(qū)多尺度土壤鹽分動態(tài)及長期均衡模型研究[D]. 武漢:武漢大學, 2020.

    SUN Guanfang. Study on soil salinity dynamics at multiple scales and long-term solute balance models in arid area[D]. Wuhan: Wuhan University, 2020.

    [17] 竇旭. 河套灌區(qū)暗管排水排鹽有效性評價與土壤水肥鹽時空變異規(guī)律研究[D]. 呼和浩特: 內(nèi)蒙古農(nóng)業(yè)大學, 2020

    DOU Xu, Study on the effectiveness evaluation of subsurface drainage and salt and spatiotemporal variation of soil water, fertilizer and salt in Hetao Irrigation District[D]. Hohhot: Inner Mongolia Agricultural University, 2020.

    [18] 竇旭, 史海濱, 苗慶豐, 等. 鹽漬化灌區(qū)土壤水鹽時空變異特征分析及地下水埋深對鹽分的影響[J]. 水土保持學報, 2019, 33(3): 246-253.

    DOU Xu, SHI Haibin, MIAO Qingfeng, et al. Temporal and spatial variability analysis of soil water and salt and the influence of groundwater depth on salt in saline irrigation area[J]. Journal of Soil and Water Conservation, 2019, 33(3): 246-253.

    [19] 楊奇勇, 楊勁松, 李曉明. 不同閾值下土壤鹽分的空間變異特征研究[J]. 土壤學報, 2011, 48(6): 1 109-1 115.

    YANG Qiyong, YANG Jingsong, LI Xiaoming. Spatial variability of soil salinity research under different thesholds[J]. Acta Pedologica Sinica, 2011, 48(6): 1 109-1 115.

    [20] 郭安廷, 崔錦霞, 許鑫, 等. 基于GIS 與地統(tǒng)計的土壤養(yǎng)分空間變異研究[J]. 中國農(nóng)學通報, 2018, 34(23): 72-79.

    GUO Anting, CUI Jinxia, XU Xin, et al. Spatial distribution of soil nutrients based on GIS and geostatistics[J]. Chinese Agricultural Science Bulletin, 2018, 34(23): 72-79.

    [21] 王瑞萍, 萬核洋, 陳帥, 等. 河套灌區(qū)烏拉特灌域春季土壤鹽堿化空間分布特征[J]. 農(nóng)業(yè)工程學報, 2021, 37(15): 105-113.

    WANG Ruiping, WAN Heyang, CHEN Shuai, et al. Spatial distribution pattern of soil salinization in farmland of Wulate sub-irrigation areas in Hetao Irrigation District in Inner Mongolia in spring[J]. Transactions of the Chinese Society of Agricultural Engineering, 2021, 37(15): 105-113.

    [22] 葛洲. 永濟灌域典型區(qū)土壤鹽堿化時空分布控制因素研究[D]. 揚州:揚州大學, 2019.

    GE Zhou. Study on controlling factors of spatial-temporal distribution of soil salinization in typical area of Yongji Irrigation area[D]. Yangzhou: Yangzhou University, 2019.

    [23] 常春龍, 楊樹青, 劉德平, 等. 河套灌區(qū)上游地下水埋深與土壤鹽分互作效應(yīng)研究[J]. 灌溉排水學報, 2014, 33(4/5): 315-319.

    CHANG Chunlong, YANG Shuqing, LIU Deping, et al. Interaction effect of groundwater depth and soil salinization in the upper of Hetao Irrigation Area[J]. Journal of Irrigation and Drainage, 2014, 33(4/5): 315-319.

    [24] 王金哲, 張光輝, 嚴明疆, 等. 環(huán)渤海平原區(qū)土壤鹽分分布特征及影響因素分析[J]. 干旱區(qū)資源與環(huán)境, 2012, 26(11): 104-109.

    WANG Jinzhe, ZHANG Guanghui, YAN Mingjiang, et al. Analysis of soil salinity distribution and influencing factors in area around Bohai sea[J]. Journal of Arid Land Resources and Environment, 2012, 26(11): 104-109.

    [25] 徐英, 葛洲, 王娟, 等. 基于指示 Kriging 法的土壤鹽漬化與地下水埋深關(guān)系研究[J]. 農(nóng)業(yè)工程學報, 2019, 35(1): 123-130.

    XU Ying, GE Zhou, WANG Juan, et al.Study on relationship between soil salinization and groundwater table depth based on indicator Kriging[J]. Transactions of the Chinese Society of Agricultural Engineering, 2019, 35(1): 123-130.

    [26] 楊會峰, 張發(fā)旺, 王貴玲, 等. 河套平原次生鹽漬化地區(qū)地下水動態(tài)調(diào)控模擬研究[J]. 南水北調(diào)與水利科技, 2011, 9(3): 63-67.

    YANG Huifeng, ZHANG Fawang, WANG Guiling, et al. Simulation study on the dynamic regulation of groundwater in secondary salinized areas of Hetao Plain[J]. South-to-North Water Diversion and Water Science and Technology, 2011, 9(3): 63-67.

    Geostatistical Analysis of Spatial Variability of Soil Water and Salt in Hetao Irrigation District

    LIU Yan’nan1, SUN Guanfang2, MAO Wei1, CHENG Xiaoyao3, ZHU Yan1, YANG Jinzhong1*

    (1.The State Key Laboratory of Water Resource and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China;2.Institute of Soil and Water Conservation, Northwest Agriculture and Forestry University, Yangling 712100, China; 3. Zhongnan Municipal Administration General Design Institute Co. Ltd., China Municipal Engineering, Wuhan 430101, China)

    【Background and objective】Soil salinity is an abiotic stress affecting agricultural production worldwide, especially in countries in arid and semi-arid regions. Soil salinity is closely correlated to soil water, both varying erratically over space and time. Understanding their spatial variability is essential to improving water use efficiency and alleviating soil salinization. This paper aims to analyze the spatial variability of soil water and salt at different soil depths, their relationship in Hetao Irrigation District, as well as their dependence on groundwater depth and salinity.【Method】The study site is at Longsheng, a typical well-canal conjunctive irrigation area in the irrigation district. Soils were sampled from the root zone (0~60 cm), subsoil (60~120 cm) and deep soil (120~180 cm), and the moisture and salt contents in each sample were measured using standard laboratory methods. Their spatial variations were analyzed using classical statistical and geostatistical methods. The relationship between groundwater and soil salinity was calculated using the Pearson correlation analysis.【Result】① The averaged soil water content at different soil depths were ranked in the order of 0~60 cm<60~120 cm<120~180 cm, while the averaged soil salt content were ranked in the order of 0~60 cm>60~120 cm>120~180 cm. The variation coefficient and nugget coefficient of soil water and salt content decreased monotonically with soil depth, while their spatial autocorrelation distance increases monotonically with the depth. Soil water and salt content in the root zone showed moderate spatial autocorrelation, while in deep soil they showed a strong spatial autocorrelation.② The variation coefficient and nugget coefficient of soil water in each soil layer are smaller than that of soil salt, while its variation range is larger than that of soil salt; there was a strong negative correlation between them.③ Soil salt in each soil layer is exponentially related to groundwater depth, but is a power-law function of groundwater salinity. When the groundwater depth is greater than 2.5 m, the soil salt content is low and the soil salinity only slightly changes with groundwater depth.【Conclusion】Salt accumulates mainly in the root zone and its spatial distribution is closely related to soil water content. Groundwater depth is the major controlling factor of soil salinity, and to reduce the risk of soil salinity, the groundwater depth should be controlled below the depth of 2.5 m.

    Hetao Irrigation District; soil water and salt; geostatistics; spatial variation; groundwater environment

    1672 - 3317(2022)09 - 0101 - 09

    P641.2;S27

    A

    10.13522/j.cnki.ggps.2021621

    劉燕楠, 孫貫芳, 毛威, 等. 基于地統(tǒng)計分析的河套灌區(qū)土壤水鹽空間變異規(guī)律分析[J]. 灌溉排水學報, 2022, 41(9): 101-109.

    LIU Yan’nan, SUN Guanfang, MAO Wei, et al. Geostatistical Analysis of Spatial Variability of Soil Water and Salt in Hetao Irrigation District[J]. Journal of Irrigation and Drainage, 2022, 41(9): 101-109.

    2021-12-13

    國家自然科學基金項目(52179041)

    劉燕楠(1998-),女。碩士研究生,主要從事土壤水地下水資源與環(huán)境方面的研究。E-mail: lyn214622@163.com

    楊金忠(1953-),男。教授,主要從事飽和-非飽和水流運動與溶質(zhì)運移方面的研究。E-mail: jzyang@whu.edu.cn

    責任編輯:趙宇龍

    猜你喜歡
    含鹽量礦化度變異
    含鹽量及含水率對鹽漬土凍脹規(guī)律影響試驗研究*
    黃河三角洲鹽漬土有機氮組成及氮有效性對土壤含鹽量的響應(yīng)*
    變異危機
    變異
    什么是水的含鹽量?
    離子組成及礦化度對低礦化度水驅(qū)采收率的影響
    巖性油氣藏(2019年5期)2019-09-11 03:54:44
    秦陵陪葬坑土遺址安全含鹽量探究
    變異的蚊子
    百科知識(2015年18期)2015-09-10 07:22:44
    利用綜合測井資料計算地下水等效NaCl溶液礦化度方法
    低礦化度水驅(qū)技術(shù)增產(chǎn)機理與適用條件
    国产一区二区在线观看av| 国产精品国产三级国产专区5o| 国产极品粉嫩免费观看在线| 亚洲欧美精品自产自拍| 视频区欧美日本亚洲| 午夜免费观看性视频| 国产亚洲精品第一综合不卡| 国产亚洲av片在线观看秒播厂| 久久人人97超碰香蕉20202| 一二三四在线观看免费中文在| 女性被躁到高潮视频| 国产老妇伦熟女老妇高清| 国产精品免费视频内射| 欧美亚洲日本最大视频资源| 欧美少妇被猛烈插入视频| 久久精品久久久久久噜噜老黄| 午夜老司机福利片| 99香蕉大伊视频| 国产欧美日韩一区二区精品| 亚洲国产av新网站| 国产激情久久老熟女| xxxhd国产人妻xxx| 一本大道久久a久久精品| 黄网站色视频无遮挡免费观看| 老汉色∧v一级毛片| 成人av一区二区三区在线看 | 在线永久观看黄色视频| 亚洲欧美色中文字幕在线| 97在线人人人人妻| 国产又爽黄色视频| 亚洲,欧美精品.| 国产91精品成人一区二区三区 | 亚洲国产精品成人久久小说| 在线 av 中文字幕| a级毛片黄视频| 国产av国产精品国产| 亚洲国产欧美网| 久久精品亚洲熟妇少妇任你| 日韩免费高清中文字幕av| 天天躁日日躁夜夜躁夜夜| 欧美 亚洲 国产 日韩一| 在线 av 中文字幕| 欧美日韩精品网址| 99国产精品一区二区三区| 狂野欧美激情性xxxx| 少妇的丰满在线观看| 国产精品秋霞免费鲁丝片| 视频区图区小说| 精品久久久久久电影网| 日日夜夜操网爽| 悠悠久久av| 99国产精品免费福利视频| 久久九九热精品免费| 成人18禁高潮啪啪吃奶动态图| 国产免费视频播放在线视频| 亚洲国产av影院在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲全国av大片| www.自偷自拍.com| 美女中出高潮动态图| 精品乱码久久久久久99久播| 午夜福利,免费看| 国产精品九九99| 黑丝袜美女国产一区| 久久精品国产亚洲av高清一级| 亚洲黑人精品在线| 这个男人来自地球电影免费观看| 精品久久久久久久毛片微露脸 | 一区二区av电影网| 一区在线观看完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av男天堂| av网站在线播放免费| 国产一区二区三区综合在线观看| 日韩中文字幕欧美一区二区| 下体分泌物呈黄色| 每晚都被弄得嗷嗷叫到高潮| 97精品久久久久久久久久精品| 亚洲国产毛片av蜜桃av| kizo精华| 麻豆av在线久日| 少妇的丰满在线观看| 精品国产国语对白av| 亚洲专区中文字幕在线| 日本精品一区二区三区蜜桃| 老鸭窝网址在线观看| 日韩欧美一区二区三区在线观看 | 啪啪无遮挡十八禁网站| 狂野欧美激情性xxxx| 日本av手机在线免费观看| 91精品伊人久久大香线蕉| 蜜桃国产av成人99| 精品国产一区二区三区四区第35| 午夜福利在线免费观看网站| 国产一区二区三区av在线| 熟女少妇亚洲综合色aaa.| 国产1区2区3区精品| 亚洲国产欧美网| 在线十欧美十亚洲十日本专区| 国产视频一区二区在线看| 蜜桃国产av成人99| 12—13女人毛片做爰片一| 热re99久久精品国产66热6| 久久久久久免费高清国产稀缺| 欧美 日韩 精品 国产| 亚洲精品美女久久av网站| 女警被强在线播放| 国产老妇伦熟女老妇高清| 少妇被粗大的猛进出69影院| 99精品久久久久人妻精品| 男男h啪啪无遮挡| 极品少妇高潮喷水抽搐| 国产精品久久久久久人妻精品电影 | 欧美xxⅹ黑人| 国产日韩一区二区三区精品不卡| 91字幕亚洲| 久久中文字幕一级| 精品第一国产精品| 国产野战对白在线观看| a在线观看视频网站| 亚洲欧美日韩高清在线视频 | 90打野战视频偷拍视频| 国产亚洲欧美精品永久| 999精品在线视频| 一个人免费在线观看的高清视频 | 黑人巨大精品欧美一区二区蜜桃| 精品免费久久久久久久清纯 | 日韩大片免费观看网站| 欧美日韩精品网址| 国产三级黄色录像| 18禁裸乳无遮挡动漫免费视频| 曰老女人黄片| 黄片小视频在线播放| 别揉我奶头~嗯~啊~动态视频 | 久久天躁狠狠躁夜夜2o2o| 久久午夜综合久久蜜桃| 极品少妇高潮喷水抽搐| 亚洲五月婷婷丁香| 一本大道久久a久久精品| 久久久精品免费免费高清| 日本av免费视频播放| 午夜日韩欧美国产| 国产主播在线观看一区二区| 视频区图区小说| 我要看黄色一级片免费的| av视频免费观看在线观看| 亚洲天堂av无毛| 久久精品国产综合久久久| 日韩电影二区| 日本一区二区免费在线视频| 下体分泌物呈黄色| 久久女婷五月综合色啪小说| 国产高清videossex| 老司机靠b影院| 精品国产一区二区三区久久久樱花| 色综合欧美亚洲国产小说| 亚洲五月色婷婷综合| 成年人免费黄色播放视频| 99久久精品国产亚洲精品| 老熟妇乱子伦视频在线观看 | 国产成人啪精品午夜网站| 黄网站色视频无遮挡免费观看| 性高湖久久久久久久久免费观看| 婷婷丁香在线五月| 欧美 日韩 精品 国产| 少妇裸体淫交视频免费看高清 | 亚洲欧美色中文字幕在线| 丝袜人妻中文字幕| 欧美人与性动交α欧美软件| 1024香蕉在线观看| 在线av久久热| 欧美 日韩 精品 国产| 国产伦人伦偷精品视频| 精品国内亚洲2022精品成人 | 精品久久久精品久久久| 国产精品二区激情视频| 久久久国产一区二区| 亚洲精品中文字幕在线视频| 天堂中文最新版在线下载| av天堂久久9| 好男人电影高清在线观看| 69av精品久久久久久 | 中文字幕制服av| 性色av乱码一区二区三区2| 91大片在线观看| 亚洲伊人色综图| 中亚洲国语对白在线视频| 丁香六月欧美| 久久亚洲国产成人精品v| 日本vs欧美在线观看视频| 下体分泌物呈黄色| 99国产综合亚洲精品| av天堂久久9| 97在线人人人人妻| 国产熟女午夜一区二区三区| 水蜜桃什么品种好| 天天操日日干夜夜撸| 久久久国产一区二区| 免费av中文字幕在线| 男人添女人高潮全过程视频| 国产精品香港三级国产av潘金莲| 免费少妇av软件| 五月开心婷婷网| av超薄肉色丝袜交足视频| 高清av免费在线| 久久精品亚洲熟妇少妇任你| 精品一区二区三卡| 午夜福利免费观看在线| 精品熟女少妇八av免费久了| 啦啦啦在线免费观看视频4| 欧美日韩亚洲高清精品| 美女高潮喷水抽搐中文字幕| 亚洲一区二区三区欧美精品| 久久久国产一区二区| 91九色精品人成在线观看| 成人av一区二区三区在线看 | 国产黄色免费在线视频| 国产精品免费大片| 99久久99久久久精品蜜桃| 亚洲专区中文字幕在线| 亚洲精品国产精品久久久不卡| 桃花免费在线播放| 欧美日韩亚洲高清精品| a级毛片在线看网站| 国产成人精品久久二区二区免费| 精品少妇黑人巨大在线播放| 久久久精品94久久精品| 国产日韩欧美在线精品| 一区二区av电影网| 蜜桃国产av成人99| 汤姆久久久久久久影院中文字幕| a级毛片黄视频| 成人免费观看视频高清| 精品国产乱子伦一区二区三区 | 国产一区二区 视频在线| 丝袜喷水一区| 人人妻人人澡人人看| 久热爱精品视频在线9| 免费高清在线观看日韩| 人人澡人人妻人| 国产高清视频在线播放一区 | 欧美日韩一级在线毛片| 窝窝影院91人妻| 91九色精品人成在线观看| 亚洲全国av大片| 叶爱在线成人免费视频播放| 中文字幕色久视频| 国产成+人综合+亚洲专区| tube8黄色片| 国产免费福利视频在线观看| 黑人猛操日本美女一级片| 久久国产精品影院| 午夜福利,免费看| 视频区图区小说| 夜夜骑夜夜射夜夜干| 丝袜喷水一区| 性色av乱码一区二区三区2| 免费不卡黄色视频| 高清欧美精品videossex| av在线老鸭窝| 91av网站免费观看| 高清av免费在线| av电影中文网址| 蜜桃国产av成人99| 91精品三级在线观看| 香蕉国产在线看| 男女国产视频网站| 亚洲自偷自拍图片 自拍| 欧美+亚洲+日韩+国产| 久久中文看片网| 肉色欧美久久久久久久蜜桃| 日韩有码中文字幕| 多毛熟女@视频| 在线 av 中文字幕| 97精品久久久久久久久久精品| 精品乱码久久久久久99久播| 亚洲激情五月婷婷啪啪| 久久久久久久久免费视频了| 欧美激情久久久久久爽电影 | 美女大奶头黄色视频| 极品人妻少妇av视频| 亚洲国产精品成人久久小说| 精品人妻在线不人妻| 十八禁人妻一区二区| 亚洲情色 制服丝袜| 丝袜在线中文字幕| 国产在线免费精品| 亚洲国产av新网站| av电影中文网址| 在线观看一区二区三区激情| 欧美日韩亚洲综合一区二区三区_| 久久久精品国产亚洲av高清涩受| 啦啦啦视频在线资源免费观看| 午夜视频精品福利| 久久精品熟女亚洲av麻豆精品| 狠狠狠狠99中文字幕| 青草久久国产| 中国国产av一级| 一二三四在线观看免费中文在| 精品一品国产午夜福利视频| 国产精品av久久久久免费| 一区二区av电影网| 亚洲av电影在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 久久精品国产亚洲av香蕉五月 | 久久精品国产a三级三级三级| 国内毛片毛片毛片毛片毛片| 国产一区二区在线观看av| 国产成人a∨麻豆精品| 国产精品99久久99久久久不卡| 嫩草影视91久久| 久久国产精品大桥未久av| 中文精品一卡2卡3卡4更新| 久久人人爽人人片av| 亚洲国产日韩一区二区| 亚洲伊人久久精品综合| 在线十欧美十亚洲十日本专区| 欧美黑人精品巨大| 中文欧美无线码| 亚洲中文日韩欧美视频| 国产一区有黄有色的免费视频| 伊人亚洲综合成人网| 两个人看的免费小视频| 各种免费的搞黄视频| 国产成人免费无遮挡视频| 丝袜喷水一区| 亚洲国产欧美一区二区综合| 精品第一国产精品| 12—13女人毛片做爰片一| 亚洲成av片中文字幕在线观看| 如日韩欧美国产精品一区二区三区| 爱豆传媒免费全集在线观看| 欧美精品啪啪一区二区三区 | 女性生殖器流出的白浆| 国产熟女午夜一区二区三区| 少妇猛男粗大的猛烈进出视频| 久久久久久久国产电影| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 国精品久久久久久国模美| 在线观看www视频免费| 精品人妻熟女毛片av久久网站| 亚洲国产av新网站| 国产亚洲一区二区精品| 十八禁人妻一区二区| 侵犯人妻中文字幕一二三四区| 久久久久久久久久久久大奶| 少妇 在线观看| 在线看a的网站| 亚洲 欧美一区二区三区| 日本五十路高清| 中亚洲国语对白在线视频| 各种免费的搞黄视频| 秋霞在线观看毛片| 国产免费现黄频在线看| 精品一区二区三区四区五区乱码| 在线十欧美十亚洲十日本专区| 人妻一区二区av| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美日韩另类电影网站| 午夜福利在线免费观看网站| 高清av免费在线| 日韩中文字幕欧美一区二区| 久久女婷五月综合色啪小说| 中文精品一卡2卡3卡4更新| 建设人人有责人人尽责人人享有的| 一个人免费在线观看的高清视频 | 亚洲中文日韩欧美视频| 精品亚洲乱码少妇综合久久| 久久免费观看电影| 真人做人爱边吃奶动态| 伊人久久大香线蕉亚洲五| 亚洲九九香蕉| 国产亚洲午夜精品一区二区久久| 久久久久久免费高清国产稀缺| 亚洲性夜色夜夜综合| 9191精品国产免费久久| 俄罗斯特黄特色一大片| 五月天丁香电影| 大型av网站在线播放| 老司机影院成人| 国产成人精品久久二区二区91| 欧美激情久久久久久爽电影 | 国产精品亚洲av一区麻豆| 一个人免费看片子| 日本猛色少妇xxxxx猛交久久| 一区福利在线观看| 制服人妻中文乱码| 日韩,欧美,国产一区二区三区| 国精品久久久久久国模美| 亚洲中文av在线| 大片免费播放器 马上看| 黑丝袜美女国产一区| 91精品伊人久久大香线蕉| 久久av网站| 中文精品一卡2卡3卡4更新| 黑人巨大精品欧美一区二区mp4| 日韩电影二区| 免费看十八禁软件| 亚洲色图综合在线观看| 一区二区日韩欧美中文字幕| 老司机午夜十八禁免费视频| 亚洲欧美日韩另类电影网站| 中文字幕人妻丝袜一区二区| 国产亚洲午夜精品一区二区久久| 亚洲精品一区蜜桃| 国产黄色免费在线视频| 国产免费一区二区三区四区乱码| 18禁国产床啪视频网站| 精品久久蜜臀av无| www.自偷自拍.com| 久久九九热精品免费| 亚洲国产精品成人久久小说| 久久久国产精品麻豆| 伦理电影免费视频| 天天影视国产精品| 精品国产乱码久久久久久小说| 久久人人爽人人片av| 成在线人永久免费视频| 午夜激情久久久久久久| 热re99久久精品国产66热6| 免费观看人在逋| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品国产av蜜桃| av网站免费在线观看视频| 久久国产亚洲av麻豆专区| 午夜福利,免费看| 99热网站在线观看| 岛国在线观看网站| 亚洲第一青青草原| 精品久久久久久电影网| 日韩 亚洲 欧美在线| 亚洲av日韩精品久久久久久密| 日本猛色少妇xxxxx猛交久久| 丝袜在线中文字幕| 亚洲av日韩在线播放| 人人妻人人添人人爽欧美一区卜| av线在线观看网站| 一级,二级,三级黄色视频| 少妇猛男粗大的猛烈进出视频| 美女午夜性视频免费| 男人操女人黄网站| 王馨瑶露胸无遮挡在线观看| 80岁老熟妇乱子伦牲交| bbb黄色大片| 男人爽女人下面视频在线观看| 美女高潮喷水抽搐中文字幕| 国产区一区二久久| 国产视频一区二区在线看| 国产精品自产拍在线观看55亚洲 | 在线观看人妻少妇| 青草久久国产| 超碰成人久久| 精品人妻熟女毛片av久久网站| 成年女人毛片免费观看观看9 | 欧美精品一区二区免费开放| 91精品伊人久久大香线蕉| 一区二区日韩欧美中文字幕| 美女主播在线视频| 手机成人av网站| 亚洲精品粉嫩美女一区| 日本猛色少妇xxxxx猛交久久| 狠狠精品人妻久久久久久综合| 大型av网站在线播放| 中文字幕av电影在线播放| 久久狼人影院| 日本一区二区免费在线视频| 国产精品免费视频内射| 大陆偷拍与自拍| 久久久久久久久免费视频了| 少妇精品久久久久久久| 亚洲天堂av无毛| 亚洲午夜精品一区,二区,三区| 国产一区二区激情短视频 | av在线老鸭窝| 欧美成狂野欧美在线观看| 欧美人与性动交α欧美软件| 91精品三级在线观看| 熟女少妇亚洲综合色aaa.| 天天躁夜夜躁狠狠躁躁| 免费看十八禁软件| 无限看片的www在线观看| 亚洲成av片中文字幕在线观看| 大码成人一级视频| 午夜视频精品福利| 丝袜人妻中文字幕| 五月天丁香电影| 午夜福利免费观看在线| 天堂8中文在线网| 久久亚洲国产成人精品v| 亚洲精品国产色婷婷电影| 午夜免费观看性视频| 国产精品免费大片| av欧美777| 国产激情久久老熟女| 男女无遮挡免费网站观看| 精品欧美一区二区三区在线| 无遮挡黄片免费观看| 女人精品久久久久毛片| 日韩欧美国产一区二区入口| 99香蕉大伊视频| 一区二区三区激情视频| cao死你这个sao货| 国产成+人综合+亚洲专区| 国产三级黄色录像| 久久久久久久久免费视频了| av视频免费观看在线观看| 色综合欧美亚洲国产小说| av在线app专区| 99热网站在线观看| 久久久久久久国产电影| 这个男人来自地球电影免费观看| 午夜福利影视在线免费观看| 欧美日本中文国产一区发布| 久久久精品94久久精品| 日韩视频在线欧美| 色婷婷av一区二区三区视频| 老汉色∧v一级毛片| 又大又爽又粗| 在线 av 中文字幕| 美女脱内裤让男人舔精品视频| 窝窝影院91人妻| 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| 在线观看一区二区三区激情| 国产一区二区在线观看av| 亚洲精品中文字幕在线视频| 亚洲精华国产精华精| 亚洲一区中文字幕在线| 亚洲人成77777在线视频| 亚洲色图综合在线观看| 精品乱码久久久久久99久播| √禁漫天堂资源中文www| 欧美av亚洲av综合av国产av| 菩萨蛮人人尽说江南好唐韦庄| 悠悠久久av| 国产一区二区三区av在线| 黑人巨大精品欧美一区二区mp4| 大型av网站在线播放| 国产在视频线精品| 久久国产精品大桥未久av| 91精品国产国语对白视频| 亚洲国产日韩一区二区| 高潮久久久久久久久久久不卡| 久久亚洲国产成人精品v| 亚洲欧洲日产国产| 在线观看免费视频网站a站| 亚洲综合色网址| 日本欧美视频一区| 日韩中文字幕视频在线看片| 国产野战对白在线观看| 亚洲国产av影院在线观看| 黑人巨大精品欧美一区二区蜜桃| 99热国产这里只有精品6| 老司机靠b影院| 亚洲色图 男人天堂 中文字幕| 久久久精品区二区三区| 国产成人影院久久av| 亚洲专区国产一区二区| av网站免费在线观看视频| 女性被躁到高潮视频| 亚洲精品中文字幕一二三四区 | 精品少妇内射三级| 亚洲精品美女久久av网站| 精品国内亚洲2022精品成人 | 热99re8久久精品国产| 狠狠精品人妻久久久久久综合| 最黄视频免费看| 成年av动漫网址| 人人妻人人澡人人看| 成人国产一区最新在线观看| 午夜老司机福利片| 国产深夜福利视频在线观看| 三级毛片av免费| tube8黄色片| 少妇的丰满在线观看| 免费少妇av软件| 91麻豆av在线| 午夜福利乱码中文字幕| 亚洲avbb在线观看| 免费少妇av软件| 大片电影免费在线观看免费| 波多野结衣av一区二区av| 少妇的丰满在线观看| 国产在线观看jvid| 中文精品一卡2卡3卡4更新| 国产淫语在线视频| 欧美国产精品一级二级三级| 热99久久久久精品小说推荐| 青春草亚洲视频在线观看| 久久亚洲国产成人精品v| 国产日韩一区二区三区精品不卡| 欧美激情 高清一区二区三区| 人人澡人人妻人| 欧美激情极品国产一区二区三区| 少妇猛男粗大的猛烈进出视频| 亚洲欧美精品自产自拍| 欧美成人午夜精品| 一区二区av电影网| 美女午夜性视频免费| 夫妻午夜视频| 欧美午夜高清在线| 亚洲av电影在线观看一区二区三区| 美女高潮喷水抽搐中文字幕| 日本av手机在线免费观看| 女性生殖器流出的白浆| 亚洲欧美一区二区三区久久| 精品国产乱码久久久久久男人| 色老头精品视频在线观看| 一级,二级,三级黄色视频| 韩国高清视频一区二区三区| 香蕉国产在线看| 亚洲 国产 在线| 成人国产av品久久久| 黄网站色视频无遮挡免费观看|