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

    中國陸地生態(tài)系統(tǒng)通量觀測站點(diǎn)空間代表性

    2013-12-20 11:27:16王紹強(qiáng)陳蝶聰何洪林閆慧敏
    生態(tài)學(xué)報(bào) 2013年24期
    關(guān)鍵詞:生態(tài)區(qū)觀測站代表性

    王紹強(qiáng),陳蝶聰,2,* ,周 蕾,何洪林,石 浩,閆慧敏,蘇 文

    (1.中國科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)網(wǎng)絡(luò)觀測與模擬重點(diǎn)實(shí)驗(yàn)室,北京 100101;2.中國科學(xué)院大學(xué),北京 100049)

    渦度相關(guān)技術(shù)是測定大氣與生態(tài)系統(tǒng)之間CO2交換、水分和能量通量最直接的方法,所提供的數(shù)據(jù)越來越多地被用作區(qū)域性和全球性陸地生態(tài)系統(tǒng)和大氣間CO2、水、能量交換的分析依據(jù)[1]。全球先后建立了500多個(gè)基于渦度相關(guān)技術(shù)的通量觀測站點(diǎn),形成全球性和區(qū)域性的通量觀測網(wǎng)絡(luò)(Fluxnet)[2],包括美國通量網(wǎng)(AmeriFLUX)、歐洲通量網(wǎng)(CarboEurope)、亞洲通量網(wǎng)(AsiaFlux)、中國通量網(wǎng)(ChinaFLUX)等共42個(gè)國家、23個(gè)區(qū)域性通量研究網(wǎng)絡(luò)[3]。由于通量站點(diǎn)構(gòu)成的網(wǎng)絡(luò)本身從未經(jīng)過正式設(shè)計(jì),新增通量觀測塔的選址并沒有嚴(yán)格地落在整個(gè)網(wǎng)絡(luò)中碳水通量空間代表性較弱的區(qū)域,而只是根據(jù)建站可行性和經(jīng)驗(yàn)性判斷,并不一定能很好地監(jiān)測特定區(qū)域范圍內(nèi)所有生態(tài)系統(tǒng)類型的通量變異特征[2]。同時(shí),通量塔的測量值僅僅代表貢獻(xiàn)區(qū)尺度的通量,還需要通過尺度擴(kuò)展方法,應(yīng)用到區(qū)域、國家、洲際或全球尺度來定量化陸地生物圈和大氣間的通量凈交換量[4]。近幾年國外在渦度相關(guān)通量尺度擴(kuò)展方面取得了較大的進(jìn)展[5-6],有7個(gè)區(qū)域性通量觀測研究網(wǎng)絡(luò)對尺度擴(kuò)展進(jìn)行了研究,研究內(nèi)容包括:(1)評價(jià)通量網(wǎng)絡(luò)的代表性[2];(2)站點(diǎn)尺度的通量如何擴(kuò)展到更大的空間尺度[7-8];(3)研究區(qū)域、洲際或全球通量的幅度、分布和年際變化[7-8];(4)評估空間異質(zhì)性和參數(shù)變異性對通量估算的影響[7,9-10]。因此,通量觀測數(shù)據(jù)從站點(diǎn)到區(qū)域的尺度擴(kuò)展,需要首先評價(jià)通量觀測站點(diǎn)的空間代表性,才能用于區(qū)域尺度的生態(tài)系統(tǒng)和大氣碳水通量交換的時(shí)空格局分析。

    國外已有科學(xué)家對美國通量網(wǎng)和歐洲通量網(wǎng)的代表性開展了定量分析評價(jià)[2,11]。例如,Hargrove等[11]研究了美國通量網(wǎng)實(shí)際設(shè)計(jì)對通量觀測數(shù)據(jù)分析的影響,確定美國通量網(wǎng)現(xiàn)有渦度通量塔對整個(gè)美國大陸通量環(huán)境的代表性程度。Sulkava等[2]提出一種網(wǎng)絡(luò)設(shè)計(jì)定量化方法(Quantitative network design,QND),用來評價(jià)現(xiàn)有網(wǎng)絡(luò)的代表性,通過增加或移除觀測站點(diǎn)來優(yōu)化網(wǎng)絡(luò),合理設(shè)計(jì)使新的通量網(wǎng)絡(luò)達(dá)到最優(yōu)化。因此,結(jié)合多元統(tǒng)計(jì)理論和空間分析技術(shù)[12],使得生成的每個(gè)生態(tài)區(qū)內(nèi)的碳通量都具有相似的影響因素,從而有助于定量表達(dá)通量塔代表的空間格局。

    中國陸地生態(tài)系統(tǒng)通量觀測研究網(wǎng)絡(luò)(ChinaFLUX)于2002年建成,最初擁有8個(gè)微氣象和16個(gè)箱式/氣相色譜法觀測站,對農(nóng)田、草地、森林和水體等典型生態(tài)系統(tǒng)與大氣間CO2、水汽、能量通量進(jìn)行長期觀測研究,帶動(dòng)了我國渦度相關(guān)通量觀測研究的迅速發(fā)展[3]。近幾年來,我國科研和教育機(jī)構(gòu)相繼在國內(nèi)建立了一批渦度相關(guān)通量觀測站點(diǎn),全國通量觀測站達(dá)到85個(gè),彌補(bǔ)了ChinaFLUX觀測站在我國生態(tài)系統(tǒng)空間分布和植被類型代表性上的不足,增強(qiáng)了我國通量觀測研究的實(shí)力。但是,中國通量觀測站點(diǎn)的建設(shè)和其他國家的網(wǎng)絡(luò)建設(shè)相似,基本上是依托已有野外臺(tái)站上進(jìn)行建設(shè)的,覆蓋了主要?dú)夂騾^(qū)內(nèi)重要的土地利用類型,然而并沒有完全考慮植物、土壤、地形和環(huán)境的多樣性、干擾和管理措施等的因素,作為一個(gè)整體來監(jiān)測區(qū)域尺度生態(tài)系統(tǒng)碳通量變異缺乏系統(tǒng)的評價(jià)。

    目前中國對通量站點(diǎn)代表性的研究更多的是利用風(fēng)浪區(qū)模型[13],結(jié)合實(shí)際觀測數(shù)據(jù),對密云[14]、黑河流域[15]等通量觀測站點(diǎn)的空間代表性進(jìn)行初步分析。由于中國地形復(fù)雜,經(jīng)緯度跨度較大,生態(tài)系統(tǒng)類型多樣,景觀破碎化較嚴(yán)重,現(xiàn)有的通量觀測站點(diǎn)僅能反映中國部分生態(tài)系統(tǒng)碳通量的空間分布特征。因此,通過對現(xiàn)有觀測站點(diǎn)空間代表性程度的定量分析評價(jià),提出我國陸地生態(tài)系統(tǒng)通量站點(diǎn)的合理布局,從而使得優(yōu)化后的通量網(wǎng)絡(luò)能夠代表大部分或主要類型的生態(tài)系統(tǒng),有利于通量數(shù)據(jù)與遙感資料有效結(jié)合,提高碳水通量觀測從站點(diǎn)擴(kuò)展到區(qū)域尺度的精度。

    1 研究方法與數(shù)據(jù)處理

    1.1 數(shù)據(jù)來源

    本文收集了中國大陸和臺(tái)灣地區(qū)的渦度相關(guān)通量觀測站點(diǎn)信息(附表1),按照站點(diǎn)的不同屬性如所在省份、所屬部門和生態(tài)系統(tǒng)類型等分類統(tǒng)計(jì),繪制現(xiàn)有通量觀測站點(diǎn)的空間分布圖,為通量生態(tài)區(qū)的劃分判定提供依據(jù)。

    通量是生態(tài)系統(tǒng)碳循環(huán)中最為重要的特征參數(shù)之一,進(jìn)行陸地生態(tài)系統(tǒng)碳通量的生態(tài)區(qū)劃時(shí),主要選擇影響植物光合作用的因素作為變量[16]。渦度相關(guān)通量塔所觀測到的是凈生態(tài)系統(tǒng) CO2交換量(Net ecosystem exchange,NEE,本文中簡稱碳通量),本文收集了11個(gè)與生態(tài)系統(tǒng)碳通量最為相關(guān)的變量,數(shù)據(jù)空間分辨率1 km×1 km,時(shí)間分辨率為年,包括氣象因素、土壤因素和地形因素的非生物因子,以及實(shí)際植被狀態(tài)、植被生產(chǎn)力的生物因子。氣象資料來自中國生態(tài)系統(tǒng)研究網(wǎng)絡(luò)(CERN)數(shù)據(jù)庫,包括1970—2000年30a平均氣溫(TEM)、降水量(PRE)、太陽輻射(RAD)和年日照時(shí)數(shù)(SUN);地形因素選擇1∶100萬的數(shù)字高程模型(DEM);土壤數(shù)據(jù)包括20世紀(jì)80年代土壤普查數(shù)據(jù)得到的土壤有機(jī)碳(SOC)和土壤氮含量(STN);實(shí)際植被狀況數(shù)據(jù)包括MODIS反演得到的2006—2008年平均葉面積指數(shù)(LAI)、2005年增強(qiáng)型植被指數(shù)(EVI)以及來自中國生態(tài)系統(tǒng)研究網(wǎng)絡(luò)(CERN)數(shù)據(jù)庫的1970—2000年平均光合有效輻射(PAR);參照歐洲通量網(wǎng)代表性評價(jià)中采用總初級生產(chǎn)力(GPP)作為評價(jià)指標(biāo)[2],本文采用BEPS模型模擬得到的1990—2000年年均GPP作為植被生產(chǎn)力變量。所選變量數(shù)據(jù)大部分為多年平均值,盡管在時(shí)間尺度上無法完全統(tǒng)一,但是也能反映碳通量影響因素的平均狀況和趨勢。

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

    1.2.1 標(biāo)準(zhǔn)化處理

    所有的輸入數(shù)據(jù)都統(tǒng)一了坐標(biāo)系和投影格式(Albers投影,南標(biāo)準(zhǔn)緯度線為25°N,北標(biāo)準(zhǔn)緯度線為47°N,中央經(jīng)線為110°E),并生成空間分辨率為1 km×1 km的柵格數(shù)據(jù)(行列數(shù)為4888×4000)。由于不同的地圖圖層測量單位的差異會(huì)直接影響聚類算法的結(jié)果[11],因此在聚類運(yùn)算之前所有的變量圖層都必須進(jìn)行標(biāo)準(zhǔn)化處理。首先統(tǒng)計(jì)每個(gè)輸入圖層所有像元的最大值和最小值,每個(gè)像元減去該變量圖層的最小值后除以最大值與最小值的差值,從而使每個(gè)變量圖層的像元值都?xì)w一化為0—1之間(圖1),所有變量圖層采用統(tǒng)一的圖例,剔除各個(gè)變量數(shù)量級的差異所導(dǎo)致的聚類誤差。

    1.2.2 相關(guān)性分析

    利用ARCGIS的空間相關(guān)性分析模塊對選取的11個(gè)輸入變量進(jìn)行兩兩相關(guān)性分析。根據(jù)相關(guān)性矩陣(表1)可以看出,土壤有機(jī)碳與土壤氮含量相關(guān)系數(shù)達(dá)90%,本文沒有隨意剔除其中一個(gè),因?yàn)橛锌赡?0%的不相關(guān)中含有額外的信息能用來區(qū)分碳通量差異顯著的兩個(gè)地區(qū)[11],所以,所選取的11個(gè)變量全都保留用于聚類運(yùn)算。

    1.2.3 主成分分析

    主成分分析(PCA)是將原來眾多的變量重新組合成一組新的互相無關(guān)的綜合變量,同時(shí)根據(jù)實(shí)際需要,從中取出幾個(gè)較少的綜合變量盡可能多地反映原來變量信息的統(tǒng)計(jì)方法。對11個(gè)已作相關(guān)性分析的聚類輸入變量進(jìn)行主成分分析,計(jì)算特征值λi(i=1,2,…,11)并按從大到小順序排列,分別計(jì)算主成分貢獻(xiàn)率和累計(jì)貢獻(xiàn)率。一般來講,為了達(dá)到降維的目的,只提取前幾個(gè)主成分,由于第一、第二、第三主成分的累計(jì)貢獻(xiàn)率達(dá)到86.68%(表2),根據(jù)累計(jì)貢獻(xiàn)率大于85%的原則[17],故選取前三個(gè)主成分,用3個(gè)新變量PC1、PC2、PC3來代替原來的11個(gè)變量。

    1.3 多元空間地理聚類

    圖1 11個(gè)經(jīng)過標(biāo)準(zhǔn)化的聚類圖層Fig.1 11 standardized clustering layers

    表1 11個(gè)輸入變量的相關(guān)性矩陣Table 1 Correlation matrix of 11 input layers

    表2 特征值及主成分貢獻(xiàn)率、累計(jì)貢獻(xiàn)率Table 2 Eigen values,principal component contribution rate and cumulative contribution rate

    空間聚類分析是將地理空間實(shí)體或地理單元集合依照某種相似性度量原則劃分為若干個(gè)類似地理空間實(shí)體或地理單元組成的類或簇的過程,類中實(shí)體或單元彼此間具有較高相似性,類與類之間的實(shí)體或單元具有較大差異性,是空間數(shù)據(jù)挖掘的重要組成部分[18-19]。本文使用ENVI軟件的非監(jiān)督分類中的K-means算法,來實(shí)現(xiàn)空間聚類(圖2)。非監(jiān)督分類即在多光譜圖像中搜尋、定義其自然相似光譜集群組的過程,其中不需要人工選擇訓(xùn)練樣本,僅需極少的人工初始輸入,計(jì)算機(jī)按照一定規(guī)則自動(dòng)地根據(jù)像元光譜或空間等特征組成集群組,然后分析者將每個(gè)組和參考數(shù)據(jù)比較,將其劃分到某一類別中去,主要有K-means算法和ISODATA 算法[20]。

    K-means算法在空間聚類各算法中一直處于核心地位,是由MacQueen于1967年提出的,在目前的聚類分析中應(yīng)用最為廣泛,具有算法簡單且收斂速度快的特點(diǎn)[21]。本研究中采用K-means算法作為多元空間地理聚類的算法,把與碳通量相關(guān)的生物與非生物因子作為空間地理聚類的輸入因子,隨機(jī)在所有圖層選擇k類中心,每個(gè)聚類中心擁有11個(gè)變量,隨后剩余的每個(gè)柵格,根據(jù)其自帶的11個(gè)變量的值計(jì)算與k類聚類中心的距離,并根據(jù)最近距離原則賦值給最近的聚類中心。聚類中心及其屬于該聚類中心的柵格重新計(jì)算新的聚類中心,不斷重復(fù)這個(gè)過程,直到準(zhǔn)則函數(shù)收斂[22]。

    K-means算法具體步驟:

    (1)選K個(gè)初始聚類中心,z1(k),z2(k),…,zk(k),其中k為迭代運(yùn)算的次序號,初始時(shí)k=1。

    (2)逐個(gè)將模式樣本x按最小距離準(zhǔn)則分配給K個(gè)聚類中心中的某一個(gè)zj(k),

    (3)計(jì)算新的聚類中心,zjk+1( ),1≤j≤K,

    式中,Nj為第j個(gè)聚類域Sj中所包含的樣本個(gè)數(shù)。

    (4)設(shè)e為迭代誤差閾值,可根據(jù)需要設(shè)置,若

    則迭代結(jié)束,否則k=k+1,返回(1)。

    1.4 通量生態(tài)區(qū)的代表性評價(jià)

    圖2 多元空間地理聚類方法示意圖Fig.2 Diagram of multivariate geographic clustering approach

    生態(tài)區(qū)劃是在對生態(tài)系統(tǒng)客觀認(rèn)識(shí)和充分研究的基礎(chǔ)上,應(yīng)用生態(tài)學(xué)原理和方法,結(jié)合地理學(xué)、氣候?qū)W、土壤學(xué)、環(huán)境科學(xué)和資源科學(xué)等多個(gè)學(xué)科的知識(shí),揭示自然生態(tài)區(qū)域的相似性和差異性規(guī)律,以及人類活動(dòng)對生態(tài)系統(tǒng)干擾的規(guī)律,從而進(jìn)行整合和分異,劃分生態(tài)環(huán)境的區(qū)域單元[23-24]。傅伯杰等[25-26]根據(jù)我國的氣候、地貌、地形、生態(tài)系統(tǒng)特點(diǎn)以及人類活動(dòng)規(guī)律等特征,將我國劃分為3個(gè)生態(tài)大區(qū)、13個(gè)生態(tài)地區(qū)、57個(gè)生態(tài)區(qū)。Williams等[27]利用多元地理聚類法,從影響通量變化的氣候、地形和土壤等16個(gè)因子中挑選出最重要的3個(gè)因子,通過主成分分析法進(jìn)行聚類運(yùn)算,繪制了美國愛荷華州的農(nóng)業(yè)生態(tài)區(qū)劃分布圖。與生態(tài)區(qū)劃相類似,通量生態(tài)區(qū)劃是根據(jù)與生態(tài)系統(tǒng)碳通量相關(guān)的環(huán)境變量的綜合影響,綜合分析各個(gè)要素或變量的空間分異特征、結(jié)構(gòu)組合和區(qū)域分布對通量的影響。通量生態(tài)區(qū)內(nèi)部的氣象因素、土壤因素和地形因素的非生物因子,以及實(shí)際植被狀態(tài)、植被生產(chǎn)力的生物因子具有高度的相似性,不同的通量生態(tài)區(qū)能夠反映區(qū)域間碳通量的差異。利用多元地理聚類分析方法對每個(gè)像元進(jìn)行空間聚類,得到不同聚類數(shù)(25類、50類、75類、100類、150類和200類)的通量生態(tài)區(qū)。根據(jù)中國現(xiàn)有通量觀測站點(diǎn)的空間分布格局,與空間聚類得到的通量生態(tài)區(qū)進(jìn)行對比疊加,綜合分析現(xiàn)有通量站點(diǎn)的代表性。

    2 結(jié)果分析

    2.1 中國通量觀測站點(diǎn)的現(xiàn)狀及分析

    中國現(xiàn)有85個(gè)通量觀測站,其中包括38個(gè)森林站、17個(gè)農(nóng)田站、16個(gè)草地站、7個(gè)荒漠站和7個(gè)濕地站(圖3)。從站點(diǎn)分布圖中可以看出,目前中國通量觀測站點(diǎn)大部分分布在中國東部地區(qū);黑河流域(甘肅、青海交界處)有8個(gè)通量站點(diǎn),分布較為密集;西部高寒區(qū)和干旱區(qū)的通量站點(diǎn)較少,新疆只有1個(gè)阜康草地站和1個(gè)塔中荒漠站,西藏北部沒有通量站點(diǎn);西南地區(qū)如廣西、貴州兩省缺少通量塔的布設(shè)。從生態(tài)系統(tǒng)類型統(tǒng)計(jì)來看,森林站點(diǎn)最多,占整個(gè)中國通量站點(diǎn)的45%,而荒漠、濕地站點(diǎn)較少。

    2.2 多元空間聚類通量生態(tài)區(qū)分析

    本文利用多元空間地理聚類方法,分別對影響碳通量的3個(gè)主成分變量圖層進(jìn)行25類、50類、75類、100類、150類和200類的空間聚類運(yùn)算,結(jié)果如圖4所示。從圖中可以看出,西北地區(qū)由于生態(tài)系統(tǒng)類型單一,環(huán)境因子相似,通量生態(tài)區(qū)劃的界限比較明顯,不同分類數(shù)量所導(dǎo)致的區(qū)劃間差異較小;南方地區(qū)生態(tài)系統(tǒng)類型錯(cuò)綜復(fù)雜,隨著聚類個(gè)數(shù)的增加,通量生態(tài)區(qū)更加破碎。因此南方地區(qū)比西北地區(qū)可能需要布設(shè)更多的通量觀測站點(diǎn),來反映復(fù)雜地形條件下生態(tài)系統(tǒng)類型的通量變化特征。

    為了與25、50、75、100、150、200 類的空間聚類運(yùn)算結(jié)果進(jìn)行對比分析,按照現(xiàn)有85個(gè)通量觀測站點(diǎn),本文還進(jìn)行了85類通量生態(tài)區(qū)的計(jì)算。隨著通量生態(tài)區(qū)數(shù)量的增加,對生態(tài)區(qū)的劃分不斷細(xì)化,通量生態(tài)區(qū)面積的平均值、最大值、最小值、標(biāo)準(zhǔn)差呈不同幅度的下降趨勢(表3)。隨著聚類數(shù)的增加,原有的比較均質(zhì)的區(qū)域由于環(huán)境因子的微小差異,被劃分成更加細(xì)致的區(qū)域(圖4)。

    圖3 中國陸地生態(tài)系統(tǒng)通量觀測站點(diǎn)的生態(tài)系統(tǒng)類型Fig.3 Ecosystem types of eddy flux tower stations of terrestrial ecosystems in China

    圖4 不同聚類數(shù)的通量生態(tài)區(qū)分布,顏色代表不同的生態(tài)區(qū)Fig.4 Flux-ecoregions based on different number of clusters

    將25類通量生態(tài)區(qū)的平均面積分別除以25、50、75、100、150、200類聚類的通量生態(tài)區(qū)相應(yīng)的平均面積,發(fā)現(xiàn)通量生態(tài)區(qū)平均面積縮小的倍數(shù)與聚類數(shù)增加的倍數(shù)一致(表3)。以表3中25類通量生態(tài)區(qū)面積的平均值、最大值、最小值、極差和標(biāo)準(zhǔn)差分別作為基準(zhǔn),通量生態(tài)區(qū)面積的平均值和標(biāo)準(zhǔn)差的縮小倍數(shù)與聚類數(shù)增加的倍數(shù)較為一致,最大值的縮小倍數(shù)落后于聚類數(shù)增加的倍數(shù),而最小值的縮小倍數(shù)大于聚類數(shù)增加的倍數(shù)。隨著通量生態(tài)區(qū)數(shù)量的增加,最小面積減小的速度明顯要快于最大面積的減少,這主要是南方地區(qū)環(huán)境因子相對復(fù)雜,分類更加細(xì)化,而北方地區(qū)環(huán)境因子單一,差異相對較小。

    表3 不同聚類數(shù)的通量生態(tài)區(qū)面積Table 3 Areas of different flux-ecoregions

    2.3 現(xiàn)有站點(diǎn)與通量生態(tài)區(qū)的綜合分析

    基于現(xiàn)有通量觀測站點(diǎn)的數(shù)量,對3個(gè)主成分變量圖層進(jìn)行了85類的聚類運(yùn)算(圖5)。將現(xiàn)有通量觀測站點(diǎn)與空間聚類得到的85類通量生態(tài)區(qū)疊加分析,在85個(gè)通量生態(tài)區(qū)中,有51個(gè)區(qū)含有通量觀測站點(diǎn),生態(tài)區(qū)覆蓋率即含有通量觀測站點(diǎn)的生態(tài)區(qū)的個(gè)數(shù)占分區(qū)個(gè)數(shù)(即聚類數(shù))的比例為60%,國土面積覆蓋率即含有通量觀測站點(diǎn)的生態(tài)區(qū)的面積占全國國土面積的比例為63.2%。對85類通量生態(tài)區(qū)中的51個(gè)含有通量觀測點(diǎn)的分區(qū)作生態(tài)區(qū)內(nèi)所含站點(diǎn)數(shù)統(tǒng)計(jì):31個(gè)區(qū)內(nèi)只含有一個(gè)站點(diǎn),12個(gè)區(qū)內(nèi)含有2個(gè)站點(diǎn),6個(gè)區(qū)內(nèi)含有3個(gè)站點(diǎn),2個(gè)區(qū)內(nèi)含有4個(gè)以上的站點(diǎn)。因此,假設(shè)全國有85類生態(tài)區(qū),目前這85個(gè)站點(diǎn)僅能代表全國60%左右的生態(tài)系統(tǒng)碳通量空間特征,還有近40%的區(qū)域需要增加通量觀測站點(diǎn)的布設(shè),即新增34個(gè)通量觀測站點(diǎn)。

    對不同聚類數(shù)的通量生態(tài)區(qū)中現(xiàn)有站點(diǎn)的覆蓋情況進(jìn)行統(tǒng)計(jì)(表4),發(fā)現(xiàn)隨著分區(qū)個(gè)數(shù)的增加,含有通量觀測站點(diǎn)的生態(tài)區(qū)的個(gè)數(shù)呈增加趨勢,但是現(xiàn)有站點(diǎn)的生態(tài)區(qū)覆蓋率和國土面積覆蓋率都在逐漸減小。85類通量生態(tài)區(qū)的國土面積覆蓋率比75類的高,是由于85類通量生態(tài)區(qū)的空間聚類是考慮了現(xiàn)有站點(diǎn)的空間分布,導(dǎo)致大部分站點(diǎn)落在面積較大的生態(tài)區(qū),沒有站點(diǎn)的生態(tài)區(qū)數(shù)量多而且面積較小,因而按照國土覆蓋率達(dá)到90%以上來說比100類通量生態(tài)區(qū)需要更多地新增站點(diǎn)。從150個(gè)通量生態(tài)區(qū)增加到200個(gè)時(shí),含有通量站點(diǎn)的生態(tài)區(qū)個(gè)數(shù)和含有1個(gè)站點(diǎn)的生態(tài)區(qū)個(gè)數(shù)都不變,意味著聚類區(qū)的增加不再對現(xiàn)有85個(gè)通量站的空間分布有明顯的區(qū)劃作用了,可以將150作為當(dāng)前通量站點(diǎn)數(shù)量增加的上限。

    現(xiàn)有的85個(gè)通量站點(diǎn)在50類和75類聚類的通量生態(tài)區(qū)中,僅能覆蓋60%—67%的國土面積,說明目前85個(gè)站點(diǎn)的分布仍然不平衡,部分地區(qū)較為集中,空間代表性需要加強(qiáng)。如果按照85類通量生態(tài)區(qū)的區(qū)劃方案,生態(tài)區(qū)數(shù)量和國土面積覆蓋率都達(dá)到90%以上的,僅僅新增加25—26個(gè)通量站點(diǎn)即可。其次,150類通量生態(tài)區(qū)中有62個(gè)區(qū)包含現(xiàn)有站點(diǎn),在此基礎(chǔ)上,至少需要增加73個(gè)站點(diǎn)才能覆蓋90%以上的通量生態(tài)區(qū)。同樣150類通量生態(tài)區(qū)中,現(xiàn)有站點(diǎn)只能覆蓋國土面積的43.5%,若要覆蓋率達(dá)到90%以上,需要新增加55個(gè)站點(diǎn),從覆蓋更多國土面積的角度考慮,能夠少增加約20個(gè)站點(diǎn),可以節(jié)省建設(shè)、運(yùn)行等費(fèi)用。

    圖5 85類通量生態(tài)區(qū)劃分Fig.5 85 clusters of flux-ecoregion regionalization

    表4 不同聚類數(shù)通量生態(tài)區(qū)內(nèi)的現(xiàn)有站點(diǎn)覆蓋情況Table 4 The covering of existing sites in different flux-ecoregions

    2.4 現(xiàn)有站點(diǎn)與自然地理區(qū)劃的綜合分析

    為了更好地分析中國通量觀測站點(diǎn)空間代表性,我們將現(xiàn)有通量觀測站點(diǎn)與自然地理區(qū)劃[28]疊加分析(圖6)。從圖6中發(fā)現(xiàn),全國55個(gè)自然地理區(qū)劃生態(tài)區(qū)中有18個(gè)生態(tài)區(qū)缺乏通量觀測站點(diǎn),主要集中在中溫帶半干旱地區(qū)、南亞熱帶濕潤地區(qū)、青藏高原寒帶干旱地區(qū)和青藏高原溫帶干旱地區(qū)。從這個(gè)意義上來講,在18個(gè)缺乏通量觀測站點(diǎn)的生態(tài)區(qū)分別增加通量站點(diǎn),也能起到增強(qiáng)空間代表性的作用。根據(jù)現(xiàn)有通量觀測站點(diǎn)與自然地理區(qū)劃生態(tài)區(qū)的疊加分析,發(fā)現(xiàn)在37個(gè)含有通量觀測點(diǎn)的自然地理區(qū)劃生態(tài)區(qū)中:13個(gè)生態(tài)區(qū)只含有1個(gè)站點(diǎn),9個(gè)生態(tài)區(qū)含有2個(gè)站點(diǎn),7個(gè)生態(tài)區(qū)含有3個(gè)站點(diǎn),4個(gè)生態(tài)區(qū)含有4個(gè)站點(diǎn),3個(gè)生態(tài)區(qū)含有5個(gè)以上的站點(diǎn)。通量觀測站分布的空間變異性較大,85個(gè)通量站點(diǎn)中,38個(gè)站點(diǎn)分布在生態(tài)區(qū)邊界,大多數(shù)站點(diǎn)位于較偏地區(qū)和接近邊界地區(qū),只有8個(gè)站點(diǎn)處于生態(tài)區(qū)中心附近,由此看來,目前通量觀測站點(diǎn)的空間代表性需要加強(qiáng)。

    圖6 通量站點(diǎn)在自然地理區(qū)劃生態(tài)區(qū)中的分布Fig.6 Spatial distribution of eddy flux tower stations in geographical regionalization

    3 討論

    國外已有科學(xué)家對美國通量網(wǎng)和歐洲通量網(wǎng)的代表性開展了定量分析評價(jià)。Hargrove等[11]首先選取調(diào)控碳通量以及影響植被光合作用與呼吸作用的氣候、土壤、干擾等25個(gè)環(huán)境因子,采用K-means聚類方法,生成9套不同的美國通量生態(tài)區(qū)劃,為美國通量網(wǎng)找到5個(gè)最優(yōu)的新增地理位置。Sulkava等[2]以歐洲通量網(wǎng)絡(luò)為例,根據(jù)三類研究問題設(shè)計(jì)8個(gè)模擬情景,將氣候要素和土壤要素或總初級生產(chǎn)力(Gross primary productivity,GPP)作為輸入變量,根據(jù)K-means++聚類分析,使同一類中數(shù)據(jù)點(diǎn)到聚類中心的平均距離平方即量化誤差最小,得到k個(gè)聚類,選取與k個(gè)聚類中心距離最近的像元作為k個(gè)樣本通量塔的位置。同時(shí)以GPP為目標(biāo)變量,根據(jù)k個(gè)樣本通量塔的數(shù)據(jù),用人工神經(jīng)網(wǎng)絡(luò)模型(Artificial neural network,ANN)進(jìn)行尺度擴(kuò)展,用GPP的均值、空間變異性、半方差函數(shù)參數(shù)(塊金值、基臺(tái)值、變程)和時(shí)空變異性這四項(xiàng)指標(biāo)來估算現(xiàn)有網(wǎng)絡(luò)的代表性。最后根據(jù)不同的設(shè)計(jì)情景,設(shè)定平均量化誤差閾值,確定現(xiàn)有站點(diǎn)中應(yīng)保留的站點(diǎn)數(shù)量及新增站點(diǎn)數(shù)量,計(jì)算尺度擴(kuò)展結(jié)果的不確定性,優(yōu)化通量網(wǎng)絡(luò)中3種植被類型的設(shè)計(jì)[2]。由此可見,多元統(tǒng)計(jì)理論和空間分析技術(shù)結(jié)合有助于優(yōu)化通量網(wǎng)絡(luò)的設(shè)計(jì),推動(dòng)通量觀測的尺度擴(kuò)展研究的進(jìn)一步深入發(fā)展。

    本研究采用多元地理變量空間聚類方法,基于分層的環(huán)境數(shù)據(jù)空間,把主要影響通量的氣候、土壤、植被特征相似的劃分為一個(gè)區(qū)域,把中國生態(tài)系統(tǒng)劃分成具有相對均質(zhì)通量特征的生態(tài)區(qū)域,確定現(xiàn)有渦度通量觀測站點(diǎn)對整個(gè)中國陸地生態(tài)系統(tǒng)的代表性程度,并根據(jù)與碳通量相關(guān)的生態(tài)系統(tǒng)特征定量化的相似性,為外推現(xiàn)有的通量觀測來估計(jì)那些沒有觀測值的生態(tài)區(qū)通量提供理論基礎(chǔ)。假定每一個(gè)通量生態(tài)區(qū)表示一類生態(tài)系統(tǒng),每一類生態(tài)系統(tǒng)至少有一個(gè)通量觀測塔,那么劃分的通量生態(tài)區(qū)越多,需要新增加的通量塔就越多。由于通量塔的建設(shè)受到觀測技術(shù)的發(fā)展程度、資金、運(yùn)行管理、電力、地形、環(huán)境等各方面條件的限制,要在目前85個(gè)站點(diǎn)基礎(chǔ)上新增加100多個(gè)站點(diǎn),從而覆蓋200個(gè)通量生態(tài)區(qū)也是不現(xiàn)實(shí)的。綜合現(xiàn)有站點(diǎn)對不同聚類數(shù)的通量生態(tài)區(qū)的覆蓋率和區(qū)內(nèi)站點(diǎn)數(shù)來看,通量生態(tài)區(qū)劃分為100—150類比較合適,需要新增的通量塔數(shù)量大約25—55個(gè),這樣既覆蓋中國生態(tài)系統(tǒng)的大部分主要類型,同時(shí)也能充分利用現(xiàn)有通量觀測站點(diǎn)網(wǎng)絡(luò)的優(yōu)勢。

    根據(jù)本文3.3節(jié)的結(jié)果,若要使國土面積覆蓋率達(dá)90%以上,對于85類通量生態(tài)區(qū),可在準(zhǔn)格爾盆地、柴達(dá)木盆地、青藏高原、阿拉善高原、呼倫貝爾高原、大興安嶺、黑龍江東北部、三江平原、云貴高原、四川盆地、東南地區(qū)等增加25個(gè)站點(diǎn);對于100類通量生態(tài)區(qū),可在準(zhǔn)格爾盆地、塔里木盆地、柴達(dá)木盆地、青藏高原、阿拉善高原、鄂爾多斯高原、大興安嶺、黑龍江東北部、三江平原、東北平原、云貴高原、四川盆地、東南地區(qū)等地區(qū)增加21個(gè)站點(diǎn);對于150類通量生態(tài)區(qū),可在全國大部分地區(qū)增加55個(gè)站點(diǎn)。

    綜合中國通量觀測站點(diǎn)的現(xiàn)狀和通量生態(tài)區(qū)的結(jié)果分析,西部高寒區(qū)和干旱區(qū)的通量站點(diǎn)較少,西南地區(qū)廣西、貴州兩省缺少通量塔的布設(shè),加之南方地區(qū)生態(tài)系統(tǒng)類型破碎復(fù)雜,未來的通量觀測站點(diǎn)布局可考慮在上述地區(qū)增加新的渦度相關(guān)通量站點(diǎn),提高中國生態(tài)系統(tǒng)通量觀測站點(diǎn)的代表性。

    4 結(jié)論

    本文通過收集主要影響碳通量的變量數(shù)據(jù)(包括氣象因素、土壤因素和地形因素的非生物因子、實(shí)際植被狀態(tài)以及植被生產(chǎn)力),用主成分分析法得到3個(gè)主成分變量,基于K-means空間聚類分析方法,分別計(jì)算出25、50、75、85、100、150、200類的通量生態(tài)區(qū)。利用中國現(xiàn)有的85個(gè)渦度相關(guān)通量觀測站點(diǎn),與75—200類的通量生態(tài)區(qū)進(jìn)行對比分析,能代表30%—60%的生態(tài)系統(tǒng)類型。綜合現(xiàn)有站點(diǎn)對不同聚類數(shù)的通量生態(tài)區(qū)的覆蓋率和區(qū)內(nèi)站點(diǎn)數(shù)來看,考慮現(xiàn)有站點(diǎn)分布格局,通量生態(tài)區(qū)劃分為100—150類比較合適;考慮到渦度相關(guān)設(shè)備運(yùn)行成本,在盡可能利用現(xiàn)有通量觀測塔的基礎(chǔ)上,通量站點(diǎn)可增加至100—150個(gè),這樣能覆蓋中國陸地生態(tài)系統(tǒng)的大部分主要類型。

    本研究采用3個(gè)主成分變量進(jìn)行通量生態(tài)區(qū)的空間聚類,由于信息量的減少會(huì)對聚類的結(jié)果和分析有一定的影響。其次,進(jìn)行K-means聚類運(yùn)算時(shí),像元變化閾值、最多迭代次數(shù)和聚類數(shù)量的設(shè)置主觀性較強(qiáng),也會(huì)影響到分類結(jié)果。第三,本文重點(diǎn)是用空間聚類方法獲得的通量生態(tài)區(qū),來評價(jià)現(xiàn)有站點(diǎn)在通量生態(tài)區(qū)的空間分布情況,尚未分析聚類中心和現(xiàn)有通量塔位置之間的關(guān)系,難以定量化通量站點(diǎn)的最優(yōu)位置。第四,由于難以獲取所有85個(gè)通量站點(diǎn)的月或年通量觀測值,本文無法對比分析通量生態(tài)區(qū)的GPP模擬值與站點(diǎn)的GPP觀測值,從而未能更加有效地定量評價(jià)新增站點(diǎn)在現(xiàn)有網(wǎng)絡(luò)中的代表性,無法確定現(xiàn)有站點(diǎn)中應(yīng)保留的站點(diǎn)數(shù)量及新增站點(diǎn)數(shù)量。由于現(xiàn)有通量觀測站點(diǎn)空間代表性有待提高,在現(xiàn)有站點(diǎn)空間格局的基礎(chǔ)上,結(jié)合中國地形的破碎化和生態(tài)系統(tǒng)類型的多樣性,需要找出更多的定量化的指標(biāo),提出我國陸地生態(tài)系統(tǒng)通量站點(diǎn)的合理布局,從而為通量觀測從站點(diǎn)擴(kuò)展到區(qū)域提供基礎(chǔ)數(shù)據(jù),使通量數(shù)據(jù)與遙感觀測資料有機(jī)地結(jié)合,更加有效地用于檢驗(yàn)過程機(jī)理模型的模擬結(jié)果。

    [1]Baldocchi D D,F(xiàn)alge E,Gu L H,Olson R,Hollinger D,Running S,Anthoni P,Bernhofer C,Davis K,Evans R,F(xiàn)uentes J,Goldstein A,Katul G,Law B,Lee X H,Malhi Y,Meyers T,Munger W,Oechel W,Paw K T,Pilegaard K,Schmid H P,Valentini R,Verma S,Vesala T,Wilson K,Wofsy S.FLUXNET:A new tool to study the temporal and spatial variability of ecosystem-Scale carbon dioxide,water vapor,and energy flux densities.American Meteorological Society,2001,82(11):2415-2434.

    [2]Sulkava M,Luyssaert S,Zaehle S,Papale D.Assessing and improving the representativeness of monitoring networks:The European flux tower network example.Journal of Geophysical Research,2011,116:G00J04,doi:10.1029/2010JG001562.

    [3]Yu G R,Sun X M.Principles of Flux Measurement in Terrestrial Ecosystems.Beijing:Higher Education Press,2006.

    [4]Xiao J F,Zhuang Q L,Baldocchi D D,Law B E,Richardson A D,Chen JQ,Oren R,Starr G,Noormets A,Ma SY,Verma SB,Wharton S,Wofsy SC,Bolstad P V,Burns SP,Cook D R,Curtis P S,Drake B G,F(xiàn)alk M,F(xiàn)ischer M L,F(xiàn)oster D R,Gu L H,Hadley JL,Hollinger D Y,Katul G G,Litvak M,Martin T A,Matamala R,McNulty S,Meyers T P,Monson R K,Munger JW,Oechel W C,Paw U K T,Schmid H P,Scott R L,Sun G,Suyker A E,Torn M S.Estimation of net ecosystem carbon exchange for the conterminous United States by combining MODIS and Ameri Flux data.Agricultural and Forest Meteorology,2008,148(11):1827-1847.

    [5]Xiao J F,Zhuang Q L,Law B E,Baldocchi D D,Chen J Q,Richardson A D,Melillo J M,Davis K J,Hollinger D Y,Wharton S,Oren R,Noormets A,F(xiàn)ischer M L,Verma SB,Cook D R,Sun G,McNulty S,Wofsy SC,Bolstad P V,Burns SP,Curtis P S,Drake B G,F(xiàn)alk M,F(xiàn)oster D R,Gu L H,Hadley JL,Katulk GG,Litvak M,Ma SY,Martinz T A,Matamala R,Meyers T P,Monson R K,Munger JW,Oechel W C,Paw U K T,Schmid H P,Scott R L,Starr G,Suyker A E,Torn M S.Assessing net ecosystem carbon exchange of U.S.terrestrial ecosystems by integrating eddy covariance flux measurements and satellite observations.Agricultural and Forest Meteorologys,2011,151(1):60-69.

    [6]Jung M,Reichstein M,Bondeau A.Towards global empirical upscaling of FLUXNET eddy covariance observations:Validation of a model tree ensemble approach using a biosphere model.Biogeosciences,2009,6(10):2001-2013.

    [7]Xiao J F,Davis K J,Urban N M,Keller K,Saliendra N Z.Upscaling carbon fluxes from towers to the regional scale:Influence of parameter variability and land cover representation on regional flux estimates.Journal of Geophysical Research,2011,116:G00J06,doi:10.1029/2010JG001568.

    [8]Jung M,Reichstein M,Margolis H A,Cescatti A,Richardson A D,Arain M A,Arneth A,Bernhofer C,Bonal D,Chen J Q,Gianelle D,Gobron N,Kiely G,Kutsch W,Lasslop G,Law B E,Lindroth A,Merbold L,Montagnani L,Moors E J,Papale D,Sottocornola M,Vaccari F,Williams C.Global patterns of land-atmosphere fluxes of carbon dioxide,latent heat,and sensible heat derived from eddy covariance,satellite,and meteorological observations.Journal of Geophysical Research,2011,116:G00J07,doi:10.1029/2010JG001566.

    [9]Desai A R.Climatic and phenological controls on coherent regional interannual variability of carbon dioxide flux in a heterogeneous landscape.Journal of Geophysical Research,2010,115:G00J02,doi:10.1029/2010JG001423.

    [10]Chasmer L,Kljun N,Hopkinson C,Brown S,Milne T,Giroux K,Barr A,Devito K,Creed I,Petrone R.Characterizing vegetation structural and topographic characteristics sampled by eddy covariance within two mature aspen stands using lidar and a flux footprint model:Scaling to MODIS.Journal of Geophysical Research,2011,116:G02026,doi:10.1029/2010JG001567.

    [11]Hargrove W W,Hoffman F M,Law B E.New analysis reveals representativeness of the AmeriFlux network.EOS,Transactions,American Geophysical Union,2003,84(48):529-544.

    [12]Hargrove W W,Hoffman F M.Potential of multivariate quantitative methods for delineation and visualization of ecoregions.Environmental Management,2005,34(S1):S39-S60.

    [13]Mi N,Yu G R,Wang P X,Wen X F,Sun X M.A preliminary study for spatial representiveness of flux observation at ChinaFLUX sites.Science in China Series D:Earth Sciences,2006,49(S2)24-35.

    [14]Gong L J,Liu SM.Flux spatial representativeness of Miyun//Chinese Meteorological Society.Proceeding of Ecological Meteorology Operation Establishment and Agricultural Meteorology Disaster Warning.Beijing:Chinese Meteorological Society,2007.

    [15]Shuang X,Liu S M,Xu Z W,Wang W Z.Investigation of spatial representativeness for surface flux measurements in the Heihe river basin.Advances in Earth Science,2009,24(7):724-733.

    [16]Song T,Wang SQ,Gao JX,Han Y W.A study on the scaling-up of CO2fluxes at Haibei station based on GISand RStechnologies.International Journal of Intelligent Information Technology Application,2010,3(2):71-78.

    [17]Zhang W T.Senior Tutorial of SPSSStatistical Analysis.Beijing:Higher Education Press,2004.

    [18]Chen SQ,Yang CC.Clusteringmethod for area geographical entities based on cluster validity function.Journal of Zhengzhou Institute of Surveying and Mapping,2006,23(1):44-47.

    [19]Liu S,Ji GL.A review of researches on spatial clustering.Journal of Nanjing Normal University:Engineering and Technology,2010,10(2):57-62.

    [20]Zhao Y S.Principles and Methods of Remote Sensing Application Analysis.Beijing:Science Press,2003.

    [21]Wang Y,Song Z B,Wu P L.A study on spatial clustering of urban function partition.Areal Research and Development,2009,28(1):27-31.

    [22]Lin D Y,Liu H P.Spatial clustering analysis in the point data distribution.Journal of Beijing Normal University:Natural Science,2006,42(4):419-423.

    [23]Liu G H,F(xiàn)u B J.The principle and characteristics of ecological regionalization.Advances in Environmental Science,1998,6(6):67-72.

    [24]Cheng Y Q,Zhang P Y.Progress on eco-geographical regionalization researches.Acta Ecologica Sinica,2006,26(10):3424-3433.

    [25]Fu B J,Chen L D,Liu G H.The objectives,tasks and characteristics of China ecological regionalization.Acta Ecologica Sinica,1999,19(5):591-595.

    [26]Fu B J,Liu G H,Chen L D,Ma K M,Li JR.Scheme of ecological regionalization in China.Acta Ecologica Sinica,2001,21(1):1-6.

    [27]Williams C L,Hargrove W W,Liebman M,James D E.Agro-ecoregionalization of Iowa using multivariate geographical clustering.Agriculture,Ecosystems and Environment,2008,123(1/3):161-174.

    [28]Huang B W.Primary scheme of comprehensive natural geographical zonation.Acta Geographica Sinica,1958,24(4):348-365.

    參考文獻(xiàn):

    [3]于貴瑞,孫曉敏.陸地生態(tài)系統(tǒng)通量觀測的原理與方法.北京:高等教育出版社,2006.

    [13]米娜,于貴瑞,溫學(xué)發(fā),孫曉敏.中國通量觀測網(wǎng)絡(luò)(ChinaFLUX)通量觀測空間代表性初步研究.中國科學(xué)(D輯:地球科學(xué)),2006,(增刊I):22-33.

    [14]宮麗娟,劉紹民.密云通量觀測的空間代表性 //中國氣象學(xué)會(huì).中國氣象學(xué)會(huì)2007年年會(huì)生態(tài)氣象業(yè)務(wù)建設(shè)與農(nóng)業(yè)氣象災(zāi)害預(yù)警分會(huì)場論文集.北京:中國氣象學(xué)會(huì),2007.

    [15]雙喜,劉紹民,徐自為,王維真.黑河流域觀測通量的空間代表性研究.地球科學(xué)進(jìn)展,2009,24(7):724-733.

    [17]張文彤.SPSS統(tǒng)計(jì)分析高級教程.北京:高等教育出版社,2004.

    [18]陳四清,楊春成.基于聚類有效性函數(shù)的面狀地理實(shí)體聚類.測繪科學(xué)技術(shù)學(xué)報(bào),2006,23(1):44-47.

    [19]柳盛,吉根林.空間聚類技術(shù)研究綜述.南京師范大學(xué)學(xué)報(bào):工程技術(shù)版,2010,10(2):57-62.

    [20]趙英時(shí).遙感應(yīng)用分析原理與方法.北京:科學(xué)出版社,2003.

    [21]王艷,宋振柏,吳佩林.城市功能分區(qū)的空間聚類方法研究及其應(yīng)用——以濟(jì)南市為例.地域研究與開發(fā),2009,28(1):27-31.

    [22]林冬云,劉慧平.應(yīng)用空間聚類進(jìn)行點(diǎn)數(shù)據(jù)分布研究.北京師范大學(xué)學(xué)報(bào):自然科學(xué)版,2006,42(4):419-423.[23]劉國華,傅伯杰.生態(tài)區(qū)劃的原則及其特征.環(huán)境科學(xué)進(jìn)展,1998,6(6):67-72.

    [24]程葉青,張平宇.生態(tài)地理區(qū)劃研究進(jìn)展.生態(tài)學(xué)報(bào),2006,26(10):3424-3433.

    [25]傅伯杰,陳利頂,劉國華.中國生態(tài)區(qū)劃的目的、任務(wù)及特點(diǎn).生態(tài)學(xué)報(bào),1999,19(5):591-595.

    [26]傅伯杰,劉國華,陳利頂,馬克明,李俊然.中國生態(tài)區(qū)劃方案.生態(tài)學(xué)報(bào),2001,21(1):1-6.

    [28]黃秉維.中國綜合自然區(qū)劃的初步草案.地理學(xué)報(bào),1958,24(4):348-365.

    附表1 中國陸地生態(tài)系統(tǒng)通量觀測站點(diǎn)信息Attached table 1 Information of eddy flux tower stations of terrestrial ecosystems in China

    續(xù)表

    猜你喜歡
    生態(tài)區(qū)觀測站代表性
    國家級非遺項(xiàng)目代表性傳承人簡介
    GPS導(dǎo)航對抗數(shù)據(jù)質(zhì)量特征實(shí)例分析
    四川省甘孜州:航拍四川稻城高海拔宇宙線觀測站
    漳州市非物質(zhì)文化遺產(chǎn)代表性項(xiàng)目代表性傳承人名錄
    閩臺(tái)地區(qū)代表性道地藥材
    大都市生態(tài)區(qū)綜合規(guī)劃思路研究
    江西建材(2018年4期)2018-04-10 12:36:58
    去中心化時(shí)差頻差直接定位方法
    非遺代表性傳承人
    ——勉沖·羅布斯達(dá)
    打造社區(qū)O2O生態(tài)區(qū) 彩生活住宅新模式問世
    世界
    亚洲一区中文字幕在线| 国产成人a∨麻豆精品| 菩萨蛮人人尽说江南好唐韦庄| 99国产综合亚洲精品| 日韩精品免费视频一区二区三区| av天堂久久9| 国产男女超爽视频在线观看| 另类精品久久| 亚洲美女黄色视频免费看| 天天影视国产精品| 99九九在线精品视频| 一本—道久久a久久精品蜜桃钙片| 性色av乱码一区二区三区2| 91麻豆精品激情在线观看国产 | 宅男免费午夜| 老司机靠b影院| 欧美日韩视频高清一区二区三区二| 久久99精品国语久久久| kizo精华| 亚洲成人免费电影在线观看 | 建设人人有责人人尽责人人享有的| 亚洲情色 制服丝袜| 亚洲人成电影观看| 在线天堂中文资源库| 老熟女久久久| 一区二区三区乱码不卡18| 男女边吃奶边做爰视频| 最新的欧美精品一区二区| 99re6热这里在线精品视频| 亚洲国产成人一精品久久久| 亚洲中文日韩欧美视频| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区久久| 宅男免费午夜| 伊人久久大香线蕉亚洲五| 一区二区三区四区激情视频| 精品久久久久久电影网| 老鸭窝网址在线观看| 性高湖久久久久久久久免费观看| 国产日韩欧美在线精品| 一二三四在线观看免费中文在| 一级毛片电影观看| 免费不卡黄色视频| 日韩av在线免费看完整版不卡| 国产在线免费精品| 80岁老熟妇乱子伦牲交| 99久久人妻综合| 国产精品99久久99久久久不卡| 国产日韩一区二区三区精品不卡| 午夜老司机福利片| 亚洲七黄色美女视频| 国产精品秋霞免费鲁丝片| 欧美日韩综合久久久久久| 人妻 亚洲 视频| 久久人人97超碰香蕉20202| 777久久人妻少妇嫩草av网站| 亚洲,一卡二卡三卡| 狂野欧美激情性xxxx| 青草久久国产| 亚洲中文av在线| 十八禁高潮呻吟视频| 亚洲欧美精品综合一区二区三区| 久久久久国产一级毛片高清牌| 亚洲一区二区三区欧美精品| 国产在线视频一区二区| 午夜福利视频在线观看免费| 亚洲成人免费电影在线观看 | 伊人久久大香线蕉亚洲五| 精品国产国语对白av| 成年女人毛片免费观看观看9 | 老司机影院成人| 久久午夜综合久久蜜桃| 又大又黄又爽视频免费| 成年美女黄网站色视频大全免费| 欧美变态另类bdsm刘玥| 嫩草影视91久久| 久久久精品免费免费高清| 中文字幕人妻丝袜一区二区| 精品免费久久久久久久清纯 | 99久久人妻综合| 欧美精品啪啪一区二区三区 | 国产有黄有色有爽视频| 国产日韩欧美视频二区| netflix在线观看网站| 少妇人妻 视频| 九草在线视频观看| 精品福利观看| 国产免费一区二区三区四区乱码| 制服人妻中文乱码| 亚洲免费av在线视频| 亚洲av欧美aⅴ国产| 久久女婷五月综合色啪小说| 99国产精品免费福利视频| 欧美在线黄色| www.自偷自拍.com| 国产精品一二三区在线看| 国产片特级美女逼逼视频| 美女福利国产在线| 国产精品一区二区免费欧美 | 国产成人精品久久久久久| 免费看不卡的av| 涩涩av久久男人的天堂| 黄网站色视频无遮挡免费观看| 高清欧美精品videossex| 黄色a级毛片大全视频| 精品国产一区二区久久| 日本午夜av视频| 交换朋友夫妻互换小说| 色网站视频免费| 1024香蕉在线观看| 咕卡用的链子| 亚洲欧美一区二区三区久久| 日韩av在线免费看完整版不卡| 欧美国产精品va在线观看不卡| 欧美日韩亚洲国产一区二区在线观看 | 十八禁高潮呻吟视频| 99久久99久久久精品蜜桃| 国产日韩一区二区三区精品不卡| 18禁黄网站禁片午夜丰满| 视频在线观看一区二区三区| 人体艺术视频欧美日本| 中文字幕人妻丝袜一区二区| 一区二区日韩欧美中文字幕| 国产野战对白在线观看| 国产日韩欧美亚洲二区| 亚洲专区中文字幕在线| 亚洲熟女毛片儿| 国产精品久久久久久精品电影小说| 黄网站色视频无遮挡免费观看| 啦啦啦在线免费观看视频4| 国产片内射在线| 午夜久久久在线观看| 亚洲精品国产一区二区精华液| 精品久久久久久电影网| 黄色片一级片一级黄色片| 中文字幕人妻熟女乱码| 免费观看a级毛片全部| 大香蕉久久成人网| 日韩av不卡免费在线播放| 欧美日韩综合久久久久久| 一本一本久久a久久精品综合妖精| 新久久久久国产一级毛片| 亚洲国产欧美日韩在线播放| a 毛片基地| 最近手机中文字幕大全| 少妇猛男粗大的猛烈进出视频| 久久热在线av| 欧美黄色淫秽网站| 国产极品粉嫩免费观看在线| 欧美日韩视频精品一区| 又大又爽又粗| 久久久久久久久久久久大奶| 中文字幕人妻丝袜制服| 日韩熟女老妇一区二区性免费视频| 欧美黄色片欧美黄色片| 久久久欧美国产精品| 国产av一区二区精品久久| 搡老岳熟女国产| 在线观看一区二区三区激情| 亚洲av在线观看美女高潮| 亚洲图色成人| 黄色 视频免费看| 老汉色∧v一级毛片| 久久久久久免费高清国产稀缺| 激情视频va一区二区三区| 日韩一区二区三区影片| 在线亚洲精品国产二区图片欧美| 中文字幕制服av| 韩国精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 韩国精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 91字幕亚洲| 在线亚洲精品国产二区图片欧美| 国产深夜福利视频在线观看| 亚洲av成人不卡在线观看播放网 | 亚洲精品国产av成人精品| 99国产精品99久久久久| 亚洲,一卡二卡三卡| 91精品伊人久久大香线蕉| 午夜日韩欧美国产| 亚洲自偷自拍图片 自拍| 亚洲av片天天在线观看| 少妇猛男粗大的猛烈进出视频| 99久久综合免费| 一区二区三区乱码不卡18| 一级毛片我不卡| 亚洲自偷自拍图片 自拍| 我要看黄色一级片免费的| 中文乱码字字幕精品一区二区三区| 在线 av 中文字幕| 国产亚洲av高清不卡| 久久久久久免费高清国产稀缺| 蜜桃国产av成人99| av天堂在线播放| 日韩电影二区| 国产精品二区激情视频| 亚洲视频免费观看视频| 国产伦理片在线播放av一区| 久久精品国产亚洲av涩爱| 成人亚洲精品一区在线观看| 国产精品av久久久久免费| 免费女性裸体啪啪无遮挡网站| 亚洲伊人久久精品综合| 一级片'在线观看视频| 黄片播放在线免费| 黄色毛片三级朝国网站| 国产精品熟女久久久久浪| 777久久人妻少妇嫩草av网站| 九色亚洲精品在线播放| 女人被躁到高潮嗷嗷叫费观| 乱人伦中国视频| 国产又爽黄色视频| 国产一卡二卡三卡精品| 大香蕉久久网| 国产成人精品久久二区二区91| 亚洲美女黄色视频免费看| 又大又黄又爽视频免费| 成人国语在线视频| 国产精品三级大全| 黄色毛片三级朝国网站| 国产片特级美女逼逼视频| 国产精品一二三区在线看| 国产深夜福利视频在线观看| 极品少妇高潮喷水抽搐| 一本一本久久a久久精品综合妖精| 性高湖久久久久久久久免费观看| 欧美激情高清一区二区三区| 欧美精品人与动牲交sv欧美| 爱豆传媒免费全集在线观看| 一边亲一边摸免费视频| 久久国产精品人妻蜜桃| 精品少妇内射三级| 亚洲国产精品成人久久小说| 蜜桃国产av成人99| 精品亚洲成国产av| 久久久欧美国产精品| 一本一本久久a久久精品综合妖精| 国产欧美亚洲国产| 成在线人永久免费视频| 日本a在线网址| 一区二区三区精品91| 久久久精品国产亚洲av高清涩受| 777米奇影视久久| 777米奇影视久久| 亚洲情色 制服丝袜| 女人爽到高潮嗷嗷叫在线视频| 777久久人妻少妇嫩草av网站| 国产一级毛片在线| 午夜影院在线不卡| 久久av网站| 99精国产麻豆久久婷婷| 亚洲精品自拍成人| 午夜老司机福利片| 国产精品欧美亚洲77777| 色视频在线一区二区三区| 国语对白做爰xxxⅹ性视频网站| 每晚都被弄得嗷嗷叫到高潮| 日本五十路高清| 久久毛片免费看一区二区三区| 视频区欧美日本亚洲| 伊人亚洲综合成人网| 日本一区二区免费在线视频| 欧美中文综合在线视频| 啦啦啦在线免费观看视频4| 高清欧美精品videossex| 日韩中文字幕欧美一区二区 | 人体艺术视频欧美日本| 免费看十八禁软件| 菩萨蛮人人尽说江南好唐韦庄| 精品久久蜜臀av无| av福利片在线| 18在线观看网站| 亚洲九九香蕉| 成年人免费黄色播放视频| 日韩 欧美 亚洲 中文字幕| 欧美日韩黄片免| 不卡av一区二区三区| 亚洲精品久久午夜乱码| 国产激情久久老熟女| 欧美日韩福利视频一区二区| 天天操日日干夜夜撸| 九草在线视频观看| 成年动漫av网址| 亚洲国产精品国产精品| 97精品久久久久久久久久精品| 亚洲中文日韩欧美视频| 亚洲av电影在线观看一区二区三区| 一本色道久久久久久精品综合| 后天国语完整版免费观看| 久久人妻熟女aⅴ| 成在线人永久免费视频| 久久久久精品人妻al黑| 国产午夜精品一二区理论片| 99国产综合亚洲精品| 男人操女人黄网站| 熟女少妇亚洲综合色aaa.| 国精品久久久久久国模美| 一级a爱视频在线免费观看| 精品人妻1区二区| 国产成人一区二区在线| cao死你这个sao货| 黄色片一级片一级黄色片| 成人18禁高潮啪啪吃奶动态图| 中文字幕高清在线视频| 麻豆乱淫一区二区| 搡老岳熟女国产| 国产精品一区二区在线观看99| 中文欧美无线码| 美女扒开内裤让男人捅视频| 啦啦啦在线免费观看视频4| 国产精品免费视频内射| videos熟女内射| 久久久国产精品麻豆| 91国产中文字幕| 美女扒开内裤让男人捅视频| 欧美乱码精品一区二区三区| 精品国产超薄肉色丝袜足j| 亚洲国产看品久久| 国产精品九九99| 香蕉国产在线看| xxx大片免费视频| 大片免费播放器 马上看| 国产精品久久久久久精品古装| 日韩电影二区| 中文字幕人妻丝袜制服| 又大又黄又爽视频免费| 50天的宝宝边吃奶边哭怎么回事| 日韩 亚洲 欧美在线| 日本91视频免费播放| 国产成人一区二区在线| 交换朋友夫妻互换小说| 免费不卡黄色视频| 精品少妇久久久久久888优播| 亚洲精品一卡2卡三卡4卡5卡 | 在线天堂中文资源库| 丝袜美腿诱惑在线| 久久性视频一级片| www.精华液| 亚洲自偷自拍图片 自拍| 亚洲国产成人一精品久久久| 日本欧美视频一区| 久久女婷五月综合色啪小说| 欧美日韩成人在线一区二区| 成人国语在线视频| 久久99热这里只频精品6学生| 亚洲成人免费电影在线观看 | 亚洲色图 男人天堂 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 热re99久久国产66热| 欧美日韩av久久| 亚洲 国产 在线| 国产精品二区激情视频| 日韩av不卡免费在线播放| 亚洲欧洲国产日韩| 久久国产精品人妻蜜桃| 亚洲专区国产一区二区| 亚洲少妇的诱惑av| 国产精品.久久久| 中文字幕人妻熟女乱码| 91麻豆av在线| 精品少妇内射三级| 日韩 亚洲 欧美在线| 亚洲中文字幕日韩| 99热国产这里只有精品6| 亚洲专区国产一区二区| 久久精品国产亚洲av高清一级| 欧美日韩国产mv在线观看视频| 国产欧美日韩综合在线一区二区| 国产亚洲欧美精品永久| 久久 成人 亚洲| 免费看av在线观看网站| 日韩 欧美 亚洲 中文字幕| 性高湖久久久久久久久免费观看| 制服诱惑二区| 亚洲精品国产色婷婷电影| 色婷婷久久久亚洲欧美| 国产精品久久久av美女十八| 久久女婷五月综合色啪小说| 久久精品亚洲av国产电影网| 只有这里有精品99| 高清av免费在线| 久久精品亚洲熟妇少妇任你| 国产亚洲精品第一综合不卡| 亚洲av日韩精品久久久久久密 | 最新在线观看一区二区三区 | 日本欧美国产在线视频| 午夜免费成人在线视频| 女人久久www免费人成看片| 欧美日韩一级在线毛片| 久久综合国产亚洲精品| 99国产精品一区二区三区| 国产日韩欧美亚洲二区| 精品国产乱码久久久久久男人| 99热网站在线观看| 王馨瑶露胸无遮挡在线观看| 成年av动漫网址| 黄色视频不卡| 国产国语露脸激情在线看| 99国产精品99久久久久| 成人国产一区最新在线观看 | 大片电影免费在线观看免费| 免费在线观看视频国产中文字幕亚洲 | 爱豆传媒免费全集在线观看| 咕卡用的链子| 在线观看人妻少妇| 天天添夜夜摸| 亚洲免费av在线视频| 9191精品国产免费久久| 亚洲图色成人| 国产熟女欧美一区二区| 亚洲国产av影院在线观看| 啦啦啦视频在线资源免费观看| 亚洲国产成人一精品久久久| 亚洲中文av在线| 国产欧美日韩综合在线一区二区| 国精品久久久久久国模美| av在线app专区| 国产熟女欧美一区二区| 中文字幕色久视频| 手机成人av网站| 欧美变态另类bdsm刘玥| 制服诱惑二区| 女性被躁到高潮视频| 18禁国产床啪视频网站| 欧美亚洲 丝袜 人妻 在线| 制服人妻中文乱码| 亚洲av综合色区一区| 精品人妻熟女毛片av久久网站| 日本午夜av视频| 一区二区三区四区激情视频| 纯流量卡能插随身wifi吗| 五月开心婷婷网| 国产麻豆69| 日本欧美国产在线视频| 欧美黑人精品巨大| 一区二区三区精品91| 亚洲 欧美一区二区三区| 在线观看人妻少妇| 国产成人影院久久av| 亚洲av电影在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 三上悠亚av全集在线观看| 狠狠精品人妻久久久久久综合| 性高湖久久久久久久久免费观看| 性少妇av在线| 国产成人a∨麻豆精品| 成年av动漫网址| 成人亚洲精品一区在线观看| 777米奇影视久久| 在线 av 中文字幕| 欧美黑人精品巨大| 一区二区三区精品91| 青青草视频在线视频观看| 丝袜在线中文字幕| 夫妻性生交免费视频一级片| 中文字幕色久视频| 手机成人av网站| 欧美老熟妇乱子伦牲交| 国产精品成人在线| netflix在线观看网站| 天天影视国产精品| 丝袜人妻中文字幕| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产看品久久| 免费看十八禁软件| 中文字幕高清在线视频| 国产一区有黄有色的免费视频| 久久久亚洲精品成人影院| 久久久久久久精品精品| 精品亚洲乱码少妇综合久久| av网站免费在线观看视频| 久久久久久久精品精品| 亚洲图色成人| 国产精品 国内视频| 男男h啪啪无遮挡| 日韩视频在线欧美| 嫁个100分男人电影在线观看 | 国产精品麻豆人妻色哟哟久久| 国产黄色免费在线视频| 热re99久久国产66热| 成年动漫av网址| 成人黄色视频免费在线看| 99国产精品一区二区三区| 免费看不卡的av| 亚洲欧洲国产日韩| 水蜜桃什么品种好| 捣出白浆h1v1| 女人精品久久久久毛片| 亚洲中文字幕日韩| 亚洲九九香蕉| 久久综合国产亚洲精品| 午夜免费成人在线视频| 一级,二级,三级黄色视频| 在线观看免费高清a一片| 国产又爽黄色视频| 午夜av观看不卡| 成年人午夜在线观看视频| 亚洲,欧美精品.| 七月丁香在线播放| 丰满饥渴人妻一区二区三| 国产精品av久久久久免费| 日本猛色少妇xxxxx猛交久久| 久久人人97超碰香蕉20202| 男的添女的下面高潮视频| 热re99久久国产66热| 国产欧美日韩一区二区三 | 欧美黄色片欧美黄色片| 国产成人欧美在线观看 | 国产精品一二三区在线看| 免费高清在线观看日韩| 国产人伦9x9x在线观看| 亚洲伊人色综图| 亚洲精品久久久久久婷婷小说| av电影中文网址| 亚洲成av片中文字幕在线观看| 日韩中文字幕视频在线看片| 国产成人欧美| 亚洲av欧美aⅴ国产| 高清不卡的av网站| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 韩国高清视频一区二区三区| 免费女性裸体啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 欧美av亚洲av综合av国产av| 亚洲天堂av无毛| 九色亚洲精品在线播放| 看十八女毛片水多多多| 午夜久久久在线观看| 国产爽快片一区二区三区| 国产成人影院久久av| 久久精品国产综合久久久| 黄片小视频在线播放| 中文精品一卡2卡3卡4更新| 欧美日韩亚洲高清精品| a级毛片在线看网站| 中文字幕精品免费在线观看视频| 欧美激情 高清一区二区三区| 精品欧美一区二区三区在线| 亚洲 国产 在线| 欧美日韩成人在线一区二区| 手机成人av网站| 久久国产精品大桥未久av| 50天的宝宝边吃奶边哭怎么回事| 国产av一区二区精品久久| 亚洲av男天堂| 国产高清视频在线播放一区 | 女人久久www免费人成看片| 黑人猛操日本美女一级片| 午夜影院在线不卡| 国产国语露脸激情在线看| 中文字幕亚洲精品专区| 亚洲天堂av无毛| 999精品在线视频| 麻豆国产av国片精品| 黄网站色视频无遮挡免费观看| 国产精品久久久久久人妻精品电影 | 亚洲色图综合在线观看| 赤兔流量卡办理| 高清av免费在线| 女人高潮潮喷娇喘18禁视频| 后天国语完整版免费观看| 亚洲欧美日韩另类电影网站| 叶爱在线成人免费视频播放| 国产成人免费无遮挡视频| 高清视频免费观看一区二区| 男女高潮啪啪啪动态图| 在线观看免费视频网站a站| 亚洲欧美成人综合另类久久久| 久久久精品区二区三区| 1024视频免费在线观看| 精品国产国语对白av| 午夜福利免费观看在线| 久久精品久久精品一区二区三区| 精品少妇黑人巨大在线播放| 多毛熟女@视频| 国产真人三级小视频在线观看| 久久久久久亚洲精品国产蜜桃av| 男人舔女人的私密视频| 最黄视频免费看| 丰满饥渴人妻一区二区三| 免费高清在线观看视频在线观看| 国产精品久久久久久精品电影小说| 狠狠婷婷综合久久久久久88av| 午夜日韩欧美国产| 美女高潮到喷水免费观看| av国产精品久久久久影院| 人人妻人人爽人人添夜夜欢视频| 精品一区二区三区av网在线观看 | 欧美少妇被猛烈插入视频| 欧美97在线视频| 精品久久蜜臀av无| 在线观看一区二区三区激情| 七月丁香在线播放| 丝袜喷水一区| 一级毛片 在线播放| videos熟女内射| 国产一区二区在线观看av| 欧美+亚洲+日韩+国产| 国产成人av教育| 丝瓜视频免费看黄片| 日韩制服骚丝袜av| 日韩制服丝袜自拍偷拍| 国产成人91sexporn| 男人舔女人的私密视频| 亚洲av成人精品一二三区| 亚洲专区中文字幕在线| bbb黄色大片| 午夜精品国产一区二区电影| 欧美日韩一级在线毛片| 亚洲图色成人| 国产一区二区在线观看av| 女人高潮潮喷娇喘18禁视频| 菩萨蛮人人尽说江南好唐韦庄|