摘 要：為探討黃河流域農(nóng)業(yè)綠色發(fā)展的影響因素和空間差異、給黃河流域農(nóng)業(yè)高質(zhì)量發(fā)展提供理論支撐和決策依據(jù)，建立了包括發(fā)展高效、生態(tài)友好和資源節(jié)約３ 個維度的農(nóng)業(yè)綠色發(fā)展水平評價指標(biāo)體系，基于黃河流域６８ 個地級市２０１１—２０２１ 年的面板數(shù)據(jù)、運用全局熵值法評估黃河流域農(nóng)業(yè)綠色發(fā)展水平，采用莫蘭指數(shù)對其空間分異格局進行分析，在此基礎(chǔ)上從經(jīng)濟、社會、自然三方面設(shè)置產(chǎn)業(yè)結(jié)構(gòu)、種植結(jié)構(gòu)、工業(yè)化水平、城鎮(zhèn)化水平、地形起伏度、年降水量６ 個指標(biāo)，采用多尺度地理加權(quán)回歸模型實證分析黃河流域農(nóng)業(yè)綠色發(fā)展的影響因素及其空間差異。結(jié)果表明：１）２０１１—２０２１ 年黃河流域農(nóng)業(yè)綠色發(fā)展水平呈上升趨勢，但研究期末整體水平仍較低，下游地區(qū)農(nóng)業(yè)綠色發(fā)展水平明顯高于中上游地區(qū)，發(fā)展高效維度在上游和下游地區(qū)上升趨勢明顯，生態(tài)友好和資源節(jié)約維度在研究期上升幅度較?。唬玻S河流域農(nóng)業(yè)綠色發(fā)展水平存在明顯的空間正相關(guān)性，高－高集聚和低－低集聚地級市相對較多；３）農(nóng)業(yè)綠色發(fā)展各影響因素的空間異質(zhì)性顯著且影響強度不同，其中產(chǎn)業(yè)結(jié)構(gòu)、種植結(jié)構(gòu)、工業(yè)化水平、地形起伏度的影響空間差異較小，城鎮(zhèn)化水平和年降水量的影響空間差異較大，產(chǎn)業(yè)結(jié)構(gòu)、工業(yè)化水平和地形起伏度有負(fù)面影響，種植結(jié)構(gòu)、城鎮(zhèn)化水平和年降水量有正面影響，影響強度大小為地形起伏度＞工業(yè)化水平＞城鎮(zhèn)化水平＞產(chǎn)業(yè)結(jié)構(gòu)＞種植結(jié)構(gòu)＞年降水量。提出了轉(zhuǎn)變農(nóng)業(yè)發(fā)展理念、加強區(qū)域協(xié)作、實施差異化農(nóng)業(yè)綠色發(fā)展策略等建議。

關(guān)鍵詞：農(nóng)業(yè)綠色發(fā)展；時空格局；影響因素；多尺度地理加權(quán)回歸模型；黃河流域

中圖分類號：Ｆ４９；Ｆ４２６．２；ＴＶ８８２．１ 文獻標(biāo)志碼：Ａ ｄｏｉ：１０．３９６９／ ｊ．ｉｓｓｎ．１０００－１３７９．２０２５．０４．００２

引用格式：周廣亮，翟曉格，張紹陽．黃河流域農(nóng)業(yè)綠色發(fā)展時空格局及影響因素研究［Ｊ］．人民黃河，２０２５，４７（４）：７－１４．

０ 引言

黨的十八大以來，我國農(nóng)業(yè)發(fā)展逐漸從數(shù)量增長型向質(zhì)量提升型轉(zhuǎn)變，生態(tài)保護、綠色發(fā)展成為農(nóng)業(yè)高質(zhì)量發(fā)展的重要特征［１－４］ ，但農(nóng)業(yè)面源污染和生態(tài)退化等問題還沒有從根本上得到遏制，農(nóng)業(yè)支持保護機制還需要進一步優(yōu)化和完善。黨的二十大報告明確提出“推動綠色發(fā)展，促進人與自然和諧共生”，２０２１ 年農(nóng)業(yè)農(nóng)村部等６ 部門制定了《“十四五”全國農(nóng)業(yè)綠色發(fā)展規(guī)劃》，使農(nóng)業(yè)綠色發(fā)展迎來新機遇。黃河流域是我國重要的農(nóng)業(yè)生產(chǎn)區(qū)和生態(tài)保護區(qū)，但農(nóng)業(yè)生產(chǎn)方式長期以高消耗、高污染的粗放型為主，農(nóng)業(yè)面源污染和生態(tài)退化等問題突出。研究黃河流域農(nóng)業(yè)綠色發(fā)展格局及其影響因素、探尋提升農(nóng)業(yè)綠色發(fā)展水平的路徑，對破解黃河流域高質(zhì)量發(fā)展的資源環(huán)境約束、鞏固生態(tài)屏障等具有重要意義。

