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

    紅樹(shù)林生長(zhǎng)氣象條件適宜度評(píng)估方法研究
    ——以廣西北海沙生紅樹(shù)林為例

    2022-03-03 09:10:40陳燕麗莫偉華
    廣西林業(yè)科學(xué) 2022年1期
    關(guān)鍵詞:氣象臺(tái)站氣象條件氣象站

    陳燕麗,張 悅,錢 拴,孫 明,莫偉華

    (1.廣西科學(xué)院廣西紅樹(shù)林研究中心 廣西紅樹(shù)林保護(hù)與利用重點(diǎn)實(shí)驗(yàn)室,廣西北海 536000;2.中國(guó)農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,北京 100094;3.國(guó)家氣象中心,北京 100081;4.廣西壯族自治區(qū)氣象科學(xué)研究所,廣西南寧 530022)

    紅樹(shù)林是生長(zhǎng)在熱帶及亞熱帶沿海潮間帶的木本常綠闊葉植物群落,是最重要的藍(lán)碳資源之一[1],生態(tài)服務(wù)功能突出[2]。紅樹(shù)林處于海洋與陸地的交錯(cuò)帶,屬于典型敏感生態(tài)脆弱帶,對(duì)氣候逆境反應(yīng)極其靈敏[3],在全球氣候變化過(guò)程中具有重要的指示作用[4]。氣候變化背景下,氣候資源時(shí)空變異的加大將增加紅樹(shù)林生態(tài)系統(tǒng)的脆弱性[5]。廣西紅樹(shù)林面積居全國(guó)第二位,紅樹(shù)林是廣西重要的濕地生態(tài)系統(tǒng)和生態(tài)保護(hù)紅線重點(diǎn)區(qū)域之一,研究氣象條件對(duì)廣西紅樹(shù)林生長(zhǎng)的影響對(duì)其保護(hù)和修復(fù)有重要意義。

    紅樹(shù)林生長(zhǎng)與氣象條件密切相關(guān)。氣溫和降水是紅樹(shù)林光合作用[6]和碳交換[7]的重要限制因子;氣溫[8-9]、降水[8]、大風(fēng)[8,10-11]和日照[8]對(duì)紅樹(shù)林生長(zhǎng)發(fā)育的影響在不同時(shí)間和地域有一定差異。氣溫、降水和海表溫度在很大程度上決定紅樹(shù)林潛在適合生長(zhǎng)區(qū)域[12-14],氣溫是人工紅樹(shù)林造林的主要限制因子[15]。紅樹(shù)林具有明顯的小氣候效應(yīng),對(duì)區(qū)域小氣候可起到明顯的調(diào)節(jié)作用[16]。現(xiàn)有的紅樹(shù)林和氣象關(guān)系研究多聚焦于兩者相互作用關(guān)系,綜合氣象條件影響分析和評(píng)估方面的報(bào)道較少。

    目前,綜合氣象條件影響分析多采用氣候適宜度模型。氣候適宜度是將溫度、光照和降水等氣象因子的數(shù)量變化定量轉(zhuǎn)化為對(duì)植物生長(zhǎng)發(fā)育、產(chǎn)量形成和品質(zhì)指標(biāo)等的適宜程度[17],評(píng)估對(duì)象主要集中在農(nóng)作物(玉米(Zea mays)[18]、小麥(Triticum aesti?vum)[19]和大豆(Glycine max)[20])和經(jīng)濟(jì)作物(橡膠[21]、煙草(Nicotiana tabacum)[22])等,在自然植被中應(yīng)用的報(bào)道較少。針對(duì)自然植被,錢拴等[23]通過(guò)年際間氣象條件指數(shù)的對(duì)比,實(shí)現(xiàn)不同時(shí)段年際間棉花生長(zhǎng)氣象條件的評(píng)價(jià);Chen 等[24]采用偏最小二乘回歸法,構(gòu)建喀斯特植被增強(qiáng)型植被指數(shù)(Enhanced veg?etation index,EVI)氣象擬合模型(Climate vegetation index,CVI),通過(guò)CVI年際差值構(gòu)建喀斯特地區(qū)植被氣象條件綜合影響評(píng)價(jià)模型。

    歸一化植被指數(shù)(Normalized difference vegeta?tion index,NDVI)是用于表征植被長(zhǎng)勢(shì)最廣泛的遙感監(jiān)測(cè)指數(shù)。本研究以廣西北海典型沙生紅樹(shù)林為研究對(duì)象,利用北海紅樹(shù)林生態(tài)氣象站和NDVI觀測(cè)數(shù)據(jù)及國(guó)家氣象臺(tái)站數(shù)據(jù),研究影響紅樹(shù)林長(zhǎng)勢(shì)的關(guān)鍵氣象因子、紅樹(shù)林小氣候效應(yīng)和氣象指標(biāo)閾值,建立紅樹(shù)林綜合氣象條件適宜度評(píng)估模型,為紅樹(shù)林保護(hù)和修復(fù)工程中開(kāi)展氣象條件影響定量評(píng)估提供技術(shù)支持。

    1 材料與方法

    1.1 研究區(qū)概況

    廣西北部灣由沿海的防城港、欽州和北海3 市組成,擁有總面積全國(guó)第二、天然林面積全國(guó)第一的紅樹(shù)林及山口和北侖河口兩個(gè)紅樹(shù)林國(guó)家級(jí)自然保護(hù)區(qū)。廣西北部灣紅樹(shù)林生態(tài)區(qū)屬南亞熱帶海洋性季風(fēng)氣候;1961—2019年,年均氣溫為22.2 ~22.9 ℃,年均降水量為1 772.3 ~2 746.4 mm,年均日照時(shí)長(zhǎng)為1 515.9 ~1 973.2 h。北海紅樹(shù)林生態(tài)氣象觀測(cè)試驗(yàn)站建成于2018年底,位于廣西北海市國(guó)家海洋科技園區(qū)的小海灣(109°18′E,21°27′N),位于廣西紅樹(shù)林生態(tài)區(qū)中部。

    1.2 數(shù)據(jù)

    國(guó)家氣象臺(tái)站數(shù)據(jù)為廣西北海國(guó)家基準(zhǔn)站觀測(cè)數(shù)據(jù),包括逐日平均氣溫、降水量、相對(duì)濕度和風(fēng)速,數(shù)據(jù)時(shí)段為1961—2019年(59年)。

    北海紅樹(shù)林生態(tài)氣象站數(shù)據(jù)來(lái)源于北海紅樹(shù)林生態(tài)氣象觀測(cè)試驗(yàn)站,包括逐日平均氣溫、最高氣溫、最低氣溫、降水量、相對(duì)濕度、風(fēng)速、地溫(5 和10 cm)、總輻射和光合有效輻射,數(shù)據(jù)時(shí)段為2018年11月8日—2019年11月16日。NDVI 數(shù)據(jù)來(lái)源于北海紅樹(shù)林生態(tài)氣象觀測(cè)試驗(yàn)站通量塔搭載的NDVI探頭,數(shù)據(jù)時(shí)段與氣象衛(wèi)星的NDVI序列相同。

    衛(wèi)星遙感數(shù)據(jù):選擇NASA MODIS 陸地產(chǎn)品組根據(jù)統(tǒng)計(jì)算法開(kāi)發(fā)的MODIS 植被指數(shù)產(chǎn)品MOD13Q1,即全球250 m分辨率、16天合成的植被指數(shù)產(chǎn)品,數(shù)據(jù)版本為V006。對(duì)所獲得的MOD13Q1遙感數(shù)據(jù)集進(jìn)行子集提取、圖像鑲嵌、數(shù)據(jù)格式轉(zhuǎn)換、投影轉(zhuǎn)換及質(zhì)量檢驗(yàn)等預(yù)處理,得到質(zhì)量可靠的NDVI 數(shù)據(jù)集,16 天合成NDVI 數(shù)據(jù),進(jìn)行最大值合成,獲取2000—2019年NDVI年最大值序列。

    1.3 方法

    1.3.1 氣象因子篩選

    利用紅樹(shù)林生態(tài)氣象站數(shù)據(jù),分析表征紅樹(shù)林長(zhǎng)勢(shì)的NDVI 與累積降水量、最高氣溫、最低氣溫、平均氣溫、相對(duì)濕度和風(fēng)速等氣象因子的相關(guān)性,結(jié)合文獻(xiàn)資料確定影響紅樹(shù)林生長(zhǎng)的關(guān)鍵氣象因子。

    1.3.2 氣象數(shù)據(jù)訂正

    紅樹(shù)林具有小氣候效應(yīng)[16],國(guó)家氣象臺(tái)站數(shù)據(jù)用于紅樹(shù)林生長(zhǎng)氣象條件評(píng)價(jià)時(shí)需進(jìn)行訂正,訂正為紅樹(shù)林生態(tài)氣象站數(shù)據(jù),用于推算歷年紅樹(shù)林生長(zhǎng)氣象指標(biāo)。利用各氣象因子國(guó)家氣象臺(tái)站與生態(tài)氣象站數(shù)據(jù)的比值,將歷年國(guó)家氣象臺(tái)站數(shù)據(jù)訂正為生態(tài)氣象站數(shù)據(jù)。計(jì)算公式如下:

    式中,Xg為國(guó)家氣象臺(tái)站數(shù)據(jù);Xs為紅樹(shù)林生態(tài)氣象站數(shù)據(jù);A為某一氣象指標(biāo)兩站數(shù)據(jù)的比值;Yn為訂正后數(shù)據(jù);Yo為原始國(guó)家氣象臺(tái)站數(shù)據(jù)。

    1.3.3 氣象指標(biāo)計(jì)算

    溫度是影響紅樹(shù)植物生長(zhǎng)的關(guān)鍵因素[25]。有效積溫是植物在某個(gè)生育期或全部生育期內(nèi)有效溫度的總和,即植物在某段時(shí)間內(nèi)日平均氣溫與生物學(xué)零度之差的總和。有效積溫可反映植物生長(zhǎng)發(fā)育對(duì)熱量的需求或衡量地區(qū)熱量資源,其剔除了生物學(xué)下限溫度以下和上限溫度以上的溫度,基本上反映了植物的生育速率與溫度的線性關(guān)系。計(jì)算公式如下[26]:

    式中,Te為評(píng)估時(shí)段的有效積溫(℃);T為日平均溫度(℃);Tb為最低溫度(℃);To為最適溫度(℃);Tm為最高溫度(℃)。研究中采用的紅樹(shù)植物3 個(gè)基點(diǎn)溫度Tb、To和Tm分別為10 ℃[8]、27 ℃[6]和38 ℃[8]。

    水分是植被生長(zhǎng)不可或缺的因子,采用生育期內(nèi)的累積降水量表示。累積降水量為在生長(zhǎng)季內(nèi)每日降水量的總和[8],計(jì)算公式如下:

    式中評(píng)估時(shí)段的累積降水量(mm);Pi為逐日降水量(mm);n為統(tǒng)計(jì)天數(shù)。

    空氣濕度對(duì)植被生長(zhǎng)有重要作用。相對(duì)濕度是指空氣中水汽壓與相同溫度下飽和水汽壓的百分比,是表征大氣干燥程度的物理量。濕度指標(biāo)采用平均相對(duì)濕度表示[27],計(jì)算公式如下:

    風(fēng)可加強(qiáng)植被和空氣的熱量交換,增加土壤蒸發(fā)和植物蒸騰,影響空氣中二氧化碳等成份的擴(kuò)散與輸送,從而影響植物生長(zhǎng)。風(fēng)速指標(biāo)采用平均風(fēng)速表示[28],計(jì)算公式如下:

    1.3.4 氣象指標(biāo)閾值確定

    依據(jù)某種氣象條件出現(xiàn)越頻繁,植被受到的鍛煉越多,其適應(yīng)能力越強(qiáng)的氣候適應(yīng)性原理,確定氣象指標(biāo)的適宜、一般適宜和不適宜閾值。設(shè)樣本依次為X1,X2,X3……Xn,計(jì)算公式如下:

    式中,Xi為某氣象指標(biāo)序列;為氣象指標(biāo)平均值;Di為氣象指標(biāo)距平絕對(duì)值;n為數(shù)據(jù)樣本總個(gè)數(shù)。對(duì)Di由小到大排序,并依據(jù)其大小,篩選對(duì)應(yīng)Xi作為某一適宜度等級(jí)閾值計(jì)算序列。其中,適宜閾值為重新排序后Di的第1 ~80個(gè)百分位值對(duì)應(yīng)的數(shù)據(jù),占總樣本數(shù)的80%;一般適宜閾值為81 ~90個(gè)百分位值對(duì)應(yīng)的數(shù)據(jù);剩余數(shù)據(jù)為不適宜閾值計(jì)算序列[29]。

    1.3.5 氣象NDVI計(jì)算

    氣候變化和人類活動(dòng)是植被變化的主要驅(qū)動(dòng)力,在非人類活動(dòng)影響下,自然植被變化的主要驅(qū)動(dòng)力為氣候變化。目前,已發(fā)展出多種方法用于分離氣候變化和人類活動(dòng)對(duì)植被的影響[30],但并未形成公認(rèn)有效的區(qū)分方法[31-32]。農(nóng)氣上常用的趨勢(shì)分離法是分離氣象和人為(生產(chǎn)力水平提高)影響的常用方法。本研究將趨勢(shì)分離法用于分析自然植被,通過(guò)分析紅樹(shù)林氣象NDVI與綜合氣象條件評(píng)分的一致性,判斷紅樹(shù)林氣象適宜度評(píng)估模型。借鑒氣象產(chǎn)量分離方法[33],將表征紅樹(shù)林長(zhǎng)勢(shì)的NDVI分解為趨勢(shì)NDVI、氣象NDVI和隨機(jī)NDVI,公式為:

    式中,y為實(shí)際NDVI;yt為趨勢(shì)NDVI;yw為氣象NDVI;Δy為隨機(jī)NDVI。隨機(jī)NDVI對(duì)實(shí)際NDVI影響較小,且無(wú)法用函數(shù)的形式表達(dá)出來(lái),一般忽略不計(jì);趨勢(shì)NDVI擬合方法采用一元線性回歸法[33]。

    2 結(jié)果與分析

    2.1 紅樹(shù)林氣象條件評(píng)估因子

    分析紅樹(shù)林冠層NDVI 與氣象因子的相關(guān)性(表1)。年尺度上,NDVI 受氣溫影響極顯著(P<0.01),受相對(duì)濕度影響顯著(P<0.05)。季節(jié)尺度上,夏季NDVI 與地溫和降水量的相關(guān)性較高,秋季NDVI 受氣溫影響極顯著(P<0.01),春季NDVI 與相對(duì)濕度的相關(guān)性最高。在年、季尺度的分析基礎(chǔ)上,進(jìn)行月尺度(表1)和日尺度(圖1)的分析。月尺度上,NDVI 在1月與風(fēng)速和輻射均呈顯著相關(guān)(P<0.05),與地溫呈極顯著相關(guān)(P<0.01);在3月與降水量和地溫均呈顯著或極顯著相關(guān)(P<0.05,P<0.01);在5月與相對(duì)濕度呈極顯著相關(guān)(P<0.01);在7月與風(fēng)速呈顯著相關(guān)(P<0.05);在8月與降水量、氣溫、相對(duì)濕度和地溫均呈極顯著相關(guān)(P<0.01);在11月與相對(duì)濕度呈極顯著相關(guān)(P<0.01)。NDVI與8月(生長(zhǎng)旺盛期)氣象因子的相關(guān)性最高。日尺度上,除降水量外,NDVI 與各氣象因子均呈顯著相關(guān)(P<0.05);與日平均氣溫、最高氣溫、最低氣溫、5 cm 地溫、總輻射和光合有效輻射均呈正相關(guān),與相對(duì)濕度、風(fēng)速和10 cm 地溫均呈負(fù)相關(guān)。降水方面的相關(guān)性不顯著可能是因?yàn)槿粘叨壬辖邓疄?的時(shí)間較多。

    圖1 2018年11月— 2019年11月紅樹(shù)林NDVI與各氣象因子日尺度上的散點(diǎn)圖Fig.1 Scatter diagrams between NDVI and daily meteorological factors in mangrove from November 2018 to November 2019

    表1 2018年11月—2019年11月紅樹(shù)林NDVI與各氣象因子的相關(guān)系數(shù)Tab.1 Correlation coefficients between NDVI and meteorological factors in mangrove from November 2018 to November 2019

    氣溫、相對(duì)濕度、地溫和輻射對(duì)紅樹(shù)林NDVI 的作用在不同時(shí)間尺度(年、季和月)上均有體現(xiàn),降水量和風(fēng)速僅在月尺度有所體現(xiàn)。通過(guò)查閱文獻(xiàn)和災(zāi)情資料,廣西北部灣紅樹(shù)林區(qū)出現(xiàn)的幾次較嚴(yán)重的病蟲(chóng)害均與高溫、連續(xù)無(wú)降水日數(shù)較多的氣候條件有關(guān)[34]。紅樹(shù)林生態(tài)氣象站觀測(cè)的地溫受潮汐影響,與普通氣象臺(tái)站地溫有較大差異,因此氣象條件評(píng)估不考慮該因子;降水量和風(fēng)速在年和季尺度上與NDVI相關(guān)程度低,但在兩個(gè)關(guān)鍵月份(3和8月)中相關(guān)性高,仍為影響紅樹(shù)林生長(zhǎng)的主要?dú)庀笠蜃?;平均氣溫和最高氣溫、最低氣溫與NDVI 相關(guān)性高,前者包含后兩者的信息,因此只選平均氣溫作為影響紅樹(shù)林生長(zhǎng)的熱量因子;考慮到多數(shù)氣象站缺少總輻射、光合有效輻射的觀測(cè),為使建立的指標(biāo)和模型適用于大多數(shù)氣象站,本研究暫不考慮輻射因子。在綜合紅樹(shù)林NDVI 與氣象因子的相關(guān)性、文獻(xiàn)資料及未來(lái)模型推廣的基礎(chǔ)上,選擇平均氣溫、降水量、相對(duì)濕度和風(fēng)速作為紅樹(shù)林生長(zhǎng)氣象條件評(píng)估因子。

    2.2 國(guó)家氣象臺(tái)站和生態(tài)氣象站數(shù)據(jù)比較

    紅樹(shù)林小氣候效應(yīng)明顯,生態(tài)氣象站只代表1 個(gè)點(diǎn),為充分利用紅樹(shù)林周圍區(qū)域更多氣象站的資料,全面評(píng)估氣象條件對(duì)紅樹(shù)林的影響,進(jìn)行國(guó)家氣象臺(tái)站和生態(tài)氣象站數(shù)據(jù)比較。各氣象因子國(guó)家氣象臺(tái)站與生態(tài)氣象站數(shù)據(jù)比值(簡(jiǎn)稱“比值”)年均值為0.8 ~1.0,年最大值為0.9 ~1.3,年最小值為0.8 ~1.0(表2)。對(duì)于年均值,降水量、平均氣溫和相對(duì)濕度國(guó)家氣象臺(tái)站與生態(tài)氣象站的數(shù)據(jù)無(wú)差異(比值1.0),風(fēng)速差異較大(比值0.8)。對(duì)于年最大值,平均氣溫和相對(duì)濕度國(guó)家氣象臺(tái)站與生態(tài)氣象站的數(shù)據(jù)無(wú)差異(比值1.0),降水量差異較大(比值1.3)。對(duì)于年最小值,降水量國(guó)家氣象臺(tái)站與生態(tài)氣象站的數(shù)據(jù)無(wú)差異(比值1.0),風(fēng)速差異較大(比值0.8)。

    表2 各氣象因子國(guó)家氣象臺(tái)站和生態(tài)氣象站數(shù)據(jù)比值Tab.2 Ratios of meteorological factors observed in national meteorological station and ecological meteorological station

    2.3 紅樹(shù)林氣象條件評(píng)估和驗(yàn)證

    利用訂正后的氣象數(shù)據(jù),計(jì)算紅樹(shù)林生育期有效積溫、累積降水量、相對(duì)濕度和風(fēng)速。依據(jù)本研究的閾值劃分方法,獲得4種評(píng)估指標(biāo)適宜、一般和不適宜3個(gè)等級(jí)區(qū)間(表3)。

    表3 單一氣象因子適宜度評(píng)估區(qū)間Tab.3 Suitability evaluation intervals of single meteorological factor

    依據(jù)單一氣象因子適宜度評(píng)估閾值,分別評(píng)價(jià)有效積溫、累積降水量、相對(duì)濕度和風(fēng)速的適宜度,將單一氣象指標(biāo)適宜、一般和不適宜數(shù)量化,評(píng)分分別為10、6 和3 分。參考紅樹(shù)林NDVI 與氣象因子相關(guān)性分析結(jié)果,采用層次分析法將有效積溫、相對(duì)濕度、累積降水量和風(fēng)速的影響權(quán)重分別設(shè)置為0.5、0.25、0.15和0.10;采用加權(quán)綜合法,建立紅樹(shù)林綜合氣象條件適宜度評(píng)估模型,計(jì)算公式如下:

    式中,T、R、H和W分別為有效積溫、累積降水量、相對(duì)濕度和風(fēng)速的適宜度評(píng)分。8 ≤S≤10 時(shí),綜合氣象條件判定為適宜;4.5 ≤S<8 時(shí),綜合氣象條件判定為一般;3 ≤S<4.5 時(shí),綜合氣象條件判定為不適宜。根據(jù)綜合氣象條件適宜度評(píng)估模型,進(jìn)行2000—2019年北海紅樹(shù)林綜合氣象條件適宜度評(píng)估(表4)。近20年,北海多數(shù)年份的綜合氣象條件適宜度為適宜,一般適宜年份有兩年(2001 和2017年),不適宜年份僅有1年(2019年)。

    表4 北海紅樹(shù)林綜合氣象條件適宜度評(píng)估Tab.4 Suitability evaluation based on multiple meteorological conditions in mangrove in Beihai

    續(xù)表4 Continued

    為驗(yàn)證綜合氣象條件適宜度評(píng)估模型的有效性,本研究選取北海市面積較大的紅樹(shù)林作為代表樣點(diǎn),獲取紅樹(shù)林MODIS NDVI年最大值序列數(shù)據(jù)并分離氣象NDVI(圖2)。紅樹(shù)林長(zhǎng)勢(shì)與綜合氣象條件評(píng)分有較高的一致性;總體來(lái)說(shuō),綜合氣象條件評(píng)分高的年份,紅樹(shù)林長(zhǎng)勢(shì)較好;綜合氣象條件評(píng)分低的年份,紅樹(shù)林長(zhǎng)勢(shì)較差。典型年份為2001、2005、2010、2017 和2019年。進(jìn)一步分析得知,2000 — 2019年綜合氣象條件評(píng)分與氣象ND?VI 全序列呈顯著相關(guān),相關(guān)系數(shù)達(dá)0.63;其中,氣象NDVI 負(fù)值年份與綜合氣象條件評(píng)分相關(guān)系數(shù)達(dá)0.80。

    圖2 北海紅樹(shù)林氣象NDVI與綜合氣象條件評(píng)分變化曲線Fig.2 Change curves of meteorological NDVI and scores of multiple meteorological conditions in mangrove in Beihai

    3 討論與結(jié)論

    3.1 討論

    已有的紅樹(shù)林和氣象相關(guān)性的研究多基于不連續(xù)生長(zhǎng)發(fā)育狀況實(shí)地調(diào)查數(shù)據(jù),已證實(shí)氣溫和降水[8]是紅樹(shù)林生長(zhǎng)發(fā)育的重要影響因子。本研究依托北海紅樹(shù)林生態(tài)氣象站NDVI 高通量觀測(cè)的優(yōu)勢(shì),從月、季和年3個(gè)尺度分析多個(gè)氣象因子與表征紅樹(shù)林長(zhǎng)勢(shì)NDVI 的相關(guān)性。結(jié)果顯示,氣溫和降水量在紅樹(shù)林生長(zhǎng)旺盛時(shí)期(8月)與NDVI 均呈極顯著相關(guān),與已有研究結(jié)論一致[8-9]。本研究發(fā)現(xiàn),在年和季尺度上,氣溫與NDVI 相關(guān)性較高,降水量與NDVI 的相關(guān)性不高,風(fēng)速與NDVI 在1 和7月呈顯著相關(guān);相對(duì)濕度與NDVI 在年尺度上呈顯著相關(guān),在5、8 和11月呈極顯著相關(guān)。本研究以氣象與NDVI 的相關(guān)性選擇氣象評(píng)估指標(biāo),理論上具有客觀性,但本研究獲取的NDVI 時(shí)間序列較短,在保障生態(tài)氣象站設(shè)備正常運(yùn)轉(zhuǎn)的前提下,獲取更長(zhǎng)時(shí)間序列的NDVI 數(shù)據(jù),可進(jìn)一步驗(yàn)證氣象指標(biāo)的合理性。

    由于紅樹(shù)林群落結(jié)構(gòu)的復(fù)雜性[35],生態(tài)氣象站群落代表性有限,評(píng)估模型在其他紅樹(shù)林生態(tài)區(qū)的適用性需進(jìn)一步驗(yàn)證。由于生態(tài)氣象站氣象數(shù)據(jù)有限,本研究?jī)H用比值法對(duì)國(guó)家氣象臺(tái)站天氣數(shù)據(jù)進(jìn)行校正,在今后的研究中,應(yīng)在更長(zhǎng)時(shí)間序列數(shù)據(jù)的基礎(chǔ)上,對(duì)訂正方法進(jìn)行改進(jìn)。同時(shí),應(yīng)結(jié)合更多年限的生育期或NDVI 數(shù)據(jù),對(duì)氣象條件適宜度評(píng)估模型進(jìn)行驗(yàn)證,提高生育期閾值評(píng)價(jià)的精度。

    氣象災(zāi)害是影響紅樹(shù)林生長(zhǎng)的關(guān)鍵因子,如高溫天氣極易誘發(fā)紅樹(shù)林病蟲(chóng)害,影響紅樹(shù)林健康生長(zhǎng)[36]。本研究建立的模型評(píng)估結(jié)果在NDVI 明顯降低年份與綜合氣象條件評(píng)分吻合度高,一定程度上表明模型對(duì)氣象災(zāi)害具有敏感性。目前,氣象災(zāi)害對(duì)紅樹(shù)林的影響多為定性研究[8],如能定量分析氣象災(zāi)害對(duì)紅樹(shù)林生長(zhǎng)的影響,可提高模型評(píng)估精度,對(duì)紅樹(shù)林應(yīng)對(duì)氣候變化和紅樹(shù)林保護(hù)修復(fù)提供更科學(xué)和合理的決策參考。

    3.2 結(jié)論

    本研究通過(guò)相關(guān)分析篩選并確定氣溫、降水量、風(fēng)速和相對(duì)濕度為影響廣西北海典型沙生紅樹(shù)林生長(zhǎng)的關(guān)鍵氣象因子;通過(guò)分析紅樹(shù)林生態(tài)氣象站與其周圍氣象站數(shù)據(jù)的差異,明確了年度氣象條件定量評(píng)估需對(duì)風(fēng)速資料進(jìn)行訂正;定量評(píng)估了影響紅樹(shù)林生長(zhǎng)的全生育期有效積溫、累積降水量、相對(duì)濕度和風(fēng)速;通過(guò)提出氣象條件影響下的NDVI 且分離出年際間氣象條件波動(dòng)影響下的氣象NDVI,消除了年際間氣象條件差異對(duì)紅樹(shù)林長(zhǎng)勢(shì)產(chǎn)生的影響;利用加權(quán)綜合法建立紅樹(shù)林綜合氣象條件適宜度評(píng)估模型,氣象條件評(píng)估結(jié)果與紅樹(shù)林長(zhǎng)勢(shì)變化有較高的一致性,建立的氣象條件適宜度評(píng)估指標(biāo)和模型可用于開(kāi)展該區(qū)域紅樹(shù)林生長(zhǎng)氣象條件影響評(píng)價(jià)。

    猜你喜歡
    氣象臺(tái)站氣象條件氣象站
    珠峰上架起世界最高氣象站
    基于氣象條件的船舶引航風(fēng)險(xiǎn)等級(jí)
    心靈氣象站
    基于基層氣象臺(tái)站預(yù)報(bào)服務(wù)的幾點(diǎn)探討
    新形勢(shì)下縣級(jí)綜合氣象業(yè)務(wù)與人才培養(yǎng)管理工作初探
    政府在基層氣象臺(tái)站保護(hù)氣象探測(cè)環(huán)境中的作用分析
    氣象條件對(duì)某新型蒸發(fā)冷卻空調(diào)的影響
    自動(dòng)氣象站應(yīng)該注意的一些防雷問(wèn)題
    縣級(jí)綜合氣象業(yè)務(wù)一體化發(fā)展的思路
    2013年十三師農(nóng)業(yè)氣象條件分析
    狂野欧美激情性xxxx| 无遮挡黄片免费观看| 婷婷六月久久综合丁香| 18禁美女被吸乳视频| 在线观看日韩欧美| 免费在线观看日本一区| 精品久久久精品久久久| 91成人精品电影| 啦啦啦 在线观看视频| 91成年电影在线观看| 香蕉国产在线看| 成熟少妇高潮喷水视频| 欧美乱码精品一区二区三区| 国产成人欧美| 无遮挡黄片免费观看| 天天添夜夜摸| 女人精品久久久久毛片| 日日夜夜操网爽| 亚洲人成网站在线播放欧美日韩| 美女午夜性视频免费| 欧美日韩中文字幕国产精品一区二区三区 | 久久久久久亚洲精品国产蜜桃av| 黄色毛片三级朝国网站| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人精品久久二区二区91| 日韩一卡2卡3卡4卡2021年| 人人澡人人妻人| 精品少妇一区二区三区视频日本电影| 亚洲精品久久成人aⅴ小说| 国内毛片毛片毛片毛片毛片| а√天堂www在线а√下载| a级毛片在线看网站| 色精品久久人妻99蜜桃| 免费在线观看完整版高清| 亚洲,欧美精品.| 亚洲熟妇中文字幕五十中出 | 91大片在线观看| 看黄色毛片网站| 侵犯人妻中文字幕一二三四区| 999久久久国产精品视频| 亚洲少妇的诱惑av| 亚洲国产看品久久| 久久久水蜜桃国产精品网| 亚洲国产精品一区二区三区在线| 超色免费av| 精品久久久久久电影网| 国产精品 欧美亚洲| 久久国产亚洲av麻豆专区| 久久国产精品影院| 国产精品影院久久| 日韩欧美三级三区| 1024视频免费在线观看| 黄色a级毛片大全视频| 天堂俺去俺来也www色官网| 亚洲激情在线av| 1024视频免费在线观看| 怎么达到女性高潮| 婷婷精品国产亚洲av在线| 视频区欧美日本亚洲| 三级毛片av免费| 一区二区三区国产精品乱码| 嫩草影视91久久| 久久久久久久午夜电影 | 99re在线观看精品视频| 亚洲五月婷婷丁香| 国产成人av激情在线播放| 激情视频va一区二区三区| 久久精品国产清高在天天线| 国产欧美日韩一区二区精品| 久久 成人 亚洲| 伊人久久大香线蕉亚洲五| 成人免费观看视频高清| 国产区一区二久久| 老司机在亚洲福利影院| 午夜福利,免费看| 久久久久九九精品影院| 窝窝影院91人妻| 国产高清视频在线播放一区| 亚洲狠狠婷婷综合久久图片| 亚洲欧美精品综合久久99| 老熟妇乱子伦视频在线观看| 亚洲成人国产一区在线观看| 99在线视频只有这里精品首页| 亚洲色图综合在线观看| 国产男靠女视频免费网站| 免费高清在线观看日韩| 精品久久久久久成人av| 亚洲精品国产区一区二| 久久久国产精品麻豆| 亚洲人成伊人成综合网2020| 欧美人与性动交α欧美软件| 久久天堂一区二区三区四区| 精品熟女少妇八av免费久了| 欧美人与性动交α欧美软件| 久久精品aⅴ一区二区三区四区| 国产黄色免费在线视频| 看片在线看免费视频| 自线自在国产av| 91成年电影在线观看| 最新在线观看一区二区三区| 久久精品亚洲熟妇少妇任你| 午夜免费观看网址| 真人一进一出gif抽搐免费| 看片在线看免费视频| 亚洲少妇的诱惑av| 亚洲精华国产精华精| 亚洲精品国产色婷婷电影| 中文亚洲av片在线观看爽| 88av欧美| 亚洲成a人片在线一区二区| 日韩免费高清中文字幕av| 成人三级黄色视频| 日韩人妻精品一区2区三区| 欧美中文综合在线视频| 欧美亚洲日本最大视频资源| 国产精品九九99| 久久这里只有精品19| 精品久久久精品久久久| 国产精品香港三级国产av潘金莲| 亚洲一区二区三区欧美精品| 久久 成人 亚洲| 精品电影一区二区在线| 狠狠狠狠99中文字幕| 一级片'在线观看视频| 国产精品一区二区精品视频观看| 国产成年人精品一区二区 | 亚洲欧美日韩无卡精品| 亚洲av熟女| 成人18禁在线播放| 日本免费一区二区三区高清不卡 | 免费在线观看亚洲国产| 岛国在线观看网站| 99精品欧美一区二区三区四区| 久久久国产成人免费| 在线国产一区二区在线| 精品人妻在线不人妻| 久久天躁狠狠躁夜夜2o2o| 亚洲情色 制服丝袜| 久久人人97超碰香蕉20202| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品国产一区二区精华液| 亚洲熟女毛片儿| 亚洲精品国产精品久久久不卡| 日本vs欧美在线观看视频| 免费日韩欧美在线观看| 99在线人妻在线中文字幕| 很黄的视频免费| 国产精品乱码一区二三区的特点 | 18禁裸乳无遮挡免费网站照片 | 99热国产这里只有精品6| 亚洲午夜理论影院| 在线观看一区二区三区| 99精品久久久久人妻精品| 午夜福利影视在线免费观看| 亚洲久久久国产精品| 精品国产一区二区久久| 黄色成人免费大全| 国产免费现黄频在线看| 黄色丝袜av网址大全| 国产无遮挡羞羞视频在线观看| 国产一区二区三区在线臀色熟女 | 亚洲av成人不卡在线观看播放网| 99热只有精品国产| 亚洲性夜色夜夜综合| 免费高清在线观看日韩| 国产极品粉嫩免费观看在线| 精品日产1卡2卡| 91精品国产国语对白视频| 一边摸一边做爽爽视频免费| 在线观看一区二区三区激情| 大陆偷拍与自拍| 久9热在线精品视频| 制服人妻中文乱码| 国产精品98久久久久久宅男小说| 人人妻,人人澡人人爽秒播| 国产精品九九99| 欧美日韩一级在线毛片| 久久青草综合色| 久久久久久人人人人人| 亚洲成人国产一区在线观看| 国产精品久久视频播放| 多毛熟女@视频| 女性被躁到高潮视频| 天天躁狠狠躁夜夜躁狠狠躁| 欧美激情高清一区二区三区| 在线av久久热| 看黄色毛片网站| 精品久久久久久,| a在线观看视频网站| 国产成人免费无遮挡视频| 欧美在线黄色| 国产精品美女特级片免费视频播放器 | 老司机亚洲免费影院| 亚洲激情在线av| 很黄的视频免费| 一二三四在线观看免费中文在| 长腿黑丝高跟| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成人免费av在线播放| 久久久久九九精品影院| 日韩免费高清中文字幕av| 在线免费观看的www视频| 很黄的视频免费| 一级黄色大片毛片| 国产精品爽爽va在线观看网站 | 色精品久久人妻99蜜桃| 国产免费男女视频| 成人免费观看视频高清| 黄色丝袜av网址大全| 波多野结衣高清无吗| 一级,二级,三级黄色视频| 成人精品一区二区免费| 国产精品国产av在线观看| 午夜老司机福利片| 久久久久久免费高清国产稀缺| 午夜福利免费观看在线| 国产精品一区二区免费欧美| 村上凉子中文字幕在线| 国产成人av教育| 亚洲专区国产一区二区| 最近最新中文字幕大全免费视频| netflix在线观看网站| 黄色女人牲交| 中出人妻视频一区二区| 亚洲专区中文字幕在线| 欧美乱妇无乱码| 99久久99久久久精品蜜桃| 国产免费av片在线观看野外av| 亚洲人成伊人成综合网2020| 每晚都被弄得嗷嗷叫到高潮| 亚洲人成电影观看| 韩国av一区二区三区四区| 午夜福利,免费看| 大码成人一级视频| 久久久久亚洲av毛片大全| 在线观看午夜福利视频| 激情在线观看视频在线高清| 久久人妻福利社区极品人妻图片| 窝窝影院91人妻| 国产免费av片在线观看野外av| aaaaa片日本免费| 日本免费一区二区三区高清不卡 | 国产精品永久免费网站| 国产av在哪里看| 久久午夜亚洲精品久久| 欧美精品啪啪一区二区三区| 丝袜美腿诱惑在线| 久久久久久久久免费视频了| 身体一侧抽搐| 悠悠久久av| 亚洲专区国产一区二区| 99国产极品粉嫩在线观看| 国产欧美日韩精品亚洲av| 亚洲人成电影观看| 亚洲情色 制服丝袜| 好男人电影高清在线观看| 19禁男女啪啪无遮挡网站| 操美女的视频在线观看| 岛国在线观看网站| 99国产极品粉嫩在线观看| 亚洲av美国av| 亚洲黑人精品在线| 国产精品 欧美亚洲| 精品国产超薄肉色丝袜足j| 久久久国产精品麻豆| 桃色一区二区三区在线观看| 久热爱精品视频在线9| 国产三级黄色录像| 最新在线观看一区二区三区| 久久国产精品影院| 淫秽高清视频在线观看| 狂野欧美激情性xxxx| 亚洲色图av天堂| 深夜精品福利| e午夜精品久久久久久久| 国产精品国产高清国产av| 一二三四社区在线视频社区8| 美女 人体艺术 gogo| 欧美乱色亚洲激情| 国产精品久久久人人做人人爽| 欧美人与性动交α欧美精品济南到| 一区二区三区国产精品乱码| 国产又爽黄色视频| 免费人成视频x8x8入口观看| 搡老乐熟女国产| 日韩免费av在线播放| 男女下面插进去视频免费观看| 法律面前人人平等表现在哪些方面| 麻豆成人av在线观看| 国产精品一区二区三区四区久久 | 国产精品一区二区免费欧美| 午夜精品国产一区二区电影| 国产精品自产拍在线观看55亚洲| 性欧美人与动物交配| 99国产精品一区二区三区| 亚洲第一av免费看| 狠狠狠狠99中文字幕| 可以免费在线观看a视频的电影网站| 热re99久久国产66热| 国产深夜福利视频在线观看| 久久中文字幕人妻熟女| 男女高潮啪啪啪动态图| 亚洲av成人不卡在线观看播放网| 大型av网站在线播放| 制服人妻中文乱码| 日本五十路高清| 亚洲国产欧美一区二区综合| 高潮久久久久久久久久久不卡| 丰满饥渴人妻一区二区三| 黄色女人牲交| 日韩大码丰满熟妇| 国产1区2区3区精品| 国产av在哪里看| 精品福利永久在线观看| 精品第一国产精品| 激情视频va一区二区三区| 色综合站精品国产| 国产成人欧美在线观看| 又黄又爽又免费观看的视频| 亚洲男人天堂网一区| 欧美一级毛片孕妇| 天天影视国产精品| 大码成人一级视频| av在线播放免费不卡| 桃色一区二区三区在线观看| 亚洲成人免费av在线播放| 精品国产一区二区久久| av天堂在线播放| 亚洲人成电影免费在线| 久久人人97超碰香蕉20202| 亚洲一卡2卡3卡4卡5卡精品中文| 999精品在线视频| 国产精品成人在线| 色哟哟哟哟哟哟| 一边摸一边抽搐一进一小说| 国产精品免费视频内射| 亚洲五月天丁香| 久久久久国产精品人妻aⅴ院| 午夜福利在线免费观看网站| 国产成人精品久久二区二区免费| 黑丝袜美女国产一区| 国产成人精品久久二区二区免费| 天天添夜夜摸| 成人永久免费在线观看视频| 精品国产乱码久久久久久男人| 欧美日韩国产mv在线观看视频| 亚洲av日韩精品久久久久久密| 国产激情欧美一区二区| 国产精品九九99| 日日爽夜夜爽网站| 十八禁网站免费在线| 夜夜躁狠狠躁天天躁| 国产精品免费一区二区三区在线| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久 | 日韩三级视频一区二区三区| 欧美大码av| 天天躁狠狠躁夜夜躁狠狠躁| 男人操女人黄网站| 精品一区二区三区av网在线观看| 我的亚洲天堂| 自线自在国产av| 十八禁人妻一区二区| 欧美人与性动交α欧美软件| 亚洲av成人av| 可以免费在线观看a视频的电影网站| 国产精品免费视频内射| 嫩草影院精品99| 亚洲欧洲精品一区二区精品久久久| 好男人电影高清在线观看| 男女下面插进去视频免费观看| 丁香六月欧美| 国产精品久久久人人做人人爽| 精品久久久久久电影网| 中文亚洲av片在线观看爽| 国产有黄有色有爽视频| 最近最新中文字幕大全电影3 | 少妇裸体淫交视频免费看高清 | 美女福利国产在线| 18禁国产床啪视频网站| 桃色一区二区三区在线观看| 亚洲精品在线美女| 欧美激情高清一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲av第一区精品v没综合| 免费在线观看视频国产中文字幕亚洲| 免费人成视频x8x8入口观看| 国产视频一区二区在线看| 超色免费av| 一二三四在线观看免费中文在| 久久精品成人免费网站| 国产免费男女视频| 欧美日韩中文字幕国产精品一区二区三区 | 成人黄色视频免费在线看| 国产成+人综合+亚洲专区| 日本 av在线| bbb黄色大片| 黄网站色视频无遮挡免费观看| а√天堂www在线а√下载| 国产高清国产精品国产三级| 欧美日本亚洲视频在线播放| 视频在线观看一区二区三区| 精品卡一卡二卡四卡免费| 他把我摸到了高潮在线观看| 99久久国产精品久久久| 成人特级黄色片久久久久久久| 免费人成视频x8x8入口观看| 中文字幕高清在线视频| 99re在线观看精品视频| 欧美日韩视频精品一区| 美国免费a级毛片| 欧美午夜高清在线| 桃色一区二区三区在线观看| 欧美日韩亚洲高清精品| 久久久久九九精品影院| 黑人猛操日本美女一级片| 亚洲 国产 在线| 男女之事视频高清在线观看| 9色porny在线观看| 精品久久久久久,| 人人妻人人澡人人看| 日韩中文字幕欧美一区二区| 亚洲欧美一区二区三区黑人| 在线观看免费高清a一片| 精品一品国产午夜福利视频| 视频区图区小说| 999精品在线视频| 国产成人影院久久av| 无限看片的www在线观看| 久久精品人人爽人人爽视色| 一级毛片女人18水好多| 另类亚洲欧美激情| 熟女少妇亚洲综合色aaa.| 美女大奶头视频| 在线观看午夜福利视频| 欧美不卡视频在线免费观看 | 国产成人精品在线电影| 制服人妻中文乱码| 亚洲精品久久成人aⅴ小说| 国产单亲对白刺激| 嫩草影视91久久| 欧美色视频一区免费| 亚洲精品成人av观看孕妇| 亚洲aⅴ乱码一区二区在线播放 | 亚洲久久久国产精品| 99国产精品99久久久久| 长腿黑丝高跟| 国产成+人综合+亚洲专区| 免费搜索国产男女视频| 日韩欧美免费精品| 欧美最黄视频在线播放免费 | 窝窝影院91人妻| 极品教师在线免费播放| 国产精品一区二区精品视频观看| 后天国语完整版免费观看| 少妇粗大呻吟视频| 少妇的丰满在线观看| 国产不卡一卡二| 91大片在线观看| av有码第一页| 国产精品久久久久成人av| 国产成人精品无人区| 亚洲中文字幕日韩| 女警被强在线播放| 自线自在国产av| 91麻豆av在线| 深夜精品福利| 午夜福利影视在线免费观看| 一夜夜www| 一级a爱视频在线免费观看| 制服诱惑二区| 亚洲成人精品中文字幕电影 | 天天影视国产精品| 黑人欧美特级aaaaaa片| 国产精品自产拍在线观看55亚洲| 九色亚洲精品在线播放| 一边摸一边抽搐一进一小说| 精品国产亚洲在线| 国产成人精品在线电影| 免费日韩欧美在线观看| 大型av网站在线播放| 侵犯人妻中文字幕一二三四区| 成人手机av| 日韩免费av在线播放| 精品国产一区二区久久| 老熟妇仑乱视频hdxx| 一级片免费观看大全| 久热爱精品视频在线9| 国产成人影院久久av| 老汉色∧v一级毛片| 国产极品粉嫩免费观看在线| 91成人精品电影| 美女扒开内裤让男人捅视频| 欧美黄色片欧美黄色片| 男男h啪啪无遮挡| 男女之事视频高清在线观看| 少妇 在线观看| 波多野结衣一区麻豆| 国产深夜福利视频在线观看| 夜夜夜夜夜久久久久| 黄色 视频免费看| 99riav亚洲国产免费| 久久久久国产精品人妻aⅴ院| 嫩草影院精品99| 可以在线观看毛片的网站| 亚洲精品久久成人aⅴ小说| bbb黄色大片| 国产精品98久久久久久宅男小说| 亚洲av电影在线进入| 精品第一国产精品| av在线播放免费不卡| 黄色女人牲交| 性少妇av在线| 国产高清视频在线播放一区| 制服诱惑二区| 女同久久另类99精品国产91| а√天堂www在线а√下载| 午夜91福利影院| 久久久国产精品麻豆| 可以在线观看毛片的网站| 91精品三级在线观看| 免费av中文字幕在线| 嫩草影视91久久| 欧美日韩视频精品一区| 丰满饥渴人妻一区二区三| 久久天躁狠狠躁夜夜2o2o| 欧美激情久久久久久爽电影 | 一级a爱片免费观看的视频| 久久香蕉激情| 亚洲 欧美 日韩 在线 免费| 色精品久久人妻99蜜桃| 黄色 视频免费看| 国产亚洲欧美98| 精品久久久久久,| 高清在线国产一区| 国产亚洲欧美在线一区二区| 免费日韩欧美在线观看| 国产欧美日韩一区二区精品| 国产成人系列免费观看| 久久国产精品人妻蜜桃| 嫩草影院精品99| 水蜜桃什么品种好| 国产色视频综合| 亚洲少妇的诱惑av| av天堂久久9| 91成年电影在线观看| 国产免费男女视频| 桃红色精品国产亚洲av| 人人妻,人人澡人人爽秒播| 免费久久久久久久精品成人欧美视频| 中文字幕人妻丝袜一区二区| 美国免费a级毛片| 久久精品成人免费网站| 欧美+亚洲+日韩+国产| 亚洲avbb在线观看| 亚洲av美国av| 国产亚洲精品一区二区www| 麻豆av在线久日| 天天躁夜夜躁狠狠躁躁| 国产亚洲欧美精品永久| 欧美日本中文国产一区发布| 日韩国内少妇激情av| 中国美女看黄片| 12—13女人毛片做爰片一| 亚洲第一欧美日韩一区二区三区| 91精品三级在线观看| 高清欧美精品videossex| 久久久久九九精品影院| 国产av在哪里看| 正在播放国产对白刺激| 国产一区二区在线av高清观看| 亚洲欧美一区二区三区久久| 18美女黄网站色大片免费观看| ponron亚洲| av网站在线播放免费| 国产成年人精品一区二区 | 欧美日韩乱码在线| 不卡一级毛片| 久久草成人影院| 国产一区二区三区综合在线观看| 热99国产精品久久久久久7| 嫁个100分男人电影在线观看| 久久人人精品亚洲av| 丰满迷人的少妇在线观看| 极品教师在线免费播放| 国产成人欧美| 精品久久久久久成人av| 大陆偷拍与自拍| 999精品在线视频| √禁漫天堂资源中文www| 嫩草影院精品99| 欧美日韩黄片免| 日韩欧美在线二视频| 国产黄a三级三级三级人| 久久国产精品影院| 狂野欧美激情性xxxx| 国产精品日韩av在线免费观看 | 欧美老熟妇乱子伦牲交| www日本在线高清视频| 美女 人体艺术 gogo| 人人妻人人爽人人添夜夜欢视频| 曰老女人黄片| 夜夜夜夜夜久久久久| 亚洲精品美女久久av网站| 一级黄色大片毛片| 久久草成人影院| 久久久国产精品麻豆| 精品一区二区三区视频在线观看免费 | 91在线观看av| 欧美色视频一区免费| 在线永久观看黄色视频| 亚洲专区字幕在线| 视频区图区小说| 日本精品一区二区三区蜜桃|