海南島農(nóng)業(yè)綠色發(fā)展指標(biāo)時(shí)空變化特征*

宋晨陽, 張建杰, 劉 玲, 馬文奇, 馬 林, 丁 尚, 趙洪偉**

海南島農(nóng)業(yè)綠色發(fā)展指標(biāo)時(shí)空變化特征*

宋晨陽1,2, 張建杰3, 劉 玲4, 馬文奇3, 馬 林4, 丁 尚2, 趙洪偉2**

(1. 海南大學(xué)林學(xué)院 海口 570228; 2. 海南大學(xué)生態(tài)與環(huán)境學(xué)院/海南省農(nóng)林環(huán)境過程與生態(tài)調(diào)控重點(diǎn)實(shí)驗(yàn)室 海口 570228; 3. 河北農(nóng)業(yè)大學(xué)資源與環(huán)境科學(xué)學(xué)院 保定 071001; 4. 中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所農(nóng)業(yè)資源研究中心 石家莊 050022)

海南自建省以來農(nóng)業(yè)生產(chǎn)規(guī)模、結(jié)構(gòu)及生產(chǎn)方式均發(fā)生了很大改變, 影響了其農(nóng)業(yè)綠色發(fā)展程度。本文借助中國農(nóng)業(yè)綠色發(fā)展指標(biāo)體系, 利用NUFER模型(NUtrient flows in Food chains, Environment and Resources use), 選取并定量計(jì)算1988—2017年海南島20項(xiàng)農(nóng)業(yè)綠色發(fā)展指標(biāo), 研究其時(shí)空變化特征, 探究制約海南島農(nóng)業(yè)綠色發(fā)展的影響因素。研究表明: 30年間海南省農(nóng)業(yè)綠色發(fā)展水平總體呈下降趨勢(shì), 較低級(jí)別(Ⅲ、Ⅳ級(jí))指標(biāo)數(shù)量由12個(gè)增加至15個(gè), 較高級(jí)別(Ⅰ、Ⅱ級(jí))指標(biāo)數(shù)量由8個(gè)減少至5個(gè)。隨著蔬菜、水果等經(jīng)濟(jì)作物種植比例增加, 畜禽養(yǎng)殖規(guī)模擴(kuò)大與集約化程度提高, 海南農(nóng)業(yè)產(chǎn)值與社會(huì)發(fā)展水平不斷升高, 農(nóng)業(yè)產(chǎn)值和農(nóng)業(yè)機(jī)械化投入分別從0.62×104￥?hm–2和3 kW?hm–2增加至16.1×104￥?hm–2和13 kW?hm–2。30年來農(nóng)業(yè)資源(氮素、磷素、農(nóng)藥和農(nóng)膜)過量投入, 其中農(nóng)藥和農(nóng)膜的使用強(qiáng)度增幅最大, 分別從8.0 kg?hm–2和0 kg?hm–2增加至41.9 kg?hm–2和34.7 kg?hm–2; 由于生產(chǎn)力水平較低, 導(dǎo)致養(yǎng)分(農(nóng)田氮素盈余量、農(nóng)田氮素徑流量、農(nóng)業(yè)源氨揮發(fā)等)環(huán)境排放量大幅增加, 農(nóng)業(yè)源氨揮發(fā)和單位播種面積農(nóng)田氮素盈余增幅最大, 分別從61.0 kg?hm–2和152.1 kg?hm–2增加至131.4 kg?hm–2和297.9 kg?hm–2。農(nóng)業(yè)資源投入過量、農(nóng)牧分離導(dǎo)致的環(huán)境污染是制約海南島農(nóng)業(yè)綠色發(fā)展的主要原因。在空間上, 由于中部為山地林區(qū), 是國家重點(diǎn)生態(tài)功能保護(hù)區(qū), 開發(fā)力度較小, 因此沿海平原的資源投入增長更快、養(yǎng)分環(huán)境損失量更多, 經(jīng)濟(jì)增長、社會(huì)發(fā)展和農(nóng)業(yè)生產(chǎn)3類指標(biāo)提升更為迅速。今后海南島應(yīng)合理優(yōu)化沿海地區(qū)種植結(jié)構(gòu), 提高生產(chǎn)管理水平, 加強(qiáng)農(nóng)牧結(jié)合減少資源損失, 實(shí)現(xiàn)海南島農(nóng)業(yè)綠色發(fā)展。

農(nóng)業(yè)綠色發(fā)展; NUFER模型; 種植業(yè)結(jié)構(gòu)變化; 農(nóng)牧系統(tǒng); 熱帶農(nóng)業(yè)

中國農(nóng)業(yè)的快速發(fā)展取得了巨大成就, 滿足了人們?nèi)找嬖鲩L的糧食需求[1]。但農(nóng)業(yè)資源消耗過高, 帶來了土壤酸化[2]、水體富營養(yǎng)化[3]、地下水硝酸鹽含量超標(biāo)[4]和生物多樣性退化[5]等一系列環(huán)境問題。我國農(nóng)業(yè)發(fā)展迫切需要提高生產(chǎn)效率、減少環(huán)境浪費(fèi), 協(xié)調(diào)其與環(huán)境保護(hù)之間的關(guān)系[6]。為了保證糧食安全、保護(hù)生態(tài)環(huán)境, 我國政府提出了農(nóng)業(yè)綠色發(fā)展道路[7]。不同學(xué)者對(duì)農(nóng)業(yè)綠色發(fā)展做出了不同的解讀, 但對(duì)其科學(xué)內(nèi)涵的解讀尚未達(dá)成共識(shí)。李谷成[8]將農(nóng)業(yè)綠色發(fā)展定義為農(nóng)業(yè)投入與產(chǎn)品生產(chǎn)、加工、流通及其消費(fèi)全鏈條的綠色化。金書秦等[9]認(rèn)為農(nóng)業(yè)綠色發(fā)展有去污、提質(zhì)、增效3個(gè)階段, 是不以環(huán)境污染為代價(jià), 通過綠色投入(包括資源、科技和經(jīng)濟(jì)的投入)、綠色消費(fèi)帶動(dòng)經(jīng)濟(jì)增長, 形成具有內(nèi)在驅(qū)動(dòng)力的良性循環(huán)。Shen等[10]提出農(nóng)業(yè)綠色發(fā)展是協(xié)調(diào)“綠色”與“發(fā)展”的關(guān)系, 把當(dāng)前資源消耗高、環(huán)境成本高的生產(chǎn)方式轉(zhuǎn)變?yōu)楦呱a(chǎn)效率、低環(huán)境排放的綠色農(nóng)業(yè)總體目標(biāo)。

農(nóng)業(yè)綠色發(fā)展是未來農(nóng)業(yè)的發(fā)展方向, 實(shí)現(xiàn)綠色發(fā)展, 首先需要評(píng)價(jià)當(dāng)前農(nóng)業(yè)發(fā)展水平, 剖析制約農(nóng)業(yè)綠色發(fā)展的關(guān)鍵問題。聯(lián)合國2015年可持續(xù)發(fā)展議程通過了17項(xiàng)可持續(xù)發(fā)展指標(biāo), 其中8項(xiàng)與農(nóng)業(yè)密切相關(guān)[11]。歐盟提出的農(nóng)業(yè)環(huán)境指數(shù)應(yīng)用也較為廣泛, 該體系可權(quán)衡農(nóng)業(yè)、經(jīng)濟(jì)、政策與環(huán)境等因素之間的相互關(guān)系[12]。魏琦等[13]以資源節(jié)約、環(huán)境友好、生態(tài)保育和質(zhì)量高效為側(cè)重點(diǎn), 構(gòu)建我國農(nóng)業(yè)綠色發(fā)展指數(shù), 指出完善數(shù)據(jù)監(jiān)測(cè)將利于未來農(nóng)業(yè)評(píng)價(jià)體系的建立, 環(huán)境污染和質(zhì)量高效是我國農(nóng)業(yè)綠色發(fā)展的突出短板。袁久和[14]從經(jīng)濟(jì)、資源環(huán)境、農(nóng)村社會(huì)發(fā)展和政策支持等方面, 評(píng)價(jià)了我國農(nóng)業(yè)發(fā)展?fàn)顩r, 并指出加強(qiáng)監(jiān)測(cè), 完善政策支持是實(shí)現(xiàn)農(nóng)業(yè)綠色發(fā)展的有效舉措。楊世琦[15]基于歐氏距離法構(gòu)建了農(nóng)業(yè)可持續(xù)發(fā)展模型, 分析影響寧夏吳忠市可持續(xù)發(fā)展的主要原因是農(nóng)業(yè)資源利用不合理、保護(hù)性農(nóng)作技術(shù)推廣不夠、土壤養(yǎng)分含量較低且農(nóng)村人均收入偏低等。Liu等[16]利用層次分析法構(gòu)建評(píng)估模型, 評(píng)價(jià)了我國1978—2017年的農(nóng)業(yè)綠色生產(chǎn)水平, 指出我國農(nóng)業(yè)綠色生產(chǎn)水平增長趨勢(shì)明顯, 但西部內(nèi)陸地區(qū)與東部地區(qū)存在較大差距, 并提出加強(qiáng)數(shù)據(jù)監(jiān)測(cè), 構(gòu)建更合理完善的農(nóng)業(yè)綠色生產(chǎn)指標(biāo)體系將有助于我國農(nóng)業(yè)現(xiàn)代化建設(shè)。

