水足跡視角下黃河流域經(jīng)濟(jì)作物水土資源匹配格局

焦士興，李玉靖，王安周，劉清春，尹義星，趙榮欽，林璐霜

（1.安陽(yáng)師范學(xué)院資源環(huán)境與旅游學(xué)院，河南安陽(yáng) 455002； 2.河南大學(xué)地理與環(huán)境學(xué)院，河南開(kāi)封 475004；3.洛陽(yáng)市第十九中學(xué)，河南洛陽(yáng) 471000； 4.山東財(cái)經(jīng)大學(xué)經(jīng)濟(jì)學(xué)院，山東濟(jì)南 250014；5.南京信息工程大學(xué)水文氣象學(xué)院，江蘇南京 210044；6.華北水利水電大學(xué)測(cè)繪與地理信息學(xué)院，河南鄭州 450046）

水土資源是農(nóng)業(yè)生產(chǎn)的核心資源，關(guān)系到糧食安全［1］。 水土資源的有效配置是農(nóng)業(yè)生產(chǎn)的前提，其合理配置有利于將資源優(yōu)勢(shì)轉(zhuǎn)變?yōu)榻?jīng)濟(jì)優(yōu)勢(shì)，促進(jìn)區(qū)域農(nóng)業(yè)持續(xù)發(fā)展［2］。 水足跡不僅可以定量反映經(jīng)濟(jì)作物生產(chǎn)中隱含的水資源量，而且可以表征經(jīng)濟(jì)作物貿(mào)易中的虛擬水資源量［3］，分析經(jīng)濟(jì)作物水足跡對(duì)于緩解農(nóng)業(yè)資源壓力、提高水土匹配度具有重要意義［4］。水土匹配度可以衡量某一區(qū)域內(nèi)水土資源的匹配狀況［5］，合理的水土匹配度可以為農(nóng)業(yè)發(fā)展提供有力支撐。 目前，主要采用基尼系數(shù)［6-7］、數(shù)據(jù)包絡(luò)分析［8-9］等方法開(kāi)展水土匹配度研究。 侯淑濤等［6］利用基尼系數(shù)，研究了黑龍江省農(nóng)業(yè)水土資源匹配協(xié)調(diào)水平及區(qū)域差異性；徐娜等［9］基于DEA 模型研究了甘肅省5個(gè)流域農(nóng)業(yè)水土資源配置效率變化趨勢(shì)及影響因素。已有研究［5-11］為水土資源匹配研究提供了理論支持和方法論基礎(chǔ)，但多以水資源總量、可利用水資源量或灌溉水量表示水資源量，而水足跡［12-13］是水資源占用的綜合表征，可以量化生產(chǎn)某種農(nóng)作物所消耗的真實(shí)用水量。

黃河流域是我國(guó)重要的工農(nóng)業(yè)生產(chǎn)核心區(qū)之一，水資源總量少，但農(nóng)用地面積大，用水需求已超出流域水資源承載能力，水資源供需矛盾突出［14］。 基于此，考慮不同農(nóng)作物需水的區(qū)域差異，采用水土匹配度、重心遷移模型分析1998—2020 年黃河流域8 個(gè)省（區(qū)）（四川省主要受長(zhǎng)江水系影響［15］，因此去除四川?。┎煌?jīng)濟(jì)作物播種面積及其水足跡，進(jìn)而探討水土資源時(shí)空匹配特征和重心軌跡演化情況。

1 數(shù)據(jù)來(lái)源和研究方法

黃河流域谷物、甜菜、茶葉、蔬菜等12 種經(jīng)濟(jì)作物生產(chǎn)總量、單產(chǎn)量以及播種面積數(shù)據(jù)來(lái)源于《中國(guó)統(tǒng)計(jì)年鑒》（1999—2021 年）。