農(nóng)業(yè)綠色發(fā)展是當(dāng)前學(xué)術(shù)界研究的熱點之一，有關(guān)學(xué)者在其內(nèi)涵闡釋、評價指標(biāo)體系及影響因素等方面開展了廣泛研究。在內(nèi)涵闡釋方面，李福奪等［５］ 從狹義角度進行研究，認(rèn)為農(nóng)業(yè)綠色發(fā)展是建立在可持續(xù)農(nóng)業(yè)實踐之上的農(nóng)業(yè)體系綠色化進程；趙會杰等［６］從廣義角度進行研究，認(rèn)為農(nóng)業(yè)綠色發(fā)展以資源和環(huán)境為載體、以制度建設(shè)和機制創(chuàng)新為保障，通過合理利用農(nóng)業(yè)資源實現(xiàn)經(jīng)濟、社會和環(huán)境的和諧統(tǒng)一；蓋美等［７］ 從系統(tǒng)論視角進行研究，認(rèn)為農(nóng)業(yè)綠色發(fā)展以綠色技術(shù)和綠色產(chǎn)品為核心，結(jié)合先進的技術(shù)設(shè)施和管理理念等，形成促進人類經(jīng)濟社會全面可持續(xù)發(fā)展的農(nóng)業(yè)發(fā)展模式。在評價指標(biāo)體系方面，王磊等［８］ 將資源節(jié)約、環(huán)境友好、產(chǎn)出高效視為衡量農(nóng)業(yè)綠色發(fā)展的關(guān)鍵指標(biāo)；樊勝岳等［９］ 在不使用外部環(huán)境變量的情況下，從農(nóng)業(yè)生產(chǎn)、生態(tài)、經(jīng)濟發(fā)展的角度評估農(nóng)業(yè)綠色發(fā)展水平；梁俊芬等［１０］ 將綠色技術(shù)、綠色產(chǎn)品的投入和產(chǎn)出引入農(nóng)業(yè)綠色發(fā)展評價指標(biāo)體系；熊素等［１１］ 認(rèn)為，在“雙碳”目標(biāo)背景下碳排放評價是農(nóng)業(yè)綠色發(fā)展評價的重要內(nèi)容；劉惠良等［１２］ 從鄉(xiāng)村振興視角進行了農(nóng)業(yè)綠色發(fā)展水平測度研究。在評價方法方面，已有研究采用的方法主要有熵權(quán)ＴＯＰＳＩＳ 法［１３］ 、主成分分析法［１４］ 、ＤＰＳＩＲ 模型［１５］ 、組合賦權(quán)法［１６］ 、層次分析法［１７］ ，一些研究采用ＤＥＡ－Ｍａｌｍｑｕｉｓｔ 模型［１８］ 、ＳＢＭ 超效率模型［１９］ 、基于混合距離函數(shù)的全局Ｍａｌｍｑｕｉｓｔ－Ｌｕｅｎｂｅｒｇｅｒ 指數(shù)［２０］ 等。在影響因素方面， 許多學(xué)者［２１－２４］ 認(rèn)為農(nóng)業(yè)產(chǎn)業(yè)集聚、居民收入水平和科技水平、財政支農(nóng)措施、合理的環(huán)境規(guī)制等對農(nóng)業(yè)綠色發(fā)展具有顯著影響。

鑒于不同地區(qū)資源稟賦的差異使得農(nóng)業(yè)綠色發(fā)展?fàn)顩r不同，即區(qū)域資源稟賦特征也是影響農(nóng)業(yè)綠色發(fā)展的重要因素，本研究基于黃河流域６８ 個地級市２０１１—２０２１ 年的面板數(shù)據(jù)、運用全局熵值法評估黃河流域農(nóng)業(yè)綠色發(fā)展水平，采用莫蘭指數(shù)對其時空分異格局進行分析，在此基礎(chǔ)上采用多尺度地理加權(quán)回歸模型探討黃河流域農(nóng)業(yè)綠色發(fā)展的影響因素及其空間差異，以期為提升黃河流域農(nóng)業(yè)綠色發(fā)展水平提供決策支持。

１ 研究設(shè)計

１．１ 樣本市選取

黃河流域是跨越我國東、中、西部地區(qū)的重要生態(tài)走廊，流域內(nèi)各地區(qū)生態(tài)環(huán)境等差異較大，自然生態(tài)系統(tǒng)復(fù)雜脆弱，流域高質(zhì)量發(fā)展面臨的生態(tài)環(huán)境約束較為突出。本研究以黃河流域各地級市為樣本單元，剔除數(shù)據(jù)缺失嚴(yán)重的部分地級市，最終確定了６８ 個樣本市。同時，以流域自然稟賦為基礎(chǔ)，借鑒相關(guān)研究，以內(nèi)蒙古托克托縣的河口鎮(zhèn)、河南鄭州市的桃花峪分別為上中游分界點、中下游分界點，把黃河流域劃分為上、中、下游３ 個區(qū)域，各區(qū)域樣本市統(tǒng)計見表１。

１．２ 研究方法

１）全局熵值法。熵值法根據(jù)各指標(biāo)的信息量來決定其權(quán)重，信息量越大熵值越小，反之，信息量越小熵值越大，是一種客觀賦權(quán)的多指標(biāo)綜合評價方法。為了比較黃河流域不同年份農(nóng)業(yè)綠色發(fā)展水平、使分析結(jié)果更加科學(xué)，借鑒趙會杰等［６］ 的研究，加入時間變量對熵值法進行改進（稱為全局熵值法），測算農(nóng)業(yè)綠色發(fā)展水平（用農(nóng)業(yè)綠色發(fā)展指數(shù)表示）。