這些研究通過不同方法、從不同尺度分析了農(nóng)業(yè)發(fā)展?fàn)顩r并提出了相應(yīng)的科學(xué)建議, 但更多的是對(duì)不同地區(qū)或不同環(huán)節(jié)發(fā)展水平的評(píng)價(jià), 缺乏從“資源-生產(chǎn)-食品加工-消費(fèi)”整個(gè)農(nóng)業(yè)生產(chǎn)和食物系統(tǒng)層面的綜合考量, 缺少對(duì)指標(biāo)間關(guān)聯(lián)的探究及其動(dòng)態(tài)響應(yīng)的定量分析。我國不同省份農(nóng)業(yè)發(fā)展特征不同, 海南島是我國唯一的“熱帶特色農(nóng)業(yè)基地”, 2018年國家確立了海南建設(shè)“全面深化改革開放試驗(yàn)區(qū)、國家生態(tài)文明試驗(yàn)區(qū)、國家重大戰(zhàn)略服務(wù)保障區(qū)、國際旅游消費(fèi)中心[17]”的戰(zhàn)略定位, 這決定了海南島未來不但要保證熱帶特色食品供應(yīng), 同時(shí)也要發(fā)展經(jīng)濟(jì)和保護(hù)生態(tài)環(huán)境。本文利用張建杰等[18]構(gòu)建的農(nóng)業(yè)綠色發(fā)展指標(biāo)體系, 研究海南島農(nóng)業(yè)綠色發(fā)展時(shí)空變化特征, 找出制約海南島農(nóng)業(yè)綠色發(fā)展的影響因素及限制因子, 以海南島為例探討大型熱帶島嶼農(nóng)業(yè)綠色發(fā)展的實(shí)現(xiàn)途徑。

1 研究區(qū)域概況與研究方法

1.1 研究區(qū)域概況

海南1988年建省, 共有4個(gè)地級(jí)市和15個(gè)縣區(qū), 全省面積3.54萬km2。本研究以海南島為研究對(duì)象, 島嶼面積3.39萬km2, 地勢(shì)中間高四周低, 以五指山、鸚哥嶺為核心向外由山地、丘陵、臺(tái)地、平原構(gòu)成環(huán)狀圈層地貌; 屬于熱帶海洋性季風(fēng)氣候, 年均氣溫22~26 ℃, 年光照達(dá)1 750~2 650 h, 年均降水量1 639 mm。2017年耕地面積43.9萬hm2, 典型土壤為赤紅壤, 復(fù)種指數(shù)200%左右。截至2017年總?cè)丝?10.4萬, 地區(qū)生產(chǎn)總值4 462.5億元, 2017年農(nóng)林牧漁產(chǎn)值占地區(qū)生產(chǎn)總值比為33.4%[19]。

海南島不同區(qū)域農(nóng)業(yè)功能劃分不同(圖1): 山區(qū)林業(yè)經(jīng)濟(jì)區(qū)發(fā)展橡膠()、特色水果和南藥等經(jīng)濟(jì)作物及特色養(yǎng)殖, 丘陵經(jīng)作養(yǎng)殖區(qū)發(fā)展經(jīng)濟(jì)作物、畜禽養(yǎng)殖綜合農(nóng)業(yè), 平原種養(yǎng)農(nóng)業(yè)區(qū)發(fā)展高效種養(yǎng)產(chǎn)業(yè), 另外中部山區(qū)熱帶雨林生態(tài)功能區(qū)是國家重點(diǎn)生態(tài)功能區(qū), 開發(fā)受限。海南省是正在創(chuàng)建的國家農(nóng)業(yè)綠色發(fā)展先行區(qū)(整省創(chuàng)建), 并于三亞市、樂東縣、陵水縣建立了集科研、生產(chǎn)、銷售、科技交流、成果轉(zhuǎn)化為一體的“南繁育種基地”。1988—2017年海南島總播種面積增加了0.6倍, 水稻()播種面積減少了22%, 熱帶水果、蔬菜和特色經(jīng)濟(jì)作物[橡膠、椰子()、檳榔()、胡椒()等]的播種面積分別增加了2.7倍、4.8倍和2.0倍, 截至2017年海南島水果、蔬菜經(jīng)濟(jì)作物總產(chǎn)量增加了8.3倍, 達(dá)756萬t; 畜禽養(yǎng)殖總數(shù)增加了5.8倍, 2017年達(dá)325萬LU(折500 kg奶牛)[19]。

圖1 海南島功能分區(qū)示意圖

1.2 數(shù)據(jù)來源

研究所需數(shù)據(jù)主要包括統(tǒng)計(jì)數(shù)據(jù)、文獻(xiàn)資料及模型原有參數(shù)[20]。統(tǒng)計(jì)數(shù)據(jù)是模型分析的基礎(chǔ), 而文獻(xiàn)資料主要為校正模型參數(shù)提供依據(jù)。

1)統(tǒng)計(jì)數(shù)據(jù): 海南島歷年農(nóng)業(yè)生產(chǎn)總值、人口數(shù)量、農(nóng)作物播種面積及產(chǎn)量、化肥、農(nóng)藥、農(nóng)膜投入量、畜禽養(yǎng)殖數(shù)量及產(chǎn)量、城市與農(nóng)村食物人均消費(fèi)水平來自于海南統(tǒng)計(jì)年鑒[19]; 由于統(tǒng)計(jì)數(shù)據(jù)中無年末蛋禽存欄量, 故通過禽蛋產(chǎn)量與FAO中國歷年蛋禽平均產(chǎn)蛋量反推蛋禽養(yǎng)殖數(shù)量。

2)文獻(xiàn)數(shù)據(jù): 1988—2017年海南島種植業(yè)不同作物養(yǎng)分含量、灌溉水氮素量[21]及氨揮發(fā)系數(shù)[22]來自文獻(xiàn), 秸稈、畜禽糞污回收利用情況來自政府報(bào)告[23-24], 生物固氮、氮沉降、徑流排放系數(shù)、動(dòng)物糞便氨揮發(fā)、N2O揮發(fā)和反硝化排放系數(shù)參考NUFER模型參數(shù)[20]。

1.3 研究方法

