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

    三江平原典型區(qū)水田時(shí)空變化及驅(qū)動因素分析

    2019-05-11 06:18:06張文琦
    關(guān)鍵詞:富錦市墾區(qū)水田

    張文琦,宋 戈

    ?

    三江平原典型區(qū)水田時(shí)空變化及驅(qū)動因素分析

    張文琦,宋 戈※

    (東北大學(xué)土地管理研究所,沈陽 110169)

    為闡明區(qū)域水田時(shí)空變化特征及其驅(qū)動因素,該文以三江平原典型地區(qū)富錦市為例,運(yùn)用網(wǎng)格單元法和地理探測器模型分析不同時(shí)段(1989—2002年、2002—2015年)水田時(shí)空變化特征及驅(qū)動因素,探討各驅(qū)動因子在水田時(shí)空變化中的作用及其相互關(guān)系。結(jié)果表明:1)1989—2015年富錦市水田擴(kuò)張明顯,新增水田34.99萬hm2,墾區(qū)水田化進(jìn)程早于農(nóng)區(qū)。1989—2002年,東部墾區(qū)水田擴(kuò)張度較高,南部次之;2002—2015年,水田擴(kuò)張度較高的地區(qū)主要分布在西南部農(nóng)區(qū)與北部墾區(qū)。2)1989—2002年,水田擴(kuò)張的主要驅(qū)動因素為政策因素和土壤類型,各因子間的交互作用大部分是雙因子增強(qiáng);2002—2015年,自然因素對水田變化影響減弱,水田擴(kuò)張的主要驅(qū)動因素是到河流距離和到鐵路距離,各因子間的交互作用大部分是非線性增強(qiáng),水田化發(fā)展到成熟階段。富錦市水田時(shí)空變化差異是由多種驅(qū)動因子共同作用的結(jié)果,研究結(jié)果可為三江平原耕地資源合理利用及水田規(guī)??茖W(xué)管控提供借鑒和參考。

    土地利用;遙感;時(shí)空變化;耕地;水田;驅(qū)動因素;地理探測器;三江平原

    0 引 言

    水田變化是人類活動影響耕地利用結(jié)構(gòu)的結(jié)果,空間表現(xiàn)為區(qū)域內(nèi)水田種植面積擴(kuò)張或縮小,水田變化不可避免地對區(qū)域農(nóng)業(yè)生產(chǎn)功能和農(nóng)業(yè)生態(tài)系統(tǒng)服務(wù)產(chǎn)生重要影響[1-2]。全球氣候變化與人類活動的雙重脅迫,直接影響地球表層土地利用結(jié)構(gòu)的穩(wěn)定性,使全球和區(qū)域尺度的耕地及耕地內(nèi)部利用結(jié)構(gòu)發(fā)生劇烈變化。三江平原地處中高緯度和歐亞大陸東端,對全球氣候變化十分敏感,地表升溫幅度和降水時(shí)空異質(zhì)性高于全球同期水平[3-5],20世紀(jì)90年代以來為應(yīng)對全球氣候變化,在高強(qiáng)度人類活動干擾下大面積種植水稻,濕地基質(zhì)逐漸變?yōu)檗r(nóng)田基質(zhì),耕地利用結(jié)構(gòu)發(fā)生明顯變化,水田擴(kuò)張?zhí)卣黠@著[6]。水田極端變化無疑會影響區(qū)域內(nèi)糧食供給和資源環(huán)境的保護(hù)[7],因此,認(rèn)識水田變化過程并揭示水田變化與驅(qū)動因子之間的相關(guān)性,對于促進(jìn)耕地利用結(jié)構(gòu)的調(diào)整優(yōu)化、區(qū)域生態(tài)環(huán)境的改善和有針對性的進(jìn)行耕地保護(hù)等具有重要的學(xué)術(shù)意義和實(shí)踐價(jià)值。

    國內(nèi)外基于耕地變化的驅(qū)動力研究已由最初的定性表達(dá),發(fā)展到如今的驅(qū)動因子定量分析。學(xué)者們通常將耕地變化的影響因素分為自然因素和社會經(jīng)濟(jì)因素兩大類[8],自然因素中的高程、坡度等被證明是影響耕地變化的主要因素,社會經(jīng)濟(jì)因素中人口與經(jīng)濟(jì)發(fā)展、交通可達(dá)性等因素也被證實(shí)是耕地變化的主要影響因素[9]。研究方法通常為自相關(guān)分析、面板數(shù)據(jù)模型、空間計(jì)量回歸模型、空間成本加權(quán)模型、地理加權(quán)回歸、空間疊加分析等[10-12]。但以往研究中多因子交互識別是一個(gè)理論難題,目前這一類研究相對較少。由王勁峰[13-14]等提出的地理探測器模型,能有效識別因子間的相互關(guān)系,探測各因子對模型的解釋力,能從龐大的空間數(shù)據(jù)中提取有用的空間關(guān)聯(lián)規(guī)則。然而,目前學(xué)者們針對該地區(qū)水田變化的研究較為少見,僅有部分研究關(guān)注水田數(shù)量變化特征和空間演化規(guī)律[15],尤其是針對水田時(shí)空變化的驅(qū)動力研究更為匱乏。

    富錦市地處三江平原腹地,普通農(nóng)區(qū)與墾區(qū)國營農(nóng)場并存,其耕地變化是整個(gè)三江平原的縮影。本文以富錦市為研究區(qū),為客觀探究水田的時(shí)空變化特征及其驅(qū)動因子,在分析富錦市水田時(shí)空變化差異的基礎(chǔ)上,引入地理探測器模型定量分析水田變化的驅(qū)動因素,揭示其驅(qū)動因子在水田擴(kuò)張過程中的作用及相互關(guān)系,確定不同時(shí)期影響水田擴(kuò)張的主導(dǎo)因子,以期為三江平原耕地保護(hù)和管理提供科學(xué)依據(jù)。

    1 研究區(qū)概況及數(shù)據(jù)來源

    1.1 研究區(qū)概況

    富錦市位于黑龍江省東北部、松花江下游南岸三江平原腹地,地理位置為46°45′35″~47°37′15″N,131°25′12″~133°26′43″E,土地總面積822 400 hm2。全境地勢地平,平均海拔60 m,地貌從西北向東南緩慢傾斜,形成西北略高,中部低平,東南稍低的沖積平原。富錦市屬中溫帶大陸性季風(fēng)氣候區(qū),四季分明,年平均氣溫2.9℃,年降水量573.9 mm,10℃以上活動積溫2 300~2 500℃,熱量狀況較好,適合農(nóng)林牧漁業(yè)生產(chǎn)。富錦市土壤主要有白漿土、草甸土、沼澤土、暗棕壤和黑土,土壤肥沃,可耕性良好。富錦市轄10個(gè)鄉(xiāng)鎮(zhèn),267個(gè)行政村,農(nóng)墾建三江分局及七星、大興、創(chuàng)業(yè)3個(gè)國營農(nóng)場在其境內(nèi),2016年末總?cè)丝?6萬人。富錦市農(nóng)區(qū)與墾區(qū)兩種體制并存,農(nóng)區(qū)耕地屬農(nóng)村集體所有,墾區(qū)耕地屬國家所有,其中富錦市所轄七星農(nóng)場、大興農(nóng)場和創(chuàng)業(yè)農(nóng)場屬于墾區(qū),其他地區(qū)屬于農(nóng)區(qū)。三江平原是全球范圍短時(shí)限內(nèi)經(jīng)歷高強(qiáng)度人類活動的區(qū)域,建國以來經(jīng)歷4次大規(guī)模土地開發(fā)高潮,農(nóng)業(yè)開發(fā)活動非常活躍,20世紀(jì)80年代末至21世紀(jì)初期是第四次墾荒高潮,當(dāng)?shù)卣扇≠Y金補(bǔ)貼等政策推行“旱改水”種植模式,大量低洼旱地改為水田,耕地利用格局變化劇烈,至2002年,富錦市農(nóng)業(yè)種植結(jié)構(gòu)調(diào)整政策基本推行完畢[16]。基于此,本研究時(shí)段確定為1989—2002年和2002—2015年,以探尋該地區(qū)水田變化的空間特征及其驅(qū)動因素。

    圖1 研究區(qū)示意圖

    1.2 數(shù)據(jù)來源與處理

    1.2.1 遙感影像數(shù)據(jù)

    本研究采用6幅Landsat遙感影像:1)獲取1989年研究區(qū)云量符合標(biāo)準(zhǔn)的耕作期Landsat5 TM影像,成像時(shí)間為5月17日和6月12日。2)獲取2002年研究區(qū)云量符合標(biāo)準(zhǔn)的Landsat7 ETM+影像,成像時(shí)間為5月17日和6月24日。3)獲取2015年研究區(qū)云量符合標(biāo)準(zhǔn)的Landsat8 OLI影像,成像時(shí)間5月21日和6月16日。依據(jù)三江平原一年一季農(nóng)作物發(fā)育時(shí)期,水稻一般4月育苗,5月和6月處于插秧和分蘗期,此時(shí)水田在遙感影像中光譜差異明顯[17],故本研究遙感影像選擇耕作初期(5月和6月)為主要數(shù)據(jù)源,分辨率為30 m×30 m,軌道號選擇114/27和115/27。

    根據(jù)2017年11月1日頒布的《土地利用現(xiàn)狀分類》(GB/T 21010—2017),結(jié)合富錦市土地利用特點(diǎn),確定土地利用分類系統(tǒng),包括6個(gè)大類:耕地、林地、草地、水域、建設(shè)用地和濕地,其中耕地分為水田和旱地,本文研究主體水田依據(jù)分類標(biāo)準(zhǔn)可界定為用于種植水稻等水生農(nóng)作物的耕地。基于ENVI 5.1軟件平臺,經(jīng)幾何校正、圖像增強(qiáng)等預(yù)處理,根據(jù)不同地類的色調(diào)、紋理等特征建立解譯標(biāo)志,采用面向?qū)ο笈c人工目視解譯相結(jié)合的方法進(jìn)行遙感影像分類,得到富錦市1989年、2002年、2015年3期土地利用現(xiàn)狀數(shù)據(jù)。采用同期土地調(diào)查資料和典型區(qū)域野外抽樣調(diào)查進(jìn)行精度驗(yàn)證,3期影像解譯總體精度分別為88.72%、91.52%和93.18%,Kappa系數(shù)分別為0.84、0.83和0.87,分類結(jié)果滿足精度要求,因此,通過該方法獲取的富錦市土地利用分類結(jié)果滿足研究所需精度。

    1.2.2 驅(qū)動因子選擇與獲取

    水田變化的空間選擇不但考慮自然因素對于水稻種植環(huán)境的適宜性和適應(yīng)性要求,而且考慮農(nóng)戶生產(chǎn)生活的經(jīng)濟(jì)性和便利性。在此基礎(chǔ)上,本研究結(jié)合已有研究成果[18]與富錦市實(shí)際情況,綜合確定高程、坡度、土壤類型、到河流距離、到鐵路距離和到城鎮(zhèn)距離和政策因素7個(gè)因素作為驅(qū)動因子。

    1)高程。高程是影響水田分布的重要指標(biāo)。DEM數(shù)字高程數(shù)據(jù)下載自地理空間數(shù)據(jù)云(http://www.gscloud.cn/),分辨率為30 m×30 m。

    2)坡度。坡度通過影響水田建設(shè)的成本和便利度來影響水田空間分布。坡度數(shù)據(jù)利用ArcGIS 10.2軟件空間分析功能坡度命令由DEM數(shù)據(jù)提取而來。

    3)土壤類型。不同土壤類型的理化性質(zhì)差異顯著,直接影響水田建設(shè)的成本和糧食產(chǎn)量。富錦市有白漿土、黑土、草甸土等7種土壤類型(不包括河流湖泊),不同土壤的理化性質(zhì)差異明顯,導(dǎo)致水田開發(fā)成本、糧食產(chǎn)量差異較大,所以土壤因子對水田變化具有一定影響。土壤數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn),通過ArcGIS 10.2以富錦市行政區(qū)劃裁剪后進(jìn)行空間柵格化處理。

    4)到河流距離。人類墾殖過程中體現(xiàn)出較強(qiáng)的親水性,河流分布是影響水田變化的重要因素。河流矢量數(shù)據(jù)來源于國家基礎(chǔ)地理信息中心(http://ngcc.sbsm.gov.cn/),距離因子利用ArcGIS 10.2軟件中Euclidean Distance工具獲取。

    5)到鐵路距離。距離因素直接影響著人類對于土地格局的干擾程度,很多研究已經(jīng)證明交通線的分布會對區(qū)域土地利用變化產(chǎn)生影響[19-20]。本研究主要考慮距離因素對土地利用變化(水田變化)的影響,自1974年三江平原福前鐵路開通以來,鐵路一直作為富錦市的主要交通線,因此到鐵路距離對當(dāng)?shù)厮镒兓哂袑?shí)際意義。數(shù)據(jù)來源與處理同上。

    6)到城鎮(zhèn)距離。綜合考慮耕作適宜性、建設(shè)經(jīng)濟(jì)性和生活便利性,到城鎮(zhèn)距離一定程度上影響研究區(qū)水田空間分布特征。數(shù)據(jù)來源與處理同上。

    7)政策因素。墾區(qū)與農(nóng)區(qū)的行政歸屬不同,有關(guān)土地管理的行政命令不同,本研究根據(jù)富錦市行政體制差異將研究區(qū)政策因素劃分為農(nóng)區(qū)和墾區(qū)2部分。在驅(qū)動力分析中政策因素按農(nóng)區(qū)與墾區(qū)的范圍分成2類,分別賦值為0與1,以區(qū)分其類別。

    對土地利用數(shù)據(jù)以及DEM、土壤、鐵路、水系和居民點(diǎn)等數(shù)據(jù)進(jìn)行地圖投影和坐標(biāo)校正,本研究在影像的空間幾何配準(zhǔn)上統(tǒng)一采用CGCS2000投影坐標(biāo)系。

    2 研究方法

    2.1 網(wǎng)格單元法

    網(wǎng)格單元(grid cell,GC)方法可以定量分析水田在時(shí)間和空間上的精細(xì)變化,并且能夠通過將相關(guān)的自然及人文因子轉(zhuǎn)入網(wǎng)格單元內(nèi),定量分析研究區(qū)域內(nèi)水田變化的驅(qū)動力[21-22]。本研究基于ArcGIS平臺,在研究區(qū)范圍內(nèi)建立矢量網(wǎng)格,綜合考慮研究區(qū)面積、研究目的以及計(jì)算機(jī)運(yùn)行效率等因素,通過反復(fù)試驗(yàn),確定網(wǎng)格單元大小為1 000 m×1 000 m,總計(jì)8 611個(gè)網(wǎng)格單元。將基于遙感數(shù)據(jù)獲得的專題圖通過ArcGIS軟件轉(zhuǎn)入對應(yīng)位置的網(wǎng)格單元內(nèi),對富錦市1989—2015年水田擴(kuò)張面積的分散程度進(jìn)行統(tǒng)計(jì)。

    在本研究中,水田擴(kuò)張度定義為每個(gè)統(tǒng)計(jì)網(wǎng)格單元內(nèi)的水田圖斑占網(wǎng)格單元總面積的比例,其動態(tài)變化量能刻畫出水田分布擴(kuò)張或消亡的時(shí)空過程。計(jì)算每個(gè)網(wǎng)格內(nèi)水田擴(kuò)張圖斑的面積,然后計(jì)算在相應(yīng)網(wǎng)格內(nèi)所占的比例,公式為

    式中為水田變化圖斑在每個(gè)網(wǎng)格內(nèi)所占的面積比例,0≤≤1;A為每個(gè)網(wǎng)格內(nèi)水田變化面積,m2;A是每個(gè)網(wǎng)格的面積,m2。

    2.2 地理探測器

    地理探測器模型是探測空間分異性,以及揭示其驅(qū)動因子的一組新興統(tǒng)計(jì)學(xué)方法,此方法核心思想是如果某個(gè)自變量對某個(gè)因變量有重要影響,那么自變量和因變量的空間分布也應(yīng)該具有相似性[23-25]。地理探測器包括4個(gè)探測器:因子探測、交互作用探測、風(fēng)險(xiǎn)區(qū)探測和生態(tài)探測,其獨(dú)特優(yōu)勢是探測兩因子交互作用于因變量,通過分別計(jì)算和比較各單因子值及兩因子疊加后值,可以判斷兩因子是否存在交互作用,以及交互作用的強(qiáng)弱、方向、線性還是非線性。本研究運(yùn)用地理探測器的因子探測器和交互探測器模型,定量分析富錦市水田變化空間差異的驅(qū)動因子。

    (1)因子探測器:探測水田變化的空間分異性,以及探測驅(qū)動因子多大程度上解釋了水田變化屬性的空間分異。用值度量,表達(dá)式為

    (2)交互探測器:定量表征2個(gè)不同驅(qū)動因子對水田分布格局的相互作用,對于影響水田化進(jìn)程的2種因子1和2,計(jì)算他們交互(疊加變量1和2兩個(gè)圖層形成新的多邊形分布)時(shí)的值,對(1)、(2)與(1∩2)進(jìn)行比較。兩因子之間關(guān)系可分為以下幾類:若(1∩2) Max((1),(2)),說明因子1和2交互后雙因子增強(qiáng);若(1∩2)=(1)+(2),說明因子1和2交互后相互獨(dú)立;若(1∩2)>(1)+(2),說明因子1和2交互后非線性增強(qiáng)。

    本研究中應(yīng)用地理探測器模型分析水田分布格局的變化與所選7種影響因子的相關(guān)性以及多種影響因子之間的交互作用。地理探測器的具體實(shí)現(xiàn)采用GeoDetector軟件(http://www.geodetector.org/),讀入離散化處理后的類型數(shù)據(jù),軟件運(yùn)行結(jié)果本研究側(cè)重分析兩部分:影響因子對水田變化格局的空間分異性的解釋力,以及這些影響因子對水田變化的空間分異性影響的交互作用。

    3 結(jié)果與分析

    3.1 水田變化空間特征分析

    1989—2015年富錦市水田擴(kuò)張明顯,新增水田34.99萬hm2?;贏rcGIS 10.2平臺,利用自然間斷點(diǎn)分級法將1989—2002年和2002—2015年水田擴(kuò)張面積進(jìn)行分級,分類數(shù)為5。與1 000 m×1 000 m網(wǎng)格層疊加,得到1989—2002年和2002—2015年富錦市1 000 m×1 000 m網(wǎng)格單元統(tǒng)計(jì)的水田擴(kuò)張度空間分布圖(圖2)。

    1989—2002年,當(dāng)?shù)卣耙缘局螡场鞭r(nóng)業(yè)種植結(jié)構(gòu)調(diào)整政策的推行,導(dǎo)致富錦市水田擴(kuò)張度在空間上呈現(xiàn)大范圍連片擴(kuò)展和零星分散。1 000 m×1 000 m格網(wǎng)內(nèi)水田擴(kuò)張度高于0.3的網(wǎng)格單元數(shù)為2 210個(gè),占總數(shù)的26%,主要聚集分布在富錦市東部墾區(qū),為水田化集中分布區(qū)。水田擴(kuò)張度在閾值0~0.102 5區(qū)間內(nèi)的網(wǎng)格單元數(shù)為5 425個(gè),占總數(shù)的63%,主要分布在富錦市中西部地區(qū)。水田擴(kuò)張度在閾值0.102 6~0.3區(qū)間內(nèi)的網(wǎng)格單元數(shù)為976個(gè),占總數(shù)的11.33%,主要零星分布在富錦市中部平原區(qū)以及沼澤濕地邊緣地帶。在規(guī)模上,東部地區(qū)的水田擴(kuò)張度明顯高于中西部地區(qū),且水田擴(kuò)張聚集密度較高,以富錦市東部墾區(qū)農(nóng)場尤為顯著。此階段,富錦市中部傳統(tǒng)農(nóng)區(qū)經(jīng)濟(jì)發(fā)展滯后,主要以墾殖沼澤濕地為主,受地形、土質(zhì)等因素影響,水田化程度較低。

    圖2 富錦市1 000 m×1 000 m網(wǎng)格單元統(tǒng)計(jì)的水田擴(kuò)張度分布

    2002—2015年,“兩江一湖工程”、“高標(biāo)準(zhǔn)基本農(nóng)田建設(shè)工程”等耕地開發(fā)整治工程實(shí)施,富錦市耕地達(dá)到較高的農(nóng)業(yè)開發(fā)水平,水田化水平大幅度提升,有條件地區(qū)基本實(shí)現(xiàn)全域水田化。該時(shí)段內(nèi),1 000 m×1 000 m格網(wǎng)內(nèi)水田擴(kuò)張度在0~0.094 3區(qū)間內(nèi)的網(wǎng)格單元數(shù)為3 718個(gè),占總網(wǎng)格單元數(shù)的36.91%,區(qū)域內(nèi)大范圍耕地基本實(shí)現(xiàn)“旱改水”。水田擴(kuò)張度在閾值0.094 4~0.266 7區(qū)間內(nèi)網(wǎng)格單元數(shù)為1 913個(gè),在閾值0.266 8~0.472 3范圍內(nèi)的網(wǎng)格單元數(shù)為1 489個(gè),2個(gè)區(qū)間內(nèi)網(wǎng)格單元數(shù)占總數(shù)的39.51%,主要分布于富錦市墾區(qū)的中南部,由于墾區(qū)20世紀(jì)90年代以來率先開展“以稻治澇”實(shí)施旱地轉(zhuǎn)水田工程,進(jìn)入水田化成熟期,2002—2015年水田擴(kuò)張度相比前一階段有所降低。1 000 m×1 000 m格網(wǎng)內(nèi)水田擴(kuò)張度高于47.24%,處于第四、五區(qū)間內(nèi)的網(wǎng)格單元數(shù)為2 031個(gè),占網(wǎng)格單元總數(shù)的23.59%,這部分區(qū)域主要位于富錦市西部農(nóng)區(qū)的錦山鎮(zhèn)、長安鎮(zhèn)和硯山鎮(zhèn)等區(qū)域,表明富錦市農(nóng)區(qū)受體制等因素制約,旱田轉(zhuǎn)水田進(jìn)程慢于東部墾區(qū)。

    3.2 水田變化空間特征驅(qū)動因素分析

    探索不同影響因素對三江平原水田化的影響及作用強(qiáng)度一直是眾多學(xué)者關(guān)注的焦點(diǎn)。本文通過定量研究的方法分析造成該區(qū)域水田化區(qū)域差異的驅(qū)動因素,并探索這些驅(qū)動因子在地理空間上對水田空間格局變化的解釋強(qiáng)度。由于地理探測器擅長識別類型數(shù)據(jù),而對于順序量、比值量或間隔量,需要進(jìn)行適當(dāng)?shù)碾x散化[13],本文將高程、坡度、到河流距離、到鐵路距離與到城鎮(zhèn)距離都用自然斷點(diǎn)分成5類,土壤類型分為8類,政策因素按照農(nóng)區(qū)與墾區(qū)行政范圍將分別賦值為0與1,以對政策因素量化和空間化表達(dá),便于在地理探測器中進(jìn)行分析(圖3)。

    圖3 水田變化驅(qū)動因子空間分布圖

    3.2.1 驅(qū)動因子探測

    運(yùn)用GeoDetector軟件進(jìn)行數(shù)據(jù)處理,借助因子探測器模塊獲得1989—2002年和2002—2015年水田變化各驅(qū)動因子的作用強(qiáng)度值(地理探測器的值),值越大表示該指標(biāo)對水田擴(kuò)張度空間分異的作用強(qiáng)度越大,且是決定水田擴(kuò)張度空間分異的主要影響因素。

    1989—2002年水田變化的驅(qū)動因子作用強(qiáng)度值從大到小排列依次為:政策因素(0.252 4)、土壤類型(0.180 1)、到河流距離(0.155 8)、到鐵路距離(0.111 4)、高程(0.100 4)、坡度(0.0277 4)、到城鎮(zhèn)距離(0.155 8)。該時(shí)期富錦市水田面積的擴(kuò)張是多因子共同作用的結(jié)果,其中政策因素起主導(dǎo)作用,土壤類型、到河流距離、到鐵路距離也具有較強(qiáng)的解釋力,而高程、坡度、到城鎮(zhèn)距離對水田面積擴(kuò)張的作用強(qiáng)度較小。表明20世紀(jì)90年代富錦市水田面積擴(kuò)張?jiān)谵r(nóng)業(yè)結(jié)構(gòu)調(diào)整的政策引領(lǐng)下,率先改造土壤類型適合種植水田、區(qū)位條件優(yōu)越的土地,導(dǎo)致富錦市東部墾區(qū)大量開發(fā)條件優(yōu)越的耕地資源被墾殖為水田。

    2002—2015年水田變化的驅(qū)動因子作用強(qiáng)度值從大到小排列依次為:到河流距離(0.041 4)、到鐵路距離(0.028 9)、高程(0.018 7)、土壤類型(0.015 3)、到城鎮(zhèn)距離(0.015 2)、政策因素(0.005 2)、坡度(0.000 7)。中國城鎮(zhèn)化進(jìn)程加快導(dǎo)致南方水稻種植面積逐漸減少,國家對糧食需求量持續(xù)增加,促使北方水稻種植面積逐漸擴(kuò)大。三江平原地區(qū)水田種植區(qū)的水稻受到消費(fèi)者認(rèn)可,導(dǎo)致市場需求量逐年上升,在市場經(jīng)濟(jì)導(dǎo)向作用下,促進(jìn)富錦市水田種植面積進(jìn)一步擴(kuò)大。另一方面,種植水稻的平均收益遠(yuǎn)高于玉米的平均收益,在比較利益的驅(qū)動下,農(nóng)民在條件適宜的地區(qū)實(shí)行旱改水工程,導(dǎo)致大面積旱地轉(zhuǎn)為水田。該時(shí)期自然因素已不再是限制水田發(fā)展關(guān)鍵,利益驅(qū)動下滿足水田改造條件的耕地都轉(zhuǎn)為水田,富錦市基本實(shí)現(xiàn)全面水田化。

    3.2.2 驅(qū)動因子交互作用探測

    富錦市水田化空間分布差異是由多種影響因子共同作用的結(jié)果,不可能存在單一因子或單一性因素影響水田的分布和變化。借助地理探測器的交互作用探測器模塊獲得1989—2002年和2002—2015年各驅(qū)動因子對水田變化的交互作用(表1、表2)。

    表1 1989—2002年富錦市水田變化各驅(qū)動因子的交互作用

    注:1、2為驅(qū)動因子,Max為(1)、(2)中最大值,Min為(1)、(2)中最小值。下同。

    Note:1 and2 are the driving factors. Max is the maximum value of(1) and(2). Min is the minimum value of(1) and(2). The same below。

    1989—2002年間各驅(qū)動因子對水田變化的交互作用大部分是雙因子增強(qiáng)和非線性增強(qiáng),不存在相互獨(dú)立起作用的因子。其中解釋力較大的交互驅(qū)動因子土壤類型/政策因素(0.311 7)、到鐵路距離/政策因素(0.298 0)、到河流距離/政策因素(0.290 2)、土壤類型/到河流距離(0.282 6)、到城鎮(zhèn)距離/政策因素(0.262 5)等均為雙因子增強(qiáng)。土壤類型/到城鎮(zhèn)距離(0.211 6)、到城鎮(zhèn)距離/到河流距離(0.185 4)、到城鎮(zhèn)距離/到鐵路距離(0.165 5)等為非線性增強(qiáng)。政策因素和土壤類型可以較大增強(qiáng)各驅(qū)動因子對水田分布空間格局變化的解釋力,說明政策引導(dǎo)對水田的空間變化有一定影響。

    2002—2015年間各因子對水田變化的交互作用大部分(1∩2)>(1)+(2),是非線性增強(qiáng)的。其中解釋力較大的交互驅(qū)動因子從大到小依次為:到河流距離/到鐵路距離(0.097 9)、土壤類型/到河流距離(0.073 1)、高程/到鐵路距離(0.071 2)、土壤類型/到鐵路距離(0.067 7)、到城鎮(zhèn)距離/到鐵路距離(0.066 8)。相互作用增強(qiáng)的主要是自然因子與社會因子交互,自然因子與社會因子非線性增強(qiáng)的貢獻(xiàn)率大于兩者之和,說明該時(shí)期水田擴(kuò)張不僅僅是單因素作用,是自然因子與社會因子綜合作用的結(jié)果,社會因子與自然作用的交互作用造成富錦市水田的大面積擴(kuò)張。

    4 討 論

    本研究對富錦市水田分布的動態(tài)變化過程和水田變化的驅(qū)動因素兩個(gè)方面進(jìn)行研究。三江平原地勢平坦、土壤肥沃、水源豐富,自然條件適宜農(nóng)業(yè)發(fā)展,自20世紀(jì)90年代以來,三江平原大面積實(shí)施“旱改水”工程[26-27]。富錦市作為三江平原的核心區(qū),水田化程度較高,由于富錦市行政管轄復(fù)雜,包含農(nóng)區(qū)與墾區(qū),體制不同導(dǎo)致富錦市水田化進(jìn)程空間差異性明顯。本研究運(yùn)用網(wǎng)格單元法明確了水田變化的大小、幅度,在驅(qū)動因素研究中運(yùn)用地理探測器模型很好的探測出1989—2015年富錦市水田變化的主導(dǎo)驅(qū)動因子和因子間的交互作用。

    從研究結(jié)果可以看出,1989—2015年富錦市水田擴(kuò)張明顯,且墾區(qū)水田化進(jìn)程早于農(nóng)區(qū)。富錦市地勢低洼,經(jīng)常發(fā)生澇災(zāi),90年代末當(dāng)?shù)貕▍^(qū)因地制宜,率先提出“以稻治澇、以稻致富”,把地勢低平易澇的不利因素變成有利因素,積極推動種植結(jié)構(gòu)向水稻調(diào)整。1989—2015年,富錦市農(nóng)區(qū)與墾區(qū)的水田變化面積占比接近。1989—2002年,墾區(qū)的水田變化面積占富錦市水田變化總面積的73.09%,遠(yuǎn)高于農(nóng)區(qū)的36.81%;2002—2015年,富錦市農(nóng)區(qū)的水田變化面積占比高于墾區(qū),說明墾區(qū)的水田化進(jìn)程早于農(nóng)區(qū)(表3)。

    表3 1989—2015年各時(shí)期水田變化面積中農(nóng)區(qū)與墾區(qū)占比

    結(jié)果表明,土壤類型對富錦市水田變化的驅(qū)動力較大。分別測算富錦市2015年水田、旱地中各土壤類型占比,水田中所占比重較高的土壤類型有草甸土、白漿土和沼澤土,旱地主要土壤類型為草甸土、沼澤土和黑土(表4)。富錦市水田中白漿土占比32.01%,遠(yuǎn)高于旱地的12.72%,白漿土質(zhì)地粘重,透水性差,低平地區(qū)的白漿土適合改種水田[28],黑土的自然肥力很高,但可耕性差,水分有效性低[29],種植水田需要的灌溉成本過高,因此富錦市旱地中黑土占比21.43%,高于水田的8.70%。

    此外,受衛(wèi)星影像分辨率和分類精度影響,結(jié)果存在一定誤差,今后應(yīng)該考慮光學(xué)影像與雷達(dá)影像相結(jié)合,并運(yùn)用Google Earth Engine云平臺等先進(jìn)手段進(jìn)行水田提取[30-31]。此外,本研究分析富錦市水田變化的空間特征及驅(qū)動因素,下一步應(yīng)結(jié)合耕地保護(hù)及農(nóng)業(yè)種植結(jié)構(gòu)調(diào)整政策,提出耕地利用的調(diào)整優(yōu)化方向,結(jié)合不同發(fā)展情景進(jìn)行模擬。

    表4 2015年水田、旱地和耕地中各土壤類型占比

    5 結(jié) 論

    本文以富錦市水田為研究對象,以遙感影像、DEM數(shù)據(jù)、土壤數(shù)據(jù)等數(shù)據(jù)為基本信息源,基于RS和GIS技術(shù),運(yùn)用地理探測器模型,對研究區(qū)1989—2015年水田變化特征及驅(qū)動因素進(jìn)行研究,揭示當(dāng)?shù)厮镒兓臻g差異的驅(qū)動因子及因子間的交互作用,主要得出以下結(jié)論:

    1)通過1 000 m×1 000 m網(wǎng)格分析富錦市1989—2002年與2002—2015年水田變化圖,分析可得,1989—2002年水田擴(kuò)張度較高的區(qū)域主要分布在富錦市東部墾區(qū),中西部農(nóng)區(qū)水田擴(kuò)張度較低;2002—2015年,富錦市大范圍耕地實(shí)施“旱改水”工程,該時(shí)期墾區(qū)在已有水田化基礎(chǔ)上繼續(xù)擴(kuò)張,西部農(nóng)區(qū)水田擴(kuò)張度高于東部墾區(qū)。

    2)1989—2002年水田擴(kuò)張主要驅(qū)動因子為:政策因素>土壤類型>到河流距離>到鐵路距離>高程>坡度>到城鎮(zhèn)距離,其中政策因素和土壤類型起主導(dǎo)作用。2002—2015年水田擴(kuò)張主要驅(qū)動因子為:到河流距離>到鐵路距離>高程>土壤類型>到城鎮(zhèn)距離>政策因素>坡度,該時(shí)期在利益驅(qū)動下富錦大部分旱地轉(zhuǎn)為水田,驅(qū)動因子作用力均衡。

    3)1989—2002年各因子間的交互作用大部分是雙因子增強(qiáng)的,解釋力較大的驅(qū)動因子有土壤類型/政策因素、到鐵路距離/政策因素、到河流距離/政策因素等;2002—2015年各因子交互作用大部分是非線性增強(qiáng)的,對水田擴(kuò)張解釋力較大的因子主要有到河流距離/到鐵路距離、土壤類型/到河流距離、高程/到鐵路距離等。

    [1] Ge D, Long H, Zhang Y, et al. Farmland transition and its influences on grain production in China[J]. Land Use Policy, 2018, 70: 94-105.

    [2] Vliet J, Groot H L F, Rietveld P, et al. Manifestations and underlying drivers of agricultural land use change in Europe[J]. Landscape and Urban Planning, 2015, 133: 24-36.

    [3] Song Jie, Cai Dongyan, Deng Jinsong, et al. Dynamics of paddy field patterns in response to urbanization: A case study of the Hang-Jia-Hu Plain[J]. Sustainability, 2015, 7: 13813-13835.

    [4] 吳莉,侯西勇,徐新良. 環(huán)渤海沿海區(qū)域耕地格局及影響因子分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2014,30 (9):1-10. Wu Li, Hou Xiyong, Xu Xinliang. Analysis of spatial pattern of farmland and its impacting factors in coastal zone of Circum Bohai[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(9): 1-10. (in Chinese with English abstract)

    [5] Yan Fengqin, Yu Lingxue, Yang Chaobin, et al. Paddy field expansion and aggregation since the mid-1950s in a cold region and its possible causes[J]. Remote sensing, 2018, 10: 384-401.

    [6] 宋戈,楊雪昕,高佳. 三江平原典型地區(qū)水田分布格局變化特征[J]. 中國土地科學(xué),2017,31(8):61-68. Song Ge, Yang Xuexin, Gao Jia. Study on the distribution patterns and characteristics of Paddy cropland in the typical area of Sanjiang Plain[J]. China Land Sciences, 2017, 31(8): 61-68. (in Chinese with English abstract)

    [7] 姜秋香,付強(qiáng),王子龍,等. 三江平原水土資源空間匹配格局[J]. 自然資源學(xué)報(bào),2011,26(2):270-277. Jiang Qiuxiang, Fu Qiang, Wang Zilong, et al. Spatial matching patterns of land and water resources in Sanjiang Plain[J]. Journal of Natural Resources, 2011, 26(2): 270-277. (in Chinese with English abstract)

    [8] Jin Jianjun, Jiang Chong, Li Lun. The economic valuation of cultivated land protection: A contingent valuation study in Wenling City, China[J]. Landscape and Urban Planning, 2013, 119: 158-164.

    [9] 郇紅艷,譚清美,朱平. 城鄉(xiāng)一體化進(jìn)程中耕地利用變化的驅(qū)動因素及區(qū)域比較[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(21):201-213. Huan Hongyan, Tan Qingmei, Zhu Ping. Driving factors of cultivated land use change in city and countryside integration process and region comparisons[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(21): 201-213. (in Chinese with English abstract)

    [10] 歐陽玲,王宗明,賈明明,等. 基于遙感的吉林省中西部耕地?cái)?shù)量和質(zhì)量空間格局變化分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(13):234-242. Ouyang Ling, Wang Zongming, Jia Mingming, et al. Spatial pattern analysis on quantity and quality of cultivated land in mid-west Jilin province based on remote sensing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(13): 234-242. (in Chinese with English abstract)

    [11] 劉敏,趙翠薇,施明輝. 貴州山區(qū)土地利用變化多尺度空間自相關(guān)分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2012,28(20):239-246. Liu Min, Zhao Cuiwei, Shi Minghui. Spatial autocorrelation analysis of multi-scale land use change at mountainous areasin Guizhou province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(20): 239-246. (in Chinese with English abstract)

    [12] 劉艷軍,于會勝,劉德剛,等. 東北地區(qū)建設(shè)用地開發(fā)強(qiáng)度格局演變的空間分異機(jī)制[J]. 地理學(xué)報(bào),2018,73(5):818-831. Liu Yanjun, Yu Huisheng, Liu Degang, et al. Spatial differentiation mechanisms of the pattern evolution of construction land development intensity in Northeast China[J]. Acta Geographica Sinica, 2018, 73(5): 818-831. (in Chinese with English abstract)

    [13] 王勁峰,徐成東. 地理探測器:原理與展望[J]. 地理學(xué)報(bào),2017,72(1):116-134. Wang Jinfeng, Xu Chengdong. Geodetector: Principle and prospective[J]. Acta Geographica Sinica, 2017, 72(1): 116-134. (in Chinese with English abstract)

    [14] Wang J F, Li X H, Christakos G, et al. Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun region, China[J]. International Journal of Geographical Information Science. 2010, 24(1): 107-127.

    [15] 李娜,雷國平,張慧,等. 水田化進(jìn)程下?lián)狭恿饔蚋貢r(shí)空變化特征[J]. 水土保持研究,2016,23(5):63-67+73. Li Na, Lei Guoping, Zhang Hui, et al. Spatial-temporal characteristics of farmland due to paddy field expansion in Naolihe River Basin[J]. Research of Soil and Water Conservation, 2016, 23(5): 63-67+73. (in Chinese with English abstract)

    [16] 朱宇. 黑龍江屯墾史[M]. 北京:社會科學(xué)文獻(xiàn)出版社,2017.

    [17] 劉煥軍,于勝男,張新樂,等. 一年一季農(nóng)作物遙感分類的時(shí)效性分析[J]. 中國農(nóng)業(yè)科學(xué),2017,50(5):830-839. Liu Huanjun, Yu Shengnan, Zhang Xinle, et al. Timeliness analysis of crop remote sensing classification one crop a year[J]. Scientia Agricultura Sinica, 2017, 50(5): 830-839. (in Chinese with English abstract)

    [18] 劉吉平,馬長迪,劉雁,等. 基于地理探測器的沼澤濕地變化驅(qū)動因子定量分析:以小三江平原為例[J]. 東北師大學(xué)報(bào):自然科學(xué)版,2017,49(2):127-135. Liu Jiping, Ma Changdi, Liu Yan, et al. Quantitative study on the driving factors of marsh change based in Geographical Dector: Case study on Small Sanjiang Plain[J]. Journal of Northeast Normal University: Natural Science Edition, 2017, 49(2):127-135. (in Chinese with English abstract)

    [19] 胡金曉,朱明,濮勵杰,等. 交通道路對區(qū)域土地利用變化的影響研究進(jìn)展[J]. 長江流域資源與環(huán)境,2018,27(1):205-214. Hu Jinxiao, Zhu Ming, Pu Lijie, et al. Progress of the research on the influence of traffic road on regional land use change[J]. Resources and Environment in the Yangtze Basin, 2018, 27(1): 205-214. (in Chinese with English abstract)

    [20] 劉超,霍永偉,許月卿,等. 生態(tài)退耕前后張家口市耕地變化及影響因素識別[J]. 自然資源學(xué)報(bào),2018,33(10):1806-1820. Liu Chao, Huo Yongwei, Xu Yueqing, et al. Changes in cultivated land and influencing factors before and after the implementation of grain for green project in Zhangjiakou City[J]. Journal of Natural Resources, 2018, 33(10): 1806-1820. (in Chinese with English abstract)

    [21] 買買提江·買提尼亞孜,阿里木江·卡斯木. 基于網(wǎng)格單元的烏魯木齊市土地覆被/利用時(shí)空變化[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(1):210-216. Maimaitijiang Maitiniyazi, Alimujiang Kasimu. Spatial-temporal change of Urumqi urban land use and land cover based on grid cell approach[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(1): 210-216. (in Chinese with English abstract)

    [22] 楊忍,劉彥隨,龍花樓,等. 基于格網(wǎng)的農(nóng)村居民點(diǎn)用地時(shí)空特征及空間指向性的地理要素識別:以環(huán)渤海地區(qū)為例[J]. 地理研究,2015,34(6):1077-1087. Yang Ren, Liu Yansui, Long Hualou. et al. patial-temporal characteristics of rural residential land usechange and spatial directivity identification based on gridin the Bohai Rim in China[J]. Geographical Research, 2015, 34(6): 1077-1087. (in Chinese with English abstract)

    [23] 趙小風(fēng),劉夢麗,王柏源. 基于地理探測器的城市建設(shè)用地規(guī)模時(shí)空差異及影響因素分析[J]. 中國土地科學(xué),2018,32(5):29-35. Zhao Xiaofeng, Liu Mengli, Wang Baiyuan. Spatio temporal disparity of urban construction land and its impact factors in China based on Geodetector[J]. China Land Sciences, 2018, 32(5): 29-35. (in Chinese with English abstract)

    [24] 李穎,馮玉,彭飛,等. 基于地理探測器的天津市生態(tài)用地格局演變[J]. 經(jīng)濟(jì)地理,2017,37(12):180-189. Li Ying, Feng Yu, Peng Fei, et al. Pattern evolvement of ecological land in Tianjin based on geodetector[J]. Economic Geography, 2017, 37(12): 180-189. (in Chinese with English abstract)

    [25] 董玉祥,徐茜,楊忍,等. 基于地理探測器的中國陸地?zé)釒П苯缣接慬J]. 地理學(xué)報(bào),2017,72(1):135-147. Dong Yuxiang, Xu Qian, Yang Ren, et al. Delineation of the northern border of the tropical zone of China's mainland using Geodetector[J]. Acta Geographica Sinica, 2007, 72(1): 135-147. (in Chinese with English abstract)

    [26] 黃妮,劉殿偉,王宗明. 1986年~2005年三江平原水田與旱地的轉(zhuǎn)化特征[J]. 資源科學(xué),2009,31(2):324-329. Hang Ni, Liu Dianwei, Wang Zongming. Study on mutual transformation characteristics between paddy field and dry land in Sanjiang Plain from 1986 to 2005[J]. Resources Science, 2009, 31(2): 324-329. (in Chinese with English abstract)

    [27] 杜國明,潘濤,尹哲睿,等. 水田化進(jìn)程中的富錦市耕地景觀格局演化規(guī)律[J]. 應(yīng)用生態(tài)學(xué)報(bào),2015,26(1):207-214. Du Guoming, Pan Tao, Yin Zherui, et al. Cultivated landscape pattern change due to the rice paddy expansion in Northeast China: A case study in Fujin[J]. Chinese Journal of Applied Ecology, 2015, 26(1): 207-214. (in Chinese with English abstract)

    [28] 王晉,莊舜堯,朱兆良. 不同種植年限水田與旱地土壤有機(jī)氮組分變化[J]. 土壤學(xué)報(bào),2014,51(2):286-294. Wang Jin, Zhuang Shunyao, Zhu Zhaoliang. Fractions of soil organic nitrogen in paddyand upland soils relative to cropping history[J]. Acta Pedologica Sinica, 2014, 51(2): 286-294. (in Chinese with English abstract)

    [29] 王嬡華,蘇以榮,李楊,等. 水田和旱地土壤有機(jī)碳周轉(zhuǎn)對水分的響應(yīng)[J]. 中國農(nóng)業(yè)科學(xué),2012,45(2):266-274. Wang Aihua, Su Yirong, Li Yang, et al. Response of the turnover of soil organic carbon to the soil moisture in paddy and uplandsoil[J]. Scientia Agricultura Sinica, 2012, 45 (2): 266-274. (in Chinese with English abstract)

    [30] Dong J, Xiao X, Menarguez M A, et al. Mapping paddy rice planting area in northeastern Asia with Landsat 8 images, phenology-based algorithm and Google Earth Engine[J]. Remote Sensing of Environment, 2016, 185: 142-154.

    [31] Yang H, Pan B, Wu W, et al. Field-based rice classification in Wuhua county through integration of multi-temporal Sentinel-1A and Landsat-8 OLI data[J]. International Journal of Applied Earth Observation & Geoinformation, 2018, 69: 226-236.

    Spatial-temporal variations and driving factor analysis of paddy fields in typical regions of Sanjiang Plain

    Zhang Wenqi, Song Ge※

    (,,110169,)

    Over the last three decades, paddy fields on the Sanjiang Plain have experienced rapid expansion. It has an important impact on regional food security, ecological environment and socio-economic development. Fujin is a typical area of Sanjiang Plain, with both common agricultural area and reclamation area, and the land use pattern is the epitome of Sanjiang Plain. In our study, used 1 000 m×1 000 m grid cell to scientifically analyze the spatial-temporal variation and driving factors of paddy fields, in Fujin from 1989 to 2002 and from 2002 to 2015. And elevation, slope, soil type, distance to railway, distance to river, distance to town and policy factors were selected as driving factors in this paper. With the help of GeoDetector, we analyzed the role of these seven factors in paddy field change and their relationship. The results indicated that such paddy fields have increased by approximately 349 900 hm2from 1989 to 2015, and the paddy field process of reclamation area was earlier than that of common agricultural area. In the early 1990s, in order to adjust the agricultural planting structure, the local government began to encourage farmers to plant paddy rice in large areas. During this period, paddy rice planting area increased rapidly. From 1989 to 2002, the areas with the high proportion of paddy field expansion were mainly located in the reclamation areas in the northeast and southeast of Fujin County, and some other agricultural areas. From 2002 to 2015, the project of cultivated land development and improvement was implemented, and the cultivated land in Fujin reached a high level of agricultural development, and all the conditional areas were basically converted to paddy fields. The areas with higher proportion of paddy field expansion were mainly located in the southwestern common agricultural area and the northern reclamation area of Fujin City. In addition, the main driving factors of paddy field expansion in 1989-2002 were policy factors and soil types. The interaction between these factors was mostly enhanced by two factors. It showed that in the 1990s, under the guidance of the policy of agricultural structure adjustment, Fujin took the lead in reforming the soil types suitable for paddy fields and the land with superior location conditions. As a result, a large number of cultivated land resources with superior development conditions in the eastern reclamation area of Fujin were reclaimed into paddy fields. From 2002 to 2015, the influence of natural factors on paddy field change becomes weaker. The main driving factors of paddy field expansion are the distance from rivers to railways, and the interaction among these factors is mostly non-linear. Natural factors are no longer the key to restrict the development of paddy fields. Driven by market interests, dryland that meets the conditions of paddy field transformation is converted to paddy fields. The spatial difference of paddy field change in Sanjiang Plain is the result of the interaction of various driving factors. The driving factors and their relationships are different in different periods. Geographic detectors can make a more comprehensive analysis of various influencing factors. The results can provide theoretical basis for the protection of cultivated land and the rational utilization of cultivated land resources in Sanjiang Plain.

    land use; remote sensing; spatial-temporal change; cultivated land; paddy field; driving factors; GeoDetector; Sanjiang plain

    2018-09-03

    2018-12-12

    國家自然科學(xué)基金資助項(xiàng)目(41571165,41071346)

    張文琦,博士生,主要研究方向?yàn)橥恋乩门c管理。Email:wenqi9094@163.com

    宋 戈,教授,博士生導(dǎo)師,主要研究方向?yàn)橥恋乩门c管理。Email:songgelaoshi@163.com

    10.11975/j.issn.1002-6819.2019.06.030

    F301.21

    A

    1002-6819(2019)-06-0244-9

    張文琦,宋 戈. 三江平原典型區(qū)水田時(shí)空變化及驅(qū)動因素分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(6):244-252. doi:10.11975/j.issn.1002-6819.2019.06.030 http://www.tcsae.org

    Zhang Wenqi, Song Ge. Spatial-temporal variations and driving factor analysis of paddy fields in typical regions of Sanjiang Plain[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(6): 244-252. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2019.06.030 http://www.tcsae.org

    猜你喜歡
    富錦市墾區(qū)水田
    黑龍江富錦:濕地風(fēng)景秀麗如畫
    先鋒廈地水田書店
    基于偏離—份額法的黑龍江省富錦市產(chǎn)業(yè)結(jié)構(gòu)變遷與競爭力分析
    基于能值理論的富錦市縣域生態(tài)經(jīng)濟(jì)系統(tǒng)研究
    哈密墾區(qū)棉花超寬膜1膜3行稀植高產(chǎn)栽培技術(shù)
    近世長三角南緣的水田經(jīng)營與環(huán)境博弈
    鳳頭豬肚豹尾說“白傳”——讀《白水田傳》
    新聞前哨(2015年2期)2015-03-11 19:29:30
    淺談富錦市寒地水稻常見病害及防治方法
    湖南省國土資源廳加強(qiáng)耕地占補(bǔ)平衡管理
    博樂墾區(qū)引進(jìn)黑綠豆種植試驗(yàn)
    亚洲色图综合在线观看| 欧美人与性动交α欧美软件| 久久久久国内视频| 十八禁高潮呻吟视频| 久久精品国产a三级三级三级| 999久久久精品免费观看国产| 亚洲熟女毛片儿| 久久精品亚洲熟妇少妇任你| 午夜日韩欧美国产| 国产精品一区二区在线不卡| av天堂久久9| 亚洲成国产人片在线观看| 免费不卡黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 两性夫妻黄色片| 午夜影院日韩av| 九色亚洲精品在线播放| 国产成+人综合+亚洲专区| 精品卡一卡二卡四卡免费| 最新在线观看一区二区三区| 一进一出抽搐动态| 在线观看午夜福利视频| 如日韩欧美国产精品一区二区三区| 最近最新免费中文字幕在线| 在线天堂中文资源库| 国产精品一区二区在线观看99| 香蕉丝袜av| 亚洲,欧美精品.| 十八禁网站免费在线| 欧美日本中文国产一区发布| 国产精品九九99| 91字幕亚洲| 午夜福利视频在线观看免费| 国产精品1区2区在线观看. | 国产成人欧美在线观看 | 亚洲美女黄片视频| 精品人妻1区二区| 国产又色又爽无遮挡免费看| 亚洲国产毛片av蜜桃av| 啦啦啦免费观看视频1| 女人被躁到高潮嗷嗷叫费观| 成年版毛片免费区| 韩国av一区二区三区四区| 看片在线看免费视频| 成人18禁在线播放| 国产成人av教育| 视频在线观看一区二区三区| 亚洲精品粉嫩美女一区| 黄片小视频在线播放| 精品国产美女av久久久久小说| 午夜免费观看网址| 亚洲av美国av| 亚洲精品美女久久av网站| 少妇被粗大的猛进出69影院| 日本五十路高清| 精品亚洲成国产av| 国产男靠女视频免费网站| 国产日韩欧美亚洲二区| 国产又爽黄色视频| 精品国产一区二区三区四区第35| 国产高清激情床上av| 波多野结衣av一区二区av| 少妇被粗大的猛进出69影院| 日本精品一区二区三区蜜桃| 欧美日韩中文字幕国产精品一区二区三区 | 日韩大码丰满熟妇| 91av网站免费观看| 精品国产美女av久久久久小说| 一二三四社区在线视频社区8| 女性生殖器流出的白浆| 国产精品自产拍在线观看55亚洲 | 精品久久久久久久久久免费视频 | 老司机午夜十八禁免费视频| 亚洲一区高清亚洲精品| 国产精品永久免费网站| 国产精品一区二区免费欧美| 国产成人免费无遮挡视频| 美女福利国产在线| 欧美性长视频在线观看| 亚洲欧美精品综合一区二区三区| 免费在线观看影片大全网站| 三上悠亚av全集在线观看| tube8黄色片| 国产精品免费一区二区三区在线 | 国产精品成人在线| 成人黄色视频免费在线看| 成在线人永久免费视频| 久久亚洲真实| 嫁个100分男人电影在线观看| 99精品久久久久人妻精品| 无人区码免费观看不卡| 婷婷丁香在线五月| 久久人妻av系列| 人妻丰满熟妇av一区二区三区 | 久久精品亚洲精品国产色婷小说| 国产精品国产高清国产av | 天天添夜夜摸| 国产成人精品无人区| 午夜福利在线观看吧| 欧美乱妇无乱码| 国产国语露脸激情在线看| 成年女人毛片免费观看观看9 | 又黄又粗又硬又大视频| 建设人人有责人人尽责人人享有的| 亚洲欧美日韩高清在线视频| 久久久久精品国产欧美久久久| 黄色a级毛片大全视频| 亚洲精品久久午夜乱码| 国产精品久久久久成人av| 国产精品久久久久久精品古装| 精品亚洲成a人片在线观看| 国产精品国产av在线观看| 黄色女人牲交| 黄色视频不卡| 日韩欧美免费精品| av超薄肉色丝袜交足视频| 久久天躁狠狠躁夜夜2o2o| 欧美亚洲 丝袜 人妻 在线| 国产亚洲欧美98| 伦理电影免费视频| 亚洲自偷自拍图片 自拍| 国产精品国产高清国产av | 两性夫妻黄色片| 在线av久久热| 黄色怎么调成土黄色| 好看av亚洲va欧美ⅴa在| 曰老女人黄片| 亚洲欧美激情在线| 色在线成人网| 久久天堂一区二区三区四区| 国产av一区二区精品久久| 成年人免费黄色播放视频| 男女之事视频高清在线观看| 后天国语完整版免费观看| 色综合欧美亚洲国产小说| 免费一级毛片在线播放高清视频 | 国产伦人伦偷精品视频| 精品国产乱码久久久久久男人| 超色免费av| 欧美黄色淫秽网站| 美女扒开内裤让男人捅视频| 欧美国产精品一级二级三级| 久久精品国产综合久久久| 少妇 在线观看| 老熟妇仑乱视频hdxx| 日韩 欧美 亚洲 中文字幕| 亚洲精品美女久久久久99蜜臀| 一夜夜www| 国产色视频综合| 精品欧美一区二区三区在线| 成人18禁高潮啪啪吃奶动态图| 精品久久久精品久久久| 黑人巨大精品欧美一区二区mp4| av天堂在线播放| 十八禁网站免费在线| 夜夜爽天天搞| 99热网站在线观看| 国产成人av激情在线播放| 免费日韩欧美在线观看| 日韩制服丝袜自拍偷拍| 黄色女人牲交| 18禁美女被吸乳视频| 一进一出抽搐动态| 男女高潮啪啪啪动态图| 亚洲av成人一区二区三| 在线观看www视频免费| 极品少妇高潮喷水抽搐| a级片在线免费高清观看视频| 天天添夜夜摸| 男人操女人黄网站| 色综合欧美亚洲国产小说| 黑人操中国人逼视频| 黄网站色视频无遮挡免费观看| 欧美精品高潮呻吟av久久| 亚洲熟妇熟女久久| 亚洲一区二区三区欧美精品| 国产精品秋霞免费鲁丝片| 成人18禁在线播放| 国产成+人综合+亚洲专区| 久久人人97超碰香蕉20202| 成年动漫av网址| 日韩欧美免费精品| 久久香蕉激情| 国产人伦9x9x在线观看| 欧美日韩瑟瑟在线播放| 悠悠久久av| a级片在线免费高清观看视频| 国产成人精品久久二区二区91| 欧美人与性动交α欧美精品济南到| 成熟少妇高潮喷水视频| 王馨瑶露胸无遮挡在线观看| 精品少妇一区二区三区视频日本电影| 午夜福利视频在线观看免费| videosex国产| 大香蕉久久成人网| 久久久久国产一级毛片高清牌| 中文字幕高清在线视频| 亚洲精品在线美女| 久久精品国产亚洲av香蕉五月 | 成人影院久久| 天天躁日日躁夜夜躁夜夜| 午夜两性在线视频| 午夜福利欧美成人| 大型av网站在线播放| e午夜精品久久久久久久| av国产精品久久久久影院| 女同久久另类99精品国产91| 精品一区二区三卡| 免费高清在线观看日韩| 女警被强在线播放| 免费在线观看日本一区| 女性被躁到高潮视频| 麻豆成人av在线观看| 又紧又爽又黄一区二区| 欧美av亚洲av综合av国产av| 90打野战视频偷拍视频| 欧美精品一区二区免费开放| 黑人猛操日本美女一级片| 一本大道久久a久久精品| 亚洲国产欧美日韩在线播放| 日韩欧美一区视频在线观看| 男女免费视频国产| 操美女的视频在线观看| 久久 成人 亚洲| 精品国产亚洲在线| 午夜福利一区二区在线看| 亚洲人成77777在线视频| 久99久视频精品免费| 50天的宝宝边吃奶边哭怎么回事| 亚洲三区欧美一区| 80岁老熟妇乱子伦牲交| 久久久久久久精品吃奶| 国产成人啪精品午夜网站| 免费日韩欧美在线观看| 国产野战对白在线观看| 人妻一区二区av| 中亚洲国语对白在线视频| 成人三级做爰电影| 一本大道久久a久久精品| 亚洲专区中文字幕在线| 中文字幕人妻丝袜一区二区| 最新的欧美精品一区二区| 国产高清国产精品国产三级| 老司机在亚洲福利影院| 国产真人三级小视频在线观看| 老司机影院毛片| 国产极品粉嫩免费观看在线| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区精品视频观看| 国产精品香港三级国产av潘金莲| 大陆偷拍与自拍| 日本wwww免费看| 国产激情欧美一区二区| 日韩免费av在线播放| 99在线人妻在线中文字幕 | 一二三四在线观看免费中文在| 亚洲成人免费av在线播放| 99国产综合亚洲精品| 人人妻人人澡人人看| 9191精品国产免费久久| 三上悠亚av全集在线观看| 捣出白浆h1v1| 亚洲熟女精品中文字幕| 免费在线观看日本一区| 丝袜在线中文字幕| 精品人妻熟女毛片av久久网站| 国产色视频综合| 怎么达到女性高潮| 不卡一级毛片| 99国产精品免费福利视频| 在线观看免费日韩欧美大片| 我的亚洲天堂| 飞空精品影院首页| 五月开心婷婷网| 亚洲熟女精品中文字幕| 午夜免费鲁丝| 久久婷婷成人综合色麻豆| 亚洲熟妇熟女久久| bbb黄色大片| 在线国产一区二区在线| 新久久久久国产一级毛片| 成人亚洲精品一区在线观看| 久9热在线精品视频| 欧美亚洲 丝袜 人妻 在线| 999久久久国产精品视频| 丰满人妻熟妇乱又伦精品不卡| 在线免费观看的www视频| 中文字幕av电影在线播放| 午夜免费鲁丝| 欧美精品啪啪一区二区三区| 久久国产乱子伦精品免费另类| 国产精品av久久久久免费| 女人被狂操c到高潮| 亚洲一区中文字幕在线| 日韩欧美一区视频在线观看| 国产精品 国内视频| 亚洲情色 制服丝袜| 国产精品综合久久久久久久免费 | 中文字幕色久视频| 亚洲专区国产一区二区| 日韩欧美三级三区| 制服诱惑二区| 中文字幕av电影在线播放| 国产在线一区二区三区精| 极品教师在线免费播放| 色综合欧美亚洲国产小说| 久热爱精品视频在线9| 99热只有精品国产| 女警被强在线播放| 好看av亚洲va欧美ⅴa在| 亚洲第一av免费看| 亚洲精品中文字幕一二三四区| 欧美中文综合在线视频| 午夜两性在线视频| 中文字幕人妻丝袜一区二区| 午夜精品在线福利| 精品卡一卡二卡四卡免费| 啦啦啦免费观看视频1| 在线观看免费日韩欧美大片| 中文字幕av电影在线播放| 无遮挡黄片免费观看| 男女下面插进去视频免费观看| 国产片内射在线| 天堂俺去俺来也www色官网| 深夜精品福利| 久久这里只有精品19| 9191精品国产免费久久| 日本五十路高清| 日本a在线网址| 色尼玛亚洲综合影院| 久久久精品区二区三区| 19禁男女啪啪无遮挡网站| e午夜精品久久久久久久| 757午夜福利合集在线观看| 国产精品久久久人人做人人爽| 久久午夜亚洲精品久久| 91精品国产国语对白视频| 丝袜在线中文字幕| 女人久久www免费人成看片| 久久精品国产99精品国产亚洲性色 | 无人区码免费观看不卡| 亚洲成人手机| 日本欧美视频一区| 高清毛片免费观看视频网站 | 亚洲人成电影观看| xxxhd国产人妻xxx| 色在线成人网| 一级片免费观看大全| 大片电影免费在线观看免费| 曰老女人黄片| 天堂动漫精品| 一进一出抽搐动态| 99国产精品一区二区三区| a级片在线免费高清观看视频| 日本一区二区免费在线视频| 亚洲va日本ⅴa欧美va伊人久久| 国产欧美日韩一区二区三| 免费在线观看黄色视频的| 精品久久久精品久久久| 一区福利在线观看| 午夜精品在线福利| 精品高清国产在线一区| 日本欧美视频一区| 大香蕉久久成人网| 香蕉久久夜色| 一区福利在线观看| 桃红色精品国产亚洲av| 啦啦啦视频在线资源免费观看| 波多野结衣av一区二区av| 日韩大码丰满熟妇| 中文字幕av电影在线播放| 巨乳人妻的诱惑在线观看| 中国美女看黄片| 怎么达到女性高潮| 久久国产乱子伦精品免费另类| av超薄肉色丝袜交足视频| 欧美最黄视频在线播放免费 | 国产亚洲精品一区二区www | 可以免费在线观看a视频的电影网站| 黑人巨大精品欧美一区二区mp4| 国产一区二区三区综合在线观看| 亚洲成国产人片在线观看| 亚洲中文字幕日韩| 成人免费观看视频高清| 欧美日韩亚洲国产一区二区在线观看 | www日本在线高清视频| 国产不卡av网站在线观看| 国产成人av教育| 久久精品亚洲精品国产色婷小说| 久久久精品免费免费高清| 久99久视频精品免费| 成人永久免费在线观看视频| 在线永久观看黄色视频| 亚洲人成伊人成综合网2020| 变态另类成人亚洲欧美熟女 | av电影中文网址| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 欧美国产精品va在线观看不卡| 国产精品二区激情视频| 国产成人免费观看mmmm| 在线观看免费视频日本深夜| 黑丝袜美女国产一区| 法律面前人人平等表现在哪些方面| 午夜视频精品福利| av网站免费在线观看视频| 免费一级毛片在线播放高清视频 | 老司机午夜十八禁免费视频| 下体分泌物呈黄色| 视频区欧美日本亚洲| 国产男女内射视频| 日韩免费av在线播放| 久久久久久久久久久久大奶| 看黄色毛片网站| av欧美777| x7x7x7水蜜桃| 黄色成人免费大全| 亚洲av电影在线进入| 国产男女内射视频| 国产成人欧美| 最近最新中文字幕大全免费视频| 人妻丰满熟妇av一区二区三区 | 国产午夜精品久久久久久| 老司机午夜福利在线观看视频| 热re99久久精品国产66热6| 俄罗斯特黄特色一大片| 亚洲 欧美一区二区三区| 午夜福利视频在线观看免费| 久久性视频一级片| 成年人午夜在线观看视频| 精品人妻熟女毛片av久久网站| 亚洲综合色网址| ponron亚洲| 母亲3免费完整高清在线观看| 久久精品国产综合久久久| 亚洲av美国av| 亚洲免费av在线视频| 国产亚洲精品久久久久久毛片 | av网站免费在线观看视频| 欧美精品一区二区免费开放| 九色亚洲精品在线播放| 少妇被粗大的猛进出69影院| 国产精品美女特级片免费视频播放器 | 久久国产亚洲av麻豆专区| 天天影视国产精品| 母亲3免费完整高清在线观看| 国产精品永久免费网站| 99国产极品粉嫩在线观看| 国产成人精品在线电影| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三区在线| 久久精品成人免费网站| 国产1区2区3区精品| 老熟妇仑乱视频hdxx| 欧美在线一区亚洲| 在线观看舔阴道视频| 老汉色av国产亚洲站长工具| 午夜久久久在线观看| av线在线观看网站| 男人操女人黄网站| 亚洲精品成人av观看孕妇| 免费女性裸体啪啪无遮挡网站| 窝窝影院91人妻| 国产欧美日韩一区二区三| 色综合欧美亚洲国产小说| 亚洲精品在线观看二区| 99香蕉大伊视频| 亚洲欧美日韩另类电影网站| 精品免费久久久久久久清纯 | 亚洲熟妇中文字幕五十中出 | 国产精品一区二区在线不卡| 午夜影院日韩av| 操出白浆在线播放| 精品一区二区三区视频在线观看免费 | 精品国产国语对白av| 99re6热这里在线精品视频| 在线播放国产精品三级| 国产又爽黄色视频| 99久久国产精品久久久| 丰满饥渴人妻一区二区三| 国产成人精品久久二区二区91| 看免费av毛片| 国产精品免费一区二区三区在线 | av不卡在线播放| 国产精品自产拍在线观看55亚洲 | 精品人妻熟女毛片av久久网站| 国产精品欧美亚洲77777| 亚洲欧美色中文字幕在线| 99香蕉大伊视频| 身体一侧抽搐| 超色免费av| 一边摸一边做爽爽视频免费| 中亚洲国语对白在线视频| 成年版毛片免费区| 男女免费视频国产| 国产免费男女视频| 成在线人永久免费视频| 少妇猛男粗大的猛烈进出视频| 日本精品一区二区三区蜜桃| 国产高清国产精品国产三级| 80岁老熟妇乱子伦牲交| 欧美一级毛片孕妇| 精品久久久久久,| 黄色片一级片一级黄色片| 老司机靠b影院| 欧美日韩视频精品一区| 欧美精品亚洲一区二区| 免费人成视频x8x8入口观看| 757午夜福利合集在线观看| 欧美日韩乱码在线| 91九色精品人成在线观看| 女性生殖器流出的白浆| 亚洲色图综合在线观看| 人人妻人人爽人人添夜夜欢视频| www.熟女人妻精品国产| 人妻 亚洲 视频| 国产精品久久电影中文字幕 | 母亲3免费完整高清在线观看| 国产精品乱码一区二三区的特点 | 91九色精品人成在线观看| 亚洲久久久国产精品| 国产精品欧美亚洲77777| 国产真人三级小视频在线观看| 首页视频小说图片口味搜索| 99久久人妻综合| 国产成人一区二区三区免费视频网站| 身体一侧抽搐| 久久久国产成人免费| 动漫黄色视频在线观看| 桃红色精品国产亚洲av| 精品国产亚洲在线| 两性夫妻黄色片| 超碰97精品在线观看| 成年人黄色毛片网站| 成人特级黄色片久久久久久久| 亚洲一区二区三区欧美精品| 又紧又爽又黄一区二区| 亚洲五月天丁香| 两性午夜刺激爽爽歪歪视频在线观看 | 黄片播放在线免费| 亚洲成人国产一区在线观看| 国产高清国产精品国产三级| 久久精品国产清高在天天线| 亚洲欧洲精品一区二区精品久久久| 国产精品国产av在线观看| 夫妻午夜视频| 黄色女人牲交| 高清在线国产一区| 午夜福利一区二区在线看| www.熟女人妻精品国产| 他把我摸到了高潮在线观看| svipshipincom国产片| av视频免费观看在线观看| 伦理电影免费视频| 亚洲精品一二三| 黄色怎么调成土黄色| 岛国毛片在线播放| 热re99久久国产66热| 老司机影院毛片| 无限看片的www在线观看| 亚洲七黄色美女视频| 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 老熟妇乱子伦视频在线观看| 大型av网站在线播放| 国产精品99久久99久久久不卡| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 大香蕉久久成人网| 免费在线观看影片大全网站| bbb黄色大片| 男女午夜视频在线观看| 狂野欧美激情性xxxx| 国产淫语在线视频| av一本久久久久| 天堂动漫精品| 黄色丝袜av网址大全| 高清视频免费观看一区二区| 丁香六月欧美| 俄罗斯特黄特色一大片| 成年人午夜在线观看视频| 在线观看午夜福利视频| 热re99久久精品国产66热6| 18禁黄网站禁片午夜丰满| 国产在线一区二区三区精| 国产精品亚洲av一区麻豆| 日韩 欧美 亚洲 中文字幕| 国产精品永久免费网站| 亚洲av美国av| 少妇猛男粗大的猛烈进出视频| 99精国产麻豆久久婷婷| 欧美国产精品va在线观看不卡| 亚洲av日韩在线播放| 亚洲精品一卡2卡三卡4卡5卡| 91国产中文字幕| 男人舔女人的私密视频| 99久久国产精品久久久| 国产一区有黄有色的免费视频| 欧洲精品卡2卡3卡4卡5卡区| 黄片小视频在线播放| 国产精品电影一区二区三区 | 丝袜美足系列| 亚洲精品中文字幕一二三四区| 国产精品免费视频内射| 高清毛片免费观看视频网站 | 午夜91福利影院| 久久久久久久精品吃奶| 国产亚洲精品久久久久久毛片 | 亚洲精品中文字幕在线视频| 国产精品美女特级片免费视频播放器 | 999久久久精品免费观看国产| 国产av又大| 18禁观看日本| 国产成人欧美在线观看 | 亚洲avbb在线观看|