周 浩，金 平，夏衛(wèi)生，雷國(guó)平

三江平原撓力河流域主要作物水分盈虧時(shí)空變化特征

周 浩1，金 平2，夏衛(wèi)生1，雷國(guó)平3※

（1. 湖南師范大學(xué)資源與環(huán)境科學(xué)學(xué)院，長(zhǎng)沙 410081；2. 湖北省城市規(guī)劃設(shè)計(jì)研究院，武漢 430071；3. 東北大學(xué)土地管理研究所，沈陽 110004）

明確中國(guó)糧食生產(chǎn)基地主要作物種植條件下的水分盈虧條件，有利于區(qū)域灌溉策略制定及保障國(guó)家糧食安全。該研究基于MODIS遙感數(shù)據(jù)和常規(guī)氣象數(shù)據(jù)，運(yùn)用Priestley-Taylor公式以及作物水分盈虧指數(shù)，以三江平原腹地的撓力河流域?yàn)檠芯繀^(qū)，揭示2000年以來該流域主要作物（中稻、春小麥和春玉米）的水分盈虧時(shí)空變化特征。結(jié)果表明：1）撓力河流域潛在蒸散量由2000年的910.25 mm增至2015年的964.04 mm，各作物需水量整體呈現(xiàn)不同幅度的上升態(tài)勢(shì)；2）各作物表現(xiàn)出不同的水分盈虧特點(diǎn)。中稻水分虧缺明顯，缺水量由東北向西南遞增，天然降水不能滿足其水分需求。春小麥水分虧缺程度低，灌溉需求量在空間上變化較小。春玉米多處于輕度缺水狀態(tài)，且西北部的缺水量小于東南部；3）流域境內(nèi)的中稻多處于輕旱狀態(tài)，2000—2015年干旱區(qū)持續(xù)擴(kuò)張，面積增長(zhǎng)幅度達(dá)到143.46%。對(duì)于旱作物而言，春小麥的缺水量尚未達(dá)到輕旱的標(biāo)準(zhǔn)，而春玉米的輕旱區(qū)面積往復(fù)變化，2000—2015年間面積減少了79.90%。作物種植結(jié)構(gòu)的快速變化，使得中稻種植區(qū)迅速擴(kuò)張，加劇了該流域的水分虧缺態(tài)勢(shì)，而春小麥和春玉米的種植區(qū)變化未對(duì)水分盈虧態(tài)勢(shì)造成顯著影響。研究可為撓力河流域的農(nóng)業(yè)結(jié)構(gòu)調(diào)整與農(nóng)田灌溉措施制定提供支持。

遙感；作物；Priestley-Taylor公式；水分盈虧；撓力河流域

0 引 言

水分是作物生長(zhǎng)發(fā)育的重要因子，而農(nóng)業(yè)旱澇的根本原因在于水分盈虧失衡[1]。對(duì)于北方干旱區(qū)而言，水分盈虧及其變化研究對(duì)于糧食生產(chǎn)與作物布局具有較強(qiáng)的現(xiàn)實(shí)意義[2]。因此，揭示不同作物種植條件下的水分盈虧條件，有利于科學(xué)引導(dǎo)農(nóng)業(yè)結(jié)構(gòu)調(diào)整及制定農(nóng)田排灌措施。

當(dāng)前，學(xué)者圍繞水分盈虧問題開展了大量研究。研究方法和手段上，?；谔镩g試驗(yàn)或傳統(tǒng)數(shù)理特征統(tǒng)計(jì)法[3-5]，來研究作物水分盈虧的分項(xiàng)特點(diǎn)，進(jìn)而分析區(qū)域水分盈虧特征。具體的評(píng)價(jià)指標(biāo)包括標(biāo)準(zhǔn)化降雨蒸散指數(shù)（Standardized Precipitation Evaporation Index，SPEI）[6]、帕爾默綜合氣象干旱指數(shù)（Palmer Drought Severity Index，PDSI）[7]、干濕指數(shù)[8]和作物水分盈虧指數(shù)（Crop Water Surplus & Deficit Index，CWSDI）[9]等，其中CWSDI指數(shù)能夠很好地反映作物水分供需矛盾大小以及作物對(duì)水分的需求程度；研究?jī)?nèi)容上，著重分析了氣象、作物和農(nóng)田土壤的水分盈虧特征，如高歌等[10]研究了中國(guó)主要流域的水分盈虧特征及其影響因素；張淑杰等[11-12]分析了東北地區(qū)氣象水分盈虧特征，并得出該地區(qū)水分盈虧量逐年下降的結(jié)論；姚曉軍等[13]和史建國(guó)等[14]分別對(duì)中國(guó)西北地區(qū)以及黃河流域的水分盈虧特征進(jìn)行了研究。上述研究對(duì)象多為較大面積的流域或地區(qū)，對(duì)于偏僻的中小尺度區(qū)域而言，由于境內(nèi)氣象水文要素?cái)?shù)據(jù)較為匱乏，不利于水分盈虧評(píng)估工作的開展。