中國農(nóng)業(yè)綠色發(fā)展指標(biāo)體系是張建杰等[18]利用NUFER模型[20]構(gòu)建的可定量計(jì)算中國及其各省農(nóng)業(yè)發(fā)展水平的指標(biāo)體系。該指標(biāo)體系從農(nóng)業(yè)生產(chǎn)、資源與環(huán)境3方面出發(fā), 研究農(nóng)產(chǎn)品“生產(chǎn)-加工-消費(fèi)”全鏈條, 經(jīng)濟(jì)增長、社會(huì)發(fā)展程度、資源投入情況、農(nóng)業(yè)生產(chǎn)水平以及生態(tài)環(huán)境排放損失量5類指標(biāo)的優(yōu)劣程度。每個(gè)類型分為直接指標(biāo)(direct)、間接指標(biāo)(indirect)和自選指標(biāo)(optional)3類, 每項(xiàng)指標(biāo)均分為Ⅰ~Ⅳ4個(gè)等級(jí), 從Ⅰ到Ⅳ依次表示農(nóng)業(yè)綠色發(fā)展程度由好到差。本研究依據(jù)系統(tǒng)性、完整性、代表性和重點(diǎn)性原則(即具體地從中國農(nóng)業(yè)綠色發(fā)展指標(biāo)體系中選出可有效解釋海南島農(nóng)業(yè)特征、突出其生產(chǎn)特色的指標(biāo)), 選取20項(xiàng)指標(biāo)定量研究海南島農(nóng)業(yè)綠色發(fā)展的時(shí)空變化特征。

經(jīng)濟(jì)增長類指標(biāo)(economy growth): 包括單位面積農(nóng)業(yè)產(chǎn)值、人均GDP、農(nóng)業(yè)產(chǎn)值占GDP比例3個(gè)指標(biāo), 反映海南島農(nóng)業(yè)經(jīng)濟(jì)效益、農(nóng)業(yè)對(duì)經(jīng)濟(jì)的貢獻(xiàn)以及農(nóng)業(yè)發(fā)展對(duì)人均生活水平的影響。

社會(huì)發(fā)展類指標(biāo)(social development): 包括人均耕地占有水平、農(nóng)田灌溉水平、農(nóng)業(yè)機(jī)械化水平3個(gè)指標(biāo), 表現(xiàn)海南島農(nóng)業(yè)生產(chǎn)條件。

資源投入類指標(biāo)(resources input): 包括氮素、磷素以及農(nóng)藥、農(nóng)膜使用強(qiáng)度4個(gè)指標(biāo), 衡量海南島當(dāng)前農(nóng)業(yè)生產(chǎn)要素投入水平。

生態(tài)環(huán)境類指標(biāo)(ecological environment): 包括農(nóng)田氮素盈余及徑流、農(nóng)業(yè)系統(tǒng)氨揮發(fā)、農(nóng)業(yè)源溫室氣體(GHG)排放強(qiáng)度、畜禽承載量5個(gè)指標(biāo), 反映海南島種養(yǎng)生產(chǎn)對(duì)水、土、氣的環(huán)境影響。

農(nóng)業(yè)生產(chǎn)類指標(biāo)(agricultural production): 包括農(nóng)田系統(tǒng)氮素利用效率、畜牧系統(tǒng)氮素利用效率、蔬菜單產(chǎn)水平、水果單產(chǎn)水平以及單位動(dòng)物蛋白產(chǎn)量5個(gè)指標(biāo), 反映海南島農(nóng)業(yè)生產(chǎn)效率及主要作物生產(chǎn)水平。

各指標(biāo)計(jì)算方法與分級(jí)標(biāo)準(zhǔn)與中國農(nóng)業(yè)綠色發(fā)展指標(biāo)體系[18]相同。

使用EXCEL 2010繪制1988—2017年海南島農(nóng)業(yè)綠色發(fā)展所有指標(biāo)結(jié)果的時(shí)間變化圖。使用ArcGIS 10.2.2繪制1988年、2017年重要指標(biāo)結(jié)果的空間分布圖。

2 結(jié)果與分析

2.1 農(nóng)業(yè)綠色發(fā)展指標(biāo)變化情況

海南島1988年和2017年農(nóng)業(yè)綠色發(fā)展指標(biāo)變化情況如表1所示。30年間綠色發(fā)展水平整體呈下降趨勢(shì), 較低級(jí)別(Ⅲ、Ⅳ級(jí))指標(biāo)數(shù)量由12個(gè)增加至15個(gè), 較高級(jí)別(Ⅰ、Ⅱ級(jí))指標(biāo)數(shù)量由8個(gè)減少至5個(gè)。經(jīng)濟(jì)增長指標(biāo)中單位面積農(nóng)業(yè)產(chǎn)值由Ⅳ級(jí)上升至Ⅱ級(jí)水平, 農(nóng)業(yè)產(chǎn)值占GDP比例從Ⅰ級(jí)水平降低至Ⅱ級(jí)水平, 人均GDP由Ⅳ級(jí)水平變?yōu)棰蚣?jí)水平。社會(huì)發(fā)展指標(biāo)中人均耕地占有水平與農(nóng)田灌溉水平指標(biāo)分別處于Ⅳ、Ⅱ級(jí)別不變, 僅農(nóng)業(yè)機(jī)械化水平由Ⅲ級(jí)提高至Ⅰ級(jí)水平。資源投入所有指標(biāo)均降低至Ⅳ級(jí)水平。生態(tài)環(huán)境指標(biāo)中, 農(nóng)田氮素盈余及畜禽承載量指標(biāo)由Ⅱ級(jí)下降至Ⅳ級(jí)水平, 農(nóng)田氮素徑流指標(biāo)由Ⅰ級(jí)水平降低至Ⅲ級(jí)水平, 農(nóng)業(yè)源GHG排放強(qiáng)度指標(biāo)30年間均處于Ⅳ級(jí)水平, 農(nóng)業(yè)系統(tǒng)氨揮發(fā)指標(biāo)由Ⅲ級(jí)降低至Ⅳ級(jí)水平。農(nóng)業(yè)生產(chǎn)指標(biāo)中水果單產(chǎn)水平、蔬菜單產(chǎn)水平和單位動(dòng)物蛋白產(chǎn)量自Ⅳ級(jí)提高至Ⅲ級(jí)水平, 氮素利用效率始終處于Ⅳ級(jí)水平。

2.2 經(jīng)濟(jì)增長指標(biāo)時(shí)空變化特征

1988—2017年海南島經(jīng)濟(jì)快速增長, 農(nóng)業(yè)對(duì)經(jīng)濟(jì)產(chǎn)值的貢獻(xiàn)顯著降低(圖2a), 單位耕地面積農(nóng)業(yè)產(chǎn)值由0.62×104￥?hm–2增加到16.1×104￥?hm–2, 提高了25.0倍; 農(nóng)業(yè)產(chǎn)值占GDP比重由74.8%降低至33.4%, 人均GDP從0.12×104￥?cap.–1增長至4.8×104￥?cap.–1。單位耕地面積農(nóng)業(yè)產(chǎn)值普遍性與代表性最強(qiáng), 因此繪制其1988年、2017年的空間分布圖(圖2b)。1988年海南島單位耕地面積農(nóng)業(yè)產(chǎn)值較低, 所有地區(qū)均處Ⅲ、Ⅳ級(jí)水平, 其中三亞市最高(1.9×104￥?hm–2), 東方市最低(0.8×104￥?hm–2); 由于南部沿海丘陵經(jīng)作養(yǎng)殖區(qū)和部分南部平原種養(yǎng)農(nóng)業(yè)區(qū)經(jīng)濟(jì)作物產(chǎn)量較高, 這些地區(qū)單位面積農(nóng)業(yè)經(jīng)濟(jì)產(chǎn)值顯著高于其他地區(qū)。2017年三亞、瓊海最高, 分別為44.4×104￥?hm–2和32.3×104￥?hm–2, 臨高縣最低, 為7.5×104￥?hm–2。

表1 1988年和2017年海南島農(nóng)業(yè)綠色發(fā)展指標(biāo)變化

圖2 1988—2017年海南島經(jīng)濟(jì)增長指標(biāo)時(shí)空變化特征(a:經(jīng)濟(jì)增長指標(biāo)時(shí)間變化趨勢(shì); b: 1988年和2017年單位面積農(nóng)業(yè)產(chǎn)值空間差異)

2.3 社會(huì)發(fā)展指標(biāo)時(shí)空變化特征