２）空間相關(guān)性分析。采用全局莫蘭指數(shù)分析黃河流域農(nóng)業(yè)綠色發(fā)展水平的空間相關(guān)性，采用局部莫蘭指數(shù)分析各樣本市農(nóng)業(yè)綠色發(fā)展的集聚特征及其分類，并運用ＡｒｃＭａｐ １０．２ 軟件對分析結(jié)果進行可視化展示。

３） 多尺度地理加權(quán)回歸模型。地理加權(quán)回歸（Ｇｅｏｇｒａｐｈｉｃａｌｌｙ Ｗｅｉｇｈｔｅｄ Ｒｅｇｒｅｓｓｉｏｎ，ＧＷＲ） 是一種局部空間回歸分析方法，基本思想是在回歸分析中引入空間權(quán)重，根據(jù)地理位置對變量值進行加權(quán)（常見的空間權(quán)函數(shù)有距離閾值法、距離反比法、Ｇａｕｓｓ 函數(shù)法等），探索研究對象在某一尺度下的空間變化及驅(qū)動因素。多尺度地理加權(quán)回歸（ＭＧＷＲ）是對ＧＷＲ 的改進，它允許不同變量具有不同的帶寬，即考慮各變量的尺度效應(yīng)，通過優(yōu)化每個特征變量的帶寬，捕捉研究對象空間關(guān)系的多尺度特性，比ＧＷＲ 更加實用，其模型形式如下：

式中：下標(biāo)ｉ 為樣本市編號、ｊ 為解釋變量編號，ｙｉ為被解釋變量（本研究中為農(nóng)業(yè)綠色發(fā)展指數(shù)），ｘｉｊ 為解釋變量， βｂｉｊ 為局域變量回歸系數(shù)（其中ｂｉｊ為第ｊ 個解釋變量的帶寬），（ｕｉ ，ｖｉ ）為樣本市空間坐標(biāo)，ｋ 為局域解釋變量的數(shù)量，εｉ為隨機擾動項。

１．３ 農(nóng)業(yè)綠色發(fā)展水平評價指標(biāo)體系構(gòu)建

農(nóng)業(yè)綠色發(fā)展既要關(guān)注農(nóng)業(yè)產(chǎn)出水平的高低，又要關(guān)注農(nóng)業(yè)生產(chǎn)過程與生態(tài)環(huán)境的協(xié)調(diào)性以及資源利用的永續(xù)性。為此，本研究借鑒查建平等［２５］ 的研究，根據(jù)農(nóng)業(yè)綠色發(fā)展的內(nèi)涵，遵循評價指標(biāo)體系構(gòu)建的系統(tǒng)性、數(shù)據(jù)可獲得性和連續(xù)性等原則，建立了包括發(fā)展高效、生態(tài)友好和資源節(jié)約３ 個維度的農(nóng)業(yè)綠色發(fā)展水平評價指標(biāo)體系（見表２）。

１）發(fā)展高效。發(fā)展高效是農(nóng)業(yè)綠色發(fā)展的根本目標(biāo)。農(nóng)業(yè)綠色發(fā)展的目的在于提高農(nóng)業(yè)生產(chǎn)效率，增加優(yōu)質(zhì)、安全、綠色的農(nóng)產(chǎn)品供給，增進人民生活福祉。因此，本研究選取單位播種面積農(nóng)機總動力、土地產(chǎn)出率、人均農(nóng)業(yè)生產(chǎn)總值、農(nóng)村居民人均可支配收入４ 個指標(biāo)來反映發(fā)展高效程度。

２）生態(tài)友好。生態(tài)友好是農(nóng)業(yè)綠色發(fā)展的基本特征。農(nóng)業(yè)綠色發(fā)展需要以良好的土地、水、生物等生態(tài)環(huán)境資源為載體，在農(nóng)業(yè)生產(chǎn)過程中應(yīng)盡量減少對生態(tài)環(huán)境的污染。因此，本研究選取化肥施用強度、農(nóng)藥施用強度、農(nóng)膜使用強度、生活垃圾無害化處理率４個指標(biāo)來反映生態(tài)友好程度。

３）資源節(jié)約。資源節(jié)約是農(nóng)業(yè)綠色發(fā)展的現(xiàn)實要求。黃河流域大部分區(qū)域位于干旱半干旱區(qū)，農(nóng)業(yè)生產(chǎn)面臨水資源緊缺等問題，為了保障農(nóng)業(yè)綠色可持續(xù)發(fā)展，應(yīng)充分提升資源節(jié)約利用水平。因此，本研究選取單位農(nóng)業(yè)產(chǎn)值耗水量、單位農(nóng)業(yè)產(chǎn)值耗電量２ 個指標(biāo)來反映資源節(jié)約程度。

１．４ 影響因素指標(biāo)設(shè)置

黃河流域農(nóng)業(yè)綠色發(fā)展受經(jīng)濟、社會、自然三方面因素的影響，各因素指標(biāo)設(shè)置見表３。