位于中國(guó)東北端的三江平原具有特有的濕地生態(tài)服務(wù)功能和糧食生產(chǎn)地位，一直是全球變化研究的重點(diǎn)地區(qū)[15-16]。自21世紀(jì)以來，在“兩江一湖”土地整治以及農(nóng)戶自發(fā)型農(nóng)業(yè)結(jié)構(gòu)調(diào)整背景下，該地區(qū)作物種植結(jié)構(gòu)持續(xù)變化[17]，對(duì)作物水分盈虧造成了顯著影響。但三江平原地處邊境地帶，境內(nèi)氣象站點(diǎn)偏少，氣象資料匱乏。因此，科學(xué)評(píng)估該資料欠缺地區(qū)作物種植條件下的水分盈虧特征，對(duì)于區(qū)域糧食生產(chǎn)安全尤為重要。近年來，隨著計(jì)算機(jī)技術(shù)的發(fā)展，水文遙感正逐漸成為資料欠缺地區(qū)獲取水分條件數(shù)據(jù)的重要手段，能夠彌補(bǔ)傳統(tǒng)監(jiān)測(cè)資料不足的劣勢(shì)[18-20]?；诖耍疚倪x取三江平原腹地的撓力河流域?yàn)檠芯繀^(qū)，采用區(qū)域遙感反演手段，來研究資料欠缺地區(qū)的主要作物（中稻、春小麥和春玉米）水分盈虧特征，明確作物種植條件與氣候濕潤(rùn)狀況的相互匹配程度，為該流域乃至三江平原的農(nóng)業(yè)結(jié)構(gòu)管理與灌溉策略制定提供依據(jù)。

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

1.1 研究區(qū)概況

撓力河流域?yàn)槿皆硟?nèi)最大的流域，總面積為2.49萬km2（圖1）。該流域地處半干旱地帶，多年平均降水量518 mm，且主要集中于6—9月份，水系自西南流向東北，地勢(shì)西南高而東北低，地貌類型是山地和平原，二者分別占流域面積的38.3%和61.7%。自建國(guó)初以來，撓力河流域歷經(jīng)多次大規(guī)模土地利用開發(fā)，導(dǎo)致濕地大面積喪失、結(jié)構(gòu)破壞和功能退化，耕地面積持續(xù)上升，并已成為三江平原主要的糧食產(chǎn)區(qū)和國(guó)家重要商品糧生產(chǎn)基地，主要作物類型為中稻、春小麥和春玉米[21]。

圖1 撓力河流域地理位置

1.2 數(shù)據(jù)來源及預(yù)處理

本文數(shù)據(jù)包括：1）MCD43B3（1 d）、MOD11A2（1 d）、MCD15A2（8 d）和LST，以上數(shù)據(jù)均來自于MODIS陸地標(biāo)準(zhǔn)3級(jí)產(chǎn)品，對(duì)數(shù)據(jù)進(jìn)行拼接、裁剪、空間范圍缺值插補(bǔ)和時(shí)間序列缺值插補(bǔ)等一系列預(yù)處理后[18]，最終得到2000、2005、2010和2015年撓力河流域的MODIS數(shù)據(jù)集，其時(shí)間分辨率和空間分辨率分別為1 d和1 km；2）逐日降水和溫度數(shù)據(jù)均來自于中國(guó)氣象數(shù)據(jù)中心“中國(guó)地面氣象資料日值數(shù)據(jù)集V3.0”的寶清站、鶴崗站、佳木斯站、勃利站、雞西站、綏芬河站、富錦站和虎林站共8個(gè)地面基準(zhǔn)氣象站點(diǎn)數(shù)據(jù)；3）土地利用分布數(shù)據(jù)通過人機(jī)交互式目視解譯方法獲得，解譯基礎(chǔ)數(shù)據(jù)源為美國(guó)陸地資源衛(wèi)星Landsat TM/OLI遙感影像，并輔以中國(guó)科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心的“中國(guó)土地利用現(xiàn)狀遙感監(jiān)測(cè)數(shù)據(jù)”；4）基礎(chǔ)地理信息數(shù)據(jù)包括SRTM的DEM數(shù)字高程模型數(shù)據(jù)（90 m）和經(jīng)、緯度空間分布數(shù)據(jù)，對(duì)基礎(chǔ)數(shù)據(jù)經(jīng)投影變化（Albers等積投影）和重采樣處理操作（1 km）；5）主要作物生長(zhǎng)周期及面積結(jié)構(gòu)系數(shù)來自于國(guó)家氣象數(shù)據(jù)中心“中國(guó)農(nóng)作物生長(zhǎng)發(fā)育和農(nóng)田土壤濕度旬值數(shù)據(jù)集”和流域境內(nèi)分縣農(nóng)業(yè)統(tǒng)計(jì)資料。

2 研究方法

2.1 潛在蒸散量

2.1.1 瞬時(shí)潛在蒸散量

蒸散量是水循環(huán)中的重要支出項(xiàng)，準(zhǔn)確估算蒸散量是流域水循環(huán)與水分盈虧研究的基礎(chǔ)。選用Priestley-Taylor[22]公式作為瞬時(shí)潛在蒸散遙感反演的基礎(chǔ)，公式如下

式中ET0為潛在蒸散量，mm；為Priestley-Taylor系數(shù)；R為地表凈輻射量，W/m2，可通過MODIS數(shù)據(jù)來推算衛(wèi)星過境時(shí)刻，并結(jié)合大氣溫度T和Emis地表發(fā)射率數(shù)據(jù)進(jìn)行計(jì)算[20]；為土壤熱通量，W/m2；為汽化潛熱，MJ/kg，由空氣定壓比熱和大氣壓（經(jīng)海拔高度進(jìn)行推算）等計(jì)算得到[23]；為飽和水汽壓-溫度曲線斜率，kPa/℃，可基于大氣溫度T的經(jīng)驗(yàn)公式進(jìn)行計(jì)算；為干濕表常數(shù)，kPa/℃，可由空氣定壓比熱，大氣壓（海拔高度推算）等計(jì)算得到。

