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

    高集約化農(nóng)區(qū)投入減量化與環(huán)境風(fēng)險降低潛勢的時空分異特征

    2017-10-13 21:02:53付永虎劉黎明王加升葉津煒
    農(nóng)業(yè)工程學(xué)報 2017年2期
    關(guān)鍵詞:灰水青浦區(qū)足跡

    付永虎,劉黎明,王加升,葉津煒,郭 赟

    ?

    高集約化農(nóng)區(qū)投入減量化與環(huán)境風(fēng)險降低潛勢的時空分異特征

    付永虎1,2,劉黎明1※,王加升1,葉津煒1,郭 赟1

    (1. 中國農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,北京 100193; 2. 淮海工學(xué)院法學(xué)院,連云港 222005)

    該文應(yīng)用氮足跡(nitrogen footprint,N footprint)、灰水足跡(grey water footprint,GWF)理論,以上海市青浦區(qū)為案例區(qū),對集約農(nóng)業(yè)土地利用系統(tǒng)的活性氮排放和非點源污染進行分析,運用非期望產(chǎn)出的SBM-Undesirable窗式分析模型,系統(tǒng)解析農(nóng)業(yè)土地利用環(huán)境效率時空分異特征,綜合評估農(nóng)業(yè)土地利用系統(tǒng)投入減量化趨勢與環(huán)境風(fēng)險降低潛勢。結(jié)果表明:1)在考慮活性氮排放和非點源污染約束下,2006-2013年青浦區(qū)農(nóng)業(yè)土地利用環(huán)境效率相對較低,2006-2013年均值僅為0.669;2)2006-2013年青浦區(qū)年均勞動力、肥料和機械動力的潛在減量比例較高,8 a間潛在年均縮減總量分別為8 104人,4 501.59 t,27 928.44 kW;非點源污染灰水足跡的潛在減排比例高于污染氮足跡,潛在減排總量年均分別為52 046.88萬m3和381.04 t?;ㄏ銟蚪值谰哂凶畲蟮臐撛谕度霚p量化與環(huán)境風(fēng)險降低比例。白鶴鎮(zhèn)、練塘鎮(zhèn)等具有較大的潛在縮減規(guī)模,應(yīng)成為青浦區(qū)農(nóng)業(yè)投入減量化和污染物減排總量控制的重點區(qū)域。該文評價結(jié)果可為制定農(nóng)業(yè)可持續(xù)發(fā)展規(guī)劃及農(nóng)業(yè)產(chǎn)業(yè)政策提供參考。

    土地利用;污染;農(nóng)業(yè);氮足跡;灰水足跡;環(huán)境效率;投入減量;環(huán)境風(fēng)險

    0 引 言

    過去30多年來,中國普遍實行了高投入-產(chǎn)出的集約農(nóng)業(yè)土地利用模式,化肥、農(nóng)業(yè)機械的廣泛使用在提高土地收益和保障糧食安全的同時,也帶來了活性氮(reactive Nitrogen,Nr)排放不斷增加、農(nóng)業(yè)非點源污染趨于嚴重等問題,集約化農(nóng)業(yè)生產(chǎn)對環(huán)境的負面影響已受到廣泛的關(guān)注。已有研究報道全球大約有70%以上的Nr排放來自于農(nóng)業(yè)生產(chǎn)[1];在中國,農(nóng)業(yè)土地利用過程導(dǎo)致的灰水足跡總量占到人類活動的60%[2]。減少Nr排放和控制農(nóng)業(yè)非點源污染風(fēng)險,已成為農(nóng)業(yè)可持續(xù)發(fā)展領(lǐng)域的重要研究方向。從農(nóng)業(yè)土地利用系統(tǒng)的角度來說,系統(tǒng)輸入和輸出的物質(zhì)量越大,人類給予生態(tài)環(huán)境的壓力就愈大。化肥、機械等投入要素的無效或低效率利用是導(dǎo)致農(nóng)業(yè)環(huán)境問題突出的主要原因之一[3],因此,通過提高農(nóng)業(yè)土地利用系統(tǒng)的環(huán)境效率,降低物質(zhì)投入與減少有害物質(zhì)的環(huán)境輸出,是實現(xiàn)農(nóng)業(yè)生態(tài)經(jīng)濟系統(tǒng)穩(wěn)步運行的關(guān)鍵。

    環(huán)境效率是衡量經(jīng)濟發(fā)展導(dǎo)致的環(huán)境代價,是可持續(xù)發(fā)展能力評價的重要指標(biāo)之一。環(huán)境效率的常用測度方法是假設(shè)一個經(jīng)濟體在等量要素或產(chǎn)出條件下,其污染排放離污染最小排放的距離[4]。土地利用系統(tǒng)環(huán)境效率分析應(yīng)既考慮農(nóng)業(yè)土地利用過程中產(chǎn)出的經(jīng)濟價值,同時也應(yīng)兼顧土地利用過程中的環(huán)境輸出,即需要考慮期望產(chǎn)出,也應(yīng)考慮非期望產(chǎn)出。有關(guān)農(nóng)業(yè)土地利用環(huán)境效率研究,不少學(xué)者從不同角度進行了深入探討并取得了諸多研究成果[5-9]。例如Marthin等[10]回顧了考慮環(huán)境效應(yīng)的農(nóng)業(yè)生產(chǎn)效率計算方法,以澳大利亞的農(nóng)業(yè)生產(chǎn)部門為案例研究,考慮環(huán)境成本的投入,對農(nóng)業(yè)生產(chǎn)效率進行了修正,研究結(jié)果表明考慮環(huán)境因素的生產(chǎn)率校準(zhǔn)方法是可持續(xù)發(fā)展研究中一個重要的生產(chǎn)率測定方法。近年來非期望產(chǎn)出的SBM方法的提出,不僅解決了投入和產(chǎn)出松弛的問題,且同時解決了非期望產(chǎn)出存在下的環(huán)境效率評價問題;是一種綜合考慮土地、勞動力、機械、肥料等投入與經(jīng)濟產(chǎn)出、環(huán)境輸出等要素相互作用的農(nóng)業(yè)環(huán)境效率的表達,能夠科學(xué)地反映各個決策單元環(huán)境效率的真實水平[11],已被眾多學(xué)者應(yīng)用于農(nóng)業(yè)環(huán)境效率和節(jié)能減排潛勢分析。Kuo等[12]以臺灣南部的58個村種植業(yè)為研究對象,應(yīng)用SBM-Undesirable模型將農(nóng)用化學(xué)品污染物輸出作為非期望產(chǎn)出,計算了各決策單元的環(huán)境效率。國內(nèi)學(xué)者分別從非點源污染、碳排放以及2者結(jié)合的視角,對國家、區(qū)域等層面進行農(nóng)業(yè)環(huán)境效率的測算,為農(nóng)業(yè)生產(chǎn)的節(jié)能減排、環(huán)境效率的提高提出了改進方向,研究成果對農(nóng)業(yè)環(huán)境管理與政策制定提供了參考依據(jù)[6,13-16]。然而,以上有關(guān)農(nóng)業(yè)環(huán)境效率研究,對非期望產(chǎn)出指標(biāo)的考查主要集中在農(nóng)業(yè)非點源污染、碳排放等因素,忽視了對農(nóng)業(yè)土地利用過程導(dǎo)致的活性氮排放的研究,以致丟失了大量的環(huán)境輸出信息,從而降低了對農(nóng)業(yè)土地利用環(huán)境效率評價的有效性和準(zhǔn)確性。其次,非點源污染的測算主要以等標(biāo)排放量為主,未能充分反映非點源污染對當(dāng)?shù)厮Y源的影響,研究成果無法科學(xué)指導(dǎo)關(guān)鍵環(huán)境風(fēng)險的縮減潛勢。作為當(dāng)前生態(tài)經(jīng)濟學(xué)和可持續(xù)發(fā)展研究領(lǐng)域的熱點之一,氮足跡和灰水足跡為評估農(nóng)業(yè)土地利用資源消耗和廢棄物排放等提供了新的理念和途徑,氮足跡和灰水足跡理論的雙重足跡分析方法應(yīng)用于農(nóng)業(yè)生產(chǎn)的環(huán)境效應(yīng)評價,能夠準(zhǔn)確量化農(nóng)業(yè)土地利用過程對環(huán)境的多重負面影響[17];同時還可為農(nóng)業(yè)土地利用環(huán)境效率評價提供全新的定量化指標(biāo)。基于此,本文主要將在以下3個方面進行拓展:1)系統(tǒng)分析農(nóng)業(yè)生產(chǎn)過程的污染氮足跡和非點源污染灰水足跡,綜合評估集約農(nóng)業(yè)土地利用的環(huán)境效應(yīng);2)將污染氮足跡和非點源污染灰水足跡納入農(nóng)業(yè)土地利用環(huán)境效率評價模型,從而更加全面、客觀地評價集約農(nóng)業(yè)土地利用系統(tǒng)的環(huán)境效率;3)通過構(gòu)建投入減量化潛力、規(guī)模,污染物減排潛力、規(guī)模模型,探討集約化農(nóng)區(qū)投入減量化與環(huán)境風(fēng)險降低潛勢。研究結(jié)果可為構(gòu)建投入減量化與低環(huán)境風(fēng)險的可持續(xù)集約農(nóng)業(yè)土地利用模式提供支撐,為政府制定農(nóng)業(yè)環(huán)境政策提供理論與方法依據(jù)。

    1 研究區(qū)概況

    本研究選取上海市青浦區(qū)作為研究單元,該區(qū)域位于上海市西郊,太湖下游,黃浦江上游。地理坐標(biāo)為30°59′~31°16′N、120°53′~121°17′E,境內(nèi)河網(wǎng)密集,淀山湖為上海市主要水源地之一。青浦區(qū)氣候溫和濕潤、日照充足,多年平均氣溫16.2 ℃,屬亞熱帶季風(fēng)氣候,多年平均降雨量1 049.1 mm。土地總面積約670.14 km2,占上海市總面積的10%,轄區(qū)地勢平坦,土地肥沃,耕地面積比重大,農(nóng)業(yè)生產(chǎn)集約化程度高,是上海市重要的糧食、蔬菜和水果生產(chǎn)基地。近年來,隨著高投入、高產(chǎn)出農(nóng)業(yè)生產(chǎn)模式的推行,Nr排放及農(nóng)業(yè)非點源污染的風(fēng)險日益增大。

    2 研究方法與參數(shù)選取

    2.1 農(nóng)業(yè)土地利用系統(tǒng)污染氮足跡和非點源污染灰水足跡的測算

    1)測算方法

    污染氮足跡指在農(nóng)業(yè)土地利用過程中所導(dǎo)致的不被作物吸收的活化氮排放[17],包含地表徑流(TN)、淋溶(TN)、氨揮發(fā)、N2O直接排放等虛擬氮足跡(virtual N footprint,NFvirtual)及農(nóng)業(yè)燃油NO和NO2排放途徑損失的農(nóng)業(yè)能源氮足跡(energy N footprint,NFenergy)。計算公式為

    式中NFpollution代表污染氮足跡,t;runoff代表徑流途徑損失的總氮,t;leach代表淋溶途徑損失的總氮,t;volatilization代表氨揮發(fā)損失的氮,t;N2O代表N2O直接排放的氮,t;fuel代表農(nóng)業(yè)燃油NO與NO2排放途徑損失的Nr,t。針對非點源污染灰水足跡的計算,本文借鑒付永虎等[18]對農(nóng)業(yè)生產(chǎn)部門灰水足跡的計算方法,計算區(qū)域農(nóng)業(yè)非點源污染灰水足跡。

    2)參數(shù)與數(shù)據(jù)來源

    表1 參數(shù)取值與數(shù)據(jù)來源

    續(xù)表

    2.2 SBM-Undesirable模型的窗式分析

    數(shù)據(jù)包絡(luò)分析方法(data envelopment analysis, DEA)是一種可測度多要素投入與多產(chǎn)出之間相對效率評價的系統(tǒng)分析方法,因其具有非主觀賦權(quán)、無需事先確定函數(shù)關(guān)系,并可對決策單元無效因素進行分析等優(yōu)點,在擁有相同結(jié)構(gòu)決策單元的復(fù)雜生產(chǎn)系統(tǒng)效率評價中具有巨大優(yōu)勢,自1978年以來已廣泛應(yīng)用于效率評價中。然而,目前DEA方法的應(yīng)用主要集中于CCR、BBC等傳統(tǒng)模型,產(chǎn)出指標(biāo)一般為期望產(chǎn)出,對污染物輸出沒有納入考慮范圍,且效率評價過程忽視了投入產(chǎn)出變量的松弛型問題,無法對無效率情況進行分析,導(dǎo)致了相對效率評價有失偏誤[27]。2001年,Tone提出了SBM(Slack-Based Measure)模型[28],有效解決了無效率情況測算的問題,并可分析松弛型問題的可改進效率值。此后,Cooper等[29]將非期望產(chǎn)出引入到SBM模型中,構(gòu)建SBM-Undesirable模型,從而實現(xiàn)了非期望產(chǎn)出條件下集約農(nóng)業(yè)土地利用環(huán)境效率評價。假設(shè)生產(chǎn)系統(tǒng)有個決策單元,每個單元有種投入、1種期望產(chǎn)出和2種非期望產(chǎn)出,針對一個包含非期望產(chǎn)出的生產(chǎn)可能性集合

    構(gòu)建SBM-Undesirable公式如下

    式中代表有界閉集,∈R,yR1與yR2;定義矩陣、Y、Y為:=[1,…,n]∈RY=[y,…,yn]∈RY=[y,…,y]∈R,>0,Y>0,Y>0。代表權(quán)重向量,∈R,≥0表示規(guī)模報酬可變(variable return to scale,簡稱VRS)。代表環(huán)境效率值,0≤≤1;0、0、0分別代表投入、期望產(chǎn)出和非期望產(chǎn)出要素;、1、2分別代表投入、期望產(chǎn)出和非期望產(chǎn)出要素的個數(shù);-、g、s分別代表投入、期望產(chǎn)出和非期望產(chǎn)出的松弛變量;λ代表權(quán)重向量。當(dāng)-=0、g=0、s=0時,=1,決策單元(decision making unit,DMU)即為有效率的。

    窗式分析由Charnes等[30]提出,通過把處于不同時間段上的同一單元看作不同的決策單元,采用特定時間段作為窗口寬度構(gòu)建生產(chǎn)前沿面。窗口隨時間序列向前滑動,每滑動一次就將最早的一個時段數(shù)據(jù)從窗口中去掉,同時增加一個新時段數(shù)據(jù),因此既能避免了傳統(tǒng)DEA分析中生產(chǎn)前沿面不連續(xù)的問題,同時也規(guī)避了跨期分析導(dǎo)致的生產(chǎn)前沿面相同的缺陷,使得對面板數(shù)據(jù)的分析更加靈活,效率計算更為準(zhǔn)確和客觀。假設(shè)有個時期(=1,…,),若窗口的寬度為(≤),則可得到=(?+1)個窗口。本文基于已有研究將窗口寬度設(shè)為=3[31],同時依據(jù)相關(guān)文獻的研究成果[13,15-16],選取影響農(nóng)業(yè)生產(chǎn)的關(guān)鍵投入和產(chǎn)出指標(biāo),構(gòu)建農(nóng)業(yè)土地利用環(huán)境效率指標(biāo)體系(表2)。農(nóng)業(yè)產(chǎn)值按照2006年不變價列入計算,非期望產(chǎn)出選取污染氮足跡和農(nóng)業(yè)非點源污染灰水足跡。將青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)農(nóng)業(yè)土地利用系統(tǒng)2006-2013年的各投入產(chǎn)出的面板數(shù)據(jù)代入SBM-Undesirable的窗式分析模型,以評價不同DMU的相對環(huán)境效率。

    表2 青浦區(qū)農(nóng)業(yè)土地利用環(huán)境效率評價指標(biāo)體系

    2.3 物質(zhì)減量及環(huán)境風(fēng)險降低潛勢評估

    根據(jù)SBM-Undesirable模型,如果<1時,DMU無效(即存在效率損失),松弛變量可反映農(nóng)業(yè)土地利用過程中環(huán)境無效率的改善途徑。因此,構(gòu)建投入減量、輸出減排比例MRP和投入減量、輸出減排規(guī)模MRS,計算公式如下

    式中MRP代表第個DMU的某一投入指標(biāo)的投入減量比例或非期望產(chǎn)出的輸出減排比例,%;S代表第個DMU的投入冗余或非期望產(chǎn)出冗余的松弛變量;AP代表第個DMU的實際投入或產(chǎn)出量;MRS代表某一指標(biāo)的投入減量或者輸出減排規(guī)模。MRP是指DMU某一指標(biāo)投入冗余或非期望產(chǎn)出冗余的松弛變量與實際投入或產(chǎn)出量的比值,反映了農(nóng)業(yè)土地利用某一投入指標(biāo)或非期望產(chǎn)出指標(biāo)的減量投入或減排的改善空間。MRS是指單個DMU的某一指標(biāo)投入冗余或非期望產(chǎn)出冗余的松弛變量與所有DMU該指標(biāo)投入冗余或非期望產(chǎn)出冗余的松弛變量總和的比值,反映了農(nóng)業(yè)土地利用某一投入指標(biāo)或者非期望產(chǎn)出指標(biāo)對研究區(qū)該指標(biāo)減投或減排的潛在規(guī)模的貢獻和影響程度。

    3 結(jié)果分析

    3.1 污染氮足跡與非點源灰水足跡分析

    3.1.1 總體分析

    農(nóng)業(yè)污染氮足跡是農(nóng)業(yè)土地利用過程對活性氮輸出的度量。2006-2013年青浦區(qū)農(nóng)業(yè)污染氮足跡在2 805.59~3 505.68 t之間波動(表3),8 a平均值為3 092.44 t。與已有研究成果相似[17],氨揮發(fā)是污染氮足跡最主要的排放途徑,2006-2013年氨揮發(fā)損失的氮占污染氮足跡的比重平均為55.95%;淋溶和徑流損失的氮是引起農(nóng)田周邊水質(zhì)惡化的重要途徑,2006-2013年年均排向水體的Nr總量在1 126.50~1 329.95 t波動,年均值1 237.69 t。N2O直接排放與燃油釋放的活性氮相對較小,年均排放總量占污染氮足跡比重的3%以下。

    2006-2013年,青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)氮肥灰水足跡均高于磷肥灰水足跡,非點源污染灰水足跡取氮肥灰水足跡計算值。8 a間農(nóng)業(yè)非點源污染灰水足跡為134 026.58~152 833.70萬m3(表3),年均值為144 969.00萬m3,單位土地利用面積灰水足跡年均5.52萬m3/hm2,2010年單位耕地面積灰水足跡是湖南省桃江縣單位農(nóng)業(yè)土地利用系統(tǒng)灰水足跡的2.30倍。青浦區(qū)多年平均補給地表水資源總量由3部分組成,其中本地徑流為2.53億m3,占1.63%,上游太湖流域來水占45.45%,潮水占52.92%,其總量達到了155.1億m3,豐富的水資源保障了社會經(jīng)濟用水需求。然而應(yīng)用水環(huán)境壓力(water environmental pressure,WEP)指標(biāo)[18],分析多年平均徑流量情景下農(nóng)業(yè)非點源污染灰水足跡的水環(huán)境壓力(指農(nóng)業(yè)非點源污染灰水足跡與研究區(qū)多年平均徑流量的比值),2006-2013年青浦區(qū)水環(huán)境受農(nóng)業(yè)生產(chǎn)的影響較大,年均WEP為5.74,農(nóng)業(yè)生產(chǎn)對稀釋水的需求已遠遠超出了青浦區(qū)天然地表徑流量。因此,如果考慮生態(tài)用水,特別是農(nóng)業(yè)生產(chǎn)的環(huán)境用水,青浦區(qū)水資源則嚴重短缺。農(nóng)業(yè)生產(chǎn)導(dǎo)致的水質(zhì)惡化引起生態(tài)用水的增大,是上海青浦區(qū)水資源出現(xiàn)水質(zhì)性缺水的重要原因之一。

    表3 上海青浦區(qū)污染氮足跡與非點源污染灰水足跡的總體分析

    3.1.2 空間分析

    青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)農(nóng)業(yè)污染氮足跡跨度較大,除白鶴鎮(zhèn)和華新鎮(zhèn)外,總體呈“西高東低”的態(tài)勢(圖1a)。2006-2013年各鄉(xiāng)鎮(zhèn)(街道)農(nóng)業(yè)污染氮足跡為25.78~1 149.62 t,年均最大值和最小值的鄉(xiāng)鎮(zhèn)(街道)分別為白鶴鎮(zhèn)(898.47 t)和盈浦街道(38.72 t)。位于青西地區(qū)的金澤鎮(zhèn)、朱家角鎮(zhèn)和練塘鎮(zhèn)農(nóng)業(yè)污染氮足跡分別高于青浦區(qū)平均水平(281.13 t)8.84%、57.50%、69.05%。其他區(qū)域則呈兩極分化的特征,位于北部的白鶴鎮(zhèn)農(nóng)業(yè)污染氮足跡遠高于區(qū)域平均值219.59%,然而位于青浦中部和東部的盈浦街道、徐涇鎮(zhèn)農(nóng)業(yè)污染氮足跡卻在60 t以下。

    青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)農(nóng)業(yè)非點源污染灰水足跡為764.68~47 769.40萬m3,2006-2013年年均最大值和最小值的鄉(xiāng)鎮(zhèn)(街道)分別為朱家角鎮(zhèn)(39 830.94萬m3)和盈浦街道(1 077.92萬m3)(圖1b)。各鄉(xiāng)鎮(zhèn)(街道)農(nóng)業(yè)非點源污染灰水足跡跨度較大,除白鶴鎮(zhèn)外,整體呈現(xiàn)青西地區(qū)高于青東地區(qū)的空間分異格局。盈浦街道、趙巷鎮(zhèn)因農(nóng)業(yè)土地利用面積較小,灰水足跡在2 000萬m3以下。青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)磷肥灰水足跡為193.86~24 627.74萬m3,2006-2013年年均最大值和最小值的鄉(xiāng)鎮(zhèn)(街道)分別為朱家角鎮(zhèn)(20 474.05萬m3)和盈浦街道(334.91萬m3)。磷肥灰水足跡的空間分布與農(nóng)業(yè)非點源污染灰水足跡的空間部分特征類似,除白鶴鎮(zhèn)外整體呈現(xiàn)青西地區(qū)高于青東地區(qū)的空間分布特征。

    3.2 農(nóng)業(yè)土地利用環(huán)境效率評估

    3.2.1 青浦區(qū)農(nóng)業(yè)土地利用環(huán)境效率的年際變化

    本文應(yīng)用SBM-Undesirable模型的窗式分析方法,基于CRS假設(shè)計算出青浦區(qū)2006-2013年農(nóng)業(yè)土地利用環(huán)境效率(圖2)??傮w來看,2006-2013年農(nóng)業(yè)土地利用環(huán)境效率呈“雙峰”型波動的態(tài)勢,在0.612~0.743之間波動,年均值僅為0.669。最高值出現(xiàn)在2008年,為0.743,處于較低水平;而2013年農(nóng)業(yè)土地利用環(huán)境效率值為0.612,這意味著如果把當(dāng)前生產(chǎn)要素的潛力全部發(fā)揮出來,可以使農(nóng)業(yè)土地利用的環(huán)境效率再增長近38.8%。

    3.2.2 空間格局演變

    為直觀地分析研究期內(nèi)青浦區(qū)農(nóng)業(yè)土地利用環(huán)境效率的空間演化格局,利用2006、2008、2011和2013年各鄉(xiāng)鎮(zhèn)(街道)截面數(shù)據(jù),分析環(huán)境效率的空間格局演變(圖3)。按照涂正革對環(huán)境技術(shù)效率的分類方法[32],將青浦區(qū)各城鎮(zhèn)分為相對高效率區(qū)(0.9-1]、較高效率區(qū)(0.8-0.9]、中等效率區(qū)(0.7-0.8]、較低效率區(qū)(0.6-0.7]和低效率區(qū)(0-0.6]。總體來看,2006-2013年青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)土地利用環(huán)境效率整體呈現(xiàn)波動變化的趨勢。2006年農(nóng)業(yè)土地利用環(huán)境高效率鄉(xiāng)鎮(zhèn)(街道)的數(shù)量為4個,均位于青浦區(qū)東部,分別為:白鶴鎮(zhèn)、重固鎮(zhèn)、華新鎮(zhèn)和徐涇鎮(zhèn);在非有效鄉(xiāng)鎮(zhèn)中,平均效率為0.418,低于平均水平0.629,其中花香橋僅為0.261,農(nóng)業(yè)土地利用極為粗放。2008年,環(huán)境高效率的DMU增加至6個,金澤鎮(zhèn)也由2006年的低效率區(qū)上升為較低效率區(qū),效率值增至0.696。2011年環(huán)境高效率的DMU數(shù)量回落至4個,金澤鎮(zhèn)土地利用環(huán)境效率上升至高效率,白鶴鎮(zhèn)、重固鎮(zhèn)和夏陽街道下降至較低效率、低效率和低效率區(qū)。2013年青浦區(qū)除盈浦街道處于高效率區(qū)外,總體呈現(xiàn)下降的態(tài)勢,總體平均環(huán)境效率也由2011年的0.625下降至2013年的0.612。

    3.3 投入減量及環(huán)境風(fēng)險降低潛勢評估

    本研究根據(jù)公式(3)和公式(4),計算2006-2013年勞動力、肥料和機械的投入減量潛力和投入減量規(guī)模,污染氮足跡和非點源污染灰水足跡的減排潛力和減排規(guī)模。

    3.3.1 投入減量與污染物減排潛力

    青浦區(qū)2006-2013年年均勞動力、肥料和機械動力的潛在減量比例分別為32.21%,25.70%和38.21%,污染氮足跡和非點源污染灰水足跡的潛在減排比例分別為12.32%和32.18%。機械動力與勞動力的投入減量潛力相對較高,非點源污染灰水足跡減排潛力較大,機械和勞動力未能有效利用,資源配置相對低效,非點源污染問題突出。從各鄉(xiāng)鎮(zhèn)(街道)的投入減量與污染物減排潛力來看(圖4),青浦區(qū)各鄉(xiāng)鎮(zhèn)(街道)普遍存在投入冗余和環(huán)境污染的現(xiàn)象,投入規(guī)模和結(jié)構(gòu)尚不盡合理、資源無效利用、非期望產(chǎn)出過多等問題。香花橋街道存在勞動力過剩,肥料施用量過高,機械投入冗余,農(nóng)業(yè)污染氮足跡和非點源污染灰水足跡排放過高的現(xiàn)象,具有最大的投入減量化與環(huán)境風(fēng)險降低潛勢,應(yīng)成為青浦區(qū)農(nóng)業(yè)土地利用系統(tǒng)設(shè)計與優(yōu)化重點關(guān)注的區(qū)域。通過提高機械利用效率,優(yōu)化勞動力配置、科學(xué)適量的化肥投入模式以及增加技術(shù)投入減少Nr和非點源污染物排放等途徑來提高該區(qū)域的農(nóng)業(yè)土地利用環(huán)境效率,進而降低環(huán)境風(fēng)險。盈浦街道農(nóng)業(yè)生產(chǎn)的勞動力、肥料、機械動力和Nr排放指標(biāo)具有較大的投入減量和污染物減排潛力,在提高勞動生產(chǎn)率、減少肥料投入和提高機械利用效率的同時,應(yīng)重點實施Nr減排的相關(guān)措施。針對非點源污染灰水足跡排放的問題,除香花橋街道外,金澤鎮(zhèn)和朱家角鎮(zhèn)具有較大的農(nóng)業(yè)非點源污染減排潛勢,應(yīng)成為非點源污染物減排重點調(diào)控的區(qū)域。華新鎮(zhèn)農(nóng)業(yè)生產(chǎn)要素配置較為合理,投入減量與污染物減排潛勢相對較低,除勞動力的潛在縮減比例為11.55%外,其他指標(biāo)的相對縮減比例均在7%以內(nèi)。

    3.3.2 投入減量與污染物減排規(guī)模

    2006-2013年青浦區(qū)農(nóng)業(yè)生產(chǎn)的勞動力、肥料和機械動力的潛在縮減總量年均分別為8 104人,4 501.59 t,27 928.44 kW;污染氮足跡和非點源污染灰水足跡的潛在減排總量分別為381.04 t和52 046.88萬m3。從各鄉(xiāng)鎮(zhèn)(街道)的潛在縮減規(guī)模來看(表4),白鶴鎮(zhèn)和香花橋街道污染氮足跡的潛在減排規(guī)模對青浦區(qū)農(nóng)業(yè)Nr減排具有較大影響,分別占區(qū)域Nr潛在減排總量的26.39%和22.18%。練塘鎮(zhèn)勞動力、肥料、機械動力和灰水足跡的潛在縮減規(guī)模對青浦區(qū)的投入減量和非點源污染減排具有較大的影響,分別占青浦區(qū)潛在縮減總量的29.30%,34.59%,30.72%和35.56%,練塘鎮(zhèn)農(nóng)業(yè)生產(chǎn)主要以水生蔬菜種植為主,其高勞動力投入、化肥過量施用等土地利用模式造成了資源的低效利用,農(nóng)業(yè)非點源污染物輸出總量較高,應(yīng)成為青浦區(qū)農(nóng)業(yè)資源節(jié)約和污染減排總量控制的重點區(qū)域。

    表4 青浦區(qū)2006-2013年關(guān)鍵指標(biāo)的潛在縮減規(guī)模

    4 討 論

    高集約化農(nóng)業(yè)生產(chǎn)方式下不可避免地伴隨著污染物的排放,糧食安全與環(huán)境安全之間日益尖銳的矛盾,使得集約農(nóng)業(yè)的可持續(xù)發(fā)展面臨著嚴峻的挑戰(zhàn)。未來十年,高集約化農(nóng)業(yè)土地利用模式是保障中國糧食安全的主要途徑,在環(huán)境約束日益增強的情況下,提高土地利用環(huán)境效率意義重大。氮足跡可從關(guān)鍵元素(氮)的角度分析農(nóng)業(yè)土地利用過程中活性氮排放的遷移與轉(zhuǎn)化途徑,重點關(guān)注活性氮的排放總量與結(jié)構(gòu)特征,分析結(jié)果可為指導(dǎo)人類生產(chǎn)方式,減少人為Nr排放提供理論和數(shù)據(jù)支撐[33-34]。灰水足跡分析方法實現(xiàn)了從水量的角度評價水污染程度的目的,直觀地反映了水污染對可用水資源量的影響[35],其計算方法通常情況下為換算成將污水稀釋至符合當(dāng)?shù)貐^(qū)域規(guī)定的水質(zhì)標(biāo)準(zhǔn)所需的淡水體積,本文灰水足跡的計算包含氮、磷等污染物進入水體的“虛擬”水量,其計算結(jié)果不僅與氮徑流、淋溶損失有關(guān),還與磷流失及當(dāng)?shù)氐乃h(huán)境質(zhì)量標(biāo)準(zhǔn)相關(guān)。根據(jù)上海市對水源地保護的要求,不同鄉(xiāng)鎮(zhèn)執(zhí)行差異化的水環(huán)境質(zhì)量標(biāo)準(zhǔn),即使氮肥流失總量相同,其灰水足跡的結(jié)果也不同。本文將高集約化農(nóng)業(yè)生產(chǎn)中活性氮排放與非點源污染物納入到農(nóng)業(yè)土地利用環(huán)境效率評價體系中,采用氮足跡、灰水足跡理論與SBM-Undesirable模型的窗式分析方法相結(jié)合,能夠更為準(zhǔn)確地評估熱點元素和關(guān)鍵環(huán)境問題下的集約農(nóng)業(yè)土地利用環(huán)境效率,繼而為降低區(qū)域農(nóng)業(yè)生產(chǎn)的環(huán)境風(fēng)險提供了具體的改進方向與途徑。

    本文應(yīng)用SBM-Undesirable模型的窗式分析方法,很好地解決了非期望產(chǎn)出存在下的環(huán)境效率的評價問題,能夠科學(xué)地反映各個決策單元環(huán)境效率的真實水平。然而,農(nóng)業(yè)土地利用系統(tǒng)是一個巨復(fù)雜的人與自然耦合系統(tǒng),部分投入、期望產(chǎn)出和非期望產(chǎn)出指標(biāo)之間具有不可分離特性,本文鑒于條件所限,未充分考慮這一點。因此,下一步可在深入分析農(nóng)業(yè)土地利用過程中投入產(chǎn)出相互關(guān)系的基礎(chǔ)上,構(gòu)建不可分離變量的非期望SBM效率評價方法,以期更為精準(zhǔn)地反映農(nóng)業(yè)生產(chǎn)的實際土地利用過程。此外,應(yīng)指出的是影響農(nóng)業(yè)土地利用環(huán)境效率的因素眾多,涉及指標(biāo)體系、計算與評價方法等諸多方面。本文根據(jù)研究區(qū)的特點及數(shù)據(jù)的可獲得性,僅選取了耕地、勞動力與化肥等5個投入產(chǎn)出指標(biāo),其中針對環(huán)境風(fēng)險的問題僅考慮了活性氮排放與非點源污染。今后,可將氮足跡、灰水足跡與生態(tài)足跡等環(huán)境類足跡整合然后與SBM相關(guān)模型結(jié)合,可綜合評估集約農(nóng)業(yè)投入減量化與環(huán)境風(fēng)險降低潛勢,進而有助于提出全面科學(xué)的環(huán)境效率改善的政策與建議,同時對構(gòu)建投入減量化與低環(huán)境風(fēng)險的集約農(nóng)業(yè)可持續(xù)土地利用模式具有重要的理論和現(xiàn)實意義。

    5 結(jié) 論

    本文以上海市青浦區(qū)為案例區(qū),將氮足跡、灰水足跡理論與SBM-Undesirable模型的窗式分析方法相結(jié)合,分析集約農(nóng)業(yè)土地利用過程對環(huán)境的負面影響,系統(tǒng)解析環(huán)境效率的時空分異特征,進而綜合評估了投入減量化趨勢與環(huán)境風(fēng)險降低潛勢。主要結(jié)論如下:

    1)青浦區(qū)污染氮足跡與非點源污染灰水足跡存在顯著的時空差異性,總體呈西高東低的態(tài)勢。2006-2013年農(nóng)業(yè)污染氮足跡與非點源污染灰水足跡年均值分別為3 092.44 t、144 969.00萬m3,農(nóng)業(yè)生產(chǎn)導(dǎo)致的水質(zhì)惡化引起生態(tài)用水的增大,是上海青浦區(qū)水資源出現(xiàn)水質(zhì)性缺水的重要原因之一。

    2)在考慮活性氮排放和非點源污染約束下,2006-2013年青浦區(qū)農(nóng)業(yè)土地利用環(huán)境效率相對較低,2006-2013年均值僅為0.669。

    3)2006-2013年青浦區(qū)年均勞動力、肥料和機械動力的潛在減量比例較高,8 a間潛在年均縮減總量分別為 8 104人,4 501.59 t,27 928.44 kW;非點源污染灰水足跡的潛在減排比例高于污染氮足跡,潛在減排總量年均分別為52 046.88萬m3和381.04 t。

    4)花香橋街道具有最大的潛在投入減量化與環(huán)境風(fēng)險降低比例。白鶴鎮(zhèn)、練塘鎮(zhèn)等具有較大的潛在縮減規(guī)模,應(yīng)成為青浦區(qū)農(nóng)業(yè)投入減量化和污染物減排總量控制的重點區(qū)域。

    [1] Chatzimpiros P, Barles S. Nitrogen food-print: N use related to meat and dairy consumption in France[J]. Biogeosciences, 2013, 10: 471-481.

    [2] 中國科學(xué)院地理科學(xué)與資源研究所世界自然基金會. 中國生態(tài)足跡報告2012[R]. 世界自然基金會(WWF)與中國科學(xué)院地理科學(xué)與資源研究所,2012.

    [3] Pang Jiaxing, Chen Xingpeng, Zhang Zilong, et al. Measuring eco-efficiency of agriculture in China[J]. Sustainability, 2016, 8(4): 398.

    [4] 張可,豐景春. 強可處置性視角下中國農(nóng)業(yè)環(huán)境效率測度及其動態(tài)演進[J]. 中國人口&資源與環(huán)境,2016,26(1):140-149.

    Zhang Ke, Feng Jingchun. Agricultural environmental efficiency and its dynamic evolution of China from the perspective of pollution's strong disposability[J]. China Population, Resources and Environment, 2016, 26(1): 140-149. (in Chinese with English abstract)

    [5] 蔡雨君. 中國農(nóng)墾區(qū)農(nóng)業(yè)環(huán)境效率及其影響因素分析[D]. 大連:東北財經(jīng)大學(xué),2015.

    Cai Yujun. The Research of Agricultural Environmental Efficiency and Impact Factors in Reclamation Area of China[D]. Dalian: Dongbei University of Finance and Economics, 2015. (in Chinese with English abstract)

    [6] 梁勁銳,席小瑾,袁宇. 基于SBM模型的農(nóng)業(yè)環(huán)境技術(shù)效率分析—以山東省為例[J]. 四川農(nóng)業(yè)大學(xué)學(xué)報,2015,33(3):338-344.

    Liang Jinrui, Xi Xiaojin, Yuan Yu. The analysis of agricultural environmental technical efficiency based on SBM model: Taking Shandong province as an example[J]. Journal of Sichuan Agricultural University, 2015, 33(3): 338-344. (in Chinese with English abstract)

    [7] 楊志海,王雅鵬. 中國農(nóng)業(yè)環(huán)境技術(shù)效率與污染減排成本分析[J]. 統(tǒng)計與信息論壇,2015,30(7):71-76.

    Yang Zhihai, Wang Yapeng. Environmentally technical efficiency and pollution abatement cost in China's agricultural sector[J]. Statistics & Information Forum, 2015, 30(7): 71-76. (in Chinese with English abstract)

    [8] 梁流濤. 考慮“非意欲”產(chǎn)出的農(nóng)業(yè)土地生產(chǎn)效率評價及其時空特征分析[J]. 資源科學(xué),2012,34(12):2249-2255.

    Liang Liutao. Agricultural land efficiency and temporal-spatial characterstics: Desired output and undesired pollution emissions[J]. Resources Science, 2012, 34(12): 2249-2255. (in Chinese with English abstract)

    [9] Aldanondo-Ochoa A M, Casasnovas-Oliva V L, Arandia- Miura A. Environmental efficiency and the impact of regulation in dryland organic vine production[J]. Land Use Policy, 2014, 36: 275-284.

    [10] Marthin N, Iain F, Ali Q, et al. Environmentally adjusted productivity measurement: An Australian case study[J]. Journal of Environmental Management, 2007, 85(2): 350-362.

    [11] Tone K. Dealing with undesirable outputs in DEA: A slacks-based measure (SBM) approach[C]//Tokyo: The Operations Research Society of Japan, 2004.

    [12] Kuo H F, Chen H L, Tsou K W, et al. Analysis of farming environmental efficiency using a dEA model with undesirable outputs[J]. Apcbee Procedia, 2014, 10(10): 154-158.

    [13] 封永剛,彭玨,鄧宗兵,等. 面源污染、碳排放雙重視角下中國耕地利用效率的時空分異[J]. 中國人口&資源與環(huán)境,2015,25(8):18-25.

    Feng Yonggang, Peng Jue, Deng Zongbing, et al. Spatial-temporal variation of cultivated land’s utilization efficiency in China based on the dual perspective of non-point source pollution and carbon emission[J]. China Population Resources and Environment, 2015, 25(8): 18-25. (in Chinese with English abstract)

    [14] 朱聰. 中國農(nóng)業(yè)環(huán)境效率研究-基于環(huán)境的三階段DEA[D]. 蚌埠:安徽財經(jīng)大學(xué),2015.

    Zhu Cong. The Analysis of Agricultural Environmental Efficiency in China Based on a Environmental Three-stage DEA[D]. Bengbu: Anhui University of Finance & Economics, 2015. (in Chinese with English abstract)

    [15] 田偉,楊璐嘉,姜靜. 低碳視角下中國農(nóng)業(yè)環(huán)境效率的測算與分析—基于非期望產(chǎn)出的SBM模型[J]. 中國農(nóng)村觀察,2014(5):59-71.

    [16] 潘丹,應(yīng)瑞瑤. 中國農(nóng)業(yè)生態(tài)效率評價方法與實證—基于非期望產(chǎn)出的SBM模型分析[J]. 生態(tài)學(xué)報,2013,33(12):3837-3845. Pan Dan, Ying Ruiyao. Agricultural eco-efficiency evaluation in China based on SBM model[J]. Ecologica Sinica, 2013, 33(12): 3837-3845. (in Chinese with English abstract)

    [17] 付永虎,劉黎明,袁承程. 農(nóng)業(yè)土地利用系統(tǒng)氮足跡與灰水足跡綜合評價[J]. 農(nóng)業(yè)工程學(xué)報,2016,32(增刊1):312-319.

    Fu Yonghu, Liu Liming, Yuan Chengcheng. Comprehensive evaluation for nitrogen footprint and grey water footprint of agricultural land use system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(Supp.1): 312-319. (in Chinese with English abstract)

    [18] 付永虎,劉黎明,起曉星,等. 基于灰水足跡的洞庭湖區(qū)糧食生產(chǎn)環(huán)境效應(yīng)評價[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(10):152-160.

    Fu Yonghu, Liu Liming, Qi Xiaoxing, et al. Environmental effects evaluation for grain production based on grey water footprint in Dongting Lake area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(10): 152-160. (in Chinese with English abstract)

    [19] 王新新. 基于農(nóng)田氮磷平衡的太湖流域環(huán)境風(fēng)險評估[D]. 北京:中國農(nóng)業(yè)大學(xué),2015.

    Wang Xinxin. The Evaluation of Environmental Risk in Tai Lake Basin Based on Cropland Nutrients Balance[D]. Beijing: China Agricultural University, 2015. (in Chinese with English abstract)

    [20] 王興仁,張福鎖,張衛(wèi)峰,等. 中國農(nóng)化服務(wù)肥料與施肥手冊[M]. 北京:中國農(nóng)業(yè)出版社,2013.

    [21] 國務(wù)院第一次全國污染源普查領(lǐng)導(dǎo)小組辦公室. 第一次全國污染源普查-農(nóng)業(yè)污染源肥料流失系數(shù)手冊[R]. 2009.

    [22] 遆超普. 不同空間尺度區(qū)域氮素收支[D]. 南京:南京農(nóng)業(yè)大學(xué),2011.

    Ti Chaopu. Regional Nitrogen Budget at Different Spatial Scales[D]. Nanjing: Nanjing Agricultural University, 2011. (in Chinese with English abstract)

    [23] 谷保靜. 人類-自然耦合系統(tǒng)氮循環(huán)研究—中國案例[D]. 杭州:浙江大學(xué),2011.

    Gu Baojing. Nitrogen Cycle of Coupled Human and Natural System-A Case Study of China[D]. Hangzhou: Zhe Jiang University, 2011. (in Chinese with English abstract)

    [24] 國家發(fā)展和改革委員應(yīng)對氣候變化司. 2005中國溫室氣體清單研究[M]. 北京:中國環(huán)境科學(xué)出版社,2014.

    [25] 周濤,王云鵬,王芳,等. 廣東省農(nóng)業(yè)氮足跡分析[J]. 中國環(huán)境科學(xué),2014,34(9):2430-2438.

    Zhou Tao, Wang Yunpeng, Wang Fang, et al. Analysis of the nitrogen footprint of agriculture in Guangdong[J]. China Environmental Science, 2014, 34(9): 2430-2438. (in Chinese with English abstract)

    [26] 郭春霞. 平原河網(wǎng)地區(qū)農(nóng)村面源污染重點源和區(qū)的識別篩選——以上海青浦區(qū)為例[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2011(8):1652-1659.

    Guo Chunxia. Identifications of critical sources and areas of rural non-point source pollution in the plain river network area: A case study of Qingpu, Shanghai, China[J]. Journal of Agro-Environment Science, 2011, 30(8): 1652-1659. (in Chinese with English abstract)

    [27] 楊清可,段學(xué)軍,葉磊,等. 基于SBM-Undesirable模型的城市土地利用效率評價——以長三角地區(qū)16城市為例[J]. 資源科學(xué),2014(4):712-721.

    Yang Qingke, Duan Xuejun, Ye Lei, et al. Efficiecy evaluation of city land utilization in the Yangtze river delta using a SBM-undesirable model[J]. Resources Science, 2014, 36(4): 712-721. (in Chinese with English abstract)

    [28] Tone K. A slacks-based measure of efficiency in data envelopment analysis[J]. European Journal of Operational Research, 2001, 130(3): 498-509.

    [29] Cooper W W, Seiford L M, Tone K. Data Envelopment Analysis (Second edition)[M]. Boston: Kluwer academic publishers, 2007.

    [30] Charnes A, Clark C T, Cooper W W A. Developmental study of data envelopment analysis in measuring the efficiency of maintenance units in the US air forces[J]. Annals of Operations Research, 1984, 2(1): 95-112.

    [31] Halkos G E, Tzeremes N G. Exploring the existence of Kuznets curve in countries’ environmental efficiency using DEA window analysis[J]. Ecological Economics, 2009, 68(7): 2168-2176.

    [32] 涂正革. 環(huán)境、資源與工業(yè)增長的協(xié)調(diào)性[J]. 經(jīng)濟研究,2008(2):93-105.

    Tu Zhengge. The coordination of industrial growth with environment and resource[J]. Economic Research Journal, 2008(2): 93-105. (in Chinese with English abstract)

    [33] 周賽軍,鄧明君,羅文兵. 氮足跡核算研究最新進展[J]. 浙江農(nóng)業(yè)科學(xué),2015,56(6):915-920.

    [34] 秦樹平,胡春勝,張玉銘,等. 氮足跡研究進展[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報,2011,19(2):462-467.

    Qin Shuping, Hu Chunsheng, Zhang Yuming, et al. Advances in nitrogen footprint research[J]. Chinese Journal of Eco-Agriculture, 2011, 19(2): 462-467. (in Chinese with English abstract)

    [35] Hoekstra A Y, Chapagain A K, Aldaya M M, 等著;劉俊國,曾昭,趙乾斌,等,譯. 水足跡評價手冊[M]. 北京:科學(xué)出版社,2012.

    Spatial-temporal variations of dematerialization of inputs and environmental risk reduction in intensive agricultural region

    Fu Yonghu1,2, Liu Liming1※, Wang Jiasheng1, Ye Jinwei1, Guo Yun1

    (1.100193,; 2.222005,)

    Evaluation of the agricultural land use environmental efficiency is an essential issue for building sustainable land use patterns with fewer agricultural inputs and lower environmental risk. In this paper, 1) nitrogen footprint (N footprint) and grey water footprint (GWF) were used to evaluate environmental effects (reactive nitrogen emission and non-point source pollution) of agricultural land use system at the regional scale; 2) the window type slacks-based measure (SBM) model based on undesirable outputs (pollution N footprint and non-point resource pollution GWF) was applied to calculate environmental efficiency of the agricultural land use system; 3) inputs and environmental risk reduction potentials were computed and evaluated based on the dual perspective of reactive nitrogen emission and non-point source pollution in an intensive agricultural area. The Qingpu District was selected as a study case for this approach. The results mainly show that: 1) There were significant spatial-temporal differences in using pollution N footprint and non-point source pollution GWF within Qingpu District. The pollution N footprint of agricultural land use system was between 2 805.59 t and 3 505.68 t in the period of 2006-2013. The annual total grey water volume for agricultural non-point source pollution ranged from 134 026.58×104m3to 152 833.70×104m3during the years 2006-2013, and its average value was 144 969.00×104m3. The negative effect of water environment under agricultural land use process was assessed by the indicator of water environmental pressure (WEP) at county level. The average value of WEP of agricultural land use system was 5.74. The high value of WEP suggested that the negative environmental effect of agricultural production increased. The spatial distribution of pollution N footprint and non-point source GWF in Qingpu District was higher value in the west but lower in the east in general. 2) According to the results of SBM-undesirable model, overall environmental efficiency of agricultural land use in Qingpu District remained at a low level, and its average efficiency was 0.669 in the period of 2006-2013. 3) The models of potential percentage reduction and potential reduction scale were built for evaluation of dematerialization and environmental risk reduction respectively based on the dual perspective of reactive nitrogen emission and non-point source pollution. The annual average percentage reduction of labor, fertilizer, mechanical power, pollution N footprint and non-point source pollution GWF were 32.21%, 25.70%, 38.21%, 12.32% and 32.18%, respectively. The average reductions of labor, fertilizer, mechanical power, pollution N footprint and non-point source pollution GWF were 8 104 people, 4 501.59 t, 27 928.44 kW, 381.04 t and 52 046.88×104m3during the years 2006-2013. Xiang Huaqiao had the highest percentage of environmental risk reduction among the 11 towns, while Baihe and Liantang had a larger scale of emission reductions at the same time. Therefore, Baihe and Liantang should be the key areas for agricultural resources conservation and pollutant emissions control within Qingpu District. The findings of this study can provide scientific supports for regional agricultural planning and policy-making.

    land use; pollution; agriculture; nitrogen footprint; grey water footprint; environmental efficiency; low input; environmental risk

    10.11975/j.issn.1002-6819.2017.02.037

    F301.2; X82

    A

    1002-6819(2017)-02-0266-10

    2016-06-19

    2016-12-17

    國家自然科學(xué)基金重點項目(41130526);淮海工學(xué)院科研基金資助項目(S201605)

    付永虎,博士生,研究方向為土地利用系統(tǒng)研究。北京 中國農(nóng)業(yè)大學(xué)土地資源管理系,100193。Email:huhu_0902@163.com

    劉黎明,教授,博士,博士生導(dǎo)師,主要研究方向為土地資源可持續(xù)利用與景觀規(guī)劃。北京 中國農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,100193。Email:liulm@cau.edu.cn

    付永虎,劉黎明,王加升,葉津煒,郭 赟. 高集約化農(nóng)區(qū)投入減量化與環(huán)境風(fēng)險降低潛勢的時空分異特征[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(2):266-275. doi:10.11975/j.issn.1002-6819.2017.02.037 http://www.tcsae.org

    Fu Yonghu, Liu Liming, Wang Jiasheng, Ye Jinwei, Guo Yun. Spatial-temporal variations of dematerialization of inputs and environmental risk reduction in intensive agricultural region[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(2): 266-275. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.02.037 http://www.tcsae.org

    猜你喜歡
    灰水青浦區(qū)足跡
    7.2.4全景聲影視廳 上海市青浦區(qū)御湖別墅香水灣
    在這里,與100個世界相遇——上海市青浦區(qū)教師進修學(xué)院附屬小學(xué)
    少先隊活動(2020年9期)2020-10-28 08:04:12
    中國灰水足跡時空格局演變及其脫鉤關(guān)系研究
    上海市青浦區(qū)大用水戶管理的應(yīng)用
    成長足跡
    足跡
    心聲歌刊(2018年6期)2018-01-24 00:56:12
    灰水回用的可行性分析
    山東化工(2017年23期)2017-05-21 03:21:23
    具有景觀效益的居民住宅區(qū)灰水回用系統(tǒng)研究
    足跡
    社會與公益(2016年2期)2016-04-13 02:49:06
    氣化裝置高壓灰水調(diào)節(jié)閥的損蝕原因分析及改造
    亚洲av不卡在线观看| 欧美一区二区精品小视频在线| 午夜福利欧美成人| 国产一区二区三区在线臀色熟女| 国产精品嫩草影院av在线观看 | 桃红色精品国产亚洲av| 婷婷丁香在线五月| 欧美激情久久久久久爽电影| 精品人妻视频免费看| 欧美日韩中文字幕国产精品一区二区三区| 精品无人区乱码1区二区| 国产探花在线观看一区二区| 一区二区三区高清视频在线| 床上黄色一级片| 久久天躁狠狠躁夜夜2o2o| 99久久精品国产亚洲精品| 99热精品在线国产| 亚洲一区二区三区色噜噜| 午夜福利在线在线| 波多野结衣高清无吗| 两人在一起打扑克的视频| 欧美日韩乱码在线| 国产91精品成人一区二区三区| 大型黄色视频在线免费观看| 在线天堂最新版资源| 日本三级黄在线观看| 亚洲欧美日韩高清专用| 国产精品嫩草影院av在线观看 | 动漫黄色视频在线观看| 不卡一级毛片| 久久精品人妻少妇| 国产精品野战在线观看| 欧美一区二区精品小视频在线| 婷婷精品国产亚洲av| 国产精品人妻久久久久久| 小蜜桃在线观看免费完整版高清| 成人国产综合亚洲| 中文字幕av成人在线电影| 婷婷精品国产亚洲av| 桃色一区二区三区在线观看| 国产精品亚洲av一区麻豆| 亚洲男人的天堂狠狠| 精品午夜福利视频在线观看一区| 丰满人妻熟妇乱又伦精品不卡| 日韩大尺度精品在线看网址| a级一级毛片免费在线观看| 91狼人影院| 国产老妇女一区| 欧美又色又爽又黄视频| 色吧在线观看| 久久精品人妻少妇| 淫妇啪啪啪对白视频| 国产人妻一区二区三区在| 97热精品久久久久久| 性插视频无遮挡在线免费观看| 深爱激情五月婷婷| 88av欧美| 色在线成人网| 91av网一区二区| 亚洲18禁久久av| 亚洲五月婷婷丁香| 国产综合懂色| 久久久久久九九精品二区国产| 国产av不卡久久| 一级作爱视频免费观看| 9191精品国产免费久久| 国产免费一级a男人的天堂| 深夜a级毛片| 国产精品久久久久久精品电影| 亚洲乱码一区二区免费版| 天天一区二区日本电影三级| 波多野结衣巨乳人妻| 日韩欧美国产在线观看| 久久精品国产亚洲av天美| 国产aⅴ精品一区二区三区波| 亚洲人成伊人成综合网2020| 99热这里只有是精品在线观看 | 精品福利观看| 日本成人三级电影网站| 99国产极品粉嫩在线观看| 亚洲综合色惰| 国产欧美日韩精品亚洲av| 一个人看视频在线观看www免费| 国内精品久久久久久久电影| 狂野欧美白嫩少妇大欣赏| 又爽又黄无遮挡网站| 日韩亚洲欧美综合| 国产成人aa在线观看| 精品日产1卡2卡| 一本精品99久久精品77| a级一级毛片免费在线观看| 成人高潮视频无遮挡免费网站| 欧美区成人在线视频| 日本免费一区二区三区高清不卡| 久久国产乱子伦精品免费另类| 91麻豆精品激情在线观看国产| 露出奶头的视频| 91在线精品国自产拍蜜月| 99riav亚洲国产免费| 亚洲成a人片在线一区二区| 别揉我奶头 嗯啊视频| 欧美乱色亚洲激情| 国产乱人视频| 久久久久久久久大av| 欧美一区二区亚洲| 啪啪无遮挡十八禁网站| 亚洲av成人av| 在线观看一区二区三区| 婷婷六月久久综合丁香| 综合色av麻豆| 最近视频中文字幕2019在线8| 51午夜福利影视在线观看| 免费观看的影片在线观看| 无人区码免费观看不卡| 欧美成人免费av一区二区三区| 亚洲三级黄色毛片| 在线十欧美十亚洲十日本专区| 1000部很黄的大片| 精品国产三级普通话版| 久久精品综合一区二区三区| 嫩草影院入口| 非洲黑人性xxxx精品又粗又长| 精品午夜福利视频在线观看一区| 波多野结衣巨乳人妻| 成人特级黄色片久久久久久久| 国产精品av视频在线免费观看| 亚洲七黄色美女视频| 精品无人区乱码1区二区| 精品一区二区三区av网在线观看| 亚洲天堂国产精品一区在线| 亚洲成av人片免费观看| 亚洲精品影视一区二区三区av| 桃红色精品国产亚洲av| 禁无遮挡网站| 亚洲精品456在线播放app | 男人狂女人下面高潮的视频| 久久久精品大字幕| 成人高潮视频无遮挡免费网站| 12—13女人毛片做爰片一| av欧美777| 欧美日韩中文字幕国产精品一区二区三区| 久久久久国内视频| 国产单亲对白刺激| 国产精品野战在线观看| 国产视频一区二区在线看| 亚洲片人在线观看| 两个人视频免费观看高清| 久久精品国产99精品国产亚洲性色| 99久久99久久久精品蜜桃| 国产乱人视频| 精品人妻1区二区| 国产高清视频在线播放一区| 欧美+亚洲+日韩+国产| 亚洲专区中文字幕在线| 国产乱人视频| 久久精品久久久久久噜噜老黄 | av天堂中文字幕网| 国产一区二区在线观看日韩| 最近中文字幕高清免费大全6 | 午夜福利在线观看吧| 搡老妇女老女人老熟妇| 国产成人福利小说| 真人一进一出gif抽搐免费| 国内久久婷婷六月综合欲色啪| 亚洲国产日韩欧美精品在线观看| 国内精品美女久久久久久| 久久人人爽人人爽人人片va | 欧美日韩中文字幕国产精品一区二区三区| 亚洲精品久久国产高清桃花| 性欧美人与动物交配| 国内精品久久久久久久电影| 国产v大片淫在线免费观看| 成熟少妇高潮喷水视频| 亚洲人成网站高清观看| 欧美日韩国产亚洲二区| 亚洲欧美日韩东京热| 91午夜精品亚洲一区二区三区 | 日韩欧美国产在线观看| 精品一区二区三区av网在线观看| 一级av片app| 又黄又爽又刺激的免费视频.| 国内精品一区二区在线观看| 成人毛片a级毛片在线播放| 91麻豆精品激情在线观看国产| 69av精品久久久久久| 亚洲人成伊人成综合网2020| 久久国产精品人妻蜜桃| 白带黄色成豆腐渣| 国产男靠女视频免费网站| 国产成年人精品一区二区| 欧美激情在线99| 成年免费大片在线观看| 看片在线看免费视频| 日韩人妻高清精品专区| 亚洲精品粉嫩美女一区| 国产成人aa在线观看| 久久久久免费精品人妻一区二区| 最好的美女福利视频网| 欧美性感艳星| 一进一出抽搐动态| 极品教师在线视频| 久久久国产成人精品二区| 亚洲精品乱码久久久v下载方式| 国产在视频线在精品| 亚洲午夜理论影院| 欧美激情久久久久久爽电影| 麻豆久久精品国产亚洲av| 国产伦精品一区二区三区视频9| 我要搜黄色片| 中文字幕精品亚洲无线码一区| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久久毛片| 欧美在线黄色| 偷拍熟女少妇极品色| 99热6这里只有精品| 怎么达到女性高潮| 亚洲电影在线观看av| 午夜a级毛片| 久久6这里有精品| 人妻久久中文字幕网| 国产探花极品一区二区| 日本与韩国留学比较| 免费看美女性在线毛片视频| 女人被狂操c到高潮| 欧美+日韩+精品| 日韩欧美精品v在线| 俄罗斯特黄特色一大片| 一级毛片久久久久久久久女| 欧美xxxx黑人xx丫x性爽| 变态另类丝袜制服| 欧美成人a在线观看| 亚洲第一欧美日韩一区二区三区| 精品国内亚洲2022精品成人| 日本黄大片高清| 精品午夜福利在线看| 午夜精品在线福利| a级一级毛片免费在线观看| 波多野结衣高清无吗| 国内毛片毛片毛片毛片毛片| 蜜桃亚洲精品一区二区三区| 网址你懂的国产日韩在线| 亚洲片人在线观看| 一个人观看的视频www高清免费观看| 1024手机看黄色片| 欧美成人性av电影在线观看| 免费人成视频x8x8入口观看| aaaaa片日本免费| 午夜精品一区二区三区免费看| 国产免费男女视频| 亚洲成人中文字幕在线播放| а√天堂www在线а√下载| 亚洲avbb在线观看| 国产老妇女一区| 日本精品一区二区三区蜜桃| 亚洲国产色片| 老女人水多毛片| 国产白丝娇喘喷水9色精品| 亚洲一区二区三区不卡视频| 国产高清视频在线播放一区| 色噜噜av男人的天堂激情| 国产亚洲欧美在线一区二区| 一个人免费在线观看电影| 麻豆av噜噜一区二区三区| 国产精品亚洲美女久久久| 亚洲成人中文字幕在线播放| 丝袜美腿在线中文| 精品人妻一区二区三区麻豆 | 成人毛片a级毛片在线播放| 久久中文看片网| 18禁黄网站禁片午夜丰满| 亚洲成av人片在线播放无| 婷婷色综合大香蕉| 亚洲精品久久国产高清桃花| 人妻丰满熟妇av一区二区三区| av天堂在线播放| 国产黄色小视频在线观看| 蜜桃亚洲精品一区二区三区| 少妇熟女aⅴ在线视频| 少妇被粗大猛烈的视频| 丝袜美腿在线中文| 亚洲欧美激情综合另类| 国产aⅴ精品一区二区三区波| 少妇熟女aⅴ在线视频| 3wmmmm亚洲av在线观看| 51午夜福利影视在线观看| 99久久精品国产亚洲精品| 国产在线男女| 亚洲欧美日韩无卡精品| 日韩av在线大香蕉| 蜜桃亚洲精品一区二区三区| 国产高清三级在线| 变态另类成人亚洲欧美熟女| 国产91精品成人一区二区三区| 天堂√8在线中文| 男女之事视频高清在线观看| 亚洲欧美激情综合另类| 757午夜福利合集在线观看| 亚洲欧美日韩无卡精品| 久久久久久久精品吃奶| 亚洲久久久久久中文字幕| 波野结衣二区三区在线| 国产成年人精品一区二区| 国产成人av教育| 午夜日韩欧美国产| 精华霜和精华液先用哪个| 欧美xxxx黑人xx丫x性爽| 免费在线观看日本一区| 精品人妻1区二区| 亚洲国产精品成人综合色| 午夜精品一区二区三区免费看| 国产野战对白在线观看| 嫩草影院入口| 脱女人内裤的视频| 精华霜和精华液先用哪个| 青草久久国产| 91九色精品人成在线观看| 丰满的人妻完整版| 精品熟女少妇八av免费久了| 内射极品少妇av片p| 欧美黑人欧美精品刺激| 桃色一区二区三区在线观看| h日本视频在线播放| 久久久久久久久大av| 日本a在线网址| 国产麻豆成人av免费视频| 精品久久久久久久久久久久久| 国产白丝娇喘喷水9色精品| 久久精品久久久久久噜噜老黄 | 欧美不卡视频在线免费观看| 国产私拍福利视频在线观看| 午夜激情欧美在线| 成年版毛片免费区| 久久精品影院6| 成年女人毛片免费观看观看9| 国产一级毛片七仙女欲春2| av在线天堂中文字幕| 女人被狂操c到高潮| 在线a可以看的网站| 在线观看一区二区三区| 一级毛片久久久久久久久女| 亚洲欧美激情综合另类| 很黄的视频免费| 亚洲三级黄色毛片| 波野结衣二区三区在线| 国产精品免费一区二区三区在线| 成人av在线播放网站| 欧美乱色亚洲激情| 九九久久精品国产亚洲av麻豆| x7x7x7水蜜桃| 国产精品人妻久久久久久| 欧美精品啪啪一区二区三区| 波多野结衣高清无吗| АⅤ资源中文在线天堂| 亚洲最大成人手机在线| 国产在线男女| 午夜免费成人在线视频| 日本免费a在线| 免费黄网站久久成人精品 | 女人十人毛片免费观看3o分钟| 精品一区二区三区av网在线观看| 色在线成人网| 国产亚洲精品久久久com| 久久久久久久久久成人| 亚洲18禁久久av| 亚洲不卡免费看| 久久热精品热| 午夜福利在线观看吧| 99久久精品一区二区三区| 精品不卡国产一区二区三区| av在线蜜桃| 久久中文看片网| 午夜老司机福利剧场| 99热精品在线国产| 国产高潮美女av| 久久精品综合一区二区三区| 亚洲成a人片在线一区二区| 日本成人三级电影网站| 黄色视频,在线免费观看| 国产在视频线在精品| 村上凉子中文字幕在线| 久久久色成人| av福利片在线观看| 最近视频中文字幕2019在线8| 国产色婷婷99| 中文字幕人成人乱码亚洲影| 无人区码免费观看不卡| 搞女人的毛片| 88av欧美| 狂野欧美白嫩少妇大欣赏| 99热这里只有精品一区| 亚洲专区国产一区二区| 国产私拍福利视频在线观看| 老司机午夜十八禁免费视频| 看十八女毛片水多多多| 精品久久久久久久久av| 色综合欧美亚洲国产小说| 夜夜看夜夜爽夜夜摸| 免费观看人在逋| 99热精品在线国产| 亚洲综合色惰| 性欧美人与动物交配| 久久久久久九九精品二区国产| 成人无遮挡网站| 特级一级黄色大片| 在线a可以看的网站| 长腿黑丝高跟| 最新在线观看一区二区三区| 亚洲最大成人中文| 久久婷婷人人爽人人干人人爱| 麻豆成人午夜福利视频| 婷婷亚洲欧美| 亚洲av日韩精品久久久久久密| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 禁无遮挡网站| 欧美绝顶高潮抽搐喷水| 麻豆国产97在线/欧美| 亚洲天堂国产精品一区在线| 亚洲精品成人久久久久久| 可以在线观看的亚洲视频| 久久香蕉精品热| 欧美zozozo另类| 国产成+人综合+亚洲专区| 深夜精品福利| 亚洲第一电影网av| 色在线成人网| av女优亚洲男人天堂| 日韩精品中文字幕看吧| 他把我摸到了高潮在线观看| 午夜福利视频1000在线观看| 小蜜桃在线观看免费完整版高清| 国产av在哪里看| 午夜福利视频1000在线观看| 伊人久久精品亚洲午夜| 国产一区二区亚洲精品在线观看| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费成人在线视频| 日日摸夜夜添夜夜添av毛片 | 成人av在线播放网站| 99久久久亚洲精品蜜臀av| 日韩 亚洲 欧美在线| 婷婷精品国产亚洲av在线| 久久草成人影院| 久久婷婷人人爽人人干人人爱| 十八禁国产超污无遮挡网站| 欧美日韩中文字幕国产精品一区二区三区| 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区免费毛片| 一区二区三区高清视频在线| 女人被狂操c到高潮| 热99在线观看视频| 亚洲在线自拍视频| 亚洲无线在线观看| 国内精品一区二区在线观看| 亚洲五月婷婷丁香| 神马国产精品三级电影在线观看| 欧美日韩亚洲国产一区二区在线观看| 我的女老师完整版在线观看| 99精品在免费线老司机午夜| 观看美女的网站| 99久久九九国产精品国产免费| 嫁个100分男人电影在线观看| 1000部很黄的大片| 性欧美人与动物交配| 国产主播在线观看一区二区| 免费在线观看日本一区| 观看美女的网站| 国产精品1区2区在线观看.| 男人和女人高潮做爰伦理| 此物有八面人人有两片| 老司机福利观看| 日韩欧美三级三区| 国产精品不卡视频一区二区 | 亚洲成人免费电影在线观看| 亚洲精华国产精华精| 亚洲精品一卡2卡三卡4卡5卡| 国产欧美日韩一区二区三| 免费观看人在逋| 国产精品亚洲美女久久久| 国产淫片久久久久久久久 | 精品一区二区三区人妻视频| 97热精品久久久久久| 99久久九九国产精品国产免费| 69人妻影院| 成人永久免费在线观看视频| 国产精品久久久久久人妻精品电影| 亚洲国产高清在线一区二区三| 无人区码免费观看不卡| 亚洲片人在线观看| 美女被艹到高潮喷水动态| 欧美日韩乱码在线| 亚洲电影在线观看av| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 窝窝影院91人妻| 国产伦人伦偷精品视频| 亚洲第一区二区三区不卡| 亚洲欧美激情综合另类| 色综合亚洲欧美另类图片| 看片在线看免费视频| 欧美zozozo另类| 国产精品亚洲一级av第二区| 午夜激情福利司机影院| 亚洲成人中文字幕在线播放| 天堂网av新在线| 超碰av人人做人人爽久久| 露出奶头的视频| 成年版毛片免费区| 国产精品美女特级片免费视频播放器| 色噜噜av男人的天堂激情| 国产色婷婷99| 国产精品精品国产色婷婷| 一级黄色大片毛片| 老司机午夜十八禁免费视频| 草草在线视频免费看| 99热精品在线国产| 欧美zozozo另类| 国产高潮美女av| 夜夜夜夜夜久久久久| 色综合亚洲欧美另类图片| 看黄色毛片网站| 亚洲五月天丁香| 日本免费一区二区三区高清不卡| 免费黄网站久久成人精品 | 天天一区二区日本电影三级| 日韩亚洲欧美综合| 国产三级在线视频| 日韩欧美国产在线观看| 欧美性感艳星| 麻豆成人av在线观看| x7x7x7水蜜桃| 国产私拍福利视频在线观看| 免费人成在线观看视频色| 免费大片18禁| 999久久久精品免费观看国产| 成人鲁丝片一二三区免费| 国产精品一及| 日韩大尺度精品在线看网址| 一二三四社区在线视频社区8| 国产久久久一区二区三区| 国产精品野战在线观看| 精品乱码久久久久久99久播| 级片在线观看| 亚洲av.av天堂| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美一区二区三区在线观看| 久久人妻av系列| 亚洲,欧美精品.| 99热精品在线国产| 村上凉子中文字幕在线| 欧美性猛交黑人性爽| 窝窝影院91人妻| 女生性感内裤真人,穿戴方法视频| 人妻丰满熟妇av一区二区三区| 精品午夜福利在线看| 99久久精品一区二区三区| 露出奶头的视频| 内射极品少妇av片p| 啪啪无遮挡十八禁网站| 永久网站在线| 老司机福利观看| 国产精品日韩av在线免费观看| 免费av不卡在线播放| 久久久久性生活片| 久久久久久久久久黄片| 亚洲欧美日韩东京热| 能在线免费观看的黄片| 黄色女人牲交| 91九色精品人成在线观看| 亚洲最大成人手机在线| 国产一区二区在线观看日韩| 精品久久久久久久人妻蜜臀av| 欧美xxxx黑人xx丫x性爽| 亚洲乱码一区二区免费版| 亚洲欧美日韩高清在线视频| 91久久精品国产一区二区成人| 亚洲 国产 在线| 国产大屁股一区二区在线视频| 一个人免费在线观看电影| 精品一区二区三区人妻视频| 免费黄网站久久成人精品 | 国产激情偷乱视频一区二区| 成人无遮挡网站| 日韩欧美精品v在线| 成人一区二区视频在线观看| 午夜福利在线观看吧| 国产午夜精品久久久久久一区二区三区 | 国产野战对白在线观看| 国内精品久久久久精免费| 久久亚洲真实| 久久人人精品亚洲av| 欧美最黄视频在线播放免费| 午夜福利在线观看吧| 日韩成人在线观看一区二区三区| 夜夜躁狠狠躁天天躁| 女人十人毛片免费观看3o分钟| 97碰自拍视频| 欧美极品一区二区三区四区| 久久久久久久久大av| 亚洲av中文字字幕乱码综合| 在线观看美女被高潮喷水网站 | 久久精品综合一区二区三区| 亚洲最大成人中文| 亚洲最大成人手机在线| 国产野战对白在线观看| 国产精品永久免费网站| 国产亚洲欧美在线一区二区| 伊人久久精品亚洲午夜| 久久性视频一级片| 国产午夜精品久久久久久一区二区三区 | 国产午夜福利久久久久久| 日韩 亚洲 欧美在线| 精品国产三级普通话版| 麻豆一二三区av精品| 色播亚洲综合网|