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

    巢湖典型流域河道沉積物磷賦存形態(tài)研究

    2021-03-22 03:34:48姜曉霞胡晨雨方喆禹

    姜曉霞,孫 玥,儲 茵,胡晨雨,王 汛,方喆禹

    巢湖典型流域河道沉積物磷賦存形態(tài)研究

    姜曉霞,孫 玥,儲 茵*,胡晨雨,王 汛,方喆禹

    (安徽農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,合肥 230036)

    河道沉積物磷的含量和形態(tài)可為流域磷的來源解析及評估其對下游受納水體的影響提供重要信息。選取巢湖典型支流豐樂河和柘皋河,在平水期對兩個流域共計27個采樣點進行了沉積物和河水采樣,沉積物采用SMT(standards measurements testing)法進行分級測定磷的不同形態(tài),包括總磷(TP)、無機磷(IP)、有機磷(OP)、鐵鋁磷(Fe/Al-P)以及鈣磷(Ca-P),水樣測定總磷(TP)和可溶性磷(DP)。結(jié)果表明,豐樂河沉積物中TP、IP和OP的平均含量分別為522.12、323.78和140.42 mg·kg-1,而柘皋河分別為1 060.70、755.74和125.31 mg·kg-1,均表現(xiàn)為TP>IP>OP。豐樂河IP中的Fe/Al-P平均含量高于Ca-P,而柘皋河則相反。柘皋河沉積物TP、IP、Fe/Al-P和Ca-P總體均高于豐樂河,但豐樂河OP含量略高于柘皋河。兩河不同采樣點沉積物不同形態(tài)磷含量均存在較大的空間差異。豐樂河和柘皋河河水TP平均值分別為0.10和0.18 mg·L-1,均以溶解態(tài)為主,易導(dǎo)致富營養(yǎng)化。兩個流域都是農(nóng)業(yè)流域,河流中磷主要來源是農(nóng)業(yè)非點源和城鎮(zhèn)生活污水的輸入。柘皋河農(nóng)田和村鎮(zhèn)都相對比較集中分布在河邊,不同來源的磷更容易輸送到達河道。柘皋河流域距入湖口近,沉積物和河水磷含量均高于豐樂河,因此柘皋河對加劇巢湖水體富營養(yǎng)化的風(fēng)險更大,但豐樂河沉積物Fe/Al-P和河水磷含量也比較高,依然存在一定的風(fēng)險。

    沉積物;磷形態(tài);河流;巢湖流域

    巢湖位于長江中下游,是我國第五大淡水湖。巢湖水體富營養(yǎng)化嚴重,經(jīng)過控制工廠及城鎮(zhèn)點源污染廢水的排放,水體中總氮總磷有下降趨勢,但仍保持較高水平。氮磷富集是導(dǎo)致水質(zhì)惡化、藍藻爆發(fā)的重要原因,其中磷是水體富營養(yǎng)化的主要限制性營養(yǎng)元素[1]。水體中的磷主要來自外源輸入,工業(yè)廢水、生活污水、農(nóng)業(yè)面源污染等通過徑流進入河流,最終匯入巢湖。另一方面,在平水期及洪水退水時期河道水體流速緩慢,從而部分沉降在河道沉積物中[2],而沉積物隨下一次洪水?dāng)y帶也將匯入巢湖。

    沉積物中磷形態(tài)主要有無機磷和有機磷,無機磷中又包括金屬氧化物結(jié)合態(tài)磷、鈣結(jié)合態(tài)磷、還原態(tài)磷、弱吸附態(tài)磷等[3]。沉積物中不同形態(tài)磷組成過程復(fù)雜,受溶解氧、pH、溫度、水力條件、沉積物組分等主要環(huán)境因素影響大[4]。沉積物中的磷不僅是水體內(nèi)源磷污染的重要來源,不同形態(tài)磷具有不同的生物有效性,也會對水體富營養(yǎng)化有潛在的風(fēng)險[5-6]。

    目前在對沉積物中不同形態(tài)磷的研究中,有Golterman等利用鰲合物EDTA提取沉積物中的鐵結(jié)合態(tài)磷、鈣結(jié)合態(tài)磷、酸可提取有機磷和堿可提取有機磷來區(qū)分不同磷灰石的形態(tài),但是DTA與EDTA對磷的提取有干擾,預(yù)處理復(fù)雜,提取過程有重復(fù)[7];李悅等提出七步連續(xù)提取法,將沉積物中的磷以吸附態(tài)磷、鋁結(jié)合磷、鐵結(jié)合磷、閉蓄態(tài)磷、自生鈣結(jié)合磷、原生碎屑磷和有機磷,區(qū)分了原生碎屑磷和自生鈣結(jié)合磷,但其主要針對海洋沉積物磷的分級測定,且步驟過于繁瑣[8];在1999年歐洲標準體系下,Ruban等研究了沉積物中磷形態(tài)的分級測定方法(SMT, standards measurements testing),提取沉積物中總磷、無機磷、有機磷、鐵鋁結(jié)合態(tài)磷和鈣結(jié)合態(tài)磷,方法簡單實用,可提供生物可利用磷的信息,且對土壤和底泥均適用[9-10]。SMT法在我國應(yīng)用廣泛,常用于水庫及湖泊沉積物磷的分級測定,如:三峽庫區(qū)[11]、福建山美水庫[12]、湖北天福廟水庫[13]、東洞庭湖[14]等;以及河道沉積物中磷的測定,如合肥市板橋河蔡田鋪污水廠尾水受納河段[15]。

    在中小流域的尺度上,采用簡單實用的磷分級測定方法來研究河道沉積物磷的不同形態(tài)和空間分布,進而探究磷的來源解析以及對下游受納水體的影響,對下游受納水體富營養(yǎng)化的治理與管理有重要意義,但這方面的研究尚比較缺乏。因此,作者選取巢湖典型支流豐樂河及柘皋河為研究對象,通過對流域分布式布點,采集河流沉積物及水樣品,再選用SMT法進行測定,分析巢湖支流水體及沉積物中不同形態(tài)磷的含量及分布特征,以期為開展小流域面源污染防治及湖泊支流統(tǒng)籌管理提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域概況

    豐樂河發(fā)源于六安市金安區(qū)張店鎮(zhèn)境內(nèi),流經(jīng)舒城縣張母橋、桃溪和豐樂三鎮(zhèn)[16],于三河鎮(zhèn)與杭埠河自西匯合后注入巢湖全長112 km。豐樂河以桃溪水文站為流域出口,劃分集水區(qū)域,流域面積為1 579 km2,耕地、林地、居民區(qū)面積分別占總面積的63%、19%和7%。河水匯入杭埠河流入巢湖,從流域出口至巢湖入口有47 km。

    柘皋河流經(jīng)柘皋鎮(zhèn)、夏閣鎮(zhèn)和中垾鎮(zhèn),全長 35 km。以柘皋大橋為流域出口劃分集水區(qū)域,流域面積為497 km2,其中耕地、林地、居民區(qū)面積分別占總面積的75%、7%和8%。河水直接注入巢湖,從流域出口至巢湖入口僅有1.6 km。豐樂河、柘皋河與巢湖的位置分布如圖1所示。

    兩河流域降水年際、年內(nèi)分布較不均勻。主要降雨過程發(fā)生在夏季,易發(fā)生暴雨及大暴雨,洪水量大,持續(xù)時間長,屬亞熱帶季風(fēng)氣候,溫和濕潤[17]。流域以丘陵為主,坡度較小,主要用地為農(nóng)業(yè)和林地用地,流域內(nèi)沒有大的工廠及城鎮(zhèn)。主要糧食作物為水稻、小麥,主要土壤類型為水稻土、黃棕壤。

    1.2 樣品的采集和處理

    通過對兩流域?qū)嵉乜疾?,基于流域?shù)字高程DEM數(shù)據(jù)利用ArcGIS對流域劃分河網(wǎng),在主要支流入河口處布設(shè)采樣點;在支流上游根據(jù)面積大小增加布設(shè)采樣點;綜合考慮流域支流長度、土地利用、人類活動及取樣條件等因素,在干流中下游增設(shè)1~2個采樣點。柘皋河中垾鎮(zhèn)至柘皋鎮(zhèn)河段周邊均是圩區(qū),無支流,故未布設(shè)采樣點。將每個采樣點作為子流域出口斷面,劃分子流域。每個采樣點水質(zhì)和沉積物均能反映所在區(qū)域河流污染的來源和輸送特征。如圖2所示,在豐樂河布設(shè)了14個采樣點,流域出口布設(shè)采樣點F1在桃溪水文站附近;柘皋河布設(shè)了13個采樣點,流域出口布設(shè)采樣點Z1在柘皋大橋。

    圖1 巢湖流域、豐樂河及柘皋河位置

    Figure 1 Chao Lake basin and location of Fengle River and Zhegao River

    本次取樣時間為2020年5月21日和22日,屬于春季取樣。根據(jù)安徽省水文局公布的2020年5月桃溪水文站實測水位數(shù)據(jù)可知,當(dāng)月16日11: 00出現(xiàn)月最高水位為9.07 m,當(dāng)月8日18:00出現(xiàn)月最低水位為8.61 m,可知當(dāng)月降雨不充沛,水位不高,屬于平水期。按照已經(jīng)布設(shè)的取樣斷面在豐樂河、柘皋河上、中、下游分別進行取樣。根據(jù)取樣斷面河道寬度和取樣難易程度,采用多點法用水質(zhì)采樣器采集河流水面下0.3 ~ 0.5 m處的表層水樣,混合水樣裝入500 mL 的聚乙烯塑料瓶中,放置低溫冷藏保溫箱中帶回;用抓斗式采泥器多點采集采樣點河道底泥,樣品封裝于聚乙烯塑料袋中帶回實驗室。水樣及時測定,底泥自然風(fēng)干,過1 mm孔徑篩后裝塑封袋標號待測。

    1.3 樣品的分析方法

    沉積物中磷分級測定采用歐洲標準測試委員會框架下的SMT法定義的總磷(TP)、無機磷(IP)、有機磷(OP)、鐵鋁磷(Fe/Al-P)以及鈣磷(Ca-P)[9-10,18-20]。該方法分3次提取測定,相比一次稱樣連續(xù)分級測定損耗更少,方法簡單實用,適合對流域尺度環(huán)境影響分析使用[9-10,21]??偭诇y定將0.2 g待測底泥樣品置于銀坩堝中高溫煅燒,使底泥中含磷礦物及有機磷化合物全部轉(zhuǎn)化為可溶性的正磷酸鹽;再取0.2 g底泥樣品用鹽酸提取,使底泥中含磷礦物全部轉(zhuǎn)化為可溶性的正磷酸鹽,用于測定無機磷含量;有機磷化合物位于無機磷提取后的殘渣中,經(jīng)高溫煅燒,將有機磷化合物全部轉(zhuǎn)化為可溶性的正磷酸鹽,再用鹽酸提?。辉偃?.2 g待測底泥樣品用氫氧化鈉溶液提取鐵鋁磷,剩余殘渣清洗后用鹽酸提取鈣磷。所有提取液均用鉬藍比色后測定吸光度。具體實驗步驟如圖3所示[22]。

    圖2 豐樂河(a)、柘皋河(b)采樣點及子流域圖

    Figure 2 Sampling sites and sub-watersheds in Fengle River (a) and Zhegao River (b)

    水體中總磷(TP)采用過硫酸鉀消解鉬酸銨分光光度法(GB/T 11 893-1 989)測定,可溶性總磷(DP)是將水樣過濾0.45 μm濾膜后,用鉬酸銨分光光度法(GB/T 11 893-1 989)測定。

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

    使用Excel 2016、SPSS 20進行數(shù)據(jù)統(tǒng)計分析,利用Origin 2018、ArcGIS 10.2進行繪圖,對兩河流不同采樣點的各形態(tài)磷進行差異分析。

    圖3 SMT磷分級測定

    Figure 3 The SMT protocol

    2 結(jié)果與分析

    2.1 河道沉積物中TP含量及其空間分布特征

    由圖4可知,豐樂河14個采樣點沉積物中TP含量介于145 ~ 1 015 mg·kg-1之間,各采樣點平均含量達到522 mg·kg-1,較巢湖春季表層沉積物TP含量平均值980 mg·kg-1相比[23],豐樂河沉積物TP含量處于中下水平,是巢湖沉積物TP含量的53%。如圖5(a)所示,不同采樣點的TP含量水平不同,Duncan分析顯示組內(nèi)水平有顯著差異。其中,F(xiàn)2采樣斷面底泥較難抓取,多處為石塊和板結(jié)巖石,經(jīng)實地調(diào)查,該處河道有過清淤,底泥沉積并不多,受來水沖擊影響大;F7采樣點出現(xiàn)高值與附近農(nóng)村居民區(qū)和大面積農(nóng)業(yè)種植有關(guān);F13采樣點TP含量最高達1 015 mg·kg-1,該采樣點附近沿河有農(nóng)戶的雞舍鴨棚,雞鴨代謝產(chǎn)生的糞便會直接排入河岸邊,降雨后會隨徑流進入水體,積累到沉積物中??梢姡浇囟▉碓磳Τ练e物TP含量的也有較大影響。

    柘皋河13個采樣點沉積物中TP含量介于534 ~1 611 mg·kg-1之間,各采樣點平均含量達到1 060 mg·kg-1,約是豐樂河表層沉積物TP含量的2倍,已達巢湖表層沉積物TP含量的108%[23],略低于杭州西湖表層沉積物TP含量[24]平均值1 209 mg·kg-1,高于太湖西北部竺山灣沉積物TP含量[25]平均值930 mg·kg-1(圖4)。柘皋河TP不同采樣點兩兩之間差異也很顯著。如圖5(b)所示,Z12采樣點的沉積物TP含量最高,達1 611.05 mg·kg-1,該處水葫蘆覆蓋水面,水體藻類懸浮物多,濁度高,透明度低,其污染源主要是河道一側(cè)前李村的生活污水和另一側(cè)農(nóng)業(yè)蔬菜大面積種植的非點源污染輸入所致;Z8、Z10和Z3均分布在城鎮(zhèn)居民區(qū),人口密度較高,受人類活動影響大;Z1采樣點的沉積物TP含量最低,為 533.88 mg·kg-1,該處為柘皋河流域的出口斷面,水位較深,河道較寬,可行船,下游是柘皋河匯入巢湖的入湖口,水面水生植物較少,河底有清淤及橋底水泥硬化,河道中央沉積物較難抓取,所以只取了岸邊底泥,代表性略低。柘皋河整體水域的沉積物磷污染較豐樂河乃至巢湖更為嚴重。

    圖 4 豐樂河、柘皋河沉積物中不同形態(tài)磷含量箱線圖

    Figure 4 Box plot of concentrations of phosphorus fractions in sediments of Fengle River and Zegao River

    2.2 河道沉積物中IP與OP含量及其空間分布特征

    根據(jù)SMT法,首先將沉積物中磷的形態(tài)分為無機和有機兩種。如圖6(a)所示,豐樂河14個采樣點沉積物中IP含量在56.94 ~ 792.34 mg·kg-1之間,各采樣點平均含量達到323.78 mg·kg-1,占TP 的比例為62%;OP含量在65.67 ~ 264.63 mg·kg-1之間,各采樣點平均含量達到140.42 mg·kg-1,占TP 的比例為27%。如圖6(b)所示,柘皋河13個采樣點沉積物中IP含量在333.7 ~ 1 349.15 mg·kg-1之間,各采樣點平均含量達到775.74 mg·kg-1,占TP 的比例為73%;OP含量在59.80 ~ 183.67 mg·kg-1之間,各采樣點平均含量為125.31 mg·kg-1,占TP 的比例為12%。說明河流表層沉積物中磷的主要成分為IP。如圖6 所示,豐樂河IP含量明顯低于柘皋河,而OP含量水平兩河差距不明顯。不同采樣點IP、OP 含量水平不同,經(jīng)組內(nèi)分析,兩兩之間隨子流域變化,IP和OP均有顯著差異,說明采樣點位置不同對磷含量有顯著影響。

    圖5 豐樂河(a)、柘皋河(b)不同采樣點TP含量

    Figure 5 Concentrations of TP in sediments at different sampling sites of Fengle River (a) and Zhegao River (b)

    圖6 豐樂河(a)、柘皋河(b)不同采樣點沉積物IP與OP含量

    Figure 6 Concentrations of IP and OP in sediments at different sampling sites in Fengle River (a) and Zhegao River (b)

    圖7 豐樂河(a)、柘皋河(b)不同采樣點Fe/Al-P與Ca-P含量

    Figure 7 Concentrations of Fe/Al-P and Ca-P in sediment at different sampling sites of Fengle River (a) and Zegao River (b)

    2.3 河道沉積物中Fe/Al-P與Ca-P含量及其空間分布特征

    IP 的主要形態(tài)為Fe/Al-P和Ca-P。如圖7(a)所示,豐樂河各采樣點間沉積物中Fe/Al-P含量為37.77 ~ 387.70 mg·kg-1,流域平均含量達到161.22 mg·kg-1,Ca-P含量為15.96 ~ 382.40 mg·kg-1,流域平均含量達到142.97 mg·kg-1,各采樣點均值Fe/Al-P占IP的比例為49.8 %,Ca-P占IP的比例為44.2%。在豐樂河多數(shù)斷面,F(xiàn)e/Al-P在IP中所占比例最高,F(xiàn)e/Al-P是指被鐵鋁氧化物及水化合物包裹的磷,不穩(wěn)定,易轉(zhuǎn)化,在還原環(huán)境下幾乎可以全部釋放,即可在一定條件下能夠被生物所利用的活性磷或潛在活性磷,沉積物中這部分磷的含量能真正反映沉積物內(nèi)源釋放能力的大小[26-27]。它的來源受人類活動影響大,主要來自生活污水、工業(yè)廢水的點源排放以及農(nóng)業(yè)面源污染隨降雨徑流的進入,F(xiàn)e/Al-P的含量可反映水體受外源污染的程度及沉積物中磷向水體內(nèi)源釋放的風(fēng)險。

    圖7(b)顯示,柘皋河各采樣點間沉積物中Fe/Al-P含量為85.60 ~ 609.67 mg·kg-1,流域平均含量達到250.36 mg·kg-1,Ca-P含量為160.25 ~ 687.07 mg·kg-1,流域平均含量達到307.89 mg·kg-1,各采樣點均值 Fe/Al-P占IP的比例為32.3%,Ca-P占IP的比例為38.7%。在柘皋河多數(shù)斷面,Ca-P在IP中所占比例最高,可能與流域內(nèi)有磷礦以及土壤類型有棕色石灰土有關(guān)。Ca-P主要是來自碳酸鈣、磷灰石以及含磷礦物相關(guān)的沉積磷,形態(tài)穩(wěn)定,不易直接釋放到水體中,一般難以被生物利用,因此可認為多數(shù)斷面水體中磷污染來自Ca-P內(nèi)源釋放的可能性并不大[28]。

    圖 8 豐樂(a)、柘皋(b)不同采樣點河水TP與DP含量

    Figure 8 Concentrations of TP and DP in water at different sites of Fengle River (a) and Zhegao River (b)

    表1 豐樂河沉積物各形態(tài)磷及水樣總磷間的相關(guān)性分析(n=14)

    注:*<0.05,**<0.01,W-TP為水樣總磷。下同。

    表2 柘皋河沉積物各形態(tài)磷及水樣總磷間的相關(guān)性分析(n=13)

    如圖7所示,不同采用點Fe/Al-P與Ca-P含量差異大,組內(nèi)分析兩兩采樣點間Fe/Al-P或Ca-P大多數(shù)都具有顯著性差異,子流域不同對均值有顯著影響。

    2.4 水中總磷、可溶性磷含量及其空間分布特征

    如圖8(a)所示,豐樂河河水TP含量介于0.03~0.15mg·L-1之間,均值約為0.10 mg·L-1,河水DP含量在0.02~0.11 mg·L-1之間,均值為0.06 mg·L-1,DP占TP的比例67.1%, F11采樣點河水總磷及可溶性磷含量均較高;由圖8(b)所示,柘皋河河水TP含量為0.07~0.59 mg·L-1之間,均值為0.18 mg·L-1,河水DP含量為0.05~0.45 mg·L-1之間,均值為0.12 mg·L-1,DP占TP的比例65.9%。湖泊水質(zhì)富營養(yǎng)化的評價標準中TP<0.10 mg·L-1為富營養(yǎng)化、0.10 mg·L-1<TP為重度富營養(yǎng)化[29]。顯然,兩河水體磷含量均已達到富營養(yǎng)化限制水平,有向下游水體輸送超標富營養(yǎng)化元素的風(fēng)險。根據(jù)《地表水環(huán)境質(zhì)量標準》(GB3838-—2002)中TP含量標準限值,豐樂河一半采樣點水質(zhì)均達Ⅱ類水,另一半采樣點水質(zhì)達Ⅲ類水;柘皋河大部分采樣點水質(zhì)達Ⅲ類水,Z12采樣點水體TP含量超Ⅴ類水0.4 mg·L-1的限值,該點沉積物中TP含量也高達1 611 mg·kg-1。柘皋河水體TP均值是豐樂河的1.8倍,DP均值是豐樂河的2倍,柘皋河河水TP、DP含量總體均高于豐樂河河水。豐樂河、柘皋河DP占TP的比例基本一致,說明平水期河流水體中磷的主要成分為可溶解性磷。

    3 討論與結(jié)論

    進一步比較兩個河流沉積物中不同形態(tài)磷之間的特征,對兩個河流沉積物中不同形態(tài)磷賦存含量及水體中總磷含量之間分別進行皮爾遜相關(guān)性分析,結(jié)果分別見表1和表2。在豐樂河中,沉積物中TP含量與IP含量、OP含量、Fe/Al-P含量和Ca-P含量均顯著正相關(guān)(表1),而沉積物中各形態(tài)磷與水體總磷均無相關(guān)性,說明了豐樂河水體TP與沉積物中各形態(tài)磷影響不大。這可能與平水期田塊沖刷較少、河水多來自地下水補給、豐樂河水體流速快與底泥交換時間相對較短等原因有關(guān)。

    由柘皋河沉積物各形態(tài)磷及水樣總磷間的皮爾遜相關(guān)性分析(表2)可得,柘皋河沉積物中TP含量與IP含量、OP含量和Fe/Al-P含量均顯著正相關(guān),Ca-P含量與TP、IP、OP和Fe/Al-P含量相關(guān)性均不明顯,柘皋河流域內(nèi)有磷礦可能影響Ca-P含量。這一結(jié)果與江立文等[26]在探究中度營養(yǎng)型湖泊與輕度營養(yǎng)型湖泊沉積物中磷形態(tài)特征及主成分分析結(jié)果是一致的。柘皋河水體總磷與沉積物中TP、IP、OP和Fe/Al-P均正相關(guān),說明了柘皋河水體TP與沉積物中各形態(tài)磷關(guān)系密切。這可能與柘皋河水體流速小、水體相對靜止、有充足的時間供底泥與水體交換且柘皋河沿河生活污水輸入較多等原因有關(guān)。

    從河流沉積物磷賦存形態(tài)的平均含量來看,豐樂河沉積物中TP、IP、OP、Fe/Al-P和Ca-P的平均含量分別為522.12、323.78、140.42、161.22和142.97 mg·kg-1,TP>IP>Fe/Al-P>Ca-P>OP,這與江立文等[26]對江西湖泊的研究結(jié)果是一致的。柘皋河沉積物中TP、IP、OP、Fe/Al-P和Ca-P的平均含量分別為1 060.70、755.74、125.31、250.36和307.89 mg·kg-1,TP>IP>Ca-P>Fe/Al-P>OP。柘皋河沉積物各形態(tài)磷含量不同于豐樂河的是Ca-P含量高于Fe/Al-P,且除OP略低于豐樂河含量外,其他含量均高于豐樂河。

    豐樂河流域內(nèi)村鎮(zhèn)分散,河道水體水深較淺,大面積農(nóng)業(yè)種植水稻。柘皋河流域面積相對較小,村鎮(zhèn)集中,河道寬而深,多圩區(qū),農(nóng)業(yè)種植多水稻、小麥和蔬菜。兩個河流磷主要來源都是農(nóng)業(yè)非點源污染以及城鎮(zhèn)生活污水的點源排放。含量差異可能與柘皋河農(nóng)田和城鎮(zhèn)都相對比較集中分布在河邊,相對強度更大,也更容易輸送至河道有關(guān)。此外,柘皋河流域有礦區(qū),部分稻田發(fā)展稻蝦結(jié)合新型農(nóng)業(yè)模式,致使礦物質(zhì)及大量餌料殘渣及蝦的代謝產(chǎn)物容易隨徑流進入水體。

    柘皋河沉積物TP、Fe/Al-P和Ca-P和河水磷含量均高于豐樂河,且距巢湖入湖口更近,對加劇巢湖水體富營養(yǎng)化的風(fēng)險更大,但豐樂河沉積物Fe/Al-P和河水中磷含量也比較高,依然存在一定的風(fēng)險。因此,治理巢湖水體富營養(yǎng)化問題要同步將支流作為治理對象,加強群眾控磷意識,減少化肥用量,生活污水要處理后排放,沿河養(yǎng)殖場的畜禽糞便要妥善處理,距離入湖口近的河段要及時清淤等。

    通過在平水期對巢湖兩個典型子流域豐樂河和柘皋河沉積物的分布式采樣,并采用SMT方法對沉積物磷的不同賦存形態(tài)進行了測定,結(jié)果表明兩個河流不同采樣點沉積物不同形態(tài)磷含量差異較大。柘皋河沉積物TP含量總體高于豐樂河。從磷的不同形態(tài)來看,沉積物中IP和OP平均含量分別是:豐樂河為323.78和140.42 mg·kg-1,而柘皋為755.74和125.31 mg·kg-1,表明兩個流域都是以IP為主。而IP中Fe/Al-P、Ca-P的平均含量分別是:豐樂河為161.22和142.97 mg·kg-1,而柘皋為250.36和307.89 mg·kg-1,豐樂河Fe/Al-P平均含量高于Ca-P,而柘皋河則相反。柘皋河沉積物IP、Fe/Al-P和Ca-P平均含量均高于豐樂河,但豐樂河OP含量略高于柘皋河。豐樂河、柘皋河河水總磷平均值分別為0.10和0.18 mg·L-1,均以溶解態(tài)為主。

    除了OP,河道沉積物其他形態(tài)磷的含量均是柘皋河高于豐樂河,柘皋河水TP和DP濃度也總體高于豐樂河。兩個流域都主要是農(nóng)業(yè)流域,這種明顯的差異可能與柘皋河農(nóng)田和城鎮(zhèn)都相對比較集中分布在河邊有關(guān),農(nóng)業(yè)非點源來源的磷,特別是城鎮(zhèn)生活污水的點源排放相對強度更大,且更容易輸送到達河道。

    柘皋河相比豐樂河距巢湖入湖口近,且沉積物TP、Fe/Al-P和Ca-P含量和河水磷含量均高于豐樂河,因此柘皋河對加劇巢湖水體富營養(yǎng)化的風(fēng)險更大,但豐樂河沉積物Fe/Al-P和河水磷含量也比較高,依然存在一定的風(fēng)險。

    [1] 蘇玉萍, 賴壽輝, 林佳, 等. 富營養(yǎng)化飲用水源地山仔水庫限制性營養(yǎng)元素研究[J]. 環(huán)境科學(xué)學(xué)報, 2015, 35(10): 3107-3113.

    [2] LIU W B, WANG S R, ZHANG L, et al. Phosphorus release characteristics of sediments in Erhai Lake and their impact on water quality[J]. Environ Earth Sci, 2015, 74(5): 3753-3766.

    [3] 葉宏萌, 楊浩, 袁旭音, 等. 基于流域沉積物氮磷形態(tài)的生態(tài)風(fēng)險評價: 以沙溪流域為例[J]. 環(huán)境化學(xué), 2020, 39(12): 3471-3479.

    [4] S?NDERGAARD M, JENSEN J P, JEPPESEN E. Internal phosphorus loading in shallow Danish lakes[J]. Hydrobiologia, 1999, 408/409: 145-152.

    [5] 萬杰, 袁旭音, 葉宏萌, 等. 洪澤湖不同入湖河流沉積物磷形態(tài)特征及生物有效性[J]. 中國環(huán)境科學(xué), 2020, 40(10): 4568-4579.

    [6] VICENTE M A F, DE MELO G V, BAPTISTA NETO J A, et al. Phosphorus fractionation distribution in Guapimirim estuary: SE Brazil[J]. Springer Plus, 2016, 5(1): 1406.

    [7] GOLTERMAN H L. Fractionation of sediment phosphate with chelating compounds: a simplification, and comparison with other methods[J]. Hydrobiologia, 1996, 335(1): 87-95

    [8] 李悅, 烏大年, 薛永先. 沉積物中不同形態(tài)磷提取方法的改進及其環(huán)境地球化學(xué)意義[J]. 海洋環(huán)境科學(xué), 1998(1): 16-21.

    [9] RUBAN V, LóPEZ-SáNCHEZ J F, PARDO P, et al. Selection and evaluation of sequential extraction procedures for the determination of phosphorus forms in lake sediment[J]. J Environ Monit, 1999, 1(1): 51-56.

    [10] RUBAN V, LóPEZ-SáNCHEZ J F, PARDO P, et al. Harmonized protocol and certified reference material for the determination of extractable contents of phosphorus in freshwater sediments: a synthesis of recent works[J]. Fresenius J Anal Chem, 2001, 370(2/3): 224-228.

    [11] 張志永, 萬成炎, 胡紅青, 等. 三峽水庫干流沉積物及消落帶土壤磷形態(tài)及其分布特征[J]. 環(huán)境科學(xué), 2018, 39(9): 4161-4168.

    [12] 路丁, 郭沛涌, 沈芳芳, 等. 福建省山美水庫入庫河道沉積物磷釋放風(fēng)險[J]. 環(huán)境化學(xué), 2015, 34(8): 1498-1505.

    [13] 劉明盟, 李永福, 葛繼穩(wěn), 等. 宜昌市天福廟水庫沉積物磷形態(tài)分布特征及其釋放通量估算[J]. 環(huán)境科學(xué)研究, 2018, 31(7): 1258-1265.

    [14] 王婷, 王坤, 姜霞. 東洞庭湖沉積物覆水后磷形態(tài)變化及其釋放量[J]. 湖泊科學(xué), 2018, 30(4): 937-947.

    [15] 湯寧, 李如忠, 王聿慶, 等. 污水廠尾水受納河段沉積物磷形態(tài)及釋放風(fēng)險效應(yīng)[J]. 環(huán)境科學(xué), 2020, 41(2): 801-808.

    [16] 儲茵, 潮洪武, 馬友華, 等. 巢湖流域豐樂河洪水事件營養(yǎng)鹽輸出動態(tài)研究[J]. 長江流域資源與環(huán)境, 2013, 22(8): 1072-1080.

    [17] 儲茵, 朱江, 夏守先, 等. 巢湖典型支流柘皋河水質(zhì)污染時空變化特征[J]. 水土保持學(xué)報, 2011, 25(4): 243-248.

    [18] YANG W Q, XIAO H, LI Y, et al. Vertical distribution and release characteristics of phosphorus forms in the sediments from the river inflow area of Dianchi Lake, China[J]. Chem Speciat Bioavailab, 2018, 30(1): 14-22.

    [19] LOH P S, YING C Y, ALNOOR H I M, et al. Comparative study on the elucidation of sedimentary phosphorus species using two methods, the SMT and SEDEX methods[J]. J Anal Methods Chem, 2020, 2020: 8548126.

    [20] 朱夢圓, 朱廣偉, 錢君龍, 等. SMT法插標分析沉積物中磷的地球化學(xué)形態(tài)[J]. 中國環(huán)境科學(xué), 2012, 32(8): 1502-1507.

    [21] 孫境蔚. 沉積物磷的分級提取方法及提取相的共性分析[J]. 環(huán)境科學(xué)與技術(shù), 2007, 30(2): 111-114,121.

    [22] PARDO P, RAURET G, LóPEZ-SáNCHEZ J F. Shortened screening method for phosphorus fractionation in sediments[J]. Anal Chimica Acta, 2004, 508(2): 201-206.

    [23] 徐康, 劉付程, 安宗勝, 等. 巢湖表層沉積物中磷賦存形態(tài)的時空變化[J]. 環(huán)境科學(xué), 2011, 32(11): 3255-3263.

    [24] 嚴攀, 徐棟, 劉子森, 等. 杭州西湖沉積物磷分析及釋放風(fēng)險[J]. 環(huán)境化學(xué), 2019, 38(7): 1479-1487.

    [25] 方家琪, 祁闖, 張新厚, 等. 太湖竺山灣沉積物碳氮磷分布特征與污染評價[J]. 環(huán)境科學(xué), 2019, 40(12): 5367-5374.

    [26] 江立文, 李銘敏, 楊銀, 等. 江西省不同營養(yǎng)類型湖泊底泥沉積物中磷的賦存形態(tài)及其差異性[J]. 環(huán)境工程, 2020, 38(9): 48-52,167.

    [27] 王叢丹, 汪金成, 楊宇, 等. 平水期洞庭湖不同形態(tài)磷賦存特征[J]. 水文, 2019, 39(6): 74-79.

    [28] 何佳, 陳春瑜, 鄧偉明, 等. 滇池水-沉積物界面磷形態(tài)分布及潛在釋放特征[J]. 湖泊科學(xué), 2015, 27(5): 799-810.

    [29] FADIRAN A O, DLAMINI S C, MAVUSO A. A comparative study of the phosphate levels in some surface and ground water bodies of Swaziland[J]. Bull Chem Soc Eth, 2008, 22(2): 197-206.

    Study on phosphorus fractions in river sediments of typical watersheds of Chao Lake basin

    JIANG Xiaoxia, SUN Yue, CHU Yin, HU Chenyu, WANG Xun, FANG Zheyu

    (School of Resources and Environment, Anhui Agricultural University, Hefei 230036)

    The contents and fractions of phosphorus in river sediment can provide important information for source analysis and the evaluation of its influence on downstream receiving water bodies. In this study, spatially distributed sampling was applied during normal-low flow period on Fengle and Zhegao, which are the two inflow rivers of Chao Lake basin. Sediment and river water samples were taken at 14 sites on Fengle river and 13 on Zhegao. The SMT(standards measurements testing) protocol was applied for the analysis of five fractions of phosphorus, including total phosphorus (TP), inorganic phosphorus (IP), organic phosphorus (OP), iron-aluminum phosphorus (Fe/Al-P) and calcium phosphorus (Ca-P); river water samples were analyzed for TP and dissolved phosphorus (DP). The results showed that the average contents of TP, IP and OP in Fengle river sediments were 522.12, 323.78 and 140.42 mg·kg-1, respectively; and which in Zhegao river sediments were 1 060.70, 755.74 and 125.31 mg·kg-1, respectively; both demonstrating an order as TP>IP>OP. For IP, the average content of Fe/Al-P in Fengle river sediment was higher than that of Ca-P, while it was the opposite for Zhegao. The average contents of TP, IP, Fe/ Al-P and Ca-P in sediments of Zhegao river were higher than those in Fengle river, but the content of OP in Fengle was slightly higher than that in Zhegao. The concentrations of phosphorus fractions in sediments demonstrated large spatial variations for both rivers. The average TP concentrations in the water of Fengle and Zhegao were 0.10 and 0.18 mg·L-1, respectively, which were mainly composed of DP, thus at a risk of eutrophication. Both Fengle and Zhegao were of typical agricultural watersheds. Agricultural nonpoint sources inflow and the input of town sewage were the main contributions of phosphorus in river sediments and water. The crop and vegetable fields, towns or villages of Zhegao watershed were mostly located along or near the riverbanks, where phosphorus could be more easily transported to river

    sediment; phosphorus fractions; river; Chao Lake basin

    X522

    A

    1672-352X (2021)06-0960-08

    channels. As Zhegao watershed was closer to Chao Lake and the phosphorus content in both sediment and river water was higher than that in Fengle, it had a greater risk of worsening the eutrophication status in Chaohu Lake. However, since Fe/Al-P in sediments and TP in river water of Fengle river were relatively high as well, there was still a risk of contributing to the eutrophication in the downstream reaches and finally the Chao Lake.

    10.13610/j.cnki.1672-352x.20220106.005

    2022-1-10 8:58:11

    [URL] https://kns.cnki.net/kcms/detail/34.1162.S.20220107.1339.003.html

    2021-02-04

    安徽省自然科學(xué)基金(1908085MD102 )和安徽省大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計劃項目(ds195209)共同資助。

    姜曉霞,碩士研究生。E-mail:jiang__xiaoxia@126.com

    通信作者:儲 茵,教授。E-mail:chuyin@ahau.edu.cn

    又爽又黄无遮挡网站| 我的老师免费观看完整版| 午夜福利高清视频| АⅤ资源中文在线天堂| 成人性生交大片免费视频hd| 高清在线国产一区| 亚洲av二区三区四区| 久久精品人妻少妇| 亚洲国产欧洲综合997久久,| 日韩欧美精品免费久久| 国产 一区 欧美 日韩| 欧美+日韩+精品| 亚洲精品亚洲一区二区| 久久人人爽人人爽人人片va| 熟女人妻精品中文字幕| 看黄色毛片网站| 免费观看在线日韩| 国产私拍福利视频在线观看| 免费观看的影片在线观看| 亚洲图色成人| 亚洲人成网站在线播放欧美日韩| 国产精品免费一区二区三区在线| 成年人黄色毛片网站| 午夜久久久久精精品| 黄色配什么色好看| 欧美bdsm另类| 午夜福利在线观看吧| 午夜福利成人在线免费观看| 精品福利观看| 国产高清激情床上av| а√天堂www在线а√下载| 久久久久性生活片| 国产精品99久久久久久久久| 少妇被粗大猛烈的视频| 亚洲精品影视一区二区三区av| 中出人妻视频一区二区| 日韩精品中文字幕看吧| 国产高清激情床上av| 免费人成视频x8x8入口观看| 别揉我奶头 嗯啊视频| 日韩精品中文字幕看吧| 久久久久久久久大av| 毛片一级片免费看久久久久 | 人妻制服诱惑在线中文字幕| 男人舔女人下体高潮全视频| videossex国产| 国产精品久久电影中文字幕| 波野结衣二区三区在线| 色吧在线观看| 国产成人影院久久av| 国内久久婷婷六月综合欲色啪| 国产毛片a区久久久久| 久久国产乱子免费精品| 欧美黑人巨大hd| 精品日产1卡2卡| 在线播放无遮挡| 老熟妇仑乱视频hdxx| 九色成人免费人妻av| 日韩亚洲欧美综合| 国产亚洲精品av在线| 男人舔奶头视频| 乱人视频在线观看| 麻豆av噜噜一区二区三区| 午夜日韩欧美国产| 少妇人妻一区二区三区视频| 欧美日韩乱码在线| 亚洲欧美精品综合久久99| 伊人久久精品亚洲午夜| 日本 av在线| 免费搜索国产男女视频| 九九爱精品视频在线观看| 欧美激情久久久久久爽电影| avwww免费| 免费人成视频x8x8入口观看| 国产在线精品亚洲第一网站| ponron亚洲| 久久久久久国产a免费观看| 国产人妻一区二区三区在| 99精品在免费线老司机午夜| 国产精品久久久久久亚洲av鲁大| 欧美一区二区亚洲| 亚洲第一电影网av| 99久久久亚洲精品蜜臀av| 变态另类丝袜制服| 中文字幕精品亚洲无线码一区| 波多野结衣巨乳人妻| 亚洲av中文av极速乱 | 黄色一级大片看看| 国产精品一及| 精品一区二区三区视频在线| 91在线观看av| 欧美bdsm另类| 午夜精品在线福利| 99久久精品一区二区三区| 日本免费一区二区三区高清不卡| 18禁在线播放成人免费| 亚洲av不卡在线观看| 成人高潮视频无遮挡免费网站| eeuss影院久久| 亚洲aⅴ乱码一区二区在线播放| 亚洲国产精品久久男人天堂| 亚洲精品一卡2卡三卡4卡5卡| avwww免费| 精品久久久久久久末码| netflix在线观看网站| 日日啪夜夜撸| 日韩av在线大香蕉| 精品一区二区三区av网在线观看| 欧美一区二区精品小视频在线| 欧美一区二区国产精品久久精品| 色av中文字幕| 两个人的视频大全免费| 免费大片18禁| 精品久久久久久久久久久久久| 全区人妻精品视频| 露出奶头的视频| 日本黄色片子视频| 又紧又爽又黄一区二区| 国产黄a三级三级三级人| 亚洲最大成人手机在线| 免费av观看视频| 国产精品福利在线免费观看| 午夜老司机福利剧场| 男女那种视频在线观看| 成人av一区二区三区在线看| 男女之事视频高清在线观看| 国产大屁股一区二区在线视频| 一区二区三区激情视频| 99在线人妻在线中文字幕| bbb黄色大片| 亚洲无线观看免费| 欧美bdsm另类| 人妻少妇偷人精品九色| 日韩在线高清观看一区二区三区 | 亚洲va在线va天堂va国产| 国产成人a区在线观看| 成人av一区二区三区在线看| 久久6这里有精品| 99久久成人亚洲精品观看| 亚洲国产高清在线一区二区三| 男女边吃奶边做爰视频| 人妻久久中文字幕网| 全区人妻精品视频| 我要搜黄色片| 精品人妻一区二区三区麻豆 | 亚洲自拍偷在线| 亚洲国产精品成人综合色| 真实男女啪啪啪动态图| 村上凉子中文字幕在线| 久久久久国内视频| 国产午夜精品论理片| 国内揄拍国产精品人妻在线| 国产真实乱freesex| 国产亚洲精品综合一区在线观看| 草草在线视频免费看| 99久久精品国产国产毛片| 亚洲av中文字字幕乱码综合| 黄片wwwwww| 在线观看午夜福利视频| 18+在线观看网站| av国产免费在线观看| 97热精品久久久久久| 可以在线观看毛片的网站| 亚洲精品亚洲一区二区| 久久久久国内视频| 亚洲,欧美,日韩| 亚洲美女视频黄频| 无人区码免费观看不卡| 欧美不卡视频在线免费观看| 亚洲熟妇熟女久久| 一卡2卡三卡四卡精品乱码亚洲| 国产视频一区二区在线看| 国产精品三级大全| 精品久久久久久成人av| 国产三级中文精品| 亚洲精品久久国产高清桃花| 尾随美女入室| 五月玫瑰六月丁香| 欧美成人免费av一区二区三区| 国产老妇女一区| 成人特级av手机在线观看| 老司机福利观看| 成人特级黄色片久久久久久久| 久久久久久国产a免费观看| 精品日产1卡2卡| 精品一区二区三区人妻视频| 久久6这里有精品| 成人综合一区亚洲| 男人舔女人下体高潮全视频| 韩国av在线不卡| 在线a可以看的网站| 亚洲精品粉嫩美女一区| 又爽又黄a免费视频| 国内少妇人妻偷人精品xxx网站| 两个人视频免费观看高清| 欧美一级a爱片免费观看看| 亚洲精品一区av在线观看| 日本 欧美在线| 啦啦啦韩国在线观看视频| 国产一区二区三区av在线 | 少妇猛男粗大的猛烈进出视频 | 观看免费一级毛片| 国产不卡一卡二| 真人一进一出gif抽搐免费| 中国美白少妇内射xxxbb| 久久久久久久久大av| 久久精品国产亚洲网站| 精品一区二区三区视频在线观看免费| а√天堂www在线а√下载| 精品久久久久久久久久免费视频| 日本五十路高清| 欧美成人一区二区免费高清观看| 日韩欧美在线乱码| 免费不卡的大黄色大毛片视频在线观看 | 亚洲熟妇熟女久久| 天堂影院成人在线观看| 人人妻,人人澡人人爽秒播| 天天躁日日操中文字幕| 国产亚洲精品综合一区在线观看| 国产亚洲91精品色在线| 欧美成人免费av一区二区三区| 悠悠久久av| 亚洲七黄色美女视频| 亚洲欧美日韩高清在线视频| 国产 一区精品| 久久久国产成人免费| 国产精品乱码一区二三区的特点| 欧美性猛交黑人性爽| av黄色大香蕉| 婷婷色综合大香蕉| 麻豆精品久久久久久蜜桃| 国产精品亚洲一级av第二区| 久久6这里有精品| 亚洲av中文字字幕乱码综合| 天堂动漫精品| 欧洲精品卡2卡3卡4卡5卡区| 搡老岳熟女国产| 成人av在线播放网站| 日韩精品中文字幕看吧| 美女xxoo啪啪120秒动态图| 级片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 熟女人妻精品中文字幕| 999久久久精品免费观看国产| 一夜夜www| 日日夜夜操网爽| avwww免费| 国产中年淑女户外野战色| 在线a可以看的网站| 免费大片18禁| 成人欧美大片| 国产精品,欧美在线| 色视频www国产| 亚洲第一区二区三区不卡| 少妇裸体淫交视频免费看高清| 精品久久久久久久久久免费视频| 亚洲综合色惰| 九九爱精品视频在线观看| 日本五十路高清| 久久久色成人| 欧美日韩综合久久久久久 | 免费观看精品视频网站| 免费电影在线观看免费观看| 老司机午夜福利在线观看视频| 22中文网久久字幕| 中国美女看黄片| 熟妇人妻久久中文字幕3abv| 草草在线视频免费看| 久久99热6这里只有精品| 亚洲一区二区三区色噜噜| 久久6这里有精品| 午夜免费男女啪啪视频观看 | 国产精华一区二区三区| 国产精品1区2区在线观看.| 欧美一区二区国产精品久久精品| 精品国内亚洲2022精品成人| 久久这里只有精品中国| 国产精品久久久久久久久免| 亚洲精品一区av在线观看| 麻豆成人午夜福利视频| 色哟哟哟哟哟哟| 欧美xxxx黑人xx丫x性爽| 国产午夜精品论理片| 有码 亚洲区| 黄色丝袜av网址大全| 村上凉子中文字幕在线| 成年免费大片在线观看| 悠悠久久av| 高清毛片免费观看视频网站| 国产乱人伦免费视频| 日本免费a在线| 小蜜桃在线观看免费完整版高清| 国产国拍精品亚洲av在线观看| 亚洲精品乱码久久久v下载方式| 免费看a级黄色片| 精品一区二区三区视频在线| 全区人妻精品视频| 最后的刺客免费高清国语| 欧美激情久久久久久爽电影| 两个人视频免费观看高清| 国产精品人妻久久久久久| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 人人妻人人澡欧美一区二区| 非洲黑人性xxxx精品又粗又长| 久久精品91蜜桃| 日本 欧美在线| 欧美精品啪啪一区二区三区| 欧美日韩乱码在线| 国产免费男女视频| 国产黄色小视频在线观看| 国产伦精品一区二区三区四那| 久久久久久久久久黄片| 一本精品99久久精品77| 精品一区二区免费观看| 在线播放国产精品三级| 亚洲美女视频黄频| 色尼玛亚洲综合影院| 黄片wwwwww| 久久久国产成人精品二区| 一区二区三区高清视频在线| 别揉我奶头~嗯~啊~动态视频| 亚洲国产高清在线一区二区三| 性插视频无遮挡在线免费观看| 欧美+日韩+精品| 久久久精品欧美日韩精品| 免费观看精品视频网站| 国产黄a三级三级三级人| 亚洲人成网站高清观看| 国产精品免费一区二区三区在线| 九色国产91popny在线| 国产伦精品一区二区三区视频9| 精品日产1卡2卡| 国产一区二区三区av在线 | 在线免费十八禁| 成人特级黄色片久久久久久久| videossex国产| 国产白丝娇喘喷水9色精品| 我要搜黄色片| 亚洲最大成人中文| 最近中文字幕高清免费大全6 | 琪琪午夜伦伦电影理论片6080| 亚洲欧美日韩高清在线视频| 变态另类丝袜制服| 天堂av国产一区二区熟女人妻| 99在线人妻在线中文字幕| 精品一区二区三区av网在线观看| 长腿黑丝高跟| a级一级毛片免费在线观看| 日本精品一区二区三区蜜桃| 天堂网av新在线| 国产私拍福利视频在线观看| 国产精品一区二区免费欧美| 一进一出好大好爽视频| 精品久久久久久久久久久久久| 草草在线视频免费看| 蜜桃久久精品国产亚洲av| 国产 一区精品| 国产欧美日韩一区二区精品| 国产伦一二天堂av在线观看| 亚洲精品色激情综合| 国产v大片淫在线免费观看| 1000部很黄的大片| av视频在线观看入口| 中文资源天堂在线| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| a级毛片免费高清观看在线播放| 老师上课跳d突然被开到最大视频| 又粗又爽又猛毛片免费看| 毛片女人毛片| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区亚洲| 波多野结衣巨乳人妻| 成人特级av手机在线观看| 中文字幕高清在线视频| 嫩草影院精品99| 中文字幕高清在线视频| 国产精品女同一区二区软件 | 一区二区三区免费毛片| 欧美三级亚洲精品| 日韩亚洲欧美综合| 嫩草影视91久久| 非洲黑人性xxxx精品又粗又长| 久久欧美精品欧美久久欧美| 99热这里只有是精品在线观看| 欧美日本亚洲视频在线播放| 亚洲av不卡在线观看| 久久精品久久久久久噜噜老黄 | 免费高清视频大片| 成人精品一区二区免费| 99九九线精品视频在线观看视频| 国产探花极品一区二区| 女人十人毛片免费观看3o分钟| 黄色配什么色好看| 成人无遮挡网站| 国产老妇女一区| 春色校园在线视频观看| 国产精品一区二区性色av| 成年人黄色毛片网站| 国产女主播在线喷水免费视频网站 | 中国美女看黄片| 国产亚洲精品久久久久久毛片| 国产精华一区二区三区| 亚洲专区中文字幕在线| 欧美色欧美亚洲另类二区| 亚洲自拍偷在线| 精品午夜福利视频在线观看一区| 一区福利在线观看| 老师上课跳d突然被开到最大视频| 亚州av有码| 黄色欧美视频在线观看| 琪琪午夜伦伦电影理论片6080| 亚洲成av人片在线播放无| 在线观看美女被高潮喷水网站| 伦精品一区二区三区| 亚洲av熟女| av天堂中文字幕网| 国产成人影院久久av| av在线蜜桃| 亚洲成人久久爱视频| 亚洲精品亚洲一区二区| 老司机午夜福利在线观看视频| avwww免费| 亚洲图色成人| 午夜a级毛片| 日本撒尿小便嘘嘘汇集6| 变态另类成人亚洲欧美熟女| av在线天堂中文字幕| 看片在线看免费视频| 国产91精品成人一区二区三区| 国产成人a区在线观看| 女生性感内裤真人,穿戴方法视频| 免费av观看视频| 亚洲av第一区精品v没综合| 精品福利观看| 亚洲欧美日韩高清专用| 日韩精品青青久久久久久| 欧美3d第一页| 国内揄拍国产精品人妻在线| 天堂动漫精品| 欧美日韩亚洲国产一区二区在线观看| 99久久中文字幕三级久久日本| 啦啦啦观看免费观看视频高清| 精品久久久久久久久av| 亚洲无线观看免费| 亚洲五月天丁香| 国产亚洲精品久久久久久毛片| 亚洲国产色片| 免费看a级黄色片| 一卡2卡三卡四卡精品乱码亚洲| 99热网站在线观看| 免费观看精品视频网站| 免费大片18禁| 又紧又爽又黄一区二区| 精品人妻熟女av久视频| 搡老岳熟女国产| 尤物成人国产欧美一区二区三区| 免费高清视频大片| 国产伦一二天堂av在线观看| 一区二区三区免费毛片| 欧美一区二区精品小视频在线| 国产亚洲欧美98| 黄片wwwwww| 男女那种视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲91精品色在线| 国产伦一二天堂av在线观看| 国模一区二区三区四区视频| 在线观看av片永久免费下载| 老司机福利观看| 日韩欧美三级三区| netflix在线观看网站| 亚洲aⅴ乱码一区二区在线播放| 全区人妻精品视频| 久久精品国产亚洲网站| 草草在线视频免费看| 国内精品久久久久精免费| 亚洲性夜色夜夜综合| 亚洲欧美日韩东京热| videossex国产| 国产爱豆传媒在线观看| 亚洲精品456在线播放app | 中文字幕高清在线视频| 欧美在线一区亚洲| 国产午夜福利久久久久久| 亚洲欧美激情综合另类| 久久精品国产亚洲av香蕉五月| 国产aⅴ精品一区二区三区波| 国产国拍精品亚洲av在线观看| 国产午夜精品论理片| 国产精品久久视频播放| 久久精品人妻少妇| 国产精品久久电影中文字幕| 九色成人免费人妻av| 内射极品少妇av片p| 日韩高清综合在线| 此物有八面人人有两片| 亚洲图色成人| 欧美性猛交黑人性爽| 国产毛片a区久久久久| 午夜精品在线福利| 国产真实伦视频高清在线观看 | 女生性感内裤真人,穿戴方法视频| 一个人看视频在线观看www免费| 国产成人影院久久av| 久久欧美精品欧美久久欧美| 欧美国产日韩亚洲一区| 哪里可以看免费的av片| 国内精品一区二区在线观看| 成年人黄色毛片网站| 午夜福利18| 观看免费一级毛片| 国内毛片毛片毛片毛片毛片| 国产熟女欧美一区二区| 国内精品久久久久久久电影| 亚洲综合色惰| 欧美日本视频| 亚洲欧美清纯卡通| 亚洲美女黄片视频| 午夜久久久久精精品| 亚洲最大成人av| 国产黄片美女视频| bbb黄色大片| 男女下面进入的视频免费午夜| 亚洲三级黄色毛片| 一夜夜www| 精品一区二区三区人妻视频| 99久国产av精品| 久久久久久久久久黄片| 搞女人的毛片| 久久国内精品自在自线图片| 免费看a级黄色片| 观看免费一级毛片| 日韩欧美免费精品| 久久香蕉精品热| x7x7x7水蜜桃| 99久久中文字幕三级久久日本| 亚洲精华国产精华液的使用体验 | 日韩欧美免费精品| 国产午夜福利久久久久久| 欧美中文日本在线观看视频| 国产精品国产三级国产av玫瑰| 国产精品不卡视频一区二区| 一本久久中文字幕| 国产一区二区三区av在线 | 亚洲第一电影网av| 舔av片在线| 麻豆成人av在线观看| 久久久久性生活片| 免费不卡的大黄色大毛片视频在线观看 | 成人高潮视频无遮挡免费网站| 好男人在线观看高清免费视频| 五月伊人婷婷丁香| 日韩大尺度精品在线看网址| 欧美激情在线99| 日韩欧美 国产精品| 精品午夜福利视频在线观看一区| 有码 亚洲区| 黄色日韩在线| 非洲黑人性xxxx精品又粗又长| 久久久久久久久大av| 91狼人影院| 在现免费观看毛片| 搡老岳熟女国产| 欧美最黄视频在线播放免费| 毛片一级片免费看久久久久 | 国产淫片久久久久久久久| 亚洲精品日韩av片在线观看| 午夜精品一区二区三区免费看| 岛国在线免费视频观看| 日韩一本色道免费dvd| 国产久久久一区二区三区| 最近视频中文字幕2019在线8| 小说图片视频综合网站| 免费看av在线观看网站| 九色成人免费人妻av| 国产精品1区2区在线观看.| 日本撒尿小便嘘嘘汇集6| 88av欧美| 欧美精品啪啪一区二区三区| 亚洲人与动物交配视频| 日日撸夜夜添| 国产探花极品一区二区| 熟女人妻精品中文字幕| 高清日韩中文字幕在线| 欧美+日韩+精品| 简卡轻食公司| 99热网站在线观看| 一个人看视频在线观看www免费| 99精品久久久久人妻精品| 日本 av在线| 中国美白少妇内射xxxbb| 久久久色成人| 男女做爰动态图高潮gif福利片| 尾随美女入室| 午夜精品久久久久久毛片777| 九九久久精品国产亚洲av麻豆| 在线播放无遮挡| 欧美色视频一区免费| 久久婷婷人人爽人人干人人爱| 国产精品不卡视频一区二区| 久久久色成人| 久久99热这里只有精品18| 伦精品一区二区三区| 国产精品一区www在线观看 | 看片在线看免费视频| 亚洲va在线va天堂va国产| 婷婷色综合大香蕉| 日韩欧美在线乱码| 男人和女人高潮做爰伦理| 亚洲美女视频黄频| 伦精品一区二区三区| 国产美女午夜福利| 高清日韩中文字幕在线|