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

    水稻灌溉用水效率和要素稟賦對化肥面源污染的影響
    ——基于分位數(shù)回歸的分析

    2017-08-11 08:46:41于偉詠漆雁斌何悅鄧鑫
    農業(yè)環(huán)境科學學報 2017年7期
    關鍵詞:面源化肥用水

    于偉詠,漆雁斌,何悅,2,鄧鑫

    (1.四川農業(yè)大學經(jīng)濟學院,成都 611130;2.長江師范學院財經(jīng)學院,重慶涪陵 408000;3.四川農業(yè)大學管理學院,成都 611130)

    水稻灌溉用水效率和要素稟賦對化肥面源污染的影響
    ——基于分位數(shù)回歸的分析

    于偉詠1,漆雁斌1,何悅1,2,鄧鑫3

    (1.四川農業(yè)大學經(jīng)濟學院,成都 611130;2.長江師范學院財經(jīng)學院,重慶涪陵 408000;3.四川農業(yè)大學管理學院,成都 611130)

    通過測度長江上游4省(市)797份樣本戶水稻灌溉用水效率和化肥面源污染量,基于分位數(shù)回歸方法,深入分析灌溉用水效率、生產(chǎn)方式和農戶要素稟賦對化肥面源污染的影響。研究表明,總體上灌溉總效率均值為0.14,純技術效率均值為0.27,規(guī)模效率均值為0.55,水稻灌溉用水效率改進空間較大,規(guī)模化程度高的平原地區(qū)用水效率更高;灌溉用水效率與化肥施用量間存在顯著替代效應,兩者呈振幅增強的波動式發(fā)展;節(jié)水灌溉能夠降低化肥面源污染,同時具有保肥效果;水旱輪作模式、耕作方式、秸稈還田方式對化肥面源污染存在影響,但未全部通過顯著性檢驗;耕地細碎化、農業(yè)補貼和商品化率會增加化肥面源污染;年齡、受教育程度減輕化肥面源污染的臨界點較高,政治參與程度對化肥面源污染影響不顯著。

    化肥面源污染;灌溉用水效率;生產(chǎn)方式;要素稟賦;分位數(shù)回歸

    農業(yè)向人類提供食物、纖維和原材料的同時,也產(chǎn)生環(huán)境污染等有害物品[1]。其中化肥的過量施用和高流失率,通過農田地表徑流、壤中流、農田排水和地下水滲漏進入水體而形成地表和地下水環(huán)境污染,從而形成由化肥造成的農業(yè)面源污染(簡稱化肥農業(yè)面源污染)[2]。我國化肥有效利用率相對發(fā)達國家偏低,有研究表明目前氮肥利用率僅為30%~35%,磷肥和鉀肥分別為10%~20%和35%~50%,平均低于發(fā)達國家15%~20%[3]。過量施肥現(xiàn)象嚴重,特別是水稻高產(chǎn)地區(qū)同時也是高污染地區(qū)[5]。據(jù)測算1990—2008年間我國因化肥施用導致的TN排放從313.27萬t增加到408.88萬t,TP排放從16.66萬t增加到25.03萬t[6]。因此,防控化肥面源污染成為緊迫的現(xiàn)實問題,同時具有潛在經(jīng)濟價值,如整個江蘇省治理農業(yè)面源污染可實現(xiàn)25.06億元收益[7],但其公共屬性決定治理需要由政府負責。

    化肥面源污染發(fā)生一是由于負外部性、“公地悲劇”、治理成本高等特征導致的“市場失靈”和政府的“規(guī)制失靈”[8],制度環(huán)境產(chǎn)生了非預期效果[9];二是從農戶生產(chǎn)行為出發(fā),生存和發(fā)展壓力,農業(yè)經(jīng)營行為短視化,環(huán)保意識淡薄,缺乏公共服務支持等原因,集中導致化肥過量不合理施用[10]。長期過度的化學投入,導致土地持續(xù)生產(chǎn)力下降,為保證產(chǎn)量,農業(yè)發(fā)展進入“化學陷阱”?;拭嬖次廴局卫硇枰喾矫嫦到y(tǒng)控制,措施之一節(jié)水灌溉在減少農田氮磷流失方面具有一定優(yōu)勢,與傳統(tǒng)灌溉相比較,節(jié)水灌溉能夠節(jié)水25%以上,滲漏量(水稻)減少31.7%,氮肥利用率增加3%~5%[11]。水資源有效利用能夠緩解農業(yè)面源污染,但目前我國水資源利用粗放、浪費嚴重,在農業(yè)灌溉過程中水的有效利用率僅為30%~40%。除了節(jié)水灌溉,文獻整理得到影響化肥面源污染的因素很多,從農戶視角可分為內部因素和外部環(huán)境因素。內部因素包括土地經(jīng)營規(guī)模、家庭勞動力、家庭收入、農戶受教育程度、環(huán)保意識、政治參與程度等[12-15];外部環(huán)境因素包括地區(qū)經(jīng)濟發(fā)展水平、農業(yè)政策、種植結構、城鄉(xiāng)二元經(jīng)濟社會結構等[16-18]。目前存在的治理障礙主要是由于農業(yè)生產(chǎn)者、環(huán)境規(guī)制部門及其他利益相關者,均存在污染治理的政治、預算和技術上的障礙[19]。本文通過測度長江上游4省(市)797份樣本戶水稻灌溉用水效率和化肥面源污染量,基于分位數(shù)回歸方法,深入探討灌溉用水效率、生產(chǎn)方式和農戶要素稟賦對化肥面源污染的影響,以期對提高農業(yè)灌溉用水效率,緩解化肥面源污染提供依據(jù)。

    1 材料與方法

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

    本研究所用數(shù)據(jù)來自于2014年對長江上游地區(qū)四川、重慶、云南、貴州四省(市)水稻種植戶的調查,共整理得到有效問卷797份,覆蓋19個縣(區(qū))。各地區(qū)調研分布情況見表1。鑒于數(shù)據(jù)的可得性和調研方便性,問卷調研主要集中在四川,其中四川、貴州、云南和重慶問卷比例分別為76.66%、9.16%、8.06%、6.02%,四川省問卷主要分布在成都平原,由于成都平原灌溉水資源豐富充沛,對于灌溉設施使用率較高。為反映調查結果的客觀真實性,作了以下驗證,調查結果統(tǒng)計戶均水稻種植面積為1813 m2,與農業(yè)部農村經(jīng)濟研究中心“我國糧食安全發(fā)展戰(zhàn)略研究”課題組2011年在農業(yè)部種植業(yè)管理司支持下開展的西南地區(qū)505個種糧調查結果中戶均1867 m2基本一致。

    1.2 研究方法

    農業(yè)灌溉用水作為農業(yè)生產(chǎn)的基本要素投入,必須和其他生產(chǎn)要素科學配合,從而帶來農業(yè)經(jīng)濟效益產(chǎn)出。農業(yè)灌溉用水效率是指,每消耗單位水資源所產(chǎn)生的糧食產(chǎn)品。具體到水資源作為投入要素的角度看,農業(yè)灌溉用水效率即為實現(xiàn)最優(yōu)產(chǎn)出下的最少農業(yè)灌溉用水投入量與實際用水量的比值[20]。測度水資源使用效率常用的方法有以計量經(jīng)濟學為基礎的參數(shù)投入法和以線性規(guī)劃為基礎的非參數(shù)方法。相比參數(shù)投入法,非參數(shù)投入方法考察被評估區(qū)域與生產(chǎn)前沿面的相對差距,測度結果為某要素的相對效率,該方法不必人為設定指標的權重。本文采用非參數(shù)估計方法的數(shù)據(jù)包絡分析方法(DEA)。

    表1 調研樣本分布情況Table 1 Distribution of research sample

    1.2.1 DEA模型

    目前學術界效率測度方法主要包括以數(shù)據(jù)包絡分析法(DEA)為主的非參數(shù)方法和以隨機前沿分析法(SFA)為主的參數(shù)方法。與SFA等參數(shù)法相比, DEA方法具有自如處理多投入多產(chǎn)出指標的復雜問題,無需對生產(chǎn)函數(shù)的形式作出假設,從而避免了由于函數(shù)設定誤差所帶來的問題。因水稻灌溉用水效率受到氣候條件、水利設施、個體行為等多方面因素影響,生產(chǎn)函數(shù)無法確定,故選擇DEA分析方法。DEA方法可分為基于投入或產(chǎn)出兩種不同方法,基于投入的DEA方法目的是為了測算生產(chǎn)單元相對給定產(chǎn)出水平下最小可能投入的效率,而基于產(chǎn)出的DEA方法則是為了度量實際產(chǎn)出與給定投入水平的最大可能產(chǎn)出差距。只有在規(guī)模收益不變的情況下,兩種方法的效率測算結果才會相等。本文側重考察作為農業(yè)投入要素的灌溉用水效率,故采用規(guī)模報酬可變假設下基于投入導向的DEA模型,具體形式如下.

    式中:n為決策單元個數(shù)(n=297);m為輸入變量個數(shù)(m=4),s為輸出變量;x為投入要素,包括土地、勞動力、技術和資本;y為產(chǎn)出收益;δ為決策單元DMU0的有效值。若δ=1或s+=s-=0,則決策單元DEA有效;若δ=1且s+≠0或s-≠0,則決策單元為弱DEA有效;0≤δ<1,則決策單元非DEA有效。即若DEA模型測度結果為1,則說明灌溉效率為DEA有效,反之則非有效。

    1.2.2 計量回歸模型

    本文研究水稻灌溉用水效率、農戶生產(chǎn)方式和要素稟賦對化肥面源污染的影響,考慮其他主要控制變量,基本回歸模型設定如下.

    式中:yi分別為化肥面源污染總氮排放量(TN)和總磷排放量(TP),xi為水稻灌溉用水效率,通過DEA方法測度而得,∑controli代表所有的控制變量之和。具體的估計策略上,本文首先采用OLS回歸方法對(2)式模型進行估計。然而,OLS回歸方法只能得到水稻灌溉效率對化肥面源污染的期望值的影響,無法分析各個因素對面源污染的分布規(guī)律的影響。Koenker和Bassett于1978年提出的分位數(shù)回歸方法(Quantile Regression,QR)可解決這個問題。從理論上講,OLS回歸是擬合因變量的條件均值與自變量之間的線性關系,而分位數(shù)回歸是通過估計因變量在0~1之間的不同分位數(shù)值,對特定分布的數(shù)據(jù)進行估計。該方法能精確解釋對于被解釋變量的變化范圍以及條件分布形狀的影響,能全面描述被解釋變量條件分布的所有情形,還可以分析各分位數(shù)條件下解釋變量對被解釋變量的作用機制。為考察不同分位數(shù)上水稻灌溉用水效率對化肥面源污染的影響因素,本文采用的分位數(shù)回歸模型形式如下.

    式中:βq為q分位數(shù)回歸系數(shù),其估計量βq可以由以下最小化問題來定義.

    1.2.3 化肥面源污染測度方法

    經(jīng)過文獻查閱和借鑒,化肥面源污染量用總氮(TN)、總磷(TP)排放量來表示,其測度依據(jù)賴斯蕓、杜鵬飛、陳吉寧[21]研究方法,每個省市的氮、磷利用系數(shù)均借鑒其測度系數(shù),化肥單元產(chǎn)污強度計算公式為.

    其中四川、重慶氮流失系數(shù)為0.10,云南、貴州氮流失系數(shù)為0.20,4省(市)磷流失系數(shù)均為0.04,復合肥氮含量65%,磷含量26%。產(chǎn)污總量即為產(chǎn)污強度與化肥使用量的乘積,表示在降水或灌溉過程中,未被利用并通過地表徑流和農田滲漏形成地表和地下水環(huán)境污染的面源污染量。

    2 結果與分析

    2.1 水稻灌溉用水效率測度與分析

    水稻灌溉用水效率受到氣候、地形、水利條件及灌溉行為等因素影響,其中土地是農業(yè)生產(chǎn)的基礎和載體,水、勞動力、資本等投入要素需要依附在土地上才能實現(xiàn)農業(yè)生產(chǎn)。土地可被看作是農業(yè)生產(chǎn)投入的固定投入要素,水等被農作物直接消耗的經(jīng)濟資源則是農業(yè)生產(chǎn)的可變投入要素,兩者作用密不可分。基于數(shù)據(jù)可獲得性和相關文獻,選擇產(chǎn)出和投入變量。參考許朗等對農業(yè)用水灌溉效率的研究[22],為研究水稻生長季的灌溉用水效率,故選擇產(chǎn)出變量為水稻產(chǎn)量。投入指標包括.①家庭水稻耕種面積。水稻種植面積為樣本戶家庭經(jīng)營水稻耕地面積,單位為m2,反映了潛在產(chǎn)出量和需水量。②水稻種植總投工。水稻種植總投工包括水稻生長過程中整地、育秧、移栽、施肥等環(huán)節(jié)總投工時,單位為日,大部分調研點水稻生產(chǎn)仍屬于勞動密集型作業(yè),單位投工會影響產(chǎn)出。③水稻種植化肥施用量?;适┯昧繛樗旧L環(huán)節(jié)所施用的化學肥料總用量,單位為kg,水稻高產(chǎn)需借助化肥等生物技術,某些生物技術甚至可減弱水稻需水程度。④灌溉用水費用。灌溉用水費用為水稻生長環(huán)節(jié)灌溉用水總花費,補貼部分不計入,單位為元,農戶會參考灌溉成本而調整灌溉行為和投入。

    運用DEAP2.0軟件得到水稻種植戶灌溉總效率、純技術效率和規(guī)模效率,效率值分布情況見表2。種植戶水稻灌溉總效率平均值為0.14,純技術效率平均值為0.27,規(guī)模效率平均值為0.55??傂仕街饕性?.2以下,說明整體上水稻灌溉用水效率較低。僅有13戶(總效率值為1)水稻種植處于農業(yè)生產(chǎn)可能集的前沿包絡面上,說明在水稻灌溉方面實現(xiàn)了水資源的最有效利用,相對有效生產(chǎn)點,其他98.37%的農戶均處于生產(chǎn)的相對無效狀態(tài),水資源利用改進空間潛力很大。以總效率平均值0.137 8為例,意味著農戶在資源稟賦不變情況下,若生產(chǎn)技術和經(jīng)營水平充分發(fā)揮平均667 m2可節(jié)省86.22%的灌溉費用。而技術效率略高于總效率,規(guī)模效率明顯高于前兩者,是總效率的4倍左右。水稻灌溉用水純技術效率中有28戶處于有效生產(chǎn)點上,無效點主要集中在0.4以下,說明不同地區(qū)、農戶擁有的水資源、土地、勞動力等資源稟賦存在差異,但灌溉技術普遍較落后、粗放,水資源浪費嚴重。以平均值0.27為例,表示農戶水稻灌溉若技術完全運用,將不會造成用水浪費72.55%,最大值與最小值技術效率差值接近于1,存在較大的差異化和兩級化,灌溉設施和技術發(fā)揮效應弱化。規(guī)模效率相對較高,效率值主要集中在0.50上下,由于平原地區(qū)因土地流轉政策推動形成了規(guī)模種植,特別是平原地區(qū)和水資源豐富地區(qū),這些地區(qū)灌溉設施和機械化水平更高,進一步提升了農業(yè)生產(chǎn)的規(guī)模效應,這一良性循環(huán)積極促使農戶為提高收益而改善農田灌溉條件,進一步提升水稻灌溉用水規(guī)模效率。

    基于不同灌溉方式的灌溉用水效率分布結果見表3。灌溉方式基于灌溉工程進行分類,分為引水灌溉、提水灌溉、井水灌溉和蓄水灌溉。可知,從灌溉方式選擇來看,目前水稻水旱輪作種植戶主要采用引水灌溉(占比64.62%),其次是提水灌溉(占比32.25%),而井水灌溉和蓄水灌溉用戶最少。從灌溉效率平均值來看,四種灌溉方式效率值均較低,且差異不大。其中蓄水灌溉效率最高,平均效率值0.195 7,其后依次是引水灌溉(0.194 3)、提水灌溉(0.193 2)、井水灌溉(0.1808)。以引水灌溉為例,其平均效率值為0.194 3,則當農戶在資源稟賦不變的情況下,若生產(chǎn)技術和經(jīng)營水平充分發(fā)揮則平均每667 m2可節(jié)約80.57%的灌溉費用,同樣提水灌溉、井水灌溉和蓄水灌溉在相同條件下平均每667 m2用水費用可分別減少80.68%、81.92%和80.43%。從效率值分布結構看,四種方式各自效率值分布主要集中在0.2以下。引水、提水、井水和蓄水灌溉效率值位于0.2以下的比例分別為71.07%、69.26%、80.00%和60.00%,可見效率值普遍偏低,存在較大的改進空間。從水資源最有效利用情況看,只有引水灌溉和提水灌溉存在效率值在農業(yè)生產(chǎn)可能集的前沿包絡面上,且占比分別僅為0.78%、2.33%。

    表2 樣本農尸灌溉效率值分布Table 2 Sample farmers′irrigation efficiency value

    表3 基于灌溉方式的效率值分布Table 3 Irrigation efficiency value based on irrigation methods

    2.2 影響因素的指標選擇與描述性統(tǒng)計

    根據(jù)公式(5)測度化肥面源污染量,結果見表4所示?;拭嬖次廴究偟?TN)平均為7.90 kg,總磷量(TP)平均為0.52 kg,水稻種植施肥主要以氮肥和復合肥為主,樣本結果顯示兩者占比93.22%,由化肥產(chǎn)污強度系數(shù)可知,主要以氮排放為主,因篇幅限制,故本文化肥面源污染以總氮量為主,即為被解釋變量。

    本文主要解釋變量為水稻灌溉用水效率、生產(chǎn)方式和農戶要素稟賦。生產(chǎn)方式具體指標包括灌溉方式、水旱輪作模式、耕作方式和秸稈還田方式。農戶要素稟賦包括土地稟賦、勞動力稟賦、資本稟賦、社會資本稟賦和技術稟賦。其中土地稟賦包括水田面積、水田塊數(shù)、地形特征,勞動力稟賦用家庭勞動力數(shù)量表示,資本稟賦具體指標包括農業(yè)收入占比、農業(yè)補貼額、商品化率,社會資本稟賦具體指標包括戶主年齡、受教育程度、是否是黨員,技術稟賦用“是否接受過農業(yè)技術培訓”表示。各指標描述性統(tǒng)計見表4。灌溉方式和水旱輪作模式序號根據(jù)求得的平均效率值進行排列,兩者均值分別為2.67、2.61。不同灌溉技術措施反映出田間不同的水分管理,尤其是在水稻分蘗后期至成熟期,對田間土壤水分控制上的差異,條件和控制包括土壤中水、肥、氣、熱狀況的改變,從而影響到需水量大小。依據(jù)工程對灌溉方式劃分,不同灌溉技術條件下的變化趨勢.深水灌溉大于淺水灌溉,淺水灌溉大于濕潤或淺濕灌溉,控制灌溉的需水量最小。不同灌溉方式在儲水、輸水過程中存在水流失情況,澆灌過程會發(fā)生田間滲漏,進而造成水資源灌溉效率存在差異。水旱輪作模式反映了水稻與其他作物品種輪作選擇,而不同輪作作物對土壤的改善效果、輪作作物對水分需求量、種植經(jīng)驗、市場等因素共同決定了輪作模式的選擇,對于提升灌溉效率、改善灌溉條件、促進農業(yè)生產(chǎn)結構轉型具有重要參考價值。

    2.3 實證結果分析

    運用Stata14.0軟件對797個樣本戶使用自助法(自助法重復400次)進行分位數(shù)回歸(見表5)。受篇幅所限,表5只列出了化肥面源污染的第25、50、75、90個分位點回歸結果。同時為進一步解釋自變量對化肥面源污染影響的完整情況,圖1列出了化肥面源污染分位數(shù)回歸的系數(shù)變化情況。從回歸結果可以看出,灌溉用水效率、農戶生產(chǎn)方式和要素稟賦在不同分位數(shù)的影響系數(shù)有明顯變化,具體體現(xiàn)在以下幾個方面.

    (1)灌溉用水效率與化肥施用量存在替代效應,兩者呈振幅增強的波動式發(fā)展。灌溉用水總效率對化肥面源污染在不同分位點上呈10%顯著性負相關。因TN是化肥施用量的正向線性函數(shù),故灌溉用水效率與化肥施用量存在替代效應,或互補效應,即隨著灌溉用水效率的提升,化肥施用量減少。也驗證了余金鳳的研究,指出節(jié)水灌溉可以提高氮肥利用率3%~5%[11]。由于淹水灌溉和干濕交替灌溉均較旱作有效地協(xié)調各時期水稻地上部、地下部生長,促進各時期氮素吸收利用,提高稻谷產(chǎn)量[23]。節(jié)水灌溉具有保肥作用,因土壤水分減少,抑制了根系水分吸收,降低騰發(fā)量,減少無效葉面騰發(fā);節(jié)水灌溉因大部分時間田面無水層,蓄雨能力較強,肥力流失則會減少。總效率對灌溉用水效率的影響系數(shù)即為彈性系數(shù),彈性系數(shù)呈增大-減小-增大趨勢,且極大值點逐步擴大,形成替代效應正累積。

    (2)節(jié)水灌溉方式能夠減少化肥面源污染。灌溉方式在0.25和0.50分位點呈顯著負相關,在0.75和0.90分位點上未通過顯著性檢驗,相關關系呈“倒U型”。由井水灌溉轉為蓄水灌溉,灌溉效率增加了3.10%,實現(xiàn)了一定程度的節(jié)水效果。由于蓄水灌溉減少儲水和輸水過程的流失,尤其是在水資源短缺地區(qū),采取井水灌溉和蓄水灌溉農戶可根據(jù)作物生長情況及時安排灌溉,灌溉管理實行非充分灌溉,有利于養(yǎng)分吸收和作物生長。由于水分狀況的改變和土壤通氣條件的改善,促進了還原物質的氧化,加速了有機質的分解和遲效養(yǎng)分的活化,使非充分灌溉稻田土壤的肥力得到充分發(fā)揮,加上滲漏大幅度減少,又相應地減少了養(yǎng)分流失。實驗發(fā)現(xiàn)與淹水灌溉相比,控制灌溉水稻節(jié)水45.9%,水分利用效率提高6.3%~79.8%,氮肥利用率增加5.2%~38.4%[24],因節(jié)水灌溉能夠明顯抑制農田排水徑流量和滲漏量[25]。節(jié)水灌溉能顯著降低稻田滲漏量,同時改善土壤的通氣狀況,水稻在一定時期內處于水分脅迫狀態(tài),限制了土壤中有害物質的產(chǎn)生,高效控制了氮、磷和COD流失關鍵時期的排水,減輕了農業(yè)面源污染程度。

    表4 變量指標的統(tǒng)計描述Table 4 Variable indicators statistical description

    (3)水旱輪作模式、耕作方式、秸稈還田方式對化肥面源污染存在影響,但未全部通過顯著性檢驗。水旱輪作模式對化肥面源污染的影響系數(shù)為負,說明輪作結構不同導致面源污染效應差異。薛利紅等基于對太湖流域連續(xù)3年田間數(shù)據(jù)比較,發(fā)現(xiàn)稻-紫云英能減少徑流總氮損失35%~40%,稻-蠶豆輪作能減少徑流氮排放25%~30%[26];洱海流域田間試驗發(fā)現(xiàn)水稻-蠶豆比水稻大蒜輪作模式減少氮素流失風險38%[27]。秸稈還田形成綠肥,可提高土壤肥力,降低稻季施肥量,減少稻季氮肥流失引起的環(huán)境風險[28]。目前堆漚還田對減少化肥面源污染效果最大,其次是直接還田,最小是焚燒還田。肖新成等[14]和吳義根等[29]也指出種植結構直接影響農業(yè)面源污染排放效率和總量。耕作方式對化肥面源污染的彈性系數(shù)為負,秸稈還田方式在0.75分位點以下為正,且在0.25分位點上5%顯著,0.90分位點上為負。實際實驗中發(fā)現(xiàn)在翻耕20 cm的基礎上進行秸稈還田或增施綠肥、畜禽糞肥則可以大幅提高水稻產(chǎn)量和氮、磷、鉀養(yǎng)分吸收量[30]。

    (4)耕地細碎化會增加化肥面源污染,一定規(guī)模下隨面積增加化肥面源污染相應增加。水田面積在各分位點上與化肥面源污染呈1%顯著性正相關,水田塊數(shù)在0.25分位點上通過顯著性檢驗。在經(jīng)營面積較大時,農戶會因勞動力不足而選擇通過增施化肥、農藥等物質生產(chǎn)資料來提升農業(yè)生產(chǎn)效益。細碎化程度越高,單位化肥投入強度越大。有研究指出適度規(guī)模能夠降低化肥投入強度,進而控制因化肥過量投入所帶來的面源污染壓力[31]。還應充分考慮地權穩(wěn)定性,若通過土地流轉擴大規(guī)模,則不一定能夠有效減少面源污染問題。龍云指出農地流轉的地權穩(wěn)定性效應、地權流動性效應和農地利用政策效應會導致轉入地農戶的農地利用方式和短期投資發(fā)生改變,會使轉入地總體面源污染水平更高[32]。結果還發(fā)現(xiàn)丘陵山區(qū)相對平原地區(qū)化肥面源污染更大,平原地區(qū)地下水充沛,田間滲漏相對更少,而丘陵山區(qū)因坡度造成水分滲漏嚴重,保肥能力較弱,農戶選擇通過增施化肥來保證產(chǎn)量。

    表5 實證回歸結果Table 5 The empirical regression results

    (5)農業(yè)補貼和商品化率會增加化肥面源污染。農業(yè)補貼在0.50分位點以下與面源污染呈正相關, 0.75以上呈負相關,說明農業(yè)補貼與化肥面源污染存在“倒U型”關系。已有研究表明農資補貼對農業(yè)化肥面源污染存在正向激發(fā)效應[33],葛繼紅的研究也證實了國家財政支農政策導致了化肥要素市場扭曲的存在[34]?;室厥袌雠で鷮е禄蔬呺H產(chǎn)出與實際價格的偏離,同時相對低廉的化肥要素價格形成了對勞動力的替代,便助漲了農戶對化肥的過量施用。商品化率越高,農戶為提高單位產(chǎn)出和收益,增加化肥投入強度,進而導致化肥面源污染增加。

    (6)年齡和受教育程度減弱化肥面源污染的臨界點較高,政治參與程度對化肥面源污染影響不顯著。年齡和受教育程度都在0.25和0.90分位點上與化肥面源污染負相關,在0.50和0.75分位點上正相關,但都未通過一致性檢驗。圖1可知,化肥面源污染隨年齡和受教育程度增加初期長時間呈緩慢提高趨勢,當達到一定程度后,出現(xiàn)明顯下降,且臨界點位置較高。政治參與程度對化肥面源污染影響系數(shù)為負,在各分位點上都不顯著,說明政治參與程度越高,獲得生產(chǎn)技術和產(chǎn)品市場信息越多,科學種植理念和環(huán)保意識越高,偏向于減少化肥投入。這是由于農業(yè)面源污染具有負外部性,大多數(shù)農戶不會關心生產(chǎn)中的環(huán)境污染問題[35],反而還會通過增施化肥來提高產(chǎn)量和節(jié)約勞動力。該結論也得到了付靜塵等的驗證[36]。

    (7)技術培訓能夠有效降低化肥面源污染。技術培訓在0.50分位點以下通過顯著性檢驗,與面源污染呈正相關關系,說明技術培訓可有效降低農戶化肥施用量。華春林等也驗證了農業(yè)教育培訓可以減少農業(yè)面源污染[37]。但應瑞瑤等指出在考慮了樣本的選擇性偏誤問題之后,農業(yè)技術培訓對農業(yè)生產(chǎn)中面源污染的減少效果并不明顯,這與中國農業(yè)技術培訓體系注重產(chǎn)品數(shù)量安全而忽視農產(chǎn)品生態(tài)環(huán)境安全以及“從上到下”的行政命令式培訓方式有關[38]。

    圖1 化肥面源污染分位數(shù)回歸的系數(shù)變化情況Figure 1 The changes of fertilizer non-point source pollution of quantile regression coefficient

    3 結論

    續(xù)圖1化肥面源污染分位數(shù)回歸的系數(shù)變化情況Continuous figure 1 The changes of fertilizer non-point source pollution of quantile regression coefficient

    (1)總體上灌溉總效率均值為0.14,純技術效率均值為0.27,規(guī)模效率均值為0.55,水資源利用改進空間較大,規(guī)模化程度高的平原地區(qū)用水效率更高。

    (2)灌溉用水效率與化肥施用量間存在顯著替代效應,兩者呈振幅增強的波動式發(fā)展。灌溉用水效率對化肥面源污染的彈性系數(shù)在不同分位點上依次為-2.38、-3.63、-2.23、-5.87,呈增大-減小-增大趨勢,且極大值逐步擴大,形成替代效應正累積。

    (3)節(jié)水灌溉能夠降低化肥面源污染,同時具有保肥效果。灌溉方式與化肥面源污染呈負相關,且在0.25和0.50分位點上通過顯著性檢驗。節(jié)水灌溉既可以提高水資源利用效率,降低稻田滲漏量,還具有控制氮、磷和COD流失,減輕面源污染的效果。

    (4)水旱輪作模式、耕作方式、秸稈還田方式對化肥面源污染存在影響,但未全部通過顯著性檢驗。水旱輪作對化肥面源污染的影響系數(shù)為負,輪作作物根系固氮能力影響了化肥氮素利用率;耕作方式的影響系數(shù)為負,但不顯著,深耕能夠提高土壤肥力、減少化肥施用;秸稈還田方式在0.25分位點上顯著正相關,對減輕化肥面源污染效果大小依次是堆漚還田、直接還田、焚燒還田。

    (5)耕地細碎化、農業(yè)補貼和商品化率會增加化肥面源污染。耕地面積在各分位點上都與化肥面源污染呈1%顯著性正相關,耕地細碎化的系數(shù)為正,在0.25分位點上1%顯著。一定規(guī)模下隨面積增加化肥面源污染會相應增加,而細碎化程度越高,單位化肥投入強度越大。

    [1]Baumg?rtner S,Dyckhoff H,Faber M.The concept of joint production and ecological economics[J].Ecological Economics,2001,36(3).365-372.

    [2]李秀芬,朱金兆,顧曉君,等.農業(yè)面源污染現(xiàn)狀與防治進展[J].中國人口·資源與環(huán)境,2010,20(4).81-84.

    LI Xiu-fen,ZHU Jin-zhao,GU Xiao-jun,et al.Current situation and control of agricultural non-point source pollution[J].China Population,Resources and Environment,2010,20(4).81-84.

    [3]張福鎖,王激清,張衛(wèi)峰,等.中國主要糧食作物肥料利用率現(xiàn)狀與提高途徑[J].土壤學報,2008,45(5).915-924.

    ZHANG Fu-suo,WANG Ji-qing,ZHANG Wei-feng,et al.Nutrient use efficiencies of major cereal crops in China and measures for improvement[J].ActaPedologicaSinica,2008,45(5).915-924.

    [4]于海蛟,趙樹君,崔遠來,等.南方典型水稻灌區(qū)肥藥及廢棄物空間變異規(guī)律[J].中國農村水利水電,2016(8).190-195.

    YU Hai-jiao,ZHAO Shu-jun,CUI Yuan-lai,et al.Spatial distribution of rice fertilizer,pesticide and waste disposal in south China’s irrigation district[J].ChinaRural Water and Hydropower,2016(8).190-195.

    [5]張郁,張崢,蘇明濤.基于化肥污染的黑龍江墾區(qū)糧食生產(chǎn)灰水足跡研究[J].干旱區(qū)資源與環(huán)境,2013,27(7).28-32.

    ZHANG Yu,ZHANG Zheng,SU Ming-tao.Research on grey water footprint based on chemical fertilizer use in the grain production in Heilongjiang reclamation area[J].Journal of Arid Land Resources and Environment,2013,27(7).28-32.

    [6]張鋒.中國化肥投入的面源污染問題研究[D].南京.南京農業(yè)大學,2011.47-49.

    ZHANG Feng.Study on the non-point pollution of the agricultural fertilizer input in China-based on the aspect of farmers’fertilizer input[D]. Nanjing.Nanjing Agricultural University,2011.47-49.

    [7]唐學玉,張海鵬,李世平.農業(yè)面源污染防控的經(jīng)濟價值.基于安全農產(chǎn)品生產(chǎn)戶視角的支付意愿分析[J].中國農村經(jīng)濟,2012(3).53-67.

    TANG Xue-yu,ZHANG Hai-peng,LI Shi-ping.The economic value of the prevention and control of agricultural surface pollution.Based on the willingness to pay for the production of safety agricultural products[J].Chinese Rural Economy,2012(3).53-67.

    [8]馬云澤.當前中國農村環(huán)境污染問題的根源及對策.基于規(guī)制經(jīng)濟學的研究視角[J].廣西民族大學學報(哲學社會科學版),2010,32 (1).18-21.

    MA Yun-ze.On the causes of and countermeasures against environmental pollution in China′s rural areas.From the research perspective of regulation economics[J].Journal of Guangxi University for Nationalities(Philosophy and Social Sciences Edition),2010,32(1).18-21.

    [9]金書秦,沈貴銀,魏珣,等.論農業(yè)面源污染的產(chǎn)生和應對[J].農業(yè)經(jīng)濟問題,2013(11).97-102.

    JIN Shu-qin,SHEN Gui-yin,WEI Xun,et al.The production and response of the contamination of agricultural surface sources[J].Issues in Agricultural Economy,2013(11).97-102.

    [10]饒靜,許翔宇,紀曉婷.我國農業(yè)面源污染現(xiàn)狀、發(fā)生機制和對策研究[J].農業(yè)經(jīng)濟問題,2011(8).81-87.

    RAO Jing,XU Xiang-yu,JI Xiao-ting.Study on the present situation, mechanism and countermeasures of the contamination of agricultural surface source in China[J].Issues in Agricultural Economy,2011(8). 81-87.

    [11]余金鳳,洪林,江洪珊.南方典型灌區(qū)節(jié)水灌溉的減污效應[J].節(jié)水灌溉,2011(8).1-4.

    YU Jin-feng,HONG Lin,JIANG Hong-shan.Effects of water-saving irrigation on pollution control in typical irrigation district in southern China[J].Water SavingIrrigation,2011(8).1-4.

    [12]汪厚安,葉慧,王雅鵬.農業(yè)面源污染與農戶經(jīng)營行為研究.對湖北農戶的實證調查與分析[J].生態(tài)經(jīng)濟,2009(9).87-91.

    WANG Hou-an,YE Hui,WANG Ya-peng.Study on the agricultural non-point source pollution and farm household′s management behavior.An empirical study and analysis of farmer households in Hubei Province[J].Ecological Economy,2009(9).87-91.

    [13]李傳桐,張廣現(xiàn).農業(yè)面源污染背后的農戶行為.基于山東省昌樂縣調查數(shù)據(jù)的面板分析[J].地域研究與開發(fā),2013,32(1).143-146,164.

    LI Chuan-tong,ZHANG Guang-xian.The farm household’s behavior behind agricultural non-point source pollution.A panel analysis on change county in Shandong Province[J].Areal Research and Development,2013,32(1).143-146,164.

    [14]肖新成,何丙輝,倪九派,等.三峽生態(tài)屏障區(qū)農業(yè)面源污染的排放效率及其影響因素[J].中國人口·資源與環(huán)境,2014,24(11).60-68.

    XIAO Xin-cheng,HE Bing-hui,NI Jiu-pai,et al.Study on emission efficiency,shadow price and motivation factors of agricultural nonpoint source pollution in ecological barrier zone of Three Gorges Reservoir Area[J].China Population,Resources and Environment,2014,24 (11).60-68.

    [15]向濤,綦勇.糧食安全與農業(yè)面源污染.以農地稟賦對化肥投入強度的影響研究[J].財經(jīng)研究,2015,41(7).132-144.

    XIANG Tao,QI Yong.Food security and agricultural non-point source pollution.Taking the impact of agricultural land endowments on fertilizer use intensity as an example[J].Journal of Finance and Economics, 2015,41(7).132-144.

    [16]Repetto R.Economic incentives for sustainable production[J].The Annals of Regional Science,1987(21).44-59.

    [17]杜江,羅珺.我國農業(yè)面源污染的經(jīng)濟成因透析[J].中國農業(yè)資源與區(qū)劃,2013,34(4).22-27.

    DU Jiang,LUO Jun.Insight to the economic reasons of agricultural non-point source pollution in China[J].Chinese Journal of Agricultural Resources and Regional Planning,2013,34(4).22-27.

    [18]梁流濤,曲福田,馮淑怡.經(jīng)濟發(fā)展與農業(yè)面源污染.分解模型與實證研究[J].長江流域資源與環(huán)境,2013(10).1369-1374.

    LIANG Liu-tao,QU Fu-tian,FENG Shu-yi.Economic development and agricultural non-point source pollution.decomposition model and empirical analysis[J].Resources and Environment in the Yangtze Basin,2013(10).1369-1374.

    [19]Dowda B M,Press D,Huertos M L.Agricultural Nonpoint Source Water Pollution Policy.the Case of California′s Central Coast[J].Agriculture,Ecosystems and Environment,2008,128(11).151-161.

    [20]王學淵,趙連閣.中國農業(yè)用水效率及影響因素.基于1997—2006年省區(qū)面板數(shù)據(jù)的SFA分析[J].農業(yè)經(jīng)濟問題,2008(3).10-20.

    WANG Xue-yuan,ZHAO Lian-ge.Agricultural water efficiency and the causal factors.A stochastic frontier analysis is based on Chinese provincial panel data 1997—2006[J].Issues in Agricultural Economy, 2008(3).10-20.

    [21]賴斯蕓,杜鵬飛,陳吉寧.基于單元分析的非點源污染調查評估方法[J].清華大學學報(自然科學版),2004,44(9).1184-1187.

    LAI Si-yun,DU Peng-fei,CHEN Ji-ning.Evaluation of non-point source pollution based on unit analysis[J].Journal of Tsinghua University(Science and Technology),2004,44(9).1184-1187.

    [22]許朗,黃鶯.農業(yè)灌溉用水效率及其影響因素分析.基于安徽省蒙城縣的實地調查[J].資源科學,2012,34(1).105-113.

    XU Lang,HUANG Ying.Measurement of irrigation water efficiencyand analysis of influential factors.An empirical study of Mengcheng county in Anhui Province[J].Resources Science,2012,34(1).105-113.

    [23]嚴奉君,孫永健,馬均,等.灌溉方式與秸稈覆蓋優(yōu)化施氮模式對秸稈腐熟特征及水稻氮素利用的影響[J].中國生態(tài)農業(yè)學報, 2016,24(11).1435-1444.

    YAN Feng-jun,SUN Yong-jian,MA Jun,et al.Effects of irrigation method and straw mulch-nitrogen management pattern on straw decomposition characteristics and nitrogen utilization of hybrid rice[J].Chinese Journal of Eco-Agriculture,2016,24(11).1435-1444.

    [24]龐桂斌,楊士紅,徐俊增.節(jié)水灌溉稻田水肥調控技術試驗研究[J].節(jié)水灌溉,2015(9).44-47,51.

    PANG Gui-bin,YANG Shi-hong,XU Jun-zeng.Field experimental study on water and fertilizer managements of rice paddy under water saving irrigation[J].Water SavingIrrigation,2015(9).44-47,51.

    [25]洪林,袁自瑛,江海濤.漳河灌區(qū)農田排水氮素流失規(guī)律研究[J].武漢大學學報(工學版),2009,42(6).681-684,690.

    HONG Lin,YUAN Zi-ying,JIANG Hai-tao.Study of laws of nitrogen losses from farmland drainage in Zhanghe irrigation district[J].EngineeringJournal of Wuhan University,2009,42(6).681-684,690.

    [26]薛利紅,楊林章,施衛(wèi)明,等.農村面源污染治理的“4R”理論與工程實踐.源頭減量技術[J].農業(yè)環(huán)境科學學報,2013,32(5).881-888.

    XUE Li-hong,YANG Lin-zhang,SHI Wei-ming,et al.Reduce-Retain-Reuse-Restore technology for controlling the agricultural nonpoint pollution in countryside in China.Source reduction technology[J].Journal of Agro-Environment Science,2013,32(5).881-888.

    [27]湯秋香,任天志,雷寶坤,等.洱海北部地區(qū)不同輪作農田氮、磷流失特性研究[J].植物營養(yǎng)與肥料學報,2011,17(3).608-615.

    TANG Qiu-xiang,REN Tian-zhi,LEI Bao-kun,et al.Characteristics of nitrogen and phosphorus loss in various crop rotation systems in northern watershed of Erhai Lake[J].Plant Nutrition and Fertilizer Science,2011,17(3).608-615.

    [28]盧萍,單玉華,楊林章,等.綠肥輪作還田堆稻田土壤溶液氮素變化及水稻產(chǎn)量的影響[J].土壤,2006,38(3).270-275.

    LU Ping,SHAN Yu-hua,YANG Lin-zhang,et al.Influence of green manure crop on nitrogen concentration in soil solution of paddy field and rice yield[J].Soils,2006,38(3).270-275.

    [29]吳義根,馮開文,李谷成.人口增長、結構調整與農業(yè)面源污染.基于空間面板STIRPAT模型的實證研究[J].農業(yè)技術經(jīng)濟,2017(3). 75-87.

    WU Yi-gen,FENG Kai-wen,LI Gu-cheng.Population growth,structural adjustment and agricultural non-point source pollution.Based on thespacepanel STIRPATmodel[J].Journal of Agrotechnical Economics, 2017(3).75-87.

    [30]吳萍萍,王家嘉,李錄久.不同耕作與施肥方式下白土的水稻產(chǎn)量及養(yǎng)分吸收量[J].植物營養(yǎng)與肥料學報,2014,20(3).754-760.

    WU Ping-ping,WANG Jia-jia,LI Lu-jiu.Rice yields and nutrient uptake under different ploughing depths and fertilizations in white soil[J].Journal of Plant Nutrition and Fertilizer,2014,20(3).754-760.

    [31]張鋒,胡浩.農戶化肥投入行為與面源污染問題研究[J].江西農業(yè)學報,2012,24(1).183-186.

    ZHANG Feng,HU Hao.Study on farmer′s chemical fertilizer input behavior and non-point source pollution[J].Acta Agriculturae Jiangxi, 2012,24(1).183-186.

    [32]龍云,任力.農地流轉對農業(yè)面源污染的影響.基于農戶行為視角[J].經(jīng)濟學家,2016(8).81-87.

    LONG Yun,REN Li.The effect of agricultural land circulation on the agricultural non-point source pollution.Based on the view of farmer′s behavior[J].Economist,2016(8).81-87.

    [33]于偉詠,漆雁斌,余華.農資補貼對化肥面源污染效應的實證研究.基于省級面板數(shù)據(jù)[J].農村經(jīng)濟,2017(2).89-94.

    YU Wei-yong,QI Yan-bin,YU Hua.An empirical study of agricultural subsidies for fertilizer-point source pollution effects.Based on provincial panel data[J].Rural Economy,2017(2).89-94.

    [34]葛繼紅,周曙東.要素市場扭曲是否激發(fā)了農業(yè)面源污染.以化肥為例[J].農業(yè)經(jīng)濟問題,2012(3).92-98.

    GE Ji-hong,ZHOU Shu-dong.Does factor market distortions stimulate the agricultural non-point source pollution?Acase study of fertilizer[J].Issues in Agricultural Economy,2012(3).92-98.

    [35]Griffin R C,Bromley D W.Agricultural Runoff as a Nonpoint Externality.A Theoretical Development[J].American Journal of Agricultural Economics,1983(70).37-49.

    [36]付靜塵,韓烈保.丹江口庫區(qū)農戶對面源污染的認知度及生產(chǎn)行為分析[J].中國人口·資源與環(huán)境,2010,20(5).70-74.

    FU Jing-chen,HAN Lie-bao.Analysis on farmers cognition degree and production behavior with non-point pollution in Danjiangkou Reservoir area[J].China Population,Resources and Environment,2010,20(5). 70-74.

    [37]華春林,陸遷,姜雅莉,等.農業(yè)教育培訓項目對減少農業(yè)面源污染的影響效果研究[J].農業(yè)技術經(jīng)濟,2013(4).83-92.

    HUA Chun-lin,LU Qian,JIANG Ya-li,et al.Study on the agricultural education training program studies the effects of reducing agricultural contamination[J].Journal of Agrotechnical Economics,2013(4).83-92.

    [38]應瑞瑤,朱勇.農業(yè)技術培訓對農業(yè)面源污染的效果評估[J].統(tǒng)計與信息論壇,2016,31(1).100-105.

    YING Rui-yao,ZHU Yong.Impact of agricultural training on reducing agricultural non-point source pollution[J].Statistics&Information Forum,2016,31(1).100-105.

    The effect of rice irrigation efficiency and related factors on fertilizer non-point source pollution based on quantile regression

    YU Wei-yong1,QI Yan-bin1,HE Yue1,2,DENG Xin3
    (1.College of Economics,Sichuan Agricultural University,Chengdu 611130,China;2.Institute of Finance and Economics,Yangtze Normal U-niversity,Chongqing Fuling 408000,China;3.College of Management,Sichuan Agricultural University,Chengdu 611130,China)

    This article measures rice irrigation water-use efficiency and fertilizer non-point source pollution of 797 households across four provinces(municipalities)in the upper reaches of Yangtze river.It analyzed the effect of rice irrigation water-use efficiency,mode of production,and effects of other factors on fertilizer non-point source pollution based on quantile regression.The results show that the average irrigation efficiency of rice is 0.14,the average scale efficiency is 0.55,the average technical efficiency is 0.27,and the rice irrigation water efficiency improvement space is larger.A substitution effect between irrigation water-use efficiency and the amount of fertilizer applied exists,which enhances the amplitude of fluctuation.Water-saving irrigation methods can reduce fertilizer non-point source pollution and preserve fertilizers.Floods and droughts(FDD)rotation mode,farming methods,and straw counters-field ways have an effect on fertilizer non-point source pollution,although not all these effects were significant.Cultivated land and finely,agricultural subsidies and commercialization rate can increase non-point source fertilizer pollution.The critical point that age and education level reduced fertilizer non-pointsource pollution is higher,but political participation has no significant effect.

    fertilizer non-point source pollution;irrigation water efficiency;production mode;factor endowments;quantile regression

    X592

    A

    1672-2043(2017)07-1274-11

    10.11654/jaes.2017-0586

    于偉詠,漆雁斌,何悅,等.水稻灌溉用水效率和要素稟賦對化肥面源污染的影響——基于分位數(shù)回歸的分析[J].農業(yè)環(huán)境科學學報,2017,36 (7).1274-1287.

    YU Wei-yong,QI Yan-bin,HE Yue,et al.The effect of rice irrigation efficiency and related factors on fertilizer non-point source pollution based on quantile regression[J].Journal of Agro-Environment Science,2017,36(7).1274-1284.

    2017-04-20

    于偉詠(1988—),男,博士研究生,四川雅安人,主要研究方向為農業(yè)資源環(huán)境。E-mail:ywy212@163.com

    *通信作者:漆雁斌E-mail:qybin@sina.com

    國家社會科學基金項目(14XGL003);四川省農村發(fā)展研究中心青年項目(CR1624,CR1719)

    Project supported:The National Social Science Foundation of China(14XGL003);The Young Fund of the Sichuan Rural Development Research Center (CR1624,CR1719)

    猜你喜歡
    面源化肥用水
    7月我國化肥進口量增加
    哪些火災不能用水撲滅?
    機電安全(2022年1期)2022-08-27 02:14:50
    農業(yè)面源污染的危害與治理
    現(xiàn)代農業(yè)離不開化肥
    澄江市農業(yè)面源污染成因及對策
    節(jié)約洗碗用水
    化肥農藥減量,怎么減?
    化肥需求增長不均 有人歡喜有人憂
    基于SWAT模型的漳河流域面源污染模擬研究
    農業(yè)面源污染對水質的影響及防治對策
    久久久久久久亚洲中文字幕| 99久久人妻综合| 最后的刺客免费高清国语| 国产精品三级大全| 国产精品秋霞免费鲁丝片| 中国三级夫妇交换| 免费观看的影片在线观看| 狠狠精品人妻久久久久久综合| 另类亚洲欧美激情| 91久久精品电影网| 大码成人一级视频| 伦理电影免费视频| 欧美成人一区二区免费高清观看| 在线免费十八禁| 天堂俺去俺来也www色官网| 一级毛片aaaaaa免费看小| 你懂的网址亚洲精品在线观看| 最新中文字幕久久久久| 亚洲婷婷狠狠爱综合网| 亚洲国产精品国产精品| 成人亚洲精品一区在线观看 | 麻豆国产97在线/欧美| 亚洲欧美清纯卡通| 成年av动漫网址| 十八禁网站网址无遮挡 | 舔av片在线| 国产 精品1| 又黄又爽又刺激的免费视频.| 亚洲成色77777| 又大又黄又爽视频免费| 王馨瑶露胸无遮挡在线观看| 中文在线观看免费www的网站| 亚洲色图av天堂| 少妇人妻一区二区三区视频| 高清不卡的av网站| av线在线观看网站| 全区人妻精品视频| 啦啦啦在线观看免费高清www| 美女内射精品一级片tv| 在线观看人妻少妇| 国产精品一区二区性色av| 色吧在线观看| 老司机影院成人| 丝袜喷水一区| 夜夜看夜夜爽夜夜摸| 精品久久久久久电影网| 天天躁日日操中文字幕| 久久久久久久久久久丰满| 大香蕉久久网| 色婷婷av一区二区三区视频| 一二三四中文在线观看免费高清| av女优亚洲男人天堂| 男女边吃奶边做爰视频| 日韩欧美一区视频在线观看 | 天天躁日日操中文字幕| 欧美成人一区二区免费高清观看| 在线观看免费高清a一片| 如何舔出高潮| 夫妻性生交免费视频一级片| 在线观看三级黄色| 久久精品夜色国产| 色5月婷婷丁香| 久久99精品国语久久久| 亚洲国产精品一区三区| 国产v大片淫在线免费观看| 99re6热这里在线精品视频| 国产一区二区三区av在线| 美女xxoo啪啪120秒动态图| 一级黄片播放器| 男人添女人高潮全过程视频| 中文在线观看免费www的网站| 成人黄色视频免费在线看| 久久97久久精品| 搡老乐熟女国产| 91在线精品国自产拍蜜月| 欧美日韩国产mv在线观看视频 | 国产精品久久久久久久久免| 女人久久www免费人成看片| 3wmmmm亚洲av在线观看| 最后的刺客免费高清国语| 观看免费一级毛片| 日韩人妻高清精品专区| 青春草亚洲视频在线观看| 少妇熟女欧美另类| 看免费成人av毛片| 色吧在线观看| 韩国av在线不卡| 欧美老熟妇乱子伦牲交| 男女下面进入的视频免费午夜| 一二三四中文在线观看免费高清| 高清在线视频一区二区三区| 日韩欧美 国产精品| 少妇的逼好多水| 成年人午夜在线观看视频| 日韩电影二区| av播播在线观看一区| 综合色丁香网| 尤物成人国产欧美一区二区三区| 久久人妻熟女aⅴ| 久久久久视频综合| 少妇熟女欧美另类| 久久女婷五月综合色啪小说| 啦啦啦在线观看免费高清www| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 久久久午夜欧美精品| 亚洲一级一片aⅴ在线观看| 国产国拍精品亚洲av在线观看| 日韩av免费高清视频| 亚洲av福利一区| 精品久久国产蜜桃| 亚洲,一卡二卡三卡| 女人十人毛片免费观看3o分钟| 欧美一区二区亚洲| 成人二区视频| 久久人妻熟女aⅴ| 亚洲欧洲日产国产| 日韩一区二区三区影片| 亚洲av二区三区四区| 少妇人妻久久综合中文| 少妇人妻精品综合一区二区| 五月开心婷婷网| 国产成人91sexporn| 国产在视频线精品| 久久精品熟女亚洲av麻豆精品| freevideosex欧美| 卡戴珊不雅视频在线播放| 久久久a久久爽久久v久久| 国产精品一区二区性色av| 噜噜噜噜噜久久久久久91| 色视频www国产| 成人高潮视频无遮挡免费网站| 日韩精品有码人妻一区| 我要看黄色一级片免费的| 成人一区二区视频在线观看| 一区二区三区精品91| 91在线精品国自产拍蜜月| 日日啪夜夜爽| 亚洲精品自拍成人| 欧美xxxx黑人xx丫x性爽| 亚洲国产毛片av蜜桃av| 国产真实伦视频高清在线观看| 啦啦啦视频在线资源免费观看| 高清在线视频一区二区三区| 只有这里有精品99| 午夜福利影视在线免费观看| 亚洲欧美成人精品一区二区| 五月伊人婷婷丁香| 下体分泌物呈黄色| 欧美激情极品国产一区二区三区 | 狂野欧美激情性bbbbbb| 亚洲精品成人av观看孕妇| 成年av动漫网址| 国产精品久久久久久久久免| 另类亚洲欧美激情| 国产一区二区三区综合在线观看 | 久久av网站| 成人免费观看视频高清| 精品国产露脸久久av麻豆| 国产黄频视频在线观看| 亚洲va在线va天堂va国产| 国产色婷婷99| 日本与韩国留学比较| 91精品国产国语对白视频| 在线观看一区二区三区激情| 欧美日韩在线观看h| 91精品伊人久久大香线蕉| tube8黄色片| 久久久久久九九精品二区国产| 午夜免费鲁丝| 国产免费又黄又爽又色| 男女下面进入的视频免费午夜| 毛片女人毛片| 免费看光身美女| 亚洲av.av天堂| 丝袜脚勾引网站| 91精品一卡2卡3卡4卡| 在线观看一区二区三区| 国产精品熟女久久久久浪| 欧美日本视频| 一本—道久久a久久精品蜜桃钙片| 水蜜桃什么品种好| 99精国产麻豆久久婷婷| 舔av片在线| 亚洲精品一区蜜桃| 一区二区三区精品91| av国产久精品久网站免费入址| 青春草视频在线免费观看| 黄片无遮挡物在线观看| 久久久色成人| 欧美激情极品国产一区二区三区 | 亚洲久久久国产精品| 久久久成人免费电影| av女优亚洲男人天堂| 免费看光身美女| 亚洲自偷自拍三级| 久久久久久人妻| 男女免费视频国产| 久久 成人 亚洲| 精品久久久噜噜| 国产精品福利在线免费观看| 下体分泌物呈黄色| 国产av码专区亚洲av| 在线免费十八禁| 免费大片黄手机在线观看| 成人国产av品久久久| 久久99热这里只频精品6学生| av在线老鸭窝| 国内揄拍国产精品人妻在线| 成年免费大片在线观看| 亚洲av日韩在线播放| 日本爱情动作片www.在线观看| 观看免费一级毛片| 久久久久精品久久久久真实原创| 国产视频内射| 人体艺术视频欧美日本| 最近中文字幕高清免费大全6| 久久国产精品男人的天堂亚洲 | tube8黄色片| 国产午夜精品久久久久久一区二区三区| 欧美xxxx黑人xx丫x性爽| 97超碰精品成人国产| 久久久成人免费电影| 久久精品国产亚洲网站| 国产免费福利视频在线观看| 毛片一级片免费看久久久久| 人人妻人人爽人人添夜夜欢视频 | 免费大片18禁| 在线观看免费视频网站a站| 亚洲欧美成人综合另类久久久| 亚洲欧美一区二区三区国产| 午夜福利在线观看免费完整高清在| 99九九线精品视频在线观看视频| 中文字幕久久专区| www.色视频.com| 国产中年淑女户外野战色| 国产黄片美女视频| 国产黄片视频在线免费观看| 国产精品国产av在线观看| 搡女人真爽免费视频火全软件| 国产中年淑女户外野战色| 高清午夜精品一区二区三区| 久久 成人 亚洲| 小蜜桃在线观看免费完整版高清| 日韩三级伦理在线观看| 特大巨黑吊av在线直播| 自拍偷自拍亚洲精品老妇| 国产在视频线精品| 三级经典国产精品| 人人妻人人添人人爽欧美一区卜 | 日本wwww免费看| 少妇 在线观看| 全区人妻精品视频| 美女福利国产在线 | 777米奇影视久久| 美女视频免费永久观看网站| 亚洲精品日韩在线中文字幕| 国产精品嫩草影院av在线观看| 国产成人freesex在线| 啦啦啦在线观看免费高清www| 免费黄网站久久成人精品| 亚洲精品乱码久久久久久按摩| 国产免费一区二区三区四区乱码| 色综合色国产| 深爱激情五月婷婷| 免费人成在线观看视频色| 久久久久精品久久久久真实原创| 嘟嘟电影网在线观看| 高清黄色对白视频在线免费看 | av黄色大香蕉| 国产精品无大码| 色婷婷久久久亚洲欧美| 色5月婷婷丁香| 亚洲怡红院男人天堂| 国产精品一区www在线观看| 一级毛片aaaaaa免费看小| 精品久久久久久久久亚洲| 日韩 亚洲 欧美在线| 日本免费在线观看一区| 18禁在线无遮挡免费观看视频| 日韩在线高清观看一区二区三区| 国产精品免费大片| 99热这里只有是精品50| 国产精品国产三级国产av玫瑰| 美女视频免费永久观看网站| 91aial.com中文字幕在线观看| 在线 av 中文字幕| 亚洲国产高清在线一区二区三| 久久 成人 亚洲| 久久亚洲国产成人精品v| xxx大片免费视频| 直男gayav资源| 国产午夜精品久久久久久一区二区三区| 国产精品一区www在线观看| 国产熟女欧美一区二区| 嫩草影院入口| 色综合色国产| 久久国产亚洲av麻豆专区| 国产精品不卡视频一区二区| 一级毛片电影观看| 国产av一区二区精品久久 | 久久精品国产a三级三级三级| 国产黄频视频在线观看| 国产视频内射| 只有这里有精品99| 麻豆精品久久久久久蜜桃| 国产免费一级a男人的天堂| 18禁在线无遮挡免费观看视频| 日本黄色日本黄色录像| 在线观看美女被高潮喷水网站| 午夜视频国产福利| 色哟哟·www| 国产色爽女视频免费观看| 九九久久精品国产亚洲av麻豆| 欧美 日韩 精品 国产| 自拍偷自拍亚洲精品老妇| 在线 av 中文字幕| 免费看光身美女| 欧美日韩视频精品一区| 欧美国产精品一级二级三级 | 亚洲高清免费不卡视频| 亚洲av不卡在线观看| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 国模一区二区三区四区视频| 国产欧美日韩一区二区三区在线 | 麻豆成人av视频| 国产亚洲午夜精品一区二区久久| 夜夜爽夜夜爽视频| 亚洲av福利一区| 国产色爽女视频免费观看| 精品一区二区免费观看| 老熟女久久久| 亚洲怡红院男人天堂| 亚洲人成网站在线播| 五月伊人婷婷丁香| 99热这里只有精品一区| 少妇丰满av| 极品教师在线视频| 婷婷色综合大香蕉| 免费黄网站久久成人精品| 婷婷色综合大香蕉| 极品教师在线视频| 26uuu在线亚洲综合色| 婷婷色av中文字幕| 少妇人妻一区二区三区视频| 晚上一个人看的免费电影| 亚洲熟女精品中文字幕| 男人狂女人下面高潮的视频| 欧美日韩在线观看h| 久久精品久久久久久久性| 亚洲精品久久久久久婷婷小说| 国产精品一区www在线观看| 两个人的视频大全免费| 国产精品99久久99久久久不卡 | 99精国产麻豆久久婷婷| 午夜福利影视在线免费观看| 黄色欧美视频在线观看| 日本猛色少妇xxxxx猛交久久| 99热这里只有是精品50| av又黄又爽大尺度在线免费看| 新久久久久国产一级毛片| 国内揄拍国产精品人妻在线| 男人狂女人下面高潮的视频| 少妇被粗大猛烈的视频| 午夜精品国产一区二区电影| 九九久久精品国产亚洲av麻豆| 欧美成人午夜免费资源| 高清日韩中文字幕在线| 国产高清三级在线| 亚洲内射少妇av| 秋霞伦理黄片| 久久久精品94久久精品| 国产精品人妻久久久久久| 国产男女内射视频| 久久久久久人妻| 国产高潮美女av| 日韩av免费高清视频| 久久精品久久久久久久性| 亚洲在久久综合| 色婷婷av一区二区三区视频| 国产人妻一区二区三区在| 高清视频免费观看一区二区| 色吧在线观看| 欧美3d第一页| 制服丝袜香蕉在线| 国产黄片视频在线免费观看| 一个人免费看片子| 一区在线观看完整版| 亚洲欧美中文字幕日韩二区| 亚洲中文av在线| 国产精品三级大全| 蜜桃久久精品国产亚洲av| 久久人妻熟女aⅴ| 久久影院123| 久久久a久久爽久久v久久| 一区在线观看完整版| 国产黄色视频一区二区在线观看| 最近中文字幕2019免费版| 亚洲国产精品一区三区| 男女无遮挡免费网站观看| 午夜激情久久久久久久| 久久这里有精品视频免费| 多毛熟女@视频| 18+在线观看网站| 成人影院久久| 伦精品一区二区三区| 天堂8中文在线网| 久久鲁丝午夜福利片| 91精品伊人久久大香线蕉| 国产高清不卡午夜福利| 中文字幕精品免费在线观看视频 | 777米奇影视久久| 美女高潮的动态| 欧美三级亚洲精品| 国产精品爽爽va在线观看网站| 国产乱人偷精品视频| 在线观看免费日韩欧美大片 | 久久久久久久久大av| 中文字幕免费在线视频6| 午夜福利在线在线| 国产精品久久久久久精品电影小说 | 高清视频免费观看一区二区| 多毛熟女@视频| 男女无遮挡免费网站观看| av又黄又爽大尺度在线免费看| 直男gayav资源| 亚洲欧洲国产日韩| 热re99久久精品国产66热6| 成人毛片a级毛片在线播放| 一区二区三区乱码不卡18| 网址你懂的国产日韩在线| 亚洲色图综合在线观看| 内射极品少妇av片p| 亚洲国产欧美人成| 成人美女网站在线观看视频| 亚洲伊人久久精品综合| 嘟嘟电影网在线观看| 久久99蜜桃精品久久| 久久久色成人| 午夜福利在线观看免费完整高清在| 一区二区三区免费毛片| 国产黄片美女视频| 久久6这里有精品| 九草在线视频观看| 亚洲色图综合在线观看| 亚洲精品国产成人久久av| 精品一区二区三区视频在线| 黄色一级大片看看| 青春草亚洲视频在线观看| kizo精华| 精品一区在线观看国产| 国产色爽女视频免费观看| 日韩精品有码人妻一区| 国产精品久久久久成人av| 亚洲精品一二三| 亚洲内射少妇av| 国产色爽女视频免费观看| 这个男人来自地球电影免费观看 | 18禁在线无遮挡免费观看视频| 国产成人午夜福利电影在线观看| 免费不卡的大黄色大毛片视频在线观看| 熟妇人妻不卡中文字幕| 哪个播放器可以免费观看大片| 日韩 亚洲 欧美在线| 少妇高潮的动态图| 久久ye,这里只有精品| 精品久久久噜噜| 亚洲精品乱久久久久久| 国产久久久一区二区三区| 日本与韩国留学比较| 国产美女午夜福利| 亚洲国产高清在线一区二区三| 高清不卡的av网站| 亚洲精品一二三| 国国产精品蜜臀av免费| 国产色婷婷99| xxx大片免费视频| 国产一区有黄有色的免费视频| 国产精品秋霞免费鲁丝片| 欧美极品一区二区三区四区| 日本av免费视频播放| 久久99热6这里只有精品| 亚洲精品视频女| 最近的中文字幕免费完整| 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看| 欧美三级亚洲精品| 中文字幕精品免费在线观看视频 | 精品一区二区免费观看| 久久久久久人妻| 又大又黄又爽视频免费| 熟妇人妻不卡中文字幕| 视频区图区小说| 九九爱精品视频在线观看| 欧美bdsm另类| 狠狠精品人妻久久久久久综合| 久久精品国产自在天天线| 久久久久久久大尺度免费视频| 最黄视频免费看| 国产爽快片一区二区三区| 亚洲国产精品专区欧美| 一个人免费看片子| 色哟哟·www| 日韩在线高清观看一区二区三区| 国产精品精品国产色婷婷| 国产探花极品一区二区| 日韩伦理黄色片| 蜜桃亚洲精品一区二区三区| 国产午夜精品久久久久久一区二区三区| 国产乱来视频区| 亚洲美女黄色视频免费看| 一级av片app| 观看av在线不卡| 美女视频免费永久观看网站| 国产永久视频网站| 99久久综合免费| 九草在线视频观看| 天堂8中文在线网| 日韩欧美一区视频在线观看 | 大片电影免费在线观看免费| 伊人久久国产一区二区| 看非洲黑人一级黄片| 欧美 日韩 精品 国产| 日韩伦理黄色片| 丝瓜视频免费看黄片| 九九久久精品国产亚洲av麻豆| 亚洲图色成人| 美女国产视频在线观看| 尤物成人国产欧美一区二区三区| 国产精品一及| 亚洲第一av免费看| 啦啦啦在线观看免费高清www| 在线观看免费视频网站a站| 一级毛片电影观看| 久久亚洲国产成人精品v| 青青草视频在线视频观看| 伦精品一区二区三区| 国产高清有码在线观看视频| 日本黄色日本黄色录像| 国产成人a∨麻豆精品| 久久久久久久久久久免费av| videossex国产| 亚洲欧美一区二区三区黑人 | 97超碰精品成人国产| 国模一区二区三区四区视频| 亚洲真实伦在线观看| 在线观看人妻少妇| 观看美女的网站| 欧美三级亚洲精品| 国产v大片淫在线免费观看| 啦啦啦在线观看免费高清www| 看免费成人av毛片| 水蜜桃什么品种好| 少妇人妻久久综合中文| 91狼人影院| 夜夜骑夜夜射夜夜干| 欧美日韩视频精品一区| 男女无遮挡免费网站观看| 亚洲aⅴ乱码一区二区在线播放| 午夜福利高清视频| 九九爱精品视频在线观看| 国产视频首页在线观看| 夜夜看夜夜爽夜夜摸| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 岛国毛片在线播放| 少妇人妻精品综合一区二区| 一级毛片 在线播放| 内射极品少妇av片p| 在线免费十八禁| 国产精品秋霞免费鲁丝片| 欧美老熟妇乱子伦牲交| 最近手机中文字幕大全| 观看av在线不卡| 久久久精品免费免费高清| 99久久综合免费| 亚洲真实伦在线观看| 成人美女网站在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 日韩av在线免费看完整版不卡| 亚洲av在线观看美女高潮| 欧美性感艳星| 午夜福利网站1000一区二区三区| 欧美极品一区二区三区四区| 国产精品爽爽va在线观看网站| 一级毛片黄色毛片免费观看视频| 欧美 日韩 精品 国产| 国国产精品蜜臀av免费| 国产有黄有色有爽视频| 妹子高潮喷水视频| 国产成人精品婷婷| 老司机影院毛片| 国产精品久久久久久精品古装| 自拍欧美九色日韩亚洲蝌蚪91 | 久久精品国产自在天天线| 99久久精品一区二区三区| 一本久久精品| 黄片wwwwww| 成人午夜精彩视频在线观看| 亚洲精品国产av成人精品| 3wmmmm亚洲av在线观看| 国产精品人妻久久久久久| 99久国产av精品国产电影| 亚洲欧洲国产日韩| 亚洲性久久影院| 人妻少妇偷人精品九色| 色吧在线观看| 大香蕉久久网| 成年免费大片在线观看| 亚洲内射少妇av| 国产熟女欧美一区二区| 午夜精品国产一区二区电影| 国产精品精品国产色婷婷| 国产高潮美女av|