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

    內蒙古草原自然資源大區(qū)植被凈初級生產力時空變化及其影響因子分析

    2022-03-27 08:49:50張子凡熊茂秋李福杰劉曉煌郝玉恒邢莉圓王新華袁鵬程
    草業(yè)科學 2022年12期
    關鍵詞:大區(qū)回歸系數(shù)時空

    張子凡,熊茂秋,李福杰,劉曉煌,郝玉恒,邢莉圓,王新華,賴 明,袁鵬程

    (1.中國地質大學(武漢)地理與信息工程學院, 湖北 武漢 430074;2.自然資源要素耦合過程與效應重點實驗室, 北京 100055;3.中國地質大學(武漢)地理與資源學院, 湖北 武漢 430074;4.中國地質調查局烏魯木齊自然資源綜合調查中心, 新疆 烏魯木齊 830000;5.中國地質調查局自然資源綜合調查指揮中心, 北京 100055;6.中國地質調查局煙臺海岸帶地質調查中心, 山東 煙臺 264000)

    內蒙古地區(qū)是我國主要的草地資源分布區(qū),也是我國的畜牧業(yè)區(qū)和糧食的主要生產地[1]。植被凈初級生產力(net primary productivity, NPP)是指綠色植物在單位時間、單位面積內所累積的有機干物質總量[2]。NPP 是生態(tài)系統(tǒng)碳匯、碳源的主要參數(shù),是陸地生態(tài)系統(tǒng)循環(huán)中必不可少的內容,也是評估生態(tài)系統(tǒng)質量與健康狀況的重要依據(jù)[3]。目前研究植被凈初級生產力的手段有實測法和模型估算法[4],實測法主要通過抽樣測量實地單位面積和時間內的生物量進行NPP 的調查和計算。隨著遙感技術的發(fā)展,利用遙感觀測數(shù)據(jù)估算區(qū)域內,長時間序列的NPP 值時空變化成為主流的研究趨勢[5]?;谶b感數(shù)據(jù)結合模型估計NPP 的方法主要有CASA模型、BEPS 模型和BIOME-BGC 模型等[6-8]。通過長時間的遙感數(shù)據(jù),結合統(tǒng)計學上的方法,對區(qū)域內NPP 值的變化趨勢和變化可持續(xù)性研究已有較多[9-12]。研究NPP 長時間序列時空變化及其影響因子,以及影響因子的時空分異性,對研究環(huán)境的變化以及不同區(qū)域碳匯的影響因素有重要意義。

    目前,利用區(qū)域的遙感觀測數(shù)據(jù),估算NPP 值,結合氣候、環(huán)境等要素,研究NPP 值變化對氣候和人類活動的響應對區(qū)域內生態(tài)環(huán)境監(jiān)測的預測有重要意義,王耠熠等[13]通過計算NPP 柵格像元點在長時間序列上的變化,利用相關性和偏相關性分析方法,研究了若爾蓋高原區(qū)NPP 值與降水量和溫度之間的關系,同時重點研究了最高溫和最低溫對NPP 的 響 應。Yin 等[14]利 用Google Earth Engine 平臺研究了森林NPP 對植被類型和氣候的響應,得出降水對NPP 的影響較強,郭睿妍等[15]利用GEE 遙感數(shù)據(jù)云平臺,采取嶺回歸分析和冗余分析的方式,研究了黃河流域森林NPP 變化及其影響因子。閆妍等[16]依據(jù)研究區(qū)內不同流域之間影響因素的不同以及CASA 模型的性質,通過從流域和省域尺度,來構建NPP 值與氣候相關要素的多元回歸模型,從時間上以月尺度來研究得出湖南省地區(qū)溫度和水汽壓對NPP 的影響較大,而降水相對較小,表明了不同區(qū)域氣候因子對NPP 值的影響不同。隨著機器學習技術的發(fā)展,在統(tǒng)計學的角度上,根據(jù)實測生物量數(shù)據(jù)和實測NPP 值結合各個因子等數(shù)據(jù)采取機器學習算法中的各類訓練的方法,對區(qū)域內NPP 和生物量進行回歸預測,可以以決策樹的方式進行非線性預測,同時可以測定各因子數(shù)據(jù)的貢獻率[17-18]。然而這些研究NPP 變化的影響因子大多是以區(qū)域每個柵格點的長時間序列的NPP 和氣候數(shù)據(jù),使用統(tǒng)計學的方法,進行相關性,偏相關性,以及整個區(qū)域的最小二乘法、嶺回歸模型和機器學習中的決策樹算法進行回歸分析。進行時間序列分析時,容易忽略對于時間變化不顯著但地理變化顯著的因子,比如地形等因子。所以只通過時間序列數(shù)據(jù)變化進行回歸分析時存在局限性,而以整個區(qū)域的柵格點為對象進行全局回歸分析,會存在空間局部不平穩(wěn)性的問題,無法反映出NPP 分布和影響因素的空間異質性。由于位于地理接近的區(qū)域,物體空間分布存在一定的自相關性和聚集性,內蒙古NPP 區(qū)域自相關性顯著[19],研究區(qū)域內NPP 時空變化的驅動因子時,進行區(qū)域全局分析多元線性回歸OLS 模型以及采用決策樹模型進行分線性回歸分析存在局部聚集等局部不平穩(wěn)性問題。而時空地理加權回歸模型,通過核函數(shù)確定時空距離的權重矩陣,可以較好地反映影響因子時空非平穩(wěn)性的問題,可以從時間和空間的角度研究驅動因素[20-21],這樣可以分析不同區(qū)域內各因子不同的回歸系數(shù)。

    本文研究基于2000-2019 年內蒙古草原自然資源大區(qū)的MOD17A3-NPP 數(shù)據(jù)、landsat7,8 的NDVI數(shù)據(jù)和國家氣象網(wǎng)站的氣象站點數(shù)據(jù)。采用統(tǒng)計分析的方法,研究該區(qū)域的NPP 時空變化及其可持續(xù)性。利用時空地理加權回歸模型,研究不同時空尺度下NPP 值的驅動因子,得出不同時空尺度的回歸系數(shù),量化不同區(qū)域NPP 變化的影響因素,同時也可為NPP 的預測提供幫助。

    1 研究區(qū)概況

    研究區(qū)域為全國自然資源綜合區(qū)劃結果中的華北草耕自然資源大區(qū)的北部區(qū)域[22-23],命名為內蒙古草原資源大區(qū),該大區(qū)是我國草地資源分布的主要區(qū)域(37°37′12″~51°02′42″ N, 106°28′40″~121°28′03″ E),主要包括河套平原和內蒙古高原、黃土高原的部分地區(qū),位于我國內蒙古自治區(qū),面積約6.51 × 108hm2,主要的自然資源為草原,約占63.42% (圖1)。

    圖1 研究區(qū)概況圖Figure 1 Overview of the study area

    2 數(shù)據(jù)源與研究方法

    2.1 數(shù)據(jù)源以及預處理

    本研究主要使用了NPP 數(shù)據(jù)、氣象數(shù)據(jù)、NDVI數(shù)據(jù)、高程數(shù)據(jù)和土地利用數(shù)據(jù)。NPP 數(shù)據(jù)來源于2000-2019 年的MODIS 陸地4 級標準數(shù)據(jù)NPP 產品數(shù)據(jù)集MOD17A3-NPP 數(shù)據(jù),來自美國航空航天局網(wǎng)站(https://ladsweb.modaps.eosdis.nasa.gov/),其空間分辨率為500 m,坐標系為WGS_1984 坐標系,投影為UTM 投影。氣象數(shù)據(jù)來自于中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn),選取了內蒙古、陜西、山西和河北省123 個氣象站點2000-2019 年的降水量、平均氣溫、相對濕度、光照時長數(shù)據(jù),采用克里金插值法進行插值,得到與NPP 數(shù)據(jù)投影坐標系相同和分辨率相同的年尺度降水、平均氣溫、相對濕度、光照時長的柵格數(shù)據(jù)。NDVI 數(shù)據(jù)來自美國陸地衛(wèi)星landsat7,8 數(shù) 據(jù) 的 反 演(https://landsat.gsfc.nasa.gov/)。高程數(shù)據(jù)和土地利用數(shù)據(jù)來自中科院資源環(huán)境數(shù)據(jù)中心(https://www.resdc.cn/),選取1 km 分辨率的土地數(shù)據(jù)和dem 數(shù)據(jù)。

    2.2 研究方法

    2.2.1 Theil-Sen Median 方法和Mann-Kendall 顯著性檢驗

    通過Theil-Sen Median 方法植被凈初級生產力的趨勢變化,利用Mann-Kendall 方法進行變化的顯著性檢驗[24]。該方法可以減少外部帶來的誤差,同時對趨勢的判斷受異常值的影響較小,有非常高的適用性,可以更準確反映長時間序列下植被凈初級生產力的變化。通過這個方法研究區(qū)域內每一個像元的時間序列變化,確定2000-2019 年NPP 的時空變化特征。Theil-Sen Median 方法公式如下:

    式中:Z為NPP變化的顯著性,n為研究的時間序列,tp為每年重復的NPP值重復的個數(shù)。定義α為顯著性檢驗水平,Z(1-α)/2為標準正態(tài)方差,當|Z| >Z(1-α)/2時,表示時間序列ET 在α水平上存在顯著變化。

    2.2.2 Hurst 指數(shù)

    Hurst 指數(shù)是用來評估時間序列信息長時間依賴的方法[25]。在本研究中,用Hurst 指數(shù)來研究植被凈初級生產力變化的可持續(xù)性,可以判斷未來變化的趨勢。具體公式和計算步驟如下:

    定義研究區(qū)間NPP 變化的時間序列:

    對Hurst 指數(shù)進擬合:

    Hurst 指數(shù)(H)在0.5 <H< 1 時,表明時間序列是一個可持續(xù)序列,即未來變化趨勢與現(xiàn)在變化趨勢一致,H越接近1,持續(xù)性變化越強。H= 0.5 則說明NPP時間序列為隨機序列,不具備時間相關性,0 <H< 0.5 則表明未來變化趨勢與現(xiàn)在有相反的趨勢,H越接近0,反持續(xù)性變化越強。

    2.2.3 時空地理加權回歸(GTWR)模型

    基于最小二乘法的多元線性回歸模型可反映全局的變量因子與因變量的關系,但對于局部聚集性比較強的數(shù)據(jù),利用全局分析無法反映因子的局部不平穩(wěn)特征。地理加權回歸模型可以反映局部空間的不平穩(wěn)性,時空地理加權回歸模型是在地理加權回歸模型基礎上,引入了時間維度[26],與地理加權回歸模型相比,能更好地揭示時空非平穩(wěn)性變化規(guī)律,不僅能夠反映空間變化趨勢,還可以提高對時間數(shù)據(jù)演變規(guī)律的認知,從而增強對時空數(shù)據(jù)的分析能力[27]。本文選擇2000-2019 年的數(shù)據(jù),研究NPP值與降水、氣溫、濕潤度、光照時長、高程、植被覆蓋度(NDVI)的相互關系,計算回歸系數(shù),同時也能根據(jù)回歸系數(shù)依據(jù)氣候NDVI 等值對未來NPP 值進行預測。時空地理加權回歸模型公式如下:

    3 結果與分析

    3.1 植被凈初級生產力的時空變化分析

    統(tǒng)計2000-2019 年各像元點NPP 均值的空間分布特征(圖2),并計算各年份區(qū)域內NPP 的均值,2000-2019 年區(qū)域NPP 平均值約為208.96 g·(m2·a)-1,NPP 值范圍為0~727.0 g·(m2·a)-1。NPP 較高的區(qū)域主要位于東部和南部,地表資源主要為森林資源和耕地資源,較低的區(qū)域主要位于西北部的荒漠地區(qū)(圖1、圖2)。從時間序列上看,2000-2019 年,內蒙古草原自然資源大區(qū)NPP 均值波動較為明顯,但整體上呈現(xiàn)較快的上升趨勢,平均每年上升約4.06 g·(m2·a)-1,NPP 均 值 在2001 年 最 低,僅 約150.55 g·(m2·a)-1,在2019 年 最 高,NPP 均 值 約 為247.79 g·(m2·a)-1。在2000-2019 年,2001-2002、2007-2008、2011-2012 和2017-2018 年,區(qū)域內NPP 均值急劇上升,在2000-2001、2006-2007 和2014-2015年下降較為明顯,其余年份變化較為平緩(圖3)。

    二是松花江黑龍江洪水量級大,持續(xù)時間長。松花江發(fā)生1998年以來最大流域性較大洪水,其中嫩江上游發(fā)生超50年一遇特大洪水,松花江上游發(fā)生超20年一遇大洪水。黑龍江發(fā)生1984年以來最大流域性大洪水,下游發(fā)生超100年一遇特大洪水。嫩江、松花江干流水位超警戒歷時46天,黑龍江干流水位超警戒歷時58天。

    圖3 2000-2019 年內蒙古草原自然資源大區(qū)植被凈初級生產力均值及變化趨勢圖Figure 3 Mean net primary productivity (NPP) values and their change trend for natural resources in the grassland region of Inner Mongolia from 2000 to 2019

    基于2000-2019 年內蒙古草原自然資源大區(qū)的NPP 柵格數(shù)據(jù),采用Theil-sen 斜率趨勢分析方法和Mann-Kendall 顯 著 性 檢 驗 方 法,取|Z| > 1.645、1.96 和2.576,表示研究序列通過了置信度為90%、95%和99%的顯著性檢驗,分別對應為弱顯著、顯著和極顯著。得到NPP 值變化趨勢以及變化顯著性的空間分布圖(圖4a),大部分區(qū)域NPP 值呈上升趨勢,約占整個區(qū)域的94.01%,僅5.99%的區(qū)域呈現(xiàn)減少趨勢。并且大部分區(qū)域NPP 值增加趨勢表現(xiàn)為極顯著增加,約占總研究區(qū)面積的60.31%,其次為顯著增加,約占研究區(qū)域的19.57%。NPP 呈減少趨勢的區(qū)域零星分布,其中不顯著減少所占的區(qū)域最多,約5.88%,極顯著、弱顯著和顯著減少的區(qū)域僅占約0.15%,零星分布于建設用地相關的區(qū)域,僅有約1.01 × 105hm2。

    圖4 2000-2019 年植被凈初級生產力值變化趨勢和變化顯著性(a)以及可持續(xù)性(b)Figure 4 Net primary productivity change trend and significance (a) and persistence of the change trend (b) from 2000 to 2019

    將NPP 的變化趨勢與Hurst 指數(shù)相疊加得出NPP 值的變化趨勢和持續(xù)形式,分為持續(xù)減少(2000-2019 年呈現(xiàn)減少趨勢,H> 0.5)、潛在增加(2000-2019 年呈減少趨勢,H< 0.5)、無變化趨勢(H= 0.5)、潛在減少(2000-2019 年呈現(xiàn)增加趨勢,H< 0.5)、持續(xù)增加(2000-2019 年呈現(xiàn)增加趨勢,H> 0.5)。由圖4b 可知,在2000-2019 年,盡管內蒙古草原資源大區(qū)94%的區(qū)域呈上升趨勢,但未來潛在減少的區(qū)域最多,約占整個區(qū)域面積的74.15%,分布于該區(qū)域西部,南部的大部分地區(qū),其次為未來持續(xù)增長的區(qū)域,約占整個區(qū)域的19.45%,主要分布于該大區(qū)的東部和東北部的有林地地區(qū)。處于不確定游離狀態(tài)的區(qū)域僅有約5.83%,主要分布在荒漠和濕地地區(qū)。未來NPP 呈現(xiàn)持續(xù)減少的區(qū)域僅有約0.23%,零星分布于該大區(qū)的中部。

    3.2 NPP 值時空變化的驅動因子分析

    通過時空地理加權回歸模型,選取2000、2005、2010、2015、2019 年的降水、年平均氣溫、光照時長、相對濕度等氣候因子以及海拔高度地形因子和NDVI 植被覆蓋度因子。通過Arcgis 漁網(wǎng)工具全局選取了1 056 個采樣點,研究區(qū)域內NPP 值時空變化與各因子之間的關系,然后進行點的插值得出植被凈初級生產力值變化和各影響要素回歸系數(shù)分布(圖5、圖6)。最終通過時空地理加權回歸模型(GTWR)得出的模擬擬合度R2= 0.925,該結果表明回歸精度較好,可以進行分析研究。

    圖6 2019 年植被凈初級生產力值變化、各影響要素回歸系數(shù)分布Figure 6 Changes in NPP of vegetation in 2019, distribution of regression coefficients of various influencing factors

    統(tǒng)計區(qū)域2000-2019 各影響因子得到平均回歸系數(shù),2000-2019 年降水、相對濕度和NDVI 變化和NPP 值相關性較強,主要為正向相關關系,氣溫、光照時長和蒸散發(fā)主要為負向相關關系(表1)。在2000 年,降水對NPP 值變化的回歸系數(shù)最高值為1.113,平均值為0.368,呈負相關的區(qū)域約占5%,在空間分布上,主要分布在西南少部分地區(qū),其他區(qū)域呈正相關,影響程度從西南到東北存在明顯的空間分異性,從西南到東北地區(qū)呈現(xiàn)線性增加的趨勢;2000 年,氣溫對NPP 影響主要為負相關,回歸系數(shù)的平均值為-0.027,最高值為0.803,最低值為-0.631,其中正向相關影響的主要分布于該區(qū)錫林郭勒盟市的西部,占整個研究區(qū)域的20%,其他區(qū)域呈負向相關關系,并且存在明顯的空間異質性,在研究區(qū)域西南影響程度較大,東北部次之,中部最小;光照時長對NPP 值的影響相對較小,在2000 年影響平均回歸系數(shù)為-0.061,最高值為-0.051,最低值為-0.066,幾乎所有區(qū)域呈負相關性影響,影響程度從西南到東北有顯著的線性減小趨勢;濕潤度對NPP 值的影響在2000 年影響平均回歸系數(shù)為0.228,最高值為0.233,最低值為0.221,幾乎所有區(qū)域呈現(xiàn)正向相關,影響程度從西南到東北呈線性增加的趨勢;2000 年蒸散發(fā)對NPP 的影響平均回歸系數(shù)為-0.020,最高值為-0.030,影響程度由西南到東北區(qū)域逐漸遞減;NDVI 值與NPP 值主要是正向相關關系,且回歸系數(shù)較高,其中在西南地區(qū)最高,由西南向東北呈遞減趨勢;2000 年海拔高度與NPP 值回歸系數(shù)的平均值為0.080,最低值為-0.032,最高值為0.177,其中海拔高度與NPP 呈現(xiàn)負相關性的區(qū)域占37.45%,主要是在該大區(qū)的西南部,其他區(qū)域呈現(xiàn)正相關,影響程度從該區(qū)西南到東北有減小的趨勢(圖5)。

    表1 2000-2019 年植被凈初級生產力影響因子回歸系數(shù)Table 1 Regression coefficients for net primary productivity influencing factors, 2000-2019

    圖5 2000 年植被凈初級生產力值變化、各影響要素回歸系數(shù)分布Figure 5 Changes in NPP of vegetation in 2000, distribution of regression coefficients of various influencing factors

    2019 年,降水、溫度、光照時長、濕度、蒸散發(fā)等氣候因子對NPP 值影響的回歸系數(shù)的平均值分別為0.257、-0.145、-0.201、0.081 和-0.001,海拔和NDVI 值的回歸系數(shù)分別為0.078 和0.559。與2000相比,降水、濕度和蒸散的影響程度在減少,光照、溫度和NDVI 值的影響程度有所增加,海拔影響基本不變(表1)。2019 年,降水對NPP 值的影響的回歸系數(shù)最高值為0.715,最低值為-0.158,其中呈負相關的區(qū)域,約占10%,主要在呼倫貝爾市西北和巴彥淖爾市東南,呈正相關的區(qū)域影響程度由西南向東北呈增加趨勢;氣溫對NPP 的影響,最高值為0.901,最低值為-1.109,大部分區(qū)域呈現(xiàn)負相關影響,約占59.56%,呈正相關的區(qū)域主要在該大區(qū)東北部和巴彥淖爾市西部;2019 年,光照時長對NPP 的影響回歸系數(shù)最高值為-0.197 4,最低值為-0.202 4,平均值為-0.201,所有區(qū)域光照時長與NPP 值呈現(xiàn)負相關影響,影響程度從西南到東北呈線性減弱;相對濕度對NPP 值的回歸系數(shù)最高為0.077,最低為0.084,影響程度較小,全部區(qū)域濕度與NPP 值呈正相關影響,與2000 年不同,影響程度由西南到東北逐漸減弱;2019 年地表蒸散對NPP 值影響的回歸系數(shù)最高值為-0.08,最低值為-0.01,影響程度變化不大,由西南到東北遞減;2019 年NDVI值和NPP 值的相關系數(shù)最低值為0.08,平均值為0.559,所有區(qū)域呈現(xiàn)正相關關系,其中在東北部和西南部的部分地區(qū),NDVI 值對NPP 的影響較小,影響程度存在顯著的聚集性;海拔高度對NPP 值的影響與2000 有所不同,最低值為0.048,最高值為0.118,全區(qū)域呈正向關,有顯著的空間異質性,由西南向東北遞增(圖6)。

    4 討論

    本研究通過對2000-2019 年NPP 值的變化及其驅動因子研究得出,降水量變化是內蒙古草原自然資源大區(qū)NPP 值增加的主要氣候因素,NDVI 值與NPP 值密切相關,但NDVI 值越高的地區(qū),NPP受NDVI 變化的影響相對較小。另外內蒙古草原自然資源大區(qū)大部分區(qū)域溫度和光照時長的變化與NPP 值呈負相關關系,同時也有學者研究指出一定程度的氣溫升高可以提高光合效率和資源利用效率從而增加NPP[29-30],可能因為在內蒙古半干旱地區(qū),由于水分較少,溫度升高會引起水分的減少,從而使植被利用水分減少反而降低了光合效率,地表水分是影響NPP 變化的主要原因。越干旱的地區(qū)溫度、光照時長和蒸散發(fā)對NPP 值的負相關性就越強(圖5、圖6)。而NPP 與溫度呈正相關的區(qū)域主要在東北部和錫林郭盟市西部植被覆蓋指數(shù)較高的地方。各類要素對NPP 值的影響從該大區(qū)西南到東北有顯著的線性變化,有明顯的空間異質性,可以為區(qū)域碳儲量和自然資源的分區(qū)研究提供幫助。

    目前本研究還存在許多不足,主要對NPP 值時空變化的趨勢進行分析。以及利用氣候、海拔、植被覆蓋度等因子對NPP 值進行時空地理加權回歸分析,通過回歸系數(shù)研究各因子對NPP 值的影響。主要還是以統(tǒng)計學以及統(tǒng)計學相關的模型進行分析,沒有深入研究內部變化機理,所得到的影響系數(shù)只是統(tǒng)計學上的相關性。由于2000 年和2019 年的氣溫和植被類型存在較大差異,有些區(qū)域的NPP值的影響回歸系數(shù)在2000 年與2019 年相差較大,甚至相關趨勢都不同,這可能是不同氣候對NPP 的影響存在一個閾值,這個閾值需要大量的時空數(shù)據(jù)分析,可能是以后需要研究的方向。另外NPP 的變化可能與人類的生產活動有很大相關性。綜合區(qū)域自然因子和工農業(yè)、畜牧業(yè)、人口密度、GDP 等社會經(jīng)濟因子研究NPP 的變化是今后需要解決的問題。

    5 結論

    在2000-2019 年間內蒙古草原自然資源區(qū)域NPP 均值呈波動上升趨勢,平均每年增加約207.72 g·(m2·a)-1。在2000-2019 年94.01%的區(qū)域NPP 值呈上升趨勢,79.88%的區(qū)域NPP 值呈顯著上升(|Z| >1.96),呈顯著下降的區(qū)域極少,僅有約0.15%;利用Hurst 指數(shù)來判斷每個像元點變化的可持續(xù)性,得出盡管2000-2019 年有94.01%的區(qū)域NPP 值呈現(xiàn)上升趨勢,但在未來會持續(xù)增加的區(qū)域僅占整個研究區(qū)域的17.45%主要分布在研究區(qū)域東部和東北部地區(qū),74.15%的區(qū)域在未來有潛在下降的趨勢。

    采用時空地理加權回歸模型,研究NPP 值變化的影響因子。得出在2000-2019 年,該研究區(qū)域內降水、相對濕度和NDVI 值的影響主要呈正相關,降水和NDVI 值對NPP 的影響較高,在2000 和2019年相關系數(shù)分別為0.368、0.257 和0.402、0.259。溫度、蒸散發(fā)和光照時間主要呈負相關關系。其中溫度和光照時間影響相對較大,2000 年和2019 年回歸系數(shù)平均值分別為-0.027、-0.145 和-0.061、-0.201。蒸散發(fā)影響較小,2000 年和2019 年回歸系數(shù)的平均值分別為-0.02 和-0.001。

    各因子對NPP 的影響存在空間異質性。降水與NPP 值的變化主要呈正相關關系,但在2000 年和2019 年仍然有一部分區(qū)域呈現(xiàn)負相關,分別占5%和10%,主要在該大區(qū)西南和東北的部分區(qū)域。溫度主要是負相關,2000 和2019 呈正相關的區(qū)域主要在錫林郭盟西部;另外光照時長、相對濕度、蒸散量和高程對NPP 值的影響在不同區(qū)域也有所不同,從該區(qū)西南到東北有顯著的線性變化趨勢。

    猜你喜歡
    大區(qū)回歸系數(shù)時空
    涪陵:工業(yè)大區(qū)打開綠色新場景
    當代黨員(2022年9期)2022-05-20 13:35:21
    跨越時空的相遇
    智利第一大區(qū)HNX礦區(qū)IOCG型銅礦床地質特征
    礦產勘查(2020年8期)2020-12-25 02:47:00
    鏡中的時空穿梭
    玩一次時空大“穿越”
    多元線性回歸的估值漂移及其判定方法
    電導法協(xié)同Logistic方程進行6種蘋果砧木抗寒性的比較
    多元線性模型中回歸系數(shù)矩陣的可估函數(shù)和協(xié)方差陣的同時Bayes估計及優(yōu)良性
    時空之門
    On Commuting Graph of Group Ring ZnS3?
    国产女主播在线喷水免费视频网站| 日本色播在线视频| 婷婷丁香在线五月| 男人爽女人下面视频在线观看| 欧美 日韩 精品 国产| 18禁国产床啪视频网站| 国产精品 国内视频| 国产成人av教育| 久久99热这里只频精品6学生| 亚洲国产精品成人久久小说| 精品亚洲成a人片在线观看| 亚洲中文av在线| 国产一区二区三区av在线| 色综合欧美亚洲国产小说| 少妇人妻 视频| 久久热在线av| 久久精品国产综合久久久| 99热全是精品| 亚洲,欧美,日韩| 久久久久久久久久久久大奶| 欧美中文综合在线视频| 99国产精品一区二区三区| 亚洲伊人色综图| 女人久久www免费人成看片| 久久99精品国语久久久| 黄色怎么调成土黄色| 午夜两性在线视频| 国产成人一区二区在线| 老鸭窝网址在线观看| 国产一区亚洲一区在线观看| www日本在线高清视频| 亚洲成人免费av在线播放| 90打野战视频偷拍视频| 久久久久久亚洲精品国产蜜桃av| 无限看片的www在线观看| 97在线人人人人妻| 成年美女黄网站色视频大全免费| 777米奇影视久久| 夜夜骑夜夜射夜夜干| 精品亚洲乱码少妇综合久久| 欧美激情高清一区二区三区| av网站在线播放免费| 国产在视频线精品| 欧美97在线视频| 日本a在线网址| 久久久久久亚洲精品国产蜜桃av| 中文乱码字字幕精品一区二区三区| 亚洲国产毛片av蜜桃av| 亚洲欧美精品自产自拍| svipshipincom国产片| 下体分泌物呈黄色| 欧美变态另类bdsm刘玥| 国产一区二区三区av在线| 精品一区二区三区四区五区乱码 | 一个人免费看片子| 精品国产乱码久久久久久男人| 99热网站在线观看| 咕卡用的链子| 欧美av亚洲av综合av国产av| 亚洲黑人精品在线| av一本久久久久| 亚洲专区国产一区二区| 一二三四社区在线视频社区8| 中文欧美无线码| 黄色a级毛片大全视频| 久久人人爽人人片av| 国产精品国产av在线观看| 大型av网站在线播放| 国产一区亚洲一区在线观看| 秋霞在线观看毛片| 久久久久精品人妻al黑| 国产成人91sexporn| 99九九在线精品视频| av在线app专区| 国产有黄有色有爽视频| 免费av中文字幕在线| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区| 亚洲国产av影院在线观看| 国产免费一区二区三区四区乱码| 亚洲少妇的诱惑av| 亚洲黑人精品在线| 精品第一国产精品| 亚洲五月色婷婷综合| 国产精品久久久人人做人人爽| 久久久久久久久久久久大奶| 亚洲欧美日韩高清在线视频 | 曰老女人黄片| 尾随美女入室| 成年人黄色毛片网站| 日本av手机在线免费观看| 老司机靠b影院| 国产一区二区激情短视频 | 中文字幕高清在线视频| 搡老岳熟女国产| 高清黄色对白视频在线免费看| 熟女少妇亚洲综合色aaa.| 好男人视频免费观看在线| 极品少妇高潮喷水抽搐| 男女下面插进去视频免费观看| 免费久久久久久久精品成人欧美视频| 妹子高潮喷水视频| 国产精品九九99| 日本wwww免费看| 色播在线永久视频| tube8黄色片| 亚洲七黄色美女视频| 高清av免费在线| 亚洲欧美一区二区三区国产| 亚洲久久久国产精品| 婷婷色综合大香蕉| 色婷婷久久久亚洲欧美| 午夜免费男女啪啪视频观看| 黄色片一级片一级黄色片| 99九九在线精品视频| 国产在线观看jvid| 亚洲国产欧美一区二区综合| 国产日韩一区二区三区精品不卡| 久久精品国产a三级三级三级| 国产有黄有色有爽视频| 青春草亚洲视频在线观看| 老汉色∧v一级毛片| av有码第一页| 欧美中文综合在线视频| 亚洲精品中文字幕在线视频| 精品国产乱码久久久久久男人| 国产亚洲欧美精品永久| 九色亚洲精品在线播放| 美女高潮到喷水免费观看| 操美女的视频在线观看| 欧美日韩精品网址| 在线精品无人区一区二区三| 亚洲免费av在线视频| 日本熟妇午夜| 91麻豆精品激情在线观看国产| 好男人在线观看高清免费视频 | 欧美在线一区亚洲| 亚洲天堂国产精品一区在线| 人人妻人人澡欧美一区二区| 久久久久九九精品影院| 好看av亚洲va欧美ⅴa在| 美女午夜性视频免费| 久久精品成人免费网站| 精品高清国产在线一区| 色综合亚洲欧美另类图片| 亚洲成av人片免费观看| 免费在线观看视频国产中文字幕亚洲| 国产精品久久久av美女十八| 午夜精品在线福利| 99久久99久久久精品蜜桃| 国产真人三级小视频在线观看| 最近最新免费中文字幕在线| 国产激情偷乱视频一区二区| 欧美在线黄色| 给我免费播放毛片高清在线观看| 日本熟妇午夜| 国产黄a三级三级三级人| 十分钟在线观看高清视频www| 欧美不卡视频在线免费观看 | 嫁个100分男人电影在线观看| 淫妇啪啪啪对白视频| 波多野结衣高清无吗| 满18在线观看网站| 国产激情偷乱视频一区二区| 一区福利在线观看| 久9热在线精品视频| 成人一区二区视频在线观看| 日韩成人在线观看一区二区三区| 1024手机看黄色片| 老熟妇乱子伦视频在线观看| 欧美性猛交黑人性爽| 成年女人毛片免费观看观看9| 又紧又爽又黄一区二区| 久久久精品欧美日韩精品| 熟女少妇亚洲综合色aaa.| 国产99白浆流出| 国产精品九九99| 久久久久国产一级毛片高清牌| 亚洲人成网站在线播放欧美日韩| 视频区欧美日本亚洲| 精品日产1卡2卡| 国产亚洲av高清不卡| 女性生殖器流出的白浆| 精品日产1卡2卡| 亚洲精品中文字幕在线视频| 日韩欧美免费精品| 久久久久久久久免费视频了| 99re在线观看精品视频| 99久久国产精品久久久| 白带黄色成豆腐渣| 久久久久久国产a免费观看| 成人三级做爰电影| 亚洲一区二区三区色噜噜| 听说在线观看完整版免费高清| 亚洲精品国产一区二区精华液| 亚洲中文字幕一区二区三区有码在线看 | 中文字幕av电影在线播放| 亚洲国产精品sss在线观看| 18美女黄网站色大片免费观看| 成人一区二区视频在线观看| 久久精品夜夜夜夜夜久久蜜豆 | 成人国语在线视频| 午夜视频精品福利| 国产精品久久视频播放| 又黄又粗又硬又大视频| 亚洲精品国产精品久久久不卡| 人妻久久中文字幕网| 欧美亚洲日本最大视频资源| 日本 欧美在线| 国产激情偷乱视频一区二区| 国产精品,欧美在线| 精品一区二区三区四区五区乱码| 精品国产亚洲在线| 99国产精品一区二区蜜桃av| 可以在线观看毛片的网站| 美女午夜性视频免费| 麻豆成人午夜福利视频| 日韩精品中文字幕看吧| 亚洲狠狠婷婷综合久久图片| 特大巨黑吊av在线直播 | 99久久精品国产亚洲精品| 一级a爱视频在线免费观看| videosex国产| 91老司机精品| 国产精品久久久人人做人人爽| 精品一区二区三区四区五区乱码| 日本成人三级电影网站| 亚洲一码二码三码区别大吗| 9191精品国产免费久久| 中文字幕人妻丝袜一区二区| 亚洲精品国产一区二区精华液| 999久久久国产精品视频| 亚洲国产欧美一区二区综合| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情极品国产一区二区三区| 最近最新中文字幕大全电影3 | 久久久国产成人免费| 女性被躁到高潮视频| 女生性感内裤真人,穿戴方法视频| 亚洲avbb在线观看| 精品人妻1区二区| 桃红色精品国产亚洲av| 超碰成人久久| 极品教师在线免费播放| 日日爽夜夜爽网站| 深夜精品福利| 久久久久精品国产欧美久久久| 久久久久亚洲av毛片大全| 丝袜人妻中文字幕| 欧美不卡视频在线免费观看 | 国产日本99.免费观看| 日本黄色视频三级网站网址| 亚洲av熟女| 精品高清国产在线一区| АⅤ资源中文在线天堂| 国产精品永久免费网站| 精品国产一区二区三区四区第35| 热99re8久久精品国产| 国产视频一区二区在线看| 日韩欧美免费精品| 一级毛片精品| 国产欧美日韩精品亚洲av| 法律面前人人平等表现在哪些方面| 日韩精品青青久久久久久| 久久久久国内视频| 国产精品亚洲美女久久久| 一二三四社区在线视频社区8| 国产精品久久电影中文字幕| 亚洲精品一区av在线观看| 国产av又大| 手机成人av网站| a级毛片a级免费在线| 黄片播放在线免费| 久久精品国产清高在天天线| 午夜久久久在线观看| 99国产极品粉嫩在线观看| 午夜福利视频1000在线观看| 少妇 在线观看| 中文字幕另类日韩欧美亚洲嫩草| 色播在线永久视频| 琪琪午夜伦伦电影理论片6080| 一本精品99久久精品77| 巨乳人妻的诱惑在线观看| 女警被强在线播放| 一区二区三区高清视频在线| 91国产中文字幕| 国产精品国产高清国产av| a级毛片a级免费在线| 人妻丰满熟妇av一区二区三区| 少妇被粗大的猛进出69影院| 校园春色视频在线观看| 成人三级黄色视频| e午夜精品久久久久久久| 少妇粗大呻吟视频| 久久久国产成人免费| 一边摸一边做爽爽视频免费| 日本 欧美在线| 欧美性长视频在线观看| 嫩草影院精品99| 国产99久久九九免费精品| 久久天堂一区二区三区四区| 成人18禁在线播放| 一区二区日韩欧美中文字幕| 99热这里只有精品一区 | 亚洲自偷自拍图片 自拍| 久久久久久久午夜电影| 国产1区2区3区精品| 亚洲激情在线av| 欧美中文日本在线观看视频| 亚洲熟女毛片儿| 免费高清在线观看日韩| 中文字幕高清在线视频| 国产av一区二区精品久久| 麻豆成人午夜福利视频| 悠悠久久av| 国产精华一区二区三区| 精品国产超薄肉色丝袜足j| 麻豆av在线久日| 岛国在线观看网站| 国产成人精品久久二区二区91| 在线观看舔阴道视频| 国产一区二区在线av高清观看| 久久热在线av| 欧美黑人欧美精品刺激| 人人澡人人妻人| 91麻豆精品激情在线观看国产| 亚洲精品国产一区二区精华液| 亚洲国产毛片av蜜桃av| 精品欧美一区二区三区在线| 精华霜和精华液先用哪个| 两人在一起打扑克的视频| 天天躁夜夜躁狠狠躁躁| 久99久视频精品免费| 不卡一级毛片| 欧美色视频一区免费| 亚洲电影在线观看av| 最近最新免费中文字幕在线| 一级黄色大片毛片| 久久精品91无色码中文字幕| 黑人巨大精品欧美一区二区mp4| 久久亚洲真实| 国产精品98久久久久久宅男小说| av视频在线观看入口| 熟妇人妻久久中文字幕3abv| a级毛片a级免费在线| 窝窝影院91人妻| 亚洲精品美女久久av网站| 亚洲avbb在线观看| 国产免费av片在线观看野外av| 色av中文字幕| 久久狼人影院| 日韩欧美在线二视频| www日本黄色视频网| 久久精品国产99精品国产亚洲性色| 女人高潮潮喷娇喘18禁视频| 国产97色在线日韩免费| 美女大奶头视频| 少妇的丰满在线观看| 国产精品影院久久| 在线av久久热| 久久国产精品男人的天堂亚洲| 精品一区二区三区四区五区乱码| 亚洲av成人不卡在线观看播放网| 99热6这里只有精品| 久久久久久人人人人人| 婷婷六月久久综合丁香| 亚洲精品国产一区二区精华液| 午夜福利欧美成人| 国产精品 欧美亚洲| 草草在线视频免费看| 亚洲avbb在线观看| 国产欧美日韩一区二区三| 正在播放国产对白刺激| 日本精品一区二区三区蜜桃| 丁香六月欧美| 免费高清在线观看日韩| 国产伦一二天堂av在线观看| av视频在线观看入口| 欧美激情久久久久久爽电影| 久久天躁狠狠躁夜夜2o2o| 久热爱精品视频在线9| 51午夜福利影视在线观看| 国产免费男女视频| 久久久久久久精品吃奶| 变态另类成人亚洲欧美熟女| 国产色视频综合| 欧美成人午夜精品| 国产色视频综合| 不卡一级毛片| 三级毛片av免费| av天堂在线播放| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 精品午夜福利视频在线观看一区| 国产精品香港三级国产av潘金莲| 高清在线国产一区| 天堂动漫精品| 九色国产91popny在线| 99热只有精品国产| 中文字幕av电影在线播放| videosex国产| 又黄又粗又硬又大视频| 欧美在线一区亚洲| 中文亚洲av片在线观看爽| 成人手机av| 亚洲中文字幕一区二区三区有码在线看 | 国产精品久久久久久精品电影 | 日韩精品免费视频一区二区三区| 在线看三级毛片| 自线自在国产av| 给我免费播放毛片高清在线观看| 亚洲av五月六月丁香网| 母亲3免费完整高清在线观看| 手机成人av网站| 18美女黄网站色大片免费观看| 亚洲三区欧美一区| 免费在线观看黄色视频的| 91字幕亚洲| 老司机午夜十八禁免费视频| 久久国产精品男人的天堂亚洲| 亚洲国产中文字幕在线视频| 18禁国产床啪视频网站| 欧美一级a爱片免费观看看 | 亚洲一卡2卡3卡4卡5卡精品中文| 人人妻,人人澡人人爽秒播| 亚洲一区高清亚洲精品| 亚洲片人在线观看| 两个人看的免费小视频| 亚洲av成人av| 精品久久久久久久人妻蜜臀av| 嫩草影院精品99| 不卡av一区二区三区| 国产高清有码在线观看视频 | 动漫黄色视频在线观看| 一级作爱视频免费观看| 自线自在国产av| e午夜精品久久久久久久| 国产一区在线观看成人免费| 成人18禁在线播放| 国产av不卡久久| 日日摸夜夜添夜夜添小说| 国产不卡一卡二| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 淫妇啪啪啪对白视频| 琪琪午夜伦伦电影理论片6080| 国产亚洲精品一区二区www| 狠狠狠狠99中文字幕| 午夜激情av网站| 久久久久九九精品影院| 999久久久国产精品视频| 精品无人区乱码1区二区| 亚洲欧美日韩无卡精品| 级片在线观看| 日韩三级视频一区二区三区| 亚洲av片天天在线观看| 男女之事视频高清在线观看| 久久国产精品影院| 国产真实乱freesex| 男女午夜视频在线观看| 午夜成年电影在线免费观看| 99精品欧美一区二区三区四区| 中文字幕高清在线视频| 久久久久久久精品吃奶| 亚洲一区二区三区色噜噜| x7x7x7水蜜桃| 日本 av在线| 一本精品99久久精品77| 久久久久久国产a免费观看| 一本大道久久a久久精品| 日韩av在线大香蕉| 不卡一级毛片| 亚洲一区高清亚洲精品| 亚洲激情在线av| 免费一级毛片在线播放高清视频| 黄频高清免费视频| 日本成人三级电影网站| 亚洲人成电影免费在线| 亚洲国产精品999在线| 成人三级做爰电影| 天堂动漫精品| 亚洲专区字幕在线| 一边摸一边抽搐一进一小说| 国产成人欧美| 免费看日本二区| 久久天堂一区二区三区四区| 一区二区三区国产精品乱码| 日韩av在线大香蕉| 成人av一区二区三区在线看| 老司机午夜福利在线观看视频| 久久天堂一区二区三区四区| 久久久久久久久中文| 99国产精品一区二区蜜桃av| 黄色 视频免费看| 久久久久久亚洲精品国产蜜桃av| 香蕉av资源在线| 一级毛片精品| 久久久久久久精品吃奶| 亚洲成人久久性| 亚洲一区高清亚洲精品| 久久中文看片网| 国产成人精品久久二区二区免费| 一级片免费观看大全| 身体一侧抽搐| 久久草成人影院| 国产精品,欧美在线| 午夜日韩欧美国产| 色精品久久人妻99蜜桃| 亚洲专区中文字幕在线| 午夜久久久在线观看| 久久久水蜜桃国产精品网| 超碰成人久久| 欧美成狂野欧美在线观看| 免费看a级黄色片| 国产高清videossex| 视频区欧美日本亚洲| 老熟妇乱子伦视频在线观看| 中文字幕人妻熟女乱码| 男人舔女人下体高潮全视频| 脱女人内裤的视频| 国产成人精品久久二区二区91| 999精品在线视频| 在线av久久热| 99精品欧美一区二区三区四区| 好看av亚洲va欧美ⅴa在| 国产成人影院久久av| 日本五十路高清| 国产av一区在线观看免费| 欧美性猛交黑人性爽| 国产黄片美女视频| 无人区码免费观看不卡| 久久精品国产亚洲av香蕉五月| 亚洲av中文字字幕乱码综合 | 亚洲中文字幕一区二区三区有码在线看 | 午夜福利免费观看在线| 老司机在亚洲福利影院| 麻豆成人午夜福利视频| 男人舔女人下体高潮全视频| 免费一级毛片在线播放高清视频| 夜夜爽天天搞| 好男人电影高清在线观看| 国产激情久久老熟女| 亚洲中文av在线| 国产黄片美女视频| 久久久精品国产亚洲av高清涩受| 丰满的人妻完整版| 色综合婷婷激情| 午夜福利在线观看吧| 欧美又色又爽又黄视频| 级片在线观看| 精品国产乱码久久久久久男人| 精品国产亚洲在线| 精品欧美一区二区三区在线| 变态另类丝袜制服| 国语自产精品视频在线第100页| 国产麻豆成人av免费视频| 国产蜜桃级精品一区二区三区| 国产亚洲精品综合一区在线观看 | 成人18禁高潮啪啪吃奶动态图| 俄罗斯特黄特色一大片| 妹子高潮喷水视频| 老司机福利观看| av有码第一页| 一级黄色大片毛片| 亚洲一区高清亚洲精品| 一二三四社区在线视频社区8| 岛国在线观看网站| 黄色片一级片一级黄色片| 色综合亚洲欧美另类图片| 91九色精品人成在线观看| 黄片大片在线免费观看| 97人妻精品一区二区三区麻豆 | 香蕉av资源在线| 美国免费a级毛片| 一夜夜www| 悠悠久久av| 国产久久久一区二区三区| 嫩草影院精品99| 亚洲精品在线美女| 婷婷丁香在线五月| 久久久久久久午夜电影| 欧美成人午夜精品| 精品欧美一区二区三区在线| 国产三级黄色录像| 国产亚洲av高清不卡| 国产单亲对白刺激| 国产精品九九99| 精品不卡国产一区二区三区| 欧美中文日本在线观看视频| 国产黄a三级三级三级人| 一二三四在线观看免费中文在| 国内精品久久久久久久电影| 美女高潮喷水抽搐中文字幕| 色老头精品视频在线观看| 美女 人体艺术 gogo| 香蕉av资源在线| 欧美性长视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲成人国产一区在线观看| 亚洲午夜理论影院| 成人手机av| 亚洲片人在线观看| 国产熟女xx| 一级作爱视频免费观看| 他把我摸到了高潮在线观看| 老司机深夜福利视频在线观看| 最近在线观看免费完整版| 午夜福利欧美成人| 91在线观看av| 99在线视频只有这里精品首页| 国产av一区二区精品久久| 一本综合久久免费| 欧美日韩福利视频一区二区| 1024香蕉在线观看| 久久久久国产精品人妻aⅴ院|