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

    紫色土坡耕地裂隙潛流的產(chǎn)流機(jī)理與膠體顆粒遷移

    2017-12-16 05:44:19唐翔宇鮮青松
    環(huán)境科學(xué)研究 2017年12期

    張 維, 唐翔宇, 鮮青松

    1.重慶工商大學(xué)旅游與國土資源學(xué)院, 重慶 400067 2.中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所, 山地表生過程與生態(tài)調(diào)控重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610041 3.中國科學(xué)院大學(xué), 北京 100049

    紫色土坡耕地裂隙潛流的產(chǎn)流機(jī)理與膠體顆粒遷移

    張 維1, 唐翔宇2*, 鮮青松2,3

    1.重慶工商大學(xué)旅游與國土資源學(xué)院, 重慶 400067 2.中國科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所, 山地表生過程與生態(tài)調(diào)控重點(diǎn)實(shí)驗(yàn)室, 四川 成都 610041 3.中國科學(xué)院大學(xué), 北京 100049

    土壤膠體是坡耕地農(nóng)化物質(zhì)遷移的主要載體. 借助18O(氧同位素)示蹤技術(shù),探索了2014年8月29日和9月10日兩場(chǎng)降雨下大型紫色土坡耕地(1 500 m2)裂隙潛流產(chǎn)流的水源來源及過程特征,并耦合了膠體顆粒釋放與遷移機(jī)理的研究. 結(jié)果表明,裂隙潛流及膠體遷移的水文過程線均總體呈快速上升和長拖尾的特征. 隨裂隙潛流產(chǎn)流開始,雨水對(duì)潛流的貢獻(xiàn)逐漸增大,并在流量快速上升段支配裂隙潛流產(chǎn)流,而潛流流量峰值前及退水階段,土壤前期可動(dòng)水是潛流的主要產(chǎn)流來源. 兩場(chǎng)降雨下裂隙潛流中膠體顆粒濃度介于0.60~6.85 mgL之間,平均值分別為1.58和2.31 mgL,水浴超聲后膠體顆粒濃度平均值分別為原樣的2.15和1.81倍. 膠體顆粒遷移速率比產(chǎn)流速率快(>30 min),表明膠體輔助坡地農(nóng)化物質(zhì)遷移的潛力較大. 對(duì)于長歷時(shí)小降雨事件,潛流中膠體的遷移動(dòng)態(tài)受潛流水化學(xué)因素〔如ρ(DOC)、ρ(Mg2+)和EC(電導(dǎo)率)〕支配,而強(qiáng)降雨事件下,潛流中膠體顆粒濃度還與潛流流量呈極顯著負(fù)相關(guān)(R2>0.5). 此外,坡地內(nèi)部產(chǎn)流方式(橫向及垂向)對(duì)裂隙潛流中膠體顆粒的遷移通量有重要影響. 研究顯示,裂隙潛流產(chǎn)流過程線結(jié)合土水勢(shì)、18O及水化學(xué)指標(biāo)的動(dòng)態(tài)變化,能夠全面揭示裂隙潛流產(chǎn)流的階段特征以及膠體顆粒釋放與遷移的機(jī)理,對(duì)于進(jìn)一步研究膠體對(duì)磷、有機(jī)農(nóng)藥等憎水性農(nóng)化物的輔助運(yùn)移特征有重要意義.

    裂隙潛流; 產(chǎn)流; 膠體; 紫色土

    坡耕地是長江上游重要的耕地資源之一,但不合理的開墾以及該地區(qū)雨季集中大雨的耦合導(dǎo)致坡耕地水土流失問題突出. 紫色土中富含的大孔隙、介孔以及泥巖中普遍發(fā)育的微細(xì)裂隙是紫色土坡耕地產(chǎn)流和污染物遷移的主要通道. 土壤大孔隙和介孔通常指的是在土水勢(shì)為0~-10及-10~-100 cm時(shí)儲(chǔ)水的孔隙[1]. 大量研究證實(shí),紫色土坡耕地是長江上游及其支流侵蝕產(chǎn)沙與農(nóng)業(yè)非點(diǎn)源污染的主要來源,對(duì)上游尤其是三峽庫區(qū)的水環(huán)境與水生態(tài)安全構(gòu)成嚴(yán)重威脅[2-4].

    已有的紫色土坡耕地產(chǎn)流的研究集中在對(duì)不同徑流結(jié)果的特征分析及影響因素的探索. 降雨坡地特性及耕作管理方式對(duì)紫色土坡耕地產(chǎn)流形式及污染物協(xié)同遷移有重要影響[5-7]. 但產(chǎn)流機(jī)理及徑流水源來源的研究則鮮有報(bào)道. 穩(wěn)定性氫氧同位素能夠反映雨水入滲及產(chǎn)流信息. 相對(duì)于小流域,坡耕地更小的空間異質(zhì)性更有利于氫氧同位素示蹤技術(shù)在坡耕地水文過程的應(yīng)用[8]. 王超[9]通過18O對(duì)紫色土坡耕地徑流水源的示蹤研究表明,紫色土坡耕地以地面下徑流為主,是雨水與土壤前期可動(dòng)水混合產(chǎn)流的結(jié)果. 但產(chǎn)流過程仍需要結(jié)合其他指標(biāo)進(jìn)行全面描述. 不同水源的化學(xué)特性及污染物含量差異較大,進(jìn)而導(dǎo)致坡耕地污染物的協(xié)同遷移輸出產(chǎn)生差異. 泥沙以及氮磷等農(nóng)業(yè)非點(diǎn)源污染物是文獻(xiàn)中廣泛報(bào)道的紫色土坡耕地輸出物質(zhì)[10-12],膠體顆粒的遷移缺乏關(guān)注. 土壤膠體是多種污染物遷移的載體,對(duì)農(nóng)化物質(zhì)如氮、磷、農(nóng)藥等具有較強(qiáng)的輔助運(yùn)移能力[13-14],同時(shí)生物膠體(細(xì)菌、病毒)本身就是污染物,對(duì)淺層地下水環(huán)境安全構(gòu)成威脅. 膠體遷移的研究集中在柱試驗(yàn)尺度人工合成膠體的遷移理論與模型擬合方面[15],坡地尺度自然膠體的遷移尤其是與產(chǎn)流過程的結(jié)合缺乏報(bào)道. 紫色土具有較高的顆粒分散性,為坡耕地提供大量的自然膠體來源. 因此,該研究以紫色土坡耕地(1 500 m2)為對(duì)象,基于18O示蹤并結(jié)合其他化學(xué)指標(biāo)的動(dòng)態(tài)變化,探索裂隙潛流產(chǎn)流機(jī)理,同時(shí)耦合產(chǎn)流過程中土壤膠體顆粒的釋放與遷移機(jī)理的研究,對(duì)于進(jìn)一步揭示膠體輔助污染物運(yùn)移機(jī)制及流域水環(huán)境保護(hù)有重要意義.

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于四川省鹽亭縣林山鄉(xiāng)截流小流域(31°16′N、105°28′E),屬中亞熱帶季風(fēng)氣候,多年(1981—2006年)平均降水量為826 mm,夏季降雨占65.5%[16]. 坡耕地是流域內(nèi)最主要的土地利用類型(>40%),以10° 以下的緩坡為主. 流域內(nèi)土壤主要包括石灰性紫色土及水稻土.

    1.2 研究方法

    1.2.1紫色土坡耕地原位監(jiān)測(cè)系統(tǒng)

    選取流域內(nèi)一塊典型的紫色土坡耕地(見圖1)為研究對(duì)象. 該坡地平均坡度為7°,與流域內(nèi)自然分布的坡耕地類似. 紫色土為泥巖快速風(fēng)化形成的幼年土,pH為8.3,w(有機(jī)質(zhì))平均值為8.75 gkg,w(砂礫)、w(粉粒)、w(黏粒)分別為27.1%、51.6%及22.3%. 紫色泥巖的容重介于2.0~3.0 gcm3之間,孔隙度>10%,微細(xì)裂隙普遍發(fā)育.

    圖1 紫色土坡耕地監(jiān)測(cè)示意Fig.1 Sketch map of the monitored sloping farmland

    坡耕地四周澆筑水泥墻隔絕外來水源,中部外圍安裝翻斗式雨量筒(分辨率0.1 mm). 坡地中部中線位置安裝張力計(jì)測(cè)定土水勢(shì),深度分別為15 cm(耕作層)、45 cm(近土巖界面). 坡地下方的橫截面設(shè)置裂隙潛流收集槽,并通過導(dǎo)流管引入定制的流量翻斗. 觀測(cè)發(fā)現(xiàn),泥巖與砂巖界面的裂隙潛流產(chǎn)流最頻繁,流量最多. 因此,重點(diǎn)關(guān)注裂隙潛流的產(chǎn)流過程及膠體遷移的特征. 雨量筒、張力計(jì)以及流量翻斗均連線至數(shù)采系統(tǒng)(CR1000,Campbell,USA)實(shí)現(xiàn)數(shù)據(jù)的自動(dòng)采集與記錄.

    1.2.2樣品采集及分析

    次降雨雨水的采集方法同文獻(xiàn)[9]. 對(duì)于裂隙潛流,從降雨開始于導(dǎo)流管末端人工采樣. 徑流水文過程上升段,采樣間隔為15 min,退水段則為30或60 min. 該研究采集了2014年8月26日及9月10日兩次降雨事件下的雨水及裂隙潛流樣品.

    雨水樣品經(jīng)0.45 μm濾膜過濾后用液態(tài)水同位素分析儀(i2120,Picarro Inc.,USA)測(cè)定18O的水平(以δ18O表示),測(cè)定精度為0.1‰. 以階段性降水量的δ18O加權(quán)均值作為次降雨的δ18O值,即:

    (1)

    式中:Ri為階段性(30 min)降水量,mm;δi為相應(yīng)雨水的δ18O值,‰;n為雨水樣品個(gè)數(shù).

    對(duì)于潛流樣品,采用復(fù)合酸度計(jì)(SensIon+MM150,USA)測(cè)定水樣pH及EC(電導(dǎo)率). 膠體顆粒濃度采用紫外-可見分光光度計(jì)(Tu-1810,China)于400 nm下比色測(cè)定. 潛流樣品水浴超聲振蕩(100 W)2 min后再一次測(cè)定膠體顆粒濃度. 膠體顆粒的數(shù)量通過式(2)估算[17]:

    Nc=6Mπd3ρ

    (2)

    式中:Nc為單位體積水樣中膠體顆粒數(shù)量,L-1;M為膠體顆粒質(zhì)量濃度,mgL;d為膠體顆粒當(dāng)量粒徑,該研究采用2 μm;ρ為膠體顆粒密度,黏土礦物為2.60 kgcm. 水樣經(jīng)0.45 μm濾膜過濾后用流動(dòng)分析儀(Auto Analyzer 3,Germany)測(cè)定ρ(DOC)(DOC為溶解性有機(jī)碳),用離子色譜(ICS-900,USA)測(cè)定ρ(Ca2+)和ρ(Mg2+).

    2 結(jié)果與分析

    2.1 裂隙潛流產(chǎn)流特征

    該研究監(jiān)測(cè)了兩場(chǎng)不同雨型(見表1)下坡耕地裂隙潛流與膠體遷移的過程變化.

    表1 降雨事件特征描述

    2014年8月26日降雨前,裂隙潛流背景流量僅為0.10 Lmin(見圖2),潛流響應(yīng)時(shí)間約為1.5 h. 此后,潛流流量逐漸上升并于降雨結(jié)束后達(dá)到峰值(2.24 Lmin). 降雨開始至Imax(最大雨強(qiáng),圖2虛線標(biāo)注),坡地中部耕作層和近土巖界面層的土水勢(shì)分別維持在-240及-780 cm,土壤大孔及介孔均處于排空狀態(tài).Imax出現(xiàn)至降雨接近結(jié)束時(shí),耕作層土水勢(shì)迅速升至-11 cm,而近土巖界面層土水勢(shì)保持穩(wěn)定,表明淺層土壤中介孔成為產(chǎn)流的主要通道,且坡地入滲以橫向的水平遷移為主. 降雨結(jié)束后潛流峰值及退水初期,耕作層土水勢(shì)逐漸降低,而近土巖界面土水勢(shì)由-780 cm快速升至-31 cm,表明此時(shí)裂隙潛流的補(bǔ)給主要為垂向的入滲水流,并且土壤介孔是主要的入滲通道.

    注: 虛線表示該日降雨Imax出現(xiàn)的時(shí)間.圖2 2014年8月26日降雨下裂隙潛流產(chǎn)流及膠體顆粒遷移特征Fig.2 Characteristics of fracture flow generation and colloid transport following the rain event on August 26th,2014

    2014年9月10日降雨是夏季典型的短歷時(shí)較強(qiáng)降雨事件(見表1). 裂隙潛流響應(yīng)時(shí)間約為30 min(見圖3),流量迅速上升至峰值(7.66 Lmin). 流量快速上升段,耕作層和近土巖界面層土水勢(shì)交替保持在-8 cm. 此時(shí),大孔隙部分充滿水,同時(shí),土壤介孔也成為產(chǎn)流通道,因而潛流流量增長速率較快. 而在退水段,兩個(gè)層次土水勢(shì)均低于-15 cm,僅介孔參與產(chǎn)流,導(dǎo)致退水過程相對(duì)較慢而呈現(xiàn)長拖尾特征. 退水段淺層和深層土水勢(shì)降低趨勢(shì)較一致,表明該次降雨事件下土層中產(chǎn)流以垂向遷移為主,可歸結(jié)為入滲水流以活塞流的形式向下趨替土壤前期可動(dòng)水的運(yùn)移,這與該階段潛流δ18O值的緩慢上升相對(duì)應(yīng). 兩次降雨下坡地產(chǎn)流方式的差異可能與降雨特性及坡地前期干旱情況有關(guān). ZHAO等[18]在該地區(qū)的的研究也證實(shí)前期濕潤以及高強(qiáng)度降雨有利于活塞流的產(chǎn)生.

    9月10日降雨Imax(圖3虛線標(biāo)注)作用下,潛流EC由621 μScm快速降至最小值(491 μScm),δ18O值也出現(xiàn)類似情況,但δ18O最小值滯后EC最小值30 min. 潛流δ18O值快速向雨水值δ18O(-16.99‰)迫近過程中,雨水對(duì)產(chǎn)流的貢獻(xiàn)越來越大并支配流量增長速率最快階段潛流的產(chǎn)流. 結(jié)合該階段土水勢(shì)的變化,可以推斷,在大孔隙部分充滿水、介孔幾乎全部充滿水且孔隙表面存在連續(xù)水膜的情況下,雨水可能繞過土壤基質(zhì),優(yōu)先通過大孔隙入滲,并經(jīng)過巖石裂隙快速形成裂隙潛流. 與前一場(chǎng)降雨響應(yīng)不同的是,潛流ρ(DOC)出現(xiàn)了兩個(gè)較明顯的濃度峰.ρ(DOC)最大值(3.45 mgL)出現(xiàn)在裂隙潛流流量快速上升段,為土壤大孔及介孔內(nèi)壁有機(jī)碳溶出的結(jié)果. 而在退水過程中段,ρ(DOC)出現(xiàn)一個(gè)次峰(2.46 mgL),可推斷為由慢速排水的介孔內(nèi)壁及基質(zhì)中有機(jī)碳溶出遷移引起. 因此,產(chǎn)流水文過程線結(jié)合土水勢(shì)、產(chǎn)流的δ18O值及水化學(xué)指標(biāo)〔如ρ(DOC)、EC〕的動(dòng)態(tài)變化,能夠全面揭示坡耕地裂隙潛流產(chǎn)流的機(jī)理,進(jìn)而有助于分析坡耕地農(nóng)化物質(zhì)的階段性遷移輸出特征.

    2.2 裂隙潛流中膠體遷移特征

    膠體顆粒濃度及數(shù)量變化如表2所示. 由表2可見,潛流中膠體顆粒濃度介于0.60~6.85 mgL之間,呈中等程度變異〔Cv(變異系數(shù))>0.5〕,與土層、巖層的空間異質(zhì)性及降雨強(qiáng)度的變化有關(guān). 膠體顆粒數(shù)量濃度介于5.51×107~6.29×108L-1之間,比文獻(xiàn)[19]中報(bào)道的其他裂隙潛流中相應(yīng)值低2~3個(gè)數(shù)量級(jí).

    注:虛線表示該日降雨Imax出現(xiàn)的時(shí)間.圖3 2014年9月10日降雨下裂隙潛流產(chǎn)流及膠體顆粒遷移特征Fig.3 Characteristics of fracture flow generation and colloid transport following the rain event on September 10th,2014

    降雨時(shí)間膠體顆粒濃度∕(mg∕L)最大值最小值平均值標(biāo)準(zhǔn)差Cv膠體顆粒數(shù)量濃度∕(107L-1)最大值最小值2014-08-264.580.601.580.8770.55342.15.512014-09-106.850.722.311.3160.56962.96.61

    從圖2、3可見,潛流產(chǎn)流開始至流量初始快速上升段,膠體顆粒濃度快速增至峰值,此后呈快速降低和長拖尾的變化特征.2014年8月26日降雨前,潛流中膠體顆粒背景值為0.65 mgL. 隨著產(chǎn)流進(jìn)行,膠體顆粒濃度迅速上升并在Imax后達(dá)到峰值(4.58 mgL),為背景濃度的7倍. 此后膠體顆粒濃度快速降低,但在降雨次峰作用下出現(xiàn)一個(gè)小峰(3.03 mgL).9月10日降雨過程兩次明顯較大降雨強(qiáng)度出現(xiàn)后,潛流中膠體顆粒濃度分別達(dá)到2.82和6.85 mgL. 因此,降雨強(qiáng)度及其分配對(duì)裂隙潛流中膠體顆粒的遷移輸出動(dòng)態(tài)有重要影響. 此外,膠體顆粒濃度的高值區(qū)主要分布在潛流初始響應(yīng)流量較低段. 兩場(chǎng)降雨潛流中膠體顆粒濃度峰均超前潛流流量峰>3 h,表明膠體顆粒遷移速率比潛流產(chǎn)流速率快. 其他學(xué)者通過柱試驗(yàn)也普遍證實(shí)了膠體顆粒比水流運(yùn)動(dòng)具有更早的穿透時(shí)間[20-22].

    裂隙潛流樣品超聲振蕩前后膠體顆粒濃度的變化如圖4所示. 超聲振蕩后膠體顆粒的濃度均有所增加,兩場(chǎng)降雨潛流樣品CaCb(超聲振蕩后的膠體顆粒濃度原樣品中膠體顆粒濃度)的平均值分別為2.15和1.81. 一方面,土壤中的膠體顆??赡芤蚰z結(jié)物(如碳酸鹽、赤鐵礦)的黏結(jié)作用以聚合物的形式存在[23]. 另一方面,潛流中ρ(Ca2+)遠(yuǎn)高于紫色土膠體的臨界絮凝濃度(24 mgL)[1],利于膠體聚合物和分散的膠體顆粒交互團(tuán)聚形成更大的顆粒而沉降[24]. 超聲振蕩可能破壞膠體間的物理黏結(jié)狀態(tài),導(dǎo)致因沉降的大顆粒的解體并重新分散懸浮在水中,進(jìn)而增加膠體顆粒的濃度.

    圖4 超聲振蕩后膠體顆粒濃度相對(duì)原樣中濃度的變化Fig.4 Relative variations of colloid concentration before and after ultrasound vibration

    3 討論

    土壤膠體的釋放主要與固液界面作用力的變化有關(guān). 界面引力包括水化力、范德華力以及化學(xué)鍵約束力等,界面斥力主要包括靜電斥力. 此外,還有滲流的水力剪切力對(duì)膠體顆粒的作用. 只有當(dāng)靜電斥力與剪切力的合力矩大于引力的合力矩時(shí),膠體顆粒才可能從土壤表面釋放進(jìn)而遷移[15,25]. 滲流的化學(xué)特征變化可能引起膠體表面化學(xué)特性的變化,進(jìn)而影響膠體釋放. 表3列出了潛流中膠體顆粒濃度與水力學(xué)及水化學(xué)指標(biāo)的相關(guān)分析結(jié)果.

    表3 裂隙潛流中膠體顆粒濃度與水力學(xué)及水化學(xué)指標(biāo)的皮爾遜相關(guān)系數(shù)

    注: ** 表示相關(guān)性在0.01水平上顯著(雙尾檢驗(yàn)).

    由表3可見,對(duì)于8月26日的長歷時(shí)小雨,潛流中膠體顆粒濃度與流量不存在相關(guān)性;而9月10日短歷時(shí)較強(qiáng)降雨下,二者呈極顯著負(fù)相關(guān). 對(duì)于水化學(xué)指標(biāo),膠體顆粒濃度與EC及ρ(Mg2+)均呈極顯著負(fù)相關(guān),而與ρ(DOC)呈極顯著正相關(guān),這與該地區(qū)井水中的研究結(jié)果[26]類似. 一方面,潛流響應(yīng)初始階段,水力剪切力增加. 另一方面,EC和ρ(Mg2+)的降低會(huì)導(dǎo)致膠體顆粒的擴(kuò)散雙電層的擴(kuò)張,進(jìn)而增加孔隙內(nèi)壁與分散的膠體顆粒間的排斥能阻,促進(jìn)膠體顆粒的釋放[15,27-28]. 同時(shí),ρ(DOC)的增加可能導(dǎo)致更多的有機(jī)物大分子吸附在膠體顆粒表面,增強(qiáng)表面負(fù)電荷[29],進(jìn)而增大固液界面的靜電斥力. 因此,水力剪切力以及靜電斥力的合力矩增大,超過引力的合力矩,從而導(dǎo)致土壤膠體的釋放. 但裂隙潛流流量對(duì)大雨強(qiáng)的響應(yīng)具有一定滯后作用. 此外,降雨前坡耕地土壤中分散的易遷移膠體顆粒是有限的[30]. 在產(chǎn)流初始階段孔隙內(nèi)壁易遷移膠體顆??焖俸谋M情況下,膠體的補(bǔ)給主要依賴于土壤基質(zhì)內(nèi)部分散的膠體顆粒向孔隙表面的慢速擴(kuò)散[31],從而導(dǎo)致流量峰值及退水段膠體顆粒濃度較低且長拖尾的現(xiàn)象. 因此,對(duì)于長歷時(shí)小降雨事件,膠體顆粒的釋放與遷移受潛流水力學(xué)特性的影響有限,而主要與水化學(xué)性質(zhì)變化有關(guān). 相反,短歷時(shí)的較強(qiáng)降雨事件下,潛流中膠體顆粒的遷移是水力學(xué)和水化學(xué)因素共同作用的結(jié)果.

    4 結(jié)論

    a) 裂隙潛流產(chǎn)流初期,雨水對(duì)潛流貢獻(xiàn)逐漸增大并在流量快速上升段支配裂隙潛流產(chǎn)流,而潛流峰值及退水段,土壤前期可動(dòng)水是潛流的主要來源.

    b) 潛流中膠體顆粒遷移速率比產(chǎn)流速率快(>30 min),膠體顆粒濃度介于0.60~6.85 mgL之間,呈中等程度變異(Cv>0.5). 水浴超聲能引起膠體顆粒濃度增加,膠體團(tuán)聚程度較高.

    c) 強(qiáng)降雨下,膠體顆粒濃度與潛流流量、EC及ρ(Mg2+)呈極顯著負(fù)相關(guān)(R2>0.5),與ρ(DOC)極顯著正相關(guān)(R2>0.6). 而長歷時(shí)小降雨事件下,膠體的遷移動(dòng)態(tài)受潛流水化學(xué)因素支配.

    d) 裂隙潛流產(chǎn)流過程線結(jié)合土水勢(shì)、產(chǎn)流的18O及水化學(xué)指標(biāo)〔如ρ(DOC)、EC〕的動(dòng)態(tài)變化,能夠全面揭示裂隙潛流產(chǎn)流的階段特征以及膠體顆粒釋放與遷移的機(jī)理,有助于進(jìn)一步研究膠體對(duì)磷、有機(jī)農(nóng)藥等憎水性農(nóng)化物的輔助運(yùn)移特征.

    [1] ZHANG Wei,TANG Xiangyu,WEISBROD N,etal.A coupled field study of subsurface fracture flow and colloid transport[J].Journal of Hydrology,2015,524:476-488.

    [2] GAO Yang,ZHU Bo,WANG Tao,etal.Seasonal change of non-point source pollution-induced bioavailable phosphorus loss:a case study of southwestern China[J].Journal of Hydrology,2012,420421:373-379.

    [3] SHEN Zhenyao,CHEN Lei,DING Xiaowen,etal.Long-term variation (1960-2003) and causal factors of non-point-source nitrogen and phosphorus in the upper reach of the Yangtze River[J].Journal of Hazardous Materials,2013,252:45-56.

    [4] TANG Jialiang,WANG Tao,ZHU Bo,etal.Tempo-spatial analysis of water quality in tributary bays of the Three Gorges Reservoir Region (China)[J].Environmental Science and Pollution Research,2015,22(21):16709-16720.

    [5] TANG Jialiang,WANG Tao,ZHU Bo,etal.Subsurface flow processes in sloping cropland of purple soil[J].Journal of Mountain Science,2012,9(1):1-9.

    [6] 鄭子成,秦鳳,李廷軒.不同坡度下紫色土地表微地形變化及其對(duì)土壤侵蝕的影響[J].農(nóng)業(yè)工程學(xué)報(bào),2015,31(8):168-175.

    ZHENG Zicheng,QIN Feng,LI Tingxuan.Changes in soil surface microrelief of purple soil under different slope gradients and its effects on soil erosion[J].Transactions of Chinese Society of Agricultural Engineering,2015,31(8):168-175.

    [7] TANG Jialiang,CHEN Xunqiang,ZHU Bo,etal.Rainfall and tillage impacts on soil erosion of sloping cropland with subtropical monsoon climate:a case study in hilly purple soil area,China[J].Journal of Mountain Science,2015,12(1):134-144.

    [8] BIRKEL C,SOULSBY C.Linking tracers,water age and conceptual models to identify dominant runoff processes in a sparsely monitored humid tropical catchment[J].Hydrological Processes,2016,30(24):4477-4493.

    [9] 王超.紫色土坡耕地徑流水源的同位素示蹤研究[D].北京:中國科學(xué)院大學(xué),2013.

    [10] ZHU Bo,WANG Tao,KUANG Fuhong,etal.Measurements of nitrate leaching from a hillslope cropland in the central Sichuan Basin,China[J].Soil Science Society of America Journal,2009,73(4):1419-1426.

    [11] GAO Yang,ZHU Bo,HE Nianping,etal.Phosphorus and carbon competitive sorption-desorption and associated non-point loss respond to natural rainfall events[J].Journal of Hydrology,2014,517:447-457.

    [12] LI Zhanbin,LI Peng,HAN Jiangang,etal.Sediment flow behavior in agro-watersheds of the purple soil region in China under different storm types and spatial scales[J].Soil & Tillage Research,2009,105(2):285-291.

    [13] GAO Changjuan,YANG Shengke,WANG Wenke,etal.Influence of humic acid colloid on adsorption of DDT in the riverbed sediments[J].Asian Journal of Chemistry,2014,26(9):2637-2642.

    [14] PANG Liping,LAFOGLER M,KNORR B,etal.Influence of colloids on the attenuation and transport of phosphorus in alluvial gravel aquifer and vadose zone media[J].Science of the Total Environment,2016,550:60-68.

    [15] ZHANG Wei,TANG Xiangyu,WEISBROD N,etal.A review of colloid transport in fractured rocks[J].Journal of Mountain Science,2012,9(6):770-787.

    [16] ZHAO Pei,TANG Xiangyu,ZHAO Peng,etal.Tracing water flow from sloping farmland to streams using oxygen-18 isotope to study a small agricultural catchment in southwest China[J].Soil & Tillage Research,2013,134:180-194.

    [17] RANI R D,SASIDHAR P.Stability assessment and characterization of colloids in coastal groundwater aquifer system at Kalpakkam[J].Environmental Earth Sciences,2011,62(2):233-243.

    [18] ZHAO Pei,TANG Xiangyu,ZHAO Peng,etal.Identifying the water source for subsurface flow with deuterium and oxygen-18 isotopes of soil water collected from tension lysimeters and cores[J].Journal of Hydrology,2013,503:1-10.

    [19] DEGUELDRE C,BAEYENS B,GOERLICH W,etal.Colloids in water from a subsurface fracture in granitic rock,Grimsel Test Site,Switzerland[J].Geochimicaet Cosmochimica Acta,1989,53(3):603-610.

    [20] TANG Xiangyu,WEISBROD N.Colloid-facilitated transport of lead in natural discrete fractures[J].Environmental Pollution,2009,157:2266-2274.

    [21] ZVIKELSKY O,WEISBROD N,DODY A.A comparison of clay colloid and artificial microsphere transport in natural discrete fractures[J].Journal of Colloid and Interface Science,2008,323(2):286-292.

    [22] 李海明,李云,翟菁,等.咸水中膠體遷移-沉積對(duì)砂介質(zhì)滲透性損失的數(shù)學(xué)模型[J].環(huán)境科學(xué)研究,2010,23(3):346-349.

    LI Haiming,LI Yun,ZHAI Jing,etal.Mathematical model of colloid transport-deposition on permeability loss of sand medium in salt water[J].Research of Environmental Sciences,2010,23(3):346-349.

    [23] WEISBROD N,DAHAN O,ADAR E M.Particle transport in unsaturated fractured chalk under arid conditions[J].Journal of Contaminant Hydrology,2002,56:117-136.

    [24] LI Haiming,WEI Jinbu,GE Yaochao,etal.The influence of Na+and Ca2+on the migration of colloids orand ammonia nitrogen in an unsaturated zone medium[J].Journal of Contaminant Hydrology,2016,194:24-29.

    [25] TORKZABAN S,KIM H N,SIMUNEK J,etal.Hysteresis of colloid retention and release in saturated porous media during transients in solution chemistry[J].Environmental Science & Technology,2010,44(5):1662-1669.

    [26] 張維,唐翔宇,鮮青松.紫色土小流域淺層井水中膠體顆粒的季節(jié)變化[J].環(huán)境科學(xué),2017,37(1):87-94.

    ZHANG Wei,TANG Xiangyu,XIAN Qingsong.Seasonal variation of colloid particles in the shallow well water of a small watershed of purple soil[J].Environmental Science,2017,37(1):87-94.

    [27] BRADFORD S A,KIM H.Implications of cation exchange on clay release and colloid-facilitated transport in porous media[J].Journal of Environmental Quality,2010,39(6):2040-2046.

    [28] 呂俊佳,許端平,李發(fā)生.不同環(huán)境因子對(duì)黑土膠體在飽和多孔介質(zhì)中運(yùn)移特性的影響[J].環(huán)境科學(xué)研究,2012,25(8):875-881.

    Lü Junjia,XU Duanping,LI Fasheng.Effects of different environmental factors on the transportation of black soil colloid in saturated porous media[J].Research of Environmental Sciences,2012,25(8):875-881.

    [29] TARASEVICH Y I,DOLENKO S A,TRIFONOVA M Y,etal.Association and colloid-chemical properties of humic acids in aqueous solutions[J].Colloid Journal,2013,75(2):207-213.

    [30] MAJDALANI S,MICHEL E,DI-PIETRO L,etal.Effects of wetting and drying cycles on in situ soil particle mobilization[J].European Journal of Soil Science,2008,59(2):147-155.

    [31] MOHANTY S K,SAIERS J E,RYAN J N.Colloid mobilization in a fractured soil:effect of pore-water exchange between preferential flow paths and soil matrix[J].Environmental Science & Technology,2016,50(5):2310-2317.

    [32] PORUBCAN A A,XU Shangbo.Colloid straining within saturated heterogeneous porous media[J].Water Research,2011,45(4):1796-1806.

    [33] KNAPPENBERGER T,ARAMRAK S,FLURY M.Transport of barrel and spherical shaped colloids in unsaturated porous media[J].Journal of Contaminant Hydrology,2015,180:69-79.

    MechanismsofFractureFlowGenerationandColloidTransportinaPurpleSoilSlopingFarmland

    ZHANG Wei1, TANG Xiangyu2*, XIAN Qingsong2,3

    1.School of Tourism and Land Resource, Chongqing Technology and Business University, Chongqing 400067, China2.Key Laboratory of Mountain Surface Processes and Ecological Regulation, Institute of Mountain Hazards and Environment, Chinese Academy of Sciences, Chengdu 610041, China3.University of Chinese Academy of Sciences, Beijing 100049, China

    Soil colloids mobilized from a sloping farmland during the flow process can act as a carrier of agricultural chemicals. On a purple soil sloping farmland (1500 m2) in Sichuan Basin, fracture flow water sources were identified using the18O tracing technique, and the dynamics of colloid transport were explored for two rain events on August 29thand September 10th, 2014. The results showed that, in response to the rain events, both hydrograph and colloid concentrations of fracture flow showed a pattern of early rapid rising to a peak followed by slow decrease with a long tail. Rainwater′s contribution to flow increased at the early stage of fracture flow, while pre-event mobile soil water appeared to be the main water source of fracture flow at peak flow discharge and the receding stage. Colloid concentration in the fracture flow varied from 0.60 to 6.85 mgL, with average concentrations of 1.58 and 2.31 mgL, respectively, during the two events. Ultrasonic dispersion treatment of the water samples resulted in an average increase of colloid concentration of 1.15 and 0.81 times for the two rain events, indicating a high degree of the aggregation of mobilized colloids. For small rain events of long duration, fracture flow chemistry (e.g.ρ(DOC),ρ(Mg2+) and EC) dominated colloid transport dynamics;however, for the heavy storms, fracture flow rate also had a strong influence on colloid transport (R2>0.5). In addition, soil water movement pattern (laterally or downward) affected the discharge of colloid from the sloping farmland. The combined use of fracture flow hydrograph′s response to soil water potential dynamics,18O and conservative flow tracers is an effective way to reveal the in-depth mechanisms of fracture flow generation and colloid transport, and to support the quantification of colloid-facilitated transport of agricultural chemicals (e.g., phosphorus and pesticides).

    fracture flow; flow generation; colloid; purple soil

    2017-03-23

    2017-09-07

    國家自然科學(xué)基金項(xiàng)目(41601539,41471268)

    張維(1987-),男,重慶人,講師,博士,主要從事環(huán)境科學(xué)與污染物水文學(xué)研究,zw512119@163.com.

    *責(zé)任作者,唐翔宇(1972-),男,福建莆田人,研究員,博士,主要從事環(huán)境科學(xué)研究,xytang@imde.ac.cn

    張維,唐翔宇,鮮青松.紫色土坡耕地裂隙潛流的產(chǎn)流機(jī)理與膠體顆粒遷移[J].環(huán)境科學(xué)研究,2017,30(12):1919-1926.

    ZHANG Wei,TANG Xiangyu,XIAN Qingsong.Mechanisms of fracture flow generation and colloid transport in a purple soil sloping farmland[J].Research of Environmental Sciences,2017,30(12):1919-1926.

    P345

    1001-6929(2017)12-1919-08

    A

    10.13198j.issn.1001-6929.2017.03.37

    三级毛片av免费| bbb黄色大片| 日本三级黄在线观看| 亚洲 欧美一区二区三区| 美国免费a级毛片| 变态另类成人亚洲欧美熟女| 久久久久久九九精品二区国产 | 精品一区二区三区av网在线观看| 国产久久久一区二区三区| 午夜久久久久精精品| 国产精品一区二区免费欧美| 中文在线观看免费www的网站 | 人人妻人人看人人澡| 香蕉久久夜色| 欧美国产精品va在线观看不卡| 国产黄片美女视频| 欧美亚洲日本最大视频资源| 日本五十路高清| 久久亚洲精品不卡| 亚洲一码二码三码区别大吗| 国产单亲对白刺激| 在线观看舔阴道视频| 亚洲国产高清在线一区二区三 | 两人在一起打扑克的视频| 国产激情久久老熟女| 亚洲欧洲精品一区二区精品久久久| 观看免费一级毛片| 成人欧美大片| 黄色a级毛片大全视频| 精品国产国语对白av| 亚洲一区中文字幕在线| 亚洲中文字幕日韩| 99国产极品粉嫩在线观看| 久久精品夜夜夜夜夜久久蜜豆 | av天堂在线播放| 国产精品香港三级国产av潘金莲| 香蕉av资源在线| 国产精品日韩av在线免费观看| 日本a在线网址| 视频在线观看一区二区三区| 亚洲精华国产精华精| 真人一进一出gif抽搐免费| 国产亚洲欧美98| 国产精品 欧美亚洲| 国产精品野战在线观看| 亚洲国产精品合色在线| 亚洲精品在线美女| 99国产极品粉嫩在线观看| 国产主播在线观看一区二区| 老汉色av国产亚洲站长工具| 制服丝袜大香蕉在线| 这个男人来自地球电影免费观看| 精品第一国产精品| 午夜两性在线视频| 在线免费观看的www视频| 又大又爽又粗| 啦啦啦韩国在线观看视频| 999精品在线视频| 男女床上黄色一级片免费看| 久99久视频精品免费| 级片在线观看| 国产亚洲精品综合一区在线观看 | 超碰成人久久| 亚洲欧洲精品一区二区精品久久久| 啦啦啦观看免费观看视频高清| 欧美乱码精品一区二区三区| 免费无遮挡裸体视频| 亚洲精品国产区一区二| 18禁观看日本| 亚洲人成77777在线视频| 久久香蕉激情| 亚洲av片天天在线观看| 亚洲国产精品999在线| 丝袜美腿诱惑在线| 特大巨黑吊av在线直播 | 一个人观看的视频www高清免费观看 | 亚洲av片天天在线观看| 国产精品久久久久久亚洲av鲁大| 久久中文字幕人妻熟女| 97超级碰碰碰精品色视频在线观看| 色播在线永久视频| 午夜福利18| 免费女性裸体啪啪无遮挡网站| 国内精品久久久久久久电影| 久久精品亚洲精品国产色婷小说| 久久精品夜夜夜夜夜久久蜜豆 | 一本一本综合久久| 国产激情偷乱视频一区二区| 女人高潮潮喷娇喘18禁视频| 成人午夜高清在线视频 | 美女高潮喷水抽搐中文字幕| 十八禁网站免费在线| 黄色视频,在线免费观看| 亚洲一卡2卡3卡4卡5卡精品中文| 麻豆一二三区av精品| 999久久久国产精品视频| 在线天堂中文资源库| 久9热在线精品视频| 免费在线观看日本一区| 精品人妻1区二区| 国产精品自产拍在线观看55亚洲| 欧美一级毛片孕妇| 日韩精品中文字幕看吧| 国产精品,欧美在线| 免费看十八禁软件| 人人澡人人妻人| 亚洲精品国产区一区二| 听说在线观看完整版免费高清| 日韩一卡2卡3卡4卡2021年| 韩国av一区二区三区四区| 999久久久精品免费观看国产| 天堂动漫精品| 成年免费大片在线观看| 老司机深夜福利视频在线观看| 91成人精品电影| 久久中文字幕一级| √禁漫天堂资源中文www| 亚洲欧美精品综合一区二区三区| 中文字幕人妻丝袜一区二区| 久久久精品国产亚洲av高清涩受| www日本在线高清视频| 久久精品夜夜夜夜夜久久蜜豆 | 欧美黑人精品巨大| 国产精品爽爽va在线观看网站 | svipshipincom国产片| 美女国产高潮福利片在线看| 久久精品aⅴ一区二区三区四区| 夜夜爽天天搞| 精品一区二区三区av网在线观看| www.自偷自拍.com| 亚洲七黄色美女视频| 精品久久久久久成人av| 亚洲人成网站在线播放欧美日韩| 免费在线观看完整版高清| 国产午夜精品久久久久久| 午夜激情福利司机影院| 激情在线观看视频在线高清| 亚洲成人国产一区在线观看| 在线观看免费午夜福利视频| 女性被躁到高潮视频| 男女做爰动态图高潮gif福利片| 欧美成人一区二区免费高清观看 | 亚洲avbb在线观看| 国产成人欧美| 制服丝袜大香蕉在线| 日韩欧美在线二视频| 一区二区三区激情视频| 久久中文字幕人妻熟女| 国产aⅴ精品一区二区三区波| 亚洲中文字幕一区二区三区有码在线看 | 精品一区二区三区视频在线观看免费| 91成年电影在线观看| 可以免费在线观看a视频的电影网站| 国产精品99久久99久久久不卡| 亚洲精品色激情综合| 麻豆国产av国片精品| 欧美日本视频| 男女床上黄色一级片免费看| 精品午夜福利视频在线观看一区| 999精品在线视频| 嫩草影院精品99| 精品午夜福利视频在线观看一区| 久久久久久人人人人人| 亚洲黑人精品在线| 一级毛片高清免费大全| 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 91九色精品人成在线观看| 亚洲片人在线观看| 国产av又大| 亚洲国产欧美网| 日本 av在线| 国产av又大| av电影中文网址| 丰满的人妻完整版| 两人在一起打扑克的视频| 国产精品一区二区精品视频观看| 琪琪午夜伦伦电影理论片6080| e午夜精品久久久久久久| 老司机在亚洲福利影院| 成在线人永久免费视频| 亚洲人成网站高清观看| 日韩精品免费视频一区二区三区| 午夜激情福利司机影院| 亚洲电影在线观看av| 国产精品爽爽va在线观看网站 | 久久国产乱子伦精品免费另类| 色精品久久人妻99蜜桃| 亚洲午夜理论影院| 欧美黑人欧美精品刺激| 午夜福利成人在线免费观看| 亚洲一区二区三区色噜噜| a级毛片a级免费在线| 欧美一级毛片孕妇| 18禁观看日本| 在线观看免费日韩欧美大片| 亚洲最大成人中文| 长腿黑丝高跟| 日韩欧美在线二视频| ponron亚洲| cao死你这个sao货| 亚洲 欧美一区二区三区| 亚洲第一av免费看| 日本一区二区免费在线视频| 婷婷精品国产亚洲av| 国产精品日韩av在线免费观看| 午夜福利成人在线免费观看| 亚洲国产欧美日韩在线播放| 18禁美女被吸乳视频| 国产aⅴ精品一区二区三区波| 99国产极品粉嫩在线观看| 欧美成人一区二区免费高清观看 | 男女那种视频在线观看| 美女午夜性视频免费| 亚洲色图 男人天堂 中文字幕| 成人永久免费在线观看视频| 国产亚洲欧美在线一区二区| 国产精品 欧美亚洲| 久久久久久九九精品二区国产 | 午夜久久久在线观看| 美女国产高潮福利片在线看| 俄罗斯特黄特色一大片| 可以在线观看毛片的网站| 一级a爱片免费观看的视频| 国产av又大| 久久久久国产精品人妻aⅴ院| 欧美黑人欧美精品刺激| 亚洲第一电影网av| 一级a爱片免费观看的视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲成av人片免费观看| 性欧美人与动物交配| 色婷婷久久久亚洲欧美| 一a级毛片在线观看| 亚洲aⅴ乱码一区二区在线播放 | 国产亚洲欧美精品永久| 美女国产高潮福利片在线看| 日本 av在线| 日韩大码丰满熟妇| 黄片播放在线免费| 999精品在线视频| 757午夜福利合集在线观看| 夜夜看夜夜爽夜夜摸| 亚洲av美国av| 99精品久久久久人妻精品| 又黄又粗又硬又大视频| 免费高清视频大片| www.熟女人妻精品国产| 久久国产精品男人的天堂亚洲| 欧美日韩福利视频一区二区| 久久国产精品影院| 狠狠狠狠99中文字幕| 国产精品电影一区二区三区| 欧美激情 高清一区二区三区| 亚洲成国产人片在线观看| 悠悠久久av| 亚洲黑人精品在线| 99国产精品一区二区三区| 国产av不卡久久| svipshipincom国产片| 最好的美女福利视频网| 亚洲人成电影免费在线| 色哟哟哟哟哟哟| 免费人成视频x8x8入口观看| 无限看片的www在线观看| 亚洲精品美女久久久久99蜜臀| 国产亚洲欧美精品永久| 91成人精品电影| 一级作爱视频免费观看| 亚洲一区二区三区色噜噜| 久久亚洲真实| 岛国在线观看网站| 男女之事视频高清在线观看| 国产亚洲精品一区二区www| 丰满的人妻完整版| 亚洲国产精品久久男人天堂| 一a级毛片在线观看| 日本撒尿小便嘘嘘汇集6| 亚洲中文字幕一区二区三区有码在线看 | 2021天堂中文幕一二区在线观 | 美女免费视频网站| av片东京热男人的天堂| www.自偷自拍.com| 国产极品粉嫩免费观看在线| 亚洲国产欧美日韩在线播放| 桃红色精品国产亚洲av| 亚洲欧美精品综合一区二区三区| 免费搜索国产男女视频| 日韩精品中文字幕看吧| 欧美成人性av电影在线观看| 无限看片的www在线观看| 国产亚洲欧美在线一区二区| 日日爽夜夜爽网站| 男女床上黄色一级片免费看| 欧美久久黑人一区二区| 欧美午夜高清在线| 精品日产1卡2卡| 亚洲av成人不卡在线观看播放网| 欧美日韩福利视频一区二区| 一进一出抽搐gif免费好疼| 黄片大片在线免费观看| 少妇 在线观看| 国产伦人伦偷精品视频| 一本综合久久免费| 丝袜美腿诱惑在线| 香蕉久久夜色| 国产精品精品国产色婷婷| 亚洲 欧美一区二区三区| 国产精品免费一区二区三区在线| 亚洲欧洲精品一区二区精品久久久| 精华霜和精华液先用哪个| xxxwww97欧美| 日本五十路高清| 精品熟女少妇八av免费久了| 欧美一区二区精品小视频在线| 亚洲五月婷婷丁香| 动漫黄色视频在线观看| 久久九九热精品免费| 精品福利观看| 母亲3免费完整高清在线观看| 欧美乱色亚洲激情| 丝袜美腿诱惑在线| 亚洲专区字幕在线| 免费看美女性在线毛片视频| 亚洲专区字幕在线| 性色av乱码一区二区三区2| av在线天堂中文字幕| 久久久久精品国产欧美久久久| 国产不卡一卡二| 首页视频小说图片口味搜索| 色在线成人网| 免费看日本二区| 午夜福利在线在线| aaaaa片日本免费| ponron亚洲| 两性夫妻黄色片| 手机成人av网站| 满18在线观看网站| 国产精品国产高清国产av| 精品国产乱码久久久久久男人| 久久久久久久久中文| 99在线人妻在线中文字幕| 精品第一国产精品| 亚洲国产欧洲综合997久久, | 国产三级黄色录像| 黄色毛片三级朝国网站| 久久久国产精品麻豆| 中文字幕av电影在线播放| 久久精品影院6| 夜夜躁狠狠躁天天躁| 欧美激情 高清一区二区三区| www.熟女人妻精品国产| 亚洲成人国产一区在线观看| 久久久久久久久免费视频了| 一二三四社区在线视频社区8| 91老司机精品| 国产成年人精品一区二区| 欧美久久黑人一区二区| 老汉色av国产亚洲站长工具| 亚洲 欧美一区二区三区| 久久精品影院6| 高潮久久久久久久久久久不卡| 岛国视频午夜一区免费看| 最新美女视频免费是黄的| 18禁国产床啪视频网站| 国产精品一区二区免费欧美| 亚洲精品美女久久av网站| 黄色a级毛片大全视频| 亚洲成av片中文字幕在线观看| 久久久久久久久免费视频了| 日韩三级视频一区二区三区| 久久久久久久久中文| 淫秽高清视频在线观看| 搡老熟女国产l中国老女人| 亚洲成人久久爱视频| 在线观看免费视频日本深夜| 国产色视频综合| 日韩中文字幕欧美一区二区| 久久午夜亚洲精品久久| 国产精品综合久久久久久久免费| 99在线视频只有这里精品首页| 黄色 视频免费看| 一区二区三区精品91| 国产视频一区二区在线看| 又大又爽又粗| 999精品在线视频| 男人操女人黄网站| 老熟妇仑乱视频hdxx| 在线观看日韩欧美| 嫩草影视91久久| 免费在线观看成人毛片| 国产亚洲精品第一综合不卡| 欧美成人免费av一区二区三区| 国产一区在线观看成人免费| 国产精品 欧美亚洲| 他把我摸到了高潮在线观看| 黄色a级毛片大全视频| 长腿黑丝高跟| 亚洲av片天天在线观看| 亚洲成人精品中文字幕电影| 色av中文字幕| 午夜免费观看网址| 久久青草综合色| 美女高潮到喷水免费观看| 精品高清国产在线一区| 成人亚洲精品一区在线观看| 又黄又粗又硬又大视频| 一个人观看的视频www高清免费观看 | 午夜福利欧美成人| 美女 人体艺术 gogo| 黑人巨大精品欧美一区二区mp4| 性色av乱码一区二区三区2| 国产99久久九九免费精品| 成人三级黄色视频| 很黄的视频免费| 一区二区三区激情视频| 狠狠狠狠99中文字幕| 亚洲国产欧美网| 美女高潮喷水抽搐中文字幕| 特大巨黑吊av在线直播 | 午夜激情福利司机影院| 亚洲国产精品999在线| 亚洲avbb在线观看| 亚洲,欧美精品.| 中文字幕久久专区| 亚洲性夜色夜夜综合| 免费在线观看视频国产中文字幕亚洲| 人成视频在线观看免费观看| 国产精品亚洲一级av第二区| 99精品久久久久人妻精品| 此物有八面人人有两片| 中出人妻视频一区二区| 午夜福利一区二区在线看| 日韩欧美在线二视频| 欧美激情 高清一区二区三区| 久久伊人香网站| 久久狼人影院| 欧美黄色淫秽网站| 亚洲成人久久性| 国产成人影院久久av| 又黄又粗又硬又大视频| 亚洲熟妇熟女久久| svipshipincom国产片| 色哟哟哟哟哟哟| 一本久久中文字幕| 欧美不卡视频在线免费观看 | 日韩高清综合在线| 亚洲精品国产一区二区精华液| 亚洲精品一区av在线观看| 一本一本综合久久| 日本成人三级电影网站| 夜夜爽天天搞| 女性被躁到高潮视频| 最近最新中文字幕大全免费视频| 我的亚洲天堂| 国产精品免费一区二区三区在线| 熟女少妇亚洲综合色aaa.| 亚洲专区字幕在线| 成人av一区二区三区在线看| 亚洲一码二码三码区别大吗| 日日摸夜夜添夜夜添小说| 国内揄拍国产精品人妻在线 | 亚洲男人的天堂狠狠| 无遮挡黄片免费观看| 久久香蕉激情| 国产99白浆流出| 一级片免费观看大全| 久久久久国产精品人妻aⅴ院| 深夜精品福利| 日韩欧美国产在线观看| 丁香欧美五月| 精品日产1卡2卡| 久久久国产欧美日韩av| 国产亚洲欧美98| 亚洲 欧美 日韩 在线 免费| 免费看十八禁软件| 一边摸一边做爽爽视频免费| 听说在线观看完整版免费高清| xxx96com| 免费高清视频大片| 99精品在免费线老司机午夜| 久久久水蜜桃国产精品网| 中文在线观看免费www的网站 | 天天添夜夜摸| 国产99久久九九免费精品| 欧美激情 高清一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 欧美在线黄色| 国产v大片淫在线免费观看| 亚洲国产精品sss在线观看| 三级毛片av免费| 天天躁夜夜躁狠狠躁躁| 欧美色视频一区免费| 熟女电影av网| 久久天躁狠狠躁夜夜2o2o| 欧美国产精品va在线观看不卡| 久久久精品欧美日韩精品| 看片在线看免费视频| 久久久水蜜桃国产精品网| 免费一级毛片在线播放高清视频| 国产成年人精品一区二区| 黄片播放在线免费| 两个人免费观看高清视频| 精品日产1卡2卡| 久久人人精品亚洲av| 久9热在线精品视频| 999精品在线视频| 麻豆国产av国片精品| 国产成人欧美| 在线观看日韩欧美| 久久久久久人人人人人| 中文字幕人妻熟女乱码| 亚洲 欧美 日韩 在线 免费| 国产成人影院久久av| 亚洲成人免费电影在线观看| www.www免费av| 亚洲中文字幕日韩| 一区二区三区精品91| 亚洲精品在线美女| 国产av又大| 一区二区三区高清视频在线| 国产欧美日韩一区二区精品| 一区二区三区高清视频在线| 欧美日韩福利视频一区二区| 韩国精品一区二区三区| x7x7x7水蜜桃| 嫩草影视91久久| 免费在线观看亚洲国产| 最好的美女福利视频网| 可以免费在线观看a视频的电影网站| 亚洲国产精品sss在线观看| 夜夜夜夜夜久久久久| 2021天堂中文幕一二区在线观 | www.自偷自拍.com| 亚洲人成伊人成综合网2020| 最好的美女福利视频网| 国产成人欧美在线观看| 18禁黄网站禁片免费观看直播| 欧美另类亚洲清纯唯美| 男女做爰动态图高潮gif福利片| 正在播放国产对白刺激| 免费在线观看完整版高清| 怎么达到女性高潮| 国产三级黄色录像| 亚洲中文日韩欧美视频| 香蕉丝袜av| 日韩精品中文字幕看吧| 精品久久久久久久末码| 亚洲中文字幕日韩| 最近在线观看免费完整版| 日韩免费av在线播放| 美女 人体艺术 gogo| 村上凉子中文字幕在线| 特大巨黑吊av在线直播 | 久久伊人香网站| 亚洲第一青青草原| 亚洲欧美日韩高清在线视频| av有码第一页| 久久久久久久久中文| 久久九九热精品免费| 一级黄色大片毛片| 国产色视频综合| 久久久久久九九精品二区国产 | 熟妇人妻久久中文字幕3abv| 欧美色视频一区免费| 男女下面进入的视频免费午夜 | 午夜免费观看网址| 黄片播放在线免费| av有码第一页| 欧美日本视频| 国产高清有码在线观看视频 | 伦理电影免费视频| 亚洲av成人av| 亚洲精品av麻豆狂野| 久久久久久国产a免费观看| 亚洲午夜理论影院| 亚洲aⅴ乱码一区二区在线播放 | 久久久精品欧美日韩精品| 亚洲中文字幕日韩| 99久久综合精品五月天人人| 国产三级黄色录像| 亚洲自偷自拍图片 自拍| 日韩av在线大香蕉| 女同久久另类99精品国产91| 欧美zozozo另类| 少妇熟女aⅴ在线视频| 免费看十八禁软件| 色哟哟哟哟哟哟| 99久久国产精品久久久| 国产成人系列免费观看| 欧美另类亚洲清纯唯美| 99热只有精品国产| 免费在线观看日本一区| 久久久久久国产a免费观看| 好看av亚洲va欧美ⅴa在| 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片| 一级作爱视频免费观看| 2021天堂中文幕一二区在线观 | 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 欧美性猛交╳xxx乱大交人| 夜夜看夜夜爽夜夜摸| 色av中文字幕| 一个人免费在线观看的高清视频| 久久精品91蜜桃| 久久草成人影院| 视频区欧美日本亚洲| 亚洲 欧美 日韩 在线 免费| 欧美日韩亚洲国产一区二区在线观看| a级毛片a级免费在线| 妹子高潮喷水视频| 欧美中文综合在线视频| 黄色视频,在线免费观看| 12—13女人毛片做爰片一|