2.1.2 瞬時(shí)潛在蒸散量的逐日尺度轉(zhuǎn)換

一天中凈輻射R呈正弦曲線變化特點(diǎn)，決定了瞬時(shí)潛在蒸散量也具有正弦曲線變化特點(diǎn)，因此可以采用正弦曲線擬合的方法實(shí)現(xiàn)瞬時(shí)潛在蒸散量的逐日尺度轉(zhuǎn)換，公式如下

式中ET0_max為日最大潛在蒸散量，mm；為模擬時(shí)刻，h；rise為日出時(shí)刻，h；set為日落時(shí)刻，h。rise和set分別對(duì)應(yīng)著一天中凈輻射值變?yōu)檎拓?fù)的時(shí)刻。

2.2 有效降水量

作物生育階段有效降水量是指作物生育階段內(nèi)降水實(shí)際補(bǔ)充到根層土壤中的凈水量，代表總降水量中的有效部分。本文利用美國(guó)農(nóng)業(yè)部土壤保持局推薦的有效降水量計(jì)算方法[24]。

式中e為逐日有效降水量，mm/d；為逐日總降水量，mm/d。

2.3 作物需水量

作物需水量為確定作物灌溉量的基礎(chǔ)，受土壤、作物、氣候等多種因素影響，主要研究方法包括田間測(cè)定和理論計(jì)算兩類，本文采用作物系數(shù)法進(jìn)行計(jì)算。公式如下[1]

式中ET為作物需水量，mm；ET0采用上文逐日潛在蒸散量的計(jì)算結(jié)果；K為作物系數(shù)，反映了作物在標(biāo)準(zhǔn)狀態(tài)下的蒸散特性。該系數(shù)受土壤、氣候、作物生長(zhǎng)狀況和管理方式等多種因素綜合影響，對(duì)于K的確定，可采用FAO推薦的標(biāo)準(zhǔn)作物系數(shù)與修正公式進(jìn)行作物系數(shù)修正[9]。公式如下

式中K為作物不同生育階段的標(biāo)準(zhǔn)作物系數(shù)；2為2 m高處的風(fēng)速，m/s；RHmin為日最低相對(duì)濕度平均值，%；為作物平均高度值，m；同時(shí)依據(jù)中稻、春小麥和春玉米的生長(zhǎng)發(fā)育特點(diǎn)，將全生育期分為為4個(gè)生育階段（初期、發(fā)育期、中期和成熟期），依據(jù)撓力河流域的作物實(shí)際生長(zhǎng)發(fā)育資料，劃分對(duì)應(yīng)的作物生育周期及確定作物系數(shù)（表1），進(jìn)而確定全生育期作物的需水量。

表1 撓力河流域主要作物生育階段劃分及其作物系數(shù)

注：K為作物系數(shù)。

Note:Kis crop coefficient.

2.4 作物水分盈虧評(píng)價(jià)

考慮到作物不同生育階段對(duì)應(yīng)著不同需水量，構(gòu)建作物水分盈虧量與作物需水量的比值CWSDI，以反映其水分盈虧程度，公式如下

式中=P-ET表示作物實(shí)際水分盈虧量的大小，mm。當(dāng)CWSDI＞0，表示水分盈余，CWSDI=0，表示水分收支平衡，CWSDI＜0，表示水分虧缺[9]。

3 結(jié)果與分析

3.1 潛在蒸散量特征

潛在蒸散量遙感反演結(jié)果顯示，撓力河流域4個(gè)研究時(shí)點(diǎn)（2000、2005、2010和2015年）的年潛在蒸散量分別為910.25、936.84、937.23和964.04 mm。在空間分布上，撓力河流域各研究時(shí)點(diǎn)的平均潛在蒸散量處于780～1 179 mm范圍，其中農(nóng)作物主要分布區(qū)的潛在蒸散值普遍高于940 mm，流域東部及南部多山，海拔較高，潛在蒸散量整體偏低。相對(duì)于流域中部的內(nèi)、外七星河地區(qū)而言，上游寶清縣地區(qū)的相對(duì)濕度低，加之受緯度帶影響，導(dǎo)致潛在蒸散量偏高（圖2）。

圖2 撓力河流域平均潛在蒸散量空間分布

3.2 有效降水量特征

撓力河流域4個(gè)研究時(shí)點(diǎn)的平均中稻全生育期有效降水量略低于春小麥和春玉米的有效降水量。中稻生長(zhǎng)期約為130 d（5月20日至9月29日），生育期內(nèi)平均有效降水量為291.81 mm，占全年總降水量的60.16%，降水較為集中。春小麥生長(zhǎng)期約為115 d（5月8日至9月5日），平均有效降水量高于中稻，為303.31 mm。春玉米的生長(zhǎng)期共135 d（4月26日至9月12日），平均有效降水量達(dá)到324.77 mm（圖3）。

3.3 作物需水量特征

