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

    呼倫湖流域磷形態(tài)分布及顆粒態(tài)磷的溯源

    2023-09-23 12:14:38蘆曉峰郭軼男王國曦
    中國環(huán)境科學 2023年9期
    關鍵詞:呼倫湖源地泥沙

    蘆曉峰,郭軼男,,王國曦,姜 霞,王 坤*

    呼倫湖流域磷形態(tài)分布及顆粒態(tài)磷的溯源

    蘆曉峰1,郭軼男1,2,王國曦2,姜 霞2,王 坤2*

    (1.沈陽農(nóng)業(yè)大學水利學院,遼寧 沈陽 110866;2.中國環(huán)境科學研究院,湖泊水污染治理與生態(tài)修復技術國家工程實驗室,國家環(huán)境保護湖泊污染控制重點實驗室,北京 100012)

    選擇呼倫湖研究區(qū)23個點位,調(diào)查呼倫湖流域水體中不同形態(tài)磷的濃度,利用復合指紋識別技術解析沉積物顆粒態(tài)磷(PP)來源及不同來源對PP相對貢獻.結果表明:呼倫湖流域水體總磷(TP)濃度范圍為0.04~0.33mg/L,平均值為0.14mg/L,水體中顆粒態(tài)磷(PP)占TP比例范圍為45.79%~92.78%.呼倫湖沉積物PP來源主要為退化草地(34.1%~100.0%)、其次為草地(0%~50.7%)和農(nóng)田(0%~15.2%);烏爾遜河沉積物PP來源主要為退化草地(49.4%~57.1%)和草地(42.9%~50.6%);新開河沉積物PP來源是退化草地(100%);引河濟湖(人工河) 沉積物PP來源主要是退化草地(83.6%).本研究為呼倫湖流域磷污染的治理及湖泊水環(huán)境的改善提供了理論依據(jù).

    呼倫湖;磷;溯源;復合指紋法

    磷是生物生長所需的關鍵營養(yǎng)元素,也是造成湖泊富營養(yǎng)化和有害藻華爆發(fā)的主要營養(yǎng)因子[1].在自然水體中,磷被分為溶解性磷(DTP)和顆粒態(tài)磷(PP),其中顆粒態(tài)磷(PP)是水體總磷(TP)的主要賦存形式,在我國近海水體及淺水湖泊中占總磷的比例可達50%以上[2].在營養(yǎng)物質(zhì)的遷移轉(zhuǎn)化過程中,部分DTP被泥沙吸附形成PP沉積在湖泊底泥中;當存在外界物理化學擾動或者在生物化學作用下,沉積物中的磷重新被釋放進入水體[3],成為磷素向水體輸出的主要途徑[4-5].研究表明,在切斷外源河流磷輸入后,湖泊底泥中的PP成為水體磷的主要來源[6].

    復合指紋圖譜技術是研究泥沙來源的有效手段,根據(jù)侵蝕地泥沙組成差異篩選出判別能力最強的指紋識別因子組合,利用復合指紋模型建立流域沉積物與潛在源地之間的“指紋”連接,定量描述各潛在源地對流域沉積物的貢獻[7].泥沙作為PP輸出的載體,通過復合指紋識別技術可以在一定程度上建立PP的輸出與來源之間的對應關系,從而探討不同源地對PP的貢獻.已有研究[8]運用復合指紋技術篩選出最佳指紋組合,分析出三峽庫區(qū)汝溪河支流消落帶泥沙中PP的主要來源是長江干流懸移質(zhì)和汝溪河上游懸移質(zhì).

    呼倫湖位于內(nèi)蒙古呼倫貝爾草原腹地,是我國內(nèi)陸高緯度地帶最大的草原型湖泊,對維系呼倫貝爾大草原的生物多樣性及豐富動植物資源起著重要作用[9].但受氣候條件的變化及周邊人類活動的影響,呼倫湖出現(xiàn)面積減小、水位下降、水質(zhì)惡化、水體富營養(yǎng)化等一系列生態(tài)環(huán)境問題[10].呼倫湖水體中的磷常年處于GB3838-2007《地表水環(huán)境質(zhì)量標準》中的V類-劣V類水平,是呼倫湖水質(zhì)不達標的關鍵限制性因子[11].

    本文以呼倫湖及其主要入湖河流作為研究對象,探明流域土壤、湖泊(河流)沉積物及水體磷元素空間分布.基于現(xiàn)場調(diào)查和文獻分析確定林地、農(nóng)田、草地和退化草地為呼倫湖流域沉積物顆粒態(tài)磷的來源,利用復合指紋識別技術估算不同土地利用類型土壤的相對貢獻,揭示了呼倫湖流域沉積物顆粒態(tài)磷的泥沙來源,為呼倫湖流域的綜合治理與生態(tài)恢復提供科學依據(jù).

    1 材料與方法

    1.1 研究區(qū)概況

    呼倫湖(48°33′~49°20′N,116°58′~117°48′E)位于內(nèi)蒙古呼倫貝爾市新巴爾虎右旗、新巴爾虎左旗和滿洲里市扎賚諾爾區(qū)之間,是內(nèi)蒙古最大的湖泊,也是中國第四大淡水湖.呼倫湖湖面面積2086.48km2,平均水深5.33m,湖長93km,蓄水量達120億m3,地處溫帶大陸性氣候地區(qū),春季炎熱短促,冬季嚴寒漫長,湖水于11月上旬封冰并于次年5月初解凍,冰封期170~180d,是我國封凍期較長的湖泊[12].呼倫湖的湖水主要來源于海拉爾河、克魯倫河,以及從東部把它與貝爾湖連接起來的烏爾遜河,東北部的新開河是一條吞吐性河流,當海拉爾河水大時,河水順著新開河匯入呼倫湖,當呼倫湖水量大時,水又順此河流入額爾古納濕地,水少時則成為內(nèi)陸湖.呼倫湖的引河濟湖(人工河)的運行加強了呼倫湖與海拉爾河的水力連通,目前也是呼倫湖水資源補給的重要通道之一[13].

    1.2 樣品采集

    在研究區(qū)域呼倫湖、克魯倫河、烏爾遜河、新開河和引河濟湖人工河共設置23個采樣點.克魯倫河布設K0、K1、K2、K3、K4、K5共6個采樣點.采樣點主要土地利用類型為草地,有著明顯牧業(yè)活動,K0為克魯倫河入湖河口;烏爾遜河布設W0、W1、W2、W3、W4、W5共6個點,主要土地利用類型為草地,部分河流沿岸區(qū)域有一定數(shù)量的防護林帶,有著明顯牧業(yè)活動,W0為烏爾遜河入湖河口;新開河布設 X0、X1、X2、X3、X4共5個點,該河流域人類活動較為明顯,附近存在居民區(qū)和農(nóng)田等,X0為新開河入湖河口;引河濟湖布設U0、U1、U2 共3個點,其主要土地利用類型為退化草地,沙化嚴重,U0為引河濟湖入湖河口;呼倫湖布設H1、H2、H3共3個點.

    本次樣品采集于2020年11月下旬~12月上旬進行,采樣點位置如圖1所示.取0.3m水面下1.5L水樣裝入塑料瓶中,并用抓泥斗采集底泥,同時采岸邊土壤樣品(林地3個、農(nóng)田3個、草地9個、退化草地8個,共23個).所采集樣品及時運回實驗室,全部保存于4℃冰箱中,分析待用.

    圖1 研究區(qū)域位置及采樣點分布

    審圖號GS(2019)1822

    1.3 研究方法

    1.3.1 水體磷測定方法 水體的總磷(TP)濃度使用采集的原水樣進行測定,測定溶解性總磷(DTP)、無機磷(IP)濃度需將原水經(jīng)0.45μm的混合纖維濾膜過濾后測定.TP和DTP均采用過硫酸鉀消解鉬銻抗分光光度法測定,對于濾后液不進行過硫酸鉀氧化分解而直接進行鉬銻抗比色法,得到IP濃度.DOP濃度利用DTP和IP差減法獲得,顆粒態(tài)磷(PP)濃度利用TP和DTP差減法獲得[14].

    1.3.2 沉積物與土壤基本理化性質(zhì)和重金屬測定 沉積物與土壤樣品置于冷凍干燥機冷干2~3d,凍干后壓碎、翻動,剔除石塊、植物根莖等雜物,用木棒研壓后通過60目土壤篩.土壤和沉積物中的總碳(TC)和總氮(TN)采用元素分析儀(Elemental Analyzer, vario EL Cube)測定;重鉻酸鉀氧化-分光光度法(HJ 615-2011)測定土壤和沉積物中的有機碳(OC);H2O2-HNO3-HF消解土壤和沉積物后,使用電感耦合等離子體原子發(fā)射光譜法(ICP-OES, PerkinElmer OPTIMA8300)測定P、K、Na、Ca、Mg、Cu、Pb、Zn、Mn等金屬元素.

    1.3.3 復合指紋分析法 (1)復合指紋因子篩選利用Kruskal-Wallis H檢驗篩選可用于研究區(qū)域泥沙來源的指紋識別因子.這種方法的零假設是所有泥沙來自于同一源地,將Hcacl值與卡方檢驗結果(Hcritical)相比,當Hcacl>Hcritical時,則零假設被拒絕,表明該因子可以區(qū)別泥沙來源,若Hcad< Hcritical時,則不能拒絕零假設,不符合指紋因子的篩選條件.

    式中:s為源地秩和;n為源地樣本個數(shù);為所有源地樣本數(shù)量總和.

    基于多元判別分析找到最佳復合指紋組合[15].其主要判別依據(jù)為Wilks'Lambda值,每步以Wilks'Lambda的概計量最小的進入判別函數(shù),λ如下式可得:

    式中:SSerror為組內(nèi)離差平方和;SSstreat為組間離差平方和.多元判別分析復合指紋因子判別率大于80%的可以作為最佳復合指紋因子組合.

    (2)源地泥沙貢獻度分析利用多元混合模型,建立不同源地土壤與沉積泥沙中最佳指紋因子的定量關系,得到各泥沙源地泥沙貢獻百分比[16].模型函數(shù)如下:

    將得出貢獻率的泥沙樣品做擬合優(yōu)度(GOF)檢驗,一般認為當GOF大于0.8時,最優(yōu)化混合模型計算結果可信[17].公式如下:

    式(3)和(4)中:es為殘差平方和;ssi為研究對象沉積泥沙中指紋識別因子的濃度;P為泥沙源地的泥沙貢獻百分比;C為泥沙源地中指紋識別因子的平均濃度;為泥沙源地數(shù)量;為指紋識別因子的數(shù)量.

    (3)源地顆粒態(tài)磷(PP)貢獻度計算查閱文獻和水文資料確定不同研究區(qū)土壤顆粒態(tài)磷(PP)輸出系數(shù).再結合之前監(jiān)測和處理的數(shù)據(jù),計算得出不同源地土壤顆粒態(tài)磷(PP)貢獻比,計算公式如下[18]:

    式中:L為泥沙源地土壤PP貢獻率;R為泥沙源地土壤PP輸出系數(shù);為泥沙源地土壤TP平均濃度;P為用地類型的泥沙貢獻率;為泥沙源地個數(shù).

    1.4 數(shù)據(jù)分析

    本文采用Excel 2010和SPSS 16.0進行數(shù)據(jù)統(tǒng)計,采樣點空間分布圖利用Arcgis 10.6進行繪制和提取,其它分布圖利用Origin 2019b進行繪制.

    2 結果與討論

    2.1 呼倫湖流域土地利用類型及土壤基本理化性質(zhì)

    如表1,土壤TP濃度與土地利用方式有關,4種土地利用類型土壤中TP濃度平均值由大到小依次為農(nóng)田(969mg/kg)、草地(656mg/kg)、林地(552mg/kg)、退化草地(444mg/kg).由于農(nóng)業(yè)生產(chǎn)過程中人為活動向土壤中輸入磷導致農(nóng)田土壤TP濃度最高,而退化草地土地沙化嚴重,對營養(yǎng)元素儲存能力較弱,容易水土流失,導致土壤TP濃度最低0.呼倫湖4種土地利用類型土壤中TN、TC、OC濃度平均值由大到小均依次為草地、林地、農(nóng)田、退化草地.草地和林地TC(35.5和24.5mg/kg)、OC(29.6和25.6mg/kg)的平均濃度顯著高于農(nóng)田和退化草地TC(22.5和18.1mg/kg)、OC(22.4和19.9mg/kg)的平均濃度;草地土壤TN濃度平均值最高,達2687mg/kg,而退化草地土壤TN濃度平均值最低,僅為1469mg/kg.土壤TC、OC以及TN的主要來源是植物地上部分的凋落物及地下根系,隨著人為開發(fā)利用和過度畜牧,土壤中的有機質(zhì)濃度逐漸降低,且導致土壤TN流失[19-20].

    表1 呼倫湖流域土地利用類型及土壤基本理化性質(zhì)

    2.2 呼倫湖流域沉積物基本理化性質(zhì)

    呼倫湖及主要出入湖河流沉積物基本理化性質(zhì)分布如表2.呼倫湖湖區(qū)及四條河流采樣點沉積物TP濃度平均值由大到小依次為呼倫湖、新開河、克魯倫河、烏爾遜河、引河.沉積物TN濃度平均值由大到小依次為新開河、烏爾遜河、克魯倫河、引河、呼倫湖;沉積物TC濃度平均值由大到小依次為新開河、克魯倫河、烏爾遜河、呼倫湖、引河;沉積物OC濃度平均值由大到小依次為新開河、克魯倫河、烏爾遜河、引河、呼倫湖.新開河作為吞吐性河流,在呼倫湖收縮期時湖水水位下降,向呼倫湖補給供水,新開河由北向南流向呼倫湖,接納了流域內(nèi)的生活污水和牲畜糞便,導致河流沉積物碳、氮濃度偏高.烏爾遜河和克魯倫河沿岸有明顯的牧業(yè)活動,沉積物碳、氮濃度也會相對較高.

    表2 呼倫湖流域沉積物基本理化性質(zhì)

    2.3 呼倫湖流域水體磷形態(tài)分布

    如圖2,呼倫湖湖區(qū)及四條河流采樣點水體總磷(WTP)濃度平均值由大到小依次為引河、克魯倫河、呼倫湖、烏爾遜河、新開河.呼倫湖各采樣點WIP、DOP、PP、WTP濃度分布如圖3.呼倫湖所有采樣點的WTP濃度的變化范圍為0.04~0.33mg/L,平均值為0.14mg/L,最大值出現(xiàn)在引河的U0點位,最小值出現(xiàn)在烏爾遜河的W4點位和新開河的X0點位,顯著超過地表水富營養(yǎng)化總磷標準(0.1mg/L)的采樣點高達15個,說明呼倫湖流域有較高的水體富營養(yǎng)化風險.本研究結果與張博等[11]2019年冬季于呼倫湖流域所測得的WTP濃度(0.09~0.36mg/L)結果類似,但明顯低于于海峰等[9]于2011~2020年的測試值(0.13~0.60mg/L),主要原因可能是由于呼倫湖水質(zhì)季節(jié)性差異顯著,TP濃度春夏高于秋冬,且于海峰等的觀測時間跨度長,流域自然與社會條件發(fā)生了變化.

    圖2 呼倫湖及主要入湖河流WTP濃度

    各采樣點的PP濃度范圍在0.02~0.24mg/L,最大值在呼倫湖湖區(qū)的W0點位,最小值在烏爾遜河的W4點位,PP在TP的占比范圍為45.79%~92.78%,平均值為72.59%,可以說明呼倫湖各采樣點中磷以顆粒態(tài)為主,溶解于水體中的磷占少部分.DOP的濃度范圍在0.003~0.3mg/L,最大值在呼倫湖的H1點位,最小值在克魯倫河的K0和新開河的X0點位.IP濃度在各個采樣點中變化很大,濃度范圍為0.004~ 0.164mg/L,IP濃度顯著超過地表水富營養(yǎng)化的標準(0.03mg/L)的點位有5個.各個采樣點中IP與DOP的比值范圍在0.03~8.5之間,表明呼倫湖水質(zhì)的磷指標變異很大,部分水體中磷的礦化程度高.

    2.4 呼倫湖流域土地利用類型對沉積物及水體磷的影響

    研究表明,呼倫湖流域沉積物磷主要通過周圍土壤顆粒物輸入[21].如圖4顯示,呼倫湖流域土壤TP與沉積物TP呈顯著的正相關關系,但整個流域無論是土壤TP還是沉積物TP與呼倫湖流域水體磷之間相關性均不顯著.以土壤和沉積物基本理化指標以及水體磷形態(tài)作為主成分分析因子,對23個采樣點進行聚類分析.結果顯示主成分因子在X軸上的占比為55.2%,在Y軸上的占比為20.7%.未受干擾土壤樣點區(qū)(草地+林地)和受干擾土壤樣點區(qū)(退化草地+農(nóng)田)在95%的置信區(qū)間內(nèi)能被顯著區(qū)分開.其中,未受干擾區(qū)域樣點的主要聚類成分因子為呼倫湖流域沉積物和土壤基本理化性質(zhì);而受干擾區(qū)域樣點的主要聚類成分因子為水體中不同磷形態(tài).相關性分析也表明,在呼倫湖流域未受干擾區(qū)域土壤TP、沉積物TP與WTP、DOP以及PP呈顯著負相關關系;而受干擾區(qū)土壤TP、沉積物TP與WTP呈顯著正相關關系.這表明土地退化是呼倫湖流域水體磷主要影響因素(圖6、圖7).部分研究顯示,呼倫湖流域土地沙化嚴重,水土流失嚴重導致大量營養(yǎng)鹽輸入水體是導致水體磷濃度偏高的主要原因[22][23].

    圖4 呼倫湖流域不同類型采樣點主成分分析

    圖5 退化草地土壤和未受干擾土壤有機碳和土壤總磷的擬合

    有機質(zhì)作為土壤重要理化指標,對沉積物和水體磷累積和形態(tài)有重要影響[24-25].在呼倫湖流域,未受干擾區(qū)域土壤有機質(zhì)與沉積物磷和有機質(zhì)呈顯著正相關關系,但和水體WTP、PP及DOP呈顯著負相關關系(圖7).表明未受干擾區(qū)土壤有機質(zhì)可能是沉積物有機質(zhì)和磷的主要來源,土壤有機質(zhì)的賦存有效減少了土壤磷向水體的輸移.已有研究證實,土壤有機質(zhì)能有效緩解土壤沙化,減少土壤磷流失,增加土壤磷賦存[26][27].通過對呼倫湖流域土壤有機碳(SOC)和土壤TP擬合發(fā)現(xiàn),在未受干擾區(qū)(草地+林地)隨著土壤有機碳減少,土壤磷流失速度逐漸加快;而在受干擾區(qū)(退化草地),SOC濃度與TP濃度呈線性擬合.主要原因可能是土壤有機質(zhì)中本身富含磷元素,同時能促進土壤團聚體的形成,有效貯存磷營養(yǎng)鹽;而當土壤有機質(zhì)流失時,土壤大團聚體分解朝沙化趨勢發(fā)展,貯存的磷營養(yǎng)鹽迅速流失,少量的磷元素主要通過有機質(zhì)貯存[28].如圖5顯示,兩條擬合線在SOC濃度為14.07mg/kg處相交,表明當呼倫湖流域SOC低于14.07mg/kg時,可以認為土壤退化,對磷營養(yǎng)鹽的貯存能力逐漸減弱.

    Fig 6 Heatmap of the correlation between soil basic properties, sediments basic properties and water phosphorus in the undisturbed area of Hulun Lake Basin

    (a)未受干擾區(qū)域 (b)受干擾區(qū)域

    2.5 呼倫湖沉積物磷來源解析

    2.5.1 最佳復合指紋因子組合的確定 本研究以Na、Al、K、Cr、Mn、Fe、Cu、Zn、Cd、Sn、Pb作為各類源地土壤待選指紋因子進行初步篩選(表S1;表S2);基于Kruskal-Wallis H[15]檢驗后,Na、Al、K、Cd的P值小于0.05,其通過檢驗,可作為初步篩選的指紋因子.經(jīng)多元判別分析,指紋因子累積判別率最高的組合為K和Cd,分別達68.2%和100.0%,因此作為識別泥沙來源的最佳指紋因子組合,并通過判別函數(shù)預測土壤樣品,全部判別正確.利用多元混合模型,建立不同源地土壤與沉積泥沙中最佳指紋因子的定量關系,得到各泥沙源地泥沙貢獻百分比[16],烏爾遜河采樣點W2、W3、W4,新開河采用點 X1、X3、X4,克魯倫河采樣點 K1、K3、K4、K5,引河采樣點U2及呼倫湖采樣點 H1、H3、W0、X0、U0的沉積物樣品源地泥沙貢獻擬合優(yōu)度(GOF)值小于0.8,未通過檢驗.

    2.5.2 呼倫湖流域不同源地土壤對沉積物磷的貢獻 對通過GOF檢驗的采樣點計算沉積物源地土壤貢獻率(表3).四種潛在源地土壤泥沙貢獻比從高到低依次為退化草地(73.2%)、草地(20.3%)、林地(3.4%)和農(nóng)田(3.1%).草地和退化草地為研究區(qū)域的主要用地類型,退化草地與其他土地用地相比極易產(chǎn)生風蝕和水蝕,從而對河流沉積物影響很大.烏爾遜河沉積物(點位W5、W1)的泥沙主要來源于草地和退化草地,草地貢獻度為34.7% ~ 42.1%,退化草地貢獻度為57.9% ~ 65.3%;呼倫湖沉積物(點位H2)泥沙來源主要是退化草地,貢獻度達100.0%,可能與當?shù)夭莸赝嘶瘜е碌乃亮魇в嘘P;新開河沉積物點位X0位于出河口,泥沙來源主要是退化草地(40.3%)和草地(42.4%),點位X2泥沙主要來源為退化草地,貢獻度為100%;引河沉積物(點位U1)沉積物泥沙主要來源于退化草地(75.8%)和農(nóng)田(18.7%).

    基于文獻和水文資料調(diào)研,呼倫湖流域林地、農(nóng)田、草地、退化草地土壤顆粒態(tài)磷輸出系數(shù)分別為54%、14%、68%及87%[29-30],對應土壤總磷平均濃度分別為552mg/kg、969mg/kg、656mg/kg和444mg/kg(表3);由公式(3)計算得出不同源地土壤對呼倫湖流域沉積物磷貢獻(表3).

    4種潛在源地土壤對沉積物磷貢獻比從高到低為退化草地(70.7%)、草地(24.7%)、農(nóng)田(1.5%)、林地(3.1%).烏爾遜河沉積物磷主要來源于草地和退化草地,草地貢獻度為42.9%~50.6%,退化草地貢獻度為49.4%~57.1%;新開河沉積物磷主要來源為退化草地,貢獻度為100%;引河沉積物磷主要來源為退化草地(83.6%)和農(nóng)田(9.2%),還有少量來源于林地(3.5%)和草地(3.8%);呼倫湖沉積物磷來源主要為退化草地,貢獻度為50.7%~100%.樊才睿等[31][32]研究表明,在相同降雨強度下,呼倫湖流域退化草地水土流失最為嚴重,流失的TP濃度要遠高于其他土地利用類型,且以PP為主,這與本文研究結果一致.

    總體上,退化草地對于河流(湖泊)沉積物以及沉積物磷的貢獻最高,分別為73.2%和70.7%;盡管農(nóng)田中土壤TP濃度最高,但由于用地類型總體占比較小,對于呼倫湖流域沉積物及沉積物磷的貢獻度最小.

    表3 呼倫湖流域不同源地土壤對沉積物及沉積物磷的貢獻

    3 結論

    3.1 呼倫湖水體總磷(WTP)濃度范圍為0.05~ 0.34mg/L,平均值為0.14mg/L,顯著超過地表水富營養(yǎng)化總磷標準(0.1mg/L),呼倫湖流域水體中的磷主要以顆粒態(tài)(PP)形式存在.呼倫湖流域水體總磷(WTP)濃度平均值由大到小依次為引河、克魯倫河、呼倫湖、烏爾遜河、新開河.

    3.2 呼倫湖流域農(nóng)田土壤總磷濃度最高,達969mg/kg,退化草地總磷濃度最低,為444mg/kg;流域沉積物總磷濃度平均值由大到小依次為呼倫湖、新開河、克魯倫河、烏爾遜河、引河.

    3.3 土壤有機質(zhì)流失導致的土地退化造成呼倫湖流域土壤磷流失,并影響水體磷濃度和形態(tài)分布,當有機碳濃度低于14.07mg/kg時,可以認為土壤退化,降低了對磷營養(yǎng)鹽的貯存能力.

    3.4 呼倫湖流域沉積物磷主要來源于周圍不同源地土壤輸入,整個流域沉積物磷和土壤磷呈顯著正相關關系.呼倫湖流域退化草地對于河流(湖泊)沉積物以及沉積物磷的貢獻率最高,分別為73.2%和70.7%,其次是草地對河流(湖泊)沉積物及沉積物磷的貢獻率為20.3%和24.7%;盡管農(nóng)田土壤磷濃度最高,但由于用地類型總體占比較小,對于呼倫湖流域沉積物磷的貢獻度最小,僅為1.5%.

    [1] 劉建康,黃祥飛.東湖生態(tài)學研究概況[J]. 環(huán)境科學, 1997. Liu J K, Huang X F. Ecological survey of East Lake [J]. Environmental Science, 1997.

    [2] Xu J, et al. Distribution, transfer process and influence factors of phosphorus at sediment-water inter face in the Huaihe River [J]. Journal of Hydrology, 2022.612.

    [3] 郭長城,王國祥,喻國華.天然泥沙對富營養(yǎng)化水體中磷的吸附特性研究[J]. 中國給水排水, 2006,(9):10-13. Guo C C, Wang G X, Yu G H. Study on adsorption of phosphorus in eutrophied water body by natural sediment [J]. China Water & Waste Water, 2006,(9):10-13.

    [4] Vanni M J, Renwick W H, Headworth J L, et al. Dissolved and particulate nutrient flux from three adjacent agricultural watersheds: Afive-yearstudy [J]. Biogeochemistry, 2001,54(1):85-114.

    [5] Foster I D L, Chapman A S, Hodgkinson R M, et al. Changing suspended sediment and particulate phosphorus loads and pathways in underdrained lowland agricultural catchments; Herefordshire and Worcestershire, U.K. [J]. Hydrobiologia, 2003,494(1-3):119-126.

    [6] 張毅敏,王 宇,楊 飛,等.太湖不同生態(tài)型湖區(qū)懸浮顆粒磷空間分布和降解速率[J]. 中國環(huán)境科學, 2016,36(7):2128-2138. Zhang Y M, Wang Y, Yang F, et al. The spatial distribution and degradation characteristic of phosphorus in suspended particulate matter among different ecological types in Taihu [J]. China Environmental Science, 2016,36(7):2128-2138.

    [7] 周子柯,王永平,滕昊蔚,等.復合指紋技術示蹤泥沙來源研究進展[J]. 泥沙研究, 2021,46(6):73-80.DOI:10.16239/j.cnki.0468-155x. 2021.06.011. Zhou Z K, Wang Y P, Teng H W, et al. Progresses of compositefingerprinting technique to trace sediment sources [J]. Journal of Sediment Research, 2021,46(06):73-80.DOI:10.16239/ j.cnki.0468-155x.2021.06.011.

    [8] 陳太麗,史忠林,王永艷,等.三峽水庫典型支流消落帶泥沙顆粒態(tài)磷復合指紋示蹤研究[J]. 農(nóng)業(yè)工程學報, 2019,35(20):118-124.Chen T L, Shi Z L, Wang Y Y, et al. Fingerprinting particulate phosphorus absorbed by sediments for riparian zone deposits in tributary of Three Gorges Reservoir [J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(20):118-124.

    [9] 于海峰,史小紅,孫 標,等.2011~2020年呼倫湖水質(zhì)及富營養(yǎng)化變化分析[J]. 干旱區(qū)研究, 2021,38(6):1534-1545.DOI:10.13866/ j.azr.2021.06.05. Yu H F, Shi X H, Sun B, et al. Analysis of water quality and eutrophication changes in Hulun Lake from 2011 to 2020 [J]. Arid ZoneResearch, 2021,38(6):1534-1545.DOI:10.13866/j.azr.2021.06.05.

    [10] 韓向紅,楊 持.呼倫湖自凈功能及其在區(qū)域環(huán)境保護中的作用分析[J]. 自然資源學報, 2002,17(6):684-690. Han X H, Yang C. An analysis of the self-purification function of Hulun Lake and its effect on regional environmental conservation [J]. Journal of Natural Resources, 2002,17(6):684-690.

    [11] 張 博,郭云艷,王書航,等.呼倫湖水體磷的時空演變及其影響因素[J]. 環(huán)境科學研究, 2021,34(4):824-830. Zhang B, Guo Y Y, Wang S H, et al. Spatial-temporal changes of phosphorus and its influential factors in Lake Hulun [J]. Research of Environmental Sciences, 2021,34(4):824-830.

    [12] 張志波.呼倫湖志 [M]. 呼和浩特:內(nèi)蒙古文化出版社, 1998. Zhang Z B. Records of Hulun Lake [M]. Hohhot: Inner Mongolia Culture Publishing House, 1998.

    [13] 王鵬飛,郭云艷,周 康,等.1961~2018年呼倫湖水面面積變化特征及其對氣候變化的響應[J]. 環(huán)境科學研究, 2021,34(4):792-800. DOI:10.13198/j.issn.1001-6929.2021.02.12.

    [14] Wang P F, Guo Y Y, Zhou K, et al. Variation of lake area of Hulun Lake during 1961~2018 and its response to climate change [J]. Research of Environmental Sciences, 2021,34(4):792-800.DOI: 10.13198/j.issn.1001-6929.2021.02.12.

    [15] Katsaounos C Z, Giokas D L, Leonardos I D, et al. Speciation of phosphorus fractionation in river sediments by explanatory data analysis [J]. Water Research, 2007,41(2):406-418.

    [16] 郝福星,黃炎和,林金石,等.指紋法研究花崗巖區(qū)典型崩崗小流域懸浮泥沙來源[J]. 水土保持學報, 2017,31(02):45-49+102.DOI: 10.13870/j.cnki.stbcxb.2017.02.009.Hao F X, Huang Y H, Lin J S, et al. Fingerprinting suspended sediment Sourcesina typical collapsing hill watershed in the granitic Region [J]. Journal of Soil and Water Conservation, 2017,31(2):45- 49+102.DOI:10.13870/j.cnki.stbcxb.2017.02.009.

    [17] Walling D E. Tracing suspended sediment sources in catchments and river systems [J]. Science of the Total Environment, 2005,344(1-3): 159-184.

    [18] Motha J A, Wallbrink P J, Hairsine P B, et al. Determining the sources of suspended sediment in a forested catchment in southeastern Australia [J]. Water Resources Research, 2003,39(3).

    [19] 許志信,李永強.草地退化與水土流失[C]//現(xiàn)代草業(yè)科學進展——中國國際草業(yè)發(fā)展大會暨中國草原學會第六屆代表大會論文集, 2002:343-346. Xu Z X, Li Y Q. Grassland degeneration and soil erosion [C]//Progress in modern prataculture science——Proceedings of the 6th congress of China international prataculture development conference and Chinese Grassland Society, 2002:343-346.

    [20] 李紹良,陳有君,關世英,等.土壤退化與草地退化關系的研究[J]. 干旱區(qū)資源與環(huán)境, 2002,16(1):92-95.DOI:10.3969/j.issn.1003- 7578.2002.01.016. Li S L, Chen Y J, Guan S Y, et al. Relationships between soil degradation and rangeland degradation [J]. Journal of Arid Land Resources and Environment, 2002,16(1):92-95.DOI:10.3969/j.issn. 1003-7578.2002.01.016.

    [21] Rossignol N, Bonis A, Jan-Bernard Bouzillé. Consequence of grazing pattern and vegetation structure on the spatial variations of net N mineralisation in a wet grassland [J]. Applied Soil Ecology, 2006, 31(1/2):62-72.

    [22] 高宏斌.基于高分辨率沉積指標重建呼倫湖近期沉積環(huán)境演變過程[D]. 呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學, 2017. Gao H B. The reconstruction of evolution processes of sedimentary environment lately in Lake Hulun based on high resolution sedimentary indicators [D]. Hohhot: Inner Mongolia Agricultural University, 2017.

    [23] 郭鴻鵬,朱靜雅,楊印生.農(nóng)業(yè)非點源污染防治技術的研究現(xiàn)狀及進展[J]. 農(nóng)業(yè)工程學報, 2008,(4):290-295. Guo H P, Zhu J Y, Yang Y S. Research status and development of technologies for controlling the agricultural non-point source pollution [J]. Transactions of the CSAE, 2008,24(4):290-295.

    [24] Nonpointpollution of surface waters with phosphorus and nitrogen [J]. Ecological Applications, 1998,8(3):559-568.

    [25] 金相燦.沉積物污染化學[M]. 北京:中國環(huán)境科學出版社, 1992. Jin X C. Sediment pollution chemistry [M]. Beijing: China Environmental Science Press, 1992.

    [26] Kastelan-Macan M, Petrovic M. The role of fulvic acids in phosphorus sorption and release from mineral particles [J]. Water Science & Technology, 1996,34(7/8):259-265.

    [27] 陶士鋒,徐曉峰,寇太記.土壤有機質(zhì)對有效磷及水提取磷濃度的影響[J]. 中國生態(tài)農(nóng)業(yè)學報, 2012,20(8):1054-1058. Tao S F, Xu X F, Kou T J. Effect of soil organic matter on the contents of available and water-extracted phosphorus [J]. Chinese Journal of Eco-Agriculture, 2012,20(8):1054-1058.

    [28] 楊春霞,王圣瑞,金相燦,等.輕組有機質(zhì)對太湖沉積物氮、磷礦化的影響[J]. 環(huán)境科學研究, 2009,22(09):1001-1007.DOI:10.13198/ j.res.2009.09.11.yangchx.003. Yang C X, Wang S R, Jin X C, et al. Effect of light fraction organic matter on mineralization of N and P in Taihu Lake sediments [J]. Research of Environmental Sciences, 2009,22(9):1001-1007.DOI: 10.13198/j.res.2009.09.11.yangchx.003.

    [29] 梁愛珍,張 延,陳學文,等.東北黑土區(qū)保護性耕作的發(fā)展現(xiàn)狀與成效研究[J]. 地理科學, 2022,42(08):1325-1335.DOI:10.13249/ j.cnki.sgs.2022.08.001. Liang A Z, Zhang Y, Chen X W, et al. Development and effects of conservation tillage in the black soil region of Northeast China [J]. Scientia Geographica Sinica, 2022,42(8):1325-1335.] doi:10.13249/j. cnki.sgs.2022.08.001.

    [30] 黃金良,洪華生,張珞平,等.基于GIS的九龍江流域農(nóng)業(yè)非點源氮磷負荷估算研究[J]. 農(nóng)業(yè)環(huán)境科學學報, 2004,(5):866-871. Huang J L, Hong H S, Zhang L P, et al. Nitrogen and phosphorus loading of agricultural non - point sources in Jiulong River watershed based on GIS [J]. Journal of Agro-Environment Science, 2004,(5): 866-871.

    [31] 張汪壽.平原區(qū)非點源污染特征及關鍵源區(qū)的識別-以武清區(qū)為例[D]. 首都師范大學, 2012. Zhang W S. The characteristic of non-point source pollution and identification of critical source areas for plain area-A case study in Wuqing District [D]. Capital Normal University, 2012.

    [32] 樊才睿,李暢游,孫 標,等.不同放牧制度對呼倫貝爾草原徑流中磷流失模擬研究[J]. 水土保持學報, 2017,31(01):17-23,29.DOI: 10.13870/j.cnki.stbcxb.2017.01.004. Fan C R, Li C Y, Sun B, et al. Study on effects so different grazing systems on phosphorus loss in surface runoff in Hulunbair grasslands [J]. Journal of Soil and Water Conservation, 2017,31(01):17-23,29. DOI:10.13870/j.cnki.stbcxb.2017.01.004.

    Distribution characteristics of various forms of phosphorus and the traceability of phosphorus in particulate form in Hulun Lake basin.

    LU Xiao-feng1, GUO Yi-nan1,2, WANG Guo-xi2, JIANG Xia2, WANG Kun2*

    (1.College of Water Conservancy, Shenyang Agricultural University, Shenyang 110866, China;2.National Engineering Laboratory for Lake Pollution Control and Ecological Restoration, State Environment Protection Key Laboratory for Lake Pollution Control, Chinese Research Academy of Environmental Sciences, Beijing 100012, China)., 2023,43(9):4810~4818

    In this paper, we selected 23 sites in the Hulun Lake area to analyze the concentration of different phosphorus forms in the sediments of Hulun Lake and its inlet and outlet rivers. The composite fingerprint identification technology was then used to explore the sources of particulate phosphorus (PP) and the relative contributions of different sources. Results showed that the total phosphorus (TP) concentration in Hulun Lake Basin ranged from 0.04 to 0.33mg/L, with an average of 0.14mg/L. The proportions of PP to TP was 45.79% ~ 92.78%. The main source of PP in Hulun Lake sediments was degraded grassland (34.1%~100.0%), followed by grassland (0%~50.7%) and farmland (0%~15.2%). Degraded grassland was also the main source of PP in the sediments of Urson River (49.4%~57.1%), Xinkai River (artificial river, 100%) and Yinhejihu River (artificial river, 83.6%). Grassland contributed an additional 42.9%~50.6% of PP in the Urson River sediments. This study provides a theoretical basis for the control of phosphorus pollution in the Hulun Lake Basin and the improvement of the lake water environment.

    Hulun Lake;phosphorus;traceability;composite fingerprinting method

    X524

    A

    1000-6923(2023)09-4810-09

    蘆曉峰(1981-),男,遼寧沈陽人,副教授,博士,主要從事水土環(huán)境與生態(tài)工程研究.發(fā)表論文40余篇.13654906041@163.com.

    蘆曉峰,郭軼男,王國曦,等.呼倫湖流域磷形態(tài)分布及顆粒態(tài)磷的溯源 [J]. 中國環(huán)境科學, 2023,43(9):4810-4818.

    Lu X F, Guo Y N, Wang G X, et al. Distribution characteristics of various forms of phosphorus and the traceability of phosphorus in particulate form in Hulun Lake basin [J]. China Environmental Science, 2023,43(9):4810-4818.

    2023-02-08

    國家重點研發(fā)計劃項目(2022YFC3204000)

    * 責任作者, 高級工程師, wangkun@craes.org.cn

    猜你喜歡
    呼倫湖源地泥沙
    泥沙做的父親
    坐在呼倫湖畔看星星
    遼河(2022年1期)2022-02-14 19:45:39
    坐在呼倫湖畔看星星
    遼河(2022年1期)2022-02-14 05:15:04
    新疆多泥沙河流水庫泥沙處理措施
    新巴爾虎寫意(組詩)
    北極光(2019年5期)2019-11-13 04:54:18
    移民與文化認同:土家族民歌《吳幺姑》探析
    北方文學(2019年21期)2019-08-18 15:23:25
    土壤團聚體對泥沙沉降速度的影響
    發(fā)源地
    氣候變化對渭河源地水文環(huán)境影響分析與探討
    不同種源地漆樹種子生物學特性研究
    青青草视频在线视频观看| 成年av动漫网址| 国产av码专区亚洲av| 永久免费av网站大全| 男女下面插进去视频免费观看| 日韩一本色道免费dvd| 国产精品不卡视频一区二区| 国产一级毛片在线| 爱豆传媒免费全集在线观看| 亚洲综合色惰| 亚洲国产欧美网| 亚洲国产色片| 国产午夜精品一二区理论片| 丝袜人妻中文字幕| 日本猛色少妇xxxxx猛交久久| 日韩免费高清中文字幕av| 老熟女久久久| videosex国产| 日本vs欧美在线观看视频| 久久精品久久精品一区二区三区| 女的被弄到高潮叫床怎么办| 亚洲三级黄色毛片| 纵有疾风起免费观看全集完整版| 欧美另类一区| 妹子高潮喷水视频| 女人久久www免费人成看片| av免费观看日本| 9色porny在线观看| 国产精品久久久久久精品古装| 少妇被粗大的猛进出69影院| 91午夜精品亚洲一区二区三区| 久久久国产一区二区| 色播在线永久视频| 日日啪夜夜爽| 国产成人精品久久久久久| 99热网站在线观看| 国产免费一区二区三区四区乱码| 精品亚洲成a人片在线观看| 看十八女毛片水多多多| 国产国语露脸激情在线看| 人人妻人人添人人爽欧美一区卜| 久热久热在线精品观看| 国产白丝娇喘喷水9色精品| 美女主播在线视频| 大片免费播放器 马上看| 久久久久久人妻| www.av在线官网国产| 五月伊人婷婷丁香| 99九九在线精品视频| 亚洲精品一区蜜桃| 午夜福利在线观看免费完整高清在| 色吧在线观看| 9热在线视频观看99| 精品视频人人做人人爽| 人成视频在线观看免费观看| 又粗又硬又长又爽又黄的视频| 国产精品久久久久成人av| 国产亚洲精品第一综合不卡| 亚洲成人av在线免费| 亚洲国产av新网站| av不卡在线播放| tube8黄色片| 欧美日韩一区二区视频在线观看视频在线| 色吧在线观看| 97精品久久久久久久久久精品| 久久这里有精品视频免费| 午夜福利网站1000一区二区三区| 男女高潮啪啪啪动态图| 如何舔出高潮| 人体艺术视频欧美日本| 一边亲一边摸免费视频| 黄网站色视频无遮挡免费观看| 另类亚洲欧美激情| 国产精品偷伦视频观看了| 国产成人av激情在线播放| 另类精品久久| 热99久久久久精品小说推荐| 亚洲成色77777| 永久网站在线| 精品人妻熟女毛片av久久网站| 国产爽快片一区二区三区| 99热全是精品| 伊人久久大香线蕉亚洲五| 精品人妻在线不人妻| 永久网站在线| 青春草亚洲视频在线观看| 日韩中文字幕欧美一区二区 | 国产亚洲一区二区精品| 成人免费观看视频高清| 啦啦啦在线观看免费高清www| 亚洲精品日本国产第一区| 成年女人在线观看亚洲视频| 亚洲精品日韩在线中文字幕| 91精品国产国语对白视频| 亚洲美女搞黄在线观看| 99精国产麻豆久久婷婷| 午夜久久久在线观看| 欧美日韩精品网址| 久久久久国产网址| 免费不卡的大黄色大毛片视频在线观看| 久久国内精品自在自线图片| 男女无遮挡免费网站观看| 麻豆av在线久日| 久久久久人妻精品一区果冻| 成人二区视频| 成人免费观看视频高清| 国产一区亚洲一区在线观看| 丰满饥渴人妻一区二区三| 久久av网站| 69精品国产乱码久久久| 亚洲av免费高清在线观看| 亚洲第一av免费看| 国产精品成人在线| 99九九在线精品视频| 在线观看一区二区三区激情| 91国产中文字幕| 91在线精品国自产拍蜜月| 男人舔女人的私密视频| 国产av码专区亚洲av| 亚洲内射少妇av| 一级a爱视频在线免费观看| 黄片小视频在线播放| 日本色播在线视频| 老司机亚洲免费影院| 一区二区日韩欧美中文字幕| 国产精品无大码| 国精品久久久久久国模美| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看| 色哟哟·www| 边亲边吃奶的免费视频| www.av在线官网国产| 午夜福利影视在线免费观看| 777久久人妻少妇嫩草av网站| 97人妻天天添夜夜摸| 晚上一个人看的免费电影| 少妇的逼水好多| 叶爱在线成人免费视频播放| 免费观看a级毛片全部| 一级毛片电影观看| 2021少妇久久久久久久久久久| 中文字幕亚洲精品专区| 久久久国产一区二区| 亚洲色图综合在线观看| 国语对白做爰xxxⅹ性视频网站| 一边摸一边做爽爽视频免费| 欧美97在线视频| 伦理电影大哥的女人| 欧美日韩视频精品一区| 2018国产大陆天天弄谢| 久久久久久人人人人人| 久久精品夜色国产| 久热这里只有精品99| 精品少妇黑人巨大在线播放| 亚洲av日韩在线播放| 日韩不卡一区二区三区视频在线| av天堂久久9| 成年人午夜在线观看视频| 久久人人97超碰香蕉20202| 欧美人与善性xxx| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品国产一区二区精华液| 人妻系列 视频| 精品酒店卫生间| 精品久久蜜臀av无| 99热全是精品| 成年美女黄网站色视频大全免费| 久久99热这里只频精品6学生| 日日啪夜夜爽| 汤姆久久久久久久影院中文字幕| 亚洲婷婷狠狠爱综合网| 国产淫语在线视频| 亚洲色图 男人天堂 中文字幕| 精品卡一卡二卡四卡免费| 久久国产精品男人的天堂亚洲| 欧美日韩成人在线一区二区| 国产精品久久久久成人av| 亚洲欧美成人精品一区二区| 成人国语在线视频| 18禁动态无遮挡网站| 国产成人精品久久久久久| 久久综合国产亚洲精品| 久久久国产欧美日韩av| 一级片免费观看大全| 日日爽夜夜爽网站| 考比视频在线观看| 69精品国产乱码久久久| 亚洲欧美精品综合一区二区三区 | 波多野结衣一区麻豆| 午夜日韩欧美国产| 国产毛片在线视频| 国产成人免费观看mmmm| 亚洲av综合色区一区| 美女脱内裤让男人舔精品视频| √禁漫天堂资源中文www| 18禁动态无遮挡网站| 久久精品国产亚洲av高清一级| 在线观看美女被高潮喷水网站| 久久久国产欧美日韩av| 少妇被粗大猛烈的视频| 秋霞在线观看毛片| 永久网站在线| 精品卡一卡二卡四卡免费| 少妇被粗大的猛进出69影院| 久久久久视频综合| 亚洲美女搞黄在线观看| 亚洲激情五月婷婷啪啪| 我要看黄色一级片免费的| 尾随美女入室| 2018国产大陆天天弄谢| 精品亚洲成国产av| 又大又黄又爽视频免费| 一边摸一边做爽爽视频免费| 免费观看av网站的网址| 久久99热这里只频精品6学生| 国产又爽黄色视频| 菩萨蛮人人尽说江南好唐韦庄| 欧美精品亚洲一区二区| 国产有黄有色有爽视频| 中文字幕av电影在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲第一区二区三区不卡| 亚洲欧美成人精品一区二区| 日韩中文字幕欧美一区二区 | 国产一区二区三区综合在线观看| 中文天堂在线官网| 免费不卡的大黄色大毛片视频在线观看| 亚洲婷婷狠狠爱综合网| 国产乱人偷精品视频| 亚洲av电影在线进入| 日韩中字成人| 人人澡人人妻人| 狠狠精品人妻久久久久久综合| 交换朋友夫妻互换小说| 亚洲,一卡二卡三卡| 国产精品香港三级国产av潘金莲 | 国产黄色免费在线视频| 免费日韩欧美在线观看| 久久久久久久久久久免费av| 国产综合精华液| 久久久久久久亚洲中文字幕| 久久精品久久久久久噜噜老黄| 97精品久久久久久久久久精品| 99国产综合亚洲精品| 国产视频首页在线观看| 亚洲三级黄色毛片| 免费不卡的大黄色大毛片视频在线观看| 国产精品久久久久久精品古装| 纵有疾风起免费观看全集完整版| 中文字幕精品免费在线观看视频| 国产探花极品一区二区| 丰满迷人的少妇在线观看| 波多野结衣av一区二区av| 寂寞人妻少妇视频99o| 久久97久久精品| 午夜免费观看性视频| 26uuu在线亚洲综合色| 亚洲精品第二区| av.在线天堂| 中文字幕av电影在线播放| 久久毛片免费看一区二区三区| 日日撸夜夜添| 国产爽快片一区二区三区| 2022亚洲国产成人精品| 下体分泌物呈黄色| 日韩一本色道免费dvd| 一级片免费观看大全| 美女高潮到喷水免费观看| 黑丝袜美女国产一区| 亚洲欧美成人综合另类久久久| 国产免费一区二区三区四区乱码| 免费黄频网站在线观看国产| 国产有黄有色有爽视频| 一级片'在线观看视频| 亚洲欧美一区二区三区久久| 国产精品久久久av美女十八| 午夜日本视频在线| 王馨瑶露胸无遮挡在线观看| 一区福利在线观看| 各种免费的搞黄视频| 18禁观看日本| 街头女战士在线观看网站| 免费在线观看视频国产中文字幕亚洲 | 亚洲av在线观看美女高潮| 日本av免费视频播放| 亚洲欧美精品综合一区二区三区 | 午夜日韩欧美国产| 日韩熟女老妇一区二区性免费视频| 中国国产av一级| 欧美在线黄色| 国产精品人妻久久久影院| 成年av动漫网址| 日本av免费视频播放| 国产精品免费视频内射| 欧美精品人与动牲交sv欧美| www日本在线高清视频| 最近中文字幕2019免费版| 久久综合国产亚洲精品| 人妻一区二区av| 欧美中文综合在线视频| 日韩电影二区| 国产成人免费观看mmmm| 97精品久久久久久久久久精品| 精品99又大又爽又粗少妇毛片| videossex国产| 久久精品久久久久久久性| 欧美黄色片欧美黄色片| 成人国产麻豆网| 又黄又粗又硬又大视频| 18禁观看日本| 少妇被粗大的猛进出69影院| 国产精品无大码| 最近中文字幕高清免费大全6| 国产深夜福利视频在线观看| 亚洲欧美精品自产自拍| 国产激情久久老熟女| 久久精品久久久久久久性| 在线天堂中文资源库| 国产精品人妻久久久影院| www.自偷自拍.com| 久久久精品国产亚洲av高清涩受| 黑人巨大精品欧美一区二区蜜桃| 最新中文字幕久久久久| 久久久久久久国产电影| 欧美日韩视频精品一区| 晚上一个人看的免费电影| 成人手机av| 韩国av在线不卡| 久久免费观看电影| 亚洲精品第二区| 久久久亚洲精品成人影院| 亚洲国产色片| a级毛片在线看网站| 成年人午夜在线观看视频| 亚洲第一区二区三区不卡| 精品久久久精品久久久| 人体艺术视频欧美日本| 高清欧美精品videossex| 亚洲av电影在线进入| 91精品国产国语对白视频| 看免费成人av毛片| 啦啦啦在线免费观看视频4| 深夜精品福利| av在线播放精品| 精品人妻熟女毛片av久久网站| 久久热在线av| 久久久久精品人妻al黑| 精品少妇久久久久久888优播| 亚洲男人天堂网一区| 国产日韩一区二区三区精品不卡| av视频免费观看在线观看| 高清av免费在线| 亚洲三级黄色毛片| 黑人猛操日本美女一级片| 日韩一区二区视频免费看| 亚洲内射少妇av| 成年女人在线观看亚洲视频| 丰满乱子伦码专区| 亚洲欧美成人精品一区二区| 综合色丁香网| 五月开心婷婷网| 日日摸夜夜添夜夜爱| 岛国毛片在线播放| 9191精品国产免费久久| 精品国产露脸久久av麻豆| 一区二区三区乱码不卡18| 久久国产精品大桥未久av| 9191精品国产免费久久| 免费观看在线日韩| 成年人午夜在线观看视频| 永久免费av网站大全| 久久久久精品性色| 亚洲av中文av极速乱| 一边摸一边做爽爽视频免费| 午夜av观看不卡| 国产一级毛片在线| 一区二区av电影网| 女人高潮潮喷娇喘18禁视频| 18禁裸乳无遮挡动漫免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 一边亲一边摸免费视频| 少妇被粗大的猛进出69影院| 久久人人爽av亚洲精品天堂| 亚洲精品美女久久久久99蜜臀 | 久久久久精品久久久久真实原创| 十八禁高潮呻吟视频| 菩萨蛮人人尽说江南好唐韦庄| 国产乱人偷精品视频| 亚洲av男天堂| 亚洲三级黄色毛片| 欧美中文综合在线视频| a级片在线免费高清观看视频| 亚洲av电影在线观看一区二区三区| 亚洲国产精品999| 国产一区二区 视频在线| 免费在线观看完整版高清| 久久热在线av| 国产精品.久久久| 99久久综合免费| 亚洲av电影在线进入| 性色avwww在线观看| 天堂中文最新版在线下载| 女性被躁到高潮视频| 电影成人av| 秋霞在线观看毛片| 99香蕉大伊视频| 欧美成人午夜免费资源| 成年女人毛片免费观看观看9 | 秋霞在线观看毛片| av免费观看日本| 欧美成人午夜免费资源| 日韩中文字幕视频在线看片| 日韩一区二区视频免费看| 欧美精品av麻豆av| 91精品三级在线观看| 国产免费视频播放在线视频| av.在线天堂| 丝袜脚勾引网站| 久久久久久久久久久久大奶| 妹子高潮喷水视频| 精品一区二区三区四区五区乱码 | av不卡在线播放| 伊人亚洲综合成人网| 少妇 在线观看| 国产极品天堂在线| 久热久热在线精品观看| 亚洲av男天堂| 久久精品aⅴ一区二区三区四区 | 精品人妻在线不人妻| 成人18禁高潮啪啪吃奶动态图| 亚洲,欧美,日韩| av国产精品久久久久影院| 亚洲人成77777在线视频| 国产精品久久久久成人av| 欧美精品人与动牲交sv欧美| 亚洲成人一二三区av| 国产国语露脸激情在线看| 成人国语在线视频| 国产精品久久久久久久久免| 一级黄片播放器| 五月开心婷婷网| 国产精品一区二区在线不卡| 精品亚洲成a人片在线观看| 天天操日日干夜夜撸| 免费观看无遮挡的男女| 久久亚洲国产成人精品v| 日韩中文字幕视频在线看片| 成年av动漫网址| 国产精品av久久久久免费| 国产精品人妻久久久影院| 少妇熟女欧美另类| 成年人免费黄色播放视频| 日日撸夜夜添| 免费黄网站久久成人精品| 两性夫妻黄色片| 少妇熟女欧美另类| 国产成人精品在线电影| 亚洲婷婷狠狠爱综合网| 天天操日日干夜夜撸| 一区二区三区激情视频| 国产一区二区在线观看av| 精品一区二区三卡| 国产一区二区 视频在线| 亚洲中文av在线| 国产在视频线精品| 日韩人妻精品一区2区三区| 韩国高清视频一区二区三区| 亚洲欧洲国产日韩| 美女中出高潮动态图| 国产精品一区二区在线不卡| 色94色欧美一区二区| 老司机亚洲免费影院| 只有这里有精品99| 亚洲国产精品国产精品| 亚洲内射少妇av| 中文字幕人妻丝袜一区二区 | 亚洲五月色婷婷综合| 男人操女人黄网站| av在线老鸭窝| 成人影院久久| 伦理电影大哥的女人| 日本欧美国产在线视频| 中文欧美无线码| 熟妇人妻不卡中文字幕| 国产男女内射视频| www.av在线官网国产| 99久国产av精品国产电影| 视频在线观看一区二区三区| 欧美成人午夜免费资源| √禁漫天堂资源中文www| 捣出白浆h1v1| 热99国产精品久久久久久7| 亚洲av国产av综合av卡| 99久久人妻综合| 亚洲伊人久久精品综合| 婷婷色av中文字幕| 丝袜喷水一区| 99热网站在线观看| 中文字幕制服av| 中文欧美无线码| 国产成人aa在线观看| 男女下面插进去视频免费观看| 女性被躁到高潮视频| 最近最新中文字幕大全免费视频 | 亚洲精品久久久久久婷婷小说| 中文乱码字字幕精品一区二区三区| 国产有黄有色有爽视频| 美女福利国产在线| 性高湖久久久久久久久免费观看| 日日摸夜夜添夜夜爱| 大香蕉久久网| 毛片一级片免费看久久久久| 国产精品一区二区在线观看99| 一二三四中文在线观看免费高清| 亚洲精品美女久久av网站| 久久ye,这里只有精品| 午夜福利,免费看| av片东京热男人的天堂| 丝袜在线中文字幕| 久久精品久久久久久噜噜老黄| 久久精品国产亚洲av天美| 七月丁香在线播放| 国产亚洲精品第一综合不卡| 免费不卡的大黄色大毛片视频在线观看| 亚洲第一av免费看| 2018国产大陆天天弄谢| 一区二区日韩欧美中文字幕| 国产片内射在线| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 一区二区三区乱码不卡18| 你懂的网址亚洲精品在线观看| 热99国产精品久久久久久7| 伊人亚洲综合成人网| 免费黄色在线免费观看| 街头女战士在线观看网站| 如何舔出高潮| 亚洲国产欧美在线一区| 久久鲁丝午夜福利片| 欧美国产精品va在线观看不卡| 成人黄色视频免费在线看| 天天躁狠狠躁夜夜躁狠狠躁| 天美传媒精品一区二区| 在线看a的网站| 婷婷色av中文字幕| 一级片'在线观看视频| 国产在线免费精品| 最近的中文字幕免费完整| 精品久久久久久电影网| 亚洲第一青青草原| 又黄又粗又硬又大视频| 亚洲精品自拍成人| 晚上一个人看的免费电影| 亚洲国产av影院在线观看| 国产精品久久久av美女十八| 久久久国产精品麻豆| 最近中文字幕高清免费大全6| 午夜影院在线不卡| 搡老乐熟女国产| 大香蕉久久网| 夜夜骑夜夜射夜夜干| 精品人妻偷拍中文字幕| 亚洲国产欧美日韩在线播放| 99热网站在线观看| 日韩欧美一区视频在线观看| 9色porny在线观看| 亚洲成色77777| 热99国产精品久久久久久7| 天天躁夜夜躁狠狠久久av| 秋霞在线观看毛片| 69精品国产乱码久久久| 日日啪夜夜爽| 男女边吃奶边做爰视频| 亚洲久久久国产精品| 好男人视频免费观看在线| 欧美精品一区二区大全| 欧美日韩亚洲高清精品| 美女国产视频在线观看| 考比视频在线观看| 久久国产亚洲av麻豆专区| 亚洲欧美精品自产自拍| 热re99久久精品国产66热6| 女人精品久久久久毛片| 久热这里只有精品99| 午夜免费鲁丝| 欧美最新免费一区二区三区| 九色亚洲精品在线播放| 观看av在线不卡| 日本91视频免费播放| 大片电影免费在线观看免费| 亚洲精品一二三| 久久青草综合色| 天堂中文最新版在线下载| 韩国av在线不卡| 国产极品粉嫩免费观看在线| 一个人免费看片子| 国产精品三级大全| 在线观看免费视频网站a站| 亚洲精品久久午夜乱码| 久久女婷五月综合色啪小说| 亚洲婷婷狠狠爱综合网| 国产高清不卡午夜福利| 亚洲av免费高清在线观看| 亚洲一区二区三区欧美精品| 国产精品久久久久久精品古装| 又大又黄又爽视频免费| 日本欧美视频一区| 成人影院久久| 久久久国产精品麻豆| 亚洲一区中文字幕在线| 汤姆久久久久久久影院中文字幕| 三上悠亚av全集在线观看| 少妇的丰满在线观看| 少妇被粗大的猛进出69影院| 国精品久久久久久国模美| 免费观看在线日韩|