海南島社會(huì)發(fā)展指標(biāo)1988—2017年時(shí)間變化結(jié)果如圖3a所示, 所有指標(biāo)均呈上升趨勢(shì)。農(nóng)田灌溉水平從31%提高到44%; 人均耕地占有水平從0.08 hm2×cap.-1提高到0.11 hm2×cap.-1, 增加38%; 農(nóng)業(yè)機(jī)械化水平從3 kW×hm-2增長至13 kW×hm-2, 增長了3.3倍。我國普遍耕地不足, 海南島農(nóng)業(yè)機(jī)械化水平30年間提升較高, 但西北地區(qū)有熱帶沙漠氣候的傾向, 因水利灌溉網(wǎng)絡(luò)覆蓋不全、溝渠堵塞等原因?qū)r(nóng)作影響較大, 因此本文關(guān)注海南島農(nóng)田灌溉水平的空間變化特征(圖3b), 1988年海南島農(nóng)田灌溉水平僅保亭縣、五指山市達(dá)到Ⅰ級(jí)水平, 約1/3地區(qū)處于Ⅲ級(jí)水平; 2017年海南島農(nóng)田灌溉水平有所提高, 2/5地區(qū)高于Ⅰ級(jí)水平, 主要集中在東南沿海丘陵經(jīng)作養(yǎng)殖區(qū)和平原種養(yǎng)農(nóng)業(yè)區(qū), 與經(jīng)濟(jì)類指標(biāo)地區(qū)分布較為一致。

圖3 1988—2017年海南島社會(huì)發(fā)展指標(biāo)時(shí)空變化特征(a:社會(huì)發(fā)展指標(biāo)時(shí)間變化趨勢(shì); b: 1988年和2017年農(nóng)業(yè)灌溉水平空間差異)

2.4 資源投入指標(biāo)時(shí)空變化特征

海南島資源投入類指標(biāo)綠色化水平下降明顯, 30年間農(nóng)業(yè)資源投入大幅增加, 2017年所有指標(biāo)均下降至Ⅳ級(jí)水平(圖4a、4c)。氮素使用強(qiáng)度由217.0 kg×hm-2增加至411.7 kg×hm-2, 增加90%; 磷素使用強(qiáng)度由13.9 kg×hm-2增加至66.2 kg×hm-2, 增加了3.8倍; 農(nóng)藥使用強(qiáng)度由8.19 kg×hm-2增加至41.9 kg×hm-2, 增加了4.1倍。農(nóng)膜自2002年投入使用, 使用強(qiáng)度由4.1 kg×hm-2增長至2017年的34.7 kg×hm-2, 增加了7.5倍。

海南島是國家農(nóng)業(yè)綠色發(fā)展先行區(qū), 而養(yǎng)分的過量投入與這一發(fā)展目標(biāo)相悖。由圖4b、4d可知, 1988年各地區(qū)資源投入較少, 氮素使用強(qiáng)度除海口、澄邁外均處于Ⅰ、Ⅱ級(jí)水平(<225 kg×hm-2),澄邁處于Ⅲ級(jí)水平(354.6 kg×hm-2), ?？谑刑幱冖艏?jí)水平(>405 kg×hm-2)(?？谑?988年與2017年地區(qū)劃分不同, 1988年農(nóng)業(yè)資源使用量過高是由于土地利用不同); 1988年各地區(qū)磷素使用強(qiáng)度均處于Ⅰ級(jí)水平(<35 kg×hm-2)。2017年平原及丘陵地區(qū)氮素投入明顯增多, 大部分地區(qū)處于Ⅲ、Ⅳ級(jí)水平, 其中瓊海和三亞地區(qū)的氮素使用強(qiáng)度最高, 分別為818.9 kg×hm-2、625.6 kg×hm-2, 僅五指山市(178.1 kg×hm-2)處于Ⅰ級(jí)水平。2017年磷素使用強(qiáng)度三亞、瓊海較高, 分別為99.6 kg×hm-2、120.2 kg×hm-2, 而五指山市使用強(qiáng)度最低, 為18.9 kg×hm-2。由于瓜果、經(jīng)濟(jì)作物等集中種植在沿海地區(qū), 而中部山區(qū)為熱帶雨林生態(tài)功能區(qū), 開發(fā)受限, 因此沿海地區(qū)相較于中部山地養(yǎng)分資源投入增加更多。

2.5 生態(tài)環(huán)境指標(biāo)時(shí)空變化特征

1988—2017年海南島生態(tài)環(huán)境類指標(biāo)綠色化水平下降趨勢(shì)明顯, 養(yǎng)分環(huán)境排放量及畜禽承載量均有增加(圖5a、5c)。農(nóng)田氮素盈余和氮素徑流分別由152.1 kg×hm-2和6.5 kg×hm-2增加至297.9 kg×hm-2和12.4 kg×hm-2, 增加近1倍。農(nóng)業(yè)系統(tǒng)氨揮發(fā)由61.0 kg×hm-2增加至131.4 kg×hm-2, 增加了1.2倍。農(nóng)業(yè)源GHG排放強(qiáng)度隨水稻播種面積減少而降低, 由5 327.3 kg(CO2eq)×hm-2降低至5 061.3kg(CO2eq)×hm-2, 減少5%。畜禽承載量由3.4LU×hm-2增加至7.4LU×hm-2。

農(nóng)田氮素盈余量可定量種植過程養(yǎng)分輸入與輸出之差, 而農(nóng)業(yè)系統(tǒng)氨揮發(fā)則可指示化肥及畜禽糞污的損失情況。圖5b表明, 養(yǎng)分環(huán)境排放地區(qū)分布與資源投入趨勢(shì)基本相同, 沿海地區(qū)較中部地區(qū)增量更高, 損失也更大。1988年農(nóng)田氮素盈余除海口市較高(249.2 kg×hm-2), 其余地區(qū)處于Ⅰ、Ⅱ級(jí)水平, 最低為保亭縣、白沙縣(75 kg×hm-2)。1988年農(nóng)業(yè)系統(tǒng)氨揮發(fā)同樣除海口市其余地區(qū)處于Ⅰ、Ⅱ級(jí)水平(圖5d), 其中五指山市、瓊中縣、保亭縣最低(15 kg×hm-2)。2017年農(nóng)田氮素盈余僅屯昌縣、瓊中市、五指山市處于Ⅱ級(jí)水平, 其中五指山市最低(87.8 kg×hm-2), 沿海臺(tái)地農(nóng)田氮素盈余較高, 瓊海市最高(649.9 kg×hm-2)。2017年農(nóng)業(yè)系統(tǒng)氨揮發(fā)全島降低至Ⅲ、Ⅳ級(jí)水平, 其中瓊海市(264.6 kg×hm-2)高于其他地區(qū)(<200.0 kg×hm-2), 五指山市最低, 為65.3 kg×hm-2, 與資源利用類指標(biāo)地區(qū)分布特征較為一致。

圖4 1988—2017年海南島資源投入指標(biāo)時(shí)空變化特征(a、c: 資源利用指標(biāo)時(shí)間變化趨勢(shì); b: 1988年和2017年氮素使用強(qiáng)度空間差異; d: 1988年和2017年磷素使用強(qiáng)度空間差異)

圖5 1988—2017年海南島生態(tài)環(huán)境指標(biāo)時(shí)空變化特征(a、c:生態(tài)環(huán)境指標(biāo)時(shí)間變化趨勢(shì); b: 1988年和2017年農(nóng)田氮素盈余量空間差異; d: 1988年和2017年農(nóng)業(yè)系統(tǒng)氨揮發(fā)空間差異)

2.6 農(nóng)業(yè)生產(chǎn)指標(biāo)時(shí)空變化特征

海南島農(nóng)業(yè)生產(chǎn)指標(biāo)綠色化水平整體呈上升趨勢(shì)(圖6a、6c), 但仍處于較低水平。1998年以前各指標(biāo)均處于Ⅰ級(jí)水平。蔬菜單產(chǎn)、水果單產(chǎn)、單位動(dòng)物蛋白產(chǎn)量由10.0 t×hm-2、9.7 t×hm-2、11.0 kg×LU-1提高至2017年的22.9 t×hm-2、20.4 t×hm-2、20.6 kg×LU-1, 分別增加了1.3倍、1.1倍和87%。農(nóng)田系統(tǒng)氮素利用效率、畜禽系統(tǒng)氮素利用效率分別由15.9%、4.0%提高至2017年的20.0%、8.3%, 分別增加了30%和1.1倍。