對(duì)于中稻而言，4個(gè)研究時(shí)點(diǎn)（2000、2005、2010和2015年）的全生育期需水量逐步上升，依次為609.63、617.02、673.23和673.65 mm。2000年和2005年，流域北端的富錦市與饒河縣交界區(qū)的需水量偏低，撓力河干流沿岸的富錦寶清段是需水高值集中區(qū)，平均值達(dá)到620 mm。至2010年，中稻需水量空間分布發(fā)生了較大變化，呈現(xiàn)出由西向東不斷遞減的特點(diǎn)，相對(duì)高值區(qū)多分布于寶清縣境內(nèi)，相對(duì)低值區(qū)多分布于北端的富錦市境內(nèi)，但面積較小。2015年，流域境內(nèi)的需水高值區(qū)分布較為零散，但由于該年份的寶清縣溫度偏低，使得區(qū)域潛在蒸散量偏低，并導(dǎo)致對(duì)應(yīng)的中稻需水量也處于較低水平。

作為旱作物的春小麥和春玉米，其需水量顯著低于中稻。4個(gè)研究時(shí)點(diǎn)的春小麥需水量均值依次為460.31、464.22、468.17和461.34 mm，且相對(duì)空間分布差異較大，2000年和2005年的相對(duì)高值區(qū)多位于內(nèi)、外七星河地區(qū)，撓力河干流沿岸的饒河寶清段和流域南端部分地區(qū)的需水量偏低，2010年，春小麥需水量在空間上整體表現(xiàn)出西高東低的特點(diǎn)，而至2015年，流域南部的需水量高于東部，北部富錦市與友誼縣交界區(qū)的需水量處于較低水平；與春小麥相比，春玉米全生育期較長(zhǎng)，其需水量高于春小麥，4個(gè)研究時(shí)點(diǎn)均值依次為510.63、502.24、522.73和518.68 mm。2000年，需水高值和低值區(qū)交錯(cuò)分布，趨勢(shì)性特征差，2005年，春玉米需水量下降幅度較小，且在空間上仍保持著高、低值混合分布的特點(diǎn)，2010年和2015年，春玉米需水量東高西低?？傮w而言，位于西部的友誼集賢地區(qū)和東部寶清縣地區(qū)的旱作物需水波動(dòng)較大（圖4）。

3.4 作物水分盈虧特征

依據(jù)CWSDI作物水分盈虧模型，評(píng)估撓力河流域主要作物的水分盈虧特征（圖5），并依據(jù)《中國(guó)農(nóng)業(yè)干旱監(jiān)測(cè)指標(biāo)方法與等級(jí)劃分》標(biāo)準(zhǔn)（無旱（?0.50

3.4.1 作物水分盈虧指數(shù)特征

各研究時(shí)點(diǎn)的中稻水分盈虧指數(shù)均小于0，反映了中稻全生育期以水分虧缺狀態(tài)為主，天然降水難以滿足該地區(qū)的灌溉需求。不同時(shí)間點(diǎn)的水分盈虧指數(shù)存在差別，對(duì)應(yīng)著中稻不同的灌溉需求度。2000年，中稻水分盈虧指數(shù)處于?0.50～?0.72范圍，平均值為?0.523，整體處于輕度干旱狀態(tài)，人工灌溉需求量由東北向西南遞增；相對(duì)于2000年，2005年的水分虧缺態(tài)勢(shì)加劇，盈虧指數(shù)均值變?yōu)?0.532，對(duì)應(yīng)的指數(shù)范圍為?0.32～?0.62；2010年水分盈虧指數(shù)處于?0.30～?0.65范圍，平均值為?0.57，盈虧指數(shù)下降幅度較大，流域水分虧缺態(tài)勢(shì)趨于惡化，而至2015年，水分盈虧指數(shù)的平均值變?yōu)?0.58。同時(shí)在空間分布上，水分盈虧指數(shù)絕對(duì)值均由東北向西南遞增，友誼集賢地區(qū)以及富錦市西南端的中稻灌溉壓力較大。

圖4 撓力河流域典型年份的作物需水量空間分布

對(duì)于春小麥而言，各研究時(shí)點(diǎn)的水分盈虧指數(shù)均值小于0，但灌溉需求量小于中稻。2000年，春小麥的水分盈虧指數(shù)處于?0.38～?0.53范圍，平均值為?0.44。依據(jù)干旱等級(jí)劃分標(biāo)準(zhǔn)，流域大部分地區(qū)的春小麥處于無旱狀態(tài)，即無需額外的人工灌溉措施來滿足春小麥的水分灌溉需求，可利用天然降水、土壤水分條件進(jìn)行正常生長(zhǎng)發(fā)育；2005年，春小麥的水分盈虧指數(shù)范圍變?yōu)?0.29～?0.45，對(duì)應(yīng)的均值為?0.359，水分虧缺量下降；2010年和2015年的水分盈虧指數(shù)均值分別為?0.362和?0.365，水分需求量增加，且在空間上表現(xiàn)出需水量差異較小的特點(diǎn)。

各研究時(shí)點(diǎn)的春玉米均為輕度缺水狀態(tài)，且缺水程度處于相似水平，西北部的友誼和集賢地區(qū)以及富錦市西南端的春玉米水分虧缺量顯著高于東南部的寶清縣地區(qū)。2000年，春玉米的水分盈虧指數(shù)范圍為?0.40～?0.61，均值達(dá)到?0.537，處于輕度干旱狀態(tài)，天然降水難以滿足其灌溉需求；相對(duì)于2000年，2005年的春玉米盈虧指數(shù)均值為?0.47，水分虧缺量下降，對(duì)應(yīng)的指數(shù)范圍為?0.42～?0.54；2010年和2015年，水分盈虧指數(shù)均值分別為?0.489和?0.484，水分灌溉需求度有所增加。

注：CWSDI為作物水分盈虧指數(shù)。

Note: CWSDI: Crop Water Surplus & Deficit Index.