1）經(jīng)濟(jì)作物水足跡。 水足跡是指人口在一定時(shí)間內(nèi)消費(fèi)的所有產(chǎn)品和服務(wù)所需要的水資源量。 作物生產(chǎn)水足跡包括藍(lán)水足跡和綠水足跡，分別指作物生產(chǎn)實(shí)際消耗的灌溉水量和有效降水量［16-17］。 目前，常采用作物蒸散量來(lái)計(jì)算有效降水量［18］，利用Penman-Monteith 公式計(jì)算蒸發(fā)蒸騰水量，并通過(guò)CROPWAT模型模擬得到藍(lán)水足跡和綠水足跡［19］。 借鑒王毅鑫等［15］2013 年研究黃河流域棉花、油料、甜菜等經(jīng)濟(jì)作物的虛擬水含量，結(jié)合不同經(jīng)濟(jì)作物產(chǎn)量，計(jì)算其水足跡。 計(jì)算公式為

式中：Wi為?。▍^(qū)）i經(jīng)濟(jì)作物水足跡，Di，c為省（區(qū)）i經(jīng)濟(jì)作物c單位質(zhì)量虛擬水含量，Bi，c為?。▍^(qū)）i經(jīng)濟(jì)作物c產(chǎn)量。

2）經(jīng)濟(jì)作物水土匹配度。 水土匹配度可以反映特定區(qū)域農(nóng)業(yè)生產(chǎn)可供水資源量與耕地資源在時(shí)空上適宜匹配的量化關(guān)系［7］，實(shí)質(zhì)上是單位面積耕地?fù)碛械乃阚E［19］。 水土匹配度越大，單位面積耕地消耗的水資源量越大，匹配性越低；水土匹配度越小，單位面積耕地消耗的水資源量越小，匹配性越高。 計(jì)算公式為

式中：Ri為?。▍^(qū)）i經(jīng)濟(jì)作物水土資源匹配度，Li為?。▍^(qū)）i經(jīng)濟(jì)作物種植面積。

3）重心遷移模型。 為了反映水土匹配度區(qū)域差異的動(dòng)態(tài)變化，引入重心遷移模型來(lái)表示時(shí)間維度的動(dòng)態(tài)演化特征［20］，重心變動(dòng)反映了水土匹配度的變化軌跡。 計(jì)算公式為［21］

式中：Xt、Yt分別為第t年重心的經(jīng)度和緯度，Xt+1、Yt+1分別為第t+1 年重心的經(jīng)度和緯度，Ri，t為?。▍^(qū)）i第t年水土匹配度，xi、yi分別為省（區(qū)）i地理中心的經(jīng)度和緯度，n為?。▍^(qū)）的數(shù)量，d為重心遷移的距離。

2 結(jié)果與分析

2.1 經(jīng)濟(jì)作物播種面積與其時(shí)空變化特征

1）黃河流域8 個(gè)省（區(qū)）不同類(lèi)型經(jīng)濟(jì)作物播種面積。 由表1 可以看出，谷物等5 類(lèi)經(jīng)濟(jì)作物播種面積逐漸上升，豆類(lèi)等6 類(lèi)經(jīng)濟(jì)作物面積不斷下降，油料播種面積保持相對(duì)穩(wěn)定。 谷物、甜菜、茶葉、蔬菜、水果5 類(lèi)經(jīng)濟(jì)作物播種面積波動(dòng)上升，其中：谷物播種面積上升的原因主要是受到惠農(nóng)政策支持；甜菜播種面積上升歸因于“十三五”以來(lái)我國(guó)開(kāi)始重視甜菜生產(chǎn)［22］；茶葉、蔬菜、水果播種面積增大的原因是，隨著人們生活質(zhì)量的提高，對(duì)健康生活的追求也提高［23］。豆類(lèi)、薯類(lèi)、棉花、麻類(lèi)、甘蔗、煙葉6 類(lèi)作物播種面積不斷下降，其中：豆類(lèi)播種面積下降主要受其經(jīng)濟(jì)效益低的影響［24］；由于育種技術(shù)進(jìn)步，薯類(lèi)單產(chǎn)提高，因此薯類(lèi)播種面積下降；棉花面積變化是種植效益和種植結(jié)構(gòu)調(diào)整的結(jié)果［25］；由于市場(chǎng)需求萎縮、效益不穩(wěn)定，因此麻類(lèi)播種面積下降［26］；甘蔗生產(chǎn)主要集中于河南，受南方甘蔗市場(chǎng)沖擊甘蔗播種面積下降；受新型電子煙快速發(fā)展以及控?zé)熣哂绊懀瑹熑~播種面積下降［27］。 油料種植面積保持相對(duì)穩(wěn)定，歸因于我國(guó)油料以進(jìn)口為主，且河南、山東為油料主產(chǎn)區(qū)［28］，政府給予優(yōu)惠政策支持。