從空間分布特征來看(圖6b), 各縣域種養(yǎng)產(chǎn)業(yè)單產(chǎn)不高。2017年蔬菜單產(chǎn)除三亞市、瓊海市處于Ⅲ級(jí)水平(35 t×hm-2), 其余地區(qū)均處于Ⅳ級(jí), 瓊中縣最低(9.7 t×hm-2); 2017年大部分市(縣)水果單產(chǎn)處于Ⅲ、Ⅳ級(jí)水平, 僅澄邁縣、臨高縣、昌江縣高達(dá)Ⅰ級(jí)水平(>29 t×hm-2), 屯昌縣最低(10.4 t×hm-2)。海南島單位動(dòng)物蛋白產(chǎn)量30年來一直處于Ⅲ、Ⅳ級(jí)水平(圖6d), 1988年中西部較低, 東部沿海地區(qū)較高, 其中五指山市最低(12.4 kg×LU-1), ?？谧罡?23.1 kg×LU-1); 2017年海南島單位動(dòng)物蛋白產(chǎn)量空間分布上呈現(xiàn)南低北高的特點(diǎn), 其中陵水縣最低(15.2 kg×LU-1), 文昌市最高(28.8 kg×LU-1)。

由此可見, 海南島農(nóng)業(yè)發(fā)展具有經(jīng)濟(jì)增長快、資源投入高、環(huán)境損失大、耕地資源緊缺、生產(chǎn)力水平低的特征; 而在空間上, 表現(xiàn)出南部地區(qū)經(jīng)濟(jì)與社會(huì)發(fā)展較好、中部山區(qū)熱帶雨林生態(tài)功能區(qū)資源投入及環(huán)境損失優(yōu)于平原種養(yǎng)農(nóng)業(yè)區(qū)及丘陵經(jīng)作養(yǎng)殖區(qū)、生產(chǎn)效率全島普遍不高的特征。

2.7 農(nóng)業(yè)綠色發(fā)展指標(biāo)間關(guān)系

海南島農(nóng)業(yè)綠色發(fā)展的經(jīng)濟(jì)增長、社會(huì)發(fā)展、資源投入、生態(tài)環(huán)境和農(nóng)業(yè)生產(chǎn)5類指標(biāo)所處等級(jí)反映其優(yōu)劣程度, 指標(biāo)間相互影響, 關(guān)系復(fù)雜。

通過對(duì)以上5類指標(biāo)的變化結(jié)果分析可以看出, 經(jīng)濟(jì)增長、資源投入和生態(tài)環(huán)境3類指標(biāo)均增加幅度最大、趨勢(shì)最明顯, 經(jīng)濟(jì)增長類指標(biāo)最能代表本地區(qū)的農(nóng)業(yè)發(fā)展?fàn)顩r, 資源投入類指標(biāo)最能反映農(nóng)業(yè)生產(chǎn)具體投入能耗的多少, 環(huán)境生產(chǎn)類指標(biāo)能表明農(nóng)業(yè)生產(chǎn)對(duì)環(huán)境的具體影響。從這3類指標(biāo)中選取最具有代表性的指標(biāo), 對(duì)其相關(guān)關(guān)系進(jìn)行分析, 得出畜禽承載量的提高明顯推動(dòng)了海南省農(nóng)業(yè)經(jīng)濟(jì)產(chǎn)值的提升, 其結(jié)果的相關(guān)性較強(qiáng)(2=0.97); 而蔬菜、水果和經(jīng)濟(jì)作物也是海南省農(nóng)業(yè)的主要產(chǎn)業(yè), 農(nóng)業(yè)經(jīng)濟(jì)產(chǎn)值隨蔬菜、水果和經(jīng)濟(jì)作物播種比例的增加大幅增長(2=0.79)。種養(yǎng)規(guī)模的擴(kuò)大勢(shì)必影響農(nóng)業(yè)資源投入及其環(huán)境損失。1988—2017年海南島熱帶水果、蔬菜和特色經(jīng)濟(jì)作物種植比例從31.2%增高至68.9%[19], 增加了1.2倍(圖7a), 與此同時(shí), 畜禽承載量由3.4 LU×hm-2增加至7.4 LU×hm-2, 增加了1.2倍(圖7b)。分析得出, 蔬菜、水果和經(jīng)濟(jì)作物播種比例的增大導(dǎo)致氮肥、農(nóng)藥使用強(qiáng)度的明顯提高以及農(nóng)田氮素盈余的增多, 該比例的增大與氮肥施用強(qiáng)度和農(nóng)田氮素盈余的增多關(guān)系最為顯著, 其2分別達(dá)0.93和0.91; 氮肥使用強(qiáng)度、農(nóng)田氮素盈余和農(nóng)業(yè)系統(tǒng)氨揮發(fā)隨畜禽承載量的增大而提高, 其中單位農(nóng)田氮素盈余及農(nóng)業(yè)系統(tǒng)氨揮發(fā)與畜禽承載量之間關(guān)系最為顯著,2分別為0.81和0.89(圖7d)。

圖6 1988—2017年海南島農(nóng)業(yè)生產(chǎn)指標(biāo)時(shí)空變化特征(a、c:農(nóng)業(yè)生產(chǎn)指標(biāo)時(shí)間變化趨勢(shì); b: 2017年單位面積蔬菜產(chǎn)量和水果產(chǎn)量空間差異; d: 1988年和2017年單位動(dòng)物蛋白產(chǎn)量空間差異)

圖7 1988—2017年海南島農(nóng)業(yè)綠色發(fā)展指標(biāo)間的關(guān)系(a:單位面積農(nóng)業(yè)產(chǎn)值與播種面積占比之間的關(guān)系; b:單位面積農(nóng)業(yè)產(chǎn)值與畜禽承載量之間的關(guān)系; c:資源投入及生態(tài)環(huán)境與播種面積占比之間的關(guān)系; d:資源投入及生態(tài)環(huán)境與畜禽承載量之間的關(guān)系)

3 討論

本研究首次定量分析了海南農(nóng)業(yè)綠色發(fā)展指標(biāo)的時(shí)空變化特征, 分析5類指標(biāo)之間的關(guān)系, 探究制約海南農(nóng)業(yè)綠色發(fā)展的關(guān)鍵問題, 得出海南島農(nóng)業(yè)發(fā)展具有以下特征:

1)農(nóng)業(yè)經(jīng)濟(jì)增長和社會(huì)發(fā)展變化趨勢(shì)一致, 1988—2017年各指標(biāo)提升明顯, 東南沿海提升速度更快。隨著社會(huì)的發(fā)展、經(jīng)濟(jì)的轉(zhuǎn)型, 農(nóng)田灌溉設(shè)備得以配備, 農(nóng)田灌溉水平隨之提高。從空間上看, 東部和南部沿海地區(qū)有效灌溉面積占比較大, 經(jīng)濟(jì)作物產(chǎn)量更高, 農(nóng)業(yè)產(chǎn)值提高更快, 說明農(nóng)業(yè)社會(huì)發(fā)展水平的提高可以促進(jìn)農(nóng)業(yè)經(jīng)濟(jì)產(chǎn)值的提升。