１）經(jīng)濟因素。采用產(chǎn)業(yè)結(jié)構(gòu)、種植結(jié)構(gòu)、工業(yè)化水平３ 個指標(biāo)來反映影響農(nóng)業(yè)綠色發(fā)展的經(jīng)濟因素。在產(chǎn)業(yè)結(jié)構(gòu)中，非農(nóng)產(chǎn)業(yè)的發(fā)展可能會帶動農(nóng)業(yè)生產(chǎn)技術(shù)和農(nóng)業(yè)社會化服務(wù)的升級，包括先進的農(nóng)業(yè)機械、高效的化肥以及科學(xué)管理方法等，但也可能導(dǎo)致很多原本投向農(nóng)業(yè)領(lǐng)域的勞動力、資金等要素轉(zhuǎn)向非農(nóng)領(lǐng)域，從而對農(nóng)業(yè)綠色發(fā)展產(chǎn)生負(fù)面影響。種植結(jié)構(gòu)也會影響農(nóng)業(yè)綠色發(fā)展，通過合理調(diào)整種植結(jié)構(gòu)，大力發(fā)展特色高效經(jīng)濟作物，可拓寬農(nóng)民增收途徑，還能緩解土地退化、生物多樣性喪失和水資源浪費等問題，從而促進農(nóng)業(yè)綠色發(fā)展。工業(yè)化有助于為延伸農(nóng)業(yè)產(chǎn)業(yè)鏈提供資金和技術(shù)支持，但也可能因侵占農(nóng)業(yè)用地、加劇生態(tài)破壞等而阻礙農(nóng)業(yè)綠色發(fā)展。

２）社會因素。采用城鎮(zhèn)化水平這一指標(biāo)來反映影響農(nóng)業(yè)綠色發(fā)展的社會因素。城鎮(zhèn)化發(fā)展可以促進農(nóng)業(yè)現(xiàn)代化和農(nóng)業(yè)生產(chǎn)方式的轉(zhuǎn)變，但城市“虹吸”效應(yīng)往往又使得農(nóng)村勞動力向城鎮(zhèn)轉(zhuǎn)移、土地利用結(jié)構(gòu)變化（如農(nóng)業(yè)用地轉(zhuǎn)為非農(nóng)業(yè)用地），對農(nóng)業(yè)綠色發(fā)展產(chǎn)生負(fù)面影響。

３）自然因素。采用地形起伏度、年降水量２ 個指標(biāo)來反映影響農(nóng)業(yè)綠色發(fā)展的自然因素。地形起伏度不僅通過坡度、坡向、海拔等因素影響氣溫和太陽光照，而且會影響農(nóng)業(yè)生產(chǎn)的規(guī)?；蜋C械化程度，進而影響農(nóng)業(yè)綠色發(fā)展。黃河流域大部分地區(qū)水資源短缺，年降水量的多少直接影響著土壤質(zhì)量、灌溉成本以及農(nóng)業(yè)生態(tài)環(huán)境等，從而對農(nóng)業(yè)綠色發(fā)展產(chǎn)生影響。

１．５ 數(shù)據(jù)來源

本研究所用面板數(shù)據(jù)主要來源于《中國城市統(tǒng)計年鑒》、各?。▍^(qū)）《統(tǒng)計年鑒》《農(nóng)村統(tǒng)計年鑒》《水資源公報》、各樣本市《統(tǒng)計年鑒》《國民經(jīng)濟和社會發(fā)展統(tǒng)計公報》等（其中涉及價格的數(shù)據(jù)均以２０１０ 年為基期剔除價格因素的影響），ＤＥＭ 數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心（ｈｔｔｐｓ：／／ ｗｗｗ．ｒｅｓｄｃ．ｃｎ／ ），對缺失的少量數(shù)據(jù)采用線性插值法或相鄰年份值替換法進行補充。

２ 農(nóng)業(yè)綠色發(fā)展水平測算結(jié)果及分析

２．１ 農(nóng)業(yè)綠色發(fā)展水平測算結(jié)果

采用全局熵值法測算６８ 個樣本市農(nóng)業(yè)綠色發(fā)展指數(shù)，結(jié)果見表４（限于篇幅，僅展示了３ 個代表年份農(nóng)業(yè)綠色發(fā)展指數(shù)大小排序前５ 位、后５ 位的樣本市測算結(jié)果及流域、各區(qū)域均值）。總體而言，２０１１—２０２１ 年黃河流域農(nóng)業(yè)綠色發(fā)展水平呈上升趨勢（流域平均農(nóng)業(yè)綠色發(fā)展指數(shù)由０．２３９ 提高到０．３９７），但研究期末農(nóng)業(yè)綠色發(fā)展水平仍較低；從各樣本市看，農(nóng)業(yè)綠色發(fā)展水平差異較大（研究期末綠色發(fā)展指數(shù)最大值０．７５６ 約為最小值０．１５１ 的５ 倍）；從不同區(qū)域看，下游地區(qū)依托優(yōu)越的地理區(qū)位和資源環(huán)境條件，農(nóng)業(yè)綠色發(fā)展水平明顯較中、上游地區(qū)高。

２．２ 農(nóng)業(yè)綠色發(fā)展水平逐年變化情況