表1 1998—2020 年黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物播種面積演變趨勢(shì)

黃河流域不同類(lèi)型經(jīng)濟(jì)作物播種面積變化趨勢(shì)上、中、下游存在差異。 黃河上游青海、甘肅、寧夏，薯類(lèi)、茶葉、蔬菜、水果4 種經(jīng)濟(jì)作物播種面積呈波動(dòng)上升趨勢(shì)，谷物、豆類(lèi)2 種經(jīng)濟(jì)作物播種面積呈波動(dòng)下降趨勢(shì)，棉花、油料、麻類(lèi)、甘蔗、甜菜、煙葉6 種經(jīng)濟(jì)作物播種面積變化不大。 黃河中游內(nèi)蒙古、陜西、山西，谷物、茶葉、蔬菜、水果4 種經(jīng)濟(jì)作物播種面積呈波動(dòng)上升趨勢(shì)；豆類(lèi)、薯類(lèi)、棉花、麻類(lèi)、甘蔗、甜菜、煙葉7 種經(jīng)濟(jì)作物播種面積呈波動(dòng)下降趨勢(shì)，油料播種面積變化不大。 黃河下游河南、山東，谷物、茶葉、蔬菜、水果4 種經(jīng)濟(jì)作物播種面積呈波動(dòng)上升趨勢(shì)，豆類(lèi)、薯類(lèi)、棉花、油料、麻類(lèi)、甘蔗、煙葉7 種經(jīng)濟(jì)作物播種面積呈波動(dòng)下降趨勢(shì)，甜菜播種面積變化不大。

2）黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物播種面積時(shí)空變化特征。 運(yùn)用ArcGIS10.2 均值分類(lèi)法，將黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物播種面積分為4 個(gè)等級(jí)，分別為Ⅰ級(jí)［486.4，4 054.0］ hm2、Ⅱ級(jí)［4 054.1，7 621.7］hm2、Ⅲ級(jí)［7 621.8，11 189.4］ hm2、Ⅳ級(jí)［11 189.5，14 757.1］ hm2，其中1998 年和2020 年播種面積空間分布見(jiàn)圖1。

圖1 1998 年和2020 年黃河流域8 個(gè)省（區(qū)）不同經(jīng)濟(jì)作物播種面積空間分布

黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物播種面積總體保持穩(wěn)定。 上游甘肅、寧夏、青海3 ?。▍^(qū)）經(jīng)濟(jì)作物播種面積保持在Ⅰ級(jí)；中游陜西播種面積保持在Ⅱ級(jí)，內(nèi)蒙古播種面積由Ⅱ級(jí)上升為Ⅲ級(jí)，山西播種面積則由Ⅱ級(jí)下降為Ⅰ級(jí)；下游河南、山東播種面積保持Ⅳ級(jí)。內(nèi)蒙古經(jīng)濟(jì)作物播種面積增加幅度最大，歸因于為追求經(jīng)濟(jì)利益而掠奪式開(kāi)發(fā)資源以及惠農(nóng)政策的實(shí)施等［29］；山西經(jīng)濟(jì)作物播種面積減小幅度最大，一方面工業(yè)快速發(fā)展導(dǎo)致資源耗竭和環(huán)境破壞，另一方面區(qū)域沙化嚴(yán)重且生態(tài)環(huán)境脆弱，從而實(shí)施了退耕還林還草工程等［30］。