2)資源投入類與生態(tài)環(huán)境類指標(biāo)表現(xiàn)為資源投入量與環(huán)境排放量大幅增多的趨勢(shì), 且沿海(如三亞、瓊海等地)較中部山地更為嚴(yán)重。其原因是由于海南熱帶果蔬、特色經(jīng)濟(jì)作物種植比例和畜禽養(yǎng)殖規(guī)模不斷擴(kuò)大導(dǎo)致的養(yǎng)分過量投入與排放增加。30年間海南果蔬、熱帶特色經(jīng)濟(jì)作物等高收益[25]、高耗肥[26-28]作物逐漸發(fā)展為種植業(yè)的主導(dǎo)產(chǎn)業(yè), 其種植比例的增加既促進(jìn)了海南島農(nóng)業(yè)經(jīng)濟(jì)的快速發(fā)展,也導(dǎo)致肥料及農(nóng)藥投入量的持續(xù)增加, 帶來了高環(huán)境排放風(fēng)險(xiǎn)。海南島農(nóng)戶生產(chǎn)1季香蕉(Lour)氮肥投入強(qiáng)度為1 300 kg×hm-2, 種植豆角()則分5次施用750 kg×hm-2氮肥, 均高于作物推薦施氮水平(150~180 kg×hm-2)[29]。同時(shí), 蔬菜水果養(yǎng)分利用效率顯著低于其他作物, 過量的氮投入會(huì)導(dǎo)致養(yǎng)分損失、土壤盈余過量、耕地退化[5], 又間接導(dǎo)致生產(chǎn)力水平降低。2017年海南島約80%市(縣)農(nóng)田氮素盈余超過240 kg×hm-2, 與丁尚等[30]研究結(jié)果較為一致, 遠(yuǎn)超荷蘭養(yǎng)分監(jiān)管系統(tǒng)(MINAS)[31]針對(duì)不同農(nóng)田核算的氮素盈余(砂質(zhì)土壤60 kg×hm-2、黏質(zhì)土壤100 kg×hm-2)。已有研究表明, 20世紀(jì)80年代以來海南島過量氮素投入和種植結(jié)構(gòu)改變使得土壤酸化進(jìn)程加快[32-33]。另外, 除了養(yǎng)分過量排放, 海南島也面臨農(nóng)藥使用超標(biāo)的問題。由于環(huán)境特殊, 農(nóng)戶為保證收益, 農(nóng)藥使用頻繁, 甚至造成食品安全問題[34]。2009年海南島農(nóng)藥使用強(qiáng)度高達(dá)56.4 kg×hm-2, 是全國平均水平(16.3 kg×hm-2)[35]的3.5倍, 2010年受“毒豇豆”事件影響, 政府出臺(tái)了蔬菜安全監(jiān)督措施, 減少高毒農(nóng)藥的使用, 農(nóng)藥使用強(qiáng)度隨之減少。同樣, 受人口增長、飲食結(jié)構(gòu)變化等因素影響, 海南島畜禽規(guī)模也在不斷擴(kuò)大, 畜禽系統(tǒng)氮素利用效率和單位動(dòng)物蛋白產(chǎn)量均隨養(yǎng)殖數(shù)量的增多而提高, 這一變化趨勢(shì)與Bai等[36]研究結(jié)論一致。與此同時(shí), 畜禽糞尿等排泄物隨之劇增, 規(guī)?；呐盼叟涮自O(shè)施沒有及時(shí)完善, 大部分糞污僅簡單處理甚至不處理直接施入周邊作物或荒地[23], 導(dǎo)致單位農(nóng)田土壤氮盈余、農(nóng)業(yè)系統(tǒng)氨揮發(fā)增加, 與馬怡斐等[37]研究結(jié)果一致。2017年瓊海市單位耕地面積畜禽承載量為14 LU×hm-2, 由于農(nóng)牧生產(chǎn)養(yǎng)分循環(huán)脫節(jié), 土壤氮盈余已達(dá)665.2 kg×hm-2, 是2014年河北欒城[37]土壤氮盈余的2.4倍(該地區(qū)單位耕地面積畜禽承載量為18 LU×hm-2)。海南沿海地區(qū)經(jīng)濟(jì)作物的種植與養(yǎng)殖更為密集, 中部山區(qū)熱帶雨林生態(tài)功能區(qū)耕作困難且開發(fā)受限, 資源投入與環(huán)境排放量水平低于沿海區(qū)域。

3)海南農(nóng)業(yè)生產(chǎn)水平雖有小幅提高, 各市縣仍處于較低水平。蔬菜水果單產(chǎn)水平較低的原因復(fù)雜, 主要有以下幾點(diǎn): 海南的酸性土壤(平均pH低于5.6[32])易導(dǎo)致肥力不足和作物養(yǎng)分利用效率不高, 蔬菜瓜果因不耐高溫不適宜夏季種植, 全年病蟲害高發(fā)且易受臺(tái)風(fēng)天氣干擾造成損失。以蔬菜為例, 2017年海南島蔬菜單產(chǎn)水平提高至22.9 kg×hm-2, 仍僅為全國水平(34.6 kg×hm-2)[35]的2/3。

4 結(jié)論

海南島農(nóng)業(yè)綠色發(fā)展水平在經(jīng)濟(jì)、社會(huì)和生產(chǎn)力方面呈上升趨勢(shì), 南部沿海地區(qū)要優(yōu)于北部地區(qū); 但就資源投入與環(huán)境排放方面來看呈下降趨勢(shì), 尤其是在沿海地區(qū), 環(huán)境損失情況較中部山區(qū)更為嚴(yán)峻。產(chǎn)生這些變化的原因主要與種植結(jié)構(gòu)與養(yǎng)殖規(guī)模的轉(zhuǎn)變有關(guān), 海南省在追求經(jīng)濟(jì)增長過程中, 蔬果、熱帶特色經(jīng)濟(jì)作物播種比例和養(yǎng)殖密度不斷增加, 且管理措施的優(yōu)化有限, 因此帶來了環(huán)境問題。中部山區(qū)作為受保護(hù)的熱帶雨林生態(tài)功能區(qū), 開發(fā)受限, 資源投入及環(huán)境污染排放水平優(yōu)于沿海地區(qū)。

資源投入過高、養(yǎng)分利用效率低和環(huán)境污染排放較大是限制海南島農(nóng)業(yè)綠色發(fā)展的主要問題。為實(shí)現(xiàn)海南地區(qū)農(nóng)業(yè)綠色發(fā)展, 今后的發(fā)展規(guī)劃中應(yīng)兼顧農(nóng)業(yè)發(fā)展與環(huán)境保護(hù)的協(xié)調(diào)統(tǒng)一。對(duì)于種植業(yè)發(fā)展, 首先需要協(xié)調(diào)經(jīng)濟(jì)作物與糧食作物種植比例, 通過減施化肥、農(nóng)藥和農(nóng)膜以減少環(huán)境風(fēng)險(xiǎn); 養(yǎng)殖業(yè)應(yīng)優(yōu)化區(qū)域布局, 合理發(fā)展畜牧業(yè)規(guī)模, 加強(qiáng)對(duì)畜禽糞污的管理, 提高資源的循環(huán)利用效率; 協(xié)調(diào)農(nóng)牧關(guān)系, 促進(jìn)農(nóng)牧可持續(xù)發(fā)展, 提高養(yǎng)分利用效率和減少環(huán)境排放。

[1] MA L, ZHANG W F, MA W Q, et al. An analysis of deve-lopments and challenges in nutrient management in China[J]. Journal of Environmental Quality, 2013, 42(4): 951–961

[2] GUO J H, LIU X J, ZHANG Y, et al. Significant acidification in major Chinese croplands[J]. Science, 2010, 327(5968): 1008–1010

[3] GIORDANO J C P, BRUSH M J, ANDERSON I C. Quantifying annual nitrogen loads to Virginia’s coastal lagoons: Sources and water quality response[J]. Estuaries and Coasts, 2011, 34(2): 297–309

[4] ZHOU J Y, GU B J, SCHLESINGER W H, et al. Significant accumulation of nitrate in Chinese semi-humid croplands[J]. Scientific Reports, 2016, 6: 25088

[5] ZHENG H, WANG L J, PENG W J, et al. Realizing the values of natural capital for inclusive, sustainable development: Informing China’s new ecological development strategy[J]. Proceedings of the National Academy of Sciences of the United States of America, 2019, 116(17): 8623–8628

[6] MA L, BAI Z H, MA W Q, et al. Exploring future food provision scenarios for China[J]. Environmental Science & Technology, 2019, 53(3): 1385–1393

