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

    浙江省銅山源水庫水體偏堿污染成因研究

    2023-08-26 00:09:07錢旭蔡梅王元元白君瑞
    人民長江 2023年8期
    關(guān)鍵詞:水質(zhì)污染

    錢旭 蔡梅 王元元 白君瑞

    摘要:

    針對中國東南丘陵地區(qū)水庫出現(xiàn)不同程度的水質(zhì)偏堿,部分水庫在局部時段內(nèi)水體pH值指標(biāo)頻繁超出地表水質(zhì)量標(biāo)準(zhǔn)的問題,以浙江衢州大(2)型水庫銅山源水庫為例,實施水質(zhì)監(jiān)測和水生態(tài)調(diào)查。通過水文、氣象、水質(zhì)、水生態(tài)指標(biāo)之間的相關(guān)性分析,揭示了自然環(huán)境本底水質(zhì)、關(guān)鍵氣象水文過程、外源污染以及藻類生長對水庫水質(zhì)偏堿污染的影響機制。結(jié)果表明:銅山源水庫水體偏堿污染具有長期性的特征,且在每年的4~8月偏堿現(xiàn)象最為嚴(yán)重。水庫周邊自然環(huán)境下的本底水質(zhì)偏堿,銅山源水庫入庫河流水體pH值高達8.02~8.47,局部時段高于地表水標(biāo)準(zhǔn)值。庫區(qū)水體pH值指標(biāo)與水溫、DO、營養(yǎng)性物質(zhì)(TP、TN)、葉綠素a等指標(biāo)在一定程度上隨著時間變化呈現(xiàn)出同步升降的正相關(guān)性。pH值變化受多種因素的綜合影響,主要表現(xiàn)為直接或間接地促進藻類光合作用增強,導(dǎo)致光合作用效應(yīng)大于水體水生動植物的呼吸作用,從而導(dǎo)致水體偏堿污染。

    關(guān) 鍵 詞:

    水庫污染; pH值; 藻類; 溶解氧; 葉綠素a; 銅山源水庫

    中圖法分類號: X131.2

    文獻標(biāo)志碼: A

    DOI:10.16232/j.cnki.1001-4179.2023.08.010

    0 引 言

    近年來,水庫逐漸成為城鎮(zhèn)生產(chǎn)、生活用水的重要來源之一,尤其是中國東南沿海省市,水庫型水源地日益成為經(jīng)濟社會可持續(xù)發(fā)展的重要保障[1]。然而,中國地表湖庫近年來多發(fā)生水體偏堿污染現(xiàn)象,北方地區(qū)天津市的于橋水庫[2-3]、山東省的大沽夾河流域[4],東部地區(qū)上海市的東風(fēng)西沙水庫[5],東南地區(qū)福建省的東牙溪水庫[6],南部地區(qū)的廣西大王灘水庫[7]、云南漁洞水庫[8]等,均出現(xiàn)過不同程度的水體pH值偏高現(xiàn)象。水庫作為一種天然或人工的湖泊,水質(zhì)特點與天然地表湖泊相似。有統(tǒng)計表明:中國湖泊的pH值具有地帶性的分布特點[9],蒙新、青藏地區(qū)絕大多數(shù)呈現(xiàn)堿性或強堿性;云貴和黃淮海地區(qū)呈弱堿性;而東北及長江中下游地區(qū)地表湖泊pH值較低,一般在6.5~8.3之間,呈現(xiàn)中性或微堿性。在21世紀(jì)前10 a,浙江省部分水庫,如千島湖(新安江)水庫、老虎潭水庫、長潭水庫等,水體長期呈現(xiàn)中性偏堿性(pH值為6.70~9.17)[10]。2016年,浙江省水庫水質(zhì)安全指標(biāo)評價以Ⅰ~Ⅲ類為主,藻類大量生長的4~9月,部分水庫存在pH值超出地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)的現(xiàn)象[11]。位于浙江省衢州市境內(nèi)的銅山源水庫作為衢江流域的大(2)型水庫,承載了浙西地區(qū)的數(shù)千萬公頃農(nóng)田的灌溉和周邊地區(qū)水力供電及漁業(yè)和觀賞旅游功能,發(fā)揮著極其重要的作用[12]。銅山源水庫國控斷面監(jiān)測數(shù)據(jù)顯示,2016~2020年銅山源水庫水體pH值在春末及夏季(4~8月)較地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)規(guī)定的pH值標(biāo)準(zhǔn)高值偏高(pH值為9)。

    湖庫水體pH值受多種因素綜合影響,是水體中二氧化碳含量、有機酸和水體污染的表征指標(biāo)[13],同時也是反映水體化學(xué)的一個重要指標(biāo),常被作為水體生態(tài)環(huán)境影響因素來考慮[14]。因此,此類酸堿平衡失調(diào)現(xiàn)象是水庫水生生態(tài)系統(tǒng)失衡的具體表現(xiàn)形式之一。目前,針對地表水體pH值超標(biāo)的研究較少,國內(nèi)一些研究認為,湖庫水體的pH值主要與水體中積累的營養(yǎng)鹽和藻類的繁殖有關(guān)[8],當(dāng)水體中藻類數(shù)量上升到一定數(shù)量級時,其生物規(guī)模、生命活動的旺盛程度等必然對水體酸堿性變化起到主導(dǎo)作用[15-16]。黃歲樑等[17-18]研究指出在水體交換強的夏季和秋季,養(yǎng)殖和非養(yǎng)殖水體中pH值與溶解氧、葉綠素a濃度存在一定程度的正相關(guān)關(guān)系,王志紅等[15]指出晝夜pH值差值與藻類的生長態(tài)勢和生命活性之間存在一定相關(guān)性。

    目前,國內(nèi)外學(xué)者多聚焦于對浮游植物本身生長習(xí)性及其與營養(yǎng)物質(zhì)等因素之間的關(guān)系研究,而浮游植物與水體pH值之間的關(guān)系尚未引起普遍關(guān)注,也未有權(quán)威的共識。針對中國東南沿海地區(qū)水庫水體偏堿污染現(xiàn)象研究也尚未全面開展,湖庫水體偏堿問題成因不明,現(xiàn)階段已嚴(yán)重阻礙國控、省控斷面水質(zhì)穩(wěn)定達標(biāo),亟需開展水體偏堿問題成因探究。此外,目前已有的研究多是基于室內(nèi)研究和原位試驗,但作為實際地表水體的水庫是一個較為復(fù)雜的綜合水生態(tài)系統(tǒng),水體酸堿性可能受到包括水文氣象因素、外源來水以及沉積物中的內(nèi)源污染等在內(nèi)的多種因素影響,更需要結(jié)合實際水質(zhì)和水生態(tài)監(jiān)測調(diào)查成果開展綜合性研究。因此,本研究針對東南沿海地區(qū)水庫普遍性偏堿問題,以浙江銅山源水庫為研究對象,采用現(xiàn)場水質(zhì)監(jiān)測和水生態(tài)調(diào)查,科學(xué)評價水庫水體的酸堿狀況,探究水庫水體偏堿污染成因,為水庫的水生態(tài)安全保障提供支撐。

    1 材料與方法

    1.1 研究對象概況

    銅山源水庫(118°15′~119°E,29°~29°15′N)位于浙江省衢州市衢江流域(圖1),集水面積180 km2,總庫容1.715億m3,年平均入庫水量2.103億m3,換水周期約為148~220 d,是一座以農(nóng)業(yè)灌溉為主,結(jié)合防洪、發(fā)電、潔水漁業(yè)、水資源保護等功能的大(2)型水庫,水質(zhì)保護目標(biāo)為Ⅱ~Ⅲ類。水庫所在區(qū)域?qū)賮啛釒Ъ撅L(fēng)性氣候區(qū),受東南季風(fēng)影響,有溫暖濕潤、四季分明、冬夏長春秋短、光照充足、降水豐沛而季節(jié)分配不均的地帶性特征。庫區(qū)雨量充沛,多年平均降雨量1 500~2 300 mm。降雨年內(nèi)分配不均,主要集中在3~7月,其降水量占全年降水量的60%以上。

    1.2 樣品采集

    1.2.1 監(jiān)測布點

    根據(jù)HJ/T 91-2002《地表水和污水監(jiān)測技術(shù)規(guī)范》中對于湖泊、水庫水質(zhì)監(jiān)測布點原則和方法的相關(guān)要求,結(jié)合水庫水深情況,在水庫進水區(qū)、出水區(qū)、深水區(qū)、淺水區(qū)、湖心區(qū)以及岸邊區(qū)相關(guān)位置布設(shè)6處(St1~St6)水質(zhì)和水生態(tài)監(jiān)測點位。其中,St6監(jiān)測點為全庫水深最深處,水深達30~40 m,且處于水庫的國控斷面[19]附近。同時,在3條入庫河道(從西到東分別為銅山溪、廟前溪、白鶴山入湖河道)入庫口上游設(shè)置St7~St8兩處水質(zhì)監(jiān)測點。監(jiān)測布點位置見圖2。

    1.2.2 樣品采集

    2021年4~8月,對St1~St8八處監(jiān)測點實施水質(zhì)加密監(jiān)測,對St1~St6六處庫區(qū)點位開展水生態(tài)調(diào)查。其中,水質(zhì)監(jiān)測頻次為1旬1次;水生態(tài)調(diào)查在3、5、7月各開展1次。

    此外,為進一步明晰庫區(qū)水質(zhì)pH值和溶解氧指標(biāo)的相關(guān)性,在4~8月每月各選取1 d開展全天(2 h間隔)的pH值和溶解氧樣本采集。為了分析庫區(qū)水體垂向pH值變化,選取深水區(qū)St6監(jiān)測點,每月開展1次垂向分層采樣。

    1.3 樣品處理與測定

    (1) 水質(zhì)指標(biāo)測定。水體pH值、水溫、溶解氧、電導(dǎo)率在現(xiàn)場采用YSI Pro plus(美國)儀器測定總氮、氨氮、總磷、高錳酸鹽指數(shù)等其他水質(zhì)指標(biāo)取水樣冷藏運輸至實驗室,分別采用堿性過硫酸鉀紫外分光光度法、奈氏試劑分光光度法、鉬酸銨分光光度法和高錳酸鹽法測定[20];葉綠素a在采樣后2 h內(nèi),抽濾處理并用液氮(-180 ℃)保存后采用丙酮-分光光度法測定。

    (2) 水生態(tài)指標(biāo)測定。水生態(tài)調(diào)查內(nèi)容主要包括浮游植物、浮游動物、底棲動物、魚類、水生植物及底質(zhì)等。浮游植物、浮游動物、底棲動物的調(diào)查指標(biāo)主要是密度、生物量、種類組成、優(yōu)勢種;魚類調(diào)查指標(biāo)主要是魚類種類組成、優(yōu)勢種;水生植物主要調(diào)查水生植物種類組成、覆蓋度、優(yōu)勢種。

    1.4 數(shù)據(jù)處理與分析

    為開展水庫水體pH值與其他水質(zhì)指標(biāo)之間的關(guān)聯(lián)性分析,引入灰色關(guān)聯(lián)度分析法(Grey Relation Analysis,GRA)。對于具備時間變化屬性的因子,在樣本數(shù)量不大的情況下,將數(shù)據(jù)進行歸一化處理、最值矯正后,通過計算數(shù)據(jù)之間的關(guān)聯(lián)度,度量兩組隨時間變化數(shù)據(jù)的變化走勢是否相近,以此判斷兩個因子之間的關(guān)聯(lián)性。灰色關(guān)聯(lián)度的取值范圍為[0,1],值越大的說明灰色關(guān)聯(lián)度越大[21]。

    本研究采用WPS office中WPS表格進行數(shù)據(jù)的處理和灰色關(guān)聯(lián)度統(tǒng)計分析;采用Origin 2017進行分析制圖。

    2 結(jié)果與討論

    2.1 水庫水體偏堿污染特征

    由2016~2020年銅山源水庫水體國控斷面每月pH值變化情況(圖3)可見,銅山源水庫水體偏堿現(xiàn)象在近5 a時間內(nèi)中長期存在。主要呈現(xiàn)出秋冬季節(jié)水體pH值較低而春夏季節(jié)較高的特征,尤其每年的4~8月為水庫水體偏堿污染較為嚴(yán)重的時期。有調(diào)查指出,浙江省內(nèi)其他水庫,如寧波市東錢湖、臺州市里石門水庫和長潭水庫等,也普遍存在著水體pH值超標(biāo)現(xiàn)象[11]。

    本研究選取銅山源水庫水體偏堿污染現(xiàn)象較為嚴(yán)重的4~8月,按旬開展水質(zhì)加密監(jiān)測,測得St1~St6各監(jiān)測點水體pH值變化情況如圖4所示。4~8月銅山源水庫庫區(qū)內(nèi)St1~St6監(jiān)測點的pH值均呈現(xiàn)偏堿性,數(shù)值在8.27~9.75之間。時間尺度上,St1~St6各監(jiān)測點在4月pH值較低(8.27~8.94),此后水體pH值呈現(xiàn)持續(xù)升高的態(tài)勢,5月上旬水庫水體pH值已高于地表水質(zhì)量標(biāo)準(zhǔn)值(9),隨后進入持續(xù)性高位波動狀態(tài),特別是6月全月和8月中下旬水體pH值普遍較高,在9.5上下波動,偏堿污染較為嚴(yán)重??臻g尺度上,St1~St6各監(jiān)測點的水體pH值變化趨勢基本一致,且同一時間各監(jiān)測點pH值相差不大??傮w而言,銅山源水庫水體在4~8月普遍呈現(xiàn)偏堿性,其中5~8月(夏秋季節(jié))超過地表水環(huán)境質(zhì)量標(biāo)準(zhǔn),產(chǎn)生偏堿污染,全庫水質(zhì)pH值的空間差異性較小。

    2.2 自然環(huán)境本底水質(zhì)特征

    對銅山源水庫的入庫河道水質(zhì)進行監(jiān)測,以研究自然環(huán)境本底水質(zhì)對水庫偏堿污染的影響,結(jié)果如圖5所示。白鶴山入庫河道(St5監(jiān)測點附近)水體4~8月的pH值變化范圍在7~8之間,總體變化幅度較小。廟前溪和銅山溪水體的pH值變化幅度較大,其中,廟前溪水體pH值變化在7~9范圍內(nèi),銅山溪水體pH值變化在7.0~9.5范圍內(nèi);且銅山溪水體pH值在7~8月間總體維持在9以上的高位,已超過地表水水質(zhì)標(biāo)準(zhǔn)。由此可知,銅山源水庫附近自然河道中流動性本底水質(zhì)已呈現(xiàn)較為明顯的偏堿性,在7~8月一度超過地表水環(huán)境標(biāo)準(zhǔn)值上限。因此,當(dāng)河道中的水體進入水庫中成為水庫水體的一部分,必然導(dǎo)致水庫水體本底pH值偏高。

    2.3 關(guān)鍵氣象水文過程對庫區(qū)水體偏堿污染的影響

    2.3.1 氣象因素特征

    氣溫和光照等氣象因素會對地表水體水質(zhì)產(chǎn)生一定間接影響。一般而言,氣溫的影響主要通過水溫反映,光照主要通過影響綠色植物的光合作用影響水質(zhì)狀況,而植物光合作用一般采用水體中的溶解氧濃度來表征。因此,本研究主要通過庫區(qū)水體的水溫以及溶解氧濃度來分析氣象因素對庫區(qū)水體偏堿污染的影響。

    由圖6可知,pH值和DO濃度曲線隨著時間變化基本呈現(xiàn)同步升降的趨勢,因此水體pH值和DO濃度兩項指標(biāo)之間存在較為明顯的正相關(guān)性?;疑P(guān)聯(lián)度分析計算表明,pH值和DO濃度之間的灰色關(guān)聯(lián)系數(shù)為0.89,具有較為顯著的關(guān)聯(lián)關(guān)系,水庫水體pH值受到DO濃度的影響。本次水生態(tài)調(diào)查結(jié)果顯示,銅山源水庫中不存在相當(dāng)規(guī)模的水生植物群落生長跡象,因此水庫水體DO濃度和pH值受到水庫中浮游植物(主要為藻類)生長的影響。

    從pH值和水溫曲線的變化趨勢可知,4~8月水體水溫呈現(xiàn)逐步上升趨勢,而同時水體pH值總體呈現(xiàn)波動上升趨勢,說明水溫和pH值具有同向變化的特征?;疑P(guān)聯(lián)度分析顯示:pH值和水溫之間的灰色關(guān)聯(lián)系數(shù)為0.83,兩者具有較為顯著的正相關(guān)關(guān)系。有研究指出,水溫變化總體可從3個方面影響水體中pH值變化:① 水溫升高一定程度上促進水生生物的呼吸作用[22],致使水體向偏酸性方向發(fā)展;② 在適當(dāng)范圍內(nèi)水溫逐步升高可促進增加水體中藻類的生物量和優(yōu)勢度,促進其光合作用,使水體向偏堿性方向發(fā)展[23~25];③ 水溫升高促進CO2氣體分子運動速率加快,容易自水面逸出,致使水體向偏堿性方向發(fā)展[26]。由此可知,水溫對于水體pH值的影響主要通過促進水體中藻類生長繁殖,影響光合作用和呼吸作用的強度來實現(xiàn)。

    從各特征時間段來看:① 4月至5月中下旬,水體水溫、DO與pH值3項指標(biāo)逐漸同步升高,表明藻類的光合作用不斷增強、水體pH值也不斷升高;此時期內(nèi)水溫處于20~25 ℃范圍內(nèi),水生態(tài)調(diào)查發(fā)現(xiàn)水庫中優(yōu)勢藻種依次為4月初的硅藻門小環(huán)藻和5月的藍藻門微囊藻,表明隨著水溫的變化,水庫中適應(yīng)性的藻類群落開始生長演替。② 6月上中旬,水溫由25 ℃逐步升高至30 ℃,適宜藍綠藻生長,同時DO濃度充足,藍綠藻類生長旺盛、光合作用較強,pH值呈現(xiàn)高位波動態(tài)勢。③ 7月中旬,水體水溫在30~32 ℃范圍內(nèi)變動,過高的水溫一定程度上抑止了藻類生長,pH值與DO出現(xiàn)同步降低趨勢。④ 7月中下旬至8月期間,水溫維持在27~30 ℃范圍內(nèi),適宜藻類生長繁殖,此時期內(nèi)水體pH值逐漸出現(xiàn)波動上升趨勢。Reynolds等人研究指出,環(huán)境因子對浮游生物類型季節(jié)性豐度和多樣性的影響差異顯著,其中,溫度、光照強度等物理因子較為重要[27]。因此,從水溫、DO與水體pH值的變化關(guān)系分析可知,水溫和DO濃度通過影響藻類的生長間接影響水體pH值變化趨勢。王斌[28]、吳阿娜[29]等亦指出,湖庫水體葉綠素a與水溫、pH值、溶解氧呈顯著正相關(guān)。

    2.3.2 關(guān)鍵水文過程特征

    降雨以及徑流等關(guān)鍵水文過程同樣可能對水體水質(zhì)產(chǎn)生影響。本研究統(tǒng)計了2016~2021年5~8月水庫降雨、出入庫水量情況,如表1所示。由表1可知,2021年5~8月水庫周邊累積降雨量為1 206.90 mm,較2016~2020年同期累計降雨量均值(920.38 mm)多。2021年5~8月累積入庫徑流量14 634.50萬m3,較2016~2020年同期入庫徑流量均值(11 168.61萬m3)大;累積出庫徑流量17 744.32萬m3,與2016~2020年同期出庫徑流量均值(17 709.59萬m3)相當(dāng)。從凈出庫水量來看,2021年5~8月累積出庫凈流量3 109.82萬m3,較2016~2020年同期累積出庫凈流量(6 540.97萬m3)減少了52.5%。因此,2021年5~8月相比于2016~2020年同期,總體存在著庫區(qū)周邊降雨較往年偏多、入庫水量較往年偏多、出庫水量與往年持平、凈出庫水量較往年大幅減少的現(xiàn)象。

    多項研究指出,出入庫徑流以及湖庫本身水動力條件顯著影響浮游植物的生長。浮游植物的豐度及結(jié)構(gòu)動態(tài)變化與溫度、光照、水體穩(wěn)定性、水力滯留時間等因素時空變化的結(jié)果相關(guān)[30-32]。也有研究指出,與銅山源水庫位于相同緯度帶的浙江省新安江水庫,浮游植物群落結(jié)構(gòu)在很大程度上受水溫、氣象條件影響[33]。水力停留時間影響浮游植物生長[1~2],在光合作用強烈的靜水或流速小于0.3 m/s的水體中,光合作用引起的溶解氧變化,必然引起水質(zhì)pH值變化[34]。銅山源水庫出庫水量偏少的現(xiàn)象容易造成水庫水體流動緩慢,延長了水庫中水流的水力停留時間,適宜的水溫、光照、溶解氧等因素為水體中藻類的生長繁殖提供有利條件,從而促進了藻類的生長和光合作用強度,影響水庫水體的酸堿性,與上述研究結(jié)果一致。

    2.4 藻類生長狀況及其對水庫水質(zhì)偏堿污染的影響

    從關(guān)鍵氣象水文因素分析可知,藻類生長繁殖所引起的光合作用強度遠超水生動植物呼吸作用強度,最終導(dǎo)致水體中DO濃度上升,可能是水庫偏堿污染現(xiàn)象產(chǎn)生的原因之一。2020年,生態(tài)環(huán)境部太湖流域東海海域生態(tài)環(huán)境監(jiān)督管理局對浙閩片的湖庫進行了水環(huán)境水生態(tài)情況調(diào)查,發(fā)現(xiàn)在藻類大量生長的4~9月,浙江省部分水庫存在pH值超標(biāo)的現(xiàn)象[11]。目前,國內(nèi)外已有多項研究發(fā)現(xiàn),水體中的pH值和藻類生長情況具有顯著相關(guān)性[10,24-25,35-36]。有多名學(xué)者提出可以采用pH值指標(biāo)來作為藻類生長和繁殖的指示和預(yù)警指標(biāo)。浮游植物主要受光照、營養(yǎng)物質(zhì)、溫度、群落結(jié)構(gòu)、生活史、分層或垂直混合和潮汐等因素的影響[37]。本研究擬從藻類生物量、水體營養(yǎng)物質(zhì)含量以及藻類光合作用強度等方面探究藻類生長對水庫水體偏堿污染的影響。

    2.4.1 藻類生物量

    4~8月,庫區(qū)水體葉綠素a(Chl.a)、藻類生物量以及pH值隨時間變化特征呈現(xiàn)較為明顯的相似性和同步性(圖7)?;疑P(guān)聯(lián)度分析計算表明,pH值和葉綠素a濃度之間的灰色關(guān)聯(lián)系數(shù)為0.63,具有一定程度的正相關(guān)關(guān)系,由此可知水體偏堿污染與水體中藻類生長狀況密切相關(guān)。4~5月中旬,葉綠素a濃度與pH值同步波動上升,此時藻類生物量為1.59 mg/L;5月下旬pH值局部降低,藻類生物量和葉綠素a同步降低為0.84 mg/L和5.5 μg/L;6~7月,水溫長期保持在25~30 ℃,藻類生物量增至3.06~3.91 mg/L,葉綠素a濃度增至7~13 μg/L,此時pH值達到9.0~9.6;8月中下旬,pH值上升至9.4~9.8,此時水體中葉綠素a濃度也同步升至12~18 μg/L,上升幅度較大。因此,銅山源水庫中藻類群落的生物量對水體pH值產(chǎn)生影響。

    此外,監(jiān)測了3條入庫河道中葉綠素a濃度變化趨勢,見圖8。由圖可知,白鶴山入庫河道葉綠素a濃度在監(jiān)測期內(nèi)一直極??;廟前溪和銅山溪水體中葉綠素a濃度在7~8月顯著升高,尤其是銅山溪水體中葉綠素a濃度在8月末達到最高值23 μg/L,說明7~8月份銅山溪和廟前溪河道水體出現(xiàn)藻類快速繁殖現(xiàn)象,在一定程度上,可能對水庫水體中藻類生長繁殖產(chǎn)生一定程度的累加作用。

    從銅山源水庫藻類群落演替過程來看,4月硅藻為優(yōu)勢種,5月水庫中藍藻為優(yōu)勢種,6月藍藻和綠藻相當(dāng),7月開始則是綠藻為優(yōu)勢種。由此可知,水庫中的綠藻在6月初步開始繁殖,至7~8月已成為能夠提供大量葉綠素a用于光合作用的藻類群體,因此也極大地促進了7~8月份水體pH值的堿性發(fā)展。

    2.4.2 營養(yǎng)物質(zhì)

    一般而言,增加的營養(yǎng)負荷被認為是藻類尤其是藍藻優(yōu)勢的驅(qū)動因素[38],一定條件下營養(yǎng)物質(zhì)相比于溫度等氣象因素,對藻類生長產(chǎn)生的影響更大[39]。藻類在湖泊中大量繁殖,導(dǎo)致水體環(huán)境中光合作用強度大于呼吸作用,促進水體偏堿發(fā)展。4~8月,TP和TN濃度與pH值變化關(guān)系曲線如圖9所示。

    由圖9可知,pH值和TP濃度隨時間的變化曲線具有較為明顯的升降趨勢同步性,同時灰色關(guān)聯(lián)度計算顯示,pH值和TP、TN濃度之間的灰色關(guān)聯(lián)系數(shù)分別為0.67和0.69,具有一定程度的正向相關(guān)關(guān)系??赏茢嗫偭诪樗畮熘械母∮沃参镉绕涫窃孱愄峁I養(yǎng)物質(zhì),間接對水庫水體的酸堿性產(chǎn)生影響。同時,庫區(qū)水體的TP濃度較低,總體維持在Ⅱ類水以內(nèi),僅有局部監(jiān)測時段在Ⅱ~Ⅲ類水之間。由此可知,水庫中TP濃度并未完全滿足藻類生長的需求,隨著TP濃度的升降,藻類生長(葉綠素a指標(biāo))表現(xiàn)出極其同步的升降變化,TP指標(biāo)可認為是銅山源水庫中藻類生長的限制性指標(biāo)。

    由圖9可知,pH值和TN濃度隨時間的變化曲線之間并未產(chǎn)生明顯相關(guān)性。王斌等人[28]亦指出,葉綠素濃度與TN濃度相關(guān)性較小。同時,水庫的TN濃度在Ⅲ類標(biāo)準(zhǔn)以內(nèi)(相比TP指標(biāo),TN含量較高),可知TN濃度對于此階段內(nèi)的藻類生長尚未產(chǎn)生限制作用,充足的氮源為藻類生長提供了豐富的養(yǎng)料。

    研究指出,降雨和徑流中的養(yǎng)分對湖泊、水庫的養(yǎng)分負荷產(chǎn)生一定影響[40~42]。銅山源水庫上游地區(qū)的外源污染主要通過入庫河道匯入,監(jiān)測3條入庫河道總磷和總氮濃度變化趨勢如圖10所示。由圖可知,廟前溪水體TP濃度遠大于銅山溪和白鶴山入庫河道。對于總氮指標(biāo)來說,銅山溪和廟前溪總體呈現(xiàn)波動趨勢,其濃度值變動范圍在0~1.43 mg/L范圍內(nèi)(Ⅰ~Ⅳ類),白鶴山入庫河道的總氮濃度在監(jiān)測時段內(nèi)普遍較高,最高時可達2.36 mg/L,屬劣Ⅴ類水質(zhì)。因此,從營養(yǎng)性物質(zhì)上來看,3條入庫河道中白鶴山的總氮污染較嚴(yán)重;廟前溪的總磷污染較嚴(yán)重。氮磷類營養(yǎng)物質(zhì)進入庫區(qū),會為水庫中藻類提供大量營養(yǎng),促進其生長和繁殖,從而導(dǎo)致水體偏堿。

    2.4.3 藻類光合作用強度

    有研究指出,在富營養(yǎng)化的地表水中,由于光合作用產(chǎn)生的O2量遠遠大于呼吸作用所需要的O2量,DO主要來源于藻類的光合作用,而光合作用引起的DO變化,必然引起水體pH值的變化[34]。有研究指出,藻類光合作用可快速減少CO2,使HCO-3的解離平衡偏移,促使水中HCO-3分解,OH-濃度或活度增加,pH值升高[4]。Trivedi[43]研究指出,湖庫水體高pH值可能是由于微觀和宏觀植被的高光合作用,使平衡向堿性方向移動。由前文分析可知,水庫水體pH值和DO總體呈現(xiàn)出消長同步的變化特征,兩者存在較為明顯的正相關(guān)關(guān)系。一般而言,在藻類生長旺盛的水體中,pH值和溶解氧因主要受藻類晝夜間交替進行的光合作用與呼吸作用的影響而呈現(xiàn)明顯的日變化規(guī)律,即pH值和溶解氧在早晨時最低,日出后整個白天內(nèi),從最小值逐漸增加,至日落前達到最大[17-18]。因此,選取水庫中3~4個監(jiān)測點開展了pH值和溶解氧兩個指標(biāo)的晝夜變化監(jiān)測,探究水庫中藻類光合作用對水體偏堿污染的影響。

    7月水庫各監(jiān)測點的pH值和溶解氧指標(biāo)晝夜變化曲線如圖11所示。雙橋入庫口(St1)、湖心區(qū)(St3)、廟前入庫口(St4)以及國控斷面(St6)的水體pH值和溶解氧指標(biāo)隨時間均呈現(xiàn)類正態(tài)分布的變化趨勢,且兩指標(biāo)的消長基本同步。其中雙橋入庫口附近水體在夜間20:00左右pH值和溶解氧指標(biāo)達到最高值,其余位置水體在傍晚16:00和17:00間達到最高,與藻類光合與呼吸作用晝夜變化的時間特征吻合程度較高。

    8月水庫各監(jiān)測點的pH值和溶解氧指標(biāo)晝夜變化曲線如圖11所示。由圖可知,8月份的pH值和溶解氧指標(biāo)晝夜變化情況與7月份極為相似,兩者都呈現(xiàn)較為明顯的先升后降情況,與藻類光合與呼吸作用晝夜變化的時間特征吻合程度較高。

    2.5 銅山源水庫水體偏堿污染形成機制

    通過對銅山源水庫的水文氣象條件、水體pH值指標(biāo)與其他水質(zhì)指標(biāo)之間的相關(guān)性分析可知,銅山源水庫水體偏堿問題是在多種因素共同影響下形成的一種綜合外在表現(xiàn)。首先,銅山源所在地區(qū)的大環(huán)境本底水質(zhì)呈現(xiàn)偏堿的現(xiàn)象,尤其是入庫河道銅山溪水質(zhì)呈現(xiàn)明顯超標(biāo)現(xiàn)象(pH>9),根據(jù)水庫多年統(tǒng)計數(shù)據(jù),銅山溪的年入庫徑流量在水庫總?cè)肓髁恐兴急戎馗哌_60%~70%,因此銅山溪河道水質(zhì)偏堿是水庫水體偏堿污染的原因之一。從2.3節(jié)和2.4節(jié)分析可知,銅山源水庫水體pH值與水溫、DO、葉綠素a、藻類生物量以及TP濃度等指標(biāo)均有著良好的正相關(guān)關(guān)系。有研究指出,水體中CO2含量會受到藻類光合作用、水生生物呼吸作用等的影響[17]。一般來說,水生生物呼吸作用產(chǎn)生CO2,而光合作用消耗水體中的CO2。光合作用所需的CO2一部分來自于大氣,一部分來自于水中HCO-3的解離(HCO-3=OH-+CO2),HCO-3在磷酸酐酶的催化作用下變化為H2CO3,然后經(jīng)脫水反應(yīng)生成CO2后被利用,在HCO-3變成H2CO3的過程中,電離出OH-。因此,光合作用使水體偏堿發(fā)展,而呼吸作用的結(jié)果則使水體向偏酸方向發(fā)展。一旦光合作用的效應(yīng)大于呼吸作用,則水體中堿性物質(zhì)的產(chǎn)生速率大于酸性物質(zhì)的產(chǎn)生速率,水體即發(fā)生偏堿污染現(xiàn)象。4~8月間銅山源水庫水體pH值、水溫、DO、TP濃度等指標(biāo)均與葉綠素a、藻類生物量呈現(xiàn)正相關(guān)關(guān)系,而且水庫中存在著較為穩(wěn)定的水動力條件,說明水庫中藻類生長繁殖較為旺盛,光合作用較強,影響了水體中pH值、DO等物化指標(biāo)的變化。這與游亮等人的研究成果一致,當(dāng)藻類數(shù)量增加到一定數(shù)量級時,其數(shù)量的多少和生命活動的旺盛程度對水體的pH值起主導(dǎo)作用[44]。監(jiān)測所得的7~8月pH值、溶解氧晝夜變化關(guān)系也恰好證明了這一點。而同時較高的pH值有利于水華藻類的生長,大量藻類繁殖又進一步提高pH值[45],形成循環(huán)反饋的機制,在一定程度上加重水質(zhì)偏堿污染。

    總體來說,銅山源水庫水體偏堿污染形成機制為:上游入庫河流本底水質(zhì)偏堿,同時夏季藻類生長、光合作用速度遠遠大于水生動植物的呼吸作用,導(dǎo)致水體pH值在夏秋季節(jié)進一步上升,使水體偏堿污染嚴(yán)重。

    3 結(jié) 論

    本文通過開展實地水質(zhì)監(jiān)測和水生態(tài)調(diào)查,揭示了自然環(huán)境本底水質(zhì)、關(guān)鍵氣象水文過程、外源污染以及藻類生長等對銅山源水庫水質(zhì)偏堿污染的影響機制,得出水庫水質(zhì)偏堿現(xiàn)象成因主要有兩個方面。

    (1) 銅山源水庫周邊自然環(huán)境內(nèi)的水體均呈現(xiàn)不同程度的堿性,本底水質(zhì)偏堿是造成水庫水質(zhì)偏堿的重要外因。

    (2) 藻類生長和繁殖所產(chǎn)生的光合作用效應(yīng)大于水體中水生動植物的呼吸作用,是造成水庫水質(zhì)偏堿的內(nèi)部原因。其中,春夏季節(jié)的水溫、光照、溶解氧、氮磷類營養(yǎng)物質(zhì)以及相對穩(wěn)定的水流條件等客觀因素,為藻類生長提供了穩(wěn)定的生長環(huán)境和充分的光合作用原料,有力促進了藻類生長繁殖,間接促進水庫水體水質(zhì)朝著偏堿的方向發(fā)展。

    近年來,中國地表湖庫出現(xiàn)不同程度的水體偏堿污染現(xiàn)象。本研究揭示的地表湖庫水體偏堿問題成因,可為出現(xiàn)此類污染現(xiàn)象的湖庫治理提供理論支撐。

    參考文獻:

    [1] 馬沛明,施練東,張俊芳,等.浙江湯浦水庫浮游植物季節(jié)演替及其影響因子分析[J].環(huán)境科學(xué),2016,37(12):4560-4569.

    [2] 朱海燕,戴學(xué)穎,王可玉.于橋水庫pH值變化原因分析及可調(diào)控措施.[J].水科學(xué)工程與技術(shù),2010(1):40-43.

    [3] 秦保平,孫韌,王德龍,等.引灤河道中水質(zhì)偏堿的起因研究[J].環(huán)境科學(xué)研究,1999,12(5):38-42.

    [4] 林影,徐家明,張廣譜,等.煙臺市夾河流域pH偏高的原因[J].山東環(huán)境,1995,2:12-13,45.

    [5] 施章薇.東風(fēng)西沙水庫原水pH值升高及出廠水鋁控制方法研究[J].建材與裝飾,2019(4):33-34.

    [6] 邊歸國.東牙溪水庫偏堿性污染的成因研究[J].福建環(huán)境,2003,20(6):15-17.

    [7] 李建民,毛小英.大王灘水庫溶解氧及pH變化規(guī)律與原因分析[J].廣東水利水電,2012(12):37-40.

    [8] 申開旭.漁洞水庫pH值堿化成因分析及控制[J].環(huán)境科學(xué)導(dǎo)刊,2011,30(1):47-49.

    [9] 王蘇民,竇鴻身.中國湖泊志[M].北京:科學(xué)出版社,1998.

    [10] 米瑋潔,俞建軍,陳暉,等.浙江飲用水源地水環(huán)境狀況及存在問題分析[J].人民長江,2012,43(13):95-98.

    [11] 生態(tài)環(huán)境部太湖流域東海海域生態(tài)環(huán)境監(jiān)督管理局.浙閩片主要湖庫基本情況調(diào)查與分析報告[R].上海:生態(tài)環(huán)境部太湖流域東海海域生態(tài)環(huán)境監(jiān)督管理局,2020.

    [12] 鐘衛(wèi)鴻,單劍鋒,薛浚,等.氮和磷對銅山源水庫優(yōu)勢藻生長影響實驗研究[J].環(huán)境污染與防治,2003,25(1):20-22.

    [13] 李洪遠.生態(tài)恢復(fù)的原理與實踐(第一版)[M].北京:化學(xué)工業(yè)出版社,2005:15-18.

    [14] 李錦秀,杜斌,孫以三.水動力條件對富營養(yǎng)化影響規(guī)律探討[J].水利水電技術(shù),2005,36(5):15-18.

    [15] 王志紅,崔福義,安全,等.pH與水庫水富營養(yǎng)化進程的相關(guān)性研究[J].給水排水,2004,30(5):37-41.

    [16] 胡艷娟.水庫藻類光合作用對pH和溶解氧的影響[J].黑龍江水利科技,2011,39(2):17-18.

    [17] 黃歲樑,臧常娟,杜勝藍,等.pH、溶解氧、葉綠素a之間相關(guān)性研究Ⅰ:養(yǎng)殖水體[J].環(huán)境工程學(xué)報,2011,5(6):1201-1208.

    [18] 黃歲樑,臧常娟,杜勝藍,等.pH、溶解氧、葉綠素a之間相關(guān)性研究Ⅱ:非養(yǎng)殖水體[J].環(huán)境工程學(xué)報,2011,5(8):1681-1688.

    [19] 國務(wù)院.水污染防治行動計劃[Z].北京:國務(wù)院,2015.

    [20] 國家環(huán)境保護總局,水和廢水監(jiān)測分析方法編委會.水和廢水監(jiān)測分析方法(第四版)[M].北京:中國環(huán)境科學(xué)出版社,2002,223-226,243-257,276-285.

    [21] 譚學(xué)瑞,鄧聚龍.灰色關(guān)聯(lián)分析:多因素統(tǒng)計分析新方法[J].統(tǒng)計研究,1995,12(3):46-48.

    [22] 祁峰,李曉東,趙艷紅,等.鹽度、光照度和溫度對漂浮剛毛藻光合作用的影響[J].大連水產(chǎn)學(xué)院學(xué)報,2008,23(5):382-386.

    [23] PEPERZAK L.Climate change and harmful algal blooms in the North Sea[J].Acta Oecologica-International Journal of Ecology,2003,24(增1):139-144.

    [24] PAERL H W,HUISMAN J.Blooms like it hot[J].Science,2008,320:57-58.

    [25] O′NEIL J,DAVIS T,BURFORD M,et al.The rise of harmful cyanobacteria blooms:the potential roles of eutrophication and climate change[J].Harmful Algae,2012,14:312-334.

    [26] 劉輝,姬泓巍,辛梅.膠州灣水體中的二氧化碳體系[J].海洋科學(xué),1998,1(6):44-47.

    [27] REYNOLDS,C S.Phytoplankton periodicity:The interactions of form,function and environmental variability[J].Freshwater Biology,1984,14(2):111-142.

    [28] 王斌,馬健,王銀亞,等.天山天池夏季葉綠素a的分布及富營養(yǎng)化特征研究[J].環(huán)境科學(xué),2015,36(7):2466-2471.

    [29] 吳阿娜,朱夢杰,湯琳,等.淀山湖藍藻水華高發(fā)期葉綠素a動態(tài)及相關(guān)環(huán)境因子分析[J].湖泊科學(xué),2011,23(1):67-72.

    [30] BECKER V,CAPUTO L,ORDóNEZ J,et al.Driving factors of the phytoplankton functional groups in a deep Mediterranean reservoir[J].Water Research,2010,44(11):3345-3354.

    [31] VOGT R J,SHARMA S,LEAVITT P R.Decadal regulation of phytoplankton abundance and water clarity in a large continental reservoir by climatic,hydrologic and trophic processes[J].Journal of Great Lakes Research,2015,41:81-90.

    [32] INEKE V G,JEROEN V W,KATLEEN V G,et al.Covariation between zooplankton community composition and cyanobacterial community dynamics in Lake Blaarmeersen(Belgium)[J].FEMS Microbiology Ecology,2008,63(2):222-237.

    [33] 笪文怡,朱廣偉,吳志旭,等.2002~2017年千島湖浮游植物群落結(jié)構(gòu)變化及其影響因素[J].湖泊科學(xué),2019,31(5):1320-1333.

    [34] 張澎浪,孫承軍.地表水體中藻類的生長對pH值及溶解氧含量的影響[J].中國環(huán)境監(jiān)測,2004,20(4):49-50.

    [35] 徐恩兵,孫健,韋雷.黑臭河道中的藻類光合作用對pH和溶解氧的影響.[J].環(huán)境工程,2018,36:78-79.

    [36] 許霄霄,王靖偉.四川省邛海水體富營養(yǎng)化評價及建議[J].水利水電快報,2021,42(8):83-87.

    [37] ALVAREZ-GNGORA C,HERRERA-SILVEIRA J A.Variations of phytoplankton community structure related to water quality trends in a tropical karstic coastal zone[J].Marine Pollution Bulletin,2005,52(1):48-60.

    [38] BROOKES J D,CAREY C C.Resilience to blooms[J].Science,2011,334:46-47.

    [39] THACKERAY S,JONES I,MABERLY S.Long-term change in the phenology of spring phytoplankton:species-specific responses to nutrient enrichment and climatic change[J].Journal of Ecology,2008,96:523-535.

    [40] JEPPESEN E,KRONVANG B,MEERHOFF M,et al.Climate change effects on runoff,catchment phosphorus loading and lake ecological state,and potential adaptations[J].Journal of Environmental Quality,2009,38:1930-1941.

    [41] JEPPESEN E,KRONVANG B,OLESEN J E,et al.Climate change effects on nitrogen loading from cultivated atchments in Europe:implications for nitrogen retention,ecological state of lakes and adaptation[J].Hydrobiologia,2011,663:1-21.

    [42] MOSS B,KOSTEN S,MEERHOF M,et al.Allied attack:climate change and eutrophication[J].Inland Waters,2011(1):101-105.

    [43] TRIVEDI,R.K.Limnology of three fresh water ponds in Manglore[C]∥National Symposium on Advances in limnology conservation of endangered fish species.Srinagar,Garhwal:Narendra Pub.House,1989,23-25.

    [44] 游亮,崔莉鳳,劉載文,等.藻類生長過程中DO、pH與葉綠素的相關(guān)性分析[J].環(huán)境科學(xué)與技術(shù),2007,30(9):42-44.

    [45] 申開旭,申時斌,宋大恩.漁洞水庫pH值異常成因分析及水資源保護建議[J].環(huán)??萍?,2014,20(1):34-38,45.

    (編輯:李 慧)

    Abstract:

    In recent years,reservoirs in the hilly areas of southeastern China have been polluted by meta-alkalization in varying degrees,and the pH values of some reservoirs have frequently exceeded the environmental quality standards of surface water during some periods of time.Taking Tongshanyuan Reservoir,a large (second) type reservoir in Quzhou City,Zhejiang Province as an example,water quality monitoring and water ecology investigation were conducted,and the correlation analysis among hydrology,meteorology,water quality and water ecology indicators were combined to reveal the influence mechanism of background water quality of the natural environment,the key meteorological and hydrological processes,the exogenous pollution and the growth of algae on the meta-alkalization pollution of the reservoir.The results showed that the meta-alkalization pollution in Tongshanyuan Reservoir had lasted for a long time and was more serious from April to August every year.The background water of the natural environment around the reservoir was alkaline.The pH values of the inflow rivers were as high as 802 to 847,which was higher than the surface water standard in some time.The pH values of water in reservoir showed a relatively positive correlation to water temperature,DO,nutrients(TP、TN),chlorophyll-a to a certain extent,which rose and fell in synchronization with the change of time.Changes in the pH value of the water body were influenced by a variety of factors,which were mainly manifested in the direct or indirect enhancement for algal photosynthesis,leading to a greater effect of photosynthesis than respiration of aquatic organisms,which resulted in meta-alkalization and pollution of the reservoir water body.

    Key words:

    reservoir pollution;pH;algae;DO;chlorophyll-a;Tongshanyuan Reservoir

    猜你喜歡
    水質(zhì)污染
    水質(zhì)抽檢豈容造假
    環(huán)境(2023年5期)2023-06-30 01:20:01
    什么是污染?
    什么是污染?
    一月冬棚養(yǎng)蝦常見水質(zhì)渾濁,要如何解決?這9大原因及處理方法你要知曉
    這條魚供不應(yīng)求!蝦蟹養(yǎng)殖戶、垂釣者的最愛,不用投喂,還能凈化水質(zhì)
    堅決打好污染防治攻堅戰(zhàn)
    堅決打好污染防治攻堅戰(zhàn)
    圖像識別在水質(zhì)檢測中的應(yīng)用
    電子制作(2018年14期)2018-08-21 01:38:16
    濟下水庫徑流水質(zhì)和垂向水質(zhì)分析及評價
    對抗塵污染,遠離“霾”伏
    都市麗人(2015年5期)2015-03-20 13:33:49
    97热精品久久久久久| 久久草成人影院| 日本色播在线视频| 久久人人精品亚洲av| 一个人看的www免费观看视频| 成人永久免费在线观看视频| 免费看日本二区| 国产真实乱freesex| 亚洲va在线va天堂va国产| 最近在线观看免费完整版| 免费搜索国产男女视频| 国模一区二区三区四区视频| 女人十人毛片免费观看3o分钟| 夜夜爽天天搞| 午夜福利在线观看吧| 99久久精品国产国产毛片| 一本久久中文字幕| 国模一区二区三区四区视频| 最后的刺客免费高清国语| 毛片女人毛片| 亚洲av第一区精品v没综合| 有码 亚洲区| 天堂动漫精品| 国产男人的电影天堂91| 91麻豆精品激情在线观看国产| 国产精品,欧美在线| 国产不卡一卡二| 久久精品国产清高在天天线| 日本精品一区二区三区蜜桃| 乱人视频在线观看| 免费看美女性在线毛片视频| 欧美日韩乱码在线| 欧美zozozo另类| 天天躁日日操中文字幕| 中国美女看黄片| 俺也久久电影网| 国产精品嫩草影院av在线观看 | 久久午夜福利片| 男人的好看免费观看在线视频| 欧美xxxx性猛交bbbb| 日韩 亚洲 欧美在线| 色视频www国产| 99热6这里只有精品| 午夜福利成人在线免费观看| 又爽又黄无遮挡网站| 欧美国产日韩亚洲一区| 亚洲人成伊人成综合网2020| 国产精品美女特级片免费视频播放器| 桃色一区二区三区在线观看| 人妻少妇偷人精品九色| 免费黄网站久久成人精品| 啦啦啦韩国在线观看视频| 久久精品国产亚洲av涩爱 | 观看美女的网站| 亚洲精品影视一区二区三区av| 欧美xxxx性猛交bbbb| 老师上课跳d突然被开到最大视频| 午夜精品一区二区三区免费看| 欧美最黄视频在线播放免费| 成人一区二区视频在线观看| 麻豆一二三区av精品| 午夜日韩欧美国产| 国产伦一二天堂av在线观看| 在线观看66精品国产| 22中文网久久字幕| 欧美bdsm另类| 99在线人妻在线中文字幕| 国产高清有码在线观看视频| 毛片女人毛片| 最新在线观看一区二区三区| 性欧美人与动物交配| 亚洲久久久久久中文字幕| 男女那种视频在线观看| 中文字幕高清在线视频| 熟妇人妻久久中文字幕3abv| 色综合色国产| 黄色配什么色好看| 黄色日韩在线| 国产乱人视频| 国产精品无大码| 美女高潮喷水抽搐中文字幕| 久久精品综合一区二区三区| 亚洲七黄色美女视频| 午夜免费激情av| 亚洲人与动物交配视频| 特大巨黑吊av在线直播| 99久久成人亚洲精品观看| 午夜免费男女啪啪视频观看 | 人妻丰满熟妇av一区二区三区| 中文字幕av在线有码专区| 在线国产一区二区在线| 婷婷精品国产亚洲av| 国产精品一区二区性色av| 亚洲真实伦在线观看| 精品国产三级普通话版| 免费高清视频大片| 免费观看人在逋| 99国产精品一区二区蜜桃av| 久久九九热精品免费| 岛国在线免费视频观看| 一个人观看的视频www高清免费观看| 亚洲av成人av| 亚洲最大成人av| 午夜福利在线观看吧| 三级毛片av免费| 精品久久久久久久久久久久久| 99热这里只有精品一区| 亚洲人成网站在线播| 99精品在免费线老司机午夜| 国产探花极品一区二区| 精品福利观看| 亚洲一区二区三区色噜噜| 亚洲成人中文字幕在线播放| 日韩欧美一区二区三区在线观看| 国产成年人精品一区二区| 校园人妻丝袜中文字幕| 男插女下体视频免费在线播放| 男人狂女人下面高潮的视频| 精品人妻一区二区三区麻豆 | 麻豆av噜噜一区二区三区| 黄色欧美视频在线观看| 国产伦一二天堂av在线观看| 国产熟女欧美一区二区| 人人妻人人澡欧美一区二区| 一进一出好大好爽视频| 五月玫瑰六月丁香| 岛国在线免费视频观看| 午夜激情欧美在线| 一级av片app| 日韩大尺度精品在线看网址| 网址你懂的国产日韩在线| 我的老师免费观看完整版| 成人三级黄色视频| 亚洲精品色激情综合| 69av精品久久久久久| 能在线免费观看的黄片| 国产高清三级在线| 国产毛片a区久久久久| 欧美又色又爽又黄视频| 亚洲av免费在线观看| 夜夜爽天天搞| 又爽又黄无遮挡网站| 亚洲午夜理论影院| 女同久久另类99精品国产91| 国产精品无大码| 男女边吃奶边做爰视频| 91久久精品国产一区二区成人| 成人av在线播放网站| 欧美成人免费av一区二区三区| 午夜视频国产福利| 永久网站在线| 亚洲内射少妇av| 69人妻影院| 久久久国产成人免费| 99久久成人亚洲精品观看| 亚洲午夜理论影院| 国产精品一区二区性色av| 日韩精品有码人妻一区| 亚洲18禁久久av| 搡老妇女老女人老熟妇| 国产欧美日韩一区二区精品| 久久精品国产鲁丝片午夜精品 | 欧美日韩中文字幕国产精品一区二区三区| 亚洲av第一区精品v没综合| 亚洲一区高清亚洲精品| 成人欧美大片| 国产老妇女一区| 国产精品无大码| 深夜精品福利| 成人午夜高清在线视频| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久免费视频| 日本熟妇午夜| 国产爱豆传媒在线观看| 久久久久久伊人网av| 狠狠狠狠99中文字幕| 国产成人av教育| 联通29元200g的流量卡| 欧美日韩综合久久久久久 | 欧美激情在线99| 亚洲欧美激情综合另类| 大又大粗又爽又黄少妇毛片口| 成人高潮视频无遮挡免费网站| 三级国产精品欧美在线观看| 久久这里只有精品中国| 精品免费久久久久久久清纯| 简卡轻食公司| 国产久久久一区二区三区| 小蜜桃在线观看免费完整版高清| 黄色日韩在线| 午夜影院日韩av| 男女边吃奶边做爰视频| 成人无遮挡网站| 啦啦啦啦在线视频资源| 免费一级毛片在线播放高清视频| 在线观看66精品国产| 九九久久精品国产亚洲av麻豆| 99九九线精品视频在线观看视频| 国产人妻一区二区三区在| 午夜福利成人在线免费观看| 小蜜桃在线观看免费完整版高清| 在线国产一区二区在线| 欧美xxxx性猛交bbbb| 美女被艹到高潮喷水动态| 成人国产综合亚洲| 国内精品美女久久久久久| 国产成人av教育| 欧美成人a在线观看| 毛片女人毛片| 18禁黄网站禁片午夜丰满| 人妻久久中文字幕网| 欧美激情国产日韩精品一区| 国产 一区精品| 日韩强制内射视频| 桃红色精品国产亚洲av| 午夜老司机福利剧场| 国产大屁股一区二区在线视频| 亚洲国产高清在线一区二区三| 免费电影在线观看免费观看| 日韩高清综合在线| 直男gayav资源| 熟女电影av网| 深夜精品福利| 久久久久九九精品影院| 午夜久久久久精精品| 亚洲中文字幕日韩| 免费电影在线观看免费观看| 中亚洲国语对白在线视频| 校园春色视频在线观看| 国产精品国产高清国产av| 热99在线观看视频| 夜夜爽天天搞| 亚洲avbb在线观看| 日韩精品中文字幕看吧| 午夜福利欧美成人| 欧美在线一区亚洲| 精品午夜福利视频在线观看一区| 韩国av一区二区三区四区| 日韩欧美国产一区二区入口| 国产精品三级大全| 简卡轻食公司| a级一级毛片免费在线观看| 国产高潮美女av| 国产精品爽爽va在线观看网站| 无人区码免费观看不卡| 日日摸夜夜添夜夜添小说| 两个人的视频大全免费| 亚洲国产欧洲综合997久久,| 1024手机看黄色片| 亚洲三级黄色毛片| 久久久久久九九精品二区国产| 少妇的逼好多水| 在线免费十八禁| 日韩一区二区视频免费看| 免费在线观看影片大全网站| 亚州av有码| 1024手机看黄色片| 国产aⅴ精品一区二区三区波| 国语自产精品视频在线第100页| 一a级毛片在线观看| 久久6这里有精品| 中文字幕免费在线视频6| 免费无遮挡裸体视频| 精品久久久久久久久亚洲 | 国产在线男女| 亚洲真实伦在线观看| 免费看美女性在线毛片视频| 最新在线观看一区二区三区| 日本 欧美在线| 欧洲精品卡2卡3卡4卡5卡区| 波野结衣二区三区在线| 香蕉av资源在线| 桃红色精品国产亚洲av| 99精品在免费线老司机午夜| 国产av在哪里看| 欧美黑人巨大hd| 午夜久久久久精精品| 日韩一本色道免费dvd| a级毛片a级免费在线| 成人二区视频| 两性午夜刺激爽爽歪歪视频在线观看| 美女免费视频网站| 嫩草影院入口| 一本精品99久久精品77| 999久久久精品免费观看国产| 国产精品人妻久久久影院| 99精品久久久久人妻精品| 免费人成在线观看视频色| 日本五十路高清| 免费搜索国产男女视频| 十八禁国产超污无遮挡网站| 999久久久精品免费观看国产| 琪琪午夜伦伦电影理论片6080| 无人区码免费观看不卡| 中文字幕av在线有码专区| 日本色播在线视频| 嫩草影院入口| 国产又黄又爽又无遮挡在线| 国产免费一级a男人的天堂| 久久精品综合一区二区三区| aaaaa片日本免费| 亚洲最大成人中文| 久久国产乱子免费精品| 亚洲黑人精品在线| 久久久久久久久久久丰满 | 国产伦精品一区二区三区视频9| 色精品久久人妻99蜜桃| 中文字幕免费在线视频6| h日本视频在线播放| 丰满的人妻完整版| 一本一本综合久久| 天堂网av新在线| 波多野结衣高清无吗| 欧美性猛交黑人性爽| 免费看日本二区| 欧美不卡视频在线免费观看| 日韩中字成人| 中国美女看黄片| 午夜精品一区二区三区免费看| 成人av一区二区三区在线看| 欧美一区二区国产精品久久精品| 一进一出抽搐动态| 亚洲无线观看免费| 黄片wwwwww| 一区二区三区四区激情视频 | 五月玫瑰六月丁香| 成年人黄色毛片网站| 国产精品野战在线观看| 欧美日本亚洲视频在线播放| 一夜夜www| 亚洲国产精品久久男人天堂| 免费观看在线日韩| 久久精品综合一区二区三区| 国产69精品久久久久777片| 亚洲专区国产一区二区| 有码 亚洲区| 久久久国产成人精品二区| 不卡视频在线观看欧美| av.在线天堂| 国产精品人妻久久久影院| 午夜爱爱视频在线播放| 欧美日韩国产亚洲二区| 我的老师免费观看完整版| 国产精品一及| 亚洲aⅴ乱码一区二区在线播放| 午夜激情欧美在线| 观看免费一级毛片| 国产午夜福利久久久久久| 亚洲av日韩精品久久久久久密| 免费看av在线观看网站| 亚洲专区国产一区二区| 国产精品自产拍在线观看55亚洲| 亚洲男人的天堂狠狠| 午夜福利高清视频| 狠狠狠狠99中文字幕| av在线老鸭窝| 亚洲av中文av极速乱 | 精品久久久久久久人妻蜜臀av| 搡老岳熟女国产| www.www免费av| 国产欧美日韩精品一区二区| 亚洲乱码一区二区免费版| 久久久久久久精品吃奶| 精品人妻一区二区三区麻豆 | 国产精品野战在线观看| 国产三级在线视频| 狠狠狠狠99中文字幕| 观看免费一级毛片| 神马国产精品三级电影在线观看| 国产精品精品国产色婷婷| 色av中文字幕| 一区二区三区高清视频在线| 国产黄色小视频在线观看| 亚洲第一电影网av| 久久久久久九九精品二区国产| 美女xxoo啪啪120秒动态图| 国产av一区在线观看免费| 亚洲专区中文字幕在线| 久久6这里有精品| 亚洲在线观看片| 国产一区二区三区av在线 | 国产精品一区二区性色av| 日韩欧美免费精品| 亚洲第一电影网av| 亚洲专区中文字幕在线| 亚洲欧美日韩卡通动漫| 啦啦啦韩国在线观看视频| 久久精品国产自在天天线| 成年人黄色毛片网站| 国产免费男女视频| 亚洲熟妇中文字幕五十中出| 亚洲欧美精品综合久久99| 99久久久亚洲精品蜜臀av| 国产伦人伦偷精品视频| 久久久精品大字幕| 国产真实伦视频高清在线观看 | 亚洲av五月六月丁香网| 精品久久久久久久久av| 老司机福利观看| 国产高清视频在线播放一区| 少妇裸体淫交视频免费看高清| 乱码一卡2卡4卡精品| 男女边吃奶边做爰视频| 亚洲成人久久爱视频| 成人国产综合亚洲| 国产国拍精品亚洲av在线观看| x7x7x7水蜜桃| 国产乱人伦免费视频| 女生性感内裤真人,穿戴方法视频| 精品99又大又爽又粗少妇毛片 | 国产精品久久视频播放| 波多野结衣巨乳人妻| 色综合亚洲欧美另类图片| x7x7x7水蜜桃| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 亚洲第一电影网av| 成人国产一区最新在线观看| 精品久久久久久久末码| 成人一区二区视频在线观看| 一级a爱片免费观看的视频| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品av在线| 国产精品人妻久久久影院| 老司机福利观看| 亚洲专区中文字幕在线| 性色avwww在线观看| 搡老熟女国产l中国老女人| 亚洲自偷自拍三级| 欧美另类亚洲清纯唯美| 91午夜精品亚洲一区二区三区 | 国产精品乱码一区二三区的特点| 88av欧美| aaaaa片日本免费| 在线观看免费视频日本深夜| 欧美性猛交╳xxx乱大交人| 久久国产精品人妻蜜桃| 色在线成人网| 给我免费播放毛片高清在线观看| 少妇熟女aⅴ在线视频| 国国产精品蜜臀av免费| 国内久久婷婷六月综合欲色啪| 国产三级在线视频| 国内精品宾馆在线| 亚洲最大成人中文| 男女边吃奶边做爰视频| 亚洲不卡免费看| 午夜福利18| 日韩亚洲欧美综合| 麻豆成人av在线观看| 日本爱情动作片www.在线观看 | 天堂网av新在线| 99久久精品热视频| 91麻豆av在线| 波多野结衣高清无吗| 亚洲四区av| 欧美激情久久久久久爽电影| 精品久久国产蜜桃| 国内久久婷婷六月综合欲色啪| www日本黄色视频网| 欧美色欧美亚洲另类二区| 免费在线观看影片大全网站| 少妇猛男粗大的猛烈进出视频 | 此物有八面人人有两片| www.www免费av| 搞女人的毛片| a级毛片a级免费在线| videossex国产| 国产欧美日韩一区二区精品| 国产亚洲精品综合一区在线观看| 国产av在哪里看| 精品一区二区三区视频在线| 熟女电影av网| 亚洲精品久久国产高清桃花| 国产白丝娇喘喷水9色精品| 在线观看美女被高潮喷水网站| 99久久精品一区二区三区| 亚洲18禁久久av| 国产 一区 欧美 日韩| 国内揄拍国产精品人妻在线| 日本黄色片子视频| 欧美又色又爽又黄视频| 久久午夜亚洲精品久久| 波野结衣二区三区在线| 岛国在线免费视频观看| 中文资源天堂在线| 国国产精品蜜臀av免费| 啦啦啦观看免费观看视频高清| 久久精品国产清高在天天线| 欧美极品一区二区三区四区| 亚洲自拍偷在线| 麻豆精品久久久久久蜜桃| 日韩人妻高清精品专区| 成年人黄色毛片网站| 国产黄片美女视频| 丝袜美腿在线中文| 中国美白少妇内射xxxbb| 国产精品三级大全| 啦啦啦观看免费观看视频高清| 亚洲自偷自拍三级| 国产 一区 欧美 日韩| 国产蜜桃级精品一区二区三区| 亚洲avbb在线观看| 午夜福利欧美成人| 国产69精品久久久久777片| 午夜福利高清视频| 成年人黄色毛片网站| 看十八女毛片水多多多| 舔av片在线| 国产精品电影一区二区三区| 国产成年人精品一区二区| 黄色欧美视频在线观看| 亚洲精品色激情综合| 国产aⅴ精品一区二区三区波| 国产高清激情床上av| 日本在线视频免费播放| 日韩精品有码人妻一区| 色av中文字幕| 少妇丰满av| 又爽又黄无遮挡网站| 亚洲美女黄片视频| 美女高潮的动态| av中文乱码字幕在线| 99热6这里只有精品| 在线观看av片永久免费下载| 国产精华一区二区三区| 小说图片视频综合网站| 欧美成人a在线观看| 黄色日韩在线| 亚洲精华国产精华精| 亚洲va日本ⅴa欧美va伊人久久| 深爱激情五月婷婷| 人人妻人人看人人澡| 欧美日韩中文字幕国产精品一区二区三区| 成人综合一区亚洲| 亚洲av成人av| 人妻丰满熟妇av一区二区三区| 色吧在线观看| 午夜久久久久精精品| 97超视频在线观看视频| 国产又黄又爽又无遮挡在线| 亚洲一级一片aⅴ在线观看| av在线老鸭窝| 免费不卡的大黄色大毛片视频在线观看 | 婷婷精品国产亚洲av| 18禁在线播放成人免费| 在线国产一区二区在线| 国产一区二区在线av高清观看| 欧美潮喷喷水| 亚洲精品粉嫩美女一区| 窝窝影院91人妻| 国产一区二区亚洲精品在线观看| 成人精品一区二区免费| 欧美性猛交黑人性爽| 日韩欧美在线乱码| 国产精品自产拍在线观看55亚洲| 在线观看舔阴道视频| 久久久精品欧美日韩精品| 国产成人一区二区在线| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 成人午夜高清在线视频| 伦理电影大哥的女人| 性插视频无遮挡在线免费观看| 国产精品一区二区三区四区免费观看 | 国产视频一区二区在线看| 午夜免费激情av| 999久久久精品免费观看国产| 久久久国产成人免费| 神马国产精品三级电影在线观看| 欧美丝袜亚洲另类 | 可以在线观看的亚洲视频| 亚洲av中文av极速乱 | 欧美日本亚洲视频在线播放| 色av中文字幕| 久久久久久久久大av| 精品久久久噜噜| 亚洲精品粉嫩美女一区| av在线亚洲专区| 亚洲国产精品成人综合色| 国产探花极品一区二区| 我的老师免费观看完整版| 人妻少妇偷人精品九色| 色在线成人网| 女人被狂操c到高潮| 国产精品免费一区二区三区在线| 久久99热6这里只有精品| 国产精品av视频在线免费观看| 久久久久精品国产欧美久久久| АⅤ资源中文在线天堂| 麻豆av噜噜一区二区三区| 99热网站在线观看| 日日夜夜操网爽| 22中文网久久字幕| 欧美另类亚洲清纯唯美| 亚洲国产欧洲综合997久久,| 老司机深夜福利视频在线观看| 久久午夜福利片| 男女啪啪激烈高潮av片| 亚洲av不卡在线观看| 亚洲精品色激情综合| 国产精品亚洲美女久久久| 校园人妻丝袜中文字幕| 全区人妻精品视频| xxxwww97欧美| 精品一区二区免费观看| 国产精品一区二区免费欧美| 亚洲最大成人av| 黄色日韩在线| 久99久视频精品免费| 五月伊人婷婷丁香| 日本一本二区三区精品| 亚洲性夜色夜夜综合| 国产乱人伦免费视频| 午夜精品久久久久久毛片777|