圖5 撓力河流域典型年份的作物水分盈虧指數(shù)空間分布特征

Fig.5 Spatial distribution of crop water surplus & deficit index in typical years in Naoli River Basin

3.4.2 作物干旱特征

依據(jù)水分盈虧指數(shù)和干旱等級(jí)劃分標(biāo)準(zhǔn)，對(duì)撓力河流域主要作物的干旱程度進(jìn)行評(píng)價(jià)，發(fā)現(xiàn)撓力河流域中稻種植區(qū)多處于輕度干旱狀態(tài)，需通過抽取地下水、渠道引水等措施，以保障該地區(qū)水稻的正常生長(zhǎng)發(fā)育。2000年，約有6.72%面積的中稻處于中度干旱狀態(tài)，且多位于富錦市東南部；2005、2010和2015年，春小麥水分盈虧態(tài)勢(shì)良好，盡管處于水分虧缺的狀態(tài)，但可通過土壤儲(chǔ)水以及適當(dāng)?shù)霓r(nóng)田灌溉措施，以保證春小麥的正常生長(zhǎng)發(fā)育；春玉米和春小麥同為旱作物，且春玉米生長(zhǎng)期的需水量顯著高于春小麥，2000—2015年，該流域春玉米表現(xiàn)出不同程度的輕度干旱特點(diǎn)，2000年的輕旱區(qū)相對(duì)面積比例達(dá)到94.43%，其余年份的面積占比依次為17.65%、31.24%和24.08%（表2）。

表2 撓力河流域作物不同干旱等級(jí)面積占比

自21世紀(jì)初以來，撓力河流域中稻種植區(qū)迅速擴(kuò)張，加劇了該流域農(nóng)田的水分虧缺情勢(shì)，而春小麥和春玉米種植結(jié)構(gòu)的變化，未對(duì)該流域的水分盈虧情勢(shì)造成顯著的負(fù)面影響。具體表現(xiàn)為：輕旱是撓力河流域中稻的主要干旱類型，2000年以來，在中稻種植面積迅速增加的背景下，輕旱區(qū)面積約增加了3 225.00 km2，16 a間擴(kuò)張幅度達(dá)到143.46%，每5 a相鄰時(shí)點(diǎn)的擴(kuò)張幅度也分別為25.71%、39.21%和39.12%；2000年，春小麥的輕旱區(qū)面積為803.00 km2，相對(duì)面積比例達(dá)到6.74%，其余年份的春小麥均處于無旱狀態(tài)；對(duì)于春玉米而言，研究期內(nèi)的輕旱區(qū)面積共減少了8 931.00 km2，變化幅度達(dá)到79.90%。

4 討 論

合理的水分條件是農(nóng)業(yè)生產(chǎn)的先決條件[16]，從柵格尺度評(píng)估區(qū)域水分盈虧特征，能夠指導(dǎo)農(nóng)田水分的科學(xué)管理。水分盈虧研究一般被分為涉及降水、蒸發(fā)等水文過程的天然狀態(tài)和額外增加人類活動(dòng)影響的人工狀態(tài)2個(gè)層次[25]。前者側(cè)重于植被水分需求量與自然供水量之間的差額問題，可從氣候、作物與農(nóng)田3個(gè)層次進(jìn)行研究[26]，被認(rèn)為是在水分管理過程中逐層增加限制性因素，其中氣候水分盈虧與作物水分盈虧是天然狀態(tài)下水分盈虧研究的初始層次，更具理論意義，農(nóng)田土壤水分盈虧是天然狀態(tài)下水分盈虧研究的實(shí)際層次，更具實(shí)踐價(jià)值[27]。本文以三江平原腹地的撓力河流域?yàn)檠芯繀^(qū)，開展作物種植條件下的水分盈虧研究，為理論意義上的水分平衡關(guān)系，后續(xù)可考慮研究農(nóng)田土壤水分盈虧關(guān)系以提高評(píng)估結(jié)果的實(shí)踐價(jià)值，或依次從作物、農(nóng)田和土壤層次來研究水分盈虧問題。

以水稻為主要種植類型的水田是撓力河流域重要的糧食生產(chǎn)載體，灌溉需水量偏大。在水田面積持續(xù)擴(kuò)張下，水田的水分虧缺態(tài)勢(shì)愈加突出，應(yīng)將水田或水稻視作水分條件管理的重點(diǎn)對(duì)象。宮興龍等[28-29]也對(duì)撓力河流域的水土資源綜合利用問題開展了系統(tǒng)性研究，并認(rèn)為水稻種植對(duì)區(qū)域干旱特征影響非常顯著，研究結(jié)論與本文基本一致。但由于作物尤其是水稻不同生長(zhǎng)階段的需水量不同，使得對(duì)應(yīng)的水分盈虧量也存在差別，后續(xù)研究應(yīng)細(xì)化作物生育階段，來揭示不同生長(zhǎng)階段的水分盈虧特點(diǎn)，并針對(duì)性地指導(dǎo)該流域農(nóng)業(yè)灌溉管理。

5 結(jié) 論

基于MODIS遙感影像和常規(guī)氣象數(shù)據(jù)，在Priestley-Taylor公式以及作物水分盈虧模型支持下，揭示2000年以來撓力河流域主要作物（中稻、春小麥和春玉米）的水分盈虧時(shí)空變化特征。主要得出以下結(jié)論：

1）撓力河流域潛在蒸散量由2000年的910.25 mm增至2015年的964.04 mm，中稻、春小麥和春玉米的全生育期有效降水量依次為291.81、303.31和324.77 mm；

