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

    1990-2021年?yáng)|太湖網(wǎng)圍養(yǎng)殖規(guī)模變化下的水質(zhì)變化特征及成因分析*

    2023-01-13 09:42:44李一凡谷孝鴻曾慶飛賈冰嬋訾鑫源陳輝輝毛志剛
    湖泊科學(xué) 2023年1期
    關(guān)鍵詞:沉水植物太湖水體

    李一凡,谷孝鴻,曾慶飛,賈冰嬋,訾鑫源,陳輝輝,毛志剛,葛 優(yōu)

    (1:中國(guó)科學(xué)院南京地理與湖泊研究所湖泊與環(huán)境國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京 210008) (2:中國(guó)科學(xué)院大學(xué),北京 100049)

    我國(guó)的湖庫(kù)網(wǎng)圍養(yǎng)殖始于1970s末。人們利用大水面自然資源,發(fā)展高效利用的漁業(yè)增養(yǎng)殖模式,為湖區(qū)增效和漁民增收致富做出了重要貢獻(xiàn)。然而,隨著集約化網(wǎng)圍養(yǎng)殖規(guī)模和密度的增加,加上流域工農(nóng)業(yè)污染的影響,水體富營(yíng)養(yǎng)化和生態(tài)退化現(xiàn)象日益嚴(yán)重。為了保護(hù)湖庫(kù)生態(tài)環(huán)境和漁業(yè)資源,各級(jí)部門(mén)有計(jì)劃地開(kāi)展了養(yǎng)殖網(wǎng)圍整治和拆除工作。東太湖是太湖東部一個(gè)典型的草型淺水湖灣,因其適宜的光熱、自然地理和水體條件發(fā)展成為我國(guó)最早開(kāi)發(fā)漁業(yè)養(yǎng)殖的湖泊之一。東太湖網(wǎng)圍養(yǎng)殖結(jié)構(gòu)和養(yǎng)殖規(guī)模經(jīng)歷了幾次變動(dòng)。首先是“以湖養(yǎng)湖”的高密度養(yǎng)魚(yú)農(nóng)業(yè)試驗(yàn)研究成功,實(shí)現(xiàn)了從粗放粗養(yǎng)的常規(guī)養(yǎng)殖模式到湖區(qū)網(wǎng)圍精養(yǎng)的轉(zhuǎn)變[1]。1990s末,隨著人們對(duì)名特優(yōu)水產(chǎn)品的大量需求,圍網(wǎng)養(yǎng)殖結(jié)構(gòu)調(diào)整,逐漸轉(zhuǎn)變?yōu)橐灾腥A絨螯蟹(簡(jiǎn)稱(chēng)河蟹)(Eriocheirsinensis)為主養(yǎng)、混養(yǎng)魚(yú)蝦的生態(tài)養(yǎng)殖模式。2000-2008年是網(wǎng)圍養(yǎng)殖發(fā)展的鼎盛時(shí)期,全湖河蟹水產(chǎn)養(yǎng)殖面積達(dá)118.56 km2。2008年底,東太湖實(shí)施了退圍還湖網(wǎng)圍整治工程,將網(wǎng)圍養(yǎng)殖水面壓縮至30.00 km2。直到2018年末,伴隨著最后一季河蟹養(yǎng)殖結(jié)束,東太湖開(kāi)始了網(wǎng)圍全面清除工作,結(jié)束了長(zhǎng)達(dá)34年的太湖網(wǎng)圍養(yǎng)殖歷史。

    按照中央環(huán)保督察要求,太湖、滆湖、長(zhǎng)蕩湖等全面拆除了湖泊養(yǎng)殖網(wǎng)圍。截至2019年,全國(guó)網(wǎng)圍養(yǎng)殖面積銳減到1.17×102km2[10]。修復(fù)退圍還湖區(qū)生態(tài)系統(tǒng)、改善水體生態(tài)環(huán)境和保護(hù)生物多樣性是湖泊治理規(guī)劃的關(guān)鍵。系統(tǒng)回顧網(wǎng)圍養(yǎng)殖水環(huán)境的長(zhǎng)期變化特征,探尋其驅(qū)動(dòng)因素、存在問(wèn)題和發(fā)展趨勢(shì),將為我國(guó)湖泊“十四五”生態(tài)修復(fù)與保護(hù)規(guī)劃提供數(shù)據(jù)支撐和借鑒。東太湖是我國(guó)網(wǎng)圍養(yǎng)殖發(fā)展最成熟、整治管理力度最大的水域,本文通過(guò)系統(tǒng)分析1990-2021年水體營(yíng)養(yǎng)鹽在不同湖區(qū)和不同網(wǎng)圍發(fā)展階段的變化規(guī)律,解析影響營(yíng)養(yǎng)鹽變動(dòng)的因素、存在問(wèn)題和發(fā)展趨勢(shì),為東太湖生態(tài)修復(fù)與保護(hù)以及全國(guó)退圍還湖區(qū)生態(tài)環(huán)境治理提供數(shù)據(jù)基礎(chǔ)和理論支持。

    1 材料與方法

    1.1 數(shù)據(jù)來(lái)源及點(diǎn)位設(shè)置

    東太湖是太湖東南部的湖灣,面積131 km2,平均水深1.3 m,是太湖主要的出水區(qū)域,承擔(dān)了重要的航運(yùn)作用,也是東部下游地區(qū)上海市等地的主要水源地[4]。本文采用的數(shù)據(jù)分為兩部分,其中1990-2005年數(shù)據(jù)來(lái)自中國(guó)科學(xué)院南京地理與湖泊研究所楊龍?jiān)毖芯繂T提供的歷史監(jiān)測(cè)數(shù)據(jù);2006-2021年數(shù)據(jù)來(lái)自本課題組長(zhǎng)期監(jiān)測(cè),監(jiān)測(cè)點(diǎn)位如圖 1所示,數(shù)據(jù)采集和分析方法均參照《湖泊富營(yíng)養(yǎng)化調(diào)查規(guī)范》(第二版)[11]的規(guī)定執(zhí)行。14個(gè)監(jiān)測(cè)點(diǎn)位包含原網(wǎng)圍區(qū)(1#、7#、13#、14#)、河口區(qū)(2#、9#、11#、12#)、湖心區(qū)(6#、8#)、瓜涇口(10#)、東茭咀(4#)和太浦河口(5#)。整個(gè)太湖的數(shù)據(jù)來(lái)自中國(guó)科學(xué)院太湖湖泊生態(tài)系統(tǒng)研究站點(diǎn)位監(jiān)測(cè)數(shù)據(jù)。

    圖1 東太湖水質(zhì)監(jiān)測(cè)調(diào)查點(diǎn)位分布(網(wǎng)格代表 2009-2018年網(wǎng)圍整治階段網(wǎng)圍分布區(qū))Fig.1 Distribution of sampling sites in East Lake Taihu (The grid represents the net-pen distribution from 2009 to 2018)

    1.2 水質(zhì)指標(biāo)的測(cè)定

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

    2006-2021年?yáng)|太湖水體TN、TP、Chl.a和CODMn濃度的年際變化數(shù)據(jù)為該年度所有點(diǎn)位和所有采樣時(shí)間的平均值±標(biāo)準(zhǔn)誤,采用R x64 4.1.1軟件中g(shù)gplot2包繪制1990-2021年水質(zhì)指標(biāo)濃度變化柱狀圖;采用AutoCAD 2010軟件繪制東太湖網(wǎng)圍規(guī)模變化的三階段水體理化指標(biāo)空間變化(2008年前網(wǎng)圍養(yǎng)殖階段、2009-2018年網(wǎng)圍整治階段和2019年后網(wǎng)圍拆除階段),數(shù)據(jù)來(lái)自該階段該點(diǎn)位所有采樣時(shí)間的平均值±標(biāo)準(zhǔn)誤,采用單因素方差分析(one-way analysis of variance)分析不同點(diǎn)位的組內(nèi)數(shù)據(jù)的顯著性差異;采用R x64 4.1.1軟件中g(shù)gplot2包對(duì)網(wǎng)圍拆除后水質(zhì)指標(biāo)的季節(jié)變化進(jìn)行分析,以保留95%的置信區(qū)間對(duì)14個(gè)采樣點(diǎn)位的數(shù)據(jù)擬合,可視化營(yíng)養(yǎng)鹽濃度的季節(jié)變化趨勢(shì);為探究東太湖30年來(lái)水體中主要污染因子和太浦河口出水量對(duì)東太湖水質(zhì)的影響,采用R x64 4.1.1軟件中g(shù)gbiplot包和factoextra包,將太浦河口出水量分別與1990-2020年?yáng)|太湖全域水質(zhì)和2008-2020年太浦河口與東茭咀兩區(qū)域水質(zhì)進(jìn)行主成分分析(principal component analysis),利用ggcorrplot包展示各指標(biāo)之間的Pearson相關(guān)系數(shù)和顯著性P值。

    選取SD、Chl.a、TN、TP和CODMn這5個(gè)與湖泊富營(yíng)養(yǎng)化狀態(tài)變化最為密切的參數(shù)開(kāi)展水體富營(yíng)養(yǎng)化狀態(tài)綜合評(píng)價(jià)。計(jì)算公式和判定標(biāo)準(zhǔn)參照《湖泊富營(yíng)養(yǎng)化調(diào)查規(guī)范》(第二版)[11]。

    2 結(jié)果與分析

    2.1 東太湖1990-2021年水質(zhì)年際變化

    東太湖水體TN濃度在1990年為0.53 mg/L,隨著網(wǎng)圍規(guī)模擴(kuò)大,到1999年TN濃度增加到1.01 mg/L。隨著養(yǎng)殖結(jié)構(gòu)改變,TN濃度有所下降,至2002年下降到0.47 mg/L。但是隨著養(yǎng)殖強(qiáng)度增加,TN濃度呈現(xiàn)增加趨勢(shì),至2008年網(wǎng)圍整治之前達(dá)到1.3 mg/L。2009-2018年間,TN濃度基本維持穩(wěn)定,平均濃度為1.33 mg/L。網(wǎng)圍拆除后,TN濃度在2020年開(kāi)始有所反彈,至2021年上升至1.87 mg/L(圖 2A)。

    1990-2021年間東太湖水體TP濃度變化趨勢(shì)和TN相似,隨著網(wǎng)圍養(yǎng)殖結(jié)構(gòu)和規(guī)模變化。1995年以前, TP平均濃度為0.015 mg/L,滿(mǎn)足地表水Ⅱ類(lèi)水質(zhì)標(biāo)準(zhǔn)。隨著網(wǎng)圍養(yǎng)殖的快速發(fā)展,TP濃度呈明顯上升趨勢(shì),到2007年網(wǎng)圍整治前達(dá)到0.06 mg/L,隨后有所降低,到 2018年網(wǎng)圍拆除前下降到0.051 mg/L。但是,近年來(lái)有上升趨勢(shì),2021年4月,TP最高達(dá)到 0.128 mg/L(圖 2B)。

    東太湖水體的CODMn濃度近30年的變化較為穩(wěn)定,平均值為4.54 mg/L,達(dá)到地表水Ⅲ類(lèi)水質(zhì)標(biāo)準(zhǔn)。2008年網(wǎng)圍整治和2018年的網(wǎng)圍全部拆除,CODMn濃度在一定時(shí)期均呈現(xiàn)了下降趨勢(shì)。2018年網(wǎng)圍拆除前CODMn濃度達(dá)到最高(6.08 mg/L)。網(wǎng)圍拆除后近兩年CODMn濃度有上升趨勢(shì),2021年達(dá)到5.72 mg/L(圖 2C)。

    東太湖水體Chl.a濃度在1990-2007年間趨于穩(wěn)定,均低于10 μg/L。2008年網(wǎng)圍整治后,2008-2010年Chl.a濃度較高,至2010年達(dá)到22.75 μg/L。2010年后呈下降趨勢(shì),至2015年降至7.21 μg/L。2016年后Chl.a濃度呈上升趨勢(shì),至2021年最高可達(dá)28.09μg/L(圖2D)。

    圖2 東太湖1990-2021年水體營(yíng)養(yǎng)鹽濃度長(zhǎng)期變化趨勢(shì)(2008年:網(wǎng)圍整治;2018年:網(wǎng)圍拆除)Fig.2 Variational trends of nutrient concentration in water of East Lake Taihu from 1990 to 2021 (2008: Net-pen culture governance, 2018: Net-pen culture demolition)

    2.2 東太湖不同水域水質(zhì)演變

    東太湖水質(zhì)在不同水域不同網(wǎng)圍整治時(shí)期變化情況如圖 3所示。東太湖原網(wǎng)圍區(qū)的水深和透明度近30年來(lái)基本維持穩(wěn)定,分別在1.35 m(P>0.05)和0.80 m(P>0.05)左右波動(dòng);湖心區(qū)的水深略高于原網(wǎng)圍區(qū),透明度近年來(lái)略有下降,由2008年之前的0.87 m下降到網(wǎng)圍拆除后的0.51 m;北部河口張家浜水體透明度近年來(lái)有所提升,網(wǎng)圍拆除后達(dá)到0.68 m;太浦河口和北部瓜涇口水深自2018年后增加明顯,分別由2018年的2.51和1.28 m提高到3.05(F=5.07,P<0.05)和3.55 m(F=79.28,P<0.01),透明度分別由2018年的0.44和0.56 m降低到網(wǎng)圍拆除后的0.29和0.46 m。東太湖入湖口的東茭咀透明度下降嚴(yán)重,從1.18 m降至0.29 m(F=15.70,P<0.01)(圖 3A)。

    原網(wǎng)圍區(qū)和湖心區(qū)水體的TN濃度基本穩(wěn)定,維持在1.10 mg/L(P>0.05)左右。位于東太湖北部的張家浜與瓜涇口的TN濃度,自2008年之后呈下降趨勢(shì),分別由2008年前的1.74和1.76 mg/L降至2019年后的1.26和1.00 mg/L。南部的東茭咀和太浦河口TN濃度上升明顯,分別由2008年前的0.68和0.75 mg/L升至2019年后的2.00 (F=7.62,P<0.01)和1.67 mg/L(F=4.44,P<0.05)。3階段內(nèi)的NH3-N濃度除了東茭咀一直呈現(xiàn)上升趨勢(shì)(F=5.43,P<0.05),其它區(qū)域均表現(xiàn)為先降低又升高,2019年之后,NH3-N濃度維持在0.22~0.32 mg/L之間 (圖3B)。

    圖3 東太湖不同網(wǎng)圍規(guī)模不同區(qū)域水質(zhì)演變 (階段1為2008年前,階段2為2009-2018年,階段3為2019年后)Fig.3 Water quality changes in different regions of East Lake Taihu

    東太湖TOC濃度在各點(diǎn)呈現(xiàn)先升高再降低的趨勢(shì),2009-2018年網(wǎng)圍整治階段TOC濃度最高,平均達(dá)到8.12 mg/L,網(wǎng)圍拆除后降低為6.00 mg/L,點(diǎn)位間差異不顯著(P>0.05)。東茭咀和太浦河口水域CODMn濃度近30年呈逐漸上升的趨勢(shì),網(wǎng)圍拆除后分別達(dá)到4.73和4.70 mg/L。瓜涇口水域的TOC和CODMn濃度近30年逐漸下降(圖3D)。

    原網(wǎng)圍區(qū)和湖心區(qū)的Chl.a濃度在2008年網(wǎng)圍整治后分別從17.92和20.77 μg/L上升至22.35和29.52 μg/L,該區(qū)域的SS濃度近30年基本維持在16 mg/L(P>0.05)左右。張家浜和瓜涇口的SS濃度在網(wǎng)圍拆除后有所下降。東茭咀和太浦河口的Chl.a濃度近30年內(nèi)均呈明顯的上升趨勢(shì),Chl.a濃度分別從2008年前的4.18和5.42 μg/L上升至2019年后的24.95 (F=14.16,P<0.01)和34.10 μg/L(F=5.96,P<0.05),東茭咀的SS濃度從1.49 mg/L顯著上升至38.38 mg/L(F=9.32,P<0.01)(圖3E)。

    2.3 東太湖網(wǎng)圍拆除后水質(zhì)季節(jié)變化

    網(wǎng)圍拆除后東太湖水體中TN和TDN濃度季節(jié)變化規(guī)律相似, 最低值均出現(xiàn)在7-10月,平均值分別為0.89和0.63 mg/L;10月至次年1月濃度較高,平均值分別為1.57和0.95 mg/L(圖4A,C)。

    圖4 東太湖網(wǎng)圍拆除后氮、磷濃度的季節(jié)變化Fig.4 Seasonal variations of nitrogen and phosphorus concentrations in East Lake Taihu after the complete renovation of net-pen

    網(wǎng)圍拆除后東太湖水體中TP和TDP的平均濃度除了2021年4月份,基本維持穩(wěn)定,分別為0.05 (P>0.05)和0.02 mg/L(P>0.05)左右。2020年和2021年1月相對(duì)較高,TP平均可達(dá)0.10 mg/L。2021年4月因采樣當(dāng)天風(fēng)浪較大,底泥再懸浮,TP(F=15.56,P<0.01)和TDP(F=133.98,P<0.01)濃度顯著升高(圖4B,D)。

    東太湖網(wǎng)圍拆除后,水體SS、CODMn、Chl.a和TOC在2019年呈現(xiàn)上升趨勢(shì),直到2020年1月分別達(dá)到38.33 mg/L、6.91 mg/L、35.30 μg/L和4.24 mg/L,沒(méi)有出現(xiàn)季節(jié)波動(dòng)。但2020-2021年,逐漸呈現(xiàn)夏秋季低、冬春季高的趨勢(shì), 2020年7-10月最低,SS、CODMn、Chl.a和TOC分別為13.29 mg/L、4.32 mg/L、14.65 μg/L和3.57 mg/L(圖5)。

    圖5 東太湖網(wǎng)圍完全拆除后SS、Chl.a、CODMn和TOC濃度變化特征Fig.5 Variations of SS, Chl.a, CODMn and DOC concentrations in East Lake Taihu after the complete renovation of net-pen

    2.4 東太湖30年不同階段水質(zhì)主成分和相關(guān)性分析

    TN是主成分分析中貢獻(xiàn)度最大的水質(zhì)指標(biāo)(圖6),表明TN是1990-2020年間東太湖主要的污染因子,TN、TP、CODMn和Chl.a濃度之間具有顯著的正相關(guān)性(P<0.05)。30年間網(wǎng)圍面積和太浦河口出水量與東太湖水體營(yíng)養(yǎng)鹽濃度變化相關(guān)性不顯著(P>0.05)。

    圖6 1990-2020年?yáng)|太湖水質(zhì)主成分分析圖和相關(guān)性熱力圖Fig.6 Principal component analysis diagram and correlation thermal map of East Lake Taihu from 1990 to 2020

    進(jìn)一步發(fā)現(xiàn)太浦河口與東茭咀水域的營(yíng)養(yǎng)鹽濃度自2008年后顯著上升,且2005年起“引江濟(jì)太”工程進(jìn)入長(zhǎng)效運(yùn)行,太浦河口出湖水量開(kāi)始增加。因此,對(duì)2008年網(wǎng)圍整治工程實(shí)施后的水質(zhì)指標(biāo)與出湖水量和降雨量進(jìn)行主成分分析,可知TN依舊是太浦河口與東茭咀的主要污染因子,TN濃度與透明度和Chl.a濃度具有顯著相關(guān)性(P<0.01)。降雨量與各水質(zhì)指標(biāo)之間并無(wú)顯著相關(guān)性(P>0.05),太浦河口出水量與Chl.a濃度有顯著正相關(guān)關(guān)系(P<0.05)(圖7)。

    圖7 2008-2020年?yáng)|太湖水質(zhì)主成分分析圖和相關(guān)性熱力圖Fig.7 Principal component analysis diagram and correlation thermal map of East Lake Taihu from 2008 to 2020

    3 討論

    3.1 網(wǎng)圍養(yǎng)殖結(jié)構(gòu)調(diào)整和規(guī)模變動(dòng)對(duì)東太湖水質(zhì)年際波動(dòng)的影響

    1990-2021年期間,東太湖TN、TP、CODMn和Chl.a濃度平均值存在年際波動(dòng),與網(wǎng)圍養(yǎng)殖結(jié)構(gòu)調(diào)整和規(guī)模變動(dòng)密切相關(guān)。東太湖是典型的草型湖泊,1980s以苦草(Vallisnerianatans(Lour.) Hara)、馬來(lái)眼子菜(PotamogetonwrightiiMorong)、微齒眼子菜(PotamogetonmaackianusA. Bennett)和輪葉黑藻(Hydrillaverticillata)占優(yōu)勢(shì)。中國(guó)科學(xué)院南京地理與湖泊研究所的一批科技人員開(kāi)創(chuàng)了我國(guó)大水面水草資源高效利用的網(wǎng)圍養(yǎng)殖模式,以草魚(yú)(Ctenopharyngodonidella)、團(tuán)頭魴(Megalobramaamblycephala)、鯉(Cyprinuscarpio)、鯽(Carassiusauratus)、青魚(yú)(Mylopharyngodonpiceus)為主養(yǎng),搭配鰱(Hypophthalmichthysmolitrix)、鳙(Aristichthysnobilis)、鱖(Sinipercachuatsi)和翹嘴紅鲌(Erythroculterilishaeformis)。到1997年,全湖網(wǎng)圍養(yǎng)魚(yú)面積達(dá)60 km2(圖 8)。隨著網(wǎng)圍養(yǎng)魚(yú)強(qiáng)度增加,水體富營(yíng)養(yǎng)化現(xiàn)象日益凸顯。據(jù)報(bào)道,水產(chǎn)品只能攝食約70%的餌料[12],多余的餌料、養(yǎng)殖藥物殘留和水產(chǎn)品的排泄物一同沉入湖底[13]。微生物經(jīng)過(guò)自身生理作用分解釋放出大量的氮、磷到水體中,增加沉積物污染風(fēng)險(xiǎn),加速湖泊富營(yíng)養(yǎng)化[14-16]。根據(jù)楊清心等的研究[17],東太湖網(wǎng)圍養(yǎng)魚(yú)產(chǎn)量為7500~11250 kg/hm2時(shí),相當(dāng)于向湖中投入超過(guò)產(chǎn)量60%的氮、磷,即每生產(chǎn)1 t的魚(yú)就會(huì)向湖內(nèi)釋放141.25 kg的氮和14.14 kg 的磷。

    圖8 東太湖網(wǎng)圍完全拆除前(1984- 2019年)的漁業(yè)養(yǎng)殖面積和河蟹養(yǎng)殖比例Fig.8 Fisheries farming area and the ratio of Chinese mitten crab farming of East Lake Taihu from 1984 to 2019

    到1990s末網(wǎng)圍養(yǎng)魚(yú)鼎盛時(shí)期,東太湖水體TN、TP和CODMn濃度達(dá)到了一個(gè)小高峰(圖 2)。因東太湖還保持較高的水草覆蓋度,Chl.a濃度變化不明顯,依舊維持較低的水平。隨后,東太湖網(wǎng)圍養(yǎng)殖模式逐漸轉(zhuǎn)向了高經(jīng)濟(jì)效益的河蟹為主,同時(shí)套養(yǎng)魚(yú)蝦的生態(tài)混養(yǎng)模式,最大程度地利用河蟹餌料,加上沉水植物伊樂(lè)藻(ElodeacanadensisMichx)和輪葉黑藻的有效管控和螺類(lèi)投放,在提高綜合效益的同時(shí)盡可能減湖泊外來(lái)物質(zhì)輸入對(duì)水質(zhì)的影響[18]。直到2000年初,水體營(yíng)養(yǎng)鹽濃度呈下降趨勢(shì)。

    隨著河蟹養(yǎng)殖強(qiáng)度增大,2008年全湖養(yǎng)殖面積達(dá)到118.56 km2,占全湖水面的90.50%。近乎全湖的網(wǎng)圍覆蓋,嚴(yán)重影響了水生植物的收割和氮、磷營(yíng)養(yǎng)鹽的輸出,從而導(dǎo)致湖泊的自?xún)裟芰ο陆?,加重了東太湖水質(zhì)的惡化[19]。Luo等也發(fā)現(xiàn)陽(yáng)澄湖中的氮、磷等營(yíng)養(yǎng)鹽濃度與網(wǎng)圍面積的大小呈正相關(guān)[20]。因此,2008年網(wǎng)圍整治前,營(yíng)養(yǎng)鹽水平達(dá)到第2個(gè)小高峰。隨后,網(wǎng)圍養(yǎng)殖整治將面積縮減為30 km2,東太湖的水質(zhì)狀況得到了改善。直到2018年養(yǎng)殖網(wǎng)圍全部拆除,2019年有個(gè)平緩的響應(yīng)期,近兩年水質(zhì)狀況不容樂(lè)觀,2020年?yáng)|太湖的綜合營(yíng)養(yǎng)狀態(tài)指數(shù)達(dá)到69,接近整個(gè)太湖水平(圖9)。

    圖9 2002-2020年太湖和東太湖 綜合營(yíng)養(yǎng)狀態(tài)指數(shù)變化趨勢(shì)Fig.9 Variations of comprehensive trophic level index of Lake Taihu and East Lake Taihu from 2002 to 2020

    3.2 網(wǎng)圍拆除后東太湖水質(zhì)變化原因分析

    影響東太湖水質(zhì)狀況的可能因素有流域河流排污和土地利用、湖體水生植物現(xiàn)存量和群落結(jié)構(gòu)改變、水文水動(dòng)力、與整個(gè)太湖的連通影響以及氣候變化和極端天氣等。

    圖10 2002-2020年太湖流域 降雨量(數(shù)據(jù)來(lái)源:水利部太湖流域管理局)Fig.10 Rainfall in Lake Taihu Basin from 1999 to 2020 (Data source: Taihu Basin Authority)

    圖11 2005-2020年太湖出入 湖水量(數(shù)據(jù)來(lái)源:太湖流域管理局)Fig.11 Volume of water input and output of Lake Taihu from 2005 to 2020 (Data source: Taihu Basin Administration Bureau)

    2)水生植物現(xiàn)存量和群落結(jié)構(gòu)改變的影響。水生植被可以固定水體中的氮、磷,并通過(guò)抑制湖泊底泥再懸浮控制沉積物中營(yíng)養(yǎng)鹽的釋放,改善水質(zhì)[24]。合理的水生植被群落結(jié)構(gòu)組成和生物量是維持健康湖泊生態(tài)系統(tǒng)的關(guān)鍵[25]。1960s,東太湖以馬來(lái)眼子菜為主,占全湖面積的90%。2002年發(fā)展為以伊樂(lè)藻+微齒眼子菜+輪葉黑藻+苦草群叢為主,共發(fā)現(xiàn)沉水植物17種,浮葉植物13種和漂浮植物9種[26]。2019年以后,由于缺少漁民對(duì)沉水植物伊樂(lè)藻、輪葉黑藻、苦草的種植和對(duì)浮葉植物的打撈,沉水植物演替為以狐尾藻(MyriophyllumverticillatumL.)、金魚(yú)藻(CeratophyllumdemersumL.)和菹草(PotamogetoncrispusL.)占優(yōu)勢(shì)。2021年沉水植物面積銳減到13.17 km2,浮葉植物野菱(Trapaincisavar.sieb.)瘋狂擴(kuò)張,達(dá)到34.97 km2,水生植物總面積占全湖總面積的45%,主要分布于原網(wǎng)圍區(qū)(表 1)。當(dāng)外界的干擾超過(guò)其耐受限度時(shí),沉水植物面積萎縮,群落結(jié)構(gòu)趨于單一化,或者消失[27-28]。當(dāng)水生植物群落結(jié)構(gòu)單一化時(shí),易形成“瘋長(zhǎng)”,腐敗壞水[29]。且浮葉植物較沉水植物有較強(qiáng)的抗風(fēng)浪能力,又缺少漁民打撈,因此網(wǎng)圍拆除后浮葉植物占優(yōu)勢(shì)[8]。傅玲研究發(fā)現(xiàn)不同生態(tài)型組合的植物群叢對(duì)氮、磷去除效果優(yōu)于單一植物群叢[30]。原網(wǎng)圍區(qū)因夏秋季金魚(yú)藻、狐尾藻、野菱生物量高于其它區(qū)域,水質(zhì)優(yōu)于其它水域,和網(wǎng)圍拆除前差異不大;但因缺少人工管護(hù),秋、冬季腐敗壞水風(fēng)險(xiǎn)加大,特別是對(duì)于生物量較高的漂浮植物野菱。沉水植物多樣性的提高及群落結(jié)構(gòu)的優(yōu)化對(duì)水生態(tài)系統(tǒng)的健康與穩(wěn)定具有重要意義。但整體來(lái)說(shuō),東太湖湖心、原網(wǎng)圍區(qū)、張家浜和瓜涇口水體營(yíng)養(yǎng)鹽相對(duì)于網(wǎng)圍拆除前變化不大, Chl.a濃度上升。水體營(yíng)養(yǎng)鹽上升最高的是太浦河口和東茭咀,整體拉升了網(wǎng)圍拆除后東太湖氮、磷營(yíng)養(yǎng)鹽的平均水平。

    表1 2017-2021年?yáng)|太湖水生植被類(lèi)群面積Tab.1 Area changes of different aquatic vegetations in East Lake Taihu from 2017 to 2021

    3)東太湖水文水動(dòng)力變化和與大太湖的連通影響。本課題組的研究發(fā)現(xiàn),網(wǎng)圍存在對(duì)表層和中層的湖流衰減影響最大,可衰減73%~96%;圍網(wǎng)拆除后,水動(dòng)力條件改變(水流運(yùn)動(dòng)對(duì)沉水植物的拉伸、攪動(dòng)、拖曳作用)對(duì)沉水植物生物量以及群落組成都有明顯影響[31-32]。狐尾藻的抗剪切力強(qiáng)于伊樂(lè)藻和輪葉黑藻,這與東太湖現(xiàn)階段沉水植物以狐尾藻占優(yōu)勢(shì)有一定關(guān)系。

    自2015年開(kāi)始,東茭咀和太浦河口處的TN、TP、Chl.a和SS濃度均出現(xiàn)明顯的上升趨勢(shì),透明度也下降到現(xiàn)在的0.30 cm左右,與整個(gè)太湖接近。楊井志成等[8]發(fā)現(xiàn)2018網(wǎng)圍拆除后,在太浦河口附近水域監(jiān)測(cè)到藍(lán)藻水華,面積為8.79 km2,占比為4.31%。研究認(rèn)為是圍網(wǎng)拆除后,水動(dòng)力增強(qiáng),原圍網(wǎng)區(qū)底泥沉積物釋放至太浦閘口附近水域,加重了該片水域的富營(yíng)養(yǎng)化程度。在本研究中,太浦河口和東茭咀水域的Chl.a濃度同太浦河口出湖水量呈顯著正相關(guān)(P<0.05)。郝文彬等[33]運(yùn)用環(huán)境流體動(dòng)力學(xué)模型也發(fā)現(xiàn),加大望虞河入湖流量可顯著改變東部湖區(qū)的水循環(huán)條件。2015和2016 年,持續(xù)兩年的特大洪水(圖,加上整個(gè)太湖東部人工水生植物刈割打撈,東茭咀附近的水生植物急劇減少[34]。網(wǎng)圍和水生植物具有降低湖水流速,阻礙沉積物中營(yíng)養(yǎng)鹽再懸浮的作用[35]。網(wǎng)圍拆除、經(jīng)由太浦河等出湖口的出湖水量增加(圖11)、水生植物消失,東茭咀附近水動(dòng)力作用加強(qiáng),加上原網(wǎng)圍區(qū)底泥營(yíng)養(yǎng)鹽釋放,綜合因素引起東太湖南部區(qū)域營(yíng)養(yǎng)鹽濃度大幅上升。在朱偉等[36]的研究中也發(fā)現(xiàn)入湖水量的增加開(kāi)始影響太湖存在的“西濁東清”水質(zhì)結(jié)構(gòu),水體交換周期縮短會(huì)使東西太湖水質(zhì)出現(xiàn)均化的現(xiàn)象,東部太湖水質(zhì)會(huì)出現(xiàn)下降的趨勢(shì)。

    4)極端天氣影響。2016年的特大洪水造成太湖出現(xiàn)磷反彈,洪水造成的年度性水體交換周期縮短已經(jīng)對(duì)東部太湖水質(zhì)產(chǎn)生了不利的影響[36],認(rèn)為2016年突然出現(xiàn)的磷反彈與這一年發(fā)生的特大洪水一次性多增加579 t磷通量有著顯著的關(guān)系。在本研究中,東太湖TN、TP、CODMn和Chl.a濃度在2016和2017年有較大提高可能與強(qiáng)降雨有關(guān),雖然相關(guān)性不顯著。

    3.3 東太湖網(wǎng)圍全拆除后的生態(tài)修復(fù)措施建議

    東太湖退圍還湖是太湖受損生態(tài)系統(tǒng)局部水域水質(zhì)改善、生態(tài)修復(fù)和多樣性保護(hù)的關(guān)鍵,也是太湖生物資源養(yǎng)護(hù)及漁業(yè)高質(zhì)量可持續(xù)發(fā)展和下游區(qū)域優(yōu)質(zhì)水源保障供給的重要基礎(chǔ)。東太湖網(wǎng)圍養(yǎng)殖拆除,目標(biāo)是重構(gòu)良好的生態(tài)系統(tǒng),但3年來(lái)的環(huán)境演變并沒(méi)有顯現(xiàn)出最初設(shè)想的效果,需要加強(qiáng)對(duì)東太湖的管理,并對(duì)近階段東太湖優(yōu)質(zhì)水生植被退化的現(xiàn)狀進(jìn)行人工干預(yù),確保恢復(fù)健康自然的草型生態(tài)系統(tǒng)格局。

    3.3.1 應(yīng)及早開(kāi)展水生植物資源調(diào)控強(qiáng)化水質(zhì)凈化 針對(duì)當(dāng)前東太湖水生植被浮葉植物菱占優(yōu)勢(shì)、覆蓋度高,沉水植物狐尾藻和金魚(yú)藻占優(yōu)勢(shì)、不易被魚(yú)類(lèi)等利用等導(dǎo)致東太湖生態(tài)系統(tǒng)服務(wù)功能降低的問(wèn)題,建議對(duì)東太湖不同區(qū)域水生植被進(jìn)行調(diào)控,組織精準(zhǔn)刈割,適度移出過(guò)量浮葉植物和沉水植物,收割控制沿湖岸蘆葦?shù)取8鶕?jù)多樣性-穩(wěn)定性理論,在目前沒(méi)有水生植物分布區(qū)域,開(kāi)展微地形地貌生態(tài)改造提高透明度恢復(fù)沉水植物,在浮葉植物優(yōu)勢(shì)區(qū)開(kāi)展風(fēng)浪消減-浮葉植物塊狀刈割-沉水植物擴(kuò)繁,在沉水植物單一區(qū)開(kāi)展沉水植物多樣性維持與調(diào)控。2021年在原網(wǎng)圍區(qū)開(kāi)展1 km2浮葉植物刈割試驗(yàn)示范,試驗(yàn)區(qū)監(jiān)測(cè)表明,沉水植物群落結(jié)構(gòu)得到改善,生物多樣性提高,生態(tài)完整性得到恢復(fù),水質(zhì)清新,顯現(xiàn)出較好的效果,為東太湖水生植被管理探索了經(jīng)驗(yàn)。

    3.3.2 加強(qiáng)對(duì)東太湖生態(tài)環(huán)境監(jiān)測(cè)并建立湖泊生態(tài)檔案 東太湖網(wǎng)圍拆除后水環(huán)境和水生態(tài)波動(dòng)較大,建立生態(tài)檔案掌握發(fā)展動(dòng)態(tài)和變化規(guī)律對(duì)于東太湖生態(tài)修復(fù)和綜合管控至關(guān)重要。因此,需要加強(qiáng)對(duì)東太湖水質(zhì)、浮游藻類(lèi)、水生植被、魚(yú)類(lèi)等生物資源,以及水位、流速和水下光照等物理環(huán)境因子的長(zhǎng)期動(dòng)態(tài)監(jiān)測(cè),合理調(diào)控東太湖的水生植物群落結(jié)構(gòu)及生物量和覆蓋度,科學(xué)評(píng)估東太湖保持Ⅲ類(lèi)以上水質(zhì)需要的物理和生物條件,確保東太湖良好生態(tài)系統(tǒng)重構(gòu)與水生態(tài)安全。

    3.3.3 創(chuàng)建東太湖漁業(yè)資源養(yǎng)護(hù)與利用新模式 東太湖草型生態(tài)系統(tǒng)是生物多樣性和漁業(yè)種質(zhì)資源保護(hù)的有利條件,為生態(tài)系統(tǒng)服務(wù)功能提升和生態(tài)漁業(yè)發(fā)展奠定了基礎(chǔ)。東太湖不能發(fā)展大規(guī)模集約化的人工養(yǎng)殖,但東太湖草型結(jié)構(gòu)及優(yōu)質(zhì)的餌料生物資源,需要在保護(hù)中利用,讓生態(tài)系統(tǒng)內(nèi)部形成良性循環(huán)。東太湖宜開(kāi)展河蟹、團(tuán)頭魴、鱖、河川沙塘鱧等太湖魚(yú)米之鄉(xiāng)特色魚(yú)類(lèi)品種的增殖放流,既可修復(fù)漁業(yè)種質(zhì)資源、維護(hù)生態(tài)系統(tǒng)完整性,又實(shí)現(xiàn)了沉水植被-底棲動(dòng)物-浮游生物-漁業(yè)的食物鏈能流耦合。通過(guò)構(gòu)建“道法自然”、“人放天養(yǎng)”的生態(tài)模式,在提升東太湖生態(tài)系統(tǒng)服務(wù)功能的同時(shí),為東太湖網(wǎng)圍拆除后獨(dú)具特色的漁業(yè)高質(zhì)量發(fā)展創(chuàng)造了新的模式,也是綠水青山就是為金山銀山生態(tài)產(chǎn)品價(jià)值轉(zhuǎn)化的最好詮釋。

    4 結(jié)論

    1)30年來(lái)東太湖的TN、TP、CODMn和Chl.a濃度整體呈上升趨勢(shì)。網(wǎng)圍養(yǎng)蟹的水質(zhì)相對(duì)于網(wǎng)圍養(yǎng)魚(yú)較優(yōu),但均隨著養(yǎng)殖強(qiáng)度的增大而變差,網(wǎng)圍規(guī)??s減對(duì)改善水質(zhì)有正面效應(yīng)。2008年后太浦河口和東茭咀透明度下降明顯,現(xiàn)在維持在30 cm左右;營(yíng)養(yǎng)鹽濃度和Chl.a濃度升高顯著。東茭咀和太浦河口水域的Chl.a濃度與太浦河口出湖水量呈顯著正相關(guān)。網(wǎng)圍拆除、出湖水量增加和水生植被消失造成的水動(dòng)力作用增強(qiáng),加上原網(wǎng)圍區(qū)底泥營(yíng)養(yǎng)鹽釋放,是引起東太湖南部區(qū)域營(yíng)養(yǎng)鹽濃度大幅上升的主要原因。

    2)東太湖網(wǎng)圍養(yǎng)殖拆除,目標(biāo)是重構(gòu)良好的生態(tài)系統(tǒng),但3年來(lái)的環(huán)境演變并沒(méi)有顯現(xiàn)出最初設(shè)想的效果,需要加強(qiáng)對(duì)東太湖的管理,及早開(kāi)展水生植物資源調(diào)控,強(qiáng)化水質(zhì)凈化,加強(qiáng)對(duì)東太湖泊湖生態(tài)環(huán)境監(jiān)測(cè)并建立湖泊生態(tài)檔案,創(chuàng)建漁業(yè)新模式,確?;謴?fù)健康自然的草型生態(tài)系統(tǒng)格局。

    致謝:感謝楊元龍副研究員提供的1990-2005年?yáng)|太湖的水質(zhì)監(jiān)測(cè)數(shù)據(jù),感謝羅菊花副研究員提供的2017-2021年?yáng)|太湖水生植被類(lèi)群面積數(shù)據(jù),感謝闞可聰、仝天衡在樣品采集上的幫助。

    猜你喜歡
    沉水植物太湖水體
    農(nóng)村黑臭水體治理和污水處理淺探
    生態(tài)修復(fù)理念在河道水體治理中的應(yīng)用
    沉水植物生長(zhǎng)影響因子研究進(jìn)展
    花卉(2020年12期)2020-01-09 07:12:52
    沉水植物種植工程實(shí)施后物種多樣性的變化
    ——以武漢紫陽(yáng)湖為例
    綠色科技(2019年2期)2019-05-21 07:43:18
    太湖思變2017
    玩具世界(2017年4期)2017-07-21 13:27:24
    廣元:治理黑臭水體 再還水清岸美
    太湖攬春
    寶藏(2017年2期)2017-03-20 13:16:42
    太湖
    中亞信息(2016年3期)2016-12-01 06:08:24
    沉水植物在修復(fù)富營(yíng)養(yǎng)湖體的研究進(jìn)展
    大科技(2016年1期)2016-07-17 02:42:07
    太湖一角
    久久久久国内视频| 身体一侧抽搐| 欧美乱色亚洲激情| 在线观看舔阴道视频| 亚洲专区字幕在线| 不卡av一区二区三区| 欧美乱色亚洲激情| 啦啦啦韩国在线观看视频| 日本三级黄在线观看| 亚洲av电影不卡..在线观看| 亚洲男人的天堂狠狠| 免费不卡黄色视频| 一区在线观看完整版| 十八禁人妻一区二区| 一本久久中文字幕| 国产亚洲av嫩草精品影院| 久久久久九九精品影院| 1024视频免费在线观看| 日本a在线网址| 国内毛片毛片毛片毛片毛片| 露出奶头的视频| 国产一区二区三区在线臀色熟女| 中文字幕人妻丝袜一区二区| 少妇粗大呻吟视频| 变态另类成人亚洲欧美熟女 | 在线国产一区二区在线| 国产亚洲精品av在线| 一本大道久久a久久精品| 国产精品乱码一区二三区的特点 | 久久久久久久午夜电影| 国产欧美日韩精品亚洲av| 深夜精品福利| 国产黄a三级三级三级人| 亚洲人成77777在线视频| 俄罗斯特黄特色一大片| 黄频高清免费视频| 两个人视频免费观看高清| 久久精品91无色码中文字幕| 一边摸一边抽搐一进一小说| 两性夫妻黄色片| 男人操女人黄网站| 禁无遮挡网站| 午夜精品久久久久久毛片777| 国产精品秋霞免费鲁丝片| 久久精品国产99精品国产亚洲性色 | 国产一区二区三区综合在线观看| 国产一区二区在线av高清观看| 国产av精品麻豆| 满18在线观看网站| 1024视频免费在线观看| 一区在线观看完整版| 女人精品久久久久毛片| 两人在一起打扑克的视频| 首页视频小说图片口味搜索| 成人国语在线视频| 十分钟在线观看高清视频www| 欧美中文综合在线视频| 成人18禁高潮啪啪吃奶动态图| 黄色a级毛片大全视频| 丁香六月欧美| av天堂久久9| 亚洲精品av麻豆狂野| 色播在线永久视频| 国产熟女xx| 国产精品乱码一区二三区的特点 | 男女之事视频高清在线观看| 中国美女看黄片| 精品欧美国产一区二区三| 99国产精品一区二区蜜桃av| 中文字幕另类日韩欧美亚洲嫩草| 国产高清激情床上av| 国产一卡二卡三卡精品| 女性生殖器流出的白浆| 欧美日韩亚洲综合一区二区三区_| 国产精品久久久久久人妻精品电影| 这个男人来自地球电影免费观看| 国产极品粉嫩免费观看在线| 亚洲免费av在线视频| 一卡2卡三卡四卡精品乱码亚洲| 九色国产91popny在线| 成人欧美大片| 在线十欧美十亚洲十日本专区| 精品不卡国产一区二区三区| 国产欧美日韩一区二区三区在线| av欧美777| 亚洲精品国产精品久久久不卡| 俄罗斯特黄特色一大片| www.www免费av| tocl精华| 成年女人毛片免费观看观看9| 操美女的视频在线观看| 性欧美人与动物交配| 国产精品永久免费网站| 一级毛片精品| 免费在线观看视频国产中文字幕亚洲| 熟妇人妻久久中文字幕3abv| 侵犯人妻中文字幕一二三四区| 亚洲精品中文字幕在线视频| 男女午夜视频在线观看| avwww免费| 黄色视频不卡| 亚洲av第一区精品v没综合| 高清毛片免费观看视频网站| 老司机午夜十八禁免费视频| 免费久久久久久久精品成人欧美视频| 两个人看的免费小视频| 真人做人爱边吃奶动态| 黄色毛片三级朝国网站| 中文字幕精品免费在线观看视频| 少妇 在线观看| 日韩免费av在线播放| 日韩 欧美 亚洲 中文字幕| 亚洲激情在线av| 免费高清在线观看日韩| 人成视频在线观看免费观看| tocl精华| 久久精品国产99精品国产亚洲性色 | aaaaa片日本免费| 国产高清videossex| 麻豆一二三区av精品| 手机成人av网站| 12—13女人毛片做爰片一| 久久国产乱子伦精品免费另类| 国产精品香港三级国产av潘金莲| 亚洲,欧美精品.| 黑丝袜美女国产一区| 成年人黄色毛片网站| 伦理电影免费视频| 成人国产综合亚洲| 亚洲国产欧美日韩在线播放| 亚洲国产毛片av蜜桃av| 国产精品,欧美在线| 亚洲精品国产精品久久久不卡| а√天堂www在线а√下载| 国产色视频综合| 中出人妻视频一区二区| 正在播放国产对白刺激| 女人爽到高潮嗷嗷叫在线视频| 免费无遮挡裸体视频| 午夜精品国产一区二区电影| 波多野结衣av一区二区av| 在线免费观看的www视频| 国产欧美日韩精品亚洲av| 一夜夜www| 久久性视频一级片| 亚洲免费av在线视频| 国内精品久久久久精免费| 曰老女人黄片| 午夜精品国产一区二区电影| 午夜福利一区二区在线看| 国产黄a三级三级三级人| 搡老妇女老女人老熟妇| 国产色视频综合| 欧美另类亚洲清纯唯美| 91九色精品人成在线观看| 亚洲 欧美 日韩 在线 免费| 美女午夜性视频免费| 欧美激情 高清一区二区三区| av在线天堂中文字幕| 久久久久国内视频| 可以免费在线观看a视频的电影网站| 久久国产精品人妻蜜桃| 岛国在线观看网站| 女人被躁到高潮嗷嗷叫费观| 国产又色又爽无遮挡免费看| АⅤ资源中文在线天堂| 国产成人精品无人区| 丝袜美腿诱惑在线| 国内毛片毛片毛片毛片毛片| 亚洲中文av在线| 日韩欧美免费精品| 欧美色视频一区免费| 大码成人一级视频| 脱女人内裤的视频| 女同久久另类99精品国产91| 久久久国产欧美日韩av| 国产亚洲欧美在线一区二区| 在线观看午夜福利视频| 女人被狂操c到高潮| 一级毛片精品| 亚洲熟女毛片儿| 日本 av在线| 极品教师在线免费播放| 免费在线观看亚洲国产| 一级毛片精品| 18禁观看日本| 亚洲第一青青草原| 男女之事视频高清在线观看| 精品一区二区三区四区五区乱码| 久久人人爽av亚洲精品天堂| 国产伦一二天堂av在线观看| 亚洲成人久久性| 男人舔女人下体高潮全视频| 久久久精品欧美日韩精品| 成人手机av| 天堂影院成人在线观看| 国产成+人综合+亚洲专区| 国产成人欧美| 午夜福利一区二区在线看| 搡老妇女老女人老熟妇| 热99re8久久精品国产| 91老司机精品| 一二三四在线观看免费中文在| 久久香蕉国产精品| 69精品国产乱码久久久| 国产精品二区激情视频| 久久人妻熟女aⅴ| 日韩一卡2卡3卡4卡2021年| 亚洲欧美日韩高清在线视频| 手机成人av网站| 搡老岳熟女国产| 免费高清在线观看日韩| 国产熟女午夜一区二区三区| 久久伊人香网站| 女性被躁到高潮视频| 国产野战对白在线观看| 亚洲av美国av| 麻豆av在线久日| 亚洲精品国产一区二区精华液| 香蕉久久夜色| 老司机靠b影院| 十八禁网站免费在线| 69av精品久久久久久| 成人国语在线视频| 国产亚洲精品av在线| 999久久久国产精品视频| 亚洲,欧美精品.| 一级片免费观看大全| 人妻久久中文字幕网| 久久久久久久久中文| 亚洲精品中文字幕在线视频| 男女床上黄色一级片免费看| 国产成人影院久久av| 亚洲第一电影网av| 国产精品国产高清国产av| 色在线成人网| 免费女性裸体啪啪无遮挡网站| 在线av久久热| 9191精品国产免费久久| 亚洲伊人色综图| 成在线人永久免费视频| 在线观看免费午夜福利视频| 可以在线观看毛片的网站| 亚洲精品国产区一区二| 91麻豆精品激情在线观看国产| 免费看美女性在线毛片视频| 三级毛片av免费| 男男h啪啪无遮挡| ponron亚洲| 久久香蕉国产精品| 亚洲一卡2卡3卡4卡5卡精品中文| 国产99久久九九免费精品| 纯流量卡能插随身wifi吗| 国产三级黄色录像| 99精品久久久久人妻精品| 脱女人内裤的视频| 亚洲 欧美 日韩 在线 免费| 一进一出好大好爽视频| 一本久久中文字幕| 亚洲avbb在线观看| 成年版毛片免费区| 男女下面进入的视频免费午夜 | 丝袜在线中文字幕| 亚洲av美国av| 黄片大片在线免费观看| 国产aⅴ精品一区二区三区波| 亚洲五月色婷婷综合| 日韩欧美一区视频在线观看| 免费搜索国产男女视频| 亚洲精品国产一区二区精华液| 欧美精品亚洲一区二区| 久久精品aⅴ一区二区三区四区| 人成视频在线观看免费观看| 国产精品永久免费网站| 少妇 在线观看| 亚洲视频免费观看视频| 一本久久中文字幕| 级片在线观看| 老汉色∧v一级毛片| 亚洲第一青青草原| 91国产中文字幕| 50天的宝宝边吃奶边哭怎么回事| 国产一级毛片七仙女欲春2 | 久久久精品国产亚洲av高清涩受| 人妻丰满熟妇av一区二区三区| 夜夜看夜夜爽夜夜摸| www.999成人在线观看| 后天国语完整版免费观看| 九色国产91popny在线| 亚洲一区二区三区色噜噜| 中文字幕人妻丝袜一区二区| 精品熟女少妇八av免费久了| 国产精品野战在线观看| 波多野结衣高清无吗| 国产精品精品国产色婷婷| 首页视频小说图片口味搜索| 欧美性长视频在线观看| 久热这里只有精品99| 久久天躁狠狠躁夜夜2o2o| 久久影院123| 午夜视频精品福利| 精品国产乱码久久久久久男人| 校园春色视频在线观看| 黄色成人免费大全| 99在线人妻在线中文字幕| 亚洲第一欧美日韩一区二区三区| 麻豆成人av在线观看| 亚洲av电影在线进入| 一卡2卡三卡四卡精品乱码亚洲| 激情视频va一区二区三区| 国产一区在线观看成人免费| 91av网站免费观看| 成熟少妇高潮喷水视频| 成人国产综合亚洲| 在线永久观看黄色视频| 亚洲黑人精品在线| 国产高清videossex| 老司机靠b影院| 精品久久久久久久毛片微露脸| aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品精品国产色婷婷| 两人在一起打扑克的视频| 国产精品乱码一区二三区的特点 | 色av中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 久久婷婷人人爽人人干人人爱 | 久久香蕉精品热| 亚洲精品粉嫩美女一区| 国产成人欧美在线观看| 一边摸一边抽搐一进一出视频| 91麻豆精品激情在线观看国产| 国产极品粉嫩免费观看在线| 9191精品国产免费久久| 亚洲美女黄片视频| 老鸭窝网址在线观看| 亚洲色图综合在线观看| av天堂在线播放| 老鸭窝网址在线观看| 久久久久久久久免费视频了| 12—13女人毛片做爰片一| 国产精品精品国产色婷婷| 一卡2卡三卡四卡精品乱码亚洲| 午夜两性在线视频| 国产一区二区激情短视频| 十分钟在线观看高清视频www| 亚洲午夜精品一区,二区,三区| 男女床上黄色一级片免费看| 日韩一卡2卡3卡4卡2021年| 美女国产高潮福利片在线看| 天堂动漫精品| 亚洲欧美激情综合另类| 欧美日韩瑟瑟在线播放| 岛国视频午夜一区免费看| 免费在线观看亚洲国产| 手机成人av网站| 一二三四社区在线视频社区8| 99久久精品国产亚洲精品| 两性午夜刺激爽爽歪歪视频在线观看 | 精品欧美国产一区二区三| 欧美乱码精品一区二区三区| av在线播放免费不卡| 精品欧美国产一区二区三| 欧美 亚洲 国产 日韩一| 国产av一区二区精品久久| 波多野结衣一区麻豆| 国产免费男女视频| 午夜免费观看网址| 三级毛片av免费| 757午夜福利合集在线观看| 亚洲av美国av| 热99re8久久精品国产| 亚洲午夜精品一区,二区,三区| 老汉色av国产亚洲站长工具| 欧美丝袜亚洲另类 | 亚洲第一欧美日韩一区二区三区| 一级黄色大片毛片| 亚洲中文av在线| 国产亚洲欧美精品永久| 日韩中文字幕欧美一区二区| 嫩草影院精品99| 大码成人一级视频| 黄色 视频免费看| 18禁美女被吸乳视频| 正在播放国产对白刺激| 亚洲专区中文字幕在线| 色综合站精品国产| 麻豆国产av国片精品| 最新在线观看一区二区三区| 国产成人欧美在线观看| 精品久久蜜臀av无| 日本一区二区免费在线视频| 90打野战视频偷拍视频| av有码第一页| 成年人黄色毛片网站| 丝袜美足系列| av欧美777| bbb黄色大片| 国产精品1区2区在线观看.| 露出奶头的视频| 午夜成年电影在线免费观看| 色哟哟哟哟哟哟| 久久人妻熟女aⅴ| 一级毛片女人18水好多| 精品国产一区二区三区四区第35| 青草久久国产| 精品久久久精品久久久| 亚洲片人在线观看| 女警被强在线播放| 黄色成人免费大全| 国产亚洲欧美在线一区二区| 成人亚洲精品av一区二区| 午夜福利高清视频| 两个人视频免费观看高清| 夜夜躁狠狠躁天天躁| 日本 欧美在线| 国产精品野战在线观看| 欧美黑人精品巨大| 国产国语露脸激情在线看| 成人国产综合亚洲| 琪琪午夜伦伦电影理论片6080| www.精华液| www.自偷自拍.com| 久久精品91无色码中文字幕| 国产成+人综合+亚洲专区| 男女床上黄色一级片免费看| 狂野欧美激情性xxxx| 999久久久国产精品视频| a级毛片在线看网站| 色播在线永久视频| 岛国视频午夜一区免费看| 99久久综合精品五月天人人| 国产成人欧美| 69av精品久久久久久| 日日夜夜操网爽| 亚洲精品久久成人aⅴ小说| 18禁黄网站禁片午夜丰满| 欧美日韩精品网址| 嫩草影院精品99| 最近最新中文字幕大全电影3 | 看片在线看免费视频| 国产蜜桃级精品一区二区三区| 欧美丝袜亚洲另类 | 国产激情欧美一区二区| avwww免费| 桃色一区二区三区在线观看| 国产99白浆流出| 母亲3免费完整高清在线观看| 国产精品亚洲美女久久久| 国产黄a三级三级三级人| www.熟女人妻精品国产| 免费看a级黄色片| 天堂影院成人在线观看| 欧美 亚洲 国产 日韩一| 19禁男女啪啪无遮挡网站| 国产av一区在线观看免费| 最好的美女福利视频网| 欧美色视频一区免费| 久久国产精品影院| 男女之事视频高清在线观看| 波多野结衣巨乳人妻| 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 人人妻,人人澡人人爽秒播| 在线观看免费日韩欧美大片| 夜夜躁狠狠躁天天躁| www.自偷自拍.com| 一本久久中文字幕| 又黄又爽又免费观看的视频| 久久久精品欧美日韩精品| 夜夜躁狠狠躁天天躁| 一进一出好大好爽视频| 国产亚洲欧美精品永久| 日本五十路高清| 国产精品美女特级片免费视频播放器 | 亚洲在线自拍视频| 国产不卡一卡二| 女人被狂操c到高潮| 韩国精品一区二区三区| 99在线视频只有这里精品首页| 国产精品亚洲一级av第二区| av网站免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲第一欧美日韩一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 在线观看免费日韩欧美大片| 黄色成人免费大全| 亚洲aⅴ乱码一区二区在线播放 | 色精品久久人妻99蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 高清在线国产一区| 亚洲国产精品999在线| 极品人妻少妇av视频| 天天躁夜夜躁狠狠躁躁| 后天国语完整版免费观看| 久久香蕉精品热| 久久精品成人免费网站| 亚洲午夜精品一区,二区,三区| 黄频高清免费视频| 国产伦一二天堂av在线观看| 国产精品爽爽va在线观看网站 | 精品国产一区二区久久| 午夜视频精品福利| 99香蕉大伊视频| 身体一侧抽搐| 给我免费播放毛片高清在线观看| 一级a爱视频在线免费观看| 欧美日韩亚洲综合一区二区三区_| 色综合婷婷激情| 欧美激情极品国产一区二区三区| 国产精品一区二区免费欧美| 久久中文看片网| 黄色丝袜av网址大全| 久久婷婷人人爽人人干人人爱 | 巨乳人妻的诱惑在线观看| 国产单亲对白刺激| 免费看a级黄色片| 亚洲人成伊人成综合网2020| 日韩欧美国产一区二区入口| 国产麻豆成人av免费视频| 久久久久久人人人人人| 两性夫妻黄色片| 999精品在线视频| 一区二区三区精品91| 精品一区二区三区视频在线观看免费| 国产精品久久视频播放| 女性被躁到高潮视频| 欧美成人性av电影在线观看| 无限看片的www在线观看| 日韩视频一区二区在线观看| 久久久久久久久中文| 亚洲精品国产精品久久久不卡| 久久久久久久久免费视频了| 1024香蕉在线观看| 午夜久久久久精精品| 国产不卡一卡二| www.精华液| 在线观看免费视频网站a站| 熟女少妇亚洲综合色aaa.| 久久狼人影院| 成熟少妇高潮喷水视频| 俄罗斯特黄特色一大片| 黄色毛片三级朝国网站| 天堂√8在线中文| 看免费av毛片| 88av欧美| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 正在播放国产对白刺激| 他把我摸到了高潮在线观看| 亚洲精品国产色婷婷电影| 国产免费av片在线观看野外av| 亚洲国产中文字幕在线视频| 日韩 欧美 亚洲 中文字幕| 狂野欧美激情性xxxx| 国产精品久久久av美女十八| 狂野欧美激情性xxxx| 久久国产乱子伦精品免费另类| x7x7x7水蜜桃| 神马国产精品三级电影在线观看 | 91成人精品电影| 日韩一卡2卡3卡4卡2021年| www.自偷自拍.com| 亚洲人成伊人成综合网2020| 精品人妻1区二区| 三级毛片av免费| 韩国精品一区二区三区| 性少妇av在线| 午夜免费成人在线视频| 性少妇av在线| 免费看a级黄色片| 国产区一区二久久| 国产男靠女视频免费网站| 老司机午夜十八禁免费视频| 日韩精品青青久久久久久| x7x7x7水蜜桃| 村上凉子中文字幕在线| 99久久综合精品五月天人人| 久久亚洲精品不卡| 成人18禁在线播放| 国产91精品成人一区二区三区| 成人18禁高潮啪啪吃奶动态图| 在线免费观看的www视频| 日韩大尺度精品在线看网址 | 久久人妻av系列| 亚洲熟妇熟女久久| 一级毛片精品| 琪琪午夜伦伦电影理论片6080| 国产成人精品久久二区二区免费| 欧美日韩乱码在线| 久久久久国产精品人妻aⅴ院| 老鸭窝网址在线观看| 亚洲三区欧美一区| 麻豆国产av国片精品| 黑人巨大精品欧美一区二区蜜桃| 夜夜夜夜夜久久久久| 国产又色又爽无遮挡免费看| 久久久久精品国产欧美久久久| 亚洲,欧美精品.| 亚洲 欧美一区二区三区| 久久精品人人爽人人爽视色| 亚洲av成人一区二区三| 精品国产乱码久久久久久男人| 超碰成人久久| 欧美激情高清一区二区三区| 精品国内亚洲2022精品成人| 国产精品,欧美在线| 看免费av毛片| 午夜免费观看网址| 色av中文字幕| 在线观看一区二区三区| 国产精品永久免费网站| 亚洲熟妇中文字幕五十中出| 黄色毛片三级朝国网站| 久久中文字幕一级|