黃河流域農(nóng)業(yè)綠色發(fā)展指數(shù)及各維度發(fā)展指數(shù)逐年演化情況見圖１。從流域整體上看，農(nóng)業(yè)綠色發(fā)展水平在２０１１—２０１５ 年緩慢上升、２０１６—２０１７ 年小幅下降、２０１８ 年后快速上升；從各區(qū)域看，上游地區(qū)農(nóng)業(yè)綠色發(fā)展水平在研究期穩(wěn)步上升，而中、下游地區(qū)農(nóng)業(yè)綠色發(fā)展水平與流域整體變化情況類似（在２０１６—２０１７ 年小幅下降、２０１８ 年后恢復(fù)上升趨勢）；從各維度來看，生態(tài)友好指數(shù)和資源節(jié)約指數(shù)在研究期保持穩(wěn)定，而發(fā)展高效指數(shù)明顯大于生態(tài)友好指數(shù)和資源節(jié)約指數(shù)且與農(nóng)業(yè)綠色發(fā)展指數(shù)變化情況一致，表明研究期黃河流域農(nóng)業(yè)綠色發(fā)展水平的變化主要受發(fā)展高效維度的影響，農(nóng)業(yè)綠色發(fā)展水平較低的主要原因是生態(tài)友好水平和資源節(jié)約水平較低，減少農(nóng)藥和化肥施用量、破解“白色污染”問題、提高資源利用效率是黃河流域各地區(qū)農(nóng)業(yè)綠色發(fā)展的關(guān)鍵。

２．３ 農(nóng)業(yè)綠色發(fā)展水平的空間分異特征

２０１１—２０２１ 年黃河流域農(nóng)業(yè)綠色發(fā)展水平全局莫蘭指數(shù)從０．５２７ 波動下降至０．４７６（見表５）、均通過了置信水平為９９％的顯著性檢驗，其值一直大于０，表明黃河流域農(nóng)業(yè)綠色發(fā)展水平存在明顯的空間正相關(guān)性，即各樣本市農(nóng)業(yè)綠色發(fā)展水平存在集聚效應(yīng)，其值在研究期呈波動下降趨勢表明集聚效應(yīng)逐漸減弱，但“強者恒強”的兩極分化現(xiàn)象仍較為明顯。

圖２ 為３ 個代表年份黃河流域農(nóng)業(yè)綠色發(fā)展水平局部空間集聚情況（圖中Ｈ－Ｈ 集聚、Ｈ－Ｌ 集聚、Ｌ－Ｈ集聚、Ｌ－Ｌ 集聚分別表示高－高集聚、高－低集聚、低－高集聚、低－低集聚）?？傮w上看，黃河流域農(nóng)業(yè)綠色發(fā)展水平局部空間集聚情況有高－高集聚和低－低集聚兩種，無高－低集聚和低－高集聚，農(nóng)業(yè)綠色發(fā)展水平空間分布相對穩(wěn)定；從各樣本市看，高－高集聚的樣本市主要分布在下游地區(qū)（有濟南、青島、淄博、煙臺、開封、濮陽等），低－低集聚的樣本市主要分布在上游（有蘭州、白銀、定西、固原、中衛(wèi)等）及中游（有天水、平?jīng)?、慶陽等）；從集聚情況演變情況看，２０１６ 年以后高－高集聚區(qū)減小并向下游山東省部分地區(qū)萎縮。

３ 農(nóng)業(yè)綠色發(fā)展水平影響因素分析

由上文可知，黃河流域農(nóng)業(yè)綠色發(fā)展水平存在明顯的空間差異。采用ＭＧＷＲ 模型從空間視角探究黃河流域農(nóng)業(yè)綠色發(fā)展水平影響因素，表６ 為ＭＧＷＲ 模型回歸結(jié)果的描述性統(tǒng)計。

從各影響因素的帶寬（作用尺度）看，產(chǎn)業(yè)結(jié)構(gòu)、種植結(jié)構(gòu)、工業(yè)化水平、地形起伏度的帶寬為６６，接近全局尺度（６８），表明這４ 個因素對農(nóng)業(yè)綠色發(fā)展影響的空間差異不大，即空間異質(zhì)性較弱；城鎮(zhèn)化水平和年降水量的帶寬分別為５６ 和４８，說明這２ 個因素對農(nóng)業(yè)綠色發(fā)展的影響具有較大空間差異，即有較強的空間異質(zhì)性。ＭＧＷＲ 模型回歸系數(shù)平均值反映了各因素對農(nóng)業(yè)綠色發(fā)展影響力的平均情況，從其屬性來看，產(chǎn)業(yè)結(jié)構(gòu)、工業(yè)化水平和地形起伏度的影響為負(fù)，其他因素的影響為正；從其絕對值來看，影響最大的因素是地形起伏度，其次為工業(yè)化水平和城鎮(zhèn)化水平，年降水量的影響最弱。

各因素回歸系數(shù)的空間差異見圖３，據(jù)此對各因素影響農(nóng)業(yè)綠色發(fā)展的空間異質(zhì)性分述如下：

１）產(chǎn)業(yè)結(jié)構(gòu)。產(chǎn)業(yè)結(jié)構(gòu)對黃河中上游地區(qū)農(nóng)業(yè)綠色發(fā)展具有顯著的負(fù)面影響，原因是產(chǎn)業(yè)結(jié)構(gòu)調(diào)整必然伴隨著各種資源的轉(zhuǎn)移，服務(wù)業(yè)等非農(nóng)產(chǎn)業(yè)的發(fā)展可能會與農(nóng)業(yè)進行土地、水資源占（利）用的競爭，影響農(nóng)業(yè)生產(chǎn)的可持續(xù)性，同時，產(chǎn)業(yè)結(jié)構(gòu)調(diào)整可能導(dǎo)致農(nóng)村勞動力向城市轉(zhuǎn)移，影響農(nóng)業(yè)生產(chǎn)的效率和規(guī)模。因此，非農(nóng)產(chǎn)業(yè)占比較高會給農(nóng)業(yè)綠色發(fā)展帶來一定壓力。

