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

    南極普里茲灣及其鄰近海域懸浮顆粒有機(jī)物的碳同位素組成及其影響因素

    2015-06-01 12:30:13任春燕陳敏高眾勇郭勞動賈仁明劉蕭邱雨生鄭敏芳
    海洋學(xué)報(bào) 2015年12期
    關(guān)鍵詞:普里南極表層

    任春燕,陳敏,2*,高眾勇,郭勞動,4,賈仁明,劉蕭,邱雨生,鄭敏芳

    (1.廈門大學(xué)海洋與地球?qū)W院,福建廈門 361102;2.廈門大學(xué)近海海洋環(huán)境科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,福建廈門 361102;3.國家海洋局第三海洋研究所大氣與全球變化重點(diǎn)實(shí)驗(yàn)室,福建廈門 361005;4.美國威斯康星大學(xué)水科學(xué)學(xué)院,密爾沃基WI 53204)

    南極普里茲灣及其鄰近海域懸浮顆粒有機(jī)物的碳同位素組成及其影響因素

    任春燕1,陳敏1,2*,高眾勇3,郭勞動1,4,賈仁明1,劉蕭1,邱雨生1,鄭敏芳1

    (1.廈門大學(xué)海洋與地球?qū)W院,福建廈門 361102;2.廈門大學(xué)近海海洋環(huán)境科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,福建廈門 361102;3.國家海洋局第三海洋研究所大氣與全球變化重點(diǎn)實(shí)驗(yàn)室,福建廈門 361005;4.美國威斯康星大學(xué)水科學(xué)學(xué)院,密爾沃基WI 53204)

    依托中國第29次南極科學(xué)考察航次開展了南大洋普里茲灣及其鄰近海域懸浮顆粒有機(jī)物碳同位素組成(δ13CPOC)的研究,結(jié)合溫度、鹽度、營養(yǎng)鹽和溶解CO2的數(shù)據(jù),揭示了影響研究海域顆粒有機(jī)物碳同位素組成的主控因素,計(jì)算出混合層中浮游植物吸收無機(jī)碳過程的碳同位素分餾因子。結(jié)果表明,普里茲灣及其鄰近海域的δ13CPOC介于-28.5‰~-21.1‰,平均值為-24.6‰,表現(xiàn)出灣內(nèi)大于灣外的特征。浮游植物同化吸收CO2過程的碳同位素分餾是影響研究海域混合層δ13CPOC的主要因素,根據(jù)δ13CPOC和1/[CO2(aq)]的線性擬合關(guān)系,計(jì)算出浮游植物同化吸收CO2過程的碳同位素分餾因子εp為23.4‰。δ13CPOC的垂直分布隨深度增加而增大,反映出顆粒有機(jī)物垂向輸送過程中顆粒有機(jī)物再礦化過程同位素分餾作用的影響。

    懸浮顆粒有機(jī)物;δ13CPOC;碳同位素分餾因子;南極普里茲灣

    1 引言

    海洋顆粒有機(jī)物(POM)是生物泵運(yùn)轉(zhuǎn)的載體,POM從產(chǎn)生到沉降遷出至海底的過程中,一方面可將部分金屬元素和污染物攜帶至沉積物中,另一方面,沉降過程中有機(jī)物的降解可將部分生源要素等釋放至水體環(huán)境中。天然碳穩(wěn)定同位素組成可用來反映海洋生物泵各種生物地球化學(xué)過程的信息,包括生物光合作用、營養(yǎng)鹽的生物吸收、有機(jī)物質(zhì)的降解、顆粒物從真光層的輸出、有機(jī)物質(zhì)的不同來源等。真光層顆粒有機(jī)物的碳同位素組成變化可提供有關(guān)海洋生物固碳過程的相關(guān)信息,有助于深入理解海洋環(huán)境中碳的生物地球化學(xué)循環(huán)。對諸多河口的研究發(fā)現(xiàn),陸源和海源δ13CPOC有不同的端元值,這為區(qū)分陸源和海源有機(jī)物的貢獻(xiàn)奠定了基礎(chǔ)[1—3]。對于開闊大洋,已發(fā)現(xiàn)不少海域δ13CPOC和溶解CO2濃度([CO2(aq)])具有負(fù)相關(guān)關(guān)系[4—9],這為古海洋學(xué)研究中通過δ13CPOC反演海水表層溶解CO2濃度的變化提供了可能[10—11]。

    南大洋是典型的高營養(yǎng)鹽、低葉綠素海區(qū)(HN-LC),南大洋約80%海域上層水體中的葉綠素a含量都小于0.5 mg/m3,生產(chǎn)力較低[12],但它是全球大氣重要的冷源和CO2匯,其吸收的CO2占全球海洋凈吸收通量的35%,因而成為研究海洋碳循環(huán)及其與全球氣候變化關(guān)系的重要區(qū)域[13]。南大洋顆粒有機(jī)物的碳同位素組成與中低緯度海域相比有明顯不同的特征,懸浮顆粒有機(jī)物和沉積有機(jī)物的δ13CPOC均明顯低于中低緯度海域的相應(yīng)值[14—15]。研究表明,南大洋浮游植物體內(nèi)有機(jī)物的δ13C值隨著緯度的增加而降低,這種降低起初被認(rèn)為是溫度變化引起光合作用過程碳同位素分餾的變化所致,但此后的研究證實(shí),溫度引起的溶解CO2濃度變化,進(jìn)而導(dǎo)致光合作用過程碳同位素分餾因子變化才是關(guān)鍵因素[4,16]。也有研究指出,南大洋較低的δ13CPOC值是由低溫、低光照、高溶解CO2濃度和低浮游植物生長速率共同所致[17]。

    普里茲灣位于南大洋印度洋扇區(qū),是南極大陸沿岸面積僅次于威德爾海和羅斯海的第三大海灣,也是我國南極科學(xué)考察的重點(diǎn)海區(qū)。由于埃默里冰架等的隔離作用,普里茲灣受陸源輸入和人類活動的影響較少。迄今為止,有關(guān)普里茲灣及其鄰近海域顆粒有機(jī)物碳同位素組成的研究還比較少,且主要集中在表層和真光層。Popp等研究了伊麗莎白公主地至普里茲灣東側(cè)一個(gè)南北向斷面懸浮顆粒有機(jī)物碳同位素的控制因素,發(fā)現(xiàn)δ13CPOC與溶解CO2濃度([CO2(aq)])之間沒有相關(guān)性,提出浮游植物種類及其生長速率的變化可能是調(diào)控δ13CPOC變化的主要因素[18]。Zhang等研究了普里茲灣上層水體碳動力學(xué)與水團(tuán)構(gòu)成的關(guān)系,發(fā)現(xiàn)海冰融化水份額與δ13CPOC之間具有顯著的正相關(guān)關(guān)系[19]。尹希杰等報(bào)道了普里茲灣及其鄰近海域表層海水δ13CPOC的空間分布特征及其影響因素,認(rèn)為研究海域δ13CPOC的分布可能主要受浮游植物吸收CO2過程的影響[20]。

    本研究利用中國第29次南極科學(xué)考察航次的機(jī)會,開展了普里茲灣及其鄰近海域整個(gè)水柱中懸浮顆粒有機(jī)物碳同位素組成的研究,結(jié)合溫度、鹽度、營養(yǎng)鹽和溶解CO2的數(shù)據(jù),探討懸浮顆粒有機(jī)物從表及底碳同位素組成空間變化的影響因素,揭示研究海域浮游植物光合作用過程的碳同位素分餾及其分餾因子。

    2 實(shí)驗(yàn)方法

    2.1 樣品采集

    懸浮顆粒物樣品采自中國第29次南極科學(xué)考察航次,該航次由“雪龍”號科考船實(shí)施,采樣時(shí)間為2013年1月31日至3月3日(南半球夏季),利用CTD采集了表層至3 500 m深度區(qū)間不同深度的海水樣品。研究海區(qū)的海冰季節(jié)變化顯著,一般而言,10月為海冰融化初期,11至12月是融化中期,1月至2月是海冰融化末期[21],本研究樣品采集時(shí)間屬于海冰融化末期。

    樣品采集區(qū)域最北到64.5°S,站位主要分布在5個(gè)經(jīng)向斷面,分別為P3(68°E)、P4(70.5°E)、P5(73° E)、P6(75°E)和P7(78°E附近)斷面,共計(jì)27個(gè)站位(圖1)。P3斷面采集了5個(gè)站位的樣品,站位北起65°S、最南至67.5°S的陸架區(qū),靠近達(dá)恩利角冰間湖。P4斷面由65°S向南延伸至陸架區(qū)弗拉姆淺灘處(67.5°S),共設(shè)有5個(gè)站位。P5斷面由65°S沿73°E向南穿過普里茲灣灣口延伸至普里茲灣灣內(nèi)68.5°S,此斷面為歷次南極科學(xué)考察的重點(diǎn)斷面,共設(shè)有6個(gè)站位。P6斷面由普里茲灣灣外66°S向南穿過四女士淺灘延伸至普里茲灣灣內(nèi)69.2°S,共設(shè)有6個(gè)站位。P7斷面由64.5°S向南穿過四女士淺灘至68°S,共設(shè)有5個(gè)站位。

    圖1 普里茲灣及其鄰近海域懸浮顆粒有機(jī)物碳同位素組成研究的采樣站位Fig.1 Sampling locations in the Prydz Bay and its adjacent areas during 2013

    2.2 方法

    海水中懸浮顆粒物樣品由CTD采集的海水樣品(5~10 dm3)經(jīng)過預(yù)先高溫灼燒的47 mm GF/F濾膜(450℃,4 h)過濾后得到。將含有顆粒物的濾膜于60℃下烘干,裝入濾膜盒,于-20℃下冷凍保存,帶回陸地實(shí)驗(yàn)室進(jìn)行POC及δ13CPOC的分析。在陸地實(shí)驗(yàn)室,濾膜干燥后用濃鹽酸酸熏48 h,去除無機(jī)碳,用Milli-Q水潤洗烘干后,包入錫舟送入元素分析儀-同位素比值質(zhì)譜儀聯(lián)機(jī)測定(EA-IRMS)(Carlo Erba NC 2500,Thermo Finnigan DELTAPLUSXP)。EAIRMS按如下條件進(jìn)行設(shè)置:氧化爐溫度為1 000℃,還原爐溫度為650℃;色譜柱(Porapak Q)溫度為50℃;He載氣流量為100 cm3/min;O2流量為150 cm3/min;參考?xì)釩O2流量為150 cm3/min。顆粒物中的有機(jī)碳通過催化氧化作用轉(zhuǎn)化為CO2,進(jìn)而測得POC含量及δ13CPOC。POC含量測定所采用的標(biāo)準(zhǔn)物質(zhì)為C6H6N2O,檢出限為0.1μmol,測量精度優(yōu)于0.2%。碳穩(wěn)定同位素的分析標(biāo)準(zhǔn)為美國南卡羅萊州白堊系皮狄組擬箭石(PDB,Peedee Belemnite)。為確保測量過程中儀器的穩(wěn)定性及所得δ13C值的準(zhǔn)確,每間隔10個(gè)待測樣品穿插一個(gè)IAEA-C8標(biāo)準(zhǔn)物質(zhì)(δ13C=-18.3‰)進(jìn)行單點(diǎn)校正,IAEA-C8標(biāo)準(zhǔn)物質(zhì)測定的相對標(biāo)準(zhǔn)偏差小于0.2‰。

    溶解CO2濃度根據(jù)DIC、p H和TA計(jì)算得到(CO2 SYS V2.1)[22],其中,DIC和TA分別采用無機(jī)碳分析儀(DIC Analyzer AS-C3)和總堿度滴定儀(Total Alkalinity Gran Titration System,AS-ALK1)進(jìn)行測量[23],DIC和TA的測定精度均優(yōu)于0.1%。p H用Thermo Orion ROSS 復(fù)合p H電極(型號8102BN)及Thermo Orion 3 star p H計(jì)測得,測量精度優(yōu)于0.005 p H單位(0.1 m V)。計(jì)算中所用碳酸離解常數(shù)K1和K2采用Mehrbach等報(bào)道[24]并經(jīng)Dickson和Millero修正的數(shù)值[25]。

    營養(yǎng)鹽數(shù)據(jù)由國家海洋局第二海洋研究所韓正兵、潘建明研究員提供。海水樣品經(jīng)0.45μm醋酸纖維膜過濾后用飽和氯化汞溶液固定,保存于潔凈的塑料瓶中,帶回陸地實(shí)驗(yàn)室。分別用銅鎘柱還原法、抗壞血酸磷鉬藍(lán)法測定硝酸鹽、活性磷酸鹽的含量(《海洋監(jiān)測規(guī)范》GB 173784-2007)。溫度、鹽度數(shù)據(jù)由SBE-911 plus CTD直接獲得,其中溫度測量的精度為±0.001℃,電導(dǎo)率的測定精度為±0.002 S/m。

    3 結(jié)果

    3.1 溫度和鹽度

    海水的表層溫度介于-1.49~0.77℃,平均值為-0.37℃。表層溫度的分布整體呈現(xiàn)東北-西南走向,灣內(nèi)溫度最低(P5-11站),由灣內(nèi)向?yàn)惩庵饾u增高,至深海區(qū)表層溫度達(dá)到最高(P7-05站)。表層水的鹽度介于32.78~34.13,平均值為33.75,最低值出現(xiàn)在P6-13站,最高值出現(xiàn)在P3-07站。表層鹽度的分布表現(xiàn)出灣內(nèi)低于灣外、灣外西側(cè)高于東側(cè)的特征。研究區(qū)域5個(gè)斷面溫鹽的垂直分布特征十分類似,故以P3和P5斷面為例加以說明(下同)。從普里茲灣灣內(nèi)外溫度、鹽度的垂直分布可以看出,研究區(qū)域從表到底主要存在南極夏季表層水(AASW),冬季殘留水(WW),繞極深層水(CDW)和南極底層水(AABW),CDW由灣外涌升至陸坡區(qū)。陸架區(qū)受海冰融化水影響,溫度、鹽度低于陸坡及深海區(qū),深海區(qū)近底層出現(xiàn)AABW的特征[26](圖2、圖3)。

    3.2 營養(yǎng)鹽

    普里茲灣及其鄰近海域表層海水硝酸鹽濃度的變化范圍為17.52~31.26μmol/dm3,平均值為27.16μmol/dm3,最低值出現(xiàn)在灣內(nèi)的P6-13站,最高值出現(xiàn)在陸坡區(qū)的P4-05站。表層海水活性磷酸鹽濃度的變化范圍為1.16~2.56μmol/dm3,平均值為2.17μmol/dm3?;钚粤姿猁}和硝酸鹽的變化趨勢一致,陸架區(qū)由于浮游植物光合作用吸收營養(yǎng)鹽更為強(qiáng)烈,硝酸鹽和磷酸鹽濃度明顯低于陸坡區(qū)和深海區(qū)。整個(gè)研究區(qū)域硝酸鹽和活性磷酸鹽的垂直變化呈隨深度增加而增加的趨勢,中深層水體變化相對較小(圖2、圖3)。

    3.3 CO2(aq)

    普里茲灣及其鄰近海域表層CO2(aq)含量介于12.6~21.8μmol/kg之間,平均值為18.5μmol/kg,最低值出現(xiàn)在灣內(nèi)的P6-13,最高值出現(xiàn)在陸坡區(qū)的P3-07站。整個(gè)研究區(qū)域CO2(aq)含量的變化趨勢呈現(xiàn)由灣內(nèi)向?yàn)惩庵饾u增加的特征;CO2(aq)的垂直分布則表現(xiàn)為隨深度增加而增大,但中深層變化較小的規(guī)律,反映出真光層CO2(aq)被浮游植物吸收利用的影響(圖2、圖3)。

    3.4 POC

    普里茲灣及其鄰近海域表層水POC含量介于0.86~10.83μmol/dm3之間,平均值為3.85μmol/dm3,變化范圍較大,最高值出現(xiàn)在營養(yǎng)鹽和CO2(aq)濃度最低的P6-13站,最低值出現(xiàn)在陸坡區(qū)弗拉姆淺灘處的P4-07站,總體低于此前中國南極科學(xué)考察航次報(bào)道的普里茲灣及其鄰近海域表層水的POC含量(1.06~93.67μmol/dm3)[27—29],可能與采樣站位的空間分布以及采樣時(shí)間的不同有關(guān)。本研究在灣內(nèi)采樣站位較少,特別是未能采集到埃默里冰架附近海域的樣品,而冰架邊緣海域的浮游植物生長一般較為旺盛,具有較高的初級生產(chǎn)力[30—31],另外,本研究樣品的采集時(shí)間處于海冰融化末期,比其他航次來得晚,此時(shí)浮游植物的藻華已處于后期。從POC的空間分布看,陸架區(qū)POC濃度高于陸坡及深海區(qū),而中深層海水POC含量變化較?。▓D2、圖3)。

    圖2 普里茲灣及其鄰近海域P3斷面溫度、鹽度、硝酸鹽濃度、磷酸鹽濃度、CO2(aq)、POC和δ13CPOC的分布Fig.2 Sectional distributions of temperature,salinity,nitrate,phosphate,CO2(aq),POC andδ13CPOCat Section P3 in the Prydz Bay and its adjacent areas

    3.5 δ13CPOC

    普里茲灣及其鄰近海域表層水δ13CPOC的變化范圍介于-28.6‰~-26.2‰,平均值為-27.4‰,最低值出現(xiàn)在弗拉姆淺灘附近的P4-07站,最高值出現(xiàn)在達(dá)恩利角附近的P3-09站。研究海域表層水的δ13CPOC明顯低于中低緯度海域的報(bào)道值,與高緯度海域的報(bào)道值相符合(表1),也與此前報(bào)道的表層水δ13CPOC隨緯度增加而降低的趨勢[5,7,9]相一致。對比普里茲灣表層海水δ13CPOC的數(shù)據(jù),本研究結(jié)果比中國第22次南極考察報(bào)道的普里茲灣表層水體δ13CPOC(-27.4‰~-19.0‰)稍低[19],可能的原因在于本研究航次未采集到埃默里冰架前沿的樣品,而海冰融化所釋放冰藻的δ13CPOC一般較高(-15‰~-8‰)[32]。本研究結(jié)果與同一航次表層δ13CPOC的報(bào)道值(變化范圍為-29.7‰~-26.3‰,平均值為-28.01‰)[20]相一致。

    圖3 普里茲灣及其鄰近海域P5斷面溫度、鹽度、硝酸鹽濃度、磷酸鹽濃度、CO2(aq)、POC和δ13CPOC的分布Fig.3 Sectional distributions of temperature,salinity,nitrate,phosphate,CO2(aq),POC andδ13CPOCat Section P5 in the Prydz Bay and its adjacent areas

    表1 海洋懸浮顆粒有機(jī)物的δ13CPOC值Tab.1 Comparisons inδ13CPOCvalues of suspended particulate organic matter from different marine environments

    普里茲灣及其鄰近海域δ13CPOC的空間分布顯示,灣內(nèi)陸架區(qū)的δ13CPOC總體高于灣外陸坡及深海區(qū),并且無論是陸架區(qū),還是陸坡區(qū)和海盆區(qū),δ13CPOC均隨著深度的增加而增加,如陸架區(qū)次表層的δ13CPOC比混合層高約2‰,而陸坡區(qū)和深海區(qū)深層水的δ13CPOC分別比次表層和混合層高約1‰和2‰(圖2、圖3)。δ13CPOC隨深度增加而增加的形成可能與生源顆粒有機(jī)物垂向輸送過程的再礦化作用有關(guān),在有機(jī)物降解過程中,會優(yōu)先降解貧13C的有機(jī)物,導(dǎo)致殘留顆粒物的δ13CPOC逐漸升高[17,33-34]。與此類似,南大洋澳大利亞扇面1 500 m以淺沉降顆粒有機(jī)物δ13CPOC隨著深度的增加而增加[17];亞北極東北太平洋也發(fā)現(xiàn)真光層以深懸浮顆粒有機(jī)物的δ13CPOC隨深度的增加而顯著增加[33]。

    4 討論

    4.1 普里茲灣混合層δ13CPOC變化的主控因素

    與中深層相比,海洋真光層中的δ13CPOC變化一般較大,其影響因素主要包括溫度、p H值、光強(qiáng)、生物利用無機(jī)碳形態(tài)、浮游植物種群、生長速率等。各因素對δ13CPOC的影響情況歸納如下:(1)溫度:研究表明,海水溫度每升高1℃,δ13CPOC增加0.3‰~0.5‰[40—41]。(2)海水p H值:海水p H值主要通過影響碳酸鹽體系無機(jī)碳的形態(tài)及其含量,從而影響浮游植物固碳底物的含量及其同位素分餾[42]。(3)光照:光強(qiáng)由于會影響初級生產(chǎn)過程,從而對δ13CPOC產(chǎn)生影響。在高緯度海域,混合層深度可作為光照強(qiáng)度變化的指標(biāo),一般而言,混合層越深,水柱平均光強(qiáng)降低,初級生產(chǎn)力下降,導(dǎo)致δ13CPOC較低;相反,混合層變淺時(shí),水柱平均光強(qiáng)增強(qiáng),初級生產(chǎn)力提高,光合作用導(dǎo)致反應(yīng)底物濃度降低,δ13CPOC增加[19,43]。(4)生物利用無機(jī)碳形態(tài):浮游植物利用的無機(jī)碳形態(tài)不同,其產(chǎn)生的δ13CPOC會產(chǎn)生差別。如果浮游植物吸收利用的是HCO-3和CO2-3,則其產(chǎn)生的δ13CPOC要高于利用CO2(aq)產(chǎn)生的δ13CPOC,這與不同固碳過程的同位素分餾因子不同有關(guān),浮游植物通過C3過程Rubisco酶吸收利用CO2(aq)的碳同位素分餾因子要大于C4過程β-羧化酶利用HCO-3和CO2-3的同位素分餾因子[44]。(5)浮游植物種群:浮游植物種群結(jié)構(gòu)的變化,也可能導(dǎo)致其產(chǎn)生的δ13CPOC不同,有研究表明,硅藻產(chǎn)生的δ13CPOC要高于甲藻[45]。(6)浮游植物生長速率:浮游植物生長速率(μ)越快,所產(chǎn)生的δ13CPOC越高[40,46]。另外,有研究表明,浮游植物脂類含量的不同也會影響δ13CPOC,脂類有機(jī)物的δ13C值一般低于氨基酸和蛋白質(zhì)的δ13C,而溫度越低,浮游植物體內(nèi)的脂類物質(zhì)含量往往越高,由此導(dǎo)致δ13CPOC的降低[17]。

    溫度、p H值、光強(qiáng)、生物利用無機(jī)碳形態(tài)、浮游植物種群、生長速率等因素引起的δ13CPOC變化本質(zhì)上都是通過影響浮游植物所吸收的無機(jī)碳含量及吸收過程的同位素分餾來實(shí)現(xiàn)的。對于海洋以C3光合作用過程為主的海洋浮游植物而言,海水中CO2(aq)的變化將對光合作用產(chǎn)生的δ13CPOC起決定性作用。已有研究發(fā)現(xiàn),很多海域的δ13CPOC與[CO2(aq)]之間具有顯著的負(fù)相關(guān)關(guān)系,如印度洋[4]、大西洋[5]、南大洋部分海域[6,8—9]等。對全球海洋表層δ13CPOC進(jìn)行的模型研究也表明,全球海洋表層水δ13CPOC的變化主要受控于[CO2(aq)]的變化[46]。δ13CPOC與[CO2(aq)]負(fù)相關(guān)性的存在主要取決于溫度或者生物活動調(diào)控的[CO2(aq)],由赤道向兩極方向,溫度逐漸降低,CO2(aq)增加,浮游植物碳吸收過程中的同位素分餾作用加強(qiáng),導(dǎo)致δ13CPOC降低[7]。

    普里茲灣及其鄰近海域遠(yuǎn)離陸地,且受埃默里冰架的隔離作用,陸源輸入很小,POC主要來自浮游植物光合作用產(chǎn)生,因而研究海域δ13CPOC的變化受陸源輸入或人類活動的影響很?。?7—29]。本研究結(jié)果顯示,研究區(qū)域混合層(50 m以淺)的δ13CPOC和[CO2(aq)]之間具有顯著的負(fù)相關(guān)關(guān)系:δ13CPOC=-0.23[CO2(aq)]-23.38(n=78,r2=0.22,p<0.000 1)(圖4),說明浮游植物同化吸收CO2過程的同位素分餾是影響普里茲灣及其鄰近海域混合層δ13CPOC變化的主要調(diào)控因素。另外,普里茲灣及其鄰近海域混合層POC含量與溶解CO2濃度之間具有顯著的負(fù)相關(guān)關(guān)系(r2=0.36,p<0.000 1)(圖5),并且溶解CO2與硝酸鹽(r2=0.68,p<0.000 1)和活性磷酸鹽(r2=0.79,p<0.000 1)之間均呈現(xiàn)顯著的正相關(guān)關(guān)系(圖6),這些都可佐證浮游生物活動影響到溶解CO2的含量,進(jìn)而通過同化吸收過程的同位素分餾影響普里茲灣及其鄰近海域δ13CPOC的變化。

    4.2 浮游植物吸收無機(jī)碳過程的碳同位素分餾因子

    圖4 普里茲灣及其鄰近海域混合層中δ13CPOC和[CO2(aq)]的關(guān)系Fig.4 Relationship betweenδ13CPOCand[CO2(aq)]in the mixed layer in the Prydz Bay and its adjacent areas

    圖5 普里茲灣及其鄰近海域混合層中[CO2(aq)]和POC的關(guān)系Fig.5 Concentrations of CO2(aq)verse particulate organic carbon(POC)in the mixed layer in the Prydz Bay and its adjacent areas

    圖6 普里茲灣及其鄰近海域混合層中[CO2(aq)]和硝酸鹽和活性磷酸鹽濃度的關(guān)系Fig.6 Correlation between[CO2(aq)]and nitrate or phosphate concentration in the mixed layer of the Prydz Bay and its adjacent areas

    盡管不同海域的δ13CPOC往往與[CO2(aq)]存在負(fù)相關(guān)關(guān)系,但不同海域二者之間的線性擬合方程往往不同,例如,南大西洋δ13CPOC與[CO2(aq)]的擬合關(guān)系為:δ13CPOC=-0.8[CO2(aq)]-12.6[16];德雷克海峽為:δ13CPOC=-0.9[CO2(aq)]-9.4[38];東北大西洋為:δ13CPOC=-1.5[CO2(aq)]-2.1[5]。不同海域δ13CPOC與[CO2(aq)]擬合方程的不同與浮游植物固碳過程中碳同位素分餾因子的變化有關(guān),生物作用改變δ13CPOC主要是通過改變碳的同位素分餾程度即分餾因子εp來進(jìn)行的。當(dāng)無機(jī)碳只是來自浮游植物通過被動擴(kuò)散所吸收的CO2(aq),考慮到浮游植物種群和生長速率,δ13CPOC與[CO2(aq)]存在如下的關(guān)系[35,37,46]:

    δ13CPOC=δ13CCO2(aq)-εp+a·μ/[CO2(aq)],(1)式中,δ13CCO2(aq)是反應(yīng)物CO2(aq)的碳同位素組成;εp是光合作用固碳過程的碳同位素分餾因子;μ代表浮游植物的生長速率,a是反映不同種群的常數(shù)。由公式可以看出,δ13CPOC受εp、μ和[CO2(aq)]的共同影響。

    本研究獲得的普里茲灣及其鄰近海域混合層δ13CPOC和1/[CO2(aq)]的線性擬合方程為:δ13CPOC=(-31.4±0.9)+69.1/[CO2(aq)](r2=0.19,p<0.000 1)(圖7)。擬合方程的截距(-31.4±0.9)即代表了(δ13CCO2(aq)-εp),而南大洋的δ13CCO2(aq)約為-8‰[47],因此,普里茲灣及其鄰近海域浮游植物同化吸收無機(jī)碳過程中的碳同位素分餾因子εp為23.4‰±0.9‰,這與文獻(xiàn)報(bào)道的浮游植物C3過程Rubsico酶吸收CO2(aq)的分餾因子(25‰~30‰)[43]一致,也再次證實(shí)浮游植物吸收溶解CO2過程是影響普里茲灣及其鄰近海域δ13CPOC變化的主要因素。

    浮游植物種群組成和生長速率的空間變化可能是影響普里茲灣及其鄰近海域δ13CPOC和[CO2(aq)]擬合相關(guān)性的因素之一。普里茲灣灣內(nèi)陸架區(qū)浮游植物優(yōu)勢種為羽紋硅藻(Nitzschiafragilariopsis)、裸鞭毛藻和個(gè)體較大的異養(yǎng)甲藻,水華發(fā)生時(shí)的優(yōu)勢種為羽紋硅藻和大的異養(yǎng)甲藻,而灣外海盆區(qū)的浮游植物優(yōu)勢種為個(gè)體較大的環(huán)紋硅藻和個(gè)體較小的甲藻,羽紋硅藻的δ13C值較高,而裸鞭毛藻和環(huán)紋硅藻的δ13C較低[43]。另外,冰融水份額也會影響δ13CPOC和[CO2(aq)]的擬合相關(guān),對普里茲灣及其鄰近海域的研究發(fā)現(xiàn),表層海水δ13CPOC與冰融水份額之間具有顯著的正相關(guān)關(guān)系,因此,海冰融化會通過影響浮游植物生長速率等影響到光合作用可利用CO2(aq)的變化,進(jìn)而導(dǎo)致δ13CPOC的空間變化[19]。除浮游植物種類和生長速率的影響外,浮游植物吸收無機(jī)碳路徑的不同也會導(dǎo)致δ13CPOC的變化,進(jìn)而影響δ13CPOC和[CO2(aq)]的相關(guān)性,目前已證實(shí)硅藻中存在的β-羧化酶可通過C4路徑吸收海水中的HCO3-[48]。

    5 結(jié)論

    圖7 普里茲灣及其鄰近海域混合層中δ13CPOC和1/[CO2(aq)]的關(guān)系Fig.7 Relationship betweenδ13CPOCand 1/[CO2(aq)]in the mixed layer of the Prydz Bay and its adjacent areas

    普里茲灣及其鄰近海域混合層中的POC含量和δ13CPOC總體表現(xiàn)出灣內(nèi)大于灣外的特征,可歸因于浮游植物同化吸收CO2過程的同位素分餾作用。根據(jù)δ13CPOC和1/[CO2(aq)]的線性擬合關(guān)系,計(jì)算出普里茲灣及其鄰近海域浮游植物同化吸收無機(jī)碳過程中的碳同位素分餾因子為23.4‰。δ13CPOC垂向分布上呈現(xiàn)隨深度增加而增大的特征,可歸因于生源顆粒有機(jī)物垂向輸送過程中再礦化所致的同位素分餾作用。

    致謝:感謝國家海洋局第二海洋研究所韓正兵、潘建明研究員提供營養(yǎng)鹽的數(shù)據(jù)。

    [1] Goni M A,Monacci N,Gisewhite R,et al.Distribution and sources of particulate organic matter in the water column and sediments of the Fly River Delta,Gulf of Papua(Papua New Guinea)[J].Estuarine,Coastal and Shelf Science,2006,69(1/2):225-245.

    [2] Guerra R,Pistocchi R,Vanucci S.Dynamics and sources of organic carbon in suspended particulate matter and sediments in Pialassa Baiona lagoon(NW Adriatic Sea,Italy)[J].Estuarine,Coastal and Shelf Science,2013,135:24-32.

    [3] Wu Y,Zhang J,Liu S M,et al.Sources and distribution of carbon within the Yangtze River system[J].Estuarine,Coastal and Shelf Science,2007,71(1/2):13-25.

    [4] Francois R,Altabet M A,Goericke R,et al.Changes in theδ13Cof surface water particulate organic matter across the subtropical convergencein the SW Indian Ocean[J].Global Biogeochemical Cycles,1993,7(3):627-644.

    [5] Rau G H,Takahashi T,Des Marais D J,et al.The relationship betweenδ13Cof organic matter and[CO2(aq)]in ocean surface water:Data from a JGOFS site in the northeast Atlantic Ocean and a model[J].Geochimica et Cosmochimica Acta,1992,56(3):1413-1419.

    [6] Fischer G.Stable carbon isotope ratios of plankton carbon and sinking organic matter from the Atlantic sector of the Southern Ocean[J].Marine Chemistry,1991,35(1/4):581-596.

    [7] Kennedy H,Robertson J.Variations in the isotopic composition of particulate organic carbon in surface waters along an 88°W transect from 67°S to 54°S[J].Deep-Sea Research PartⅡ,1995,42(4/5):1109-1122.

    [8] Lourey M J,Trull T W,Tilbrook B.Sensitivity ofδ13Cof Southern Ocean suspended and sinking organic matter to temperature,nutrient utilization,and atmospheric CO2[J].Deep-Sea Research PartⅠ:Oceanographic Research Papers,2004,51(2):281-305.

    [9] Rau G H,Sullivan C W,Gordon I L.δ13Candδ15N variations in Weddell Sea particulate organic matter[J].Marine Chemistry,1991,35(1/4):355-369.

    [10] Bentaleb I,F(xiàn)ontugne M,Descolas-Gros C,et al.Organic carbon isotopic composition of phytoplankton and sea-surface pCO2reconstructions in the Southern Indian Ocean during the last 50,000 yr[J].Organic Geochemistry,1996,24(4):399-410.

    [11] Freeman K H,Hayes J M.Fractionation of carbon isotopes by phytoplankton and estimates of ancient CO2levels[J].Global Biogeochemical Cycles,1992,6(2):185-198.

    [12] Fukuchi M.Phytoplankton chlorophyll stocks in the Antarctic Ocean[J].Journal of the Oceanographical Society of Japan,1980,36(1):73-84.

    [13] Arrigo K R,Van Dijken G,Long M.Coastal Southern Ocean:a strong anthropogenic CO2sink[J].Geophysical Research Letters,2008,35(21),doi:10.1029/2008GL035624.

    [14] Cozzi S,Cantoni C.Stable isotope(δ13C andδ15N)composition of particulate organic matter,nutrients and dissolved organic matter during spring ice retreat at Terra Nova Bay[J].Antarctic Science,2010,23(1):43-56.

    [15] Dehairs F,Kopczynska E,Nielsen P,et al.δ13C of Southern Ocean suspended organic matter during spring and early summer:regional and temporal variability[J].Deep-Sea Research PartⅡ,1997,44(1/2):129-142.

    [16] Rau G H,Takahashi T,Des Marais D J.Latitudinal variations in planktonδ13C:implications for CO2and productivity in past oceans[J].Nature,1989,341(6242):516-518.

    [17] Wada E,Terazaki M,Kabaya Y,et al.15N and13C abundances in the Antartic Ocean with emphasis on the biogeochemical structure of the food web[J].Deep-Sea Research Part A.Oceanographic Research Papers,1987,34(5/6):829-841.

    [18] Popp B N,Trull T,Kenig F,et al.Controls on the carbon isotopic composition of southern ocean phytoplankton[J].Global Biogeochemical Cycles,1999,13(4):827-843.

    [19] Zhang R,Zheng M F,Chen M,et al.An isotopic perspective on the correlation of surface ocean carbon dynamics and sea ice melting in Prydz Bay(Antarctica)during austral summer[J].Deep-Sea Research PartⅠ:Oceanographic Research Papers,2014,83:24-33.

    [20] 尹希杰,李云海,喬磊,等.南極普里茲灣海域夏季表層水體顆粒有機(jī)碳及其同位素分布特征[J].極地研究,2014,26(1):159-166.

    Yin Xijie,Li Yunhai,Qiao Lei,et al.Distribution of particulate organic carbon(POC)andδ13CPOCin surface waters in summer in the Prydz Bay,Antarctica[J].Chinese Journal of Polar Research,2014,26(1):159-166.

    [21] 鄭少軍,史久新.南極普里茲灣鄰近海域海冰生消發(fā)展特征分析[J].中國海洋大學(xué)學(xué)報(bào),2011,41(7/8):9-16.

    Zheng Shaojun,Shi Jiuxin.The characteristic of sea ice growth and melt in the Prydz Bay region,Antarctica[J].Periodical of Ocean University of China,2011,41(7/8):9-16.

    [22] Lewis E,Wallace D W R.Program Developed for CO2System Calculations.ORNL/CDIAC-105[Z].Oak Ridge,Tennessee:Carbon Dioxide Information Analysis Center,Oak Ridge National Laboratory,U.S.Department of Energy,1998.

    [23] 張凡.南極普里茲灣碳酸鹽體系與碳通量研究[D].廈門:國家海洋局第三海洋研究所,2013.

    Zhang Fan.Carbonate system and air-sea CO2fluxes in Prydz Bay,the Antarctica[D].Xiamen:The Third Institute of Oceanography,State Oceanic Administration,2013.

    [24] Mehrbach C,Culberson C H,Hawley J E,et al.Measurement of the apparent dissociation constants of carbonic acid in seawater at atmospheric pressure[J].Limnology and Oceanography,1973,18(6):897-907.

    [25] Dickson A G,Millero F J.A comparison of the equilibrium constants for the dissociation of carbonic acid in seawater media[J].Deep-Sea Research Part A:Oceanographic Research Papers,1987,34(10):1733-1743.

    [26] 高郭平,董兆乾,侍茂崇,等.南極普里茲灣關(guān)鍵物理海洋學(xué)問題研究進(jìn)展及未來趨勢[J].上海海洋大學(xué)學(xué)報(bào),2013,22(2):313-320.

    Gao Guoping,Dong Zhaoqian,Shi Maochong,et al.Advances of physical oceanographic study on Prydz Bay and adjacent region,Antarctica[J].Journal of Shanghai Ocean University,2013,22(2):313-320.

    [27] 扈傳昱,張海生,潘建明,等.夏季南極普里茲灣碳的生物地球化學(xué)循環(huán)Ⅱ:POC的分布特征[J].極地研究,2001,13(3):195-204.

    Hu Chuanyu,Zhang Haisheng,Pan Jianming,et al.The biogeochemistry of carbon cycle in summer of the Prydz Bay,AntarcticaⅡ:characteristics of POC distribution[J].Chinese Journal of Polar Research,2001,13(3):195-204.

    [28] 于培松,扈傳昱,朱根海,等.南極普里茲灣海域顆粒有機(jī)碳的變化特征[J].海洋學(xué)報(bào),2011,33(4):181-186.

    Yu Peisong,Hu Chuanyu,Zhu Genhai,et al.Characteristics of particulate organic carbon in the Prydz Bay of Antarctica[J].Haiyang Xuebao,2011,33(4):181-186.

    [29] 韓正兵,扈傳昱,薛斌,等.2007/2008年和2008/2009年夏季南大洋以及普里茲灣POC的分布與變化[J].極地研究,2011,23(1):11-18.

    Han Zhengbin,Hu Chuanyu,Xue Bin,et al.Particulate organic carbon in the surface water of South Ocean and Prydz Bay during the austral summer of 2007/2008 and 2008/2009[J].Chinese Journal of Polar Research,2011,23(1):11-18.

    [30] 邱雨生,黃奕普,劉廣山,等.南極普里茲灣及鄰近海域初級生產(chǎn)力的時(shí)空變異[J].廈門大學(xué)學(xué)報(bào)(自然科學(xué)版),2004,43(5):676-681.

    Qiu Yusheng,Huang Yipu,Liu Guangshan,et al.Spatial and temporal variations of primary productivity in Prydz Bay and its adjacent sea area,Antarctica[J].Journal of Xiamen University(Natural Science),2004,43(5):676-681.

    [31] 劉子琳,蔡昱明,寧修仁,等.1999/2000年夏季南極普里茲灣及灣口區(qū)葉綠素a和初級生產(chǎn)力[J].極地研究,2001,13(1):1-12.

    Liu Zilin,Cai Yuming,Ning Xiuren,et al.Primary productivity and standing stock of phytoplankton in the Prydz Bay and the adjacent northern sea area during the austral summer of 1999/2000[J].Chinese Journal of Polar Research,2001,13(1):1-12.

    [32] Gibson J A E,Trull T W,Nichols P D.Sedimentation of13C-rich organic matter from Antarctic sea-ice algae:A potential indicator of past sea-ice extent[J].Geology,1999,27:331-334.

    [33] Wu Jinping,Calvert S E,Wong C S,et al.Carbon and nitrogen isotopic composition of sedimenting particulate material at station papa in the subarctic northeast Pacific[J].Deep-Sea Research PartⅡ:Topical Studies in Oceanography,1999,46(11/12):2793-2832.

    [34] Guo Laodong,Tanaka T,Wang Deli,et al.Distributions,speciation and stableisotope composition of organic matter in the southeastern Bering Sea[J].Marine Chemistry,2004,91(1/4):211-226.

    [35] Laws E A,Popp B N,Bidigare R R,et al.Dependence of phytoplankton carbon isotopic composition on growth rate and[CO2]aq:theoretical considerations and experimental results[J].Geochimica et Cosmochimica Acta,1995,59(6):1131-1138.

    [36] Liu K-K,Kao S-J,Hu H-C,et al.Carbon isotopic composition of suspended and sinking particulate organic matter in the northern South China Seafrom production to deposition[J].Deep-Sea Research PartⅡ:Topical Studies in Oceanography,2007,54(14/15):1504-1527.

    [37] Zhang Run,Chen Min,Guo Laodong,et al.Variations in the isotopic composition of particulate organic carbon and their relation with carbon dynamics in the western Arctic Ocean[J].Deep-Sea Research PartⅡ:Topical Studies in Oceanography,2012,81-84:72-78.

    [38] Rau G H,Takahashi T,Des Marais D J,et al.Particulate organic matterδ13C variations across the Drake Passage[J].Journal of Geophysical Research:Oceans,1991,96(C8):15131-15135.

    [39] Kopczyńska E E,Goeyens L,Semeneh M,et al.Phytoplankton composition and cell carbon distribution in Prydz Bay,Antarctica:relation to organic particulate matter and itsδ13C values[J].Journal of Plankton Research,1995,17(4):685-707.

    [40] Burkhardt S,Riebesell U,Zondervan I.Effects of growth rate,CO2concentration,and cell size on the stable carbon isotope fractionation in marine phytoplankton[J].Geochimica et Cosmochimica Acta,1999,63(22):3729-3741.

    [41] Hinga K R,Arthur M A,Pilson M E Q,et al.Carbon isotope fractionation by marine phytoplankton in culture:the effects of CO2concentration,p H,temperature,and species[J].Global Biogeochemical Cycles,1994,8(1):91-102.

    [42] Wefer G,Killingley J S.Carbon isotopes in organic matter from a benthic algaHalimedaincrassata(Bermuda):effects of light intensity[J].Chemical Geology,1986,59:321-326.

    [43] Drake B L.Using models of carbon isotope fractionation during photosynthesis to understand the natural fractionation ratio[J].Radiocarbon,2014,56(1):29-38.

    [44] Falkowski P G.Species variability in the fractionation of13C and12C by marine phytoplankton[J].Journal of Plankton Research,1991,13(S1):21-28.

    [45] Laws E A,Bidigare R R,Popp B N.Effect of growth rate and CO2concentration on carbon isotopic fractionation by the marine diatomPhaeodactylumtricornutum[J].Limnology and Oceanography,1997,42(7):1552-1560.

    [46] Hofmann M,Wolf-Gladrow D A,Takahashi T,et al.Stable carbon isotope distribution of particulate organic matter in the ocean:a model study[J].Marine Chemistry,2000,72(2/4):131-150.

    [47] Henley S F,Annett A L,Ganeshram R S,et al.Factors influencing the stable carbon isotopic composition of suspended and sinking organic matter in the coastal Antarctic sea ice environment[J].Biogeosciences,2012,9(3):1137-1157.

    [48] Tortell P D,Payne C,Gueguen C,et al.Inorganic carbon uptake by Southern Ocean phytoplankton[J].Limnology and Oceanography,2008,53(4):1266-1278.

    Stable carbon isotopic composition of suspended particulate organic matter in the Prydz Bay and its adjacent areas

    Ren Chunyan1,Chen Min1,2,Gao Zhongyong3,Guo Laodong1,4,Jia Renming1,Liu Xiao1,Qiu Yusheng1,Zheng Minfang1

    (1.College of Ocean and Earth Sciences,Xiamen University,Xiamen 361102,China;2.State Key Laboratory of Marine Environmental Science,Xiamen University,Xiamen 361102,China;3.Key Lab of Global Change and Marine-Atmospheric Chemistry,Third Institute of Oceanography,State Oceanic Administration,Xiamen 361005,China;4.School of Freshwater Sciences,University of Wisconsin-Milwaukee,WI 53204,USA)

    To better understand carbon and particle dynamics in the Prydz Bay and its adjacent areas,stable carbon isotopic composition in suspended particulate organic matter(POM)along with aqueous CO2,nutrients and hydrographic data were measured to examine the spatial and vertical distributions ofδ13CPOCand their controlling factors during the 29thChinese Antarctic Research Expedition.Furthermore,carbon isotope fractionation factor during the assimilation of dissolved inorganic carbon(DIC)by phytoplankton in the mixed layer was estimated.Our resultsshowed thatδ13CPOCranged from-28.5‰to-21.1‰with an average of-24.6‰,showing an increase trend from the inner bay to open waters outside the Bay.Carbon isotope fractionation factor during phytoplankton DIC assimilation in the mixed layer was estimated to be 23.4‰based on the relationship betweenδ13CPOCand 1/[CO2(aq)],and was the major controlling factor forδ13CPOCvalues in the mixed layer in the Prydz Bay and its adjacent areas.Values ofδ13CPOCincreased in general with increasing depth in the study areas,ascribing to the isotope fractionation during POM remineralization.

    suspended particulate organic matter;δ13CPOC;isotope fractionation factor;Prydz Bay;Antarctic

    X145

    A

    0253-4193(2015)12-0074-11

    任春燕,陳敏,高眾勇,等.南極普里茲灣及其鄰近海域懸浮顆粒有機(jī)物的碳同位素組成及其影響因素[J].海洋學(xué)報(bào),2015,37(12):74—84,

    10.3969/j.issn.0253-4193.2015.12.008

    Ren Chunyan,Chen Min,Gao Zhongyong,et al.Stable carbon isotopic composition of suspended particulate organic matter in the Prydz Bay and its adjacent areas[J].Haiyang Xuebao,2015,37(12):74—84,doi:10.3969/j.issn.0253-4193.2015.12.008

    2015-09-30;

    2015-10-18。

    南北極環(huán)境綜合考察與評估專項(xiàng)(CHINARE01-04-03,04-01-06,04-04,04-03-05);國家自然科學(xué)基金杰出青年基金項(xiàng)目(41125020)。

    任春燕(1985—),女,天津市人,博士研究生,從事同位素海洋化學(xué)研究。E-mail:wounderful007@163.com

    *通信作者:陳敏(1970—),男,廣東省興寧市人,教授,博士,從事同位素海洋化學(xué)研究。E-mail:mchen@xmu.edu.cn

    猜你喜歡
    普里南極表層
    半潛式平臺表層卡套管處理與認(rèn)識
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對磷的吸收及釋放研究進(jìn)展
    我去南極能住哪兒
    來到南極要補(bǔ)課
    南極大逃亡
    難忘的南極之旅
    聯(lián)合國助理秘書長兼婦女署副執(zhí)行主任普里來訪
    氬弧熔覆原位合成Ti(C,N)-WC增強(qiáng)鎳基表層復(fù)合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    超聲波光整強(qiáng)化40Cr表層顯微硬度研究
    国产成人精品久久二区二区91 | 日韩欧美精品免费久久| 韩国高清视频一区二区三区| 中文乱码字字幕精品一区二区三区| 免费在线观看黄色视频的| 婷婷色综合www| 天美传媒精品一区二区| 久久久久人妻精品一区果冻| 国产人伦9x9x在线观看 | kizo精华| 国产成人一区二区在线| 精品酒店卫生间| 精品国产乱码久久久久久男人| 午夜福利在线观看免费完整高清在| av网站免费在线观看视频| 最新的欧美精品一区二区| 国产精品一国产av| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久精品电影小说| av线在线观看网站| 国产男女超爽视频在线观看| 女人精品久久久久毛片| av不卡在线播放| 十八禁高潮呻吟视频| 国产一区二区激情短视频 | 伊人久久大香线蕉亚洲五| 人人妻人人爽人人添夜夜欢视频| 国产精品偷伦视频观看了| 丰满乱子伦码专区| 少妇 在线观看| 成年动漫av网址| 大话2 男鬼变身卡| 色哟哟·www| 国产精品熟女久久久久浪| 国产亚洲精品第一综合不卡| 日本-黄色视频高清免费观看| 亚洲伊人久久精品综合| 国产精品久久久久久av不卡| 女人高潮潮喷娇喘18禁视频| 亚洲中文av在线| 婷婷色综合大香蕉| 国产精品久久久久久av不卡| 久久精品国产亚洲av天美| 久久久久国产网址| av福利片在线| 精品卡一卡二卡四卡免费| 纯流量卡能插随身wifi吗| 中文字幕人妻丝袜一区二区 | 国产日韩一区二区三区精品不卡| 免费看av在线观看网站| 99re6热这里在线精品视频| 最新的欧美精品一区二区| 男女高潮啪啪啪动态图| 免费看av在线观看网站| 国产在线视频一区二区| 亚洲欧美成人综合另类久久久| 亚洲精品av麻豆狂野| 国产精品久久久久久久久免| 国产成人精品在线电影| 精品久久蜜臀av无| 久久精品夜色国产| 91精品三级在线观看| 亚洲美女黄色视频免费看| 成人漫画全彩无遮挡| 丰满饥渴人妻一区二区三| 欧美变态另类bdsm刘玥| 亚洲熟女精品中文字幕| 精品国产国语对白av| www.自偷自拍.com| 黄频高清免费视频| 制服人妻中文乱码| 啦啦啦啦在线视频资源| 国产白丝娇喘喷水9色精品| 免费播放大片免费观看视频在线观看| 综合色丁香网| 久久精品国产a三级三级三级| 久久久国产精品麻豆| 国产一区二区在线观看av| 咕卡用的链子| 午夜福利乱码中文字幕| 青草久久国产| 王馨瑶露胸无遮挡在线观看| 青春草视频在线免费观看| 久久久久久久久久久久大奶| 两个人看的免费小视频| 久久人人97超碰香蕉20202| 国产精品.久久久| 久久久国产一区二区| 日本猛色少妇xxxxx猛交久久| 一区二区三区乱码不卡18| 波多野结衣av一区二区av| 欧美日韩综合久久久久久| 18禁裸乳无遮挡动漫免费视频| 久久精品国产综合久久久| 精品一区二区三区四区五区乱码 | 一区二区三区精品91| 国产成人免费观看mmmm| 亚洲四区av| 伊人久久国产一区二区| 男女免费视频国产| 青春草亚洲视频在线观看| 久久久精品区二区三区| 精品酒店卫生间| 亚洲综合色惰| 最近2019中文字幕mv第一页| 黄片小视频在线播放| 香蕉国产在线看| 中文字幕av电影在线播放| 亚洲美女搞黄在线观看| 在线天堂最新版资源| 777米奇影视久久| 下体分泌物呈黄色| 久久久久国产网址| xxx大片免费视频| 亚洲国产色片| 岛国毛片在线播放| 久久久国产一区二区| 侵犯人妻中文字幕一二三四区| 国产又色又爽无遮挡免| 国产国语露脸激情在线看| 亚洲国产最新在线播放| 国产成人精品一,二区| 电影成人av| 国产一区二区 视频在线| 如何舔出高潮| 狠狠婷婷综合久久久久久88av| 国产精品免费大片| 久久久久久人人人人人| 大话2 男鬼变身卡| 久久精品久久精品一区二区三区| 欧美av亚洲av综合av国产av | 国产精品免费大片| 中文字幕制服av| 美女脱内裤让男人舔精品视频| 国产亚洲欧美精品永久| 蜜桃国产av成人99| 欧美最新免费一区二区三区| 久久久亚洲精品成人影院| 免费播放大片免费观看视频在线观看| 国产一区二区 视频在线| 视频在线观看一区二区三区| 亚洲经典国产精华液单| 伊人久久大香线蕉亚洲五| 男人添女人高潮全过程视频| 中文字幕人妻丝袜制服| 国产熟女午夜一区二区三区| 伦理电影免费视频| 99热全是精品| 国产精品国产av在线观看| 女性生殖器流出的白浆| 亚洲国产精品一区二区三区在线| 亚洲美女黄色视频免费看| 一边亲一边摸免费视频| 精品国产国语对白av| 国产精品蜜桃在线观看| 亚洲精品成人av观看孕妇| 欧美精品亚洲一区二区| 成年女人在线观看亚洲视频| 久久久久久久久久人人人人人人| 欧美亚洲 丝袜 人妻 在线| 亚洲内射少妇av| 黄色视频在线播放观看不卡| 久久这里有精品视频免费| 午夜福利在线观看免费完整高清在| 日韩中文字幕视频在线看片| 黄片播放在线免费| 国产精品国产三级国产专区5o| www.自偷自拍.com| 黄色 视频免费看| 99久国产av精品国产电影| 色视频在线一区二区三区| 在线观看美女被高潮喷水网站| 97精品久久久久久久久久精品| 亚洲少妇的诱惑av| 91精品三级在线观看| 欧美精品一区二区大全| 亚洲精华国产精华液的使用体验| 久久久久久免费高清国产稀缺| 中文字幕av电影在线播放| 七月丁香在线播放| 天堂中文最新版在线下载| 日韩电影二区| 亚洲一区中文字幕在线| 色哟哟·www| 在线观看免费日韩欧美大片| 久久精品国产亚洲av涩爱| 久久人人97超碰香蕉20202| 性色avwww在线观看| 成年女人毛片免费观看观看9 | 欧美日韩亚洲国产一区二区在线观看 | 纯流量卡能插随身wifi吗| 好男人视频免费观看在线| 爱豆传媒免费全集在线观看| 80岁老熟妇乱子伦牲交| 日本欧美国产在线视频| 日本猛色少妇xxxxx猛交久久| 国产极品天堂在线| 亚洲精品美女久久av网站| 国产深夜福利视频在线观看| 国语对白做爰xxxⅹ性视频网站| 又粗又硬又长又爽又黄的视频| 电影成人av| 久久精品国产亚洲av涩爱| 久久精品国产a三级三级三级| 亚洲av免费高清在线观看| 国产视频首页在线观看| 婷婷色综合大香蕉| 精品国产露脸久久av麻豆| 国产淫语在线视频| 国产乱来视频区| 日日爽夜夜爽网站| 一级毛片 在线播放| 老鸭窝网址在线观看| 精品少妇黑人巨大在线播放| 韩国精品一区二区三区| 国产精品久久久久久精品电影小说| 两个人看的免费小视频| 亚洲av成人精品一二三区| 国产免费福利视频在线观看| 久久亚洲国产成人精品v| 久久久久久久国产电影| 婷婷色综合www| 久久影院123| 99热网站在线观看| 亚洲av日韩在线播放| 18禁国产床啪视频网站| 亚洲国产精品国产精品| 人妻少妇偷人精品九色| 久久久久久久精品精品| 一本—道久久a久久精品蜜桃钙片| 亚洲av国产av综合av卡| 黄色一级大片看看| 国产野战对白在线观看| 少妇人妻精品综合一区二区| 久久国产亚洲av麻豆专区| av网站在线播放免费| 亚洲熟女精品中文字幕| 国产成人精品久久久久久| 国产精品二区激情视频| 国产野战对白在线观看| 韩国高清视频一区二区三区| 一本久久精品| 黄色怎么调成土黄色| 中文字幕色久视频| 国产黄频视频在线观看| 精品国产露脸久久av麻豆| 1024香蕉在线观看| av线在线观看网站| 麻豆精品久久久久久蜜桃| 日韩免费高清中文字幕av| 青青草视频在线视频观看| 一本色道久久久久久精品综合| 亚洲美女黄色视频免费看| 午夜福利网站1000一区二区三区| 日本黄色日本黄色录像| 一区二区日韩欧美中文字幕| 伊人久久大香线蕉亚洲五| 免费黄色在线免费观看| 精品久久久精品久久久| 人人妻人人添人人爽欧美一区卜| 午夜激情久久久久久久| 少妇被粗大猛烈的视频| 久久久亚洲精品成人影院| av.在线天堂| 免费久久久久久久精品成人欧美视频| 精品少妇一区二区三区视频日本电影 | 成人毛片60女人毛片免费| 欧美亚洲 丝袜 人妻 在线| 美女国产视频在线观看| 丝袜美足系列| av卡一久久| 国产av国产精品国产| 亚洲国产精品一区三区| 免费av中文字幕在线| 爱豆传媒免费全集在线观看| 成人手机av| 久久99一区二区三区| 叶爱在线成人免费视频播放| 国产成人精品福利久久| 中文字幕人妻丝袜一区二区 | 免费观看av网站的网址| 亚洲精品美女久久av网站| 一本色道久久久久久精品综合| 最近的中文字幕免费完整| 亚洲人成网站在线观看播放| 色网站视频免费| 男女啪啪激烈高潮av片| 搡女人真爽免费视频火全软件| 中文精品一卡2卡3卡4更新| 亚洲色图综合在线观看| 久久久久久免费高清国产稀缺| 超碰97精品在线观看| 免费看av在线观看网站| 秋霞在线观看毛片| 精品少妇一区二区三区视频日本电影 | 日本91视频免费播放| 久久久久久久大尺度免费视频| 三级国产精品片| 久久久久国产精品人妻一区二区| 老司机影院成人| 色婷婷久久久亚洲欧美| 超碰成人久久| 超碰成人久久| 2021少妇久久久久久久久久久| 亚洲欧美一区二区三区久久| 国产免费视频播放在线视频| 久久精品久久精品一区二区三区| 美女xxoo啪啪120秒动态图| 国产视频首页在线观看| 男女国产视频网站| 久久ye,这里只有精品| 国产视频首页在线观看| 久久ye,这里只有精品| 男女国产视频网站| 成人影院久久| 精品第一国产精品| 黑人欧美特级aaaaaa片| 中国国产av一级| 两个人看的免费小视频| 青青草视频在线视频观看| 国产综合精华液| 精品少妇久久久久久888优播| 国产精品女同一区二区软件| 欧美日韩视频高清一区二区三区二| 国产成人av激情在线播放| 久久人人爽av亚洲精品天堂| 久久人人爽av亚洲精品天堂| 我的亚洲天堂| 日日摸夜夜添夜夜爱| 国产1区2区3区精品| 国产免费一区二区三区四区乱码| 亚洲国产欧美网| 精品少妇内射三级| 精品国产国语对白av| 亚洲综合精品二区| 国语对白做爰xxxⅹ性视频网站| 婷婷色综合www| 丝袜脚勾引网站| 最近的中文字幕免费完整| 三级国产精品片| 女的被弄到高潮叫床怎么办| 伊人久久大香线蕉亚洲五| 十八禁高潮呻吟视频| 亚洲第一av免费看| 久久久精品免费免费高清| 久久亚洲国产成人精品v| 亚洲精品美女久久久久99蜜臀 | 日本wwww免费看| 99热网站在线观看| 亚洲一级一片aⅴ在线观看| 精品人妻偷拍中文字幕| 亚洲第一区二区三区不卡| 99国产精品免费福利视频| av有码第一页| 中文字幕另类日韩欧美亚洲嫩草| 成人毛片60女人毛片免费| 最新的欧美精品一区二区| 成人黄色视频免费在线看| 国产精品秋霞免费鲁丝片| 免费不卡的大黄色大毛片视频在线观看| 国产男人的电影天堂91| 久久ye,这里只有精品| 啦啦啦在线免费观看视频4| 一区二区三区乱码不卡18| 黄色视频在线播放观看不卡| 妹子高潮喷水视频| www.自偷自拍.com| 日韩电影二区| 十八禁网站网址无遮挡| 99国产综合亚洲精品| 伊人亚洲综合成人网| av免费观看日本| 亚洲精品乱久久久久久| 国产人伦9x9x在线观看 | 日韩精品免费视频一区二区三区| 国产精品不卡视频一区二区| 女人久久www免费人成看片| 在线看a的网站| 久久99精品国语久久久| 又大又黄又爽视频免费| 大片电影免费在线观看免费| 久久久精品免费免费高清| 人妻一区二区av| 国产精品av久久久久免费| 视频在线观看一区二区三区| 久久精品久久久久久噜噜老黄| 一本色道久久久久久精品综合| 亚洲欧洲精品一区二区精品久久久 | 国产亚洲一区二区精品| 三上悠亚av全集在线观看| 久久狼人影院| 久久鲁丝午夜福利片| 亚洲av日韩在线播放| 一区在线观看完整版| 国产一区亚洲一区在线观看| 亚洲av电影在线观看一区二区三区| 精品国产国语对白av| 麻豆乱淫一区二区| videos熟女内射| 18在线观看网站| 99久久精品国产国产毛片| 免费在线观看完整版高清| 欧美日韩精品成人综合77777| 国产免费视频播放在线视频| 日日啪夜夜爽| 国产成人精品久久久久久| 亚洲欧洲日产国产| 少妇 在线观看| 日韩电影二区| 熟女电影av网| 精品国产乱码久久久久久小说| 一级毛片电影观看| 2021少妇久久久久久久久久久| 欧美 日韩 精品 国产| 一区福利在线观看| 色网站视频免费| 哪个播放器可以免费观看大片| 日韩av免费高清视频| 午夜91福利影院| 老汉色∧v一级毛片| 成人午夜精彩视频在线观看| 熟女电影av网| 国产极品粉嫩免费观看在线| 丰满乱子伦码专区| 中文乱码字字幕精品一区二区三区| 黄色一级大片看看| 久久久国产一区二区| 人妻系列 视频| 在线天堂中文资源库| 美女国产视频在线观看| a级毛片黄视频| 亚洲国产精品999| 国产 一区精品| 精品国产国语对白av| 久久这里有精品视频免费| 男女免费视频国产| 亚洲精品乱久久久久久| 大片免费播放器 马上看| 久久久久久久大尺度免费视频| 国产精品不卡视频一区二区| 亚洲av欧美aⅴ国产| 777米奇影视久久| 亚洲综合精品二区| 免费黄网站久久成人精品| 亚洲成av片中文字幕在线观看 | 丝袜喷水一区| 久久久欧美国产精品| 亚洲精品aⅴ在线观看| 精品一区二区免费观看| 亚洲精品av麻豆狂野| 婷婷色av中文字幕| 日韩在线高清观看一区二区三区| 国产成人精品久久久久久| 一级毛片 在线播放| 国产免费一区二区三区四区乱码| 国产欧美日韩一区二区三区在线| 欧美+日韩+精品| 大香蕉久久成人网| www.熟女人妻精品国产| a 毛片基地| 日韩精品免费视频一区二区三区| 男女免费视频国产| 可以免费在线观看a视频的电影网站 | 伦理电影大哥的女人| 婷婷色综合大香蕉| xxxhd国产人妻xxx| 久久久久久免费高清国产稀缺| 欧美老熟妇乱子伦牲交| 在线观看免费日韩欧美大片| 久久影院123| 国产亚洲一区二区精品| 欧美老熟妇乱子伦牲交| 日产精品乱码卡一卡2卡三| 男女国产视频网站| 秋霞在线观看毛片| 国产精品不卡视频一区二区| 欧美日韩成人在线一区二区| 青春草亚洲视频在线观看| 青春草国产在线视频| 人妻一区二区av| 日本色播在线视频| 久久精品久久精品一区二区三区| 国产一级毛片在线| 国产在线免费精品| 精品国产一区二区久久| 丰满少妇做爰视频| videos熟女内射| 日本-黄色视频高清免费观看| 日韩中文字幕视频在线看片| 亚洲欧美一区二区三区国产| 国产精品久久久久久精品古装| 男女边摸边吃奶| 97精品久久久久久久久久精品| 18在线观看网站| 久久久精品区二区三区| av在线观看视频网站免费| 性色av一级| 日韩熟女老妇一区二区性免费视频| 国产精品成人在线| 久久久久久久久免费视频了| 欧美成人午夜精品| 亚洲精品美女久久久久99蜜臀 | 中国国产av一级| 欧美 日韩 精品 国产| 波野结衣二区三区在线| 爱豆传媒免费全集在线观看| 9热在线视频观看99| 尾随美女入室| 国产欧美日韩综合在线一区二区| 国产成人a∨麻豆精品| 久久久a久久爽久久v久久| 交换朋友夫妻互换小说| 成人18禁高潮啪啪吃奶动态图| 久热久热在线精品观看| 亚洲人成77777在线视频| 爱豆传媒免费全集在线观看| 最近2019中文字幕mv第一页| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 制服诱惑二区| 久久午夜福利片| 日韩视频在线欧美| 国产黄频视频在线观看| 高清不卡的av网站| 一级爰片在线观看| 午夜日本视频在线| 黄色毛片三级朝国网站| 丁香六月天网| 日韩av不卡免费在线播放| 欧美日韩视频高清一区二区三区二| 九九爱精品视频在线观看| 久久精品国产亚洲av高清一级| 欧美精品av麻豆av| www.自偷自拍.com| 边亲边吃奶的免费视频| 人人澡人人妻人| 国产野战对白在线观看| 一级片免费观看大全| 日本爱情动作片www.在线观看| 涩涩av久久男人的天堂| 亚洲一区二区三区欧美精品| 久久 成人 亚洲| 国产xxxxx性猛交| 男人操女人黄网站| 亚洲精品国产av成人精品| 蜜桃在线观看..| 日韩制服骚丝袜av| 免费黄色在线免费观看| 女人久久www免费人成看片| 国产成人精品婷婷| 欧美成人精品欧美一级黄| 欧美xxⅹ黑人| av国产精品久久久久影院| 国产xxxxx性猛交| 亚洲天堂av无毛| 老司机亚洲免费影院| 在线天堂中文资源库| 国产亚洲欧美精品永久| 女的被弄到高潮叫床怎么办| 满18在线观看网站| 免费观看无遮挡的男女| 男人添女人高潮全过程视频| 国产欧美日韩综合在线一区二区| 久久久久精品性色| videosex国产| 在线观看国产h片| av在线老鸭窝| 最黄视频免费看| 亚洲欧美一区二区三区国产| 波野结衣二区三区在线| 在线观看国产h片| av在线老鸭窝| 欧美日韩亚洲国产一区二区在线观看 | 国产精品欧美亚洲77777| videos熟女内射| 99久久中文字幕三级久久日本| 国产野战对白在线观看| 亚洲欧美成人综合另类久久久| 少妇 在线观看| 亚洲欧美精品自产自拍| 在线观看免费高清a一片| 18禁国产床啪视频网站| 啦啦啦在线观看免费高清www| 亚洲美女搞黄在线观看| 中文欧美无线码| 国产精品免费视频内射| 欧美xxⅹ黑人| 午夜91福利影院| 大香蕉久久网| 久久热在线av| 精品一品国产午夜福利视频| 亚洲av.av天堂| 男人舔女人的私密视频| 街头女战士在线观看网站| 下体分泌物呈黄色| 肉色欧美久久久久久久蜜桃| 少妇猛男粗大的猛烈进出视频| 熟女av电影| 色哟哟·www| 国产免费又黄又爽又色| 国产成人91sexporn| 极品少妇高潮喷水抽搐| 国产成人精品婷婷| 国产精品不卡视频一区二区| 成人漫画全彩无遮挡| 91精品三级在线观看| 国产av国产精品国产| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 欧美日韩综合久久久久久| 夜夜骑夜夜射夜夜干| 最近最新中文字幕大全免费视频 | 亚洲av国产av综合av卡| 国产免费一区二区三区四区乱码| 伦理电影大哥的女人| 啦啦啦在线免费观看视频4|