[7] 韓長賦. 大力推進(jìn)農(nóng)業(yè)綠色發(fā)展[J]. 農(nóng)村實(shí)用技術(shù), 2018, (3): 12–13 HAN C F. Vigorously promote green development of agriculture[J]. Applicable Technologies for Rural Areas, 2018, (3): 12–13

[8] 李谷成. 中國農(nóng)業(yè)的綠色生產(chǎn)率革命: 1978—2008年[J]. 經(jīng)濟(jì)學(xué), 2014, 13(2): 537–558 LI G C. The green productivity revolution of agriculture in China from 1978 to 2008[J]. China Economic Quarterly, 2014, 13(2): 537–558

[9] 金書秦, 牛坤玉, 韓冬梅. 農(nóng)業(yè)綠色發(fā)展路徑及其“十四五”取向[J]. 改革, 2020, 33(2): 30–39 JIN S Q, NIU K Y, HAN D M. The path of agricultural green development and its orientation in the 14th five-year plan period[J]. Reform, 2020, 33(2): 30–39

[10] SHEN J B, ZHU Q C, JIAO X Q, et al. Agriculture green development: A model for China and the world[J]. Frontiers of Agricultural Science and Engineering, 2020, 7(1): 5–13

[11] PéTRY F. Sustainability Issues in Agricultural and Rural Development Policies: Vol. 2: Traineer’s-Kit[M]. Roma: FAO, 1995

[12] 高薇, 汪嵐, 欒勝基. 歐盟農(nóng)業(yè)環(huán)境足跡指數(shù)法述評(píng)[J]. 生態(tài)經(jīng)濟(jì), 2015, 31(9): 176–180 GAO W, WANG L, LUAN S J. Review of the EU agri-environmental footprint index method[J]. Ecological Economy, 2015, 31(9): 176–180

