• <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)機理與適用條件
    精品国内亚洲2022精品成人| 女生性感内裤真人,穿戴方法视频| 啦啦啦观看免费观看视频高清| 一本一本综合久久| 成人免费观看视频高清| 国产私拍福利视频在线观看| 日韩成人在线观看一区二区三区| 成人特级黄色片久久久久久久| 岛国在线观看网站| 国产精品亚洲一级av第二区| 欧美av亚洲av综合av国产av| 侵犯人妻中文字幕一二三四区| 国产亚洲欧美精品永久| 亚洲成av人片免费观看| 男女之事视频高清在线观看| 亚洲精品美女久久av网站| 亚洲片人在线观看| 国产野战对白在线观看| 午夜福利高清视频| 日韩中文字幕欧美一区二区| 免费在线观看视频国产中文字幕亚洲| 一级片免费观看大全| 午夜福利在线在线| 怎么达到女性高潮| 丁香六月欧美| 午夜成年电影在线免费观看| videosex国产| 每晚都被弄得嗷嗷叫到高潮| 男男h啪啪无遮挡| xxx96com| av福利片在线| 免费搜索国产男女视频| 久久中文字幕一级| 久久久久久久午夜电影| 亚洲,欧美精品.| 草草在线视频免费看| 一级a爱片免费观看的视频| 久久欧美精品欧美久久欧美| 久久久久国产精品人妻aⅴ院| 久久伊人香网站| 男人的好看免费观看在线视频 | 两个人看的免费小视频| 国产精华一区二区三区| 亚洲无线在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 九色国产91popny在线| e午夜精品久久久久久久| 麻豆成人午夜福利视频| 老鸭窝网址在线观看| 久久久久久久午夜电影| 久久中文看片网| 精品国产一区二区三区四区第35| 日本成人三级电影网站| 在线十欧美十亚洲十日本专区| 一级毛片精品| 国产激情欧美一区二区| 怎么达到女性高潮| 亚洲 欧美一区二区三区| 成人18禁在线播放| 亚洲五月婷婷丁香| 黄色视频,在线免费观看| 午夜免费观看网址| 搡老岳熟女国产| 亚洲真实伦在线观看| 亚洲一区二区三区色噜噜| 欧美国产精品va在线观看不卡| 亚洲国产欧美网| 丁香六月欧美| 成人欧美大片| 久久午夜综合久久蜜桃| √禁漫天堂资源中文www| 女人被狂操c到高潮| 亚洲欧美激情综合另类| 啦啦啦 在线观看视频| 日韩欧美一区视频在线观看| 久久久久国产精品人妻aⅴ院| 亚洲精品国产一区二区精华液| bbb黄色大片| 欧美在线一区亚洲| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲一区高清亚洲精品| 欧美日韩福利视频一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 99在线视频只有这里精品首页| 亚洲色图 男人天堂 中文字幕| 一区二区三区精品91| 99riav亚洲国产免费| 一区二区三区精品91| 好男人在线观看高清免费视频 | 窝窝影院91人妻| 日韩欧美在线二视频| 别揉我奶头~嗯~啊~动态视频| 18禁美女被吸乳视频| 午夜福利欧美成人| 亚洲 国产 在线| 少妇粗大呻吟视频| 曰老女人黄片| 日韩国内少妇激情av| 99久久99久久久精品蜜桃| 中文字幕av电影在线播放| 91老司机精品| 丰满人妻熟妇乱又伦精品不卡| 老司机福利观看| e午夜精品久久久久久久| 1024视频免费在线观看| 91国产中文字幕| 久久亚洲真实| 男女视频在线观看网站免费 | 久久久久久久久久黄片| 久久久国产精品麻豆| 51午夜福利影视在线观看| 久久久久久久精品吃奶| 男人的好看免费观看在线视频 | 变态另类丝袜制服| 亚洲久久久国产精品| 午夜免费观看网址| 国产一区二区激情短视频| videosex国产| 日韩三级视频一区二区三区| 亚洲国产欧洲综合997久久, | 中出人妻视频一区二区| 久久伊人香网站| 久久青草综合色| 热re99久久国产66热| 99国产精品一区二区三区| 日韩国内少妇激情av| 精品福利观看| 动漫黄色视频在线观看| 亚洲第一电影网av| 99热6这里只有精品| 午夜福利视频1000在线观看| 19禁男女啪啪无遮挡网站| 久久国产精品男人的天堂亚洲| 成人国产一区最新在线观看| www日本在线高清视频| 亚洲欧洲精品一区二区精品久久久| 国产亚洲欧美98| 欧美激情高清一区二区三区| 黄色视频不卡| 无遮挡黄片免费观看| 国产精品久久久久久亚洲av鲁大| 欧美一区二区精品小视频在线| 欧美亚洲日本最大视频资源| 久久久久久九九精品二区国产 | 久久久久精品国产欧美久久久| 两人在一起打扑克的视频| 可以在线观看的亚洲视频| 午夜精品在线福利| 91国产中文字幕| 国产色视频综合| 国产熟女午夜一区二区三区| 欧美人与性动交α欧美精品济南到| 欧美人与性动交α欧美精品济南到| 在线观看日韩欧美| 97超级碰碰碰精品色视频在线观看| 国产视频内射| 欧美成人免费av一区二区三区| 久久国产精品人妻蜜桃| 国产在线精品亚洲第一网站| 淫妇啪啪啪对白视频| 人人澡人人妻人| 精品一区二区三区av网在线观看| 可以在线观看毛片的网站| 色播在线永久视频| 久久久久免费精品人妻一区二区 | 又黄又粗又硬又大视频| 国产精品亚洲一级av第二区| 999久久久国产精品视频| 国产真人三级小视频在线观看| 精品乱码久久久久久99久播| 久久久久久大精品| 日韩视频一区二区在线观看| 亚洲av电影在线进入| 男女视频在线观看网站免费 | 午夜福利免费观看在线| 色综合站精品国产| 亚洲狠狠婷婷综合久久图片| 国产高清有码在线观看视频 | 久久亚洲精品不卡| 免费电影在线观看免费观看| 国产一区二区三区在线臀色熟女| 母亲3免费完整高清在线观看| 国产熟女午夜一区二区三区| 人妻久久中文字幕网| 免费观看精品视频网站| www国产在线视频色| x7x7x7水蜜桃| 国产精品久久电影中文字幕| 99久久精品国产亚洲精品| 国产亚洲欧美98| 制服诱惑二区| 国产精品 国内视频| 男女视频在线观看网站免费 | 国产精品日韩av在线免费观看| 亚洲av电影在线进入| 日本 欧美在线| 中文字幕人妻熟女乱码| 中出人妻视频一区二区| 黄片播放在线免费| 精品卡一卡二卡四卡免费| 亚洲午夜精品一区,二区,三区| 精品电影一区二区在线| 无遮挡黄片免费观看| 日韩欧美一区二区三区在线观看| 国产精品日韩av在线免费观看| 午夜福利免费观看在线| 精品免费久久久久久久清纯| 久久亚洲真实| 在线观看免费日韩欧美大片| 欧美色视频一区免费| 亚洲人成伊人成综合网2020| 免费看日本二区| 在线播放国产精品三级| 99国产精品99久久久久| 手机成人av网站| 老汉色∧v一级毛片| 人人妻人人澡欧美一区二区| 国产精品香港三级国产av潘金莲| 宅男免费午夜| 欧美日韩中文字幕国产精品一区二区三区| 大型av网站在线播放| 亚洲色图 男人天堂 中文字幕| 免费在线观看成人毛片| 午夜久久久在线观看| 成人精品一区二区免费| 亚洲一区中文字幕在线| 国内少妇人妻偷人精品xxx网站 | 高清毛片免费观看视频网站| 欧美中文日本在线观看视频| 久久久久久亚洲精品国产蜜桃av| 精品午夜福利视频在线观看一区| 亚洲av电影不卡..在线观看| 中文字幕高清在线视频| 熟女电影av网| 国产真人三级小视频在线观看| 亚洲av日韩精品久久久久久密| 69av精品久久久久久| 国产欧美日韩一区二区精品| 丰满人妻熟妇乱又伦精品不卡| 女人爽到高潮嗷嗷叫在线视频| 亚洲电影在线观看av| 天天躁狠狠躁夜夜躁狠狠躁| 久久久水蜜桃国产精品网| 母亲3免费完整高清在线观看| x7x7x7水蜜桃| 黄色a级毛片大全视频| 免费无遮挡裸体视频| 久久国产精品男人的天堂亚洲| av天堂在线播放| 男女下面进入的视频免费午夜 | 色综合亚洲欧美另类图片| 亚洲第一青青草原| 亚洲激情在线av| 精品久久久久久久毛片微露脸| 久久久国产精品麻豆| 亚洲性夜色夜夜综合| 国产精华一区二区三区| 可以免费在线观看a视频的电影网站| 精品久久久久久久久久免费视频| 国产精品国产高清国产av| 国产av在哪里看| 啦啦啦 在线观看视频| 桃红色精品国产亚洲av| 色综合欧美亚洲国产小说| 国产亚洲精品一区二区www| 大型黄色视频在线免费观看| 国产精品免费一区二区三区在线| 免费搜索国产男女视频| 国产97色在线日韩免费| 国产成人啪精品午夜网站| 成人精品一区二区免费| 亚洲人成网站在线播放欧美日韩| 两个人的视频大全免费| 精品久久久久久久久亚洲| 久久人人爽人人爽人人片va| 卡戴珊不雅视频在线播放| 三级毛片av免费| 中文字幕人妻熟人妻熟丝袜美| 国产精品不卡视频一区二区| 亚洲国产高清在线一区二区三| 看非洲黑人一级黄片| a级毛片a级免费在线| 99久国产av精品国产电影| 亚洲美女搞黄在线观看 | 久久亚洲精品不卡| 激情 狠狠 欧美| 午夜爱爱视频在线播放| 在线播放国产精品三级| 日韩精品中文字幕看吧| 波多野结衣高清无吗| 亚洲精品色激情综合| 插阴视频在线观看视频| 最新在线观看一区二区三区| 真人做人爱边吃奶动态| 久久6这里有精品| 天天躁日日操中文字幕| 国产v大片淫在线免费观看| 国产探花极品一区二区| 又粗又爽又猛毛片免费看| 国产一区亚洲一区在线观看| 男女视频在线观看网站免费| 色视频www国产| 国产在线男女| 91久久精品电影网| 精品久久国产蜜桃| 亚洲国产色片| 观看免费一级毛片| 亚洲成人av在线免费| 一a级毛片在线观看| 一个人免费在线观看电影| 日本撒尿小便嘘嘘汇集6| 天堂av国产一区二区熟女人妻| 日本黄大片高清| 亚洲自偷自拍三级| 欧美色视频一区免费| 最近最新中文字幕大全电影3| 国产乱人视频| 日日干狠狠操夜夜爽| 成人美女网站在线观看视频| av在线观看视频网站免费| 国产精品女同一区二区软件| 精品无人区乱码1区二区| 午夜免费男女啪啪视频观看 | 日本精品一区二区三区蜜桃| 国产一区二区三区在线臀色熟女| 一区二区三区免费毛片| 国产真实乱freesex| 一级av片app| 18+在线观看网站| 亚洲色图av天堂| 精品人妻视频免费看| 亚洲一级一片aⅴ在线观看| 成人漫画全彩无遮挡| 亚洲av免费高清在线观看| 美女被艹到高潮喷水动态| 身体一侧抽搐| 国产成人福利小说| 乱系列少妇在线播放| 在线国产一区二区在线| 亚洲精品亚洲一区二区| 亚洲精品国产成人久久av| 97超级碰碰碰精品色视频在线观看| 91在线观看av| 精品人妻视频免费看| 狠狠狠狠99中文字幕| 极品教师在线视频| 国产伦在线观看视频一区| 久久久成人免费电影| 午夜激情欧美在线| 日日摸夜夜添夜夜添av毛片| 日本五十路高清| 国产精品爽爽va在线观看网站| 国产黄色视频一区二区在线观看 | 成人av一区二区三区在线看| 菩萨蛮人人尽说江南好唐韦庄 | 欧美成人免费av一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成人精品中文字幕电影| 18禁在线播放成人免费| 国产精品一二三区在线看| 国产精品乱码一区二三区的特点| 自拍偷自拍亚洲精品老妇| 综合色av麻豆| 久久久久性生活片| 18+在线观看网站| 无遮挡黄片免费观看| 一夜夜www| 女人被狂操c到高潮| 亚洲最大成人中文| 日韩精品有码人妻一区| 一个人观看的视频www高清免费观看| 中文字幕免费在线视频6| 亚洲成人av在线免费| 欧美精品国产亚洲| 最近2019中文字幕mv第一页| 搞女人的毛片| 高清日韩中文字幕在线| 麻豆国产av国片精品| 国产精品野战在线观看| 久久久久久九九精品二区国产| 麻豆国产97在线/欧美| 久久人人爽人人片av| 亚洲无线观看免费| 精品久久久久久久久亚洲| 性插视频无遮挡在线免费观看| 真实男女啪啪啪动态图| 国内精品美女久久久久久| 女生性感内裤真人,穿戴方法视频| 麻豆乱淫一区二区| 日本熟妇午夜| 久久久久久久久大av| 国产淫片久久久久久久久| 亚洲欧美成人精品一区二区| av国产免费在线观看| 久久午夜福利片| 精品一区二区三区av网在线观看| 男女啪啪激烈高潮av片| а√天堂www在线а√下载| 亚洲国产精品国产精品| av在线蜜桃| 国产精品av视频在线免费观看| 搡老妇女老女人老熟妇| 国产视频内射| 18+在线观看网站| 精品人妻偷拍中文字幕| 久久久久久久久久黄片| 日韩欧美精品免费久久| 精品日产1卡2卡| 亚洲最大成人中文| 亚洲成人中文字幕在线播放| 少妇裸体淫交视频免费看高清| 亚洲成a人片在线一区二区| 亚洲经典国产精华液单| 国产亚洲91精品色在线| 国产日本99.免费观看| 久久鲁丝午夜福利片| 日韩欧美三级三区| 久久精品夜色国产| 人妻制服诱惑在线中文字幕| 久久精品国产清高在天天线| 午夜爱爱视频在线播放| 国产探花极品一区二区| 久久久久久久久大av| 黄色一级大片看看| 婷婷六月久久综合丁香| 久久久色成人| 精品少妇黑人巨大在线播放 | 亚洲国产精品成人久久小说 | 精品一区二区三区av网在线观看| 免费在线观看影片大全网站| 国产精品精品国产色婷婷| 高清毛片免费看| 成人毛片a级毛片在线播放| 国产伦一二天堂av在线观看| 国产高清不卡午夜福利| 天堂√8在线中文| 欧美3d第一页| 婷婷精品国产亚洲av| 狂野欧美激情性xxxx在线观看| 欧美另类亚洲清纯唯美| 人妻夜夜爽99麻豆av| 久久综合国产亚洲精品| 午夜久久久久精精品| 中文字幕精品亚洲无线码一区| 亚洲天堂国产精品一区在线| 午夜影院日韩av| 日本a在线网址| 校园春色视频在线观看| 变态另类成人亚洲欧美熟女| 日日撸夜夜添| 在线免费十八禁| 波多野结衣巨乳人妻| 久久久成人免费电影| 午夜免费激情av| 亚洲av成人精品一区久久| 国产高潮美女av| 激情 狠狠 欧美| 晚上一个人看的免费电影| 三级经典国产精品| 欧美一区二区亚洲| 亚洲三级黄色毛片| 日韩制服骚丝袜av| 国产美女午夜福利| 三级男女做爰猛烈吃奶摸视频| 久久久久久大精品| 又爽又黄a免费视频| 最好的美女福利视频网| 毛片女人毛片| 国产精品国产高清国产av| 男女边吃奶边做爰视频| 精品99又大又爽又粗少妇毛片| 男女啪啪激烈高潮av片| 日韩高清综合在线| 日日摸夜夜添夜夜添小说| 成人高潮视频无遮挡免费网站| 亚洲中文日韩欧美视频| 自拍偷自拍亚洲精品老妇| 亚洲精品国产av成人精品 | 亚洲美女黄片视频| 俺也久久电影网| 免费av观看视频| 色吧在线观看| 日日摸夜夜添夜夜添av毛片| 日韩制服骚丝袜av| 老司机午夜福利在线观看视频| 国产成人精品久久久久久| 日日摸夜夜添夜夜爱| 一级毛片aaaaaa免费看小| 精品国产三级普通话版| 久久99热6这里只有精品| 寂寞人妻少妇视频99o| 日日摸夜夜添夜夜爱| 在线播放无遮挡| 99久久中文字幕三级久久日本| 性欧美人与动物交配| 永久网站在线| 久久午夜福利片| 国产成人freesex在线 | 女人被狂操c到高潮| 两个人的视频大全免费| 男女视频在线观看网站免费| 18+在线观看网站| 亚洲自偷自拍三级| 欧美高清性xxxxhd video| 国产午夜精品久久久久久一区二区三区 | 人妻夜夜爽99麻豆av| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看 | 国产精品,欧美在线| 美女 人体艺术 gogo| 小说图片视频综合网站| 搡老妇女老女人老熟妇| 看片在线看免费视频| 国产精品人妻久久久影院| 夜夜夜夜夜久久久久| 在线天堂最新版资源| 一区二区三区四区激情视频 | 国产久久久一区二区三区| 深夜a级毛片| 少妇人妻精品综合一区二区 | 国产 一区 欧美 日韩| 久久这里只有精品中国| 女的被弄到高潮叫床怎么办| 成人av一区二区三区在线看| 91久久精品国产一区二区成人| 国产伦精品一区二区三区四那| 国内揄拍国产精品人妻在线| 亚洲人与动物交配视频| 午夜福利高清视频| 六月丁香七月| 国产午夜精品论理片| 亚洲性夜色夜夜综合| 最近最新中文字幕大全电影3| 国产男靠女视频免费网站| 欧美色视频一区免费| 又黄又爽又免费观看的视频| 97在线视频观看| 最近最新中文字幕大全电影3| 少妇人妻一区二区三区视频| 美女高潮的动态| 亚洲真实伦在线观看| 欧美中文日本在线观看视频| 男人舔奶头视频| 少妇人妻精品综合一区二区 | 国产在视频线在精品| 国产不卡一卡二| 啦啦啦啦在线视频资源| 亚洲av中文av极速乱| 最后的刺客免费高清国语| 精品福利观看| 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验 | 国产精品人妻久久久影院| 老女人水多毛片| 国产精品久久久久久精品电影| 久久精品国产亚洲av天美| 成人国产麻豆网| 亚洲av电影不卡..在线观看| 亚洲成人精品中文字幕电影| 久久久久性生活片| 又粗又爽又猛毛片免费看| 人妻少妇偷人精品九色| 欧美三级亚洲精品| 亚洲精品国产av成人精品 | 黄色视频,在线免费观看| 亚洲av中文字字幕乱码综合| 亚洲精品粉嫩美女一区| 亚洲三级黄色毛片| 五月玫瑰六月丁香| 精品久久久久久久久久免费视频| 我的女老师完整版在线观看| 老师上课跳d突然被开到最大视频| 免费av毛片视频| 国产精品永久免费网站| 久久婷婷人人爽人人干人人爱| 亚洲不卡免费看| 久久婷婷人人爽人人干人人爱| 国产91av在线免费观看| 久久亚洲精品不卡| 99热这里只有精品一区| 看非洲黑人一级黄片| 欧美成人精品欧美一级黄| 搡老岳熟女国产| 成人亚洲精品av一区二区| 精华霜和精华液先用哪个| 极品教师在线视频| 最新中文字幕久久久久| 高清毛片免费看| 草草在线视频免费看| 亚洲人成网站在线播| 男人舔女人下体高潮全视频| 午夜福利成人在线免费观看| 国产精品久久久久久精品电影| 中文字幕人妻熟人妻熟丝袜美| 久久久久久久久大av| 真实男女啪啪啪动态图| 国产精品人妻久久久影院| 免费观看精品视频网站| 国产一区二区激情短视频| 亚洲av中文字字幕乱码综合| 淫秽高清视频在线观看| 久久国内精品自在自线图片| 高清午夜精品一区二区三区 | 亚洲七黄色美女视频| 国产精品一区二区三区四区免费观看 | 99riav亚洲国产免费| 日本黄色视频三级网站网址| 亚洲av免费高清在线观看| 日韩欧美 国产精品| 久久久国产成人免费| 国产三级中文精品| 久久99热6这里只有精品|