2）撓力河流域的中稻以水分虧缺特點(diǎn)為主，缺水量由東北向西南遞增，天然降水難以滿足該流域水田的灌溉需求；春小麥的水分盈虧態(tài)勢(shì)良好，各年份間的空間需水量處于同一水平；春玉米的水分虧缺態(tài)勢(shì)較為嚴(yán)峻，且虧缺量由東南部向西北部緩慢增加；

3）撓力河流域中稻多處于輕度干旱狀態(tài)，2000—2015年面積持續(xù)增長(zhǎng)了143.46%。流域的春小麥多未達(dá)到缺水的干旱標(biāo)準(zhǔn)，而春玉米的輕旱區(qū)面積減少了79.90%。中稻面積的迅速擴(kuò)張，加劇了該流域的水分虧缺態(tài)勢(shì)，春小麥和春玉米種植結(jié)構(gòu)的改變，未對(duì)水分盈虧造成顯著的負(fù)面影響。

[1] 高超，李學(xué)文，孫艷偉，等. 淮河流域夏玉米生育階段需水量及農(nóng)業(yè)干旱時(shí)空特征[J]. 作物學(xué)報(bào)，2019，45(2)：297-309. Gao Chao, Li Xuewen, Sun Yanwei, et al. Spatio temporal characteristics of water requirement and agricultural drought during summer maize season in Huaihe River Basin[J]. Acta Agronomica Sinica, 2019, 45(2): 297-309. (in Chinese with English abstract)

[2] 賈艷青，張勃. 1960—2016年中國(guó)北方地區(qū)極端干濕事件演變特征[J]. 自然資源學(xué)報(bào)，2019，34(7)：1543-1554. Jia Yanqing, Zhangbo. Spatio-temporal changes of the extreme drought and wet events in Northern China from 1960 to 2016[J]. Journal of Natural Resources, 2019, 34(7): 1543-1554. (in Chinese with English abstract)

[3] 陳東東，王曉東，王森，等. 四川省潛在蒸散量變化及其氣候影響因素分析[J]. 中國(guó)農(nóng)業(yè)氣象，2017，38(9)：548-557. Chen Dongdong, Wang Xiaodong, Wang Sen, et al. Potential evapotranspiration changes and its effects of meteorological factors across Sichuan Province[J]. Chinese Journal of Agrometeorology, 2017, 38(9): 548-557. (in Chinese with English abstract)

[4] 南紀(jì)琴，王景雷，秦安振，等. 中國(guó)西北旱區(qū)農(nóng)業(yè)水土資源利用情景潛力研究[J]. 自然資源學(xué)報(bào)，2017，32(2)：292-300. Nan Jiqin, Wang Jinglei, Qin Anzhen, et al. Study on utilization potential of agricultural soil and water resources’ in Northwest Arid Area[J]. Journal of Natural Resources, 2017, 32(2): 292-300. (in Chinese with English abstract)

