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

    淤泥質(zhì)潮灘重金屬和溶解性有機(jī)質(zhì)的潮汐變化特征及其環(huán)境影響評價

    2022-12-21 06:28:26張麗聰肖凱張鵬李海龍王俊堅李鎮(zhèn)揚王芬芳徐華林郭躍華
    生態(tài)環(huán)境學(xué)報 2022年11期
    關(guān)鍵詞:潮灘潮位潮汐

    張麗聰,肖凱*,張鵬,李海龍,王俊堅,李鎮(zhèn)揚,王芬芳,徐華林,郭躍華

    1. 南方科技大學(xué)環(huán)境科學(xué)與工程學(xué)院,廣東 深圳 518055;2. 中國地質(zhì)大學(xué)(北京)水資源與環(huán)境學(xué)院,北京 100083;3. 廈門大學(xué)環(huán)境與生態(tài)學(xué)院,福建 廈門 361102;4. 廣東內(nèi)伶仃福田國家級自然保護(hù)區(qū)管理局,廣東 深圳 518040;5. 廣東省中交一航局生態(tài)工程有限公司,廣東 深圳 518107

    淤泥質(zhì)潮灘指富含泥質(zhì)沉積物的潮灘,位于平均大潮高潮線與平均大潮低潮線之間的地帶,可分為光灘和沼澤帶(楊世倫,2003)。在世界范圍內(nèi),淤泥質(zhì)潮灘廣泛分布于緯度 73°以下的一些海灣和河口等處。中國潮灘岸線總長約為4000 km(王穎,1990),具有規(guī)模大、分布廣、類型多的特點(陳思明,2018)。潮灘具有多種生態(tài)價值,如過濾污染物、為動植物提供棲息地和提高海岸的防護(hù)能力等(Goodwin,2001;Wang et al.,2012;Leonardi et al.,2016)。但近幾十年來,隨著沿海經(jīng)濟(jì)的發(fā)展,海岸帶出現(xiàn)了不同程度的重金屬污染問題。

    泥質(zhì)潮灘作為海陸交互帶,它們是環(huán)境相關(guān)溶質(zhì)(如重金屬、碳和營養(yǎng)鹽)活躍的生物地球化學(xué)轉(zhuǎn)化區(qū),由于涉及含水層,其溶質(zhì)成分復(fù)雜。水體中重金屬濃度和溶解態(tài)有機(jī)質(zhì)(Dissolved Organic Matter,DOM)特征受到物理驅(qū)動、一系列復(fù)雜的水巖相互作用和生物地球化學(xué)作用而表現(xiàn)出復(fù)雜的變化(O’Connor et al.,2018;Liu et al.,2019;Zhang et al.,2022)。這些作用包括載體上離子的吸附和解吸、有機(jī)物質(zhì)的再礦化、氧化還原反應(yīng)產(chǎn)生或消耗金屬氧化物等(李海龍等,2015)。潮灘水體中重金屬如果過度富集不僅會影響動植物的生長,污染食物鏈,而且會影響海洋重金屬的均衡(Marchand et al.,2016;Trezzi et al.,2016)。DOM是全球碳循環(huán)的重要組成部分,在水生生態(tài)系統(tǒng)的生物地球化學(xué)過程發(fā)揮著重要的作用,但同樣具有生物毒性和遺傳毒性(He et al.,2015;Tang et al.,2020)。有色溶解有機(jī)質(zhì)(Colored Dissolved Organic Matter,CDOM)是指在紫外-可見波長范圍內(nèi)可以吸收和重新發(fā)射光的DOM部分,可有效反映DOM的來源、分子特性(芳香性和分子量)或微生物成分等(Zhou et al.,2010;Pain et al.,2019)。已有大量的研究表明,有機(jī)質(zhì)的再礦化作用與絡(luò)合作用均會對重金屬產(chǎn)生影響(Roy et al.,2011;Liu et al.,2019)。

    潮灘沉積物和地下水中的重金屬和有機(jī)碳可以通過地表水-地下水直接的對流與擴(kuò)散作用,及底棲動物洞穴沖刷輸送到地表水中(Burdige,2006;Kalnejais et al.,2015;Xiao et al.,2021;Pan et al.,2022)。海底地下水排泄(Submarine Groundwater Discharge,SGD)已被認(rèn)為是碳、重金屬、營養(yǎng)鹽向海排泄的重要途徑。大量研究表明,SGD攜帶的溶解態(tài)重金屬和有機(jī)碳在量級上和河流輸送的通量相當(dāng),有些甚至遠(yuǎn)超過河流通量(Chen et al.,2018;Wang et al.,2019)。當(dāng)具有不同組分的地下水和海水混合時,潮灘內(nèi)的地表水-地下水對流和擴(kuò)散過程可能會增加或減少重金屬和DOM濃度,進(jìn)而影響污染物在沉積物中的遷移轉(zhuǎn)化過程(Beck et al.,2008;Seidel et al.,2014;O′Connor et al.,2018)。但是鮮有研究同時關(guān)注重金屬和DOM特征在潮汐尺度上的動態(tài)特征及相互作用機(jī)制。另外,目前對于DOM的研究主要集中于地表水體中或者地表河口區(qū)域,對于潮灘淺層地下水的研究也較少(Pain et al.,2019)。

    綜上所述,選取漳江口的紅樹林濕地自然保護(hù)區(qū)的淤泥質(zhì)光灘作為研究場地,在潮汐尺度內(nèi)收集地下水和地表水樣本,以研究重金屬(Cr、Cu、Zn、As、Ba、Pb和Hg)濃度和有機(jī)質(zhì)特征隨潮汐的變化規(guī)律和驅(qū)動機(jī)制,定量評估重金屬和有機(jī)質(zhì)在沉積物-水界面的對流和擴(kuò)散通量,評估其生態(tài)環(huán)境效應(yīng)和污染風(fēng)險。研究成果有助于增強(qiáng)對淤泥質(zhì)潮灘水動力過程及其對地下水重金屬和有機(jī)質(zhì)循環(huán)影響的了解,有助于潮灘特殊生境的保護(hù)和有效管理。

    1 材料與方法

    1.1 研究區(qū)概況

    研究區(qū)位于福建省云霄縣漳江口紅樹林國家自然保護(hù)區(qū)(117°24'07″—117°30'00″E,23°53'45″—23°56'00″N),具體位置如圖1所示。該保護(hù)區(qū)由河口水域、紅樹林濕地、光灘和鹽沼組成,總面積為2360.0 hm2。該區(qū)域?qū)賮啛釒ШQ笮约撅L(fēng)氣候,夏熱冬溫,年均降雨量為1714.5 mm。潮汐類型屬于不規(guī)則半日潮,平均潮差為2.32 m(Wang et al.,2022)。監(jiān)測點位于自然保護(hù)區(qū)的淤泥質(zhì)潮灘,與潮溝相鄰。沉積物以淤泥和砂質(zhì)淤泥為主,厚度超過2 m,其中含有大量植物殘體、鈣質(zhì)、軟體動物殘體和石灰物質(zhì),有機(jī)質(zhì)含量大多在2.5%—10.0%范圍內(nèi)(薛博,2007)。

    圖1 漳江口的地理位置和監(jiān)測點的分布示意圖(改自Wang et al., 2022)Figure 1 Locations of the Zhangjiang Estuary and monitoring point (modified from Wang et al., 2022)

    1.2 野外工作方法

    1.2.1 地下水動態(tài)監(jiān)測

    監(jiān)測點安裝有“對井”裝置,可以實現(xiàn)對兩個不同深度地下水的水壓、電導(dǎo)率和溫度的自動監(jiān)測(Hou et al.,2016)。圖2為對井裝置結(jié)構(gòu),其由上下兩個相同的井筒和一根固定的鏈接桿組成,井筒有內(nèi)外套筒,外套筒筒壁上均勻分布著直徑為 5 mm的透水孔,內(nèi)筒中放置了傳感器(CTD-Divers)。監(jiān)測時間為2020年10月11日13:00—20:00,頻率為60分鐘/次。

    圖2 對井裝置的組成Figure 2 The schematic structure diagrams showing

    1.2.2 水樣采集

    從低潮到高潮期間,地下水每隔2小時采集1次,共采集3次;地表水每隔1小時采集1次,共采集8次(圖3)。為減小采樣活動對監(jiān)測的地下水位產(chǎn)生影響,采樣點選取在“對井”裝置2 m附近。地下水采集是使用挖坑滲水法收集 0.5 m深的水樣;地表水樣品使用Niskin采水器采集(水面以下50 cm)(Xiao et al.,2022)。水樣的理化參數(shù)包括鹽度、pH和溶解氧(DO),使用便攜式多參儀(Hanna,HI 9829)進(jìn)行現(xiàn)場測量。水樣過0.45 μm醋酸纖維濾膜保存于預(yù)先經(jīng)過特殊處理的樣品瓶中。重金屬樣品加入HNO3酸化至pH為1.5—2.0。有機(jī)質(zhì)樣品置于棕色玻璃瓶中,4 ℃下保存,待測。

    圖3 地下水監(jiān)測周期與采樣頻率Figure 3 Groundwater monitoring period and sampling frequency

    1.3 室內(nèi)實驗

    水樣中的Cr、Cu、Zn、As、Ba和Pb使用電感耦合等離子體質(zhì)譜儀(ICP-MS,Thermo Fisher Scientific,USA)測量,Hg采用原子熒光光譜儀測量。為減小測量誤差,每份樣品均分別3次測樣后取平均質(zhì)量濃度。水樣中溶解性有機(jī)碳質(zhì)量濃度(DOC)通過總有機(jī)碳分析儀(TOC-LCSH/CSN)配備ASI-L(Shimadzu,Japan)自動進(jìn)樣裝置檢測,樣品檢測前,通過滴入鹽酸將樣品酸化至pH=2左右,去除水樣中無機(jī)碳干擾。水樣中DOM的光譜特征采用紫外-可見光分光光度計(Shimadzu,Japan)測定,以Milli-Q超純水作為空白水樣。掃描波長范圍為200—800 nm,波段步長為1 nm。本文共考慮了 SUVA254、E2/E3、a254、a350和a3655種DOM光譜特征,具體計算方法和含義描述見表1。

    表1 溶解性有機(jī)質(zhì)光譜特征指數(shù)含義Table 1 Meaning of spectral characteristic index of chromophoric dissolved organic matter

    1.4 計算方法

    1.4.1 沉積物-水界面的交換通量

    本研究計算的重金屬和 DOM 在沉積物-水界面的交換通量包括水力梯度驅(qū)動的對流通量和濃度梯度驅(qū)動的擴(kuò)散通量。對流通量采用“對井法”(Hou et al.,2016)計算,根據(jù)廣義達(dá)西定律,首先計算地下水-地表水交換速率qn:

    式中:

    Kv——沉積物垂向滲透系數(shù),利用原位降水頭法測定(Wang et al.,2014),其值為 1.2×10-7m·s-1;

    J——上下兩個單井的水力梯度,計算方法為:

    式中:

    hup和hlow——上下井的等效淡水水頭;

    ΔL——上下井的間距,值為0.4 m;

    cup和clow——上、下井處地下水的鹽度,利用監(jiān)測的電導(dǎo)率和溫度數(shù)據(jù)獲得(Cox et al.,1967);

    φ——淡水與海水的動態(tài)黏度之比;

    ε——常數(shù),為 7.143×10-4m3·kg-1,用于描述密度和鹽度間的線性關(guān)系。

    若qn>0,表示為地下水排泄;qn<0,則表示地表水入滲到含水層中。

    對流攜帶的重金屬和DOM通量F可按如下公式估算(Qu et al.,2017):

    式中:

    此外,本研究還考慮了垂向分子擴(kuò)散通量,根據(jù)菲克第一定律,目標(biāo)溶質(zhì)在地下水與地表水之間的擴(kuò)散通量(F0)可表示為(Fick,1855):

    式中:

    φ0——沉積物在深度為1—10 cm和20—30 cm孔隙度的平均值,為0.4。

    Ds——溶質(zhì)的分子擴(kuò)散系數(shù),當(dāng)沉積物的孔隙度φ0<0.7 時為(Ullman et al.,1987):

    式中:

    D0——目標(biāo)溶質(zhì)在無限稀釋條件下的擴(kuò)散系數(shù)。不同金屬(Cr、Cu、Zn、As、Ba、Pb 和 Hg)和 DOC 的 D0分別為 7.95×10-6、6.59×10-6、6.86×10-6、 6.56×10-6、 1.84×10-6、 8.41×10-6、14.00×10-6、100.00×10-6cm2·s-1(Poisson et al.,1983;Dharmaprakash et al.,1989;Kuss et al.,2009;Wang et al.,2016;Yao et al.,2020)。

    1.4.2 重金屬污染評價

    采用單因子污染指數(shù)法(Pi)和綜合污染指數(shù)法(Psyn)(Hakanson,1980;Sun et al.,2020)評估研究區(qū)地下水和地表水中重金屬的污染狀況,具體公式如下:

    式中:

    Ci——第i種金屬的平均質(zhì)量濃度;

    Si——第i種金屬評價標(biāo)準(zhǔn),地表水采用《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB 3838—2002)中的Ⅰ類標(biāo)準(zhǔn),對地下水采用國土資源部修訂的《地下水質(zhì)量標(biāo)準(zhǔn)》(GB/T 14848—2017)的Ⅰ類標(biāo)準(zhǔn)。根據(jù) Pi值可將水質(zhì)可分為四類:清潔(Pi<1,Ⅰ類)、輕度污染(1<Pi<3,Ⅱ類)、中度污染(3<Pi<6,Ⅲ類),和重度污染(Pi>6,Ⅳ類)(Luo et al.,2022)。同樣,根據(jù) Psyn值可將水質(zhì)分為清潔(Psyn≤1)、輕度污染(1<Psyn≤2)、中度污染(2<Psyn≤3)和重度污染(Psyn>3)4個等級(Sun et al.,2020)。

    1.4.3 統(tǒng)計分析

    采用單因素方差分析(ANOVA)和 Tukey事后檢驗研究地表水和地下水中相關(guān)物理化學(xué)參數(shù)的差異性,當(dāng)顯著性水平P<0.01時為差異極顯著;P<0.05時為顯著。同時,對于不滿足單因素方差分析數(shù)據(jù),需通過對數(shù)轉(zhuǎn)換以滿足方差齊性要求。采用Mantel test分析和Pearson相關(guān)分析研究相關(guān)參數(shù)的相關(guān)性。統(tǒng)計分析及數(shù)據(jù)可視化利用 R語言(R3.6.0軟件)完成,示意圖和統(tǒng)計圖分別在CorelDRAW 2018和Sigmaplot 14.0軟件中完成。

    2 結(jié)果

    2.1 理化參數(shù)、重金屬和有機(jī)質(zhì)的潮汐變化

    在取樣期間,地表水中的鹽度、pH和DO總體呈現(xiàn)先下降后上升的變化趨勢,在低潮時達(dá)到最低值,在漲潮過程中呈增大趨勢(圖4)。在地下水中,鹽度和pH較穩(wěn)定,DO隨潮水位的增加而迅速降低,從 7.01 mg·L-1減少到 2.82 mg·L-1。此外,地表水中的pH和DO均高于地下水,且存在顯著差異(pH:P=0.0004;DO:P=0.002),而鹽度差異性不顯著(P=0.720)。

    圖4 潮灘地表水和地下水中理化參數(shù)(鹽度,pH和DO)的潮汐變化Figure 4 Tidal variations of physicochemicla parameters (salinity, pH and DO) in the surface water and groundwater in the tidal mudflat

    地表水中重金屬平均質(zhì)量濃度的大小順序為:Ba>Zn>Pb>As>Cu>Cr>Hg;地下水中的為:Ba>Zn>As>Cu>Cr>Pb>Hg。在均值水平上,地表水中 Cr、Zn、As和Hg與地下水中的質(zhì)量濃度很接近,但地表水中的Cu、Ba和Pb分別是地下水中3.4、2.5和9.2倍。在取樣期間,地表水和地下水中重金屬的質(zhì)量濃度隨潮汐的變化較復(fù)雜(圖5)。除As和Pb,地表水中的重金屬均在低潮位時達(dá)到最大值。除Pb外,地表水中的重金屬質(zhì)量濃度隨著潮位的增加總體呈現(xiàn)先上升后下降的變化趨勢。地下水中Ba和Cu的變化較穩(wěn)定,而Cr、As和Hg表現(xiàn)出先增加后減少的趨勢,Zn表現(xiàn)出持續(xù)增加的趨勢。根據(jù)方差分析,地表水和地下水中的As和Ba存在顯著差異(As:P=0.049;Ba:P=0.038),而 Cr、Cu、Zn、Pb和 Hg不存在顯著差異(Cr:P=0.569;Cu:P=0.332;Zn:P=0.900;Pb:P=0.563;Hg:P=0.925)。

    圖5 潮灘地表水和地下水中重金屬質(zhì)量濃度的潮汐變化Figure 5 Tidal variations of heavy metals in surface water and groundwater in the tidal mudflat

    在取樣期間,地表水中的 DOC質(zhì)量濃度、CDOM的光譜特征指數(shù)隨潮汐的變化明顯,通常在低潮時達(dá)到最大值,并隨著潮汐水位的增加而降低。對比,在地下水中這些指標(biāo)隨著潮汐水位的增加而呈現(xiàn)出增加的趨勢(圖6)。地表水和地下水中DOC的平均值分別為 (3.11±0.29) mg·L-1和 (3.32±0.16) mg·L-1。根據(jù)方差分析,地表水和地下水中DOC、CDOM紫外光譜特征(SUVA254,a350)均存在顯著差異(SUVA254:P=0.014;a350:P=0.003),而DOC質(zhì)量濃度和溶解有機(jī)質(zhì)相對分子量(E2/E3)不存在顯著差異(DOC:P=0.275;E2/E3:P=0.132)。

    圖6 潮灘地表水和地下水中DOC和CDOM光譜特征指數(shù)的潮汐變化Figure 6 Tidal variations of the DOC and spectral characteristic indexes of CDOM in the surface water and groundwater in the tidal mudflat

    2.2 沉積物-水界面水交換速率

    根據(jù)公式(1),計算得到監(jiān)測點處地下水-地表水交換速率,結(jié)果如圖7所示。地下水-地表水交換速率具有明顯的周期性變化,但由于水力梯度和滲透系數(shù)較小,導(dǎo)致整體交換速率較低,平均地下水排泄速率和地表水入滲速率分別為 3.22×10-5m·h-1和2.10×10-5m·h-1,凈排泄速率為1.12×10-5m·h-1。在取樣期間,地下水的排泄速率從4.75×10-5m·h-1(13:00)減小到0(17:00),之后隨著潮位的增加,交換方式轉(zhuǎn)變?yōu)楹K霛B(2.16×10-6—5.40×10-5m·h-1)。

    圖7 監(jiān)測期間地下水位和地下水-地表水交換速率Figure 7 Time series of groundwater level and groundwater-surface water exchange rates during the monitoring period

    2.3 沉積物-水界面重金屬和DOC的交換通量

    根據(jù)公式(3)、(4),計算得到的重金屬和DOC在監(jiān)測期內(nèi)通過沉積物-水界面的交換通量如圖 8所示。地下水交換攜帶的重金屬排泄量的順序為:Ba>Zn>As>Cu>Cr>Pb>Hg,其中地下水中大部分重金屬表現(xiàn)為向地表水凈排泄,即潮灘沉積物為重金屬的源;而Cu和Pb的凈排泄量小于0,表現(xiàn)為匯。擴(kuò)散通量大小順序為:Ba>Zn>Pb>As>Hg>Cr>Cu,其中Cr、Cu、Hg表現(xiàn)為由地表水向地下水?dāng)U散,而其他重金屬表現(xiàn)出向地表水?dāng)U散。除了Pb的擴(kuò)散通量約為地下水排泄攜帶的通量的20%外,其他重金屬通過地下水排泄攜帶的對流通量是擴(kuò)散通量的20—200倍。因此,除了Pb外,其他重金屬的擴(kuò)散通量相比于對流通量可以忽略不計。此外,地下水中 DOC 的 SGD 通量為 0.11 mg·m-2·h-1,擴(kuò)散通量為 3.45×10-3mg·m-2·h-1,相差 31 倍,同樣可以忽略。

    圖8 重金屬和DOC在沉積物-水界面的對流和擴(kuò)散通量Figure 8 Calculated connective and diffusive fluxes for heavy metals and DOC

    2.4 重金屬污染評價

    根據(jù)單因子污染指數(shù)法(Pi),漳江口的潮灘水體現(xiàn)出了一定程度的污染(除 Cr、Cu和 Pb外)(圖9)。由于缺乏Ba相應(yīng)的水質(zhì)標(biāo)準(zhǔn),且其生態(tài)危害性較小,因此不考慮Ba的污染狀況。其中,Hg的污染最嚴(yán)重,地表水和地下水中Hg的平均Pi值分別達(dá)到了中度和輕度污染水平,地表水中的最大值甚至達(dá)到了重度污染水平。此外,Zn和As也表現(xiàn)出了潛在的污染風(fēng)險。除As和Hg在地表水和地下水中污染情況差距較大外,其他重金屬在兩種水體中的污染水平相當(dāng)。但依據(jù)綜合污染指數(shù)法(Psyn),地表水和地下水中的 Psyn分別為 3.90和1.98,屬于重度污染和輕度污染水平??傮w上,地表水中7種重金屬的污染程度大于地下水,說明潮灘具有潛在的凈化重金屬污染物的作用。

    圖9 地下水和地表水中重金屬的污染指數(shù)Figure 9 The pollution index of heavy metals in the groundwater and surface water

    3 討論

    3.1 重金屬的潮汐變化的影響因素

    本研究基于Mantel測試檢驗了潮汐因素1(潮位和鹽度)和化學(xué)因素2(pH和DO)對水體中重金屬和有機(jī)質(zhì)的影響(圖10)。結(jié)果表明:除因素2與As具有顯著相關(guān)性外(P<0.05),因素1和2對其他重金屬的影響均較小。但是,隨著潮位的升高,地表水和地下水中的大部分重金屬均表現(xiàn)出顯著的變化(圖5),這是因為在地表水與地下水交換過程中,其中的重金屬與有機(jī)質(zhì)含量的變化并非受單因素控制,而是物理混合、一系列復(fù)雜的水巖相互作用和生物地球化學(xué)反應(yīng)共同作用的結(jié)果(李海龍等,2015)。

    圖10 環(huán)境因子的相關(guān)性Figure 10 Correlations between environmental factors

    在地表水中,Cr、Cu、Zn、Ba和Hg的質(zhì)量濃度最大值通常出現(xiàn)在低潮期間,這表明其可能來源于地下水的排泄或者上游河水的輸入。地表水中的Cu和Ba的質(zhì)量濃度明顯大于地下水,但監(jiān)測期間發(fā)生海水-地下水交換時地下水中這兩種重金屬的變化較小,且Cu和Ba與潮位、鹽度、pH和DO均不相關(guān),這可能是因為地表水與地下水交換所憑借的環(huán)境介質(zhì)——潮灘沉積物具有一定程度的阻隔地表水中的Cu和Ba侵入地下水的能力。先前研究表明,沉積物具有接受、儲存和循環(huán)水、營養(yǎng)物質(zhì)及能量的能力(Anderson et al.,1984),此外,還具有保留、分散和轉(zhuǎn)化化學(xué)物質(zhì)和生物的能力,從而起到環(huán)境過濾器或緩沖區(qū)的作用(Gregorich et al.,1994)。潮灘是由陸地生態(tài)系統(tǒng)向海洋生態(tài)系統(tǒng)的陸海交界過渡地帶,可作為陸源性物質(zhì)進(jìn)入海洋之前的過濾器(Ferreira et al.,2017)。因此,地表水與地下水中溶質(zhì)質(zhì)量濃度差別較大,且隨著地表水和地下水混合其含量變化不明顯時,地表水與地下水交換的環(huán)境介質(zhì)——潮灘沉積物在其中發(fā)揮了重要的作用。此外,Cr和 Cu呈極顯著正相關(guān)(r=0.775),表明二者可能為同一來源。特別是上游紅樹林地下水中重金屬的輸入,Cr、Zn、Cu和Pb在上游的紅樹林區(qū)具有較高的含量(Hong et al.,2021;Wang et al.,2022)。因此,上游河水的輸入將對Cr、Zn、Cu、Pb和Ba的含量產(chǎn)生重要的影響。在相同的取樣時刻,特別是在最低潮的取樣,地下水中的Zn、As和Hg通常高于地表水中的質(zhì)量濃度。因此,地下水排放攜帶的Zn、As和Hg將對地表水中的重金屬動態(tài)產(chǎn)生重要影響。地表水中Pb在高潮時達(dá)到最大值,也可能受到海水輸入的影響。

    隨著潮位的增加,地表水中的重金屬(除 Pb外)與地下水中Cr、As和Hg的質(zhì)量濃度呈現(xiàn)先增大后減小的趨勢,這與漲潮期間先發(fā)生地下水排泄后發(fā)生地表水入滲的兩個時期(圖7)相對應(yīng)。在地下水排泄時期,地下水中的Cr、As和Hg的質(zhì)量濃度增大可能是因為地下水的擾動使沉積動力環(huán)境更活躍,使沉積物中的Cr、As和Hg暴露在氧化或次氧化環(huán)境中,從而促進(jìn)其釋放到水體中,質(zhì)量濃度增大。而在海水入滲期間,金屬質(zhì)量濃度減小則主要是因為稀釋作用。圖4中DO的質(zhì)量濃度變化趨勢證實了這一點,隨著潮位的升高,地表水中的DO質(zhì)量濃度增大,但地下水中的DO持續(xù)且迅速下降,表明地下水中具有較多的氧化還原反應(yīng)、有機(jī)物質(zhì)再礦化和金屬氧化物等消耗DO的反應(yīng)發(fā)生。已有大量研究表明,潮間帶地下水中Fe和Mn的氧化還原循環(huán)會對其他重金屬產(chǎn)生重要的影響(O′Connor et al.,2018;Liu et al.,2019)。特別是漳江口沉積物中Fe和Mn元素含量遠(yuǎn)高于其他重金屬(Yang et al.,2020)。例如,鐵(氫)氧化物可以有效地清除地下水中的Cr、As和Hg(Bone et al.,2006;Richard et al.,2016)。當(dāng)?shù)叵滤優(yōu)檫€原環(huán)境時,鐵(氫)氧化物溶解釋放沉積物附著的重金屬,這可能是地下水中Cr、As和Hg質(zhì)量濃度增加的原因。在15:00之后,隨著地表水的入滲,Cr、Cu、As和Hg的質(zhì)量濃度開始減少,這主要可能是地表水的稀釋引起的(圖5)。

    3.2 有機(jī)質(zhì)的潮汐變化及影響因素

    如圖10所示,根據(jù)Mantel’s測試的結(jié)果,DOC和DOM的光譜特征參數(shù)受潮汐因素1(鹽度和潮位)的影響較低,但是DOM的光譜特征參數(shù)受化學(xué)因素2(pH和DO)的影響較明顯,這表明水體中的DOM含量受生物和化學(xué)反應(yīng)作用(因素2)的影響更強(qiáng)烈,物理作用(因素1)影響較小。

    地表水中DOC在最低潮時達(dá)到最大值,并隨著潮位的增加而逐漸降低(圖6)。這表明地下水排泄是地表水中 DOC的重要來源。此外,DOC與CDOM、CDOM與SUVA254具有相似的變化趨勢,且具有顯著的正相關(guān)(DOC與CDOM:r=0.646;CDOM與SUVA254:r=0.914),表明地下水中高芳香性有機(jī)質(zhì)在有機(jī)質(zhì)中占比較高。但隨著潮位的上升,DOM 的分子量趨于穩(wěn)定,這是因為海水中含有較多的小分子量有機(jī)質(zhì)(Richard et al.,2016)。在地下水中,DOC、CDOM、SUVA254和E2/E3隨著潮位的上升表現(xiàn)出增加的趨勢。pH、DO與CDOM和SUVA254表現(xiàn)出了顯著的負(fù)相關(guān)(pH與CDOM:r= -0.603;pH與SUVA254:r= -0.692;DO與CDOM:r= -0.696;DO與CDOM:r= -0.801;),這表明高芳香性DOM主要來源于地下水,并隨著咸淡水混合、潮汐波流等因素使沉積物頻繁地再懸浮、移動和干擾生物活動等,已沉積的有機(jī)質(zhì)會不斷地從沉積物釋放到水體中。新產(chǎn)生的DOM具有高芳香化程度和低分子量的特點,與低潮時地表水中 DOM的特征相同,更加證明了地下水的排泄對地表水中DOM的影響。

    3.3 重金屬和DOC的排泄通量和環(huán)境效應(yīng)

    相較于國內(nèi)外其他的研究場地,漳江口泥質(zhì)潮灘 SGD攜帶的重金屬通量通常處于較低的水平,但對于毒性較強(qiáng)的As和Pb,其通量分別是中國膠州灣鹽沼和河口濕地的34倍和21倍,而Zn和Pb通量最高分別可達(dá)地中海西北部的100倍和20倍(表2)。此外,Cu和Pb的凈通量小于0,表明其會在沉積物中不斷累積。因此,漳江口紅樹林濕地淤泥質(zhì)潮灘是Cu和Pb的匯,是其他重金屬和有機(jī)質(zhì)的源。地下水中的DOC質(zhì)量濃度(3.32 mg·L-1)明顯高于茅尾海(1.44 mg·L-1)和珍珠灣紅樹林生態(tài)系統(tǒng)(0.72—2.28 mg·L-1)中的,但SGD輸送的DOC通量卻遠(yuǎn)遠(yuǎn)小于這兩個研究場地(王亞麗,2020)。較小的重金屬和DOC通量主要是由于較小的SGD速率引起的。

    表2 世界范圍內(nèi)對流和擴(kuò)散輸送的重金屬通量的對比Table 2 Comparisons of heavy metal fluxes transported by convection and diffusion globally

    從污染水平看,漳江口淤泥質(zhì)潮灘地表水和地下水分別處于重度污染和輕度污染水平,其中 Hg的污染最嚴(yán)重,Zn和As表現(xiàn)出潛在的污染風(fēng)險。并且,上游SGD輸送的重金屬及DOM可能對潮灘和海洋中重金屬及碳均衡產(chǎn)生重要的影響,在相同的時間,Wang et al.,(2022)計算了泥質(zhì)潮灘上游的紅樹林和鹽沼區(qū)的SGD速率分別為0.18 m·d-1和0.04 m·d-1,分別是本研究的23、5倍,這是因為淤泥質(zhì)潮灘的滲透系數(shù)(Kv)遠(yuǎn)小于紅樹林和鹽沼區(qū)的滲透系數(shù)。淤泥質(zhì)潮灘滲透性差,同時受不同植被類型的影響。此外,紅樹林區(qū)具有富含有機(jī)質(zhì)的特征,重金屬Cr、Zn、Cu和Pb的含量也最高(Wang et al.,2022)。

    4 結(jié)論

    本研究詳細(xì)調(diào)查和分析了漳江口淤泥質(zhì)潮灘地下水和地表水中重金屬和有機(jī)質(zhì)在潮汐尺度的變化特征、重金屬和有機(jī)質(zhì)在沉積物-水界面的交換通量、重金屬污染評價,得到的主要結(jié)論如下:

    (1)地表水中的重金屬最大值通常出現(xiàn)在低潮期間,并隨著潮位的增加而減小。其動態(tài)主要受上游河水輸入或地下水排泄控制。地下水中Cr、Zn、As和 Hg受生物地球化學(xué)反應(yīng)和地表水稀釋的共同影響。地表水中DOC和DOM特征參數(shù)隨潮位的增加而減少,在地下水中相反。地下水中的DOM,具有大分子量和高芳香性的特點,主要受pH和DO的影響。統(tǒng)計分析表明,由于潮灘復(fù)雜的環(huán)境和水文過程,地表水和地下水中有機(jī)質(zhì)與重金屬的相關(guān)性較低。

    (2)相較于對流通量,重金屬和DOC的擴(kuò)散通量可以忽略不計。SGD攜帶的重金屬和DOC的通量范圍分別是 0.19—155 μg·m-2·d-1和 0—2.56 mg·m-2·d-1,而擴(kuò)散通量范圍分別為-8.30×10-4—1.72 μg·m-2·d-1和 0—8.30×10-2mg·m-2·d-1。除 As和Pb外,其余重金屬和DOC的交換通量較小,但這主要是由于較低的SGD速率引起的。

    (3)重金屬污染評價結(jié)果表明,漳江口淤泥質(zhì)潮灘地表水和地下水中重金屬的綜合污染程度分別處于重度污染和輕度污染水平。其中,地表水中的Hg處于中度污染水平,是地表水污染的最主要因素;地下水中的Hg和As處于輕度污染水平,是地下水污染的主要因素;其余重金屬均處于清潔水平;As表現(xiàn)出潛在的污染風(fēng)險。

    猜你喜歡
    潮灘潮位潮汐
    杭州灣北岸臨港新片區(qū)岸段潮灘沉積特征
    海洋通報(2022年5期)2022-11-30 12:03:48
    潮汐與戰(zhàn)爭(上)
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    崇明東北潮灘泥沙輸運過程的研究
    海洋通報(2021年1期)2021-07-23 01:55:40
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    潮灘干濕轉(zhuǎn)換的地貌發(fā)育物理模型及動力機(jī)制
    海洋通報(2020年3期)2020-10-20 09:26:46
    絕美海灘
    多潮位站海道地形測量潮位控制方法研究
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測研究
    潮汐式灌溉控制系統(tǒng)的設(shè)計及應(yīng)用
    電子制作(2017年9期)2017-04-17 03:00:56
    91久久精品国产一区二区三区| 又大又黄又爽视频免费| 国产 一区精品| 一本—道久久a久久精品蜜桃钙片| 色视频在线一区二区三区| 日韩亚洲欧美综合| 亚洲av成人精品一区久久| 亚洲人成77777在线视频| 午夜日本视频在线| 亚洲欧洲精品一区二区精品久久久 | 久久久久视频综合| 亚洲精品亚洲一区二区| 插阴视频在线观看视频| 这个男人来自地球电影免费观看 | 午夜福利视频在线观看免费| 成人二区视频| 女性被躁到高潮视频| 99久久精品国产国产毛片| 日韩 亚洲 欧美在线| 久久国内精品自在自线图片| 能在线免费看毛片的网站| 午夜久久久在线观看| 亚洲国产精品一区二区三区在线| 大陆偷拍与自拍| 色5月婷婷丁香| 少妇的逼水好多| 日日啪夜夜爽| 夜夜骑夜夜射夜夜干| 22中文网久久字幕| 欧美日本中文国产一区发布| 一级毛片aaaaaa免费看小| 精品99又大又爽又粗少妇毛片| a级毛片黄视频| 最黄视频免费看| 尾随美女入室| 国产69精品久久久久777片| 国产伦精品一区二区三区视频9| 只有这里有精品99| 免费看不卡的av| 亚洲精品中文字幕在线视频| 插阴视频在线观看视频| 青春草亚洲视频在线观看| 看非洲黑人一级黄片| 天堂中文最新版在线下载| 国产又色又爽无遮挡免| 国产色婷婷99| 18禁观看日本| 91精品国产九色| 国产免费一区二区三区四区乱码| 久久99蜜桃精品久久| 在线观看人妻少妇| 久久久久久久国产电影| 黄色怎么调成土黄色| 日韩强制内射视频| 五月伊人婷婷丁香| 精品少妇黑人巨大在线播放| 一区二区三区乱码不卡18| 亚洲精品久久午夜乱码| 97超视频在线观看视频| 人妻夜夜爽99麻豆av| 一个人看视频在线观看www免费| 亚洲精品国产av蜜桃| 搡老乐熟女国产| 中文天堂在线官网| 两个人免费观看高清视频| 熟女人妻精品中文字幕| 青春草亚洲视频在线观看| 一级,二级,三级黄色视频| 色网站视频免费| 久久精品人人爽人人爽视色| 久久毛片免费看一区二区三区| 一个人免费看片子| 亚洲少妇的诱惑av| 天天影视国产精品| 夜夜爽夜夜爽视频| 国产成人午夜福利电影在线观看| 韩国av在线不卡| 男人添女人高潮全过程视频| 91久久精品电影网| 中国国产av一级| 亚洲精品成人av观看孕妇| 精品一品国产午夜福利视频| 这个男人来自地球电影免费观看 | 老女人水多毛片| 国产精品.久久久| 久久精品久久久久久噜噜老黄| 校园人妻丝袜中文字幕| 自线自在国产av| 人妻制服诱惑在线中文字幕| 日本午夜av视频| 国产成人av激情在线播放 | av网站免费在线观看视频| 国产一区亚洲一区在线观看| 欧美精品亚洲一区二区| 蜜桃在线观看..| 成人毛片a级毛片在线播放| 亚洲综合色网址| 日韩亚洲欧美综合| 午夜免费观看性视频| 欧美精品一区二区大全| 亚洲国产欧美在线一区| 日韩亚洲欧美综合| 国产成人午夜福利电影在线观看| 91国产中文字幕| 欧美少妇被猛烈插入视频| 色婷婷久久久亚洲欧美| 中国美白少妇内射xxxbb| 亚洲精品,欧美精品| 99九九线精品视频在线观看视频| 晚上一个人看的免费电影| 婷婷成人精品国产| 最近手机中文字幕大全| 热99国产精品久久久久久7| 日韩中文字幕视频在线看片| 多毛熟女@视频| 少妇精品久久久久久久| 丰满迷人的少妇在线观看| 精品视频人人做人人爽| 久久久欧美国产精品| 成人免费观看视频高清| 交换朋友夫妻互换小说| 国产免费一级a男人的天堂| 一级爰片在线观看| 人人妻人人爽人人添夜夜欢视频| 女人久久www免费人成看片| 99国产精品免费福利视频| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| 亚洲内射少妇av| 哪个播放器可以免费观看大片| 国产精品国产三级专区第一集| 在线看a的网站| 日本与韩国留学比较| 这个男人来自地球电影免费观看 | 亚洲伊人久久精品综合| 色哟哟·www| 国产老妇伦熟女老妇高清| 国产一区有黄有色的免费视频| 黄色配什么色好看| 亚洲av中文av极速乱| 精品久久蜜臀av无| 欧美另类一区| a 毛片基地| 国产视频内射| 久久精品国产亚洲网站| 日本91视频免费播放| 黄色配什么色好看| 国产有黄有色有爽视频| 亚洲精品456在线播放app| 午夜激情福利司机影院| 亚洲国产色片| 性色av一级| 狠狠精品人妻久久久久久综合| 自线自在国产av| 一级二级三级毛片免费看| 国产成人精品婷婷| 亚洲国产av影院在线观看| 欧美人与性动交α欧美精品济南到 | 成人午夜精彩视频在线观看| 99热这里只有是精品在线观看| 成人亚洲欧美一区二区av| 97在线人人人人妻| 韩国av在线不卡| 欧美成人精品欧美一级黄| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 日韩三级伦理在线观看| 免费av不卡在线播放| 国产白丝娇喘喷水9色精品| 我的老师免费观看完整版| 成人毛片a级毛片在线播放| 少妇丰满av| 卡戴珊不雅视频在线播放| 少妇 在线观看| 另类精品久久| 一级爰片在线观看| 国产精品一区二区在线观看99| 九草在线视频观看| 精品熟女少妇av免费看| 最新的欧美精品一区二区| 美女脱内裤让男人舔精品视频| 熟女电影av网| 满18在线观看网站| 99热国产这里只有精品6| 亚洲欧洲日产国产| 免费av中文字幕在线| 又粗又硬又长又爽又黄的视频| 人人妻人人澡人人看| 蜜臀久久99精品久久宅男| 精品视频人人做人人爽| 五月开心婷婷网| 免费观看a级毛片全部| 国产片内射在线| 一本一本综合久久| 新久久久久国产一级毛片| 三上悠亚av全集在线观看| 久久久国产欧美日韩av| 欧美激情国产日韩精品一区| 亚洲av欧美aⅴ国产| 国产免费一级a男人的天堂| 一区在线观看完整版| 欧美少妇被猛烈插入视频| 狂野欧美白嫩少妇大欣赏| 又黄又爽又刺激的免费视频.| 欧美日韩一区二区视频在线观看视频在线| 亚洲一级一片aⅴ在线观看| 久久影院123| 免费观看性生交大片5| 亚洲欧美清纯卡通| 日韩一区二区三区影片| 日韩视频在线欧美| 人成视频在线观看免费观看| 国产成人精品久久久久久| 最近最新中文字幕免费大全7| 日本av手机在线免费观看| 晚上一个人看的免费电影| 午夜福利视频在线观看免费| 欧美激情国产日韩精品一区| 2022亚洲国产成人精品| 中国三级夫妇交换| 搡老乐熟女国产| 春色校园在线视频观看| 18禁观看日本| 九九爱精品视频在线观看| 国产精品熟女久久久久浪| 免费观看在线日韩| 69精品国产乱码久久久| 一级片'在线观看视频| 国产日韩欧美视频二区| 最近手机中文字幕大全| 亚洲成人一二三区av| 久热久热在线精品观看| av黄色大香蕉| 在线精品无人区一区二区三| 国产成人精品无人区| 国产欧美日韩综合在线一区二区| 七月丁香在线播放| 18禁裸乳无遮挡动漫免费视频| 人人妻人人添人人爽欧美一区卜| 日本与韩国留学比较| 99re6热这里在线精品视频| 夜夜看夜夜爽夜夜摸| 看非洲黑人一级黄片| 插阴视频在线观看视频| 日韩电影二区| 欧美日韩在线观看h| 国产免费一区二区三区四区乱码| 国产色爽女视频免费观看| 男的添女的下面高潮视频| av天堂久久9| 97精品久久久久久久久久精品| 成人18禁高潮啪啪吃奶动态图 | 久久精品夜色国产| 三上悠亚av全集在线观看| 如何舔出高潮| 99re6热这里在线精品视频| 边亲边吃奶的免费视频| 少妇熟女欧美另类| 久久97久久精品| 在现免费观看毛片| 美女视频免费永久观看网站| 午夜免费观看性视频| 汤姆久久久久久久影院中文字幕| 视频中文字幕在线观看| 内地一区二区视频在线| 日本-黄色视频高清免费观看| 免费久久久久久久精品成人欧美视频 | a 毛片基地| 在现免费观看毛片| 特大巨黑吊av在线直播| 国产欧美亚洲国产| 爱豆传媒免费全集在线观看| 成年人午夜在线观看视频| 国产日韩欧美在线精品| 人成视频在线观看免费观看| 亚洲国产欧美在线一区| 成人无遮挡网站| 亚洲色图综合在线观看| 99国产综合亚洲精品| 亚洲av.av天堂| 2021少妇久久久久久久久久久| 少妇人妻精品综合一区二区| 美女中出高潮动态图| 91久久精品电影网| 五月玫瑰六月丁香| 人人澡人人妻人| 午夜激情福利司机影院| 中文字幕免费在线视频6| 成年人免费黄色播放视频| 久久人妻熟女aⅴ| 久久久久精品久久久久真实原创| av天堂久久9| 女性被躁到高潮视频| 精品久久久噜噜| 国产一区二区三区综合在线观看 | 另类精品久久| 老女人水多毛片| 国产av精品麻豆| 日韩不卡一区二区三区视频在线| 美女xxoo啪啪120秒动态图| 久久久久国产精品人妻一区二区| 精品人妻在线不人妻| 亚洲av欧美aⅴ国产| 亚洲综合色网址| 女性被躁到高潮视频| 日本欧美视频一区| av.在线天堂| 成人手机av| 永久网站在线| av网站免费在线观看视频| 免费av中文字幕在线| 免费看av在线观看网站| 热99久久久久精品小说推荐| 国产熟女欧美一区二区| 日产精品乱码卡一卡2卡三| 如日韩欧美国产精品一区二区三区 | 亚洲性久久影院| 日本午夜av视频| 国模一区二区三区四区视频| 国产高清三级在线| 天天操日日干夜夜撸| 国产视频首页在线观看| 女性被躁到高潮视频| 香蕉精品网在线| 亚洲欧美中文字幕日韩二区| 嫩草影院入口| 午夜福利视频在线观看免费| 久久国产精品男人的天堂亚洲 | 看免费成人av毛片| 午夜激情久久久久久久| 久久久精品区二区三区| 国产高清不卡午夜福利| 91久久精品国产一区二区三区| 久热久热在线精品观看| 蜜桃在线观看..| 精品久久久久久电影网| 国产亚洲精品久久久com| 日韩在线高清观看一区二区三区| 看非洲黑人一级黄片| 免费观看的影片在线观看| 成人综合一区亚洲| 少妇高潮的动态图| av免费在线看不卡| 精品国产国语对白av| 91精品国产九色| 人妻一区二区av| av播播在线观看一区| 精品一区在线观看国产| 春色校园在线视频观看| 天天影视国产精品| 国产成人精品在线电影| 婷婷成人精品国产| 热re99久久精品国产66热6| 精品一品国产午夜福利视频| 国产一区二区三区av在线| 热re99久久国产66热| 亚洲国产成人一精品久久久| 国产成人免费无遮挡视频| 2021少妇久久久久久久久久久| 精品一区二区三卡| 久久久亚洲精品成人影院| 国产综合精华液| 亚洲精品久久成人aⅴ小说 | 能在线免费看毛片的网站| 国产色爽女视频免费观看| 男人操女人黄网站| 国产av国产精品国产| 亚洲国产毛片av蜜桃av| 亚洲精品aⅴ在线观看| 伦理电影免费视频| 欧美少妇被猛烈插入视频| 国产69精品久久久久777片| 99re6热这里在线精品视频| 亚洲国产精品专区欧美| 国产日韩欧美视频二区| 亚洲第一区二区三区不卡| 色视频在线一区二区三区| 亚洲少妇的诱惑av| 婷婷成人精品国产| 国产日韩欧美在线精品| 最黄视频免费看| 国产精品一区www在线观看| 国产伦理片在线播放av一区| 一本—道久久a久久精品蜜桃钙片| 久久99热6这里只有精品| 国产av码专区亚洲av| 观看美女的网站| 精品人妻一区二区三区麻豆| 美女大奶头黄色视频| 久久久久人妻精品一区果冻| 欧美亚洲日本最大视频资源| 亚洲成色77777| 一级片'在线观看视频| 另类精品久久| 国产精品99久久久久久久久| 美女主播在线视频| 91精品国产国语对白视频| 久久久精品免费免费高清| 精品国产一区二区久久| 一区二区三区乱码不卡18| 尾随美女入室| 一边亲一边摸免费视频| 香蕉精品网在线| 女性生殖器流出的白浆| 国产男女内射视频| 熟妇人妻不卡中文字幕| 亚洲欧美精品自产自拍| 国产av国产精品国产| av在线app专区| 久久热精品热| 婷婷色综合大香蕉| 国产精品 国内视频| av在线播放精品| av黄色大香蕉| 日韩人妻高清精品专区| 美女xxoo啪啪120秒动态图| 久久精品国产鲁丝片午夜精品| 制服丝袜香蕉在线| 亚洲国产精品一区三区| 久热久热在线精品观看| 精品国产露脸久久av麻豆| 男女国产视频网站| 建设人人有责人人尽责人人享有的| 久久国产精品大桥未久av| 免费久久久久久久精品成人欧美视频 | 亚州av有码| 99久国产av精品国产电影| 精品人妻在线不人妻| 亚洲精品乱久久久久久| 狠狠精品人妻久久久久久综合| 汤姆久久久久久久影院中文字幕| 久久久久网色| 黑丝袜美女国产一区| 久久精品久久精品一区二区三区| 另类亚洲欧美激情| 一级毛片 在线播放| 韩国高清视频一区二区三区| 亚洲av日韩在线播放| 国语对白做爰xxxⅹ性视频网站| 性高湖久久久久久久久免费观看| 青春草亚洲视频在线观看| 日韩欧美一区视频在线观看| 日本vs欧美在线观看视频| 在线观看www视频免费| 日韩,欧美,国产一区二区三区| xxx大片免费视频| 在线观看www视频免费| 午夜福利在线观看免费完整高清在| 国产国拍精品亚洲av在线观看| 亚洲一级一片aⅴ在线观看| 亚洲精品日本国产第一区| 亚洲av在线观看美女高潮| 秋霞伦理黄片| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久久久免| 亚洲av国产av综合av卡| 免费日韩欧美在线观看| 99精国产麻豆久久婷婷| 卡戴珊不雅视频在线播放| 精品少妇黑人巨大在线播放| 少妇的逼好多水| 亚洲精品456在线播放app| 午夜福利在线观看免费完整高清在| 亚洲美女搞黄在线观看| 这个男人来自地球电影免费观看 | 国产精品久久久久成人av| 国产精品一区二区在线观看99| 亚洲精品av麻豆狂野| 人人妻人人爽人人添夜夜欢视频| 国产片内射在线| 成年女人在线观看亚洲视频| 搡女人真爽免费视频火全软件| 最后的刺客免费高清国语| 插阴视频在线观看视频| 另类亚洲欧美激情| 大片电影免费在线观看免费| 欧美成人午夜免费资源| 最新中文字幕久久久久| 女的被弄到高潮叫床怎么办| 日本午夜av视频| 婷婷色综合www| 色94色欧美一区二区| 黑人欧美特级aaaaaa片| 在线观看免费日韩欧美大片 | 亚洲av在线观看美女高潮| 久久免费观看电影| 美女cb高潮喷水在线观看| 最近中文字幕2019免费版| 日本爱情动作片www.在线观看| 尾随美女入室| 激情五月婷婷亚洲| 国产一区有黄有色的免费视频| 最黄视频免费看| 超碰97精品在线观看| 久久99蜜桃精品久久| 看非洲黑人一级黄片| 最近最新中文字幕免费大全7| 99热6这里只有精品| 母亲3免费完整高清在线观看 | 高清毛片免费看| 日本av免费视频播放| 秋霞在线观看毛片| 成人亚洲欧美一区二区av| 久久人人爽av亚洲精品天堂| 午夜福利视频在线观看免费| 亚洲国产日韩一区二区| 亚洲欧洲日产国产| 十八禁高潮呻吟视频| 国产精品一国产av| 国产片特级美女逼逼视频| 天堂中文最新版在线下载| 春色校园在线视频观看| 美女脱内裤让男人舔精品视频| 国产亚洲精品久久久com| 制服丝袜香蕉在线| 一本一本综合久久| 欧美三级亚洲精品| 下体分泌物呈黄色| 亚洲av中文av极速乱| 97精品久久久久久久久久精品| 日产精品乱码卡一卡2卡三| 晚上一个人看的免费电影| 成人手机av| 99热这里只有精品一区| 国产成人免费观看mmmm| 精品亚洲成a人片在线观看| 久久久久久久久大av| 亚洲av在线观看美女高潮| 十分钟在线观看高清视频www| 又黄又爽又刺激的免费视频.| 高清av免费在线| 亚洲精品一二三| 午夜免费鲁丝| 少妇高潮的动态图| 精品卡一卡二卡四卡免费| 亚洲成人av在线免费| 日本av免费视频播放| 国产爽快片一区二区三区| 国产精品一国产av| 国产精品久久久久久久久免| 精品久久久久久久久亚洲| 黑人巨大精品欧美一区二区蜜桃 | 91精品三级在线观看| 亚洲欧美清纯卡通| 午夜日本视频在线| 色网站视频免费| av福利片在线| 2021少妇久久久久久久久久久| 一区二区日韩欧美中文字幕 | 欧美bdsm另类| 日韩亚洲欧美综合| 最近中文字幕2019免费版| 国产视频首页在线观看| 十八禁高潮呻吟视频| 中文欧美无线码| 久久午夜福利片| 国产午夜精品一二区理论片| 这个男人来自地球电影免费观看 | 成人毛片a级毛片在线播放| 大陆偷拍与自拍| 性高湖久久久久久久久免费观看| 亚洲欧洲国产日韩| 久久久久网色| 免费看光身美女| 自拍欧美九色日韩亚洲蝌蚪91| 在线观看www视频免费| 岛国毛片在线播放| 精品久久久久久久久亚洲| 日韩免费高清中文字幕av| 免费黄色在线免费观看| 亚洲欧美日韩卡通动漫| 国产深夜福利视频在线观看| 久久精品国产自在天天线| 18禁观看日本| 2022亚洲国产成人精品| 国产视频内射| 成人毛片a级毛片在线播放| 51国产日韩欧美| 国产永久视频网站| 国产高清国产精品国产三级| 香蕉精品网在线| 国产精品一区二区在线观看99| 日韩一本色道免费dvd| 99九九在线精品视频| 久久久欧美国产精品| 亚洲经典国产精华液单| 日韩伦理黄色片| 欧美日韩亚洲高清精品| 午夜福利网站1000一区二区三区| 91在线精品国自产拍蜜月| 精品一品国产午夜福利视频| 欧美精品国产亚洲| 亚洲精品456在线播放app| 18在线观看网站| 在现免费观看毛片| 国产精品一区二区在线观看99| 成人18禁高潮啪啪吃奶动态图 | 91精品国产九色| 天堂中文最新版在线下载| 大片电影免费在线观看免费| 在线亚洲精品国产二区图片欧美 | 日韩免费高清中文字幕av| 伦理电影大哥的女人| 18禁在线播放成人免费| 亚洲精品乱码久久久v下载方式| 成人手机av| 少妇 在线观看| 久久99一区二区三区| 国产av国产精品国产| 啦啦啦在线观看免费高清www| 91精品国产九色| 99热国产这里只有精品6| 最后的刺客免费高清国语| 色网站视频免费| 最近中文字幕高清免费大全6| 日韩一区二区视频免费看|