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

    近30年來東海長江口泥質(zhì)區(qū)浮游植物生產(chǎn)力與群落結構變化的生物標志物記錄?

    2020-12-18 03:19:30郭新宇張海龍
    關鍵詞:甲藻陸源硅藻

    郭新宇, 張海龍, 李 莉, 畢 蓉

    (中國海洋大學海洋化學理論與工程技術教育部重點實驗室, 山東 青島 266100)

    東海是我國重要的邊緣海,具有廣闊的大陸架,處于東亞季風的影響區(qū)域。東海環(huán)流體系復雜,主要包括臺灣暖流、長江沖淡水、東海沿岸流、黑潮等。其中,在海底斜坡地貌與東亞季風協(xié)同作用下,東海沿岸流、臺灣暖流等在東海近岸相互作用,形成沿岸上升流[1]。同時,長江三角洲等河口地區(qū)受人類活動影響明顯,長江沖淡水攜帶大量陸源物質(zhì)進入東海。因此,東海的生態(tài)環(huán)境受到人類活動與環(huán)境因素的雙重影響。

    近年來,關于東海近岸泥質(zhì)區(qū)的研究逐漸增多,主要集中于浮游植物生產(chǎn)力與群落結構變化,以及不同來源有機質(zhì)比例變化兩個方面。東海近岸存在著長江口泥質(zhì)區(qū)、閩浙泥質(zhì)區(qū)以及濟州島西南泥質(zhì)區(qū)[2],其中,長江口泥質(zhì)區(qū)位于陸-海作用交界區(qū),受長江輸入的影響明顯。Zhao等[3]通過沉積物中多種色素濃度的變化重建了百年尺度長江口海域浮游植物生產(chǎn)力的變化,表明1980年代以來長江口附近海域浮游植物生產(chǎn)力明顯升高,這與人類活動導致長江輸入的營養(yǎng)鹽濃度増加有關。王江濤等[4]研究了近50年營養(yǎng)鹽變化對長江口海域浮游植物群落結構的影響,結果表明,1980年代以來,長江口硅藻比例降低,甲藻比例升高;引發(fā)這一變化的原因是人類活動導致氮營養(yǎng)鹽過量輸入以及N/P比升髙。Duan等[5]通過研究沉積物中多種生物標志物重建了閩浙泥質(zhì)區(qū)百年尺度浮游植物群落結構的變化,結果表明,1950年代以來,該區(qū)域硅藻比例顯著升高,主要是因為沿岸上升流提供了豐富的硅酸鹽。Cao等[6]運用總有機質(zhì)與生物標志物分析的方法,重建了百年尺度東海陸架區(qū)不同來源有機質(zhì)比值變化,結果表明,1980年代以來,輸入到東海的陸源物質(zhì)逐年降低,主要是因為各種大壩,特別是三峽大壩的攔截作用。綜上可得,目前已有的研究多是聚焦于單一人類活動或自然因素對生態(tài)環(huán)境的影響,對多種人類活動和自然因素對浮游植物生產(chǎn)力、群落結構和陸源輸入相對貢獻綜合影響的研究還較少。

    本文運用總有機質(zhì)與多種生物標志物分析的手段,選擇位于東海長江口泥質(zhì)區(qū)的柱狀樣DH3-1進行研究。通過探討浮游植物生物標志物(菜子甾醇、甲藻甾醇、C37長鏈烯酮)的含量及相對比例、陸源貢獻指標(δ13CTOC、TMBR(Terrestrial and Marine Biomarker Ratio)、BIT(Branched and Isoprenoid Tetraether))的年代記錄,探究DH3-1站位浮游植物生產(chǎn)力、群落結構和陸源輸入相對貢獻的時間變化格局;并結合前人的研究成果,分析人類活動與自然因素產(chǎn)生的綜合影響。

    1 實驗樣品與方法

    1.1 實驗樣品的采集

    本研究于2010年6月搭載國家自然科學基金東海春季共享航次,由科學三號調(diào)查船通過箱式采樣器采集,插管獲得柱狀樣DH3-1。此柱狀樣長50 cm,每隔2 cm取樣,共得25個樣品,樣品置于-20℃保存。采樣站位如圖1所示。

    (陰影部分為泥質(zhì)區(qū),站位DH3-1;CDW:長江沖淡水;KC:黑潮;TWWC:臺灣暖流;ECSCC:東海沿岸流。The shaded areas represent mud depositions. CDW: Changjiang Diluted Water; KC: Kuroshio Current; TWWC: Taiwan Warm Current; ECSCC: the East China Sea Coastal Current.)

    圖1 采樣站位(DH3-1)及東海主要流系(夏季表層)示意圖[7]
    Fig.1 Schematic illustration of the core site (DH3-1) and the average position of currents in the East China Sea (surface layer in summer)[7]

    1.2 實驗方法

    1.2.1210Pb 定年 柱狀樣DH3-1通過210Pb定年,在臺灣中央研究院地球科學研究所完成,用γ分析方法直接測定210Pb及其214Pb(226Ra)。分析儀器為美國ORTEC生產(chǎn)的高純鍺探測器(Ortec HPGe GEM/Lo-Ax/GMX)。利用作圖法,即CIC(Constant Initial210Pb Concentration)模式計算該站位的沉積速率,約為1.8 cm·a-1,選取2010年為測年零點,時間跨度約為27年(1983—2010年)。

    1.2.2 TOC及δ13CTOC分析 沉積物樣品經(jīng)冷凍干燥、研磨后,稱取1 g左右加入玻璃瓶中,滴加4 mol/L 的HCl溶液,去除無機碳;震蕩、超聲、離心后去除上清液,加入5~7 mL Milli-Q水。重復數(shù)次直至檢測上清液為近中性。將洗為近中性的樣品置于55 ℃的烘箱中烘干,烘干后放入干燥器中平衡24 h,之后研磨使其均質(zhì)化。

    用錫杯準確稱取適量去除無機碳的樣品,通過CHN元素分析儀(Thermo Flash 2000)測定總有機碳(TOC)含量(TOC測定的標準偏差為±0.02 wt%(n=6));用錫杯準確稱取適量已去除無機碳的樣品,通過同位素質(zhì)譜儀(Thermo Delta V mass spectrometer),在連續(xù)流模式下測定δ13CTOC(δ13CTOC測定的標準偏差為±0.2‰(n=6))。

    1.2.3 生物標志物分析 生物標志物的分析主要分為樣品預處理及儀器分析兩個步驟。

    (1)樣品預處理:沉積物樣品經(jīng)冷凍干燥、研磨后,稱取5 g左右加入Teflon瓶中,加入CH2Cl2/CH3OH(v/v=3∶1)萃取液和三種內(nèi)標(19醇24氘烷,n-C46GDGT與19酸)。震蕩,超聲,離心,重復萃取4次后N2吹干富集。加入6% KOH的CH3OH溶液進行過夜堿水解。加入正己烷萃取4次后N2吹干富集。然后進行硅膠柱分離,用正己烷淋洗得到非極性組分(烷烴),再用CH2Cl2/CH3OH(v/v =95∶5)混合試劑淋洗得到中性組分(醇類、烯酮和GDGTs)。烷烴組分吹干、加異辛烷定容后上機分析。中性組分平均分成兩份,一份吹干后,加CH2Cl2和BSTFA試劑進行衍生化反應(70 ℃,1 h)后上機分析,用于測定甾醇與烯酮組分。另一份過0.45 μm PTFE濾膜,吹干、定容后上機分析測定GDGTs組分。堿水解后的正己烷萃取剩余液中加入過量6 mol/L HCl溶液,再加入正己烷萃取4次后N2吹干富集;加入HCl/CH3OH(v/v =5∶95)溶液進行酯交換反應(70 ℃,12 h)。加入約2 mL Milli-Q水,用正己烷萃取4次后N2吹干富集。吹干后加異辛烷定容后上機分析,用于測定脂肪酸組分。

    (2)儀器分析:烷烴、甾醇、烯酮、脂肪酸組分使用氣相色譜儀(GC,Agilent 7890A)進行定量分析。使用HP-1型毛細管柱(50 m×0.32 mm×0.17 μm),不分流進樣;載氣為H2,流速為1.3 mL/min。FID檢測器溫度、進樣口溫度和柱箱初始爐溫分別設定為300、300和80 ℃。

    GDGTs組分使用高效液相色譜-質(zhì)譜聯(lián)用儀(HPLC-MS, Agilent 1200/Waters-Quattro Ultima)進行定量分析。液相部分,使用正相鍵合色譜柱(150 mm×2.1 mm),柱溫為30 ℃,流動相為正己烷與異丙醇的混合溶液。采用梯度洗脫,設置流速為0.2 mL/min。質(zhì)譜部分,使用正相色譜柱與大氣壓化學電離源(APCI,atmospheric pressure chemical ionization),采用單離子掃描法。參數(shù)設定條件如下:噴霧器電壓為60 psi,溫度為400 ℃;霧化氮氣流速為600 L/h;錐孔器80 L/h;APCI源溫度為95 ℃,探針溫度為550 ℃。

    2 結果與討論

    2.1 海源生物標志物含量與浮游植物生產(chǎn)力的變化

    圖2(a)展示了過去近30年間(約1983—2010年),DH3-1站位TOC含量與各海源生物標志物含量及其總量的變化。TOC含量的變化范圍是0.24%~0.75%,平均值為0.47%;約1992年以前,TOC含量無明顯變化;約1992年以后開始增加,出現(xiàn)明顯波動。海源生物標志物中,菜子甾醇含量的變化范圍是99~655 ng·g-1,平均值為282 ng·g-1;甲藻甾醇含量的變化范圍是167~929 ng·g-1,平均值為549 ng·g-1;約1983—2003年,長鏈烯酮含量的變化范圍是11~84 ng·g-1,平均值為35 ng·g-1;約2004—2010年,長鏈烯酮含量低于儀器檢出限。三種浮游植物生物標志物含量之和的變化范圍是279~1 533 ng·g-1,平均值為857 ng·g-1。泉古菌醇含量的變化范圍是294~1 076 ng ·g-1,平均值為552 ng·g-1。

    ((1)TOC:總有機碳;(2)B:菜子甾醇含量;(3)D:甲藻甾醇含量;(4)A:C37長鏈烯酮含量;(5)∑(B+D+A):菜子甾醇+甲藻甾醇+C37長鏈烯酮總量;(6)Crenarcheaol:泉古菌醇含量;(7)B/TOC:菜子甾醇含量與TOC比值;(8)D/TOC:甲藻甾醇含量與TOC比值;(9)A/TOC:C37長鏈烯酮含量與TOC比值;(10)∑(B+D+A)/TOC:菜子甾醇+甲藻甾醇+C37長鏈烯酮總量與TOC比值;(11)Crenarcheaol /TOC:泉古菌醇含量與TOC比值。虛線往右代表發(fā)生明顯變化的時間。藍色箭頭代表變化的趨勢。(1) TOC: total organic carbon; (2) B: brassicasterol content; (3) D: dinosterol content; (4) A: C37alkenones content; (5) ∑(B+D+A): brassicasterol + dinosterol + C37alkenones contents; (6) Crenarcheaol: crenarcheaol content; (7) B/TOC: TOC-normalized brassicasterol content; (8) D/TOC: TOC-normalized dinosterol content; (9) A/TOC: TOC-normalized C37alkenones content; (10) ∑(B+D+A)/TOC: TOC-normalized brassicasterol + dinosterol + C37alkenones contents; (11) Crenarcheaol/TOC: TOC-normalized crenarcheaol content. Dashed lines represent the time at which major changes of biomarkers started. Blue arrows represent the trend of changes.)

    圖2 柱狀樣DH3-1海源生物標志物含量(ng·g-1)(a)及與其TOC比值(ng·g-1·%-1)(b)的變化
    Fig.2 Marine biomarker contents (ng·g-1) (a) and TOC-normalized biomarker contents (ng·g-1·%-1) (b) in the core DH3-1

    生物標志物含量的變化會受到生產(chǎn)力變化、降解作用、沉積速率等因素的影響[8]。將各種生物標志物的含量與TOC做比值,可以消除成巖作用、沉積速率與粒度的影響(見圖2(b))。從圖中可以看出,三種浮游植物生物標志物含量、生物標志物含量/TOC比值的變化趨勢基本一致。由此可以看出,在過去近30年間,浮游植物生物標志物含量的變化主要受生產(chǎn)力變化影響。

    雖然生物標志物含量并不能嚴格代表其來源藻類的絕對數(shù)量,但是其變化可以代表對應藻類數(shù)量的相對變化[9],這是用生物標志物含量的變化重建生產(chǎn)力變化的基礎。菜子甾醇、甲藻甾醇、C37長鏈烯酮分別主要由硅藻、甲藻、顆石藻產(chǎn)生,因此,其含量的變化可以分別用來指示硅藻、甲藻、顆石藻生產(chǎn)力的變化[10];總量的變化可以用來指示浮游植物總生產(chǎn)力的變化。從圖2中可以看出,菜子甾醇、甲藻甾醇、C37長鏈烯酮含量及其含量總和均呈現(xiàn)增加趨勢,這表明在過去近30年間,DH3-1站位的硅藻生產(chǎn)力、甲藻生產(chǎn)力、顆石藻生產(chǎn)力與浮游植物總生產(chǎn)力均呈現(xiàn)上升趨勢。DH3-1站位生產(chǎn)力的重建結果與前人在該海域的研究結果具有較好的一致性。Xing等[11]運用生物標志物法重建了位于東海近岸泥質(zhì)區(qū)的五個站位(CJ43,DH5-1,DH6-3,32#,34#)浮游植物生產(chǎn)力在百年尺度上的變化;結果表明:近百年尤其是1980年代以來,五個站位浮游植物生產(chǎn)力的變化均呈現(xiàn)上升趨勢。

    2.2 浮游植物生物標志物比例與浮游植物群落結構的變化

    圖3表示了DH3-1站位過去近30年間浮游植物生物標志物相對比例(每種生物標志物與總生物標志物含量之比)的變化。菜子甾醇相對比例的變化范圍是0.19~0.43,平均值為0.32;甲藻甾醇相對比例的變化范圍是0.57~0.76,平均值為0.64;C37長鏈烯酮相對比例的變化范圍是0.03~0.09,平均值為0.05。

    由于生物標志物母源細胞的產(chǎn)生數(shù)量及其在環(huán)境中的留存效率不確定[17],所以生物標志物的比例不能嚴格代表其來源藻類的群落占比,但是其比例的變化可以指示其來源藻類群落占比的變化。如圖3所示,約1992年以前,硅藻、甲藻的相對比例沒有明顯變化。約1992—2001年,硅藻相對比例降低,甲藻相對比例升高;約2001年以后,兩者呈現(xiàn)出相反的變化趨勢,即硅藻相對比例升高,甲藻相對比例降低。顆石藻相對比例一直在相對低值波動較低且無明顯變化規(guī)律。

    DH3-1站位的群落結構重建結果與前人在該海域的研究結果具有較好的一致性。Xing等[11]通過生物標志物法重建了東海近岸泥質(zhì)區(qū)的五個站位(CJ43,DH5-1,DH6-3,32#,34#)百年尺度上浮游植物群落結構的變化,結果表明:東海不同位置硅藻與甲藻相對比例的變化趨勢大致可以分為三種類型:①硅藻比例降低,甲藻比例升高;②硅藻比例升高,甲藻比例降低;③硅藻甲藻比例交替增減。DH3-1站位群落結構的重建結果符合“硅藻甲藻比例交替增減”的特征,即,硅藻比例先降低后升高,甲藻比例先升高后降低。對比DH3-1站位與CJ43等五個站位的位置與群落結構的變化,可以初步得出:在東海不同區(qū)域,硅藻與甲藻相對比例的變化明顯不同(見圖4)。在泥質(zhì)區(qū)北部及外側,呈現(xiàn)硅藻比例降低,甲藻比例升高(站位CJ43,DH5-1)和硅藻甲藻比例交替增減(站位DH3-1,DH6-3)的變化;在泥質(zhì)區(qū)內(nèi)側,呈現(xiàn)硅藻比例升高,甲藻比例降低(站位32#,34#)的變化。又如,王云龍等[18]直接統(tǒng)計了長江口及附近海域2000—2003年20個監(jiān)測點春季浮游植物的數(shù)量與種類的水平分布與變化情況。結果顯示:從2000—2003年,硅藻比例依次為:78.78%,75%,75.61%,76.09%;甲藻比例依次為:12.12%,22.73%,12.20%,10.87%,即,以2001年為分界,呈現(xiàn)硅藻比例先下降后上升、甲藻比例先上升后下降的變化特征。這一結果表明,長江口鄰近海域群落組成的鏡檢結果與本文的重建結果具有較為一致的變化特征。

    ((1)B/∑MB:菜子甾醇比例;(2)D/∑MB:甲藻甾醇比例;(3)A/∑MB:C37長鏈烯酮比例;B:菜子甾醇含量;D:甲藻甾醇含量;A:C37長鏈烯酮含量;∑MB:浮游植物生物標志物總量(菜子甾醇+甲藻甾醇+C37長鏈烯酮)。虛線往右代表發(fā)生明顯變化的時間。藍色箭頭代表變化的趨勢。(1) B/∑MB: the ratio between brassicasterol and the sum of three biomarkers; (2) D/∑MB: the ratio between dinosterol and the sum of three biomarkers; (3) A/∑MB: the ratio between C37alkenones and the sum of three biomarkers; B: brassicasterol content; D: dinosterol content; A: C37alkenones content; ∑MB: the sum of three phytoplankton biomarkers (brassicasterol + dinosterol + C37alkenones). Dashed lines represent the time at which major changes of biomarkers started. Blue arrows represent the trend of changes.)

    圖3 柱狀樣DH3-1浮游植物生物標志物與總生物標志物含量之比的變化
    Fig.3 Marine phytoplankton biomarker ratios in the core DH3-1

    (站位CJ43,DH5-1,DH6-3,32#,34#信息引用自文獻[11]。Site information for CJ43,DH5-1,DH6-3,32#,34# is cited from reference[11].)

    圖4 東海硅藻與甲藻相對比例呈現(xiàn)不同變化特征的站位分布
    Fig.4 The distribution of sites with different ratios of diatoms and dinoflagellates

    前人研究表明,營養(yǎng)鹽結構的改變會對浮游植物群落結構產(chǎn)生顯著影響[19]。Si、P含量高的情況下,Si∶N增大而N∶P減小,硅藻更具競爭優(yōu)勢[20];N含量尤其是DIN含量高的情況下,Si∶N減小而N∶P增大,甲藻更具競爭優(yōu)勢[21]。根據(jù)楊穎[15]等統(tǒng)計的1984—2014年30年期間長江口鄰近海域中營養(yǎng)鹽含量的數(shù)據(jù),可以發(fā)現(xiàn):過去近30年期間,長江口海域水質(zhì)中N∶P 比值在不同年份存在明顯的轉折變化。其中,約1992年以前,N∶P基本無明顯的變化趨勢,而2000年以后N∶P總體低于1987—2000年間的值。N∶P的這一變化較好地解釋了本研究浮游植物群落結構的演替規(guī)律,即約1992年以前,硅藻、甲藻的相對比例沒有明顯變化,而約2001年以后,硅藻相對比例升高,甲藻相對比例降低。

    營養(yǎng)鹽結構的上述變化則可能是由于不同年份,驅(qū)動其變化的機制不同。長江每年都向東海輸送大量的營養(yǎng)鹽,是東海營養(yǎng)鹽的主要來源。長江流域地區(qū)復雜多變的人類活動使長江輸入東海的營養(yǎng)鹽含量明顯改變,致使東海營養(yǎng)鹽結構發(fā)生變化,從而導致東海浮游植物群落結構改變。特別是1980年以來,長江流域大規(guī)模修建水壩、開展水土保持工程等人類活動則導致東海Si輸入量明顯減少[22]。由于水壩的攔截作用會導致長江入海的沉積物通量明顯減少[23],從而導致Si的入海通量減少[24],而N的入海通量卻并沒有顯著減少。這是因為長江流域水庫大多屬于“攔沙不攔水”型水庫[25],對P、Si(多以顆粒態(tài)存在)的攔截較為明顯,對N(多以溶解態(tài)存在)的攔截作用較弱。除了人類活動,東海沿岸上升流也會對營養(yǎng)鹽結構產(chǎn)生顯著影響。上升流可以促進深層水向上涌動,為上層水體帶來更多有機物[26],尤其是輸入了高含量的Si、P營養(yǎng)鹽,從而導致營養(yǎng)鹽結構呈現(xiàn)Si∶N增大,N∶P減小的變化,導致群落結構向硅藻比例增大的趨勢演變。Chen等[27]研究了渤、黃、東海不同季節(jié)表層海水溫度、鹽度、氮磷硅營養(yǎng)鹽濃度分布,結果表明:東海沿岸存在的高硅酸鹽濃度區(qū),與沿岸上升流的分布區(qū)域基本一致。硅營養(yǎng)鹽濃度高值區(qū)的存在同時也證明了上升流的存在。由于上升流帶來的深層水溫度較低,因此,可以將海水次表層溫度作為判定上升流是否存在的指標。重建DH3-1站位次表層海水溫度(TEX86-T)的變化,以此判定DH3-1站位在過去約30年間是否存在上升流(見圖5)。重建結果顯示,約1999年以后,DH3-1站位的海水次表層溫度呈現(xiàn)降低趨勢,表明約1999年以后研究區(qū)域有明顯的上升流存在。

    (虛線代表發(fā)生明顯變化的時間。The dashed line represents the time at which major changes started.)

    圖5 DH3-1站位次表層海水溫度變化
    Fig.5 TEX86-T of core DH3-1

    綜上可得,人類活動和上升流(自然因素)會對營養(yǎng)鹽結構產(chǎn)生截然不同的影響。人類活動加劇會導致營養(yǎng)鹽結構呈現(xiàn)Si∶N減小,N∶P增大的變化;上升流則會導致營養(yǎng)鹽結構呈現(xiàn)Si∶N增大,N∶P減小的變化。因此,DH3-1站位硅藻與甲藻比例交替變化,可能是因為不同年份人類活動與自然因素的作用強度不同。1992—2001年,DH3-1站位硅藻比例下降、甲藻比例上升,可能是因為人類活動的影響大于自然因素;2001年以后,DH3-1站位硅藻比例上升、甲藻比例下降,可能是因為自然因素的影響大于人類活動。通過對比Xing[11]等提出的東海近岸不同區(qū)域硅藻與甲藻比例變化的驅(qū)動機制與DH3-1站位硅藻與甲藻比例變化的影響因素,可以發(fā)現(xiàn):在長江口附近,硅藻與甲藻比例變化主要受人類活動驅(qū)動;但是在該海域的小范圍內(nèi),某些年份自然因素仍可以產(chǎn)生顯著影響。

    從圖2中可以看出,DH3-1站位的顆石藻比例一直處于低值。前人研究表明,長江口及其鄰近海域的浮游植物群落主要由沿岸廣溫型的硅藻和甲藻組成[28],這與本文的研究結果一致。造成這種群落結構分布主要是由于顆石藻生長的環(huán)境與硅藻、甲藻明顯不同,高鹽、高溫、低營養(yǎng)鹽環(huán)境更利于顆石藻生長[29]。而且,現(xiàn)代調(diào)查顯示,東海近岸的顆石藻比例隨著離岸距離和水深的增加而增大[30]。近30年來,人類活動加劇導致東海營養(yǎng)鹽水平顯著提升,且DH3-1站位位于東海近岸,這些因素都導致DH3-1站位顆石藻相對比例一直較低。

    2.3 陸源生物標志物含量與陸源輸入相對貢獻的變化

    圖6(a)表示了DH3-1站位過去近30年間TOC含量與各陸源生物標志物含量及其總量的變化。其中,C27+C29+C31長鏈正構烷烴含量的變化范圍是327~1 164 ng·g-1,平均值為531 ng·g-1。C24+C26+C28+C30一元正構飽和直鏈脂肪酸含量的變化范圍是0.95~3.03 μg·g-1,平均值為1.81 μg·g-1。支鏈GDGTs含量的變化范圍是27~92 ng·g-1,平均值為56 ng·g-1。將各種陸源生物標志物的含量與TOC做比值,可以消除成巖作用、沉積速率與粒度的影響,如圖6(b)所示。三種陸源生物標志物含量、生物標志物含量/TOC比值的變化趨勢基本一致。由此可以看出,在過去近30年間,三種陸源生物標志物含量的改變主要是由于受物質(zhì)輸入的影響,而不受降解等其他因素的影響。

    長鏈正構烷烴以及一元正構飽和直鏈脂肪酸主要產(chǎn)生于陸源高等植物表層葉蠟[31],可以通過風塵及河流顆粒物輸入海洋。沉積物中長鏈正構烷烴(以C27、C29、C31為主)記錄可以反演陸源高等植被有機質(zhì)的輸入[32];C24~C36一元正構偶數(shù)碳飽和直鏈脂肪酸同樣可以指示陸源有機質(zhì)輸入[33]。支鏈GDGTs主要由土壤厭氧細菌產(chǎn)生,通過河流輸入海洋,其含量可以用來指示陸源土壤有機質(zhì)輸入[34]。

    (1)∑(C27+C29+C31n-alkanes):C27+C29+C31長鏈正構烷烴含量;(2)Branched-GDGTs:支鏈GDGTs含量;(3)∑(C24+C26+C28+C30n-SFAs):C24+C26+C28+C30一元正構飽和直鏈脂肪酸含量;(4)∑(C27+C29+C31n-alkanes)/TOC:C27+C29+C31長鏈正構烷烴含量與TOC 比值;(5)Branched-GDGTs/TOC:支鏈GDGTs含量與TOC比值;(6)∑(C24+C26+C28+C30n-SFAs)/TOC:C24+C26+C28+C30一元正構飽和直鏈脂肪酸含量與TOC比值。虛線往右代表發(fā)生明顯變化的時間。(1) C27+C29+C31long-chainn-alkanes contents; (2) Branched-GDGTs contents; (3) C24+C26+C28+C30n-saturated fatty acids contents; (4) TOC-normalized C27+C29+C31n-alkanes contents; (5) TOC-normalized branched-GDGTs contents; (6) TOC-normalized C24+C26+C28+C30n-SFAs contents. Dashed lines represent the time of at which major changes started. Blue arrows represent the trend of changes.)

    圖6 柱狀樣DH3-1陸源生物標志物含量(ng·g-1或μg·g-1) (a) 及與TOC比值 (ng·g-1·%-1或μg·g-1·%-1)(b)的變化
    Fig.6 Terrestrial biomarker contents (ng·g-1or μg·g-1)(a) and TOC-normalized biomarker contents (ng·g-1·%-1or μg·g-1·%-1)(b) in the core DH3-1

    從圖6中可以看出,三種陸源生物標志物含量具有相似的變化特點,即約1992年以前,無明顯變化;約1992—2004年,波動明顯,稍有增加;約2004年以后,明顯減少。這表明約1992年以前,陸源有機質(zhì)輸入量沒有明顯變化;約1992—2004年,稍有增加。約2004年以后明顯減少。

    為了區(qū)分海源與陸源有機質(zhì)的相對比例,指示陸源輸入的相對貢獻,利用不同來源生物標志物含量的比值構建 TMBR指標[35]和BIT指標[34]。TMBR指標是將指示陸源有機質(zhì)的長鏈正構烷烴含量(C27+C29+C31n-alkanes)與指示海源有機質(zhì)的浮游植物生物標志物含量(菜子甾醇+甲藻甾醇+長鏈烯酮)作比值,計算公式為:TMBR = ∑(C27+C29+C31n-alkanes)/[∑(C27+C29+C31n-alkanes)+ ∑(菜子甾醇+甲藻甾醇+長鏈烯酮)]。TMBR的取值范圍為0~1。TMBR為0時,有機質(zhì)全部為海源輸入;TMBR為1時,有機質(zhì)全部為陸源輸入。BIT指標是將指示陸源土壤有機質(zhì)的三種支鏈 GDGTs含量(I+II+III)與指示海源有機質(zhì)的泉古菌醇(IV)含量作比值,計算公式為:BIT=(I+II+III)/(I+II+III+IV)。BIT的取值范圍為0~1。BIT為0時,有機質(zhì)全部為海源輸入;BIT為1時,有機質(zhì)全部為陸源輸入。

    除了TMBR與BIT指標,δ13CTOC也可以用來評估陸源有機質(zhì)的相對貢獻,通常而言,陸源有機質(zhì)δ13CTOC比海源有機質(zhì)δ13CTOC更負[36]。

    圖7表示了DH3-1站位過去近30年間陸源輸入相對貢獻的變化。其中,δ13CTOC的變化范圍是-23.17‰~-22.71‰,平均值為-22.80‰。TMBR的變化范圍是0.21~0.77,平均值為0.45。BIT的變化范圍是0.07~0.15,平均值為0.09。從圖中可以看出,過去近30年間,δ13CTOC整體變正;TMBR、BIT整體減小;這表明在過去近30年間,不同方式輸入的陸源有機質(zhì)的相對貢獻均呈現(xiàn)降低趨勢。

    ((1)δ13CTOC;(2)TMBR(=∑(C27+C29+C31n-alkanes)/[∑(C27+C29+C31n-alkanes)+ ∑(菜子甾醇+甲藻甾醇+長鏈烯酮)];(3)BIT=支鏈GDGTs/(支鏈GDGTs +泉古菌醇)。藍色箭頭代表變化的趨勢。(1) δ13CTOC; (2) TMBR(=∑(C27+C29+C31n-alkanes)/[∑(C27+C29+C31n-alkanes)+∑(brassicasterol + dinosterol + C37alkenones)]; (3)BIT=branched-GDGTs/(branched-GDGTs + crenarcheaol). Blue arrows represent the trend of changes.)

    圖7 柱狀樣DH3-1陸源輸入相對貢獻指標的變化
    Fig.7 Terrestrial relative contribution index in the core DH3-1

    DH3-1站位陸源輸入相對貢獻的重建結果與前人在該海域的研究結果具有較好的一致性。Cao等[6]通過研究總有機質(zhì)與生物標志物重建了東海近岸三個站位(DH5-1,DH6-2,DH6-3)百年尺度上陸源輸入相對貢獻的變化情況;重建結果表明:近百年來,尤其是1980年代以后,三個站位的陸源輸入相對貢獻均明顯降低。

    分析導致陸源輸入相對貢獻變化的原因,需要綜合考慮影響陸源有機質(zhì)與海源有機質(zhì)輸入的原因。一是海源有機質(zhì)輸入方面:近30年以來,工農(nóng)業(yè)迅速發(fā)展導致氮磷類化肥使用量劇增,使得東海營養(yǎng)鹽含量迅速增加,加之沿岸上升流的影響,導致海洋初級生產(chǎn)力提高,海源有機質(zhì)輸入量增加。二是陸源有機質(zhì)輸入方面:1980年代以來,長江流域大規(guī)模修建水壩與開展水土保持工程等人類活動,造成長江入海的沉積物通量明顯減少,導致通過長江向東海輸入的陸源有機質(zhì)明顯減少[37]。但是,由于近百年來東亞冬季風增強[38],搬運了更多的陸源物質(zhì)輸入海洋[39],導致通過大氣向東海輸入的陸源有機質(zhì)明顯增加。人類活動和季風搬運(自然因素)會對陸源有機質(zhì)輸入量產(chǎn)生截然不同的影響,即人類活動加劇會導致陸源有機質(zhì)輸入量減少,而季風搬運增強則會導致陸源有機質(zhì)輸入量增加。

    從圖7中可以看出,DH3-1站位陸源生物標志物的含量在約2004年前后交替變化,即約1992—2004年稍有增加,約2004年以后明顯減少。這可能是因為不同年份人類活動與自然因素的作用強度不同,從而導致陸源有機質(zhì)輸入量呈現(xiàn)不同的變化。約1992—2004年,陸源有機質(zhì)輸入量稍有增加,可能是因為自然因素的影響大于人類活動;約2004年以后,陸源有機質(zhì)輸入量明顯減少,可能是因為人類活動的影響大于自然因素。

    圖2的結果則表明,過去近30年間,DH3-1站位海源有機質(zhì)輸入量整體增加,這主要是因為近30年以來,工農(nóng)業(yè)迅速發(fā)展導致氮磷類化肥使用量劇增,使得東海營養(yǎng)鹽含量迅速增加,加之沿岸上升流的影響,導致DH3-1站位的初級生產(chǎn)力提高。

    綜上可知,在約2004年以前,陸源有機質(zhì)與海源有機質(zhì)含量均呈現(xiàn)增加的趨勢,但是由于兩者的控制因素存在差異,導致陸、海源有機質(zhì)增加的程度不一樣。海源有機質(zhì)增加的比重大于陸源有機質(zhì),導致陸源有機質(zhì)輸入的相對貢獻減少。約2004年以后,海源有機質(zhì)含量呈現(xiàn)增加的趨勢,陸源有機質(zhì)含量呈現(xiàn)減少的趨勢,導致陸源有機質(zhì)輸入的相對貢獻減少。

    3 結論

    (1)菜子甾醇、甲藻甾醇、長鏈烯酮含量變化的記錄顯示,在過去的近30年間,DH3-1站位浮游植物生產(chǎn)力呈現(xiàn)上升趨勢。

    (2)約1992年以來,DH3-1站位的群落結構發(fā)生明顯改變。約1992—2001年期間硅藻比例下降、甲藻比例上升;約2001年以后,硅藻比例上升、甲藻比例下降。

    (3)約1992年以來,DH3-1站位的陸源輸入相對貢獻發(fā)生明顯改變。δ13CTOC、TMBR、BIT三個指標指示的陸源輸入相對貢獻均呈現(xiàn)減小趨勢。長鏈正構烷烴、一元正構飽和直鏈脂肪酸、支鏈 GDGTs含量變化的記錄顯示,約1992—2004年期間,陸源有機質(zhì)輸入量增加,與此同時,海源有機質(zhì)輸入量也增加且增加比重大于陸源有機質(zhì),所以陸源輸入相對貢獻減??;約2004年以后,陸源有機質(zhì)輸入量減少,海源有機質(zhì)輸入量依然增加,所以陸源有機質(zhì)輸入的相對貢獻減小。

    致謝:本文所使用的210Pb定年數(shù)據(jù)由臺灣中央研究院地球科學研究所陳蔚芳老師提供,作者對此表示誠摯謝意。

    猜你喜歡
    甲藻陸源硅藻
    甲藻甾醇在海洋和湖泊古環(huán)境重建中的應用
    硅藻——海里的“生物寶石”
    大自然探索(2023年7期)2023-11-14 13:07:36
    硅藻實驗室“打工記”
    放學
    北京文學(2023年3期)2023-03-13 01:34:00
    硅藻指數(shù)在水環(huán)境監(jiān)測與評價中的應用分析
    養(yǎng)只恐龍當寵物
    學與玩(2019年8期)2019-10-29 03:32:16
    疑似甲藻水華監(jiān)測及應急管理探究
    陸源有機碳對萊州灣浮游動物能量貢獻的初步研究
    控制陸源污染,保護海洋環(huán)境
    硅藻泥涂料的現(xiàn)狀和發(fā)展趨勢
    上海建材(2018年1期)2018-04-18 12:15:16
    久久热精品热| 中文字幕人妻熟人妻熟丝袜美| 久久久色成人| 精品国内亚洲2022精品成人| 欧美日韩一区二区视频在线观看视频在线 | 国产精品蜜桃在线观看| 国产大屁股一区二区在线视频| 久久久久免费精品人妻一区二区| 亚洲欧美精品专区久久| 一级二级三级毛片免费看| 2021少妇久久久久久久久久久| 国产久久久一区二区三区| 91精品一卡2卡3卡4卡| 中文天堂在线官网| 亚洲欧洲日产国产| 韩国av在线不卡| 国产精品电影一区二区三区| 两个人的视频大全免费| a级毛片免费高清观看在线播放| 好男人在线观看高清免费视频| 日本三级黄在线观看| 免费看光身美女| av在线观看视频网站免费| 国产成年人精品一区二区| 最后的刺客免费高清国语| av线在线观看网站| 波多野结衣巨乳人妻| 1000部很黄的大片| 久久久久久大精品| 桃色一区二区三区在线观看| 亚洲性久久影院| 蜜臀久久99精品久久宅男| 亚洲欧洲日产国产| 精品无人区乱码1区二区| 亚洲丝袜综合中文字幕| 亚洲精品久久久久久婷婷小说 | 精品人妻一区二区三区麻豆| 亚洲电影在线观看av| 国产精品一及| 亚洲欧美日韩高清专用| 可以在线观看毛片的网站| 国产免费福利视频在线观看| 亚洲国产欧洲综合997久久,| 能在线免费观看的黄片| 国产精品不卡视频一区二区| 国产极品天堂在线| 美女xxoo啪啪120秒动态图| 少妇人妻精品综合一区二区| 熟妇人妻久久中文字幕3abv| 久久99蜜桃精品久久| 乱人视频在线观看| 在现免费观看毛片| 亚洲精品456在线播放app| 成人一区二区视频在线观看| 国产精品熟女久久久久浪| 最近手机中文字幕大全| 成人av在线播放网站| 免费黄色在线免费观看| 女人十人毛片免费观看3o分钟| or卡值多少钱| 国产精品熟女久久久久浪| 午夜久久久久精精品| 熟女电影av网| 麻豆国产97在线/欧美| 99在线视频只有这里精品首页| 国产免费又黄又爽又色| 日韩av在线免费看完整版不卡| 亚洲国产欧美人成| 在线免费观看不下载黄p国产| 久久久久性生活片| 亚洲图色成人| 亚洲真实伦在线观看| 国产高清视频在线观看网站| 最近视频中文字幕2019在线8| 97热精品久久久久久| 99热这里只有是精品在线观看| 少妇人妻精品综合一区二区| 麻豆精品久久久久久蜜桃| 好男人在线观看高清免费视频| 嘟嘟电影网在线观看| 成人美女网站在线观看视频| 久久精品国产亚洲av涩爱| 精品一区二区三区人妻视频| eeuss影院久久| 久久人人爽人人爽人人片va| 欧美zozozo另类| 国产黄色视频一区二区在线观看 | 九草在线视频观看| 国产成人精品久久久久久| 日韩人妻高清精品专区| 久久这里只有精品中国| 成人漫画全彩无遮挡| 岛国在线免费视频观看| 中文乱码字字幕精品一区二区三区 | 国产高清视频在线观看网站| 一级二级三级毛片免费看| 一个人免费在线观看电影| 国产精品熟女久久久久浪| 国产精品伦人一区二区| 国产一区有黄有色的免费视频 | 欧美高清性xxxxhd video| 中文字幕久久专区| 亚洲三级黄色毛片| 日本熟妇午夜| 精品一区二区三区人妻视频| 国产精品福利在线免费观看| 久久久久久伊人网av| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩卡通动漫| 好男人在线观看高清免费视频| 久久精品国产99精品国产亚洲性色| 高清在线视频一区二区三区 | 欧美精品国产亚洲| 嫩草影院新地址| 日韩亚洲欧美综合| 久久久国产成人精品二区| 国产日韩欧美在线精品| 亚洲性久久影院| 国产极品精品免费视频能看的| 日产精品乱码卡一卡2卡三| 国产一区二区在线观看日韩| 亚洲熟妇中文字幕五十中出| 国产黄色小视频在线观看| 日本一二三区视频观看| 最近手机中文字幕大全| 麻豆乱淫一区二区| 亚洲天堂国产精品一区在线| 天堂网av新在线| 人人妻人人看人人澡| 白带黄色成豆腐渣| 日韩一区二区三区影片| 一级毛片aaaaaa免费看小| 99久久九九国产精品国产免费| 成人毛片60女人毛片免费| av.在线天堂| 黄片无遮挡物在线观看| 亚州av有码| 美女xxoo啪啪120秒动态图| 欧美成人一区二区免费高清观看| 三级毛片av免费| 别揉我奶头 嗯啊视频| 久久亚洲国产成人精品v| 国产av码专区亚洲av| 国产极品天堂在线| 色吧在线观看| 淫秽高清视频在线观看| 久久午夜福利片| 国产精品美女特级片免费视频播放器| 一卡2卡三卡四卡精品乱码亚洲| 一级爰片在线观看| 国产伦理片在线播放av一区| 日本黄大片高清| 久久精品国产自在天天线| 99久久九九国产精品国产免费| 两个人的视频大全免费| 久久这里只有精品中国| 亚洲自偷自拍三级| av视频在线观看入口| 国产人妻一区二区三区在| 蜜桃久久精品国产亚洲av| 欧美高清成人免费视频www| 99视频精品全部免费 在线| 自拍偷自拍亚洲精品老妇| 亚洲av福利一区| 午夜福利在线观看吧| 如何舔出高潮| 午夜免费男女啪啪视频观看| 国语自产精品视频在线第100页| 男的添女的下面高潮视频| 亚洲欧洲国产日韩| 偷拍熟女少妇极品色| 久久久久九九精品影院| 蜜桃亚洲精品一区二区三区| 赤兔流量卡办理| 九九在线视频观看精品| 欧美日韩在线观看h| 国产黄a三级三级三级人| 亚洲aⅴ乱码一区二区在线播放| 欧美xxxx性猛交bbbb| 国产在线一区二区三区精 | 神马国产精品三级电影在线观看| 国产精品久久久久久精品电影| 久久精品久久精品一区二区三区| 日韩一区二区三区影片| 99热精品在线国产| av.在线天堂| 美女黄网站色视频| 97热精品久久久久久| 欧美成人精品欧美一级黄| 18+在线观看网站| 国产免费福利视频在线观看| 国产成人91sexporn| 校园人妻丝袜中文字幕| 久久精品夜色国产| 观看免费一级毛片| 国产精品一区二区三区四区免费观看| 亚洲av中文av极速乱| 天堂√8在线中文| 国内精品宾馆在线| 直男gayav资源| 国内精品一区二区在线观看| 亚洲成人精品中文字幕电影| 狂野欧美激情性xxxx在线观看| 国产精品美女特级片免费视频播放器| 国产精品人妻久久久久久| 国产三级在线视频| 乱人视频在线观看| 秋霞在线观看毛片| av专区在线播放| 中文天堂在线官网| 国产亚洲91精品色在线| 国产精品嫩草影院av在线观看| 久久久a久久爽久久v久久| 熟女电影av网| 久久99热这里只有精品18| 亚洲欧美日韩卡通动漫| 国产精品久久久久久久久免| 欧美人与善性xxx| 午夜福利在线在线| 国产亚洲午夜精品一区二区久久 | 日日啪夜夜撸| 精品久久久久久久久av| 2021少妇久久久久久久久久久| 国产久久久一区二区三区| 久久鲁丝午夜福利片| 只有这里有精品99| 噜噜噜噜噜久久久久久91| av女优亚洲男人天堂| 黄色欧美视频在线观看| 91午夜精品亚洲一区二区三区| 成年女人永久免费观看视频| 非洲黑人性xxxx精品又粗又长| 精品国内亚洲2022精品成人| 插逼视频在线观看| 免费电影在线观看免费观看| 日韩高清综合在线| 国产成人精品一,二区| 精品国产露脸久久av麻豆 | 在线天堂最新版资源| 七月丁香在线播放| 精品久久久久久久末码| 中文字幕亚洲精品专区| 色综合色国产| 久久婷婷人人爽人人干人人爱| 91精品伊人久久大香线蕉| 高清在线视频一区二区三区 | 亚洲av福利一区| 日本欧美国产在线视频| 成人美女网站在线观看视频| 永久免费av网站大全| 国产av不卡久久| 国产片特级美女逼逼视频| 欧美97在线视频| 久久精品人妻少妇| 久久久久久大精品| 少妇熟女欧美另类| 亚洲精华国产精华液的使用体验| 亚洲精品自拍成人| 一级毛片久久久久久久久女| 亚洲精品456在线播放app| 亚洲经典国产精华液单| 日本wwww免费看| 亚洲婷婷狠狠爱综合网| 天堂影院成人在线观看| 插逼视频在线观看| 少妇的逼好多水| 三级毛片av免费| 如何舔出高潮| 国产黄片美女视频| www.色视频.com| 国产精华一区二区三区| 欧美一级a爱片免费观看看| 国产女主播在线喷水免费视频网站 | 少妇高潮的动态图| 中文字幕熟女人妻在线| 亚洲欧美中文字幕日韩二区| 亚洲人与动物交配视频| 国产探花在线观看一区二区| 在线天堂最新版资源| 亚洲综合色惰| 久久韩国三级中文字幕| 国产中年淑女户外野战色| 国产激情偷乱视频一区二区| 欧美高清性xxxxhd video| 国产伦在线观看视频一区| 亚洲精品,欧美精品| 99九九线精品视频在线观看视频| 精品无人区乱码1区二区| 国产激情偷乱视频一区二区| 亚洲精品aⅴ在线观看| 色吧在线观看| 一本久久精品| 亚洲av电影不卡..在线观看| 伊人久久精品亚洲午夜| 亚洲精品日韩av片在线观看| 舔av片在线| 人妻少妇偷人精品九色| 联通29元200g的流量卡| 国产极品精品免费视频能看的| 最新中文字幕久久久久| 亚洲成人精品中文字幕电影| 婷婷色综合大香蕉| 国产成人免费观看mmmm| 九色成人免费人妻av| 久久久久久大精品| 亚洲av电影在线观看一区二区三区 | 欧美xxxx性猛交bbbb| 青春草国产在线视频| 在线播放国产精品三级| 国产高清三级在线| 国产亚洲精品av在线| av福利片在线观看| 边亲边吃奶的免费视频| 午夜福利在线在线| 国产私拍福利视频在线观看| 乱人视频在线观看| 久久韩国三级中文字幕| 亚洲国产精品成人久久小说| 国产久久久一区二区三区| 看黄色毛片网站| 日本黄色视频三级网站网址| 久久欧美精品欧美久久欧美| 国产 一区精品| 色吧在线观看| 久久6这里有精品| 国产精品久久久久久精品电影小说 | 国产真实伦视频高清在线观看| 丰满人妻一区二区三区视频av| 日韩三级伦理在线观看| 深爱激情五月婷婷| 男插女下体视频免费在线播放| 晚上一个人看的免费电影| 真实男女啪啪啪动态图| 熟妇人妻久久中文字幕3abv| 日韩欧美 国产精品| 亚洲自拍偷在线| 亚洲欧美成人综合另类久久久 | 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦韩国在线观看视频| 91精品一卡2卡3卡4卡| 亚洲色图av天堂| 国内精品美女久久久久久| 男女啪啪激烈高潮av片| av国产免费在线观看| 久久人人爽人人爽人人片va| 免费观看性生交大片5| 男人和女人高潮做爰伦理| 51国产日韩欧美| 国产在视频线在精品| 日韩,欧美,国产一区二区三区 | 女人十人毛片免费观看3o分钟| 亚洲欧美中文字幕日韩二区| 亚洲精品456在线播放app| 又黄又爽又刺激的免费视频.| 免费一级毛片在线播放高清视频| 欧美日韩国产亚洲二区| 熟女电影av网| 久久国产乱子免费精品| 中国国产av一级| 日韩亚洲欧美综合| 久久草成人影院| 少妇的逼水好多| 欧美一区二区精品小视频在线| 亚洲综合色惰| 久久久久久国产a免费观看| 一区二区三区乱码不卡18| 精品久久久久久成人av| 少妇的逼水好多| 日韩一本色道免费dvd| www.色视频.com| 两个人视频免费观看高清| 午夜免费激情av| av免费观看日本| 九九久久精品国产亚洲av麻豆| АⅤ资源中文在线天堂| 少妇猛男粗大的猛烈进出视频 | 国产伦一二天堂av在线观看| 你懂的网址亚洲精品在线观看 | 成人午夜高清在线视频| 国产精华一区二区三区| 18禁在线播放成人免费| 少妇熟女欧美另类| 天堂网av新在线| 白带黄色成豆腐渣| 麻豆成人av视频| av免费观看日本| 精品欧美国产一区二区三| 成人综合一区亚洲| 免费看a级黄色片| 视频中文字幕在线观看| 国产91av在线免费观看| 久久精品国产亚洲av涩爱| 亚洲精品久久久久久婷婷小说 | 男女国产视频网站| 最近的中文字幕免费完整| 精品国内亚洲2022精品成人| 国产一级毛片七仙女欲春2| 美女国产视频在线观看| 久久人人爽人人片av| 日本免费在线观看一区| 69人妻影院| 少妇的逼水好多| 免费不卡的大黄色大毛片视频在线观看 | 日韩av在线大香蕉| kizo精华| 一边亲一边摸免费视频| 亚洲精品乱码久久久v下载方式| 神马国产精品三级电影在线观看| 亚洲国产最新在线播放| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 久久精品国产鲁丝片午夜精品| 欧美日韩在线观看h| 国产成人aa在线观看| 精品国产露脸久久av麻豆 | 九九热线精品视视频播放| 中文在线观看免费www的网站| 美女大奶头视频| 国产美女午夜福利| 天堂中文最新版在线下载 | 亚洲色图av天堂| 欧美一区二区亚洲| 免费观看精品视频网站| 久久久久久久久大av| 亚洲国产成人一精品久久久| 日产精品乱码卡一卡2卡三| 亚洲国产精品成人久久小说| 全区人妻精品视频| 国产免费视频播放在线视频 | 久久人人爽人人爽人人片va| 少妇熟女欧美另类| 永久网站在线| 欧美zozozo另类| av在线亚洲专区| 国产精品三级大全| 天堂影院成人在线观看| 26uuu在线亚洲综合色| 日本一二三区视频观看| 亚洲最大成人中文| 国产亚洲最大av| 国产成人精品一,二区| av国产免费在线观看| 最近的中文字幕免费完整| 日本免费在线观看一区| 日日干狠狠操夜夜爽| 欧美高清成人免费视频www| 国产黄a三级三级三级人| 亚洲人成网站高清观看| 搡女人真爽免费视频火全软件| 美女高潮的动态| 日本五十路高清| 好男人在线观看高清免费视频| 色网站视频免费| 青春草亚洲视频在线观看| 国产淫语在线视频| 国产精品国产三级国产av玫瑰| 干丝袜人妻中文字幕| 久久精品人妻少妇| 久久久久网色| 精品国产露脸久久av麻豆 | 九九久久精品国产亚洲av麻豆| 国产乱来视频区| 午夜福利网站1000一区二区三区| 男插女下体视频免费在线播放| 小说图片视频综合网站| 91在线精品国自产拍蜜月| 在线播放无遮挡| 色哟哟·www| 午夜精品一区二区三区免费看| 国产精品一及| 国产亚洲av片在线观看秒播厂 | 汤姆久久久久久久影院中文字幕 | 一二三四中文在线观看免费高清| 久久鲁丝午夜福利片| 国产精品女同一区二区软件| 免费一级毛片在线播放高清视频| 国产高清视频在线观看网站| 亚洲天堂国产精品一区在线| www.色视频.com| 99久久成人亚洲精品观看| 国产伦精品一区二区三区视频9| 午夜精品国产一区二区电影 | 国产视频内射| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 日韩欧美 国产精品| 日韩亚洲欧美综合| 夫妻性生交免费视频一级片| 亚洲欧美日韩卡通动漫| 色5月婷婷丁香| 91久久精品国产一区二区三区| 成人一区二区视频在线观看| 亚洲内射少妇av| 免费观看精品视频网站| 国产成人精品婷婷| 直男gayav资源| 26uuu在线亚洲综合色| 嘟嘟电影网在线观看| 蜜臀久久99精品久久宅男| 国产乱来视频区| 村上凉子中文字幕在线| 热99re8久久精品国产| 免费一级毛片在线播放高清视频| 久久久久久久久久黄片| 丝袜喷水一区| 在线免费观看的www视频| 观看免费一级毛片| 免费av毛片视频| 极品教师在线视频| 国产亚洲5aaaaa淫片| 韩国av在线不卡| 汤姆久久久久久久影院中文字幕 | 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美日韩东京热| 岛国在线免费视频观看| 国产精品乱码一区二三区的特点| 中国国产av一级| 中国美白少妇内射xxxbb| 日本爱情动作片www.在线观看| 日韩,欧美,国产一区二区三区 | 亚洲色图av天堂| 一区二区三区免费毛片| 夫妻性生交免费视频一级片| 精品久久久久久久末码| 精品国产露脸久久av麻豆 | 搡老妇女老女人老熟妇| 国产探花极品一区二区| 两个人视频免费观看高清| 久久亚洲国产成人精品v| 91aial.com中文字幕在线观看| 最近最新中文字幕免费大全7| 午夜激情福利司机影院| 亚洲av二区三区四区| 日韩av在线大香蕉| 国产精品国产高清国产av| av卡一久久| 99视频精品全部免费 在线| 性色avwww在线观看| 国产一级毛片七仙女欲春2| 成人av在线播放网站| 久久热精品热| 桃色一区二区三区在线观看| 欧美日韩在线观看h| 成人亚洲精品av一区二区| 久久亚洲精品不卡| 国产 一区 欧美 日韩| 日韩三级伦理在线观看| 成人欧美大片| 乱人视频在线观看| 日本熟妇午夜| 色播亚洲综合网| 国产熟女欧美一区二区| 国内精品宾馆在线| 国产乱人视频| 欧美3d第一页| eeuss影院久久| 国产精品永久免费网站| 丰满人妻一区二区三区视频av| 国产淫语在线视频| 激情 狠狠 欧美| 伊人久久精品亚洲午夜| 性插视频无遮挡在线免费观看| 小说图片视频综合网站| 精品欧美国产一区二区三| 午夜日本视频在线| 久久久国产成人精品二区| 色噜噜av男人的天堂激情| 中文字幕人妻熟人妻熟丝袜美| 国产成人精品婷婷| 国产成人午夜福利电影在线观看| 校园人妻丝袜中文字幕| 亚洲经典国产精华液单| 男人的好看免费观看在线视频| 91久久精品电影网| 国产亚洲一区二区精品| 男人舔女人下体高潮全视频| 国产 一区 欧美 日韩| 国产一区二区在线av高清观看| 亚洲国产日韩欧美精品在线观看| 免费观看人在逋| 国产精品1区2区在线观看.| 亚洲av免费高清在线观看| 亚洲精品日韩在线中文字幕| 一级av片app| 久久久久久久久大av| 亚洲丝袜综合中文字幕| 丰满乱子伦码专区| 天天一区二区日本电影三级| 一边摸一边抽搐一进一小说| 一区二区三区高清视频在线| 久久精品夜色国产| 精华霜和精华液先用哪个| 级片在线观看| 精品久久久久久久久av| 欧美激情国产日韩精品一区| 亚洲精品影视一区二区三区av| 寂寞人妻少妇视频99o| 99热全是精品| 26uuu在线亚洲综合色| 十八禁国产超污无遮挡网站| 久久久精品94久久精品| 内射极品少妇av片p| 一区二区三区免费毛片| 国产女主播在线喷水免费视频网站 | 国产av不卡久久| 亚洲欧美成人精品一区二区| 国产高清有码在线观看视频| 国产老妇女一区| 亚洲欧美一区二区三区国产| 亚洲欧美日韩无卡精品| 亚洲精品一区蜜桃| 变态另类丝袜制服| av在线亚洲专区| 国模一区二区三区四区视频| 欧美一区二区精品小视频在线|