黃河流域經(jīng)濟(jì)作物播種面積呈現(xiàn)東高西低的分布格局。 從播種面積均值來(lái)看，河南（1 396.26 萬(wàn)hm2）＞山東（1 116.03 萬(wàn)hm2）＞內(nèi)蒙古（644.18 萬(wàn)hm2）＞陜西（508.88 萬(wàn)hm2）＞山西（393.10 萬(wàn)hm2）＞甘肅（374.89萬(wàn)hm2）＞寧夏（104.66 萬(wàn)hm2）＞青海（48.47 萬(wàn)hm2）。河南經(jīng)濟(jì)作物播種面積均值最大，歸因于平原面積廣闊且土壤肥沃［31］，有利于農(nóng)業(yè)生產(chǎn)和耕作；青海均值最小，主要原因是地處世界屋脊青藏高原，自然環(huán)境惡劣，不利于作物生長(zhǎng)且品種單一［32］。

2.2 經(jīng)濟(jì)作物水足跡演變趨勢(shì)及時(shí)空變化特征

1）經(jīng)濟(jì)作物水足跡演變趨勢(shì)。 利用式（1）計(jì)算黃河流域12 種經(jīng)濟(jì)作物水足跡，結(jié)果見(jiàn)表2。 黃河流域8 個(gè)?。▍^(qū)）谷物、油料、甜菜、茶葉、蔬菜、水果6 種經(jīng)濟(jì)作物水足跡波動(dòng)上升，豆類(lèi)、薯類(lèi)、棉花、麻類(lèi)、甘蔗、煙葉6 種經(jīng)濟(jì)作物水足跡波動(dòng)下降。 經(jīng)濟(jì)作物水足跡與其播種面積變化趨勢(shì)基本一致（油料除外），原因是，經(jīng)濟(jì)作物水足跡受其播種面積變化影響明顯。

表2 1998—2020 年黃河流域8 省（區(qū)）經(jīng)濟(jì)作物水足跡演變趨勢(shì)

黃河流域經(jīng)濟(jì)作物水足跡變化趨勢(shì)上、中、下游存在差異。 黃河上游青海、甘肅、寧夏，谷物、薯類(lèi)、蔬菜、水果4 種經(jīng)濟(jì)作物水足跡波動(dòng)上升，豆類(lèi)、甜菜、棉花3 種經(jīng)濟(jì)作物水足跡波動(dòng)下降，油料、麻類(lèi)、甘蔗、煙葉、茶葉5 種經(jīng)濟(jì)作物水足跡變化不大。 中游內(nèi)蒙古、陜西、山西，谷物、油料、蔬菜、水果、茶葉5 種經(jīng)濟(jì)作物水足跡波動(dòng)上升，棉花、煙葉2 種經(jīng)濟(jì)作物水足跡波動(dòng)下降，豆類(lèi)、薯類(lèi)、麻類(lèi)、甘蔗、甜菜5 種經(jīng)濟(jì)作物水足跡變化相對(duì)平緩。 下游河南、山東，谷物、茶葉、蔬菜、水果4 種經(jīng)濟(jì)作物水足跡波動(dòng)上升，豆類(lèi)、薯類(lèi)、棉花、麻類(lèi)、甘蔗、煙葉6 種經(jīng)濟(jì)作物水足跡波動(dòng)下降，油料、甜菜2 種經(jīng)濟(jì)作物水足跡變化相對(duì)平緩。

2） 經(jīng)濟(jì)作物水足跡時(shí)空變化特征。 運(yùn)用ArcGIS10.2 均值分類(lèi)法，將黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水足跡分為4 個(gè)等級(jí)，分別為Ⅰ級(jí)［317.8，3 316.2］億m3、Ⅱ級(jí)［3 316.3，6 314.6］億m3、Ⅲ級(jí)［6 314.7，9 313.0］億m3、Ⅳ級(jí)［9 313.1，12 311.4］億m3，水足跡等級(jí)越高表明經(jīng)濟(jì)作物在生長(zhǎng)過(guò)程中所消耗水資源量越多，不利于水資源持續(xù)利用。 1998 年、2020 年黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水足跡空間分布見(jiàn)圖2。

圖2 1998 年、2020 年黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物水足跡空間分布

黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水足跡總體呈增大趨勢(shì)。 從1998 年到2020 年，黃河上游寧夏、青海經(jīng)濟(jì)作物水足跡保持在Ⅰ級(jí)，甘肅經(jīng)濟(jì)作物水足跡由Ⅰ級(jí)上升為Ⅱ級(jí)；中游山西經(jīng)濟(jì)作物水足跡保持在Ⅰ級(jí)，陜西、內(nèi)蒙古經(jīng)濟(jì)作物水足跡則由Ⅰ級(jí)分別上升為Ⅱ級(jí)、Ⅲ級(jí)；黃河下游山東經(jīng)濟(jì)作物水足跡保持在Ⅲ級(jí)，河南經(jīng)濟(jì)作物水足跡則由Ⅲ級(jí)上升為Ⅳ級(jí)。 從黃河流域來(lái)看，經(jīng)濟(jì)作物水足跡呈增大趨勢(shì)，表明經(jīng)濟(jì)作物用水增多，主要原因是水資源利用效率下降［33］；同時(shí)經(jīng)濟(jì)作物播種面積由1998 年的4 475.69 萬(wàn)hm2增加到2020 年的4 817.50 萬(wàn)hm2，增長(zhǎng)率達(dá)到7.6%，水足跡隨之增大。 中游內(nèi)蒙古經(jīng)濟(jì)作物水足跡增幅最大，主要原因是受政策影響播種面積由1998 年的590.86 萬(wàn)hm2增加到2020 年的816.59 萬(wàn)hm2，增長(zhǎng)率達(dá)到了38.2%，播種面積增大導(dǎo)致水足跡增大。

黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水足跡呈現(xiàn)東部和中部高、西部低的分布格局。 從經(jīng)濟(jì)作物水足跡均值來(lái)看，河南（10 049.7 億m3）＞山東（7 894.4 億m3）＞內(nèi)蒙古（4 029.0 億m3）＞陜西（3 373.4 億m3）＞甘肅（2 677.0億m3）＞山西（2 651.5 億m3）＞青海（308.1 億m3）＞寧夏（268.0 億m3）。 河南經(jīng)濟(jì)作物水足跡均值最大，歸因于經(jīng)濟(jì)作物播種面積最大，且種植結(jié)構(gòu)調(diào)整導(dǎo)致種植面積增大，從而造成經(jīng)濟(jì)作物水足跡增大，并導(dǎo)致了非常嚴(yán)重的水生態(tài)赤字［34］；寧夏經(jīng)濟(jì)作物水足跡均值最小，由于寧夏氣候干燥且年均降水量?jī)H166～647 mm，自然生態(tài)環(huán)境脆弱，不利于耕作，因此經(jīng)濟(jì)作物水足跡最?。?5］。 從水資源安全考慮，應(yīng)該發(fā)展高效節(jié)水農(nóng)業(yè)，因地制宜發(fā)展蔬菜、水果等低耗水經(jīng)濟(jì)作物，適量進(jìn)口水密集型產(chǎn)品，實(shí)現(xiàn)水資源合理利用。

2.3 經(jīng)濟(jì)作物水土資源匹配度時(shí)空分析

利用式（2）計(jì)算黃河流域8 ?。▍^(qū)）經(jīng)濟(jì)作物水土資源匹配度，運(yùn)用ArcGIS10.2 均值分類(lèi)法劃分為4 個(gè)等級(jí)，分別為Ⅰ級(jí)［0.48，0.60］萬(wàn)m3／hm2、Ⅱ級(jí)［0.61，0.72］萬(wàn)m3／hm2、Ⅲ級(jí)［0.73，0.84］萬(wàn)m3／hm2、Ⅳ級(jí)［0.85，0.96］萬(wàn)m3／hm2，進(jìn)而分析其空間格局變化，見(jiàn)圖3。 水土匹配度是指區(qū)域內(nèi)單位面積耕地可擁有的經(jīng)濟(jì)作物水足跡，水土匹配度越小，表明單位面積耕地經(jīng)濟(jì)作物生產(chǎn)所消耗的水資源量越少，農(nóng)業(yè)生產(chǎn)水土資源可持續(xù)利用越好；水土匹配度越大，匹配程度越差［36］。

圖3 1998 年、2020 年黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物水土匹配度

黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物水土匹配度呈增大趨勢(shì)。 黃河上游寧夏、甘肅經(jīng)濟(jì)作物水土匹配度分別由Ⅰ級(jí)、Ⅱ級(jí)增大至Ⅳ級(jí)，青海經(jīng)濟(jì)作物水土匹配度由Ⅰ級(jí)增大至Ⅱ級(jí)；中游內(nèi)蒙古、山西經(jīng)濟(jì)作物水土匹配度分別由Ⅰ級(jí)、Ⅱ級(jí)增大至Ⅲ級(jí)，陜西經(jīng)濟(jì)作物水土匹配度由Ⅰ級(jí)增大為Ⅳ級(jí)；下游河南、山東經(jīng)濟(jì)作物水土匹配度由Ⅰ級(jí)增大至Ⅲ級(jí)。 黃河流域水土匹配度均有不同程度增大，說(shuō)明單位耕地面積經(jīng)濟(jì)作物生產(chǎn)消耗的水資源量增多，主要原因是農(nóng)業(yè)水資源生產(chǎn)效率不高，且受技術(shù)效率影響較大，而較低的技術(shù)效率限制了水資源生產(chǎn)效率的提高［37］。 其中寧夏和陜西經(jīng)濟(jì)作物水土匹配度變化最大，主要原因是，農(nóng)業(yè)是寧夏主導(dǎo)產(chǎn)業(yè)，農(nóng)業(yè)種植以高耗水作物為主；陜西經(jīng)濟(jì)作物主要集中在南方山區(qū)，農(nóng)田水利建設(shè)難以實(shí)施，且農(nóng)業(yè)用水占60%以上，限制了用水效率的提高。

黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水土匹配度呈現(xiàn)中南部高、北部和西部低的分布格局。 從均值來(lái)看，寧夏（0.734 萬(wàn)m3／hm2）＞河南（0.717 萬(wàn)m3／hm2）＞甘肅（0.712 萬(wàn)m3／hm2）＞山東（0.708 萬(wàn)m3／hm2）＞山西（0.676萬(wàn)m3／hm2）＞陜西（0.662 萬(wàn)m3／hm2）＞青海（0.638萬(wàn)m3／hm2）＞內(nèi)蒙古（0.613 萬(wàn)m3／hm2）。 寧夏水土匹配度均值最大，主要原因是寧夏氣候干燥且年均降水量?jī)H166～647 mm，區(qū)域水資源整體稀缺，水生態(tài)環(huán)境脆弱，而寧夏平原是重要灌溉農(nóng)業(yè)區(qū)，水資源生產(chǎn)率低，經(jīng)濟(jì)效益不高［34，38］；內(nèi)蒙古水土匹配度均值最小，主要原因是內(nèi)蒙古地區(qū)為以雨養(yǎng)農(nóng)業(yè)或牧草業(yè)為主的旱作區(qū)［35］。

2.4 黃河流域經(jīng)濟(jì)作物水土匹配度重心軌跡演化

利用式（3）、式（4）計(jì)算黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水土匹配度重心，以探討其空間演變特征，見(jiàn)圖4。 黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物水土匹配度重心軌跡呈現(xiàn)由東北向西南移動(dòng)的趨勢(shì)，1998—2020 年年均移動(dòng)距離為0.54 km，累計(jì)移動(dòng)307.12 km，表明西南部經(jīng)濟(jì)作物水土匹配度增大，這與高蕓等［5］研究結(jié)果一致。 不同年份經(jīng)濟(jì)作物水土匹配度空間變化具有不平衡性，其中2018 年經(jīng)濟(jì)作物水土匹配度重心向東偏南方向移動(dòng)了約48.65 km，移動(dòng)距離最大；2007 年經(jīng)濟(jì)作物水土匹配度重心移動(dòng)距離最小，向南偏西方向移動(dòng)了約1.53 km。

圖4 1998—2020 年黃河流域8 個(gè)省（區(qū)）經(jīng)濟(jì)作物匹配度重心軌跡

黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水土匹配度重心呈現(xiàn)4 個(gè)階段的變化特征，第一階段為1998—2009年，經(jīng)濟(jì)作物水土匹配度重心位于（109°2'38″E，37°20'38″N），穩(wěn)定在黃河中游范圍；第二階段為2009—2013 年，經(jīng)濟(jì)作物水土匹配度重心向西偏南方向遷移，偏移距離約為43.51 km，主要原因是，黃河上游經(jīng)濟(jì)作物水土匹配度增大；第三階段為2013—2018年，經(jīng)濟(jì)作物水土匹配度重心開(kāi)始向東偏南移動(dòng)，偏移距離約為96.38 km，主要原因是，黃河下游河南、山東經(jīng)濟(jì)作物水土匹配度增大；第四階段為2018—2020 年，經(jīng)濟(jì)作物水土匹配度重心向西偏北方向遷移，偏移距離約為63.63 km，歸因于黃河上游經(jīng)濟(jì)作物水土匹配度增大。

3 結(jié)論

基于水足跡視角，運(yùn)用水土匹配度、重心遷移模型，結(jié)合ArcGIS10.2 軟件，分析了黃河流域8 個(gè)?。▍^(qū)）1998—2020 年不同類(lèi)型經(jīng)濟(jì)作物播種面積及其水足跡，探討了經(jīng)濟(jì)作物水土資源時(shí)空匹配格局及重心演化特征。

1）黃河流域谷物等5 類(lèi)經(jīng)濟(jì)作物播種面積逐漸上升，豆類(lèi)等6 類(lèi)經(jīng)濟(jì)作物播種面積不斷下降，油料播種面積則保持相對(duì)穩(wěn)定；不同類(lèi)型經(jīng)濟(jì)作物播種面積上、中、下游存在差異。 黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物播種面積總體保持穩(wěn)定，呈現(xiàn)東高西低分布格局。為了做到以水定城、以水定地、以水定人、以水定產(chǎn)，把水資源作為最大的剛性約束，要不斷優(yōu)化種植結(jié)構(gòu)，抑制大面積種植高耗水作物，不斷發(fā)展節(jié)水技術(shù)，實(shí)施全流域節(jié)水行動(dòng)，推動(dòng)用水方式由粗放向節(jié)約集約轉(zhuǎn)變。

2）谷物等6 類(lèi)經(jīng)濟(jì)作物水足跡波動(dòng)上升，豆類(lèi)等6 類(lèi)經(jīng)濟(jì)作物水足跡波動(dòng)下降；不同經(jīng)濟(jì)作物水足跡上、中、下游存在差異。 黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水足跡總體呈現(xiàn)增大趨勢(shì)，呈現(xiàn)東部和中部高、西部低的分布格局。 為緩解黃河流域水足跡增大趨勢(shì)，要大力推廣資源節(jié)約型農(nóng)業(yè)生產(chǎn)技術(shù)，提高農(nóng)機(jī)和水資源利用效率；加強(qiáng)科技支撐，提高水資源生產(chǎn)效率，促進(jìn)農(nóng)業(yè)高效發(fā)展；推廣節(jié)水灌溉新技術(shù)，改善灌溉條件，發(fā)展綠色、高效節(jié)水農(nóng)業(yè)。

3）黃河流域8 個(gè)?。▍^(qū)）經(jīng)濟(jì)作物水土匹配度有增大趨勢(shì)，呈現(xiàn)中南部高、北部和西部低的分布格局；水土匹配度重心由東北向西南方向移動(dòng)，表現(xiàn)為1998—2009 年處于相對(duì)穩(wěn)定階段、2009—2013 年向西偏南方向移動(dòng)、2013—2018 年向東偏南方向移動(dòng)、2018—2020 年向西偏北方向移動(dòng)。 黃河上、中、下游要發(fā)揮各自?xún)?yōu)勢(shì)，推動(dòng)流域協(xié)同合作，不斷優(yōu)化水土配置，逐步縮小區(qū)域差異；黃河流域上游地區(qū)要注重生態(tài)環(huán)境保護(hù)，提高用水效率；中游地區(qū)面對(duì)水土流失的困境，要加強(qiáng)水土保持工作，著力發(fā)展環(huán)境友好型農(nóng)業(yè)；下游地區(qū)要推廣農(nóng)業(yè)節(jié)水技術(shù)；同時(shí)建議實(shí)施水系聯(lián)通、跨流域調(diào)水等工程，以改善因水土不平衡造成的結(jié)構(gòu)性缺水［5］。