２）種植結(jié)構(gòu)。種植結(jié)構(gòu)對黃河流域農(nóng)業(yè)綠色發(fā)展的影響呈現(xiàn)自北向南逐漸減弱趨勢，具體而言，與青海、甘肅、河南、山東等地相比，內(nèi)蒙古、山西和陜西農(nóng)業(yè)種植結(jié)構(gòu)的調(diào)整對農(nóng)業(yè)綠色發(fā)展水平的提升效果更為顯著。內(nèi)蒙古因地制宜發(fā)展區(qū)域特色農(nóng)產(chǎn)品，不僅提升了農(nóng)業(yè)產(chǎn)出效益，而且減少了農(nóng)藥、化肥等的污染；山西和陜西水資源相對匱乏，調(diào)整種植結(jié)構(gòu)可以有效改善生態(tài)環(huán)境、提高水資源利用效率，使農(nóng)業(yè)綠色發(fā)展水平顯著提升。而青海和甘肅的部分地區(qū)農(nóng)業(yè)種植規(guī)模較小、種植結(jié)構(gòu)單一，加之受自然資源限制，難以滿足經(jīng)濟價值較高農(nóng)作物的種植要求，種植結(jié)構(gòu)調(diào)整對農(nóng)業(yè)綠色發(fā)展影響較弱。

３）工業(yè)化水平。工業(yè)化水平對黃河流域農(nóng)業(yè)綠色發(fā)展的影響強度呈現(xiàn)由東北向西南逐漸減弱的趨勢，其中對中下游地區(qū)農(nóng)業(yè)綠色發(fā)展具有顯著的負(fù)面影響，這與工業(yè)化發(fā)展的區(qū)域不平衡性有關(guān)，位于黃河中下游的山西、河南等省份工業(yè)化水平的提高給環(huán)境和資源帶來巨大壓力，對農(nóng)業(yè)綠色發(fā)展的負(fù)面影響較強，而黃河上游各省份工業(yè)化水平相對較低，對農(nóng)業(yè)綠色發(fā)展的影響不顯著。

４）城鎮(zhèn)化水平。城鎮(zhèn)化水平對農(nóng)業(yè)綠色發(fā)展的影響在黃河上游地區(qū)和下游地區(qū)相對較強、在中游地區(qū)相對較弱，其原因是：城鎮(zhèn)化水平的提升為黃河上、下游地區(qū)農(nóng)業(yè)高質(zhì)量發(fā)展在資金、技術(shù)、專業(yè)人才等方面注入了新的活力，但黃河中游地區(qū)人口數(shù)量較多、農(nóng)業(yè)基礎(chǔ)龐大，城鎮(zhèn)化發(fā)展可能造成農(nóng)村社會結(jié)構(gòu)的改變，在土地、水資源利用方面與農(nóng)業(yè)產(chǎn)生沖突等，因此與上、下游地區(qū)相比，中游地區(qū)城鎮(zhèn)化發(fā)展對農(nóng)業(yè)綠色發(fā)展的推動作用相對較小。

５）地形起伏度。地形起伏度對農(nóng)業(yè)綠色發(fā)展具有顯著的負(fù)面影響，影響強度呈現(xiàn)由西向東（從上游向下游）逐漸遞減的趨勢。黃河中上游地區(qū)的地形起伏度較大，農(nóng)業(yè)生產(chǎn)機械化難以大規(guī)模實施，且存在灌溉條件較差、水土流失嚴(yán)重等問題，對農(nóng)業(yè)綠色發(fā)展產(chǎn)生較大影響。

６）年降水量。年降水量對農(nóng)業(yè)綠色發(fā)展的影響在中上游地區(qū)的甘肅省和陜西省部分地區(qū)、山西省呂梁等地區(qū)較強，在其他地區(qū)較弱，原因是甘肅和陜西年降水量較少、水資源匱乏，影響農(nóng)作物的生長和灌溉用水，在長期干旱條件下土地鹽堿化加劇、肥力降低，使得農(nóng)田生態(tài)系統(tǒng)退化、災(zāi)害風(fēng)險增加，不利于農(nóng)業(yè)綠色可持續(xù)發(fā)展。

４ 結(jié)論與建議

４．１ 結(jié)論

１）２０１１—２０２１ 年黃河流域農(nóng)業(yè)綠色發(fā)展水平呈上升趨勢，但研究期末整體水平仍較低；從各區(qū)域來看，下游地區(qū)農(nóng)業(yè)綠色發(fā)展水平明顯高于中上游地區(qū)；從農(nóng)業(yè)綠色發(fā)展的３ 個維度看，發(fā)展高效維度在上游和下游地區(qū)上升趨勢明顯，生態(tài)友好和資源節(jié)約維度在研究期上升幅度較小。

２）黃河流域農(nóng)業(yè)綠色發(fā)展水平存在明顯的空間正相關(guān)性，高－高集聚和低－低集聚地級市相對較多，無高－低集聚和低－高集聚情況。

