• <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
    日韩,欧美,国产一区二区三区| 九色成人免费人妻av| 日本av手机在线免费观看| 国产淫片久久久久久久久| 国产精品蜜桃在线观看| 久久精品熟女亚洲av麻豆精品| 看十八女毛片水多多多| 看非洲黑人一级黄片| 深爱激情五月婷婷| 成年免费大片在线观看| 青春草国产在线视频| 国产免费视频播放在线视频| 亚洲第一区二区三区不卡| 亚洲最大成人中文| 国产高清不卡午夜福利| 欧美三级亚洲精品| 国产精品国产三级专区第一集| 免费高清在线观看视频在线观看| 一本—道久久a久久精品蜜桃钙片| 精品久久久精品久久久| 国内揄拍国产精品人妻在线| 日韩av不卡免费在线播放| 18+在线观看网站| 晚上一个人看的免费电影| 最近2019中文字幕mv第一页| 免费黄色在线免费观看| 欧美3d第一页| 久久久久久久久久久丰满| 久久毛片免费看一区二区三区| 国产在线视频一区二区| 日韩伦理黄色片| 校园人妻丝袜中文字幕| 欧美高清性xxxxhd video| 日本黄色日本黄色录像| 一区二区av电影网| 99久久精品国产国产毛片| 最近最新中文字幕大全电影3| 国产精品.久久久| 精华霜和精华液先用哪个| 久久久久久久精品精品| 蜜桃久久精品国产亚洲av| 久久久久国产网址| 久久久成人免费电影| 午夜福利视频精品| 国产亚洲av片在线观看秒播厂| 夫妻性生交免费视频一级片| 最近最新中文字幕免费大全7| 能在线免费看毛片的网站| 国产在线一区二区三区精| 黄色日韩在线| 十分钟在线观看高清视频www | 少妇裸体淫交视频免费看高清| 国产视频内射| 国产成人freesex在线| 日产精品乱码卡一卡2卡三| 亚洲精品国产色婷婷电影| 狂野欧美白嫩少妇大欣赏| 久久婷婷青草| 婷婷色麻豆天堂久久| 国产精品爽爽va在线观看网站| 狂野欧美激情性bbbbbb| 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 超碰av人人做人人爽久久| 国产亚洲一区二区精品| 日韩,欧美,国产一区二区三区| 婷婷色综合大香蕉| av国产久精品久网站免费入址| 直男gayav资源| 黑丝袜美女国产一区| www.色视频.com| 欧美xxxx黑人xx丫x性爽| 尤物成人国产欧美一区二区三区| 国产免费一级a男人的天堂| 国产亚洲91精品色在线| 久久人人爽人人片av| 亚洲欧美成人精品一区二区| 精品久久久久久久久亚洲| 久久6这里有精品| 久久久久久久国产电影| 波野结衣二区三区在线| 国产亚洲91精品色在线| 国产无遮挡羞羞视频在线观看| 99re6热这里在线精品视频| 成人黄色视频免费在线看| 校园人妻丝袜中文字幕| 在线观看人妻少妇| 日本-黄色视频高清免费观看| av又黄又爽大尺度在线免费看| 久久久色成人| 亚洲精品一二三| 只有这里有精品99| 国产精品久久久久久av不卡| 国产深夜福利视频在线观看| 国产免费一级a男人的天堂| 多毛熟女@视频| 韩国av在线不卡| 在线播放无遮挡| 亚洲精品日韩av片在线观看| 一区二区三区乱码不卡18| 国产一区有黄有色的免费视频| 国产成人精品久久久久久| 黑人高潮一二区| 亚洲国产精品成人久久小说| 伊人久久精品亚洲午夜| 七月丁香在线播放| 国产精品嫩草影院av在线观看| 中文字幕久久专区| 91aial.com中文字幕在线观看| 精品少妇久久久久久888优播| 女的被弄到高潮叫床怎么办| 亚洲电影在线观看av| 国产亚洲91精品色在线| 韩国高清视频一区二区三区| 久久久久久久久久人人人人人人| 国产男女内射视频| 亚洲精品456在线播放app| 女人久久www免费人成看片| 久久久久久伊人网av| 99久国产av精品国产电影| 成人一区二区视频在线观看| 国产亚洲5aaaaa淫片| 亚洲精品久久午夜乱码| 亚洲人成网站高清观看| av播播在线观看一区| 国产高清有码在线观看视频| 亚洲人成网站高清观看| 99视频精品全部免费 在线| 国产高清不卡午夜福利| 亚洲人成网站在线播| 18禁动态无遮挡网站| 少妇精品久久久久久久| 欧美+日韩+精品| 美女脱内裤让男人舔精品视频| 一区在线观看完整版| 一级a做视频免费观看| 丰满少妇做爰视频| 晚上一个人看的免费电影| 国产精品人妻久久久久久| 夜夜看夜夜爽夜夜摸| 精华霜和精华液先用哪个| 搡老乐熟女国产| 成年人午夜在线观看视频| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 一区二区三区免费毛片| 久久久久久人妻| av专区在线播放| 黄色一级大片看看| 街头女战士在线观看网站| 亚洲av福利一区| 女人十人毛片免费观看3o分钟| 国产 精品1| 国产深夜福利视频在线观看| 欧美老熟妇乱子伦牲交| 纵有疾风起免费观看全集完整版| 搡女人真爽免费视频火全软件| 久久人人爽av亚洲精品天堂 | 国产精品不卡视频一区二区| 三级国产精品欧美在线观看| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| 欧美丝袜亚洲另类| 小蜜桃在线观看免费完整版高清| 日本黄大片高清| 国产精品免费大片| 99热这里只有精品一区| av网站免费在线观看视频| 综合色丁香网| 精品人妻一区二区三区麻豆| 日本av免费视频播放| 六月丁香七月| 熟女av电影| 亚洲国产精品专区欧美| 国产乱来视频区| 国产黄片美女视频| 久久影院123| 欧美成人午夜免费资源| 亚洲真实伦在线观看| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 26uuu在线亚洲综合色| 亚洲av免费高清在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产亚洲午夜精品一区二区久久| 久久精品国产鲁丝片午夜精品| 国产永久视频网站| 久久午夜福利片| 国产伦精品一区二区三区四那| 久久精品国产自在天天线| 热99国产精品久久久久久7| 国产伦在线观看视频一区| 色吧在线观看| a级毛色黄片| av福利片在线观看| 免费大片18禁| 在线观看人妻少妇| 少妇被粗大猛烈的视频| 丝袜脚勾引网站| 看免费成人av毛片| 色哟哟·www| 日韩一区二区三区影片| 国产成人精品久久久久久| 久久久久久伊人网av| 99久久精品国产国产毛片| 成人一区二区视频在线观看| 亚洲激情五月婷婷啪啪| 插阴视频在线观看视频| 久久亚洲国产成人精品v| 26uuu在线亚洲综合色| 简卡轻食公司| 高清午夜精品一区二区三区| 亚洲色图综合在线观看| 赤兔流量卡办理| 在线观看一区二区三区| 少妇 在线观看| 夫妻性生交免费视频一级片| 偷拍熟女少妇极品色| videos熟女内射| 亚洲激情五月婷婷啪啪| 一本一本综合久久| 人人妻人人澡人人爽人人夜夜| 久久久精品94久久精品| 尤物成人国产欧美一区二区三区| 美女主播在线视频| 免费黄色在线免费观看| 高清黄色对白视频在线免费看 | 特大巨黑吊av在线直播| 国内少妇人妻偷人精品xxx网站| 国精品久久久久久国模美| 国产欧美日韩精品一区二区| 久久久久视频综合| 国产深夜福利视频在线观看| 日本色播在线视频| 在线免费观看不下载黄p国产| 一个人看的www免费观看视频| 激情五月婷婷亚洲| 夜夜爽夜夜爽视频| 七月丁香在线播放| 九九爱精品视频在线观看| 国产黄色免费在线视频| 亚洲欧美中文字幕日韩二区| 亚洲第一区二区三区不卡| 免费观看性生交大片5| 美女中出高潮动态图| 国内揄拍国产精品人妻在线| 亚洲av综合色区一区| 纵有疾风起免费观看全集完整版| 春色校园在线视频观看| 久久久久久久久久久免费av| 最近中文字幕高清免费大全6| 免费人成在线观看视频色| 视频区图区小说| a级毛色黄片| 日本黄大片高清| 91精品伊人久久大香线蕉| 午夜免费男女啪啪视频观看| 有码 亚洲区| 国产日韩欧美亚洲二区| 男女下面进入的视频免费午夜| 久久精品夜色国产| 亚洲av欧美aⅴ国产| av卡一久久| 在线播放无遮挡| 日韩中文字幕视频在线看片 | 在线观看三级黄色| 免费黄网站久久成人精品| 久久久午夜欧美精品| 大香蕉97超碰在线| 女人久久www免费人成看片| 久久久久人妻精品一区果冻| av专区在线播放| 99久久中文字幕三级久久日本| 香蕉精品网在线| 亚洲精品aⅴ在线观看| 97热精品久久久久久| 99re6热这里在线精品视频| 多毛熟女@视频| 观看免费一级毛片| 欧美bdsm另类| 日本黄色片子视频| 2022亚洲国产成人精品| 国产老妇伦熟女老妇高清| 成人亚洲欧美一区二区av| 一级毛片电影观看| 国产 一区 欧美 日韩| 街头女战士在线观看网站| 乱系列少妇在线播放| 肉色欧美久久久久久久蜜桃| 午夜老司机福利剧场| 国产在线视频一区二区| 人人妻人人添人人爽欧美一区卜 | 一本色道久久久久久精品综合| 国产91av在线免费观看| 亚洲精品日韩av片在线观看| 婷婷色av中文字幕| 国产亚洲5aaaaa淫片| a 毛片基地| 最近中文字幕2019免费版| 深爱激情五月婷婷| 赤兔流量卡办理| 国产亚洲91精品色在线| 成人综合一区亚洲| 高清日韩中文字幕在线| 久久久久久人妻| 少妇人妻一区二区三区视频| 久久久久性生活片| 男女边摸边吃奶| 日韩亚洲欧美综合| 欧美日韩综合久久久久久| 91在线精品国自产拍蜜月| 欧美高清成人免费视频www| 国产男女内射视频| av专区在线播放| 国产精品人妻久久久影院| 日日啪夜夜爽| 亚洲伊人久久精品综合| 女的被弄到高潮叫床怎么办| av播播在线观看一区| 国产中年淑女户外野战色| 亚洲av日韩在线播放| 欧美日韩国产mv在线观看视频 | 国产淫片久久久久久久久| 夜夜看夜夜爽夜夜摸| 欧美xxxx黑人xx丫x性爽| 亚洲中文av在线| 午夜免费鲁丝| 国产精品99久久久久久久久| 一级毛片久久久久久久久女| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 久久这里有精品视频免费| 又粗又硬又长又爽又黄的视频| 最近最新中文字幕大全电影3| 久久99热这里只有精品18| 在线亚洲精品国产二区图片欧美 | 男女边摸边吃奶| 丝袜脚勾引网站| 国产女主播在线喷水免费视频网站| 99久久精品国产国产毛片| 99视频精品全部免费 在线| 美女福利国产在线 | 99九九线精品视频在线观看视频| 人妻一区二区av| 久久久久人妻精品一区果冻| 麻豆成人av视频| 激情五月婷婷亚洲| 日韩在线高清观看一区二区三区| 欧美一区二区亚洲| 免费看光身美女| 美女xxoo啪啪120秒动态图| 另类亚洲欧美激情| 久久ye,这里只有精品| 日韩一区二区视频免费看| 中文精品一卡2卡3卡4更新| 老女人水多毛片| 99久久中文字幕三级久久日本| 日日啪夜夜爽| 亚洲一级一片aⅴ在线观看| av在线app专区| 日本vs欧美在线观看视频 | 国产黄片美女视频| 能在线免费看毛片的网站| av天堂中文字幕网| 夜夜骑夜夜射夜夜干| 最近中文字幕2019免费版| 丰满少妇做爰视频| 亚洲欧美精品自产自拍| 黄色视频在线播放观看不卡| 亚洲精品视频女| 久久久久国产精品人妻一区二区| 亚洲av国产av综合av卡| 色视频www国产| 中文欧美无线码| 亚洲成人一二三区av| 少妇人妻一区二区三区视频| 久久久久久久大尺度免费视频| 成人黄色视频免费在线看| 蜜桃亚洲精品一区二区三区| 精品国产乱码久久久久久小说| 中文字幕免费在线视频6| 国产一区二区在线观看日韩| 久久精品国产a三级三级三级| 国产精品av视频在线免费观看| 免费黄网站久久成人精品| 一区二区av电影网| 国产色爽女视频免费观看| 国产欧美另类精品又又久久亚洲欧美| 成人一区二区视频在线观看| 搡老乐熟女国产| 久久av网站| 人人妻人人添人人爽欧美一区卜 | 成人18禁高潮啪啪吃奶动态图 | 如何舔出高潮| 看非洲黑人一级黄片| 精品99又大又爽又粗少妇毛片| 亚洲欧洲国产日韩| 涩涩av久久男人的天堂| 五月玫瑰六月丁香| 99久久精品一区二区三区| 伦理电影免费视频| 在线精品无人区一区二区三 | 肉色欧美久久久久久久蜜桃| 女性被躁到高潮视频| 国产爱豆传媒在线观看| 性色av一级| 亚洲熟女精品中文字幕| 免费黄网站久久成人精品| 少妇人妻久久综合中文| 欧美97在线视频| 日本猛色少妇xxxxx猛交久久| 国产在线男女| 午夜免费鲁丝| 黑丝袜美女国产一区| 性色avwww在线观看| 午夜福利在线观看免费完整高清在| 色婷婷久久久亚洲欧美| 久久精品人妻少妇| 一级av片app| 视频区图区小说| 精品亚洲成a人片在线观看 | 亚洲性久久影院| 一区二区三区四区激情视频| 久久久久久久久久人人人人人人| 女性被躁到高潮视频| 男女边摸边吃奶| 欧美丝袜亚洲另类| 免费播放大片免费观看视频在线观看| 乱码一卡2卡4卡精品| 国产成人精品婷婷| 久久毛片免费看一区二区三区| 青青草视频在线视频观看| 日韩,欧美,国产一区二区三区| 嫩草影院入口| 两个人的视频大全免费| 99热全是精品| 又大又黄又爽视频免费| 亚洲人成网站高清观看| 国产视频首页在线观看| 免费看光身美女| 欧美xxxx黑人xx丫x性爽| 国产美女午夜福利| 永久免费av网站大全| 国产美女午夜福利| 国产真实伦视频高清在线观看| 18禁动态无遮挡网站| www.色视频.com| 成人美女网站在线观看视频| 舔av片在线| 国产在线免费精品| 美女cb高潮喷水在线观看| 纵有疾风起免费观看全集完整版| 久久99热6这里只有精品| 欧美bdsm另类| 最近最新中文字幕大全电影3| 亚洲内射少妇av| 国产伦理片在线播放av一区| 色哟哟·www| 欧美人与善性xxx| 欧美成人a在线观看| av.在线天堂| 最近的中文字幕免费完整| 国产欧美亚洲国产| 18禁裸乳无遮挡免费网站照片| 国产女主播在线喷水免费视频网站| 日韩一区二区三区影片| 日韩中字成人| 亚洲人成网站高清观看| 狂野欧美激情性xxxx在线观看| 国产在线视频一区二区| 我要看黄色一级片免费的| 欧美精品一区二区大全| 建设人人有责人人尽责人人享有的 | 中文欧美无线码| av线在线观看网站| 国产男女内射视频| 亚洲欧美日韩东京热| av卡一久久| 午夜福利网站1000一区二区三区| 天堂俺去俺来也www色官网| 在线观看一区二区三区| 亚洲三级黄色毛片| 亚洲精品久久午夜乱码| 亚洲电影在线观看av| 99热6这里只有精品| av国产免费在线观看| 国产大屁股一区二区在线视频| 国产精品欧美亚洲77777| 国产一区二区三区综合在线观看 | 国产永久视频网站| 中文字幕亚洲精品专区| 国产精品av视频在线免费观看| 欧美97在线视频| 欧美日韩视频高清一区二区三区二| 久久国产精品大桥未久av | 80岁老熟妇乱子伦牲交| av国产精品久久久久影院| av.在线天堂| 99热全是精品| 亚洲国产精品专区欧美| 另类亚洲欧美激情| av卡一久久| 22中文网久久字幕| 91精品国产九色| 97超视频在线观看视频| 国产毛片在线视频| 午夜福利影视在线免费观看| 免费高清在线观看视频在线观看| 97在线人人人人妻| 噜噜噜噜噜久久久久久91| 岛国毛片在线播放| 国产精品久久久久久精品电影小说 | 亚洲av免费高清在线观看| 自拍偷自拍亚洲精品老妇| 蜜桃久久精品国产亚洲av| 九色成人免费人妻av| 日韩,欧美,国产一区二区三区| 毛片女人毛片| 欧美日韩视频精品一区| 又黄又爽又刺激的免费视频.| 91久久精品国产一区二区三区| 亚洲经典国产精华液单| 丝袜脚勾引网站| 黄片wwwwww| 国产淫语在线视频| 一区二区三区四区激情视频| 1000部很黄的大片| 中文字幕制服av| 国产国拍精品亚洲av在线观看| 不卡视频在线观看欧美| 男人爽女人下面视频在线观看| 久久国产乱子免费精品| 99热这里只有是精品50| 久久99热这里只频精品6学生| 亚洲欧美一区二区三区黑人 | 亚洲中文av在线| 亚洲国产精品成人久久小说| 免费人妻精品一区二区三区视频| 欧美丝袜亚洲另类| 在线观看免费日韩欧美大片 | 国产在视频线精品| 自拍偷自拍亚洲精品老妇| 老女人水多毛片| 国产综合精华液| 美女脱内裤让男人舔精品视频| 久久久久性生活片| 一本—道久久a久久精品蜜桃钙片| 国产精品人妻久久久久久| videossex国产| 亚洲真实伦在线观看| 最近最新中文字幕大全电影3| 成人漫画全彩无遮挡| 免费观看av网站的网址| 国产v大片淫在线免费观看| 亚洲国产欧美在线一区| 黄色一级大片看看| 欧美xxxx性猛交bbbb| 亚洲怡红院男人天堂| 97超碰精品成人国产| 国产在视频线精品| 午夜福利在线在线| 老师上课跳d突然被开到最大视频| 久久综合国产亚洲精品| 欧美xxⅹ黑人| 男人狂女人下面高潮的视频| 又粗又硬又长又爽又黄的视频| 大香蕉久久网| 久久精品夜色国产| 99久久中文字幕三级久久日本| 亚洲,欧美,日韩| 女性被躁到高潮视频| 国内揄拍国产精品人妻在线| 人人妻人人看人人澡| 一区二区三区乱码不卡18| 亚洲欧美日韩无卡精品| 久久婷婷青草| 99久久综合免费| 日日啪夜夜撸| 国产成人精品一,二区| 欧美zozozo另类| 日韩av在线免费看完整版不卡| 尤物成人国产欧美一区二区三区| 免费大片18禁| 1000部很黄的大片| 国产久久久一区二区三区| 亚洲av电影在线观看一区二区三区| 美女脱内裤让男人舔精品视频| 一级av片app| 亚洲欧洲日产国产| 99视频精品全部免费 在线| 超碰97精品在线观看| 免费久久久久久久精品成人欧美视频 | 亚洲精品亚洲一区二区| 乱码一卡2卡4卡精品| 欧美一区二区亚洲| 国产精品一区二区性色av| 熟妇人妻不卡中文字幕| 激情 狠狠 欧美| 亚洲欧美精品自产自拍| 熟妇人妻不卡中文字幕| 我要看日韩黄色一级片| 免费观看在线日韩| 久久久久久久国产电影| 久久影院123| 1000部很黄的大片| 国产在线一区二区三区精| 青青草视频在线视频观看| 欧美精品人与动牲交sv欧美| 中文字幕人妻熟人妻熟丝袜美| 啦啦啦中文免费视频观看日本| 男女啪啪激烈高潮av片| 又黄又爽又刺激的免费视频.| 久久国产亚洲av麻豆专区| 日韩人妻高清精品专区| 国产色爽女视频免费观看|