[5] 趙俊芳，郭建平，徐精文，等. 基于濕潤(rùn)指數(shù)的中國(guó)干濕狀況變化趨勢(shì)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(8)：18-24. Zhao Junfang, Guo Jianping, Xu Jingwen, et al. Trends of Chinese dry-wet condition based on wetness index[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(8): 18-24. (in Chinese with English abstract)

[6] Vicente-Serrano S M, Beguería S, López-Morenoj J L. A multi-scalar drought index sensitive to global warming: The standardized precipitation evapotranspiration index[J]. Journal of Climate, 2010, 23(7): 1696-1718.

[7] Rayne S, Forest K. Evidence for increasingly variable Palmer Drought Severity Index in the United States since 1895[J]. Science of the Total Environment, 2016, 544: 792-796.

[8] 任志遠(yuǎn)，李冬玉，劉憲鋒. 1960—2011年秦嶺南北地表干濕變化時(shí)空特征[J]. 陜西師范大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014，42(1)：81-85. Ren Zhiyuan, Li Dongyu, Liu Xianfeng. Spatial and temporal characteristics of surface dry-wet variations in Qinling Mountain from 1960 to 2011[J]. Journal of Shaanxi Normal University: Natural Science Edition, 2014, 42(1): 81-85. (in Chinese with English abstract)

[9] 魏新光，王鐵良，李波，等. 遼寧省玉米地水分盈虧時(shí)空分布特征及灌溉模式分區(qū)研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2018，34(23)：119-126. Wei Xinguang, Wang Tieliang, Li Bo, et al. Temporal and spatial distribution characteristics of maize water surplus deficit and irrigation mode partition in Liaoning province[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(23): 119-126. (in Chinese with English abstract)

[10] 高歌，許崇育. 1961—2010年中國(guó)十大流域水分盈虧量時(shí)空變化特征[J]. 地理學(xué)報(bào)，2015，70(3)：380-391. Gao Ge, Xu Chongyu. Characteristics of water surplus and deficit change in 10 major river basins in China during 1961-2010[J]. Acta Geographica Sinica, 2015, 70(3): 380-391. (in Chinese with English abstract)

[11] 張淑杰，張玉書，蔡福，等. 1961年至2007年東北地區(qū)水分盈虧變化趨勢(shì)及其影響與響應(yīng)[J]. 資源科學(xué)，2010，32(4)：663-670. Zhang Shujie, Zhang Yushu, Cai Fu, et al. Changes of water budget and their effects and responses in northeastern China from 1961 to 2007[J]. Resources Science, 2010, 32(4): 663-670. (in Chinese with English abstract)

[12] 曾麗紅，宋開山，張柏. 1951—2008 年東北地區(qū)水分盈虧量時(shí)空格局[J]. 自然資源學(xué)報(bào)，2011，26(5)：858-870. Zeng Lihong, Song Kaishan, Zhang Bai. Spatial and temporal structure of water deficit over northeast China during the period of 1951-2008[J]. Journal of Natural Resources, 2011, 26(5): 858-870. (in Chinese with English abstract)

[13] 姚曉軍，張曉，孫美平，等. 1960—2010 年中國(guó)西北地區(qū)水分盈虧量時(shí)空特征[J]. 地理研究，2013，32(4)：607-616. Yao Xiaojun, Zhang Xiao, Sun Meiping, et al. Spatial-temporal characteristics of water deficit in Northwest China from 1960 to 2010[J]. Geographical Research, 2013, 32(4): 607-616. (in Chinese with English abstract)

[14] 史建國(guó)，嚴(yán)昌榮，何文清，等. 黃河流域水分虧缺時(shí)空格局變化研究[J]. 自然資源學(xué)報(bào)，2008，23(1)：113-119. Shi Jianguo, Yan Changrong, He Wenqing, et al. Study on spatial and temporal variation of water surplus and deficiency in Yellow River basin[J]. Journal of Natural Resources, 2008, 23(1): 113-119. (in Chinese with English abstract)

[15] 吳海燕，孫甜田，范作偉，等. 東北地區(qū)主要糧食作物對(duì)氣候變化的響應(yīng)及其產(chǎn)量效應(yīng)[J]. 農(nóng)業(yè)資源與環(huán)境學(xué)報(bào)，2014，31(4)：299-307. Wu Haiyan, Sun Tiantian, Fan Zuowei, et al. The major food crops in response to climate change and its yield effect in Northeast of China[J]. Journal of Agricultural Resources and Environment, 2014, 31(4): 299-307. (in Chinese with English abstract)

[16] 劉彥隨，甘紅，張富剛. 中國(guó)東北地區(qū)農(nóng)業(yè)水土資源匹配格局[J]. 地理學(xué)報(bào)，2006，61(8)：847-854. Liu Yansui, Gan Hong, Zhang Fugang. Analysis of the matching patterns of land and water resources in Northeast China[J]. Acta Geographica Sinica, 2006, 61(8): 847-854. (in Chinese with English abstract)

[17] 杜國(guó)明，李昀，于鳳榮，等. 基于遙感的2000—2009年三江平原北部耕地變化特征分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2012，28(1)：225-229. Du Guoming, Li Yun, Yu Fengrong et al. Change characteristics analysis of farmland in Northern Sanjiang Plain in 2000—2009 based on remote sensing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(1): 225-229. (in Chinese with English abstract)

[18] 楊勝天，趙長(zhǎng)森. 遙感水文[M]. 北京：科學(xué)出版社，2015.

[19] 徐宗學(xué)，趙捷. 生態(tài)水文模型開發(fā)和應(yīng)用：回顧與展望[J]. 水利學(xué)報(bào)，2016，47(3)：346-354. Xu Zongxue, Zhao Jie. Development and applications of eco-hydrological models: past and future[J]. Journal of Hydraulic Engineering, 2016, 47(3): 346-354. (in Chinese with English abstract)

[20] 劉昌明，白鵬，王中根，等. 稀缺資料流域水文計(jì)算若干研究：以青藏高原為例[J]. 水利學(xué)報(bào)，2016，47(3)：272-282. Liu Changming, Bai Peng, Wang Zhonggen, et al. Study on prediction of ungagged basins：A case study on the Tibetan Plateau[J]. Journal of Hydraulic Engineering, 2016, 47(3): 272-282. (in Chinese with English abstract)

[21] 周浩，雷國(guó)平，楊雪昕，等. RCPs氣候情景下?lián)狭恿饔蚋厮临Y源平衡效應(yīng)模擬[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2018，29(4)：1302-1312. Zhou Hao Lei Guoping, Yang Xuexin, et al. Simulation for balanced effect of soil and water resources on cultivated land in Naoli River Basin, Northeast China under the RCPs climate scene[J]. Chinese Journal of Applied Ecology, 2018, 29(4): 1302-1312. (in Chinese with English abstract)

[22] Priestley C H B, Taylor R J. On the assessment of surface heat flux and evaporation using large-scale parameters[J]. Monthly Weather Review, 1972, 100(2): 81-92.

[23] Liang S. Narrowband to broadband conversions of land surface albedo I Algorithms[J]. Remote Sensing of Environment, 2001, 76(2): 213-238.

[24] 羅那那，巴特爾·巴克，吳燕鋒. 基于標(biāo)準(zhǔn)化降水指數(shù)的北疆地區(qū)近52年旱澇變化特征[J]. 水土保持研究，2017，24(2)：293-299. Luo Nana, Batur Bake, Wu Yanfeng. Characteristics of drought and flood based on standardized precipitation index in Northern Xinjiang in recent 52 years [J]. Research of Soil and Water Conservation, 2017, 24(2): 293-299. (in Chinese with English abstract)

[25] 任志遠(yuǎn)，郭彩玲. 區(qū)域水土資源平衡與灌溉優(yōu)化模型研究：以陜西關(guān)中灌區(qū)為例[J]. 干旱區(qū)地理，2000(3)：264-268. Ren Zhiyuan, Guo Cailing. Regional water and soil resources balance & optimizing model research：taking Guanzhong Irrigation Region in Shanxi as an Example[J]. Arid Land Geography, 2000(3): 264-268. (in Chinese with English abstract)

[26] 賈紹鳳，封志明，李麗娟，等. 資源地理與水土資源研究成果與展望[J]. 自然資源學(xué)報(bào)，2010，25(9)：1445-1457. Jia Shaofeng, Feng Zhiming, LI Lijuan, et al. Research results and prospects of resource geography and water & land resources[J]. Journal of Natural Resources, 2010, 25(9): 1445-1457. (in Chinese with English abstract)

[27] 楊艷昭，張偉科，劉登偉. 內(nèi)蒙古水土資源平衡及其水資源承載能力[J]. 干旱區(qū)地理，2008(3)：436-441. Yang Yanzhao, Zhang Weike, Liu Dengwei. Ater and land balance and water resources carrying capacity in Inner Mongolia[J]. Arid Land Geography, 2008(3): 436-441. (in Chinese with English abstract)

[28] 宮興龍，付強(qiáng)，關(guān)英紅，等. 土地結(jié)構(gòu)改變對(duì)水土資源平衡的影響分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2018，49(4)：319-329. Gong Xinglong, Fu Qiang, Guan Yinghong, et al. Influence of land structure change on balance of soil and water resources[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(4): 319-329. (in Chinese with English abstract)

[29] 張瑩，雷國(guó)平，張弘強(qiáng)，等. 微觀尺度分析撓力河流域耕地利用水土資源匹配時(shí)空動(dòng)態(tài)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(8)：185-194. Zhang Ying, Lei Guoping, Zhang Hongqiang, et al. Spatiotemporal dynamics of land and water resources matching of cultivated land use based on micro scale in Naoli River Basin[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(8): 185-194. (in Chinese with English abstract)

Spatial-temporal characteristics of water surplus and deficit for the main crops in Naoli River Basin of Sanjiang Plain, China

Zhou Hao1, Jin Ping2, Xia Weisheng1, Lei Guoping3※

(1.,,410081,; 2.,430071,;3.,,110004,)

Naoli River Basin is located within the Sanjiang Plain in Heilongjiang Province, China. The agricultural structure of Sanjiang Plain has been greatly adjusted with a large amount of dry land crops into wetland rice since the early 2000s, thereby resulting in a series of problems in the utilization of water resources in Naoli River Basin. Therefore, it is necessary to explore the water profit and loss situation of main crops, in order to implement regional irrigation strategies, and thereby to serve as the national food security. In this study, the Priestley Taylor formula, and Crop Water Surplus and Deficit Index (CWSDI) were used to reveal the water profit and loss situation of the middle rice, spring wheat and spring maize, in Naoli River Basin. The fundamental resources included the MODIS image, the long-term sequenced meteorological and the DEM data sets during the processing. Results showed that the Potential Evapotranspiration (ET0) increased from 910.25 mm in 2000 to 964.04 mm in 2015 in this river basin. The water demand of main crops, including middle crop, spring wheat and spring maize, also increased at different ranges as a whole. There was a decrease trend in the CWSDI from 2000 to 2015 in Naoli River Basin, indicating the natural precipitations can not meet the irrigation demand of middle rice, while the absolute value of index increased from the northeast to southwest in the spatial distribution. The CWSDI of spring wheat was larger than that of the middle rice, showing similar characteristics in the whole spatial distribution for each year. The CWSDI of spring maize indicated that the spatial characteristics of northwest were slightly lower than that of the southeast, indicating both regions in the condition of mild water shortage in Naoli River Basin. Drought classification was also evaluated using the CWSDI. The middle rice, accounting for the area proportion of 143.46%, was basically in mild drought state in the research basin. Despite at the state of water shortage, spring wheat has not reached the drought standard. The cultivated area of spring maize has decreased by 79.90% in the mild drought state in Naoli River Basin. The reason can be that the spring wheat and spring maize have been partly substituted for the middle rice with a rapid increasing trend since 2000, and thereby the planting of wetland crops, such as rice, led to a severe shortage of water. The spring wheat and spring maize demonstrated an extremely rare negative impact on the water surplus and deficit situation. These findings can provide the promising references when implementing the farmland irrigation schemes and agricultural structure adjustment in Naoli River Basin, China.

remote sensing; crops; Priestley-Taylor formula; water surplus & deficit; Naoli River Basin.

周浩，金平，夏衛(wèi)生，等. 三江平原撓力河流域主要作物水分盈虧時(shí)空變化特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(14)：159-166.doi：10.11975/j.issn.1002-6819.2020.14.019 http://www.tcsae.org

Zhou Hao, Jin Ping, Xia Weisheng, et al. Spatial-temporal characteristics of water surplus and deficit for the main crops in Naoli River Basin of Sanjiang Plain, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 159-166. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.14.019 http://www.tcsae.org

2019-11-12

2020-06-28

國(guó)家自然科學(xué)基金項(xiàng)目（41671520）；國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目（2016YFC0502406）

周浩，博士，講師，主要從事耕地及國(guó)土資源綜合利用研究。Email：zhouhao7404@ 163.com

雷國(guó)平，博士，教授，主要從事土地利用規(guī)劃與管理研究。Email：guopinglei@126.com

10.11975/j.issn.1002-6819.2020.14.019

F323.2；F301.21

A

1002-6819(2020)-14-0159-08