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

    廣西欽州灣百年來紅樹林演變的有機碳同位素和孢粉示蹤及其影響因素

    2015-06-01 12:30:12夏鵬孟憲偉平愛平李珍
    海洋學(xué)報 2015年3期
    關(guān)鍵詞:孢粉紅樹欽州

    夏鵬,孟憲偉*,平愛平,李珍

    (1.國家海洋局第一海洋研究所,山東青島 266061;2.華東師范大學(xué)河口海岸學(xué)國家重點實驗室,上海 200062)

    廣西欽州灣百年來紅樹林演變的有機碳同位素和孢粉示蹤及其影響因素

    夏鵬1,孟憲偉1*,平愛平1,李珍2

    (1.國家海洋局第一海洋研究所,山東青島 266061;2.華東師范大學(xué)河口海岸學(xué)國家重點實驗室,上海 200062)

    氣候變化和人類活動制約下的紅樹林演變是一種長時間尺度效應(yīng),而沉積物則是記錄這種響應(yīng)的最佳檔案。采用古生態(tài)學(xué)研究思路,選擇有效的示蹤參數(shù)是解讀這一檔案的有效途徑。本文通過廣西欽州灣紅樹林區(qū)1根柱狀沉積物中有機碳同位素(δ13C)、C/N分析和孢粉鑒定,以沉積物中的紅樹林源有機碳貢獻(xiàn)和紅樹植物孢粉組合為示蹤參數(shù),在210Pb年齡框架構(gòu)建的基礎(chǔ)上,通過對比研究紅樹林源有機碳貢獻(xiàn)和紅樹孢粉組合特征,揭示百年來欽州灣紅樹林興衰和群落演替規(guī)律;進而重塑其紅樹林演變歷史:興盛期(1864-1918年)、衰退期(1918-1968年)和低谷期(1968-2007年),這與土地替代、遙感分析和現(xiàn)場勘測數(shù)據(jù)較為吻合。結(jié)合氣候變化和人類活動資料,發(fā)現(xiàn)欽州灣紅樹林的近期衰退主要源于人類活動的影響,尤其是蝦塘圍墾;而與氣候變化關(guān)系不大。

    紅樹林;有機碳同位素;孢粉;氣候變化;人類活動

    1 引言

    高生產(chǎn)力的紅樹林生態(tài)系統(tǒng)[平均生產(chǎn)力(以碳計)2 500 mg/(m2·d)]覆蓋了60%~75%的熱帶海岸線[1—2],占陸源碳輸出入海量的11%和海洋埋藏總有機碳的15%[3—4],在近岸環(huán)境生物地球化學(xué)循環(huán)過程中起到了重要作用。此外,紅樹林可以攔截進入沿海的陸源營養(yǎng)物質(zhì)、污染物和沉積物,減輕高負(fù)荷的人為輸入[5]。紅樹林生態(tài)系統(tǒng)對氣候變化、人類活動非常敏感[6—7],為了解當(dāng)前乃至將來紅樹林生態(tài)系統(tǒng)對外界環(huán)境變化的響應(yīng),掌握歷史時期紅樹林生態(tài)系統(tǒng)的興衰演變規(guī)律是至關(guān)重要的[8]。高生產(chǎn)力和低泥沙呼吸速率相結(jié)合,使紅樹林沉積物有著長期封存有機碳的潛力[9]。埋藏的紅樹林有機遺骸化石,如花、果、葉和孢粉等,是用于紅樹林古生態(tài)重建的重要化石證據(jù)[10—11]。潮濕的紅樹林區(qū)沉積物易于形成厭氧環(huán)境,利于這些化石記錄的長期保存[11]。

    δ13C和C/N比值再輔以其他示蹤劑(如:類脂化合物、氨基酸、木質(zhì)素和長鏈烷烴等)已被廣泛應(yīng)用于近岸和河口沉積物中有機質(zhì)來源的辨識[12—13],也被用于示蹤食物鏈中紅樹林源有機質(zhì)的轉(zhuǎn)移路徑[7,10]。然而,有機質(zhì)來源的復(fù)雜性和早期成巖作用,致使不同端元的有機質(zhì)同位素信號發(fā)生不同程度的改變[12—14]。沉積物中有機質(zhì)來源的重塑僅反映被保存下來的殘存有機組分,而非當(dāng)初有機質(zhì)的實際埋藏通量;例如,紅樹林碳比浮游藻類碳頑固,將會被優(yōu)先保存[15—16]。而紅樹孢粉是示蹤紅樹林沼澤沉積環(huán)境的獨特指標(biāo),可以識別不同紅樹種群的群落演替特征[11,17]。因此,孢粉圖譜有助于彌補有機碳同位素的某些限制[17—18],增加紅樹林古生態(tài)恢復(fù)的可靠性。但是,在紅樹林古生態(tài)恢復(fù)的國內(nèi)外文獻(xiàn)中,有機碳同位素和孢粉數(shù)據(jù)相互印證的研究卻鮮有報道。

    本文通過廣西欽州灣紅樹林濕地的1根柱狀沉積物中δ13C、TOC、TN分析以及孢粉鑒定,以沉積物中紅樹林源有機碳貢獻(xiàn)和紅樹孢粉組合為示蹤參數(shù),在210Pb年齡框架構(gòu)建的基礎(chǔ)上,通過紅樹林源有機碳貢獻(xiàn)和紅樹孢粉組合特征的對比研究,揭示百年來欽州灣紅樹林興衰和群落演替規(guī)律,進而重塑其紅樹林演變歷史;結(jié)合欽州灣氣候變化和人類活動資料,建立百年來廣西欽州灣紅樹林演變對氣候變化和人類活動的響應(yīng)關(guān)系。

    2 樣品采集與分析

    2.1 樣品采集

    2007年4月27日,趁低潮時在欽州灣茅嶺江口南側(cè)的紅樹林沼澤外緣潮灘插管獲取了3根80~90 cm長的柱狀巖心,從中選取1根90 cm長且沉積連續(xù)的柱狀巖心Q37作為研究對象(圖1),其沉積類型以粉砂質(zhì)砂(TS)為主。利用美國Trimble公司的DSM 212H雙信標(biāo)DGPS進行定位,平面定位精度優(yōu)于1 m。在樣品采集和運輸?shù)倪^程中盡量保持沉積物的原始狀態(tài)。實驗室內(nèi)按2 cm間隔進行切割,子樣品用于210Pb測年、有機碳同位素和孢粉等分析。同步采集Q37柱周邊不同種類的紅樹植物葉片,用于紅樹林端元有機碳同位素的測定。

    圖1 區(qū)域背景和Q37柱的站位分布圖Fig.1 Location of Maowei Sea and sampling site

    2.2 分析方法

    210Pb測年:對入管蠟封并放置3周的5 g干樣進行210Pb、226Ra和137Cs的放射性比活度測試,采用EG&G Ortec公司生產(chǎn)的高純鍺低本底γ能譜儀完成[19];但由于沉積物中137Cs含量過低,均未超過檢出限。分析工作在華東師大河口海岸國家重點實驗室完成。李貞等[20]基于CIC模型計算了Q37柱的沉積速率為0.63 cm/a,時間跨度為143 a;由于巖芯頂部的年代為樣品采集時間2007年,底層年代為1864年。

    孢粉分析:每個樣品取樣5 g,經(jīng)鹽酸處理去碳酸鈣后,用重液分離法使孢粉富集、制片。Leitz顯微鏡下放大250倍觀察統(tǒng)計2個20 mm×20 mm蓋玻片下的孢粉數(shù)量,作為該樣品的化石豐度。對豐度在50粒/片以上的進行統(tǒng)計,計算屬種個數(shù)和百分含量。參考圖譜分別為《中國熱帶亞熱帶被子植物花粉形態(tài)》和《中國植物花粉形態(tài)》[21]。分析工作在華東師大河口海岸國家重點實驗室完成。運用程序CONISS[22]進行聚類分析,并參考孢粉成分的相對變化來確定孢粉的組合帶[20]。

    TOC和TN分析:TOC(有機碳)和TN(總氮)分別利用重鉻酸鉀—硫酸氧化—硫酸亞鐵法滴定法和重鉻酸鉀—硫酸消化—凱氏定氮法測試[23]。分析工作在國家海洋局第一海洋研究所海洋沉積與環(huán)境地質(zhì)實驗室完成。有機碳埋藏通量BForg的計算公式為:

    式中,BForg為有機碳埋藏通量[單位:mg/(cm·a)];Ci為有機碳含量(單位:mg/g);Si為沉積物累積速率(單位:cm/a);ρi為沉積物的干密度(單位:g/cm);i為沉積物某一層的深度(單位:cm)。

    有機碳同位素分析:樣品用清水洗凈、晾干,先用1.5 mol/L鹽酸處理以去除無機碳酸鹽,置于烘箱中烘干(65℃),研磨過40目篩。研磨后的樣品在Delta-PlusXP質(zhì)譜儀(Thermo Scientific,Bremen,Germany)上進行碳同位素值測定,帶元素分析儀(Flash EA1112,CE Instruments,Wigan,UK),燃燒溫度為1 020℃。分析工作在中國農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院穩(wěn)定同位素實驗室完成,碳同位素值用國際PDB標(biāo)準(zhǔn)表示,平行樣分析誤差小于0.2‰。

    3 結(jié)果

    3.1 孢粉圖譜

    Q37柱中孢粉濃度變化范圍介于4 536~17 076粒/g,表層孢粉濃度總體略低于底層。鑒定植物種類主要包括:(1)陸生木本植物(Pinus、Alnus、Morus、Cedrus、Meliaceae、Syzygium、Eucalyptus、Anacardiaceae、Engelhardia、Sapindus、Eurya、Cratoxylum、Castanopsis、Tetrameles、Mallotus、Rutaceae、Quercus、Liquidambar、Lonicera、Quercus、Hamamelidaceae);(2)半紅樹植物(Arenga、Malvaceae、Casuarina、Trachycarpus、Quisqualis);(3)紅樹植物(Bruguiera、Aegiceras、Rhizophora);(4)中生草本植物(Daphne、Loranthus、Poaceae、Piper、Euphorbia、Orchidaceae);(5)耐旱草本植物(Artemisia、Compositae、Chenopodium);(6)濕生/水生草本植物(Lilium、Cyperaceae、Acorus);(7)蕨類孢子(Gleicheniaceae、Cyatheaceae、Polypodiaceae、Pteridiaceae)。其中,草本植物類型的孢粉含量最高,介于36.2%~54.5%之間;其次是木本植物,介于21.5%~45.5%;蕨類孢子含量最少。孢粉組合的變化主要表現(xiàn)在針葉、熱帶雨林以及紅樹植物類型的變化,結(jié)合聚類分析結(jié)果,可劃分為3個孢粉組合帶(圖2)。

    圖2 Q37柱孢粉組合特征及其CONISS分帶Fig.2 Pollen diagram and its three zones based on CONISS cluster analysis from Core Q37

    孢粉Ⅰ帶(88~60 cm,1864-1911年):該帶以孢粉組合Gramineae-Cyatheaceae-Rhizophora-Eurya-Piper為主要特征,主要孢粉類型為木本植物,介于25.7%~53.6%;其中,紅樹花粉(10%~22%)和熱帶喬木花粉占主導(dǎo)位置,紅樹植物主要為R.stylosa(紅海欖,3.6%~12.7%)和A.corniculatum(桐花樹,4.5%~8.4%),B.gymnorrhiza(木欖)花粉主要出現(xiàn)在底部(圖3)。熱帶喬木類花粉以Eurya(最大13.9%)、Mallotus(0~5.7%)和Piper(2.4%~9.5%)為主,Pinus的含量極低。草本植物類型中以禾本科為主,其花粉含量與木本相近,但是二者的變化趨勢相反。耐鹽的濱藜及蒿等也有少量呈現(xiàn),表明沉積物固結(jié)程度較高。孢粉總濃度介于74 000~102 000粒/g,但蕨類孢子含量較低,反映了較弱的水動力環(huán)境。參照現(xiàn)代紅樹孢粉譜特征,表明該期區(qū)域植被以熱帶雨林為主;采樣點應(yīng)處于以紅海欖為主的群落邊緣,外緣有桐花樹分布。

    圖3 紅樹植物孢粉圖式Fig.3 Vertical distribution of mangrove species

    孢粉Ⅱ帶(60~30 cm,1911-1960年):該帶以木本植物花粉含量降低、紅樹植物花粉銳減為主要特征,且草本植物花粉保持穩(wěn)定高值;其孢粉組合為Gramineae-Cyatheaceae-Cyperaceae-Rhizophora-Polypodiaceae。孢粉Ⅱ帶濃度比Ⅰ帶略有增加,保持在較高水平;蕨類孢子也略有增加(尤其是Polypodiaceae和Cyatheaceae),反映了較弱的水動力環(huán)境。熱帶植物花粉類型多樣性變化不大,但含量明顯降低;針葉植物成分增加,其中Pinus占6.0%~7.0%。紅樹植物與半紅樹植物成分均相應(yīng)減少,紅海欖含量不足10%。而以Gramineae為主的中生草本植物含量增大,濕生/淡水生植物Cyperaceae增加至5.4%~8.3%。

    孢粉Ⅲ帶(30~0 cm,1960-2007年):孢粉組合特征為Gramineae-Pinus-Cyatheaceae-Cyperaceae-Polypodiaceae。孢粉總濃度銳減到24 000~52 000粒/g,而Pinus含量遞增至18.4%~20.8%,蕨類孢子含量也達(dá)到最大,表明水動力強度明顯增強。紅樹植物花粉降到最低值,伴隨著桐花樹增加和紅海欖減少;但二者含量均不足5%,表明均由異地攜入,非原生紅樹植物產(chǎn)生。與人類活動密切相關(guān)的Casuarina和Morus等人工防護林類型明顯增加,而半紅樹種類和熱帶植物花粉類型均相應(yīng)減少。草本植物中仍以Gramineae為主要成分,但含量略有降低,濕生/淡水生植物中的Cyperaceae、Lilium和Acorus明顯增加。紅樹植物與半紅樹植物成分均減少,桐花樹含量不足10%;而中生草本植物孢粉含量增大,增加至5.4%~8.3%。該孢粉組合反映了區(qū)域熱帶雨林減少和人類種植及次生林的增加,該區(qū)以紅樹林消退、水動力環(huán)境增強為主要變化特征。

    綜上所述,Q37柱總體反映了3段不同的沉積環(huán)境:(1)1864-1911年間,該區(qū)為紅樹林灘,群落優(yōu)勢種為紅海欖和桐花樹,且處于紅樹林群落內(nèi)緣;(2)1911-1960年間,紅樹林灘呈現(xiàn)衰退趨勢,伴隨著草本植物的大量增加及區(qū)域植物類型的多樣性增強,均反映出受河流影響增強的趨勢[24];(3)1960-2007年間,孢粉總濃度迅速減少,而蕨類孢子和Pinus明顯增加,反映出水動力環(huán)境有所增強[24—25],應(yīng)與灘面上覆紅樹植物的減少有關(guān)。

    3.2 有機碳同位素和C/N比

    3.2.1 端元值的確定

    欽州灣潮間帶沉積物中埋藏的有機碳具有多源性,如陸源、紅樹林源和海洋生物源;因此選擇合適的參數(shù)(δ13C和C/N)定量分離出沉積物中紅樹林源有機碳的貢獻(xiàn)成為利用紅樹林源有機碳貢獻(xiàn)示蹤紅樹林演變的關(guān)鍵。而δ15N同位素不適于示蹤有機質(zhì)來源的主要原因是15N易于沿食物鏈富集[26-27]以及富15N農(nóng)業(yè)和城鎮(zhèn)廢水的排放[28-29]。3個不同端元有機質(zhì)之間的δ13C和C/N比值存在著明顯差異(見圖4a),這是定量分離沉積物中有機質(zhì)來源的先決條件之一。對不同紅樹植物葉片(即桐花樹、秋茄、木欖和紅海欖)的δ13C或C/N進行兩兩T檢驗,p均小于0.05,表明不同紅樹植物的δ13C或C/N不存在顯著差異;可以組合在一起計算紅樹林源有機碳的端元值。

    圖4 三端元有機碳的δ13C和C/N散點分布(a)和Q37柱中有機碳的三端元解析圖(b)Fig.4 δ13C and C/N for three end member samples(a),andδ13C and C/N in ternary mixing diagrams for mangrove,terrestrial OM and phytoplankton(b)

    紅樹植物葉片的δ13C(-28.75‰±0.89‰)和C/N(40.15±11.18)值與肯尼亞[31]和墨西哥灣[30]的紅樹植物葉片δ13C和C/N較為一致。然而,相比紅海欖(δ13C:-29.14‰±0.46‰,C/N:40.15±0.15)、桐花樹(δ13C:-29.13‰±0.90‰,C/N∶40.13± 7.44)和秋茄(δ13C:-28.35‰±0.85‰,C/N:37.38 ±13.11),木欖有著較高的δ13C和C/N值(δ13C:-28.39‰±0.88‰,C/N:45.75±18.12;圖4a)。陸源有機碳的平均δ13C(-24.12‰±0.60‰)和C/N比(12.76±2.11)與珠江上游表層沉積物(δ13C:-24.9‰±1.0‰,C/N:13.3±1.4)[32]的較為接近,源于河床沉積物相似的物源背景。本次選取南海北部陸坡浮游植物的平均δ13C(-16.05‰±0.82‰)和C/N比(6.46±0.10)作為對海洋生物源有機碳的端元值[32—33]。

    3.2.2 有機碳來源解析

    Q37柱中δ13C、TOC、TN、C/N摩爾比和BForg的垂向分布特征如圖5a~d所示。TOC、TN、C/N和BForg的變化趨勢較為一致,兩兩之間保持高度正相關(guān)關(guān)系;其中,TOC-TN、TOC-C/N和TOC-BForg的相關(guān)系數(shù)分別為0.85(P<0.01)、0.83(P<0.01)和0.99(P<0.01)。而δ13C與TOC、TN、C/N、BForg的變化趨勢恰好相反,即較重的δ13C值對應(yīng)較低的TOC含量;其中,δ13C-TOC和δ13C-C/N的相關(guān)系數(shù)分別為-0.84(P<0.01)和-0.76(P<0.01)。

    基于Q37柱中δ13C-C/N的三端元散點分布特征(圖4b),78%的樣品點落于10%耐受水平下的三角區(qū)域內(nèi),表明可用三端元計算公式[34]來解釋。其余落于三角區(qū)域外的10個樣品點多位于深層沉積物中,其埋藏有機質(zhì)可能受早期成巖作用的影響發(fā)生輕微改變[12,14,35];可用鄰近兩端元來近似解釋[34]。結(jié)合210Pb年代框架,可將紅樹林源有機碳貢獻(xiàn)大致劃分為3個階段:

    階段Ⅰ(89~56 cm,1864-1918年):沉積物具有較高的TOC和TN含量,均值分別為1.51%和0.06%,且其含量自下往上遞減趨勢明顯;而δ13C的表現(xiàn)相對穩(wěn)定,介于-27.89‰~-25.59‰(均值-27.34‰)。該段有機碳的埋藏通量較高,平均有機碳埋藏通量為12.15 mg/(cm2·a)[6.76~22.14 mg/(cm2·a)]。紅樹林源有機碳貢獻(xiàn)處于高值區(qū),平均貢獻(xiàn)為81.45%,其次為陸源;即高含量的有機碳主要源于紅樹林的貢獻(xiàn)。

    階段Ⅱ(56~24 cm,1918-1965年):TOC和TN含量明顯遞減,均值分別為0.66%和0.03%,總體呈現(xiàn)“先減后增”的輕微變化趨勢;而δ13C遞增趨勢明顯,介于-26.72‰~-23.89‰(均值-25.05‰)。該段有機碳的平均埋藏通量為5.88 mg/(cm2·a)[2.62~10.15 mg/(cm2·a)],埋藏通量較之階段Ⅰ減少了50%以上。該段陸源、紅樹林源和海源有機碳的平均貢獻(xiàn)分別為50.00%、30.92%和19.08%;由于陸源有機碳的入海量年際變化不大,有機碳埋藏通量的大幅減少應(yīng)與紅樹林源有機碳的減少有關(guān)。

    階段Ⅲ(24~0 cm,1965-2007年):TOC和TN含量保持在較低水平上,均值分別為0.42%和0.03%;而δ13C介于-25.42‰~-23.75‰(均值-24.33‰)。該段有機碳的平均埋藏通量為3.95 mg/(cm2·a)[2.62~7.92mg/(cm2·a)],埋藏通量較之階段Ⅱ減少了30%以上。該段陸源、紅樹林源和海源有機碳的平均貢獻(xiàn)量分別為63.93%、12.85%和23.22%;有機碳埋藏通量的持續(xù)減少主要源于紅樹林源有機碳的減少。

    圖5 Q37柱中δ13C(a),TOC(實心)和TN(空心)(b),C/N摩爾比(c),BForg(d),有機碳貢獻(xiàn)(e)的垂向分布特征Fig.5 Vertical distributions ofδ13C(a),TOC(filled circles)and TN(empty circles)(b),C/N(c),BForg(d)and% OM source and depth(e)

    4 討論

    4.1 基于紅樹孢粉和有機碳同位素綜合示蹤紅樹林演變

    由于有機碳同位素復(fù)雜的物源特征及其可能受早期成巖作用或有機質(zhì)降解作用的影響,在單獨示蹤紅樹林演變時存在弊端;需要與更直觀指示紅樹林生長狀態(tài)的紅樹花粉含量進行相互印證[11]。為了研究紅樹林源有機碳貢獻(xiàn)(見圖6a)與紅樹植物花粉變化趨勢是否一致(見圖6b),基于最小二乘法的線性回歸被廣泛應(yīng)用于調(diào)查兩個變量間的關(guān)系(見圖6c)。研究表明:二者之間呈顯著正相關(guān)(r=0.80,P<0.01),且散點均落于95%預(yù)測區(qū)間內(nèi);表明有機碳同位素在示蹤歷史時期紅樹林演變階段時取得與紅樹花粉一致的結(jié)論。因此,紅樹林生長演化可以劃分為3個階段:興盛期(1864-1918年)、衰退期(1918-1968年)和低谷期(1968-2007年)。這與土地替代、遙感分析和現(xiàn)場實測同期數(shù)據(jù)較為吻合:廣西海岸曾經(jīng)有23 904 hm2紅樹林,是現(xiàn)存面積的近3倍[36];但到解放后的1955年減少了61%,僅剩下9 351.2 hm2[37];自1955年以后紅樹林面積急速下降,到1988年廣西紅樹林面積下降到1955年的50%,即4 671.4 hm2[37]。

    4.2 與氣候變化和人類活動的關(guān)系

    近年來,紅樹林生態(tài)系統(tǒng)受氣候變化和人類活動的影響嚴(yán)重[6];為了查清紅樹林衰退(見圖7a~b)與二者之間的相互關(guān)系,本文選用了氣溫、降雨量和日照時數(shù)(見圖7e~g)作為氣候因子以及圍填海面積和海堤長度(見圖7c~d)作為人類活動因子進行對比研究。研究發(fā)現(xiàn):近50年來的氣溫具有十年際變化趨勢,一個周期約15 a,與太平洋十年濤動(PDO)[38]類似,且呈逐年遞增的趨勢;日照時數(shù)卻呈逐年減少的趨勢;而降雨量年際增減趨勢不明顯。顯然,氣溫和日照時數(shù)的輕微變化不是造成紅樹林大幅衰退的主要原因。而近期記錄的欽州灣圍填海面積和海堤長度總體呈快速遞增趨勢,尤其是蝦塘面積;這與范航清[36]記錄的紅樹林區(qū)蝦塘圍墾是導(dǎo)致近期紅樹林生態(tài)系統(tǒng)快速衰退的主要原因相吻合。

    圖6 紅樹林源有機碳貢獻(xiàn)(a)和紅樹花粉的年際變化(b)及其散點關(guān)系圖(c)Fig.6 Mangrove-derived OM contribution(a),mangrove pollen content(b),and its scatter plots(c)

    圖7 紅樹林演化(a、b)與人類活動(c、d)和氣候變化(e~g)的關(guān)系Fig.7 Mangrove development(a,b)and its relationship with anthropogenic activity(c,d)and climate change(e,g)

    5 結(jié)論

    由于孢粉數(shù)據(jù)難于準(zhǔn)確定量和有機碳同位素復(fù)雜的物源特征及其受早期成巖作用的影響,它們在單獨示蹤近期紅樹林演化狀況時互有利弊。因此,本文以沉積物中紅樹林源有機碳貢獻(xiàn)和紅樹植物孢粉組合為示蹤參數(shù),在210Pb年代框架和有機碳端元值構(gòu)建的基礎(chǔ)上,通過紅樹林源有機碳貢獻(xiàn)和紅樹植物孢粉組合變化的對比研究,相互印證并揭示了近150年來欽州灣紅樹林興衰和群落演替規(guī)律。進而將欽州灣紅樹林的演化劃分為3個階段:興盛期(1864-1918年)、衰退期(1918-1968年)和低谷期(1968-2007年),這與土地替代、遙感分析和現(xiàn)場實測數(shù)據(jù)相吻合。并結(jié)合氣候變化和人類活動資料,發(fā)現(xiàn)近期欽州灣紅樹林衰退主要源于人類活動的影響,尤其是蝦塘圍墾;而與氣候變化關(guān)系不大。

    致謝:中國農(nóng)業(yè)大學(xué)王國安教授、華東師范大學(xué)李貞、國家海洋局第一海洋研究所王湘芹、楊剛、朱志偉等參加了現(xiàn)場觀測、樣品采集和分析工作,在此一并感謝!

    [1] Clough B.Mangrove forest productivity and biomass accumulation in Hinchinbrook Channel,Australia[J].Mangroves and Salt Marshes,1998,2:191-198.

    [2] Jennerjahn T C,Ittekkot V.Relevance of mangroves for the production and deposition of organic matter along tropical continental margins[J].Naturwissenschaften,2002,89:23-30.

    [3] Alongi D M.The dynamics of benthic nutrient pools and fluxes in tropical mangrove forests[J].Journal of Marine Research,1996,54:123-148.

    [4] Lee S Y.Mangrove outwelling:a review[J].Hydrobiologia,1995,295:203-212.

    [5] Valiela I,Cole M L.Comparative evidence that salt marshes and mangroves may protect seagrass meadows from land-derived nitrogen loads[J].E-cosystems,2002,5:92-102.

    [6] Alongi D M.Mangrove forests:resilience,protection from tsunamis,and responses to global climate change[J].Estuarine,Coastal and Shelf Science,2008,76:1-13.

    [7] Chong V C,Low C B,Ichikawa T.Contribution of mangrove detritus to juvenile prawn nutrition:a dual stable isotope study in a Malaysian mangrove forest[J].Marine Biology,2001,138:77-86.

    [8] Valiela I,Bowen J L,York J K.Mangrove forests:one of the world's threatened major tropical environments[J].Bioscience,2001,51:807-815.

    [9] Kristensen E,Bouillon S,Dittmar T,et al.Organic carbon dynamics in mangrove ecosystems:A review[J].Aquatic Botany,2008,89:201-219.

    [10] Bouillon S,Connolly R M,Lee S Y.Organic matter exchange and cycling in mangrove ecosystems:recent insights from stable isotope studies[J].Journal of Sea Research,2008,59:44-58.

    [11] Ellison JC.Long-term retrospection on mangrove development using sediment cores and pollen analysis:A review[J].Aquatic Botany,2008,89:93-104.

    [12] Graham M C,Eaves M A,F(xiàn)armer JG,et al.A study of carbon and nitrogen stable isotope and elemental ratios as potential indicators of source and fate of organic matter in sediments of the Forth Estuary,Scotland[J].Estuarine,Coastal and Shelf Science,2001,52:375-380.

    [13] Yamamuro M.Chemical tracers of sediment organic matter origins in two coastal lagoons[J].Journal of Marine Systems,2000,26:127-134.

    [14] Fourqurean J W,Schrlau J E.Changes in nutrient content and stable isotope ratios of C and N during decomposition of seagrasses and mangrove leaves along a nutrient availability gradient in Florida Bay,USA[J].Chemistry and Ecology,2003,19,373-390.

    [15] De Boer W.Biomass dynamics of seagrasses and the role of mangrove and seagrass vegetation as different nutrient sources for an intertidal ecosystem[J].Aquatic Botany,2000,66:225-239.

    [16] Holmer M,Olsen A B.Role of decomposition of mangrove and seagrass detritus in sediment carbon and nitrogen cycling in a tropical mangrove forest[J].Marine Ecology Progress Series,2002,230:87-101.

    [17] Wooller M J,Swain D L,F(xiàn)icken K J,et al.Late Quaternary vegetation changes around Lake Rutundu,Mount Kenya,East Africa:evidence from grass cuticles,pollen and stable carbon isotopes[J].Journal of Quaternary Science,2003,18(1):3-15.

    [18] Ficken K J,Wooller MJ,Swain D L,et al.Reconstruction of a sub-alpine grassdominated ecosystem,Lake,Rutundu,Mount Kenya:a novel multiproxy approach[J].Palaeogeography,Palaeoclimatology,Palaeoecology,2002,177:137-150.

    [19] San Miguel E G,Bolívar J P,García-tenorio R.Vertical distribution of Th-isotope ratios,210Pb,226Ra and137Cs in sediment cores from an estuary affected by anthropogenic releases[J].Science of the Total Environment,2004,318:143-157.

    [20] 李貞,李珍,張衛(wèi)國,等.廣西欽州灣海岸帶孢粉組合和沉積環(huán)境演變[J].第四紀(jì)研究,2010,30(3):598-608.

    Li Zhen,Li Zhen,Zhang Weiguo,et al.Pollen assemblage and sedimentary environment evolution in the coastal zone in Qinzhou Bay,Guangxi province[J].Quaternary Sciences,2010,30:598-608.

    [21] 王伏雄,錢南芬,張玉龍,等.中國植物花粉形態(tài)[M].北京:科學(xué)出版社,1995.

    Wang Fuxiong,Qian Nanfen,Zhang Yulong,et al.Pollen morphology of Chinese plant[M].Beijing:Science Press,1995.

    [22] Grimm E C.Coniss:A Fortran 77 program forstratigra phically constrained cluster analysis by method of incremental sum of squares[J].Computers&Geosciences,1987,13:13-35.

    [23] 劉建軍,吳敬祿.太湖大浦湖區(qū)近百年來湖泊記錄的環(huán)境信息[J].古地理學(xué)報,2006,8(4):559-564.

    Liu Jianjun,Wu Jinglu.Environmental information of recent 100 years recorded sediments of Dapu area in Taihu Lake[J].Journal of palaeogeography,2006,8:559-564.

    [24] Van Der Kaars S.Pollen distribution in marine sediments from the south-eastern Indonesian waters[J].Palaeogeography,Palaeoclimatology,Palaeoecology,2001,171:341-361.

    [25] Hooghiemstra H,Lézine A M,Leroy S A G,et al.Late Quaternary palynology in marine sediments:A synthesis of the understanding of pollen distribution patterns in the NW African setting[J].Quaternary International,2006,148:29-44.

    [26] Cole M L,Valiela I,Kroeger K D,et al.Assessment of aδ15N isotopic method to indicate anthropogenic eutrophication in aquatic ecosystems[J].Journal of Environmental Quality,2004,33:124-132.

    [27] Kwak T J,Zedler JB.Food web analysis of southern California coastal wetlands using multiple stableisotopes[J].Oecologia,1997,110(2):262-277.

    [28] Mutchler T,Dunton K H,Townsend-small A,et al.Isotopic and elemental indicators of nutrient sources and status of coastal habitats in the Caribbean Sea,Yucatan Peninsula,Mexico[J].Estuarine Coastal and Shelf Sciences,2007,74:449-457.

    [29] Savage C.Tracing the influence of sewage nitrogen in a coastal ecosystem using stable nitrogen isotopes[J].Ambio,2005,34(2):145-150.

    [30] Gonneeaa M E,Paytana A,Herrera-silveira J A.Tracing organic matter sources and carbon burial in mangrove sediments over the past 160 years[J].Estuarine,Coastal and Shelf Science,2004,61:211-227.

    [31] Rao R G,Woitchik A F,Goeyens L,et al.Carbon,nitrogen contents and stable carbon isotope abundancein mangrove leaves from an East-African coastal lagoon(Keyna)[J].Aquatic Botany,1994,7:175-183.

    [32] He B,Dai M,Huang W,et al.Sources and accumulation of organic carbon in the Pearl River Estuary surface sediment as indicated by elemental,stable carbon isotopic,and carbohydrate compositions[J].Biogeosciences,2010,7:3343-3362.

    [33] 侯秀瓊.流沙灣海洋生物穩(wěn)定碳氮同位素應(yīng)用的初步研究[D].湛江:廣東海洋大學(xué),2009:16-17.

    Hou Xiuqiong.Preliminary study on marine species in Liusha Bay by using stable carbon and nitrogen isotopes[D].Zhanjing:Guangdong Ocean University,2009,16-17.

    [34] Dittmar T,Lara R J,Kattner G.River or mangrove?Tracing major organic matter sources in tropical Brazilian coastal waters[J].Marine Chemistry,2001,73:253-271.

    [35] Yamamuro M,Kayanne H.Rapid direct determination of organic carbon and nitrogen in carbonate-bearing sediments with a Yanaco MT-5 CHN analyzer[J].Limnology and Oceanography,1995,40(5):1001-1005.

    [36] 范航清.廣西海岸紅樹林現(xiàn)狀及人為干擾[M]//范航清,梁士楚.中國紅樹林研究與管理.北京:科學(xué)出版社,1995.

    Fan Hangqing.Mangroves status and factitious disturbance in the coast of Guangxi[M]//Mangroves Study and Management in China.Beijing:Science Press,1995,189-202.

    [37] 黃鵠,戴志軍,胡自寧,等.廣西海岸環(huán)境脆弱性研究[M].北京:海洋出版社,2005.

    Huang Hao,Dai Zhijun,Hu Zining,et al.Study on the coastal environmental vulnerability in Guangxi[M].Beijing:China Ocean Press,2005.

    [38] Mantua N J,Hare S R,Zhang Y,et al.A Pacific decadal climate oscillation with impacts on salmon[J].Bulletin of the American Meteorological Society,1997,78:1069-1079.

    Historical retrospection on mangrove development using stable carbon isotopes and pollen analysis,and its response to climate change and human activity

    Xia Peng1,Meng Xianwei1,F(xiàn)eng Aiping1,Li Zhen2

    (1.First Institute of Oceanography,State Oceanic Administration,Qingdao 266061,China;2.State Key Laboratory of Estuarine and Coastal Research,East China Normal University,Shanghai 200062,China)

    Mangrove evolution is the long time-scale process under the condition of climate change and human activity,and is clearly recorded by the surrounding sediments.It is a new palaeoecological approach to trace mangrove evolution using organic carbon isotopes and mangrove pollen from sedimentary records.In the study,a sediment core was collected from a mangrove swamp of the Qinzhou Bay,Guangxi.Based on the three end-number model,δ13C and C/N analysis were used to calculate the contribution of mangrove-derived organic carbon.Compared with mangrove pollen assemblage and210Pb dating,it revealed the temporal evolution of mangrove development and its community succession in the recent 150 years.The mangrove evolution could be divided into three major stages:flourish(1864-1918 A.D.),degradation(1918-1968 A.D.)and low-level period(1968-2007 A.D.),which was consistent with the previous reports.An significantly degradation in mangrove ecosystem recently corresponded to the increase of reclamation area(especially for the shrimp-pond reclamation)and artificial seawalls,rather than the climate change as recorded in the region.

    mangrove;carbon;isotope;pollen;climate change

    P736.4

    A

    0253-4193(2015)03-0077-09

    夏鵬,孟憲偉,豐愛平,等.廣西欽州灣百年來紅樹林演變的有機碳同位素和孢粉示蹤及其影響因素[J].海洋學(xué)報,2015,37(3):77—85,

    10.3969/j.issn.0253-4193.2015.03.008

    Xia Peng,Meng Xianwei,F(xiàn)eng Aiping,et al.Historical retrospection on mangrove development using stable carbon isotopes and pollen analysis,and its response to climate change and human activity[J].Haiyang Xuebao,2015,37(3):77—85,doi:10.3969/j.issn.0253-4193.2015.03.008

    2014-02-26;

    2014-06-29。

    國家自然科學(xué)基金(41206057,41376075);973計劃(2010CB951203);中央級公益性科研院所基本科研業(yè)務(wù)費專項資金(2012G18);國家海洋局青年基金(2013326)。

    夏鵬(1982—),男,山東省濰坊市人,副研,博士,從事海洋地球化學(xué)研究。E-mail:pengxia@fio.org.cn

    *通信作者:孟憲偉(1963—),男,遼寧省遼中縣人,研究員,博士,從事同位素地球化學(xué)研究。E-mail:mxw@fio.org.cn

    猜你喜歡
    孢粉紅樹欽州
    東北地區(qū)末次盛冰期以來孢粉數(shù)據(jù)庫構(gòu)建技術(shù)與應(yīng)用
    北部灣港欽州自動化集裝箱碼頭/獲獎證書
    水運工程(2022年10期)2022-11-02 05:26:48
    紅樹林誕生記
    西部陸海新通道背景下加強欽州水上消防建設(shè)的思考
    水上消防(2021年4期)2021-11-05 08:51:34
    紅樹植物無瓣海桑中重金屬元素的分布與富集特征
    海洋通報(2020年5期)2021-01-14 09:27:04
    孢粉在流域生態(tài)環(huán)境中的指示作用及研究進展
    人民珠江(2019年10期)2019-11-11 03:22:22
    西藏仲巴縣隆格爾地區(qū)漸新統(tǒng)日貢拉組孢粉組合的發(fā)現(xiàn)及其意義
    欽州坭興陶藝造型設(shè)計初探
    “欽州大蠔”豪氣崛起
    上海師范大學(xué)揭示花粉壁孢粉素合成的調(diào)控機制
    蔬菜(2018年8期)2018-01-16 15:39:27
    亚洲最大成人中文| 亚洲综合色惰| 少妇被粗大猛烈的视频| 久久精品久久久久久噜噜老黄| xxx大片免费视频| 亚洲欧美日韩卡通动漫| 一级av片app| 波多野结衣巨乳人妻| 日本免费在线观看一区| 精品一区二区免费观看| 好男人视频免费观看在线| 免费看日本二区| 看黄色毛片网站| av在线蜜桃| 我的女老师完整版在线观看| 女人十人毛片免费观看3o分钟| 人妻系列 视频| 成人一区二区视频在线观看| 亚洲欧美精品自产自拍| 日韩不卡一区二区三区视频在线| 啦啦啦中文免费视频观看日本| 51国产日韩欧美| 国产日韩欧美在线精品| 国产免费一级a男人的天堂| 亚洲av男天堂| 91久久精品国产一区二区三区| 欧美三级亚洲精品| 搞女人的毛片| 午夜久久久久精精品| 午夜福利视频精品| 久久精品久久精品一区二区三区| 高清欧美精品videossex| av卡一久久| 精品久久久精品久久久| 日本欧美国产在线视频| 精品久久久精品久久久| 亚洲成人久久爱视频| 韩国av在线不卡| 又爽又黄无遮挡网站| 久久精品国产鲁丝片午夜精品| 直男gayav资源| 成人综合一区亚洲| 久久久久免费精品人妻一区二区| 久久亚洲国产成人精品v| 亚洲欧美精品自产自拍| 欧美另类一区| a级毛色黄片| 欧美变态另类bdsm刘玥| 国产精品无大码| 一级二级三级毛片免费看| 三级国产精品片| 午夜免费男女啪啪视频观看| 亚洲欧美精品专区久久| 狂野欧美激情性xxxx在线观看| av女优亚洲男人天堂| 联通29元200g的流量卡| 2022亚洲国产成人精品| 搡老妇女老女人老熟妇| 精品99又大又爽又粗少妇毛片| 国产一区二区三区综合在线观看 | 亚洲欧美一区二区三区国产| 亚洲欧美成人精品一区二区| 2021天堂中文幕一二区在线观| 91久久精品国产一区二区三区| 少妇人妻精品综合一区二区| 日日摸夜夜添夜夜添av毛片| 日韩一本色道免费dvd| 少妇人妻一区二区三区视频| 久久久久久伊人网av| a级毛片免费高清观看在线播放| 成人av在线播放网站| 日韩电影二区| 久久热精品热| 一级毛片久久久久久久久女| 嘟嘟电影网在线观看| 精品少妇黑人巨大在线播放| 日本wwww免费看| 一个人观看的视频www高清免费观看| 狠狠精品人妻久久久久久综合| 免费大片黄手机在线观看| 精品久久久久久久久亚洲| 国产av国产精品国产| 九九爱精品视频在线观看| 不卡视频在线观看欧美| 热99在线观看视频| 一区二区三区免费毛片| 精品国内亚洲2022精品成人| 久久99精品国语久久久| 日本免费a在线| 国产v大片淫在线免费观看| 99久久九九国产精品国产免费| 最新中文字幕久久久久| 99久久精品国产国产毛片| 国产成人a区在线观看| 国产一级毛片七仙女欲春2| 国产亚洲91精品色在线| 久久精品夜色国产| 欧美人与善性xxx| 免费观看a级毛片全部| 成人午夜高清在线视频| 久久久a久久爽久久v久久| 国产一级毛片七仙女欲春2| 一边亲一边摸免费视频| 久久精品久久久久久久性| 国产探花极品一区二区| 亚洲av男天堂| 免费在线观看成人毛片| 亚洲色图av天堂| av免费在线看不卡| 熟女人妻精品中文字幕| 欧美激情久久久久久爽电影| 亚洲自偷自拍三级| a级毛色黄片| 天天一区二区日本电影三级| 日韩三级伦理在线观看| 久久久久性生活片| 黄色日韩在线| 色哟哟·www| 韩国av在线不卡| 免费看光身美女| 亚洲欧美成人综合另类久久久| 久久午夜福利片| 久久久久久久亚洲中文字幕| 国产亚洲5aaaaa淫片| 搡老妇女老女人老熟妇| 日日啪夜夜爽| 人妻夜夜爽99麻豆av| 狂野欧美激情性xxxx在线观看| 国产精品伦人一区二区| 日本免费在线观看一区| 免费看美女性在线毛片视频| 乱码一卡2卡4卡精品| 午夜福利网站1000一区二区三区| 综合色丁香网| 亚洲在线观看片| 日韩精品青青久久久久久| 午夜亚洲福利在线播放| 国产大屁股一区二区在线视频| 国产午夜福利久久久久久| 精品久久国产蜜桃| 内地一区二区视频在线| 精品人妻偷拍中文字幕| 免费大片18禁| 亚洲国产精品国产精品| 久久久亚洲精品成人影院| 偷拍熟女少妇极品色| 特大巨黑吊av在线直播| 人妻夜夜爽99麻豆av| 乱人视频在线观看| 亚洲图色成人| 亚洲av中文av极速乱| 国产黄片视频在线免费观看| 精品一区二区三卡| 亚洲成人精品中文字幕电影| 亚洲婷婷狠狠爱综合网| 久久这里有精品视频免费| 中文在线观看免费www的网站| 欧美三级亚洲精品| 91在线精品国自产拍蜜月| 又爽又黄a免费视频| 成人毛片60女人毛片免费| 一级毛片电影观看| 国产亚洲一区二区精品| 精品国产露脸久久av麻豆 | 国产久久久一区二区三区| 天堂中文最新版在线下载 | 欧美激情在线99| 男女那种视频在线观看| 亚洲熟女精品中文字幕| 在线 av 中文字幕| 国产成人精品福利久久| 国国产精品蜜臀av免费| 高清在线视频一区二区三区| 亚洲成人精品中文字幕电影| 麻豆成人午夜福利视频| 人人妻人人澡人人爽人人夜夜 | 日本三级黄在线观看| 麻豆久久精品国产亚洲av| 亚洲av一区综合| 中文字幕免费在线视频6| av.在线天堂| 搡老乐熟女国产| 成人亚洲精品av一区二区| 人体艺术视频欧美日本| 亚洲av一区综合| 国产大屁股一区二区在线视频| 男女边摸边吃奶| 亚洲国产精品成人综合色| 国产色爽女视频免费观看| 国产精品国产三级专区第一集| 成人午夜高清在线视频| 免费看不卡的av| 最新中文字幕久久久久| 水蜜桃什么品种好| 老司机影院毛片| 国产精品三级大全| videos熟女内射| 午夜日本视频在线| 亚洲精品国产成人久久av| 男人爽女人下面视频在线观看| 久久99热6这里只有精品| 国产 一区 欧美 日韩| 日韩国内少妇激情av| 永久免费av网站大全| 日日撸夜夜添| 91在线精品国自产拍蜜月| 久久精品久久久久久久性| 亚洲人成网站在线观看播放| 午夜福利在线观看免费完整高清在| 日韩制服骚丝袜av| 国产av国产精品国产| 99热这里只有是精品在线观看| 亚洲三级黄色毛片| 国产精品一区二区性色av| 色吧在线观看| 亚洲人与动物交配视频| 精品不卡国产一区二区三区| 最近中文字幕高清免费大全6| 欧美一级a爱片免费观看看| 少妇人妻一区二区三区视频| 51国产日韩欧美| 国产精品一区二区在线观看99 | 在线观看人妻少妇| 91aial.com中文字幕在线观看| 99久久精品国产国产毛片| 99re6热这里在线精品视频| 国产精品国产三级国产av玫瑰| 国产69精品久久久久777片| 在线天堂最新版资源| 国产综合精华液| 天堂网av新在线| 日本午夜av视频| 亚洲国产精品成人久久小说| 大陆偷拍与自拍| 搡老乐熟女国产| 丝袜喷水一区| 搡女人真爽免费视频火全软件| 免费黄色在线免费观看| 欧美三级亚洲精品| 精品久久久久久久末码| 不卡视频在线观看欧美| 久久久精品欧美日韩精品| videos熟女内射| 国产真实伦视频高清在线观看| 男人和女人高潮做爰伦理| www.av在线官网国产| 97在线视频观看| 国产在线一区二区三区精| 汤姆久久久久久久影院中文字幕 | 白带黄色成豆腐渣| 搡老妇女老女人老熟妇| 在线播放无遮挡| 亚洲自偷自拍三级| 又爽又黄a免费视频| 中文欧美无线码| 亚洲四区av| 黄色一级大片看看| 天美传媒精品一区二区| 亚洲欧美成人精品一区二区| 国产一区二区亚洲精品在线观看| 国产免费又黄又爽又色| 亚洲av.av天堂| 免费看av在线观看网站| 国产极品天堂在线| 久久久久精品性色| 天堂中文最新版在线下载 | 亚洲精品日本国产第一区| 免费大片黄手机在线观看| 男女国产视频网站| 午夜免费观看性视频| 日韩制服骚丝袜av| 久久精品久久久久久久性| 亚洲精品aⅴ在线观看| 嫩草影院新地址| 日韩一区二区视频免费看| 国产精品一及| 成年av动漫网址| 久久99蜜桃精品久久| 老司机影院成人| h日本视频在线播放| 国产精品女同一区二区软件| 大又大粗又爽又黄少妇毛片口| 国产精品av视频在线免费观看| 日韩一区二区三区影片| 国产永久视频网站| 久久久久网色| 国产精品一区二区性色av| 国模一区二区三区四区视频| 三级经典国产精品| 欧美zozozo另类| 尾随美女入室| 搡老乐熟女国产| 啦啦啦韩国在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 超碰97精品在线观看| 男女啪啪激烈高潮av片| 国产精品美女特级片免费视频播放器| 国产高潮美女av| 亚洲国产最新在线播放| 男女那种视频在线观看| 男的添女的下面高潮视频| 两个人视频免费观看高清| 国产色婷婷99| 在线观看一区二区三区| 天堂俺去俺来也www色官网 | 99久国产av精品国产电影| 亚洲成色77777| 综合色丁香网| 国产成人一区二区在线| 亚洲精品,欧美精品| 99久国产av精品国产电影| 欧美日韩综合久久久久久| 日韩成人av中文字幕在线观看| 久久久欧美国产精品| 五月玫瑰六月丁香| 老师上课跳d突然被开到最大视频| 免费高清在线观看视频在线观看| 欧美xxxx性猛交bbbb| 欧美日韩一区二区视频在线观看视频在线 | 身体一侧抽搐| 99久久精品国产国产毛片| .国产精品久久| 女人十人毛片免费观看3o分钟| 国产伦一二天堂av在线观看| 乱人视频在线观看| 亚洲无线观看免费| 欧美性感艳星| 高清视频免费观看一区二区 | 网址你懂的国产日韩在线| 观看美女的网站| 可以在线观看毛片的网站| 日韩大片免费观看网站| 中文在线观看免费www的网站| 亚洲av成人精品一二三区| 国产精品蜜桃在线观看| 精华霜和精华液先用哪个| 啦啦啦啦在线视频资源| 久久国产乱子免费精品| 日本黄色片子视频| 美女大奶头视频| a级毛色黄片| 国产有黄有色有爽视频| 精品一区在线观看国产| 色尼玛亚洲综合影院| 精品久久久久久久久亚洲| 又爽又黄无遮挡网站| 国产三级在线视频| 久久精品久久久久久久性| 国产av在哪里看| 久99久视频精品免费| 1000部很黄的大片| 日韩 亚洲 欧美在线| 国产精品一区二区性色av| 国内精品宾馆在线| 永久免费av网站大全| 老司机影院成人| 欧美日韩一区二区视频在线观看视频在线 | 国产一区二区亚洲精品在线观看| 国产午夜精品一二区理论片| 综合色av麻豆| 日韩在线高清观看一区二区三区| 噜噜噜噜噜久久久久久91| 精品国内亚洲2022精品成人| 一边亲一边摸免费视频| 天堂av国产一区二区熟女人妻| 日日撸夜夜添| 成人美女网站在线观看视频| 高清av免费在线| 久久精品熟女亚洲av麻豆精品 | 国产精品久久久久久久久免| 婷婷色综合大香蕉| 亚洲图色成人| 久久久久久伊人网av| 极品少妇高潮喷水抽搐| 亚洲综合色惰| 在线观看人妻少妇| 亚洲av成人av| 成人综合一区亚洲| 久久久久久久久久黄片| 黄色欧美视频在线观看| 日韩av免费高清视频| 精品国内亚洲2022精品成人| 亚洲精品国产成人久久av| 国产精品人妻久久久久久| 嫩草影院精品99| 国产伦理片在线播放av一区| 嫩草影院新地址| 亚洲最大成人中文| 好男人视频免费观看在线| 在线a可以看的网站| 嫩草影院入口| av.在线天堂| 男人和女人高潮做爰伦理| 国产女主播在线喷水免费视频网站 | 国产精品99久久久久久久久| 日韩av在线大香蕉| 熟女电影av网| 成人性生交大片免费视频hd| 男女边摸边吃奶| 国产乱来视频区| av在线观看视频网站免费| 国产精品国产三级国产专区5o| 国产三级在线视频| 久久久久久久国产电影| 七月丁香在线播放| 成年女人在线观看亚洲视频 | av国产久精品久网站免费入址| 小蜜桃在线观看免费完整版高清| 久久久久久久久大av| 午夜精品国产一区二区电影 | 一级片'在线观看视频| 一级爰片在线观看| 国产成人精品婷婷| 国产在线男女| 午夜亚洲福利在线播放| av福利片在线观看| 亚洲最大成人手机在线| 久久久久久九九精品二区国产| 国产中年淑女户外野战色| 日日啪夜夜撸| 午夜福利成人在线免费观看| 色吧在线观看| 午夜视频国产福利| 午夜亚洲福利在线播放| av福利片在线观看| 麻豆成人午夜福利视频| 自拍偷自拍亚洲精品老妇| 国产亚洲一区二区精品| eeuss影院久久| 亚洲国产精品国产精品| 国产国拍精品亚洲av在线观看| 80岁老熟妇乱子伦牲交| 视频中文字幕在线观看| 国产精品1区2区在线观看.| 成人无遮挡网站| 永久网站在线| 97精品久久久久久久久久精品| 国产极品天堂在线| 99热这里只有精品一区| 久久这里只有精品中国| 亚洲人成网站高清观看| or卡值多少钱| 少妇熟女aⅴ在线视频| 久久午夜福利片| 寂寞人妻少妇视频99o| 日韩 亚洲 欧美在线| 91aial.com中文字幕在线观看| 国产91av在线免费观看| 亚洲精品色激情综合| 可以在线观看毛片的网站| av在线天堂中文字幕| 少妇高潮的动态图| 国产精品熟女久久久久浪| 亚洲精品视频女| 免费看日本二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲在线自拍视频| 听说在线观看完整版免费高清| 久久精品国产亚洲av涩爱| 91久久精品国产一区二区成人| 欧美另类一区| 国产伦一二天堂av在线观看| 精品久久久久久久末码| 成人毛片a级毛片在线播放| 久久人人爽人人片av| 3wmmmm亚洲av在线观看| 日本一本二区三区精品| 青春草亚洲视频在线观看| 免费看光身美女| 一级毛片我不卡| 22中文网久久字幕| 一级av片app| 美女国产视频在线观看| 成人一区二区视频在线观看| 日韩成人av中文字幕在线观看| 夫妻性生交免费视频一级片| 国产美女午夜福利| 日韩伦理黄色片| 亚洲av成人精品一区久久| 亚洲最大成人av| 天天一区二区日本电影三级| 99久久精品国产国产毛片| 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 黑人高潮一二区| 小蜜桃在线观看免费完整版高清| 久久6这里有精品| 熟女电影av网| 日韩一区二区三区影片| 免费播放大片免费观看视频在线观看| 免费少妇av软件| 99九九线精品视频在线观看视频| .国产精品久久| 成人性生交大片免费视频hd| 3wmmmm亚洲av在线观看| 欧美一级a爱片免费观看看| 亚洲av电影在线观看一区二区三区 | 免费观看在线日韩| 三级国产精品片| 亚洲怡红院男人天堂| 国产精品福利在线免费观看| 中文在线观看免费www的网站| 久久久精品94久久精品| 看十八女毛片水多多多| 成人av在线播放网站| 97人妻精品一区二区三区麻豆| 老女人水多毛片| 我要看日韩黄色一级片| 三级经典国产精品| 国内精品宾馆在线| 亚洲欧美精品专区久久| 十八禁国产超污无遮挡网站| 嫩草影院精品99| 老司机影院毛片| 最近的中文字幕免费完整| 又大又黄又爽视频免费| 国产探花在线观看一区二区| 国产黄片视频在线免费观看| 亚洲av日韩在线播放| 国产精品伦人一区二区| av又黄又爽大尺度在线免费看| 日韩在线高清观看一区二区三区| 欧美成人一区二区免费高清观看| 亚洲av中文av极速乱| 亚洲不卡免费看| 亚洲欧洲日产国产| 欧美不卡视频在线免费观看| 能在线免费观看的黄片| 成人国产麻豆网| 国产伦精品一区二区三区四那| 久久精品久久精品一区二区三区| 1000部很黄的大片| 国产美女午夜福利| 日产精品乱码卡一卡2卡三| 麻豆乱淫一区二区| 97人妻精品一区二区三区麻豆| 国产成人精品一,二区| 一级毛片久久久久久久久女| 日本-黄色视频高清免费观看| a级毛色黄片| 我的女老师完整版在线观看| 色5月婷婷丁香| 亚洲综合精品二区| 国产黄频视频在线观看| 国产精品麻豆人妻色哟哟久久 | 夫妻性生交免费视频一级片| 秋霞在线观看毛片| 99久久中文字幕三级久久日本| 亚洲无线观看免费| 亚洲最大成人手机在线| 五月伊人婷婷丁香| 国产欧美日韩精品一区二区| 日韩av在线大香蕉| 亚洲精品亚洲一区二区| av黄色大香蕉| 亚洲av电影在线观看一区二区三区 | 纵有疾风起免费观看全集完整版 | 自拍偷自拍亚洲精品老妇| 真实男女啪啪啪动态图| 一个人免费在线观看电影| 国产免费一级a男人的天堂| 午夜精品国产一区二区电影 | 亚洲国产精品成人综合色| 国产乱人视频| av线在线观看网站| 国产精品精品国产色婷婷| 夜夜看夜夜爽夜夜摸| 亚洲一区高清亚洲精品| 又粗又硬又长又爽又黄的视频| 欧美激情国产日韩精品一区| 亚洲av免费高清在线观看| 亚洲成人中文字幕在线播放| 免费人成在线观看视频色| 国产精品久久久久久av不卡| 汤姆久久久久久久影院中文字幕 | 亚洲av成人精品一二三区| 国产精品人妻久久久久久| 久久久久九九精品影院| 欧美日韩在线观看h| 亚洲三级黄色毛片| av卡一久久| 久久99精品国语久久久| 午夜精品一区二区三区免费看| 丝袜美腿在线中文| 国内少妇人妻偷人精品xxx网站| 亚洲国产精品国产精品| 亚洲av一区综合| 91久久精品电影网| 亚洲av不卡在线观看| 在线免费观看的www视频| 赤兔流量卡办理| 天天躁日日操中文字幕| 亚洲av国产av综合av卡| 丝瓜视频免费看黄片| 国产成人精品婷婷| 色5月婷婷丁香| 亚洲最大成人手机在线| 日韩精品有码人妻一区| 人人妻人人看人人澡| 亚洲va在线va天堂va国产| 日韩欧美精品v在线| 成人欧美大片| 一区二区三区乱码不卡18| 国产真实伦视频高清在线观看| 秋霞伦理黄片| 久久久久九九精品影院| 少妇被粗大猛烈的视频| 免费av观看视频| 老司机影院成人| 日韩制服骚丝袜av| 99re6热这里在线精品视频| 国产精品一区www在线观看| 亚洲综合色惰| 成人高潮视频无遮挡免费网站| 亚洲精品第二区|