[13] 魏琦, 張斌, 金書秦. 中國農(nóng)業(yè)綠色發(fā)展指數(shù)構(gòu)建及區(qū)域比較研究[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2018, (11): 11–20 WEI Q, ZHANG B, JIN S Q. A study on construction and regional comparison of agricultural green development index in China[J]. Issues in Agricultural Economy, 2018, (11): 11–20

[14] 袁久和. 我國農(nóng)村綠色發(fā)展水平與影響因素的實(shí)證分析[J]. 山西農(nóng)業(yè)大學(xué)學(xué)報(bào): 社會(huì)科學(xué)版, 2019, 18(6): 46–53 YUAN J H. Empirical analysis of green development level and influencing factors in rural areas of China[J]. Journal of Shanxi Agricultural University: Social Science Edition, 2019, 18(6): 46–53

[15] 楊世琦. 基于歐氏距離的農(nóng)業(yè)可持續(xù)發(fā)展評(píng)價(jià)理論構(gòu)建與實(shí)例驗(yàn)證[J]. 生態(tài)學(xué)報(bào), 2017, 37(11): 3840–3848 YANG S Q. A case study of a novel sustainable agricultural development evaluation method based on Euclidean distance theory[J]. Acta Ecologica Sinica, 2017, 37(11): 3840–3848

[16] LIU Y F, SUN D S, WANG H J, et al. An evaluation of China’s agricultural green production: 1978–2017[J]. Journal of Cleaner Production, 2020, 243: 118483

[17] 新華社. 中共中央國務(wù)院關(guān)于支持海南全面深化改革開放的指導(dǎo)意見[EB/OL]. 北京: 新華社, (2018-04-14). http://www.gov.cn/zhengce/2018-04/14/content_5282456.htm Xinhua News. Guiding opinions of the CPC Central Committee and the State Council on supporting Hainan’s comprehensive deepening reform and opening-up[EB/OL]. Beijing: Xinhua News, (2018-04-14). http://www.gov.cn/zhengce/ 2018-04/14/content_5282456.htm

[18] 張建杰, 崔石磊, 馬林, 等. 中國農(nóng)業(yè)綠色發(fā)展指標(biāo)體系的構(gòu)建與例證[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2020, DOI: 10.13930/j.cnki.cjea.200069 ZHANG J J, CUI S L, MA L, et al. Construction of a green development index system for agriculture in China and examples[J]. Chinese Journal of Eco-Agriculture, 2020, DOI: 10.13930/j.cnki.cjea.200069

[19] 海南省統(tǒng)計(jì)局, 國家統(tǒng)計(jì)局海南調(diào)查總隊(duì). 海南統(tǒng)計(jì)年鑒[M]. 北京: 中國統(tǒng)計(jì)出版社, 2018 Hainan Provincial Bureau of Statistics, Hainan Survey Team of National Bureau of Statistics. Hainan Statistical Yearbook 2018[M]. Beijing: China Statistics Press, 2018

[20] MA L, MA W Q, VELTHOF G F, et al. Modeling nutrient flows in the food chain of China[J]. Journal of Environmental Quality, 2010, 39(4): 1279–1289

[21] 劉曉利. 我國“農(nóng)田-畜牧-營養(yǎng)-環(huán)境”體系氮素養(yǎng)分循環(huán)與平衡[D]. 保定: 河北農(nóng)業(yè)大學(xué), 2005 LIU X L. Nitrogen cycling and balance in “agriculture-livestock-nutrition-environment” system in China[D]. Baoding: Hebei Agricultural University, 2005

[22] 阮云澤, 張茂星, 陳鵬, 等. 熱帶人工橡膠林地磚紅壤中氨揮發(fā)規(guī)律的研究[J]. 土壤, 2014, 46(3): 466–469 RUAN Y Z, ZHANG M X, CHEN P, et al. Characterization of ammonia volatilization from latosol soils in rubber plantation in tropical region[J]. Soils, 2014, 46(3): 466–469

[23] 海南省農(nóng)業(yè)農(nóng)村廳. 關(guān)于加大支持儋州市推進(jìn)畜禽養(yǎng)殖糞污資源化利用的建議[EB/OL]. (2018-04-09). http://www. hainan.gov.cn/rdjydata-6527.html Hainan Provincial Department of Agriculture. Suggestions on strengthening the support of Danzhou to promote the utilization of manure from livestock and poultry[EB/OL]. (2018-04-09). http://www.hainan.gov.cn/rdjydata-6527.html

[24] 海南省農(nóng)業(yè)廳.《海南省“十三五”秸稈綜合利用實(shí)施方案》的公告[EB/OL]. (2017-12-19). https://www. nongjitong. com/ news/ 2017/426707.htmlHainan Provincial Department of Agriculture. Implementation plan of comprehensive utilization of straw in the “Thirteenth Five-Year Plan” in Hainan Province[EB/OL]. (2017-12-19). https://www.nongjitong.com/news/2017/426707.html

[25] 國家發(fā)展和改革委員會(huì)價(jià)格司. 全國農(nóng)產(chǎn)品成本收益資料匯編2018[M]. 北京: 中國統(tǒng)計(jì)出版社, 2018 National Development and Reform Commission of the People’ Repbulic of China. Compilation of National Agricultural Product Cost-Benefit Data[M]. Beijing: China Statistics Press, 2018

[26] 冀宏杰, 張認(rèn)連, 武淑霞, 等. 太湖流域農(nóng)田肥料投入與養(yǎng)分平衡狀況分析[J]. 中國土壤與肥料, 2008, (5): 70–75 JI H J, ZHANG R L, WU S X, et al. Analysis of fertilizer input and nutrient balance of farmland in Taihu watershed[J]. Soil and Fertilizer Sciences in China, 2008, (5): 70–75

[27] YAN Z J, LIU P P, LI Y H, et al. Phosphorus in China’s intensive vegetable production systems: Overfertilization, soil enrichment, and environmental implications[J]. Journal of Environmental Quality, 2013, 42(4): 982–989

[28] 武良, 張衛(wèi)峰, 陳新平, 等. 中國農(nóng)田氮肥投入和生產(chǎn)效率[J]. 中國土壤與肥料, 2016, (4): 76–83 WU L, ZHANG W F, CHEN X P, et al. Nitrogen fertilizer input and nitrogen use efficiency in Chinese farmland[J]. Soil and Fertilizer Sciences in China, 2016, (4): 76–83

[29] 朱兆良. 中國土壤氮素研究[J]. 土壤學(xué)報(bào), 2008, 45(5): 778–783 ZHU Z L. Research on soil nitrogen in China[J]. Acta Pedologica Sinica, 2008, 45(5): 778–783

[30] 丁尚, 郭浩浩, 程鳴宇, 等. 海南農(nóng)牧生產(chǎn)體系氮素養(yǎng)分流動(dòng)時(shí)空變化特征[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2019, 27(2): 246–256 DING S, GUO H H, CHENG M Y, et al. Spatial and temporal characteristics of nitrogen flows in farming and animal production system in Hainan Province[J]. Chinese Journal of Eco-Agriculture, 2019, 27(2): 246–256

[31] OENEMA O, KROS H, DE VRIES W. Approaches and uncertainties in nutrient budgets: Implications for nutrient management and environmental policies[J]. European Journal of Agronomy, 2003, 20(1/2): 3–16

[32] 鄧萬剛. 海南島土壤酸度數(shù)據(jù)庫的建立及土壤酸度變異研究[D]. 儋州: 華南熱帶農(nóng)業(yè)大學(xué), 2005 DENG W G. Foundation of soil acidity database and study on variation of soil acidity in Hainan island[D]. Danzhou: South China University of Tropical Agriculture, 2005

[33] 趙青云, 王輝, 莊輝發(fā), 等. 海南香草蘭園土壤酸化現(xiàn)狀及酸化原因分析[J]. 熱帶農(nóng)業(yè)科技, 2014, 37(1): 12–13 ZHAO Q Y, WANG H, ZHUANG H F, et al. Investigation and analysis of the soil acidification of vanilla orchard in Hainan[J]. Tropical Agricultural Science & Technology, 2014, 37(1): 12–13

[34] 梁偉紅, 羅微, 劉燕群, 等. 海南冬季瓜菜產(chǎn)業(yè)風(fēng)險(xiǎn)及可持續(xù)發(fā)展對(duì)策研究[J]. 江蘇農(nóng)業(yè)科學(xué), 2013, 41(10): 431–433 LIANG W H, LUO W, LIU Y Q, et al. Study on risks and sustainable development countermeasures of melon-vegetable industry in Hainan Province[J]. Jiangsu Agricultural Sciences, 2013, 41(10): 431–433

[35] 中華人民共和國國家統(tǒng)計(jì)局. 中國統(tǒng)計(jì)年鑒[M]. 北京: 中國統(tǒng)計(jì)出版社, 2019 National Bureau of Statistics of the People’s Republic of China. China Statistical Yearbook[M]. Beijing: China Statistics Press, 2019

[36] BAI Z H, MA L, QIN W, et al. Changes in pig production in China and their effects on nitrogen and phosphorus use and losses[J]. Environmental Science & Technology, 2014, 48(21): 12742–12749

[37] 馬怡斐, 柏兆海, 馬林, 等. 欒城城郊型農(nóng)牧系統(tǒng)養(yǎng)分流動(dòng)與環(huán)境排放時(shí)空特征[J]. 中國農(nóng)業(yè)科學(xué), 2018, 51(3): 493–506MA Y F, BAI Z H, MA L, et al. Temporal and spatial changes of nutrient flows and losses in the Peri-Urban crop-livestock system in Luancheng[J]. Scientia Agricultura Sinica, 2018, 51(3): 493–506

Spatial and temporal characteristics of agricultural green development indicators in Hainan Island*

SONG Chenyang1,2,ZHANG Jianjie3,LIU Ling4,MA Wenqi3,MA Lin4, DING Shang2, ZHAO Hongwei2**

(1. College of Forestry, Hainan University, Haikou 570228, China; 2. College of Ecology and Environment, Hainan University / Key Laboratory of A&F Environmental Processes and Ecological Regulation of Hainan Province, Haikou 570228, China; 3. College of Resources and Environmental Sciences, Hebei Agricultural University, Baoding 071001, China; 4. Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences, Shijiazhuang 050022, China)

Since the establishment of Hainan Province, the scale, structure, and methods of agricultural production have undergone significant changes, which are critical with respect to the progress of agricultural green development (AGD). In this study we adopted China’s AGD indicator system with 20 indicators, and using the NUFER (NUtrient flows in Food chains, Environment and Resources use) model, we quantitatively calculated the spatial and temporal changes in Hainan Island AGD indicators from 1988 to 2017. Moreover, we discussed the factors restricting AGD in Hainan Island. The results revealed that the overall level of Hainan Island’s AGD had shown a declining trend over the past 30 years. The numbers of indicators at lower levels (Ⅲ and Ⅳ) had increased from 12 to 15, whereas those at higher levels (Ⅰ and Ⅱ) had decreased from 8 to 5. With increase in the proportion of economic crops such as vegetables and fruits, the scale of livestock and poultry farming and the degree of intensification, Hainan’s agricultural output value and social development level had continued to increase. Furthermore, the agricultural output value per unit cultivated land area and level of agricultural mechanization had increased from 6 200 ￥?hm–2and 3 kW?hm–2to 161 000 ￥?hm–2and 13 kW?hm–2, respectively. Over the past 30 years, there had been an increasingly excessive input of agricultural resources (nitrogen, phosphorus, pesticides, and agricultural film), among which, the intensities of pesticide and agricultural film usage had shown the highest increases, rising from 8.0 kg?hm–2and 0 kg?hm–2to 41.9 kg?hm–2and 34.7 kg?hm–2, respectively. As a consequence of the low productivity levels, the environmental emissions of applied nutrients (e.g., nitrogen surplus in farmland, nitrogen runoff in farmland, and NH3volatilization from agricultural system) had increased significantly, with NH3volatilization and nitrogen surplus per unit sown area showing the largest increases, rising from 61.0 kg?hm–2and 152.1 kg?hm–2to 131.4 kg?hm–2and 297.9 kg?hm–2, respectively. Environmental pollution caused by an excessive input of agricultural resources and the separation of agriculture and animal husbandry had been the main factors contributing to the restriction of AGD in Hainan Island. Given that the central part of the island is covered in montane forest and has been designated as a national key ecological function protection area, the scale of development in this region had been relatively limited. Consequently, the coastal plains had experience more rapid growth in resource input. Although this had a detrimental impact on the environment, the economy growth, social development, and agricultural productivity indicators for these areas had improved to a greater extent than those for mountain forest area. In the future, Hainan Island should rationally optimize the planting structure in coastal areas, improve the level of production management, and strengthen the combination of agriculture and animal husbandry to reduce resource losses and achieve the island’s requisite AGD goals.

Agricultural green development; NUFER model; Change in planting structure; Agriculture and animal husbandry system; Tropical agriculture

S19; S158.5

10.13930/j.cnki.cjea.200011

宋晨陽, 張建杰, 劉玲, 馬文奇, 馬林, 丁尚, 趙洪偉. 海南島農(nóng)業(yè)綠色發(fā)展指標(biāo)時(shí)空變化特征[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2020, 28(8): 1156-1167

SONG C Y, ZHANG J J, LIU L, MA W Q, MA L, DING S, ZHAO H W. Spatial and temporal characteristics of agricultural green development indicators in Hainan Island[J]. Chinese Journal of Eco-Agriculture, 2020, 28(8): 1156-1167

* 海南省重大科技項(xiàng)目(ZDKJ2017002)資助

趙洪偉, 研究方向?yàn)槲廴旧鷳B(tài)學(xué)。E-mail: hwzhao@hainu.edu.cn

宋晨陽, 研究方向?yàn)閰^(qū)域養(yǎng)分資源管理。E-mail: Cysong1995@163.com

2020-01-06

2020-05-15

* This study was supported by the Major Science and Technology Program of Hainan Province (ZDKJ2017002).

, E-mail: hwzhao@hainu.edu.cn

Jan. 6, 2020;

May 15, 2020