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

    會(huì)仙喀斯特濕地三種典型植物葉片碳同位素(δ13C)特征及其指示意義

    2021-07-19 01:43沈育伊張德楠徐廣平滕秋梅周龍武黃科朝牟芝熠孫英杰
    廣西植物 2021年5期
    關(guān)鍵詞:水葫蘆喀斯特同位素

    沈育伊 張德楠 徐廣平 滕秋梅 周龍武 黃科朝 牟芝熠 孫英杰

    摘 要: ?為探討會(huì)仙喀斯特濕地不同生長環(huán)境下植物葉片碳(C)、氮(N)、磷(P)、穩(wěn)定碳同位素(δ13C)特征及其生態(tài)學(xué)指示意義,該文以挺水植物蘆葦、浮水植物水葫蘆和沉水植物金魚藻為研究對(duì)象,分析了三種典型不同生活型植物葉片的δ13 C特征及種間和微生境的差異,并基于植物碳同位素與碳酸酐酶顯著正相關(guān)的二端元模型,估算了植物利用光合作用固定的碳酸氫根離子(HCO-3)的碳量。結(jié)果表明:(1)三種植物葉片δ13C的變化范圍為-28.47‰~-21.69‰,平均值為-24.83‰,不同生活型植物間δ13C存在差異,金魚藻>水葫蘆>蘆葦。(2)植物δ13C值與葉片C、N和P元素含量間呈顯著正相關(guān),與C/N、C/P和N/P呈顯著負(fù)相關(guān)關(guān)系,與底泥的有機(jī)質(zhì)、速效氮、總氮、速效磷和總磷含量呈顯著正相關(guān)。(3)會(huì)仙喀斯特濕地三種不同生活型植物葉片N/P 平均值為10.34,表現(xiàn)出植物受N、P共同影響的特征。(4)δ13C的變化特征,揭示了三種水生植物可能通過增加磷利用效率來促進(jìn)低水分利用率環(huán)境下的碳的合成,通過提高植物水分利用效率的策略來代償較低的氮素利用效率。(5)蘆葦光合作用固定HCO-3碳量為159.60 t·a-1·km-2,水葫蘆為10.80 t·a-1·km-2,金魚藻為9.24 t·a-1·km-2,平均值為59.88 t·a-1·km-2。會(huì)仙喀斯特濕地植物的不同生活型、光合作用途徑和生長微環(huán)境,是影響葉片δ13C變化的主要因素。

    關(guān)鍵詞: 喀斯特濕地, 水生植物, 營養(yǎng)元素, 穩(wěn)定碳同位素(δ13C), 碳匯

    中圖分類號(hào): ?Q178.1

    文獻(xiàn)標(biāo)識(shí)碼: ?A

    文章編號(hào): ?1000-3142(2021)05-0769-11

    Abstract: ?In order to understand characteristics and significance of carbon (C), nitrogen (N), phosphorus (P) and stable carbon isotope composition (δ13C) in plant leaves of three typical aquatic plants under different growth environments, three life forms of aquatic plants including emerging plant (Phragmites australis), emersion plant (Eichhornia crassipes) and submerged plant (Ceratophyllum demersumin) of the Huixian karst wetland were selected as study objects. The differences of carbon isotope compositions of interspecies and different habitat conditions were studied, and the percentage of HCO-3 uptaken from water by aquatic photosynthesis was also estimated using the double-meta model. The results were as follows: (1) The foliar δ13C of three diffirent plant life forms ranged from -28.47‰~-21.69‰ with an average of -24.83‰. There were significant differences among species, and the sequence of relative value was Ceratophyllum demersum>Eichhornia crassipes>Phragmites australis, and among the three types of plants,Phragmites australis and Ceratophyllum demersum show the lowest and highest δ13 C, respectively. (2) For three types of halophytes,δ13C showed significantly positive correlations with foliar C, N and P. The δ13C was also positively correlated with the contents of organic matter, available nitrogen, total nitrogen, available phosphorus and total phosphorus of sediment. All types of halophytes showed a negative correlation among δ13C, C/N, N/P and C/P. (3) The foliar N/P ratio of halophytes was 10.34, indicating a co-limitation by N and P for halophyte growth in the Huixian karst wetland. (4) These results suggested that three halophytes were adapted to environment through a high water use efficiency at the expense of decreased nitrogen use efficiency, and might improve C sequestration by increasing phosphorus use efficiency in the habitat of low water availability. (5) The value of HCO-3 carbon sequestration by photosynhesis was 159.60 t·a-1·km-2 for Phragmites australis, 10.80 t·a-1·km-2 for Eichhornia crassipes, 9.24 t·a-1·km-2 for Ceratophyllum demersum, with a mean value of 59.88 t·a-1·km-2. Different plant life forms, pathways of carbon fixation in photosynthesis and micro-environment were the important factors on affecting plant foliar δ13C in the Huixian karst wetland.

    Key words: karst wetland, aquatic macrophytes, nutrient elements, stable carbon isotopic δ13C, carbon sink

    碳(C)是植物生長中重要的元素,水生植物的穩(wěn)定碳同位素(δ13C)記錄了植物生長過程中外界環(huán)境的信息變化,可揭示碳的生物地球化學(xué)循環(huán)過程(Van et al.,2018)。穩(wěn)定碳同位素已被廣泛應(yīng)用于研究水生生態(tài)系統(tǒng)演變和營養(yǎng)傳遞規(guī)律等方面(Dijkstra et al.,2016;韓瓊和王寶利,2017;劉金亮等,2020)。會(huì)仙喀斯特濕地與漓江水系構(gòu)成了“桂林之腎”,對(duì)維護(hù)漓江大水系體系——“六江四庫四湖一濕地”的生態(tài)安全意義重大,發(fā)揮著重要的生態(tài)功能,研究濕地水生植物葉片的穩(wěn)定碳同位素特征,對(duì)深入揭示喀斯特濕地水生植物與環(huán)境間的相互關(guān)系具有重要指示意義。

    近年來,利用穩(wěn)定碳同位素技術(shù),開展植物個(gè)體與環(huán)境間關(guān)系的研究逐漸增多,如林清和王紹令(2001)對(duì)沉水植物龍須眼子菜與環(huán)境之間的關(guān)系進(jìn)行了分析,黃亮等(2003)對(duì)長江中游若干湖泊水生植物體內(nèi)碳氮磷及δ13C的分布進(jìn)行研究,章程等(2013)研究了桂林會(huì)仙水生植物δ13C特征并對(duì)固碳量進(jìn)行了估算,陳清等(2015)對(duì)天津沼澤濕地蘆葦葉片δ13C特征及其環(huán)境影響因素進(jìn)行了研究,以上研究主要集中于碳同位素在各濕地中的應(yīng)用及其環(huán)境指示意義。而植物葉片δ13C值在不同環(huán)境條件下會(huì)發(fā)生變化,并可能在植物種內(nèi)和種間呈現(xiàn)差異性(Sprent et al.,1996)。δ13C與氮(N)、磷(P)元素化學(xué)計(jì)量學(xué)的關(guān)系,也可能會(huì)因不同物種和不同區(qū)域環(huán)境而呈現(xiàn)不同(Zhou et al.,2016;Ma et al.,2018)。

    隨著人類活動(dòng)的增加,會(huì)仙喀斯特濕地水體富營養(yǎng)化、濕地面積縮小等問題突出,影響著濕地生態(tài)系統(tǒng)的功能(徐廣平等,2019;滕秋梅等,2021)。一般沉水植物的葉子生長在水面之下,而浮水植物和挺水植物的葉子生長在水面之上,因而不同生活型植物生長的微環(huán)境有所不同。目前會(huì)仙喀斯特濕地不同生活型植物葉片δ13C值指示的植物水分利用效率(water use efficiency,WUE)與氮磷養(yǎng)分利用效率(nitrogen use efficiency,NUE;phosphorus use efficiency,PUE)間的權(quán)衡關(guān)系尚不明確。根據(jù)水生植物的δ13C值,亦可指示水生植物光合作用固定碳酸氫根離子(HCO-3)的比例和固碳量。鑒于此,本研究以會(huì)仙喀斯特濕地挺水植物蘆葦(Phragmites australis)、浮水植物水葫蘆(Eichhornia crassipes)和沉水植物金魚藻(Ceratophyllum demersum)為研究對(duì)象,探討不同生長環(huán)境下典型水生植物葉片δ13C特征與微環(huán)境間的關(guān)系、δ13C的季節(jié)性變化趨勢(shì),及其光合作用固定碳酸氫根離子(HCO-3)的比例,旨在為豐富桂林會(huì)仙喀斯特濕地生態(tài)學(xué)和生態(tài)系統(tǒng)服務(wù)功能等的研究提供基礎(chǔ)數(shù)據(jù),為廣西桂林會(huì)仙喀斯特國家濕地公園保護(hù)和建設(shè)提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    會(huì)仙喀斯特濕地位于廣西桂林市臨桂縣會(huì)仙鎮(zhèn),東至雁山區(qū),西至四塘鄉(xiāng),25°1′30″—25°11′15″ N、110°8′15″—110°18′00″ E,總面積約為120 km2。氣候?qū)賮啛釒Ъ撅L(fēng)氣候,年平均氣溫為16.5~20.5 ℃,年降水量為1 890.4 mm。土壤以紅黃壤和紅壤為主,濕地植物種類較多,且生長茂盛,主要建群種有華克拉莎(Cladium chinense)、長苞香蒲群落(Typha angustata)、金魚藻等(孫英杰等,2018)。從濕地沿岸到湖心區(qū)域分布有挺水植物蘆葦、浮水植物水葫蘆和沉水植物金魚藻三種代表性的優(yōu)勢(shì)水生植物。濕地試驗(yàn)采樣區(qū)水環(huán)境特征見表1。

    2 材料與方法

    2.1 樣地設(shè)置和樣品的采集

    2013年9月,沿水位梯度不同生長環(huán)境,在會(huì)仙喀斯特濕地從河岸附近到湖心深水區(qū)設(shè)置采樣區(qū),按生活型選取代表性的挺水植物蘆葦、浮水植物水葫蘆和沉水植物金魚藻三個(gè)不同植物群落,每個(gè)群落類型各選擇 3 塊20 m × 20 m 的樣地。在樣地中采集生長狀況良好、個(gè)體大小接近的蘆葦、水葫蘆和金魚藻三種植物各10株,將每個(gè)植株按東西南北四個(gè)方向的健康葉混勻作為一個(gè)樣品。采集的樣品清洗干凈后放入無菌密封袋中帶回實(shí)驗(yàn)室處理。在各樣地中,采用直徑5cm的土壤取樣器(土鉆),隨機(jī)采集3個(gè)濕地底泥沉積物樣品混勻?yàn)橐粋€(gè)樣品,裝入無菌密封袋后帶回實(shí)驗(yàn)室。δ13C同位素季節(jié)變化的樣品分別于2013年3月(春季)、6月(夏季)、9月(秋季)和12月(冬季)采集。

    2.2元素含量及δ13C同位素的分析

    土壤樣品分析參考孫英杰等(2018)的方法,有機(jī)質(zhì)(OM)采用重鉻酸鉀氧化法,全氮(TN)通過元素分析儀(德國Vario ELIII)測(cè)定,全磷(TP)用濃硫酸-高氯酸消煮,鉬銻抗比色法(美國Agilent8453紫外-可見分光光度計(jì)),速效氮(AN)用堿解擴(kuò)散法,速效磷(AP)用碳酸氫鈉浸提,鉬銻抗比色法。植物整株生物量的計(jì)算,隨機(jī)采集水生植物,將實(shí)測(cè)樣方面積的生物量,轉(zhuǎn)換為單位面積的生物量(g·m-2)(滕秋梅等,2021)。將帶回的植物樣品洗凈后,用自來水沖洗干凈,再用去離子水淋洗2~3次,自然晾干,105 ℃殺青30 min 后,50 ℃ 烘箱至恒重,分別粉碎后過80 目尼龍篩,混合均勻,裝入無菌密封袋中備分析用。葉片碳和氮含量的測(cè)定:采用無機(jī)酸淋洗的方法對(duì)葉片無機(jī)碳進(jìn)行清除,然后在德國Elmentar VarioEL Ⅲ元素分析儀上測(cè)定。葉片磷含量的測(cè)定:采用HNO3-HClO4消化,在ICP-AES(IRIS IntrepidⅡ XSP)上測(cè)定,回收率在91%以上。葉片鈣(Ca)含量的測(cè)定:樣品經(jīng)HNO3-HF-HClO4微波輔助消解后(意大利MileStone ETHOSE型微波消解儀),利用PE Elan6000型電感耦合等離子體質(zhì)譜儀(ICP-MS,Perkin-Elmer,USA)測(cè)定。δ13C同位素采用DELTAplus/ XL 型穩(wěn)定同位素比值質(zhì)譜計(jì)(Finnigan MAT公司)進(jìn)行測(cè)定,碳同位素參考PDB國際標(biāo)準(zhǔn),測(cè)定精度± 0.2 ‰。

    2.3 植物利用HCO-3比例的估算

    參考吳沿友等(2011)和李瑞等(2015)利用水生植物δ13C值與碳酸酐酶間的關(guān)系建立的二端元模型,來估算植物利用HCO-3作為無機(jī)碳源進(jìn)行光合作用大小比例,計(jì)算公式如下:

    fB=(δT-δA)/(δB-δA)

    式中:δT為試驗(yàn)區(qū)水生植物葉片的δ13C值;δA為基本上不利用HCO-3作無機(jī)碳源、碳酸酐酶活力極低的植物葉片的δ13C值;δB為極少利用二氧化碳作碳源而以 HCO-3為主要無機(jī)碳源的δ13C值;fB則為植物利用的HCO-3占無機(jī)碳源比例。

    基于試驗(yàn)區(qū)水生植物δ13C值的特征,通過植物碳同位素與碳酸酐酶顯著正相關(guān)的二端元模型(吳沿友等,2011),結(jié)合參考章程等(2013)的估算方法,本研究選擇不利用碳酸氫根的懸鈴木葉片δ13C(-31.56‰)作為無機(jī)碳源的δA,選擇三種植物葉片中的最大值-19.19‰作為δB端元,因此,fB的計(jì)算公式修正為fB=(δT+31.56)/(-19.19+31.56)。通過fB 和水生植物的固碳量,計(jì)算出水生植物通過光合作用固定的HCO-3的碳量。

    2.4 數(shù)據(jù)處理

    數(shù)據(jù)統(tǒng)計(jì)在Mcrosoft Excel 2007上進(jìn)行,采用SPSS 19.0統(tǒng)計(jì)軟件進(jìn)行不同植物間元素含量的單因素方差分析。

    3 結(jié)果與分析

    3.1 三種典型濕地植物葉片C、N、P、Ca與δ13C含量特征

    C、N、P是植物正常生長發(fā)育需要的三種重要元素,其含量特征能反映植物的生長特點(diǎn)及其與周圍環(huán)境間的關(guān)系。從表2可以看出,濕地蘆葦葉片的碳含量最高,平均為(35.50±0.51)%,其次是金魚藻,平均為(33.80±0.51)%,水葫蘆的碳含量最低(34.62±0.64)%。三種植物葉片氮元素含量表現(xiàn)為蘆葦>水葫蘆>金魚藻,其中蘆葦葉片最高,為(3.64±0.30)%。葉片磷元素含量則表現(xiàn)為金魚藻>水葫蘆>蘆葦。鈣元素含量與碳元素含量表現(xiàn)出類似規(guī)律:蘆葦>金魚藻>水葫蘆,蘆葦葉片的鈣含量達(dá)到(2.24±0.15)%,而水葫蘆最小,為(0.94±0.05)%。三種不同生活型植物δ13C含量均值的大小范圍在-28.47‰~-21.69‰,平均值為-24.83‰,其中金魚藻的δ13C值最大,蘆葦?shù)摩?3C值最小。除了葉片C之外,葉片N、P、Ca和δ13C含量在不同生活型植物間差異顯著(P<0.05)。

    從表3可以看出,挺水植物蘆葦葉片δ13C的季節(jié)變化表現(xiàn)為秋季>夏季>冬季>春季,浮水植物水葫蘆葉片δ13C的季節(jié)變化為秋季>冬季>夏季>春季,蘆葦和水葫蘆的δ13C值在季節(jié)間無顯著差異(P>0.05)。沉水植物金魚藻葉片δ13C的季節(jié)變化則表現(xiàn)為秋季>夏季>春季>冬季,秋夏季分別與春冬季之間有顯著差異(P<0.05)。不同季節(jié)植物葉片δ13C值的大小在不同生活型植物間差異顯著,一致表現(xiàn)為金魚藻>水葫蘆>蘆葦。

    表4是與其他地區(qū)植物葉片δ13C值的比較,從中可以看出,本試驗(yàn)區(qū)典型植物的δ13C值接近于貴州茂蘭板寨水生植物(李瑞等,2015)、貴州喀斯特山區(qū)植物(楊成等,2007)和江蘇鹽城濱海鹽生植物的δ13C值。在不同地區(qū),因陸生植物,水生植物和地區(qū)環(huán)境間的差異,植物葉片δ13C值的大小范圍表現(xiàn)出不同。

    3.2 濕地植物葉片δ13C與C、N、P、C/N、C/P、N/P的關(guān)系

    對(duì)會(huì)仙喀斯特濕地不同生活型植物葉片δ13C與C、N、P、C/N、C/P、N/P的相關(guān)性進(jìn)行分析(表5),結(jié)果表明,三種濕地植物葉片的δ13C值與C、N和P元素含量間呈現(xiàn)顯著和極顯著正相關(guān)關(guān)系(P<0.05,P<0.01),相關(guān)系數(shù)分別為0.563、0.458和0.765。δ13C與C/N、C/P和N/P間呈現(xiàn)顯著和極顯著負(fù)相關(guān)關(guān)系(P<0.05,P<0.01),相關(guān)系數(shù)分別為-0.617、-0.395和-0.488。

    從表6可以看出,會(huì)仙喀斯特濕地底泥的養(yǎng)分特征顯著影響了三種植物葉片的δ13C值,簡(jiǎn)單相關(guān)分析結(jié)果表明,植物葉片δ13C值與底泥的有機(jī)質(zhì)、速效氮、總氮、速效磷和總磷含量間呈現(xiàn)顯著正相關(guān)(P<0.05,P<0.01)。由于底泥養(yǎng)分元素之間的可能存在共變性和相互影響,進(jìn)一步采用偏相關(guān)方法進(jìn)行了補(bǔ)充分析,結(jié)果表明葉片δ13C值也與底泥的有機(jī)質(zhì)、速效氮、總氮和總磷含量含量呈顯著正相關(guān)(P<0.05),與速效磷含量沒有顯著性相關(guān)。

    3.3 植物光合作用固定HCO-3碳量的估算

    從表7可以看出,三種水生植物的生物量和固碳量的大小關(guān)系依次為蘆葦>水葫蘆>金魚藻,總體上挺水植物蘆葦?shù)墓烫剂勘瘸了参锝痿~藻固碳量高。金魚藻利用HCO-3作為光合作用碳源的比例最大,f B數(shù)值為20.15%,大小關(guān)系表現(xiàn)為金魚藻>蘆葦>水葫蘆。蘆葦通過光合作用固定的HCO-3的碳量最大,為159.60 t·a-1·km-2。三種水生植物吸收利用HCO-3占無機(jī)碳源的比例(fB)平均值為18.19%,由光合作用固定的HCO-3碳量的平均值為59.88 t·a-1·km-2。

    4 討論

    4.1 三種水生植物葉片δ13C特征與養(yǎng)分元素含量的關(guān)系

    在會(huì)仙喀斯特濕地該試驗(yàn)區(qū)域中,三種典型濕地水生植物葉片碳含量平均值為34.64 %,低于中國東黃海海島濕地植物葉片的碳含量(郭超等,2018),也低于其他陸地植物葉片的碳含量(Elser et al.,2000;任書杰等,2012),說明會(huì)仙喀斯特濕地這三種代表性水生植物的碳含量相對(duì)較低,這與會(huì)仙喀斯特濕地外界環(huán)境的人為擾動(dòng)較大有關(guān),這類似于濱海濕地(李征等,2012)。會(huì)仙喀斯特濕地植物葉片C/N、C/P 平均值低于Elser et al.(2000)報(bào)道的全球植物葉片的平均值(C/N:23.8,C/P :300.9)及中國部分研究區(qū)域(任書杰等,2012;郭超等,2018),表明會(huì)仙喀斯特濕地水生植物植株單位養(yǎng)分供應(yīng)的能力較低,這與其特殊的喀斯特生境有關(guān),會(huì)仙喀斯特濕地的人為干擾較嚴(yán)重(徐廣平等,2019;滕秋梅等,2021),對(duì)植被生產(chǎn)力的提高產(chǎn)生了影響。δ13C值的大小關(guān)系表現(xiàn)為沉水植物>浮水植物>挺水植物,這與貴州茂蘭(李瑞等,2015)報(bào)道的挺水植物>浮水植物>沉水植物有所不同,也與草海(張強(qiáng)等,2012)、會(huì)仙(章程等,2013)、長江中游(黃亮等,2003)等地報(bào)道的沉水植物>挺水植物>浮水植物略有差異,尤其與會(huì)仙(章程等,2013)亦有差異,原因可能是與不同采樣區(qū)域間存在的異質(zhì)性有關(guān)。吳沿友等(2011)的研究表明,植物葉片δ13C值具有較大的變異性,也與碳酸酐酶的活力大小有很大關(guān)系。可見,由于植物外界生長微環(huán)境的不同,對(duì)會(huì)仙喀斯特濕地植物穩(wěn)定碳同位素分餾有較大的影響,浮水植物水葫蘆在大氣CO2等碳源不充足時(shí),可能還會(huì)利用一部分水體中的HCO-3,因此其δ13C值介于沉水植物金魚藻和挺水植物蘆葦之間。

    在會(huì)仙喀斯特濕地,三種典型水生植物δ13C平均值的變化范圍在-28.47‰~-21.69‰之間,平均值-24.83‰,與Keely et al.(1992)報(bào)道的水生植物的δ13C值波動(dòng)范圍為-50‰~-11‰相符,與長江中游湖泊水生植物的δ13C值(黃亮等,2003)和長江口潮灘蘆葦δ13C的平均值(王偉偉,2011)較為接近,高于王文文(2012)報(bào)道的黃河三角洲蘆葦δ13C值-27.7‰ ,該差異可能主要是由于區(qū)域環(huán)境條件的不同而引起的。本研究中,沉水植物金魚藻的δ13C最高,黃亮等(2003)的研究結(jié)果也表明,由于水體中無機(jī)碳δ13C組成較高,使得沉水植物葉子的δ13C大于挺水植物和浮水植物。因此,會(huì)仙喀斯特濕地水生植物的不同生活型和不同植物種間的差異,是導(dǎo)致植物葉片δ13C值呈現(xiàn)差異性的原因。

    類似陸地植物,水生植物中同樣分布著C3、C4 和 CAM 植物,但C3 型水生植物居多(Keeley,1998)。挺水植物可以直接利用大氣中的CO2進(jìn)行光合作用,其δ13C分布在-30‰~-24‰之間,與C3 植物的δ13C變化范圍相重合,相似光合作用途徑的水生植物種類還包括浮葉植物和漂浮植物(Keeley,1998)。挺水植物、浮葉植物和漂浮植物的葉片暴露在大氣中,δ13C主要受到大氣中CO2的影響,這使得這三種類型植物的δ13C變化相對(duì)較?。▌⒔鹆恋?,2020)。相反,沉水植物由于水中CO2供給量的不足,可以利用溶解的HCO-3作為碳源,導(dǎo)致其δ13C值較C3 植物偏重(Raven,1972),并且與C4 植物的δ13C變化范圍部分有所重疊。在本研究中,沉水植物金魚藻的碳同位素高于浮水植物水葫蘆和挺水植物蘆葦,這與其他研究結(jié)論相一致(Raven,1972;Madsen & Breinholt,1995),主要原因可能是由于金魚藻葉片在光合作用中,利用了水體中的 CO2、HCO-3 等多個(gè)碳源,前人已有的研究也證明HCO-3比CO2更能富集δ13C(Mook et al.,1974)。

    在本研究中,金魚藻能夠利用溶解在水中的二氧化碳,還可能利用水體中的HCO-3,水葫蘆主要利用大氣中的二氧化碳,選擇性的利用一部分水體中的HCO-3,而蘆葦則主要利用大氣中的二氧化碳。由于浮水植物和挺水植物的碳源多為空氣中的 CO2氣體,空氣中CO2的碳同位素值在各個(gè)季節(jié)一般無差異,在各自的光合作用中沒產(chǎn)生顯著的同位素效應(yīng)(劉勇麗,2015),因此,本研究中的水葫蘆和蘆葦葉片δ13C值的季節(jié)變化間無顯著差異性(P>0.05)。前人研究表明,在溫度較高的夏季和秋季,馬來眼子菜(Potamogeton wrightii)、狐尾藻(Myriophyllum verticillatum)等沉水植物與金魚藻相互競(jìng)爭(zhēng)碳源,引起水體中的碳源不足,碳同位素的分餾受到抑制,金魚藻從而富集較多δ13C(Maberly & Madsen,1998;劉勇麗,2015),與上述研究結(jié)果類似,本該研究中沉水植物金魚藻葉片δ13C值季節(jié)之間有明顯的變化,主要在秋夏兩季分別與春冬兩季之間有顯著差異(P<0.05),秋季和夏季金魚藻葉片的δ13C值較高。

    一般C3植物葉片氮含量與葉片δ13C值呈正相關(guān)(Cernusak et al.,2013),與C/N(長期氮素利用效率)之間存在負(fù)相關(guān)(Gong et al.,2011),但在外界環(huán)境因素干擾下會(huì)發(fā)生變化(Sheng et al.,2011)。δ13C 隨葉片N濃度的增加而增加,原因是植物光合能力隨葉片 N 濃度的增加而增強(qiáng),由于δ13C與光合能力呈正比,所以葉片N與δ13C呈正相關(guān)(Zhou et al.,2016;Ma et al.,2018);馬劍英等(2008)研究表明紅砂植物葉片P與δ13C正相關(guān);Zhou et al.(2016)的研究結(jié)果表明,δ13C與P存在負(fù)相關(guān)關(guān)系。在本研究中,三種植物葉片δ13C與N、P含量顯著正相關(guān)(表4),表明三種植物可能通過采取N、P光合相關(guān)型的策略來適應(yīng)喀斯特濕地的環(huán)境。在土壤養(yǎng)分高的環(huán)境中,植物光合速率較快,植物葉片的δ13C值較大(Domingues et al.,2010),本研究結(jié)果也支持此結(jié)論,而三種植物葉片的δ13C相對(duì)不是很高,這也說明會(huì)仙喀斯特濕地養(yǎng)分相對(duì)較貧瘠。三種植物葉片δ13C值與底泥的有機(jī)質(zhì)、速效氮、總氮、速效磷和總磷含量呈顯著正相關(guān),表明底泥營養(yǎng)元素是影響三種植物葉片δ13C變化的重要因素。由于不同植物對(duì)外界環(huán)境的適應(yīng)性不同,因此其δ13C值與養(yǎng)分元素的相關(guān)性存在強(qiáng)弱的差異性(Wei et al.,2008)。

    4.2 三種植物葉片葉片

    δ13C與C、N、P元素化學(xué)計(jì)量學(xué)的關(guān)系

    植物δ13C能夠反映與植物生長過程相關(guān)聯(lián)的環(huán)境變化的信息,對(duì)環(huán)境變化有一定的指示意義(劉賢趙等,2014)。葉片氮含量可以反映植物碳、氮分配的經(jīng)濟(jì)性,碳氮比可以表征植物的長期氮素利用效率(展小云等,2012)?;谫|(zhì)量的植物葉片氮含量與葉片δ13C值間的關(guān)系具有不確定性,有報(bào)道指出二者呈現(xiàn)顯著正相關(guān)性(Zhang et al.,2015),也有研究結(jié)果表明有弱的負(fù)相關(guān)性(Hultine & Marshall,2000)。本研究結(jié)果表明,植物δ13C值與N元素含量間呈顯著正相關(guān)關(guān)系(P<0.05)。一般而言,植物對(duì)水分與氮素的利用效率存在優(yōu)化適應(yīng)策略,植物葉片的δ13C(長期水分利用效率)與葉碳氮比(長期氮素利用效率)之間呈現(xiàn)負(fù)相關(guān)性,表明擁有較高水分利用效率的植物多采用低氮素利用效率的適應(yīng)策略(陳清等,2015)。這是由于植物固定單位碳一般需要獲取更多氮,如采用低氮素利用效率的適應(yīng)策略,暗示植物則必須增大其水分利用效率(WUE)以代償減少的N利用率(Zhou et al.,2016;Wu et al.,2020)。本研究的結(jié)論也支持上述理論,會(huì)仙喀斯特濕地三種典型植物的δ13C與C/N呈極顯著負(fù)相關(guān)(P<0.01),說明會(huì)仙喀斯特濕地在外界人為干擾的影響下,這三種典型濕地植物為了保持其較高的水分利用效率,可能采取了低氮素利用效率的適應(yīng)策略。

    N/P常用來表征植物養(yǎng)分虧缺狀況,植物N/P<14表明植物生長受N限制,植物N/P>16,表明植物生長受P限制(Cernusak et al.,2013;Jia et al.,2018)。本研究結(jié)果說明,蘆葦N/P為22.21,可能主要受P限制;水葫蘆和金魚藻N/P<14,則可能主要受到N的限制;而這三種濕地植物葉片N/P平均值為10.34,低于全球陸生植物平均值13.8(Elser et al.,2000),說明N、P可能是共同制約植物生長的因素,也說明這三種典型植物的N、P養(yǎng)分虧缺成為植被生態(tài)恢復(fù)的影響因子。水生植物穩(wěn)定碳同位素(δ13C)的分餾過程較復(fù)雜,水生植物的δ13C主要與其吸收利用的碳源和水體環(huán)境要素等密切相關(guān)(劉金亮等,2020)。三種植物葉片δ13C與N/P含量顯著正相關(guān),進(jìn)一步表明三種植物通過與N、P光合相關(guān)的策略來適應(yīng)外界環(huán)境,與C/P成反比,說明三種植物可能主要通過提高P的利用率來促進(jìn)在低水分利用率環(huán)境下的碳的合成,暗示出外界環(huán)境壓力對(duì)不同生長環(huán)境下植物的氮、磷的適應(yīng)策略產(chǎn)生了較大的影響。會(huì)仙喀斯特濕地三種典型植物葉片δ13C特征與C/N、N/P和C/P的相互關(guān)系,與濱海鹽生植物的響應(yīng)特征相類似(夏大娟等,2020)。

    4.3 三種植物光合作用固定HCO-3碳量的估算

    水生植物利用的碳主要源于大氣、水體和沉積物的孔隙水中(Milligan et al.,2010),因不同碳源具有不同的碳同位素信息,不同類型的水生植物利用無機(jī)碳的方式也不一樣(劉金亮等,2020)。水生植物在進(jìn)行光合作用時(shí),多具有利用HCO-3作為碳源的能力(Prins & Elzenga,1989)。水生植物通過碳酸酐酶來利用HCO-3進(jìn)行光合作用,與Ca2+、HCO-3、pH等發(fā)生化學(xué)反應(yīng)(陳波等,2014;張?zhí)盏龋?014),引起了植物的穩(wěn)定碳同位素分餾(吳沿友等,2011),通過穩(wěn)定碳同位素與碳酸酐酶間的關(guān)系建立的二端元模型,可以估算水生植物光合作用利用HCO-3比例和固碳量(章程等,2013;李瑞等,2015)??λ固氐貐^(qū)淡水的pH 值多高于7,溶解的無機(jī)碳形式主要以HCO-3為主(Madsen,1983)。在本研究試驗(yàn)區(qū)域中,pH值在7.13~7.36之間,三種典型水生植物利用HCO-3作為光合作用碳源的比例的平均值為28.19%,低于貴州茂蘭(李瑞等,2015)的平均值47.84%,也低于章程等(2013)在會(huì)仙濕地其他區(qū)域計(jì)算的46.59%,究其原因,這可能與同一種植物生長在不同的水環(huán)境中有不同的δ13 C有關(guān)(林清和王紹令,2001),水生植物的δ13 C主要與其吸收利用的碳源和水體環(huán)境要素等密切相關(guān)(劉金亮等,2020)。

    本研究試驗(yàn)區(qū)中,挺水植物蘆葦?shù)纳锪亢凸烫剂孔畲?,其通過光合作用固定的HCO-3的碳量也最大,盡管沉水植物金魚藻的fB數(shù)值最大,但由于其生物量相對(duì)不高,固碳量相對(duì)較低,因而其光合作用固定的HCO-3的碳量最小,但其仍表現(xiàn)出一定的HCO-3固碳量。這從側(cè)面說明,由于人為對(duì)會(huì)仙濕地沉水植物金魚藻等的刈割和破壞,降低了金魚藻的生物量和HCO-3固碳量。三種典型水生植物固定的HCO-3碳量的平均值為59.88 t·a-1·km-2,低于章程等(2013)在會(huì)仙濕地其他區(qū)域計(jì)算的76.74 t·a-1·km-2,高于貴州茂蘭(李瑞等,2015)的平均值19.52 t·a-1·km-2??梢娍λ固貪竦厮参锿ㄟ^光合作用,可以將碳酸鹽巖風(fēng)化形成的HCO-3轉(zhuǎn)化為有機(jī)碳并且儲(chǔ)存在植物體內(nèi),水生植物通過光合作用表現(xiàn)出一定的碳匯效應(yīng)(章程等,2013)。同時(shí),本研究也表明,由于受到植被生活型、環(huán)境調(diào)控等因素的影響,加之植物生長微環(huán)境的差異,使不同植物在光合作用中對(duì)不同碳源的利用程度有所不同,產(chǎn)生同位素分餾,即使同一種植物在同一大的區(qū)域中也會(huì)表現(xiàn)出一定的同位素差異。另外,由于會(huì)仙喀斯特濕地水體中動(dòng)物可能會(huì)吞食部分水生植物,因而估算的這三種植物HCO-3的固碳量可能會(huì)偏小。水體環(huán)境因素(可溶性無機(jī)碳濃度、pH值、溫度和營養(yǎng)鹽含量等)、植物自身的生理學(xué)特征等對(duì)植物葉片δ13C的影響,需要后續(xù)進(jìn)一步深入的探討。

    5 結(jié)論

    會(huì)仙喀斯特濕地三種不同生活型水生植物δ13C值的大小范圍在-28.47‰~-21.69‰之間,平均值為-24.83‰,其中金魚藻的δ13C值最大,蘆葦?shù)摩?3C值最小。三種濕地植物葉片N/P 平均值為10.34,表現(xiàn)出其受N、P共同限制的特征。植物δ13C值與C、N和P元素含量間均呈顯著正相關(guān)關(guān)系(P<0.05),與C/N、C/P和N/P均呈顯著負(fù)相關(guān)關(guān)系(P<0.05,P<0.01),表明三種水生植物可能通過增加磷利用效率來促進(jìn)低水分利用率環(huán)境下的碳的合成,采用提高水分利用效率的策略來代償減小的氮素利用效率。基于測(cè)定的δ13C值,推測(cè)會(huì)仙喀斯特濕地中挺水植物蘆葦和漂浮植物水葫蘆是C3型植物,沉水植物金魚藻選擇擬C4的光合作用途徑。會(huì)仙喀斯特濕地三種植物的不同生活型、光合作用的途徑、水體底泥養(yǎng)分元素含量和生長微環(huán)境,共同促使植物同位素特征表現(xiàn)出一定的差異性,三種典型水生植物具有一定的碳匯效應(yīng)。

    致謝 感謝于倩倩、田壘、劉建春和程桂霞等在試驗(yàn)樣品分析方面提供的幫助!

    參考文獻(xiàn):

    CERNUSAK LA, UBIERNA N, WINTER K, et al., 2013. Environmental and physiological determinants of carbon isotope discrimination in terrestrial plants [J]. New Phytol, 200(4): 950-965.

    CHEN B, YANG R, LIU ZH, et al., 2014. Effects of aquatic phototrophs on diurnal hydrochemical and δ13CDIC variations in an epikarst spring and two spring-fed ponds of Laqiao, Maolan, SW China [J]. Geochim, 43(4): 375-385. ?[陳波, 楊睿, 劉再華, 等, 2014. 水生光合生物對(duì)茂蘭拉橋泉及其下游水化學(xué)和δ13CDIC晝夜變化的影響 [J]. 地球化學(xué), 43(4): 375-385.]

    CHEN Q, WANG YD, GUO CC, et al., 2015. Foliar stable carbon isotope ratios of Phragmites australis and the relevant environmental factors in marsh wetlands in Tianjin [J]. Chin J Plan Ecol, 39(11): 1044-1052. ?[陳清, 王義東, 郭長城, 等, 2015. 天津沼澤濕地蘆葦葉片的碳穩(wěn)定同位素比值分布特征及其環(huán)境影響因素 [J]. 植物生態(tài)學(xué)報(bào), 39(11): 1044-1052.]

    DIJKSTRA FA, CARRILLO Y, ASPINWALL MJ, et al., 2016. Waters, nitrogen and phosphorus use efficiencies of four tree species in response to variable water and nutrient supply [J]. Plant Soil, 406(1-2): 187-199.

    DOMINGUES TF, MEIR P, FELDPAUSCH TR, et al., 2010. Colimitation of photosynthetic capacity by nitrogen and phosphorus in West Africa woodlands [J]. Plant Cell Environ, 33, 959-980.

    ELSER JJ, FAGAN WF, DENNO RF, et al., 2000. Nutritional constraints in terrestrial and freshwater food webs [J]. Nature, 408(6812): 578-580.

    GONG XY, CHEN Q, LIN S, et al., 2011. Tradeoffs between nitrogen and water use efficiency in dominant species of the semiarid steppe of Inner Mongolia [J]. Plant Soil, 340: 227-238.

    GUO C, TUO B, SU T, et al., 2018. Carbon, nitrogen and phosphorus stoichiometry of five common herbaceous species across islands in the Yellow Sea and the East China Sea [J]. Chin J App Ecol, 29(2): 380-388. ?[郭超, 妥彬, 蘇田, 等, 2018. 中國東黃海海島 5 種常見草本的碳氮磷化學(xué)計(jì)量特征 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 29(2): 380-388.]

    HAN Q, WANG BL, 2017. Mechanisms of stable carbon isotope fractionation by phytoplankton and its environmental applications [J]. Chin J Ecol, 36(5): 1436-1443. ?[韓瓊, 王寶利, 2017. 浮游植物穩(wěn)定碳同位素分餾機(jī)制及環(huán)境應(yīng)用 [J]. 生態(tài)學(xué)雜志, 36(5): 1436-1443.]

    HUANG L, WU Y, ZHANG J, et al., 2003. Distribution of C, N, P and δ13C in aquatic plants of some lakes in the Middle Yangtze Valley [J]. Acta Geosci Sin, 24(6): 515-518. ?[黃亮, 吳瑩, 張經(jīng), 等, 2003. 長江中游若干湖泊水生植物體內(nèi)C、N、P及δ13C分布 [J]. 地球?qū)W報(bào), 24(6): 515-518.]

    HULTINE KR, MARSHALL JD, 2000. Altitude trends in conifer leaf morphology and stable carbon isotope composition [J]. Oecologia, 123: 32-40.

    JIA J, BAI JH, WANG W, et al., 2018. Changes of biogenic elements in Phragmites australis and Suaeda salsa from salt marshes in Yellow River Delta,China [J]. Chin Geogr Sci, 28(3): 411-419.

    KEELY JE, 1998. CAM photosynthesis in submerged aquatic plants [J]. Bot Rev, 64(2): 121-175.

    KEELY JE, SADQUIST DR, 1992. Carbon: frest water plants [J]. Plant Cell Environ, 15: 1021-1035.

    LI R, YU S, SUN PA, et al., 2015. Characteristics of δ13C in typical aquatic plants and carbon sequestration by plant photosynthesis in the Banzhai catchment, Maolan of Guizhou Province [J]. Carsol Sin, 34(1): 9-16. ?[李瑞, 于奭, 孫平安, 等, 2015. 貴州茂蘭板寨水域水生植物δ13C特征及光合作用固定HCO-3碳量估算 [J]. 中國巖溶, 34(1): 9-16.]

    LI Z, HAN L, LIU YH, et al., 2012. C, N and P stoichiometric characteristics in leaves of Suaeda salsa during different growth phase in coastal wetlands of China [J]. Chin J Plan Ecol, 36(10): 1054-1061. ?[李征, 韓琳, 劉玉虹, 等, 2012. 濱海鹽地堿蓬不同生長階段葉片 C、N、P 化學(xué)計(jì)量特征 [J]. 植物生態(tài)學(xué)報(bào), 36(10): ?1054-1061.]

    LIN Q, WANG SL, 2001. The composition of stable carbon isotope and some influencing factors of submerged plant [J]. Acta Ecol Sin, 21(5): 806-809. ?[林清, 王紹令, 2001. 沉水植物穩(wěn)定碳同位素組成及影響因素分析 [J]. 生態(tài)學(xué)報(bào), 21(5): 806-809.]

    LIU JL, XUE B, YAO SC, et al., 2020. Mechanisms of stable carbon isotope fractionation of aquatic plants and the research advances of application [J]. Acta Ecol Sin, 40(8): 2533-2544. ?[劉金亮, 薛濱, 姚書春, 等, 2020. 湖泊水生植物穩(wěn)定碳同位素分餾機(jī)制與應(yīng)用研究進(jìn)展 [J]. 生態(tài)學(xué)報(bào), 40(8): 2533-2544.]

    LIU XZ, ZHANG Y, SU Q, et al., 2014. Research progress in responses of modern terrestrial plant carbon isotope composition to climate change [J]. Adv Earth Sci, 29(12): 1341-1354. ?[劉賢趙, 張勇, 宿慶, 等, 2014. 現(xiàn)代陸生植物碳同位素組成對(duì)氣候變化的響應(yīng)研究進(jìn)展 [J]. 地球科學(xué)進(jìn)展, 29(12): 1341-1354.]

    LIU YL, 2015. On the stable carbon and nitrogen isotope ecology of aquatic macrophytes in Lake Taihu watershed [D]. Xiangtan: Xiangtan University. ?[劉勇麗, 2015. 太湖流域水生植物穩(wěn)定碳氮同位素的生態(tài)學(xué)研究 [D]. 湘潭: 湘潭大學(xué).]

    MA F, LIANG WY, ZHOU ZN, et al., 2018. Spatial ariation in leaf stable carbon isotope composition of three caragana species in northern China [J]. Forests, 9(6): 297.

    MA JY, CHEN FH, XIA DS, et al., 2008. Correlations between leaf δ13C and physiological parameters of desert plant Reaumuria soongorica [J]. Chin J App Ecol, 19(5): 1166-1171. ?[馬劍英, 陳發(fā)虎, 夏敦勝, 等, 2008. 荒漠植物紅砂葉片δ13 C值與生理指標(biāo)的關(guān)系 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 19(5): 1166-1171.]

    MABERLY SC, MADSEN TV, 1998. Affinity for CO2 in relation to the ability of freshwater macrophytes to use HCO-3 [J]. Funct Ecol, 12(1): 99-106.

    MADSEN TV, 1983. Growth and photosynthetic acclimation by Ranunculus aqutuatilis L.in response to inorgnic carbon availability [J]. New Phytol, 125(4): 707-715.

    MADSEN TV, BREINHOLT M, 1995. Effects of air contact on growth, inorganic carbon sources, and nitrogen uptake by an amphibious freshwater macrophyte [J]. Plant Physiol, 107(1): 149-154.

    MILLIGAN HE, PRETZLAW TD, HUMPHRIES MM, 2010. Stable isotope differentiation of freshwater and terrestrial vascular plants in two subarctic regions [J]. Acoscience, 17(3): 265-275.

    MOOK WG, BOMMERSON JC, STAVERMAN WH, 1974. Carbon isotope fraction between dissolved bicarbonate and gaseous carbon dioxide [J]. Earth Planet Sci Lett, 22(2): 169-176.

    PRINS HBA, ELZENGA JTM, 1989. Bicarbonate utilization: function and mechanism [J]. Aquat Bot, 34(1-3): 59-83.

    QU CM, GAN XG, SU B, et al., 2001. The characteristics of foliar δ13C values of plants and plant water use efficiency indicated by δ13C values in two fragmented rainforests in Xishuangbanna, Yunnan [J]. Acta Bot Sin, ?43(2): 186-192. ?[渠春梅, 韓興國, 蘇波, 等, 2001. 云南西雙版納片斷化熱帶雨林植物葉片δ13 C值的特點(diǎn)及其對(duì)水分利用效率的指示 [J]. 植物學(xué)報(bào), 43(2): 186-192.]

    RAVEN JA, 1972. Endogenous inorganic carbon sources in plant photosynthesis [J]. New Phytol, 71(2): 227-247.

    REN SJ, YU GR, JIANG CM, et al., 2012. Stoichiometric characteristics of leaf carbon, nitrogen, and phosphorus of 102 dominant species in forest ecosystems along the North-South Transect of East China [J]. Chin J App Ecol, 23(3): 581-586. ?[任書杰, 于貴瑞, 姜春明, 等, 2012. 中國東部南北樣帶森林生態(tài)系統(tǒng)102個(gè)優(yōu)勢(shì)種葉片碳氮磷化學(xué)計(jì)量學(xué)統(tǒng)計(jì)特征 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 23(3): 581-586.]

    SHENG WP, REN SJ, YU GR, et al., 2011. Patterns and driving factors of WUE and NUE in natural forest ecosystems along the Northsouth Transect of Eastern China [J]. J Geogr Sci, 21: 651-665.

    SPRENT JI, GEOGHEGAN IE, WHITTY PW, et al., 1996. Natural abundance of 15N and 13C in nodulated legumes and other plants in the cerrado and neighbouring regions of Brazil [J]. Oecology, 105(4): 440-446.

    SUN YJ, XU GP, SHEN YY, et al., 2018. Soil enzyme activities of phragmites australis community area in Huixian karst wetland, Guilin [J]. Wetland Sci, 16(2): 196-203. ?[孫英杰, 徐廣平, 沈育伊, 等, 2018. 桂林會(huì)仙喀斯特濕地蘆葦群落區(qū)土壤酶活性 [J]. 濕地科學(xué), 16(2): 196-203.]

    TENG QM, SUN YJ, SHEN YY, et al., 2021. Growth and phenotypic plasticity variability of Eichhornia crassipes in response to different eutrophic water in karst wetland [J]. J Lake Sci, 33(1): 123-137. ?[滕秋梅, 孫英杰, 沈育伊, 等, 2021. 喀斯特濕地水葫蘆(Eichhornia crassipes)生長及表型可塑性特征對(duì)不同富營養(yǎng)化水體的響應(yīng) [J]. 湖泊科學(xué), 33(1): 123-137.]

    VAN HM, CHAKRABORTY A, DAVIES KL, et al., 2018. The stable isotope composition of organic and inorganic fossils in lake sediment records: Current understanding, challenges, and future directions [J]. Quaternary Sci Rev, 196: 154-176.

    WANG WW, 2011. Nutrient dynamic and stable isotope indicators in tidal flats of the Yangtze Estuary [D]. Shanghai: East China Normal University. ?[王偉偉, 2011. 長江口潮灘營養(yǎng)動(dòng)態(tài)與穩(wěn)定同位素指示研究 [D]. 上海: 華東師范大學(xué).]

    WANG WW, 2012. Response of δ13C values of wetland ecosystem plants to soil salinity in Yellow River Delta [D]. Yantai: Ludong University. ?[王文文, 2012.黃河三角洲濕地生態(tài)系統(tǒng)植物δ13C值對(duì)土壤鹽分的響應(yīng) [D]. 煙臺(tái): 魯東大學(xué).]

    WEI L, YAN CL, GUO XY, et al., 2008. Variation in the δ13C of two mangrove plants is correlated with stomatal response to salinity [J]. J Plant Growth Regul, 27: 263-269.

    WU HJ, YANG T, LIU X, et al., 2020. Towards an integrated nutrient management in crop species to improve nitrogen and phosphorus use efficiencies of Chaohu Watershed [J]. J Clean Prod, 272: 122765.

    WU YY, XING DK, LIU Y, 2011. The characteristics of bicarbonate used by plants [J]. Earth Environ, 39(2): 273-277. ?[吳沿友, 邢德科, 劉瑩, 2011. 植物利用碳酸氫根離子的特征分析 [J]. 地球與環(huán)境, 39(2): 273-277.]

    XIA DJ, LIU QR, LILIANN Z, et al., 2020. Foliar δ13C correlates with elemental stoichiometry in halophytes of coastal wetlands [J]. Acta Ecol Sin, 40(7): 2215-2224. ?[夏大娟, 劉秋蓉, Liliann Zou, 等, 2020. 濱海鹽生植物葉片δ13 C與主要養(yǎng)分元素化學(xué)計(jì)量的關(guān)系 [J]. 生態(tài)學(xué)報(bào), 40(7): 2215-2224.]

    XU GP, LI YQ, SHEN YY, et al., 2019. Soil organic carbon distribution and components in different plant communities along a water table gradient in the Huixian karst wetland in Guilin [J]. Chin J Environ Sci, 40(3): 1491-1503. ?[徐廣平, 李艷瓊, 沈育伊, 等, 2019. 桂林會(huì)仙喀斯特濕地水位梯度下不同植物群落土壤有機(jī)碳及其組分特征 [J]. 環(huán)境科學(xué), 40(3): 1491-1503.]

    YANG C, LIU CQ, SONG ZL, et al., 2007. Carbon isotope compositions of plant leaves from kasrst mountainous area, Guizhou, China [J]. Earth Environ, 35(1): 33-38. ?[楊成, 劉叢強(qiáng), 宋照亮, 等, 2007. 貴州喀斯特山區(qū)植物葉片碳同位素組成研究 [J]. 地球與環(huán)境, 35(1): 33-38.]

    ZHAN XY, YU GR, SHENG WP, et al., 2012. Foliar water use efficiency and nitrogen use efficiency of dominant plant species in main forests along the Northsouth Transect of East China [J]. Chin J App Ecol, 23(3): 587-594. ?[展小云, 于貴瑞, 盛文萍, 等, 2012. 中國東部南北樣帶森林優(yōu)勢(shì)植物葉片的水分利用效率和氮素利用效率 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 23(3): 587-594.]

    ZHANG C, XIE YQ, NING LD, et al., 2013. Characteristics of δ13C in typical aquatic plants and carbon sequestration in Huixian karst wetland, Guilin [J]. Carsol Sin, 32(3): 247-251. ?[章程, 謝運(yùn)球, 寧良丹, 等, 2013. 桂林會(huì)仙巖溶濕地典型水生植物δ13C特征與固碳量估算 [J]. 中國巖溶, 32(3): 247-251.]

    ZHANG JH, GU L, BAO FC, et al., 2015. Nitrogen control of 13C enrichment in heterotrophic organs relative to leaves in a landscape building desert plant species [J]. Biogeosciences, 12: 15-27.

    ZHANG Q, 2012. The stability of carbon sink effect related to carbonate rock dissolution: A case study of the Caohai lake geological carbon sink [J]. Acta Geosci Sin, 33(6): 947-952. ?[張強(qiáng), 2012. 巖溶地質(zhì)碳匯的穩(wěn)定性-以貴州草海地質(zhì)碳匯為例 [J]. 地球?qū)W報(bào), 33(6): 947-952.]

    ZHANG T, PU JB, YUAN DX, et al., 2014. Diel variations of hydrochemistry and influencing factors in a surface stream in subtropical karst area, SW China [J]. Chin J Environ Sci, 35(8): 102-109. ?[張?zhí)眨?蒲俊兵, 袁道先, 等, 2014. 亞熱帶典型巖溶區(qū)地表溪流水文地球化學(xué)晝夜變化及其影響因素研究 [J]. 環(huán)境科學(xué), 35(8): 102-109.]

    ZHOU YC, CHENG XL, FAN JW, et al., 2016. Relationships between foliar carbon isotope composition and elements of C3 species in grasslands of Inner Mongolia, China [J]. Plant Ecol, 217(7): 883-897.

    (責(zé)任編輯 李 莉)

    猜你喜歡
    水葫蘆喀斯特同位素
    水葫蘆
    大自然的鬼斧神工:“中國南方喀斯特”
    桃子湖的水葫蘆
    元素含量及穩(wěn)定同位素比值用于網(wǎng)銷帶魚產(chǎn)地溯源
    水葫蘆
    2017第四屆全國穩(wěn)定同位素制備與應(yīng)用技術(shù)交流會(huì)
    2016年生物高考有關(guān)同位素標(biāo)記法的實(shí)驗(yàn)考查
    熱脫附—單體穩(wěn)定同位素技術(shù)測(cè)定揮發(fā)性有機(jī)物
    水葫蘆的秘密
    穿越喀斯特
    亚洲中文av在线| 侵犯人妻中文字幕一二三四区| 久久天堂一区二区三区四区| 国产欧美日韩一区二区三| 欧美黑人精品巨大| 日韩欧美一区二区三区在线观看| 国产亚洲精品第一综合不卡| 黄色 视频免费看| 999久久久精品免费观看国产| 免费在线观看影片大全网站| 黑人巨大精品欧美一区二区蜜桃| 村上凉子中文字幕在线| 亚洲少妇的诱惑av| 一进一出抽搐gif免费好疼 | 黑人巨大精品欧美一区二区mp4| 色婷婷久久久亚洲欧美| 老司机亚洲免费影院| 99国产精品一区二区三区| 老司机深夜福利视频在线观看| 欧美乱码精品一区二区三区| xxxhd国产人妻xxx| 久久久久久久久免费视频了| 高清av免费在线| 高清av免费在线| 91麻豆av在线| 亚洲人成伊人成综合网2020| 久久久精品国产亚洲av高清涩受| 老司机午夜十八禁免费视频| 最新在线观看一区二区三区| 一级片'在线观看视频| 在线播放国产精品三级| 日日夜夜操网爽| 热re99久久国产66热| 亚洲一区中文字幕在线| 成年版毛片免费区| 亚洲va日本ⅴa欧美va伊人久久| 国产精品免费视频内射| 变态另类成人亚洲欧美熟女 | 日韩精品免费视频一区二区三区| 身体一侧抽搐| 老熟妇仑乱视频hdxx| 国产精品综合久久久久久久免费 | 窝窝影院91人妻| 久久 成人 亚洲| 中文字幕精品免费在线观看视频| 9热在线视频观看99| 免费在线观看视频国产中文字幕亚洲| 亚洲欧美一区二区三区久久| 国产成+人综合+亚洲专区| 国产亚洲精品一区二区www| 女性被躁到高潮视频| 欧美人与性动交α欧美软件| 中文字幕av电影在线播放| 高清在线国产一区| 久久精品国产99精品国产亚洲性色 | 香蕉丝袜av| 国产一区二区激情短视频| 国产精品日韩av在线免费观看 | 欧美日本中文国产一区发布| 欧美午夜高清在线| avwww免费| 久久精品aⅴ一区二区三区四区| 婷婷六月久久综合丁香| 亚洲国产看品久久| 80岁老熟妇乱子伦牲交| 久久中文看片网| 亚洲成av片中文字幕在线观看| 国产av一区二区精品久久| 岛国视频午夜一区免费看| 亚洲中文av在线| 91av网站免费观看| 中亚洲国语对白在线视频| 超色免费av| 后天国语完整版免费观看| 一区在线观看完整版| 欧美黄色淫秽网站| 夜夜夜夜夜久久久久| 亚洲人成伊人成综合网2020| 欧美一级毛片孕妇| 欧美中文综合在线视频| 亚洲欧洲精品一区二区精品久久久| 欧美大码av| 国产1区2区3区精品| 国产精品久久久久成人av| 老司机深夜福利视频在线观看| 欧美不卡视频在线免费观看 | 一级黄色大片毛片| 亚洲成人免费电影在线观看| 亚洲一区二区三区色噜噜 | 一个人免费在线观看的高清视频| 亚洲中文av在线| 51午夜福利影视在线观看| 国产黄色免费在线视频| 亚洲人成伊人成综合网2020| 亚洲精品国产一区二区精华液| 国产片内射在线| 亚洲午夜理论影院| 亚洲五月色婷婷综合| 久久人人97超碰香蕉20202| 他把我摸到了高潮在线观看| 国产一区二区激情短视频| 男女下面插进去视频免费观看| 三级毛片av免费| 成人影院久久| 亚洲av片天天在线观看| 身体一侧抽搐| 在线视频色国产色| 深夜精品福利| 亚洲七黄色美女视频| 国产成人av教育| 成人亚洲精品av一区二区 | 日韩欧美一区二区三区在线观看| 国产午夜精品久久久久久| 满18在线观看网站| bbb黄色大片| 高清av免费在线| 日日夜夜操网爽| 亚洲一区中文字幕在线| 日韩精品青青久久久久久| 日韩欧美三级三区| 69av精品久久久久久| 亚洲人成电影免费在线| 欧美日韩亚洲国产一区二区在线观看| 99国产精品99久久久久| 麻豆成人av在线观看| 神马国产精品三级电影在线观看 | 免费观看人在逋| 国产精品偷伦视频观看了| 极品人妻少妇av视频| 国产成人av激情在线播放| 757午夜福利合集在线观看| 午夜免费成人在线视频| 亚洲精品美女久久久久99蜜臀| 精品熟女少妇八av免费久了| 首页视频小说图片口味搜索| 高潮久久久久久久久久久不卡| 国产精品99久久99久久久不卡| www.自偷自拍.com| 免费少妇av软件| 久久久久久久午夜电影 | 国产1区2区3区精品| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 国产欧美日韩精品亚洲av| 大香蕉久久成人网| 亚洲国产中文字幕在线视频| 亚洲国产欧美网| 成人av一区二区三区在线看| a级片在线免费高清观看视频| 欧美成人性av电影在线观看| 成人国语在线视频| 国产99久久九九免费精品| 亚洲五月色婷婷综合| 色在线成人网| 美女高潮到喷水免费观看| 嫁个100分男人电影在线观看| 狂野欧美激情性xxxx| 亚洲av成人一区二区三| 黄色 视频免费看| 男女床上黄色一级片免费看| 1024视频免费在线观看| 国产97色在线日韩免费| 男男h啪啪无遮挡| 老司机靠b影院| 日韩欧美三级三区| 五月开心婷婷网| 高清毛片免费观看视频网站 | 亚洲精品久久午夜乱码| 一边摸一边做爽爽视频免费| 成人国产一区最新在线观看| 亚洲精品美女久久av网站| 亚洲人成电影观看| 国产欧美日韩一区二区精品| 高清在线国产一区| 变态另类成人亚洲欧美熟女 | 在线天堂中文资源库| 在线观看免费日韩欧美大片| 老司机福利观看| 亚洲国产欧美网| 日本免费一区二区三区高清不卡 | 日韩欧美免费精品| 99精品久久久久人妻精品| 人妻丰满熟妇av一区二区三区| 97人妻天天添夜夜摸| 精品一品国产午夜福利视频| 久久中文字幕一级| 三级毛片av免费| 日日摸夜夜添夜夜添小说| 色哟哟哟哟哟哟| 美女高潮喷水抽搐中文字幕| 又大又爽又粗| 桃色一区二区三区在线观看| 看片在线看免费视频| 黑人巨大精品欧美一区二区mp4| 日韩精品青青久久久久久| 国产亚洲欧美在线一区二区| 成人黄色视频免费在线看| 日本三级黄在线观看| 欧美精品啪啪一区二区三区| 久久人人精品亚洲av| 亚洲自拍偷在线| 久久精品人人爽人人爽视色| 很黄的视频免费| 国产精品一区二区免费欧美| 97人妻天天添夜夜摸| 亚洲专区中文字幕在线| 99在线人妻在线中文字幕| 成人手机av| 亚洲色图综合在线观看| 久久精品91蜜桃| 日本 av在线| 免费搜索国产男女视频| 成年女人毛片免费观看观看9| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 在线观看免费视频网站a站| 亚洲七黄色美女视频| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| 超碰97精品在线观看| 国产蜜桃级精品一区二区三区| cao死你这个sao货| 成人三级做爰电影| 咕卡用的链子| 久久欧美精品欧美久久欧美| 国产精品二区激情视频| 十八禁人妻一区二区| 在线播放国产精品三级| 女性被躁到高潮视频| 亚洲av成人不卡在线观看播放网| 国产精品日韩av在线免费观看 | 色婷婷av一区二区三区视频| 看片在线看免费视频| 亚洲精品国产区一区二| 亚洲熟妇中文字幕五十中出 | 成年版毛片免费区| 99re在线观看精品视频| 极品教师在线免费播放| 99久久精品国产亚洲精品| 91九色精品人成在线观看| 搡老岳熟女国产| 亚洲欧美一区二区三区黑人| 久久精品国产99精品国产亚洲性色 | 午夜免费观看网址| 成人手机av| 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 亚洲精品中文字幕一二三四区| 女人精品久久久久毛片| 国产伦一二天堂av在线观看| 亚洲国产欧美网| 波多野结衣高清无吗| 欧美日韩乱码在线| 亚洲avbb在线观看| 亚洲欧美日韩另类电影网站| 男女下面进入的视频免费午夜 | 亚洲免费av在线视频| 国产av一区在线观看免费| 老司机在亚洲福利影院| 国产免费av片在线观看野外av| 桃红色精品国产亚洲av| 成人18禁高潮啪啪吃奶动态图| 久热爱精品视频在线9| 亚洲成av片中文字幕在线观看| 国产成人一区二区三区免费视频网站| 亚洲中文日韩欧美视频| 久久久久久人人人人人| 美女高潮到喷水免费观看| 69精品国产乱码久久久| 高清在线国产一区| 在线观看66精品国产| 他把我摸到了高潮在线观看| 免费看十八禁软件| 精品午夜福利视频在线观看一区| 美女国产高潮福利片在线看| 精品人妻在线不人妻| 国产成人影院久久av| 黄片播放在线免费| 身体一侧抽搐| 天堂√8在线中文| 国产亚洲精品久久久久5区| 91麻豆av在线| 女同久久另类99精品国产91| 亚洲自拍偷在线| 国内毛片毛片毛片毛片毛片| 久久精品国产99精品国产亚洲性色 | 国产精品久久久久久人妻精品电影| 日本三级黄在线观看| 超碰97精品在线观看| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 亚洲中文字幕日韩| 宅男免费午夜| 日韩有码中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人精品中文字幕电影 | 人妻久久中文字幕网| 国产视频一区二区在线看| 国产欧美日韩精品亚洲av| 操出白浆在线播放| 女人爽到高潮嗷嗷叫在线视频| 在线视频色国产色| 一级片免费观看大全| 黑人欧美特级aaaaaa片| 少妇粗大呻吟视频| 99国产综合亚洲精品| 精品日产1卡2卡| 精品国产超薄肉色丝袜足j| 亚洲av成人av| 国产精品美女特级片免费视频播放器 | 黄片小视频在线播放| 人人妻人人澡人人看| 午夜福利影视在线免费观看| 日韩欧美一区视频在线观看| 成年人黄色毛片网站| 国产成人精品久久二区二区免费| 午夜免费成人在线视频| 免费人成视频x8x8入口观看| 国产aⅴ精品一区二区三区波| 午夜视频精品福利| 日韩欧美免费精品| 久久国产精品人妻蜜桃| 老熟妇仑乱视频hdxx| 狠狠狠狠99中文字幕| 国产亚洲欧美精品永久| 亚洲专区国产一区二区| 日本黄色日本黄色录像| 久久久久久亚洲精品国产蜜桃av| 欧美日韩精品网址| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 国产精品一区二区精品视频观看| 久久久水蜜桃国产精品网| 精品日产1卡2卡| 亚洲三区欧美一区| 亚洲av日韩精品久久久久久密| 成人三级做爰电影| 欧美 亚洲 国产 日韩一| 在线观看一区二区三区激情| 久久久久国内视频| 亚洲在线自拍视频| 午夜影院日韩av| 亚洲五月天丁香| 精品高清国产在线一区| 免费少妇av软件| 在线观看免费日韩欧美大片| 丁香六月欧美| 久久狼人影院| 久久九九热精品免费| 免费在线观看视频国产中文字幕亚洲| 欧美午夜高清在线| 国产极品粉嫩免费观看在线| 两性午夜刺激爽爽歪歪视频在线观看 | 在线观看66精品国产| 国产午夜精品久久久久久| 老司机福利观看| 亚洲欧美日韩另类电影网站| 五月开心婷婷网| 亚洲第一欧美日韩一区二区三区| 在线观看www视频免费| 亚洲国产精品一区二区三区在线| 国产日韩一区二区三区精品不卡| 亚洲专区中文字幕在线| 精品福利永久在线观看| 国产精品一区二区免费欧美| 女人高潮潮喷娇喘18禁视频| 视频在线观看一区二区三区| 午夜91福利影院| 他把我摸到了高潮在线观看| 成人手机av| 亚洲专区字幕在线| 18禁国产床啪视频网站| 日本wwww免费看| 一边摸一边做爽爽视频免费| 国产有黄有色有爽视频| xxx96com| 久热这里只有精品99| 久久久精品国产亚洲av高清涩受| 色婷婷久久久亚洲欧美| 亚洲成av片中文字幕在线观看| 久久久国产成人精品二区 | 欧美精品一区二区免费开放| 久久精品亚洲熟妇少妇任你| 久久精品影院6| 黑人猛操日本美女一级片| 国产深夜福利视频在线观看| av网站免费在线观看视频| 可以在线观看毛片的网站| 国产欧美日韩一区二区三| 亚洲免费av在线视频| 少妇 在线观看| 天天添夜夜摸| 成人18禁高潮啪啪吃奶动态图| 在线观看一区二区三区激情| 电影成人av| 亚洲国产欧美网| 午夜免费观看网址| 欧美中文日本在线观看视频| 窝窝影院91人妻| 丝袜在线中文字幕| 亚洲精品中文字幕一二三四区| 久久午夜综合久久蜜桃| 纯流量卡能插随身wifi吗| av超薄肉色丝袜交足视频| 日日夜夜操网爽| 一a级毛片在线观看| 真人一进一出gif抽搐免费| 亚洲欧美一区二区三区久久| 岛国视频午夜一区免费看| 久99久视频精品免费| 久久精品亚洲熟妇少妇任你| 99精品欧美一区二区三区四区| 国产单亲对白刺激| 黄网站色视频无遮挡免费观看| 国产亚洲欧美在线一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 每晚都被弄得嗷嗷叫到高潮| 母亲3免费完整高清在线观看| 国产精品电影一区二区三区| 精品国产一区二区三区四区第35| 久久这里只有精品19| 亚洲一区二区三区色噜噜 | 亚洲欧美精品综合久久99| 久久人人97超碰香蕉20202| 丝袜人妻中文字幕| 91在线观看av| a级片在线免费高清观看视频| 欧美丝袜亚洲另类 | 少妇被粗大的猛进出69影院| 久久人人爽av亚洲精品天堂| 在线观看舔阴道视频| 在线永久观看黄色视频| 亚洲成a人片在线一区二区| 侵犯人妻中文字幕一二三四区| 成在线人永久免费视频| 9色porny在线观看| 波多野结衣av一区二区av| 99riav亚洲国产免费| 国产亚洲精品一区二区www| 免费在线观看完整版高清| 在线av久久热| 亚洲精品中文字幕在线视频| 母亲3免费完整高清在线观看| 国产亚洲精品综合一区在线观看 | 亚洲九九香蕉| xxxhd国产人妻xxx| 婷婷六月久久综合丁香| 十分钟在线观看高清视频www| 亚洲精品国产区一区二| 久久中文字幕一级| 亚洲在线自拍视频| 日韩 欧美 亚洲 中文字幕| 精品无人区乱码1区二区| 不卡一级毛片| 亚洲精华国产精华精| 国产精品久久久久成人av| 久久香蕉国产精品| 国产高清videossex| 波多野结衣av一区二区av| 欧美日韩福利视频一区二区| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 亚洲午夜理论影院| 在线天堂中文资源库| 国产精品 欧美亚洲| 亚洲成人精品中文字幕电影 | 中文字幕av电影在线播放| 亚洲,欧美精品.| 亚洲精品国产一区二区精华液| 免费看a级黄色片| 制服诱惑二区| 水蜜桃什么品种好| 日本a在线网址| 午夜91福利影院| 欧美激情极品国产一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 婷婷精品国产亚洲av在线| 99精品久久久久人妻精品| 亚洲成人免费电影在线观看| 人成视频在线观看免费观看| 日韩免费高清中文字幕av| 视频区欧美日本亚洲| 亚洲情色 制服丝袜| 男女下面进入的视频免费午夜 | 国产伦一二天堂av在线观看| 精品人妻在线不人妻| 两性夫妻黄色片| www国产在线视频色| 久久久久九九精品影院| 99精品久久久久人妻精品| 在线观看日韩欧美| 法律面前人人平等表现在哪些方面| 欧美日韩精品网址| 美女午夜性视频免费| 美女 人体艺术 gogo| 黄色毛片三级朝国网站| 亚洲精品久久成人aⅴ小说| 亚洲第一欧美日韩一区二区三区| 乱人伦中国视频| 亚洲人成伊人成综合网2020| 欧美+亚洲+日韩+国产| 69精品国产乱码久久久| 国产99白浆流出| 欧美乱妇无乱码| 亚洲中文字幕日韩| 亚洲五月天丁香| 一边摸一边抽搐一进一出视频| 日本免费a在线| 丝袜美足系列| 五月开心婷婷网| 一级毛片精品| 十八禁网站免费在线| 午夜免费鲁丝| 深夜精品福利| 两个人看的免费小视频| 国产91精品成人一区二区三区| 国产片内射在线| 免费av毛片视频| 成人影院久久| 亚洲激情在线av| 亚洲欧洲精品一区二区精品久久久| av欧美777| √禁漫天堂资源中文www| 美女午夜性视频免费| av网站在线播放免费| 精品一区二区三卡| 精品欧美一区二区三区在线| 久久国产乱子伦精品免费另类| av片东京热男人的天堂| 国产av在哪里看| 日本vs欧美在线观看视频| 伊人久久大香线蕉亚洲五| 99热国产这里只有精品6| 91字幕亚洲| 老司机午夜十八禁免费视频| 好男人电影高清在线观看| 视频区图区小说| 乱人伦中国视频| 免费不卡黄色视频| 男女午夜视频在线观看| 最新在线观看一区二区三区| 村上凉子中文字幕在线| 性欧美人与动物交配| 国产三级在线视频| 欧美中文日本在线观看视频| 精品久久久久久,| 国产成人欧美在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲成a人片在线一区二区| 村上凉子中文字幕在线| 青草久久国产| 麻豆久久精品国产亚洲av | 国产精品99久久99久久久不卡| 亚洲黑人精品在线| 热99国产精品久久久久久7| 国产成人精品无人区| 欧美精品亚洲一区二区| 中国美女看黄片| 女同久久另类99精品国产91| 桃红色精品国产亚洲av| 丝袜在线中文字幕| 久久香蕉激情| 波多野结衣高清无吗| 首页视频小说图片口味搜索| 无遮挡黄片免费观看| av网站在线播放免费| 啦啦啦免费观看视频1| 亚洲avbb在线观看| 丰满的人妻完整版| 天堂影院成人在线观看| 亚洲男人天堂网一区| 亚洲精品av麻豆狂野| 久久精品人人爽人人爽视色| 精品久久久精品久久久| 51午夜福利影视在线观看| 欧美日韩精品网址| 久久精品亚洲av国产电影网| 身体一侧抽搐| 精品久久久久久电影网| 色在线成人网| 校园春色视频在线观看| 亚洲精品一二三| 50天的宝宝边吃奶边哭怎么回事| 久久久水蜜桃国产精品网| 免费一级毛片在线播放高清视频 | а√天堂www在线а√下载| 黄色成人免费大全| 久久九九热精品免费| 欧美+亚洲+日韩+国产| 国产精品亚洲av一区麻豆| av在线播放免费不卡| 国产精品一区二区在线不卡| 欧美黑人欧美精品刺激| 国产激情欧美一区二区| 俄罗斯特黄特色一大片| 免费日韩欧美在线观看| 一级片免费观看大全| 天堂俺去俺来也www色官网| 麻豆国产av国片精品| 午夜免费成人在线视频| 少妇的丰满在线观看| av视频免费观看在线观看| 高潮久久久久久久久久久不卡| 国产精品日韩av在线免费观看 | 久久精品国产综合久久久| 日日干狠狠操夜夜爽| 国产精品亚洲一级av第二区| 久9热在线精品视频| 日韩精品青青久久久久久| 少妇的丰满在线观看| 久久亚洲精品不卡| 91老司机精品| 精品久久久久久久毛片微露脸| 99国产精品一区二区三区| 免费看十八禁软件|