３）農(nóng)業(yè)綠色發(fā)展各影響因素的空間異質(zhì)性顯著且影響強度不同，其中產(chǎn)業(yè)結(jié)構(gòu)、種植結(jié)構(gòu)、工業(yè)化水平、地形起伏度的影響空間差異較小，城鎮(zhèn)化水平和年降水量的影響空間差異較大，產(chǎn)業(yè)結(jié)構(gòu)、工業(yè)化水平和地形起伏度有負(fù)面影響，種植結(jié)構(gòu)、城鎮(zhèn)化水平和年降水量有正面影響，影響強度大小為地形起伏度＞工業(yè)化水平＞城鎮(zhèn)化水平＞產(chǎn)業(yè)結(jié)構(gòu)＞種植結(jié)構(gòu)＞年降水量。

４．２ 建議

１）轉(zhuǎn)變農(nóng)業(yè)發(fā)展理念，優(yōu)化農(nóng)業(yè)生產(chǎn)模式。一是將農(nóng)業(yè)與綠色發(fā)展有效結(jié)合，重視環(huán)境質(zhì)量、生態(tài)保護和農(nóng)業(yè)資源的保護，確保農(nóng)業(yè)發(fā)展與生態(tài)環(huán)境的和諧共生；二是改變農(nóng)業(yè)生產(chǎn)經(jīng)營模式，建立以綠色技術(shù)為核心的農(nóng)業(yè)生產(chǎn)體系，推動綠色技術(shù)在農(nóng)業(yè)生產(chǎn)、加工和流通等各個環(huán)節(jié)的應(yīng)用；三是調(diào)整農(nóng)業(yè)種植結(jié)構(gòu)，加強綠色農(nóng)產(chǎn)品供應(yīng)與產(chǎn)業(yè)鏈建設(shè)，提高農(nóng)業(yè)的質(zhì)量、效益和市場競爭力，進一步激發(fā)農(nóng)業(yè)綠色發(fā)展的內(nèi)在驅(qū)動力。

２）加強區(qū)域協(xié)作，形成黃河流域農(nóng)業(yè)綠色發(fā)展合力。發(fā)揮黃河流域各地區(qū)的農(nóng)業(yè)綠色發(fā)展優(yōu)勢，加強區(qū)域間協(xié)同合作。中、下游地區(qū)經(jīng)濟發(fā)展和生產(chǎn)技術(shù)水平較高，應(yīng)發(fā)揮輻射帶動作用，向上游地區(qū)輸送先進農(nóng)業(yè)管理理念和農(nóng)業(yè)生產(chǎn)技術(shù)；上游地區(qū)可利用自身優(yōu)勢，為中下游地區(qū)提供優(yōu)質(zhì)的初級農(nóng)產(chǎn)品等，拓展黃河流域農(nóng)業(yè)綠色發(fā)展的產(chǎn)業(yè)鏈。

３）因地制宜，實施差異化農(nóng)業(yè)綠色發(fā)展策略。黃河上游地區(qū)可兼顧生態(tài)保護與農(nóng)業(yè)效益提升，通過加強技術(shù)研發(fā)與培訓(xùn)，開發(fā)更多適宜的農(nóng)作物品種；中游地區(qū)面臨資源環(huán)境的約束，應(yīng)當(dāng)注重農(nóng)業(yè)資源高效節(jié)約利用，大力發(fā)展生態(tài)友好與資源節(jié)約型農(nóng)業(yè)；下游地區(qū)應(yīng)持續(xù)推進農(nóng)業(yè)節(jié)水、減肥、減藥，加快環(huán)境治理和生態(tài)修復(fù)，降低農(nóng)業(yè)面源污染，不斷提升農(nóng)業(yè)綠色發(fā)展水平。

參考文獻：

［１］ 金書秦，沈貴銀．中國農(nóng)業(yè)面源污染的困境擺脫與綠色轉(zhuǎn)型［Ｊ］．改革，２０１３（５）：７９－８７．

［２］ 魏琦，張斌，金書秦．中國農(nóng)業(yè)綠色發(fā)展指數(shù)構(gòu)建及區(qū)域比較研究［Ｊ］．農(nóng)業(yè)經(jīng)濟問題，２０１８（１１）：１１－２０．

［３］ 辜勝阻，曹冬梅，韓龍艷．“十三五”中國城鎮(zhèn)化六大轉(zhuǎn)型與健康發(fā)展［Ｊ］．中國人口·資源與環(huán)境，２０１７，２７（４）：６－１５．

［４］ 黃群慧，魏后凱，劉元春，等．學(xué)習(xí)黨的二十大精神筆談［Ｊ］．財貿(mào)經(jīng)濟，２０２３，４４（２）：５－２３．

［５］ 李福奪，楊鵬，尹昌斌．我國農(nóng)業(yè)綠色發(fā)展的基本理論與研究展望［Ｊ］．中國農(nóng)業(yè)資源與區(qū)劃，２０２０，４１（１０）：１－７．

［６］ 趙會杰，于法穩(wěn)．基于熵值法的糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色發(fā)展水平評價［Ｊ］．改革，２０１９（１１）：１３６－１４６．

［７］ 蓋美，楊苘菲，何亞寧．東北糧食主產(chǎn)區(qū)農(nóng)業(yè)綠色發(fā)展水平時空演化及其影響因素［Ｊ］．資源科學(xué)，２０２２，４４（５）：９２７－９４２．

［８］ 王磊，馬金銘．?dāng)?shù)字普惠金融影響農(nóng)業(yè)綠色發(fā)展的機制與效應(yīng)［Ｊ］．華南農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），２０２３，２２（６）：１４－２７．

［９］ 樊勝岳，李耀龍，馬曉杰，等．?dāng)?shù)字化水平對農(nóng)業(yè)綠色發(fā)展影響的實證研究：基于中國３０ 個省份的面板數(shù)據(jù)［Ｊ］．世界農(nóng)業(yè)，２０２１（１２）：４－１６．

［１０］ 梁俊芬，方偉，萬忠，等．中國省域農(nóng)業(yè)科技創(chuàng)新能力評價：基于綠色發(fā)展視角［Ｊ］．科技管理研究，２０２０，４０（９）：６０－６７．

［１１］ 熊素，羅蓉．“雙碳”目標(biāo)下中國農(nóng)業(yè)綠色發(fā)展：理論框架、困境審視及破局之道［Ｊ］．農(nóng)村經(jīng)濟，２０２３（２）：１０６－１１５．

［１２］ 劉惠良，劉紅峰，尹少華．基于鄉(xiāng)村振興視域下農(nóng)業(yè)綠色發(fā)展水平測度研究［Ｊ／ ＯＬ］．經(jīng)濟地理（２０２２－１０－２８）［２０２５－０２－２６］．ｈｔｔｐｓ：／／ ｎａｖｉ．ｃｎｋｉ．ｎｅｔ／ ｋｎａｖｉ／ ．

［１３］ 余永琦，王長松，彭柳林，等．基于熵權(quán)ＴＯＰＳＩＳ 模型的農(nóng)業(yè)綠色發(fā)展水平評價與障礙因素分析：以江西省為例［Ｊ］．中國農(nóng)業(yè)資源與區(qū)劃，２０２２，４３（２）：１８７－１９６．

［１４］ 明翠琴．中國農(nóng)業(yè)綠色增長評價指標(biāo)體系的構(gòu)建及實證［Ｊ］．技術(shù)經(jīng)濟與管理研究，２０２１（９）：１０８－１１３．

［１５］ 金賽美．中國省際農(nóng)業(yè)綠色發(fā)展水平及區(qū)域差異評價［Ｊ］．求索，２０１９（２）：８９－９５．

［１６］ 喻保華，王肖楊，宋春曉，等．中國農(nóng)業(yè)綠色發(fā)展時空演化及耦合協(xié)調(diào)研究［Ｊ］．生態(tài)經(jīng)濟，２０２３，３９（５）：１３２－１３９．

［１７］ 段呈，余成群，李少偉，等．西藏高原縣域尺度農(nóng)業(yè)綠色發(fā)展水平綜合評價：以白朗縣為例［Ｊ］．中國農(nóng)業(yè)資源與區(qū)劃，２０２３，４４（５）：８２－９０．

［１８］ 鄢曹政，殷旅江，何波．物流業(yè)集聚、空間溢出效應(yīng)與農(nóng)業(yè)綠色全要素生產(chǎn)率：基于省域數(shù)據(jù)的實證分析［Ｊ］．中國流通經(jīng)濟，２０２２，３６（９）：３－１６．

［１９］ 馬國群，譚硯文．環(huán)境規(guī)制對農(nóng)業(yè)綠色全要素生產(chǎn)率的影響研究：基于面板門檻模型的分析［Ｊ］．農(nóng)業(yè)技術(shù)經(jīng)濟，２０２１（５）：７７－９２．

［２０］ 李谷成，李欠男．“兩型社會”試驗區(qū)的設(shè)立促進了農(nóng)業(yè)綠色發(fā)展嗎？：基于ＰＳＭ－ＤＩＤ 模型的實證［Ｊ］．農(nóng)林經(jīng)濟管理學(xué)報，２０２２，２１（２）：１２７－１３５．。

［２１］ 田彩紅，李琳，廖斌．農(nóng)村一二三產(chǎn)業(yè)融合能否促進農(nóng)業(yè)綠色發(fā)展？：以長江經(jīng)濟帶５７９ 個縣域為例［Ｊ］．自然資源學(xué)報，２０２４，３９（３）：６０１－６１９．

［２２］ 唐一帆，吳波．財政支農(nóng)促進了農(nóng)業(yè)綠色發(fā)展嗎？：基于ＰＶＡＲ 模型的實證檢驗［Ｊ］．湖南農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），２０２２，２３（６）：４６－５４．

［２３］ 臺德進，蔡榮．糧食價格支持政策與農(nóng)業(yè)綠色增長：促進或抑制［Ｊ］．貴州財經(jīng)大學(xué)學(xué)報，２０２４（２）：９１－１００．

［２４］ 郭海紅，李樹超．環(huán)境規(guī)制、空間效應(yīng)與農(nóng)業(yè)綠色發(fā)展［Ｊ］．研究與發(fā)展管理，２０２２，３４（２）：５４－６７．

［２５］ 查建平，周霞，周玉璽．黃河流域農(nóng)業(yè)綠色發(fā)展水平綜合評價分析［Ｊ］． 中國農(nóng)業(yè)資源與區(qū)劃，２０２２，４３（１）：１８－２８．

【責(zé)任編輯 張智民】

基金項目：國家社會科學(xué)基金資助項目（２３ＢＪＹ２４９）；河南省高校哲學(xué)社會科學(xué)創(chuàng)新團隊項目（２０２４－ＣＸＴＤ－０４）