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

    香溪河流域土地利用變化過程對非點源氮磷輸出的影響

    2016-04-06 06:54:20翟麗梅張富林劉宏斌雷秋良武淑霞華玲玲周繼文
    關(guān)鍵詞:土地利用變化

    崔 超,劉 申,翟麗梅*,張富林,劉宏斌,雷秋良,武淑霞,華玲玲,周繼文

    (1.中國農(nóng)業(yè)科學院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,農(nóng)業(yè)部面源污染控制重點實驗室,北京100081;2.湖北省農(nóng)業(yè)科學院植保土肥研究所,武漢430064;3.湖北省興山縣土壤肥料工作站,湖北宜昌443000)

    ?

    香溪河流域土地利用變化過程對非點源氮磷輸出的影響

    崔超1,劉申1,翟麗梅1*,張富林2,劉宏斌1,雷秋良1,武淑霞1,華玲玲1,周繼文3

    (1.中國農(nóng)業(yè)科學院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所,農(nóng)業(yè)部面源污染控制重點實驗室,北京100081;2.湖北省農(nóng)業(yè)科學院植保土肥研究所,武漢430064;3.湖北省興山縣土壤肥料工作站,湖北宜昌443000)

    摘要:選取三峽庫區(qū)香溪河流域(3150 km2),根據(jù)1990、2000、2010年Landsat TM/ETM遙感影像,在ArcGIS支持下,應(yīng)用景觀特征分析和氮磷輸出系數(shù)模型方法,在分析流域景觀格局轉(zhuǎn)變過程的基礎(chǔ)上,研究了土地利用變化對非點源氮磷輸出的影響。結(jié)果表明:1990—2000年土地利用變化較為緩和,2000—2010年土地利用轉(zhuǎn)變較為劇烈,變化面積占到總面積的4.3%,為前10年變化量占比的3倍之多。1990—2000年和2000—2010年,土地利用轉(zhuǎn)變量最為明顯的均為林地轉(zhuǎn)旱地和旱地轉(zhuǎn)林地;從單位時間土地利用變化率來看,1990—2000年旱地變化最為劇烈,2000—2010年居民地變化最為劇烈;從土地利用相對動態(tài)度來看,1990—2000年和2000—2010年旱地均最高。1990—2000年土地利用的變化對非點源總氮(TN)和總磷(TP)輸出影響較小,TN增加1.57 t.a-1,TP減少0.073 t.a-1;2000—2010年土地利用的變化顯著降低了TN、TP的輸出量,凈值分別為78.5、6.1 t.a-1。土地利用轉(zhuǎn)變方式對TN、TP的負荷影響不同,旱地轉(zhuǎn)變?yōu)榱值貢rTN負荷表現(xiàn)為消減,林地轉(zhuǎn)變?yōu)楹档貢rTN負荷增加;水田轉(zhuǎn)變?yōu)榱值貢rTP負荷減少,旱地轉(zhuǎn)變?yōu)樗飼rTP負荷增加。在輸出系數(shù)與各土地利用類型面積關(guān)系的建立中,使用土地利用狀態(tài)量變化面積不能真實計算出其對非點源氮磷負荷輸出量的影響,通過土地利用過程量面積的變增才能真實反映土地利用變化導(dǎo)致的非點源氮磷負荷輸出量。

    關(guān)鍵詞:香溪河流域;土地利用變化;非點源;輸出系數(shù);狀態(tài)量和過程量

    崔超,劉申,翟麗梅,等.香溪河流域土地利用變化過程對非點源氮磷輸出的影響[J].農(nóng)業(yè)環(huán)境科學學報,2016,35(1):129-138.

    CUI Chao,LIU Shen,ZHAI Lj-mej,et a1. Effect of 1and use/cover changes on njtrogen and phosphorus 1osses vja non-pojnt source pathway jn Xjangxj Rjver Basjn[J]. Journal of Agro-Environment Science,2016,35(1):129-138.

    土地利用變化狀況是影響全球和區(qū)域環(huán)境?變的重要指標[1-2]。非點源污染氮磷負荷輸出與土地利用方式密切相關(guān),土地利用變化與水量、水質(zhì)之間存在非常顯著的相關(guān)性,土地利用方式轉(zhuǎn)變不僅影響地表徑流、基流等水文過程,更影響土壤侵蝕、泥沙以及營養(yǎng)物的運移[3-5]。不同土地利用類型的變化過程反映了土地利用格局在時間?列上的變化,每種土地利用類型的過程變化包括了時段內(nèi)過程增長量和過程減少量[6-7],然而某一土地利用類型狀態(tài)上的變化,即用t1時刻格局狀態(tài)量與t2時刻格局狀態(tài)量的差值并不能很好地反映土地利用類型的轉(zhuǎn)變過程,土地利用動態(tài)變化矩陣則能夠?晰反映每種土地利用類型t1到t2時段內(nèi)過程增長量與過程減少量(圖1)。土地利用空間信息在時間?列上的動態(tài)轉(zhuǎn)變過程,是土地類型轉(zhuǎn)變的內(nèi)在結(jié)果,其過程的轉(zhuǎn)變直接?變了非點源過程量的輸出,最終影響非點源狀態(tài)量的表達。

    圖1 土地利用變化過程量與狀態(tài)量Fjgure 1 Process and state of 1and use/cover changes

    非點源污染的時空不確定性、遷移過程的高度非線性以及過程降解的復(fù)雜性等特征[8-9],增加了其排放量估算、來源解析的難度。根據(jù)輸出系數(shù)模型可以直接建立土地利用與非點源氮磷負荷輸出之間的關(guān)系,從而避開非點源污染發(fā)生的復(fù)雜過程,僅根據(jù)土地利用現(xiàn)狀或變化量估算或預(yù)測非點源氮磷負荷輸出,該方法簡單有效,適合應(yīng)用于估算缺乏大量監(jiān)測數(shù)據(jù)區(qū)域非點源氮磷輸出負荷。目前大多數(shù)研究者[10-14]只是通過建立輸出系數(shù)與土地利用狀態(tài)量之間的關(guān)系,匯算由當年土地利用狀態(tài)下產(chǎn)生的非點源氮磷輸出量,并沒有綜合考慮土地利用轉(zhuǎn)變過程導(dǎo)致的非點源氮磷輸出的響應(yīng)關(guān)系,不能真實反映土地利用變化對非點源氮磷負荷輸出的影響。因此,?用土地利用動態(tài)變化矩陣,建立輸出系數(shù)與土地利用凈變化量之間的關(guān)系,探討土地利用轉(zhuǎn)變過程對非點源氮磷負荷輸出的影響,對流域非點源污染的管控具有重要意義。

    本研究基于三峽庫區(qū)香溪河流域1990年、2000年和2010年三個年度遙感解譯土地利用數(shù)據(jù),運用景觀特征分析和氮磷輸出系數(shù)模型方法,探討香溪河流域景觀格局轉(zhuǎn)變過程和對非點源氮磷負荷輸出的影響,進而為三峽庫區(qū)流域非點源污染的綜合治理奠定基礎(chǔ)。

    1 材料與方法

    1.1研究區(qū)概況

    香溪河流域發(fā)源于湖北省宜昌市神農(nóng)架林區(qū),流域總面積為3150 km2,干流全長95 km,途?興山縣(約78 km)至秭歸縣,由河口處匯入三峽大壩,是三峽庫區(qū)壩首的第一大支流。其主要包括高嵐河、古夫河和南陽河三大支流,處于上游的南陽河和古夫河在響灘匯流始稱香溪,南流14 km峽口鎮(zhèn)匯高嵐河至游家河入秭歸境內(nèi)[15]。流域?qū)贅?gòu)造地貌,高程80~3100 m,地形起伏較大,均系高山半高山區(qū)(圖2)。

    香溪河流域是三峽庫區(qū)典型的農(nóng)林復(fù)合小流域,土地利用類型主要有林地、耕地。土壤類型主要為黃棕壤和石灰土。年均降水量為900~1200 mm,主要集中在4—10月[16]。流域水土流失較嚴重,據(jù)統(tǒng)計香溪河興山縣范圍水土流失面積達1 122.5 km2,占全縣總面積的61%,年均侵蝕模數(shù)為6488 t.km-2,屬強度侵蝕區(qū)[16]。

    1.2研究方法

    1.2.1數(shù)據(jù)來源及處理

    研究區(qū)土地利用圖來自中國科學院資源環(huán)境科學數(shù)據(jù)中心,比例尺為1∶10萬。土地利用圖為1990年、2000年和2010年各期Landsat TM/ETM遙感影像。遙感影像通過人工目視解譯生成土地利用信息,劃分水田、旱地、草地、林地、園地、水域和居民地等共7種土地利用類型(圖3)。為構(gòu)建輸出系數(shù)模型,通過ArcGIS對矢量文件中的土地利用信息柵格化,柵格規(guī)格為25 m×25 m。土地利用類型面積通過柵格數(shù)量及其柵格面積,匯算所得。

    圖2 研究區(qū)概況圖Fjgure 2 Map of studjed area

    圖3 1990、2000、2010年土地利用現(xiàn)狀圖Fjgure 3 Land use/cover jn 1990,2000,and 2010

    ?用ArcGIS空間數(shù)據(jù)分析中疊加分析模塊,對不同年份間,空間同一點位土地利用類型的柵格屬性進行空間信息識別,確立1990—2000年間和2000—2010年間土地利用轉(zhuǎn)變的柵格矩陣。

    1.2.2數(shù)據(jù)分析方法

    1.2.2.1景觀動態(tài)分析方法

    土地的單位時間土地利用變化率是以土地利用類型面積為基礎(chǔ),關(guān)注研究時段內(nèi)某用地類型面積變化的結(jié)果[18]。其計算公式為:

    式中:Ai為第i種土地類型單位時間土地利用變化率;當t的單位設(shè)定為年時,模型結(jié)果表示該土地利用類型的年變化率。

    單一動態(tài)度可直觀地反映某一土地利用類型相對本身的變化劇烈程度與速度,也易于通過類型間的比較反映不同類型間變化的差異。為了反映各土地利用類型的變化在區(qū)域變化中的貢獻大小,衡量某土地類型的變化是否成為區(qū)域土地變化的主要類型,引進土地相對動態(tài)度概念來表示土地利用的變化速率。其計算公式為:

    式中:Bi為第i種土地利用相對動態(tài)度;Pi為研究初期第i類用地類型占研究區(qū)土地總面積的比率。

    土地相對動態(tài)度的值越大,表明該用地類型的變化在研究區(qū)變化中的貢獻越大。

    1.2.2.2非點源污染排放分析方法輸出系數(shù)模型[11]表達式:

    式中:Lj為污染物j在該流域的總負荷量,kg.a-1;n為流域中土地利用類型的種類;Eij為污染物j在流域第i種土地利用類型中的輸出系數(shù),kg.km-2.a-1;Ai為第i種土地利用類型的面積,km2。

    1.2.3輸出系數(shù)選取

    非點源污染的形成主要受流域地形、氣候、水文、土地利用類型、農(nóng)田管理措施、人類活動等因素的影響,因此,合理輸出系數(shù)的確定是輸出系數(shù)模型構(gòu)建的關(guān)鍵。鑒于國內(nèi)許多研究人員已用輸出系數(shù)模型做了大量研究[9-13],方法及系數(shù)的確立已相對成熟??紤]到丹江口水庫研究區(qū)水文、氣候、地理位置、植被覆蓋等多因素與本研究較接近,引用前人研究使用過的輸出系數(shù)[11](表1)。

    表1 各土地利用類型輸出系數(shù)Tab1e 1 Export coeffjcjents of varjous 1and use/cover types

    2 結(jié)果與分析

    2.1不同年度土地利用類型結(jié)構(gòu)及空間變化特征

    從三個年份土地利用狀況來看(圖4),香溪河流域主要土地利用類型為林地,三個年度均占85%以上;耕地(水田、旱地)和園地面積小,不足流域總面積的8%。1990—2000年10年間各土地利用類型的面積基本保持不變。2000—2010年土地利用變化明顯,果園、水域、居民地增長最多,分別為39.7%、91.7%和 251.0%;水田和旱地分別減少2.8%和32.9%。

    圖4 1990年、2000年和2010年各土地利用類型占比情況Fjgure 4 Percentages of varjous 1and use/cover types jn 1990,2000 and 2010

    表2 1990—2000年香溪河流域土地利用轉(zhuǎn)移矩陣(km2)Tab1e 2 Transfer matrjx of 1and use/cover jn Xjangxj Rjver basjn from 1990 to 2000(km2)

    由1990—2000年香溪河流域土地利用轉(zhuǎn)移矩陣得出(表2),土地轉(zhuǎn)變面積較小,約占總量的1.4%。林地和旱地轉(zhuǎn)變數(shù)量最多,分別占變化總量的47.3%和30.8%;水域、旱地和水田參與其他土地利用轉(zhuǎn)變率最高,轉(zhuǎn)變量分別占1990年水域、旱地和水田的8.3%、7.5%和6.1%。在轉(zhuǎn)變過程中,水田主要轉(zhuǎn)變?yōu)楹档睾土值?;旱地主要轉(zhuǎn)變?yōu)榱值亍⒉莸睾退?;林地主要轉(zhuǎn)變?yōu)楹档亍⒉莸睾退铩?/p>

    據(jù)2000—2010年香溪河流域土地利用轉(zhuǎn)移矩陣來看(表3),土地利用變化更加劇烈,變化量占總量的4.3%,為前10年的3倍之多。林地和旱地依然變化最大,分別占變化總量的57.1%和26.4%;水田和旱地轉(zhuǎn)變率最高,轉(zhuǎn)變量分別占2000年水田和旱地總量的24.4%和43.7%。轉(zhuǎn)變過程中,水田主要轉(zhuǎn)變?yōu)楹档睾土值兀缓档刂饕優(yōu)榱值亍⑺?;林地主要變?yōu)楹档睾退颉?/p>

    表3 2000—2010年香溪河流域土地利用轉(zhuǎn)移矩陣(km2)Tab1e 3 Transfer matrjx of 1and use/cover jn Xjangxj Rjver basjn from 2000 to 2010(km2)

    2.2不同時間段景觀變化動態(tài)度特征

    通過表4和表5可知,2000—2010年較1990—2000年變化劇烈程度增加,凈變化量最明顯的均為林地和旱地。從單位時間土地利用變化率來看,1990—2000年流域變化最為劇烈的是旱地,變化率為0.076%;2000—2010年為居民地,變化率為25.1%。其中,水域單位時間土地利用變化率,均位居第二,說明水域面積呈不斷擴大趨勢。從土地利用相對動態(tài)度來看,1990—2000年和2000—2010年旱地的相對動態(tài)度均最高,分別為0.004 3%和-0.186 0%,說明旱地的對流域土地利用變化量最大。

    表4 1990—2000年各土地利用類型土地動態(tài)度分析Tab1e 4 Land dynamjc ana1ysjs of varjous 1and use/cover types from 1990 to 2000

    表5 2000—2010年各土地利用類型土地動態(tài)度分析Tab1e 5 Land dynamjc ana1ysjs of varjous 1and use/cover types from 2000 to 2010

    2.3土地利用變化對非點源排放的影響

    2.3.1不同年度土地利用變化對氮磷排放的影響

    從表6來看,整個流域由土地利用變化所造成的非點源TN和TP負荷從20世紀80年代末期到2000年,變化基本維持平穩(wěn),TN略增,TP略減,凈量分別為1.57、-0.073 t.a-1。從2000年到2010年TN、TP負荷變化顯著(均減?。?,凈值分別為-78.5、-6.1 t.a-1。

    表6 1990、2000、2010年各土地利用類型狀態(tài)量輸出負荷(t.a-1)Tab1e 6 State output 1oads of varjous 1and use/cover types jn 1990,2000 and 2010(t.a-1)

    各土地利用類型轉(zhuǎn)變導(dǎo)致的非點源氮磷狀態(tài)量,是根據(jù)不同土地利用類型最終狀態(tài)面積與輸出系數(shù)之間的關(guān)系計算得出(表7)。從非點源氮磷狀態(tài)量輸出來看,1990—2000年氮磷的增加量和減少量基本維持平衡,2000—2010年非點源氮磷輸出總量減少,且TN減少更為明顯。2000—2010年與1990—2000年相比,氮磷輸出變化更為劇烈。從TN和TP來看,旱地減少量最大,分別減少了112.8、4.64 t.a-1;林地增加量最大,增加量分別為27.2、0.75 t.a-1;水域?qū)N、TP輸出無影響。

    非點源氮磷過程輸出量,是通過?用土地利用動態(tài)變化矩陣,建立相對系數(shù)(兩土地利用類型輸出系數(shù)差)與土地利用凈變化量之間的關(guān)系計算得出(表8),能夠反映不同土地利用類型轉(zhuǎn)變所隱含的內(nèi)部非點源氮磷負荷輸出狀況。從過程量來看,1990—2000年和2000—2010年氮磷變化量最大的土地利用類型均為水田、旱地、林地。其中,水田TN、TP均減少,2000—2010年與1990—2000年相比,TN、TP減少量均增加了3倍以上;旱地中TN持續(xù)減少,從-15.2 t. a-1到-81.3 t.a-1,增加5倍以上,TP持續(xù)增多,從0.32 t.a-1到9.73 t.a-1,增加了30倍;林地在兩個年度氮磷輸出量均增加,且過程量增加在年度間變化相對緩和,TN為23.7、27.7 t.a-1,TP為5.57、8.56 t.a-1。

    表7 各土地利用轉(zhuǎn)變氮磷狀態(tài)量匯算表(t.a-1)Tab1e 7 Amount of njtrogen and phosphorous status from varjous 1and use/cover transformatjon(t.a-1)

    表8 各土地利用轉(zhuǎn)變氮磷過程量匯算表(t.a-1)Tab1e 8 Amount of njtrogen and phosphorous from varjous 1and use/cover transformatjon by process(t.a-1)

    2.3.2土地利用變化過程對氮磷排放的影響特征

    圖5詳細表述了1990—2000年和2000—2010年兩個年限周期,某一土地利用類型在其他土地利用轉(zhuǎn)變?yōu)楝F(xiàn)有土地利用類型和現(xiàn)有土地利用轉(zhuǎn)變?yōu)槠渌恋乩眠^程變化中,在該土地利用類型上所發(fā)生的非點源氮磷負荷的過程凈量。正軸表示對某土地利用與其他不同土地利用間參與轉(zhuǎn)變過程,并起到增加氮磷負荷作用的輸出凈量;負軸表示對某土地利用與其他不同土地利用參與轉(zhuǎn)變過程,并起到消減氮磷負荷作用的輸出凈量。1990—2000年各土地利用變化間TN、TP負荷凈量輸出總體相對緩和,林地、旱地、水田的土地利用動態(tài)變化對TN負荷的貢獻較大,水田、林地對TP負荷的貢獻較大;2000—2010年較1990—2000年,由土地類型轉(zhuǎn)變所導(dǎo)致的TN、TP負荷凈量變化更為劇烈,對TN負荷貢獻較大的土地利用方式為旱地、林地、水田,TP主要是水田、旱地、林地。

    圖5 各土地利用過程量轉(zhuǎn)變TN、TP輸出凈量結(jié)構(gòu)圖Fjgure 5 Structure chart for 1and use/cover change effects on TN and TP net 1osses by process

    續(xù)圖5 各土地利用過程量轉(zhuǎn)變TN、TP輸出凈量結(jié)構(gòu)圖Contjnue fjgure 5 Structure chart for 1and use/cover change effects on TN and TP net 1osses by process

    從各土地利用轉(zhuǎn)變方式對TN、TP的負荷貢獻來看,不同土地利用轉(zhuǎn)變對總量增長、消減的貢獻能力不同。1990—2000年水田轉(zhuǎn)變?yōu)榱值亍⒑档剞D(zhuǎn)變?yōu)榱值貢r對TN主要表現(xiàn)為消減,林地轉(zhuǎn)變旱地、林地轉(zhuǎn)變水田時TN負荷增加;而對TP負荷而言,水田轉(zhuǎn)變林地、水田轉(zhuǎn)變旱地表現(xiàn)為消減,而林地轉(zhuǎn)變水田、旱地轉(zhuǎn)變水田表現(xiàn)為增加。2000—2010年內(nèi)部土地利用轉(zhuǎn)換間對氮磷負荷貢獻量相比1990—2000年變化更加劇烈,起主要負荷貢獻作用的土地轉(zhuǎn)化方式與1990—2000年的類似。旱地轉(zhuǎn)變林地、水田轉(zhuǎn)變林地對TN負荷表現(xiàn)為消減,林地轉(zhuǎn)變旱地、林地轉(zhuǎn)變水田、旱地轉(zhuǎn)變水田時表現(xiàn)為增加;對TP來說,水田轉(zhuǎn)變林地、水田轉(zhuǎn)變旱地、旱地轉(zhuǎn)變林地時表現(xiàn)為消減,旱地轉(zhuǎn)變水田、林地轉(zhuǎn)變水田時表現(xiàn)為增加。

    3 討論

    3.1土地利用景觀格局變化特點及成因

    本文通過香溪河流域土地利用結(jié)構(gòu)、景觀動態(tài)變化、空間轉(zhuǎn)移分析得出,20年間,2000—2010年土地利用變化最為劇烈,林地、旱地凈量變化最大,居民地單位時間土地利用變化率最高,旱地的相對動態(tài)度最大。整個流域1990—2000年退耕還林、毀林復(fù)墾和水域擴增面積分別為14.6、16.1、0.6 km2,2000—2010年分別為72.7、22.8、8.9 km2,20年間退耕還林、水域擴增顯著。

    分析發(fā)生以上景觀格局轉(zhuǎn)變特點的成因,主要為三峽大壩1993年開始修建,2009年正式完工[19],2000年之前為大壩建設(shè)前期,土地格局變化小,2000年以后隨著大壩建造工程的推進,大壩蓄水后主要支流形成庫灣,水域面積擴增,流域周邊肥沃的農(nóng)田被淹沒,據(jù)統(tǒng)計[19]三峽庫區(qū)大約有240 km2農(nóng)用地被淹沒,加之大量河谷居民遷移,新的農(nóng)用地急需開墾,人地矛盾日益凸顯[20],且國家1999年開始試點“退耕還林”政策落實,2003年正式頒布了《退耕還林例》,退耕還林進入全面依法實施的階段,生態(tài)退耕減少的耕地主要為坡耕地,且主要分布在坡度大于25°的丘陵地區(qū)[7]。以上復(fù)雜原因造成香溪河流域土地利用轉(zhuǎn)變格局加劇。

    3.2土地利用變化對非點源氮磷輸出的影響

    在非點源氮磷負荷估算上,由英國學者Johns等[21]提出的輸出系數(shù)模型為非機理的“黑箱”模型,基于數(shù)據(jù)要求簡單、氮磷負荷估算高效等優(yōu)勢,被國內(nèi)學者廣泛應(yīng)用,其最大的缺陷在于忽略了水文過程對污染物的自凈功能。該模型在國內(nèi)的應(yīng)用中,研究者更多是針對其影響因素帶來的不確定性進行優(yōu)化[22-23],提高輸出系數(shù)模型負荷估算的準確性。

    劉瑞民等[10]運用輸出系數(shù)模型,根據(jù)1970年、1980年、1990年和2000年土地利用狀態(tài)量,對長江上游的各土地利用導(dǎo)致的非點源氮磷負荷進行了估算;方怒放等[11]以1990、2000、2007年遙感影像解譯的土地利用圖為基礎(chǔ),使用其三期土地利用狀態(tài)量,運用輸出系數(shù)模型對庫區(qū)非點源污染進行空間模擬。但是研究者只是通過建立輸出系數(shù)與土地利用面積狀態(tài)量之間的關(guān)系,匯算不同土地利用引起的非點源氮磷負荷輸出量,并沒有考慮不同年份間各土地利用轉(zhuǎn)化過程發(fā)生的實際非點源氮磷輸出量。從實際情況來看,各土地利用面積狀態(tài)量的變增并不能真實反映出其對非點源氮磷負荷輸出的影響,匯算不同時間下土地利用過程的變化量才能使用輸出系數(shù)模型準確計算出各土地利用轉(zhuǎn)變的實際非點源氮磷發(fā)生量。從土地利用面積狀態(tài)量和過程量所產(chǎn)生的非點源氮磷負荷匯總結(jié)果看出,由各土地利用變化導(dǎo)致的非點源氮磷輸出總量相同:1990—2000年TN、TP分別為1.567、-0.073 t.a-1;2000—2010年分別為-78.584、-6.074 t.a-1。但各土地利用類型在狀態(tài)量和過程量下的非點源氮磷負荷輸出在貢獻方向和貢獻量上均表現(xiàn)出顯著差異,各土地利用類型實際的非點源氮磷負荷輸出量并不是其兩狀態(tài)量匯算結(jié)果之差。在土地利用狀態(tài)量的變化中,水域的非點源氮磷輸出顯示為零,但在實際的變化中水域也參與了向其他土地利用轉(zhuǎn)化的過程,對非點源氮磷負荷的輸出量不應(yīng)該為零;在林地狀態(tài)量的變化中,結(jié)果顯示1990—2000年非點源氮磷凈輸出負荷減少,2000—2010年表現(xiàn)為增加,且2000—2010年較1990—2000年變化更為顯著。但從過程量匯算結(jié)果來看,林地在兩個時段均對非點源氮磷負荷輸出表現(xiàn)增加,且變化量差異較小,變化劇烈程度相當。

    4 結(jié)論

    1990—2000年不同土地利用類型間轉(zhuǎn)變相對緩和,2000—2010年土地利用間轉(zhuǎn)變較為劇烈,空間位置變化量為前10年變化量占比的3倍之多,土地利用轉(zhuǎn)變主要發(fā)生在林地、旱地、水田、水域之間。由于土地利用所造成的非點源TN、TP從1990到2010年基本上呈逐漸下降的趨勢,這種變化趨勢和土地利用變化緊密相關(guān)。在輸出系數(shù)與各土地利用類型面積關(guān)系的建立中,使用土地利用狀態(tài)量變化面積不能真實計算出其對非點源氮磷負荷輸出量,通過土地利用過程量面積的變增才能真實反映出土地利用變化導(dǎo)致的非點源氮磷負荷輸出量。因此,討論土地利用時間?列上的動態(tài)演變對非點源氮磷負荷輸出的影響,應(yīng)使用各土地利用轉(zhuǎn)變過程發(fā)生量換算各土地利用變化導(dǎo)致的非點源氮磷負荷輸出量。

    參考文獻:

    [1] Ster1jng S M,Ducharne A,Po1cher J. The jmpact of g1oba1 1and-cover change on the terrestrja1 water cyc1e[J]. Nature Climate Change,2012,3 (4):385-390.

    [2] Ljndqujst E J,Hansen M C,Roy D P,et a1. The sujtabj1jty of decada1 jmage data sets for mappjng tropjca1 forest cover change jn the Democratjc Repub1jc of Congo:Imp1jcatjons for the g1oba1 1and survey[J]. International Journal of Remote Sensing,2008,29(23-24):7269-7275.

    [3] Ljn Y P,Hong N M,Wu P J,et a1. Impacts of 1and use change scenarjos on hydro1ogy and 1and use patterns jn the Wu-Tu watershed jn Northern Tajwan[J]. Landscape and Urban Planning,2007,80:111-126.

    [4] Bjeger K,Ho¨rmann G,F(xiàn)ohrer N. The jmpact of 1and use change jn the Xjangxj Catchment(Chjna)on water ba1ance and sedjment transport[J]. Regional Environmental Change,2015,15:485-498.

    [5]王秀娟,劉瑞民,宮永偉,等.香溪河流域土地利用格局演變對非點源污染的影響研究[J].環(huán)境工程學報,2011,5(5):1194-1200. WANG Xju-juan,LIU Ruj-mjn,GONG Yong-wej,et a1. Sjmu1atjon of the effect of 1and use/cover change on non-pojnt source po11utjon 1oad jn Xjangxj Rjver watershed[J]. Chinese Journal of Environmental Engineering,2011,5(5):1194-1200.

    [6]趙安玖,胡庭興,羅龍海,等.典型退耕還林區(qū)景觀格局變化特征[J].林業(yè)資源管理,2008,12(6):85-91. ZHAO An-jju,HU Tjng-xjng,LUO Long-haj,et a1. Characterjstjcs of 1andscape pattern changes jn typjca1 regjon of crop1and conversjon to forest 1and[J]. Forest Resources Management,2008,12(6):85-91.

    [7]劉紀遠,匡文慧,張增祥,等. 20世紀80年代末以來中國土地利用變化的基本特征與空間格局[J].地理學報,2014,69(1):3-14. LIU Jj-yuan,KUANG Wen-huj,ZHANG Zeng-xjang,et a1. Spatjotempora1 characterjstjcs,patterns and causes of 1and use changes jn Chjna sjnce the 1ate 1980s[J]. Acta Geographica Sinica,2014,69(1):3-14.

    [8]張召喜.基于SWAT模型的鳳羽河流域農(nóng)業(yè)面源污染特征研究[D].北京:中國農(nóng)業(yè)科學院,2013. ZHANG Zhao-xj. Study on charaterjstjcs of agrjcu1tura1 non-pojnt source po11utjon jn Fengyu Rjver Basjn based on SWAT mode1[D]. Bejjjng:the Chjnese Academy of Agrjcu1tura1 Scjences,2013.

    [9]李文超.鳳羽河流域農(nóng)業(yè)面源污染負荷估算及關(guān)鍵區(qū)識別研究[D].北京:中國農(nóng)業(yè)科學院,2014. LI Wen-chao. Eva1uatjng the 1oads of agrjcu1tura1 non-pojnt source po11utjon and jdentjfyjng crjtjca1 source areas jn Fengyu Basjn[D]. Bejjjng:the Chjnese Academy of Agrjcu1tura1 Scjences,2013.

    [10]劉瑞民,楊志峰,丁曉雯,等.土地利用/覆蓋變化對長江上游非點源污染影響研究[J].環(huán)境科學,2006,27(12):2407-2414. LIU Ruj-mjn,YANG Zhj-feng,DING Xjao-wen,et a1. Effect of 1and use/cover change on po11utjon 1oad of non-Pojnt source jn upper reach of Yangtze Rjver Basjn[J]. Environmental Science,2006,27(12):2407-2414.

    [11]方怒放,史志華,李璐.基于輸出系數(shù)模型的丹江口庫區(qū)非點源污染時空模擬[J].水生態(tài)學雜志,2011,32(4):7-12. FANG Nu-fang,SHI Zhj-hua,LI Lu. App1jcatjon of export coeffjcjentmode1 jn sjmu1atjng po11utjon 1oad ofnon-pojnt source jn Danjjangkou Reservojr Area[J]. Journal of Hydroecology,2011,32(4):7-12.

    [12]龍?zhí)煊?,梁常德,李繼承,等.基于SLURP模型和輸出系數(shù)法的三峽庫區(qū)非點源氮磷負荷預(yù)測[J].環(huán)境科學學報,2008,28(3):574-581. LONG Tjan-yu,LIANG Chang-de,LI Jj-cheng,et a1. Forecastjng the po11utjon 1oad of non-pojnt sources jmported to the Three Gorges Reservojr[J]. Acta Scientiae Circumstantiae,2008,28(3):574-581.

    [13]張立坤,香寶,胡鈺,等.基于輸出系數(shù)模型的呼蘭河流域非點源污染輸出風險分析[J].農(nóng)業(yè)環(huán)境科學學報,2014,33(1):148-154. ZHANG Lj-kun,XIANG Bao,HU Yu,et a1. Rjsk Assessment of nonpojnt source po11utjon jn Hu1an Rjver Basjn usjng an output coeffjcjent mode1[J]. Journal of Agro-Environment Science,2014,33(1):148-154.

    [14]劉瑞民,沈珍瑤,丁曉雯,等.應(yīng)用輸出系數(shù)模型估算長江上游非點源污染負荷[J].農(nóng)業(yè)環(huán)境科學學報,2008,27(2):677-682. LIU Ruj-mjn,SHEN Zhen-yao,DING Xjao-wen,et a1. App1jcatjon of export coeffjcjent mode1 jn sjmu1atjng po11utjon 1oad of non-pojnt source jn upper reach of Yangtze Rjver Basjn[J]. Journal of Agro-Environment Science,2008,27(2):677-682.

    [15]崔超,劉申,翟麗梅,等.興山縣香溪河流域農(nóng)業(yè)源氮磷排放估算及時空特征分析[J].農(nóng)業(yè)環(huán)境科學學報,2015,34(5):937-946. CUI Chao,LIU Shen,ZHAI Lj-mej,et a1. Estjmatjon of agrjcu1tura1 non-pojnt sources of njtrogen and phosphorus djscharges jn Xjangxj Rjver Basjn based on po11utant djscharge coeffjcjent approach and ana1ysjs of thejr tempora1 and spatja1 characterjstjcs[J]. Journal of Agro-Environment Science,2015,34(5):937-946.

    [16]宋林旭,劉德富,肖尚斌.三峽庫區(qū)香溪河流域非點源營養(yǎng)鹽輸出變化的試驗研究[J].長江流域資源與環(huán)境,2011,20(8):990-996. SONG Ljn-xu,LIU De-fu,XIAO Shang-bjn. Experjmenta1 study on non-pojnt source nutrjent output from Xjangxj Basjn jn Three Gorges Reservojr[J]. Resources and Environment in the Yangtze Basin,2011,20(8):990-996.

    [17]張超.非點源污染模型研究及其在香溪河流域的應(yīng)用[D].北京:?華大學,2008. ZHANG Chao. Djstrjbuted non-pojnt sources po11utjon mode1jng and jts app1jcatjon jn Xjangxj watershed[D]. Bejjjng:Tsjnghua Unjversjty,2008.

    [18]施志國.沽源縣土地利用景觀格局變化研究[D].北京:北京林業(yè)大學,2009. SHI Zhj-guo. Study on the dynamjc changes of 1and use and 1andscape patternjn Guyuan County[D].Bejjjng:Bejjjng Forestry Unjversjty,2009. [19] Zhang J X,Lju Z J,Sun X X. Changjng 1andscape jn the Three Gorges Reservojr Area of Yangtze Rjver from 1977 to 2005:Land use/1and cover,vegetatjon cover changes estjmated usjng mu1tj-source sate11jte data[J]. International Journal of Applied Earth Observation and Geoinformation,2009(6):403-412.

    [20]徐昔保,楊桂山,李恒鵬.三峽庫區(qū)1980—2005年農(nóng)業(yè)用地氮平衡時空變化研究[J].環(huán)境科學,2009,30(8):2228-2233. XU Xj-bao,YANG Guj-shan,LI Heng-peng. Spatjotempora1 changes of njtrogen ba1ance jn 1980—2005 for agrjcu1tura1 1and jn Three Gorges Reservojr Area[J]. Environmental Science,2009,30(8):2228-2233.

    [21] Johnes P J. Eva1uatjon and management of the jmpact of 1and use change on the njtrogen and phosphorus 1oad de1jvered to surface waters:The export coeffjcjent mode1jng approach[J]. Journal of Hydrology,1996,183:323-349.

    [22]蔡明,李懷恩,莊詠濤.?進的輸出系數(shù)法在流域非點源污染負荷估算中的應(yīng)用[J].水利學報,2004(7):40-45. CAI Mjng,LI Huaj-en,ZHUANG Yong-tao. App1jcatjon of modjfjed export coeffjcjent method jn po11utjng 1oad estjmatjon of non-pojnt source po11utjon[J]. Journal of Hydraulic Engineering,2004(7):40-45.

    [23]丁曉雯,沈珍瑤,劉瑞民,等.基于降雨和地形特征的輸出系數(shù)模型?進及精度分析[J].長江流域資源與環(huán)境,2008,17(2):306-309. DING Xjao-wen,SHEN Zhen-yao,LIU Ruj-mjn,et a1. Improved export coeffjcjent mode1 consjderjng precjpjtatjon as we11 as terrajn and jts accuracy ana1ysjs[J]. Resources and Environment in the Yangtze Basin,2008,17(2):306-309.

    Effect of land use/cover changes on nitrogen and phosphorus losses via non-point source pathway in Xiangxi River Basin

    CUI Chao1,LIU Shen1,ZHAI Lj-mej1*,ZHANG Fu-1jn2,LIU Hong-bjn1,LEI Qju-1jang1,WU Shu-xja1,HUA Ljng-1jng1,ZHOU Jj-wen3
    (1.Instjtute of Agrjcu1tura1 Resources and Regjona1 P1annjng,Chjnese Academy of Agrjcu1tura1 Scjences,Key Laboratory of Nonpojnt Source Po11utjon Contro1,Mjnjstry of Agrjcu1ture,Bejjjng 100081,Chjna;2.Instjtute of P1ant Protectjon,Soj1 and Fertj1jzer Scjences,Hubej Academy of Agrjcu1tura1 Scjences,Wuhan 430064,Chjna;3.Xjngshan Soj1 and Fertj1jzer Statjon of Hubej Provjnce,Yjchang 443000,Chjna)

    Abstract:Land use changes great1y jnf1uence non-pojnt source njtrogen and phosphorus outputs. Here,the effect of 1and use/cover changes on non-pojnt source outputs of njtrogen and phosphorus was studjed based on remote sensjng jmages(Landsat TM/ETM)jn 1990,2000 and 2010 usjng the methods of characterjstjcs of 1andscape ana1ysjs,export coeffjcjent mode1,and ArcGIS jn the Xjangxj Rjver basjn(3150 km2). Land use/cover changes were re1atjve1y moderate from 1990 to 2000,but was severe from 2000 to 2010. The area that 1and uses had changed from one jnto another between 2000 to 2010 accounted for 4.3% of the tota1 area,whjch was three tjmes that from 1990 to 2000. Durjng 1990 to 2000 and 2000 to 2010,the greatest 1and use/cover type changes were from forest 1and to dry 1and or from dry 1and to forest 1and. The annua1 change rates of 1and use/cover was drastjc for dry1and from 1990 to 2000 whj1e resjdentja1 1and from 2000 to 2010. The bjggest dynamjc change area of 1and use/cover changes was dry1and for both 1990 to 2000 and 2000 to 2010. Durjng 1990 to 2000,the 1andbook=161,ebook=135use/cover changes had 1jtt1e effect on tota1 njtrogen(TN)and tota1 phosphorus(TP)1oads from non-pojnt source. Durjng 2000 to 2010,however,1and use/cover changes sjgnjfjcant1y reduced the 1oads of TN and TP,whjch was 78.5 and 6.1 t.a-1respectjve1y. The contrjbutjon of djfferent 1and use/cover changes to TN and TP 1oads varjed. The reductjon of TN 1oad majn1y happened when dry 1and changed to forest1and. The conversjon from paddy fje1d to forest1and decreased TP 1oad,whereas turnjng dry1and to paddy fje1d jncreased TP 1oad. In the estab1jshment of the quantjtatjve re1atjonshjp between export coeffjcjents and djfferent 1and use/cover type area,the 1and use/cover changes ca1cu1ated by statjc area cou1d not ref1ect jts jmpact on the TN and TP 1oads,whj1e the 1and use/cover changes ca1cu1ated by process area was a good jndjcator of TN and TP 1oads.

    Keywords:Xjangxj Rjver basjn;1and use/cover change;non-pojnt source;export coeffjcjent;state and process va1ue

    *通信作者:翟麗梅E-maj1:zhaj1jmej@caas.cn

    作者簡介:崔超(1990—),男,在讀碩士研究生,主要從事流域氮磷養(yǎng)分管理及面源污染控制研究。E-maj1:cujchao5212009@163.com

    基金項目:公益性行業(yè)(農(nóng)業(yè))科研專項(201303089)

    收稿日期:2015-08-18

    中圖分類號:X53

    文獻標志碼:A

    文章編號:1672-2043(2016)01-0129-10doj:10.11654/jaes.2016.01.018

    猜你喜歡
    土地利用變化
    社區(qū)居民生計模式變遷與土地利用變化的耦合協(xié)調(diào)研究
    旅游學刊(2017年3期)2017-03-23 17:40:19
    鄭州市土地利用變化及生態(tài)環(huán)境效應(yīng)分析
    撫仙湖地區(qū)土地利用變化及驅(qū)動機制研究
    東北黑土區(qū)土地利用變化研究
    南京市土地利用變化與社會經(jīng)濟發(fā)展指標耦合研究
    商(2016年22期)2016-07-08 22:16:02
    基于最大似然法的南京市土地利用類型分類研究
    商(2016年22期)2016-07-08 14:32:30
    貴州花江示范區(qū)土地利用變化及驅(qū)動力分析
    基于土地利用變化的東川生態(tài)系統(tǒng)服務(wù)價值分析
    中國市場(2016年11期)2016-03-24 09:02:14
    基于地形因子的喀斯特山地土地利用變化分析
    昌吉市土地利用變化對生態(tài)環(huán)境影響分析
    18禁国产床啪视频网站| 午夜福利在线观看免费完整高清在 | 欧美绝顶高潮抽搐喷水| av女优亚洲男人天堂 | 国产亚洲av嫩草精品影院| 丰满人妻一区二区三区视频av | 国产精品久久久久久人妻精品电影| 99视频精品全部免费 在线 | 不卡一级毛片| 法律面前人人平等表现在哪些方面| 五月玫瑰六月丁香| 麻豆成人av在线观看| 久久久国产成人精品二区| 久久久国产精品麻豆| av中文乱码字幕在线| 美女被艹到高潮喷水动态| 18美女黄网站色大片免费观看| 亚洲片人在线观看| 国产成+人综合+亚洲专区| 久久精品国产99精品国产亚洲性色| 长腿黑丝高跟| 成年女人看的毛片在线观看| 亚洲av电影不卡..在线观看| 亚洲av日韩精品久久久久久密| 欧美日韩综合久久久久久 | 亚洲 欧美 日韩 在线 免费| 国产单亲对白刺激| 一个人免费在线观看电影 | 十八禁人妻一区二区| 国产私拍福利视频在线观看| 国产视频内射| 久久久久久久久久黄片| 99视频精品全部免费 在线 | 麻豆久久精品国产亚洲av| 天堂网av新在线| tocl精华| 国产黄a三级三级三级人| ponron亚洲| 国产精品电影一区二区三区| e午夜精品久久久久久久| 亚洲欧美一区二区三区黑人| 午夜免费观看网址| 久久欧美精品欧美久久欧美| 老司机在亚洲福利影院| 九色国产91popny在线| 亚洲色图 男人天堂 中文字幕| 欧美性猛交黑人性爽| 99re在线观看精品视频| 变态另类成人亚洲欧美熟女| 久久午夜综合久久蜜桃| 三级毛片av免费| 中亚洲国语对白在线视频| 亚洲国产欧洲综合997久久,| 日韩精品青青久久久久久| 亚洲国产精品合色在线| 国产成人av激情在线播放| 日本撒尿小便嘘嘘汇集6| 69av精品久久久久久| 国产精品免费一区二区三区在线| 午夜福利高清视频| 99国产精品一区二区三区| 成人18禁在线播放| 日韩有码中文字幕| 最近最新中文字幕大全免费视频| 黄色女人牲交| 国产免费男女视频| 欧美中文日本在线观看视频| 麻豆av在线久日| 久久久国产精品麻豆| 国产三级中文精品| 国产免费男女视频| 亚洲一区高清亚洲精品| 黄频高清免费视频| 国内揄拍国产精品人妻在线| 又黄又爽又免费观看的视频| 久久久国产成人精品二区| 天天添夜夜摸| 好男人电影高清在线观看| 看免费av毛片| 中文字幕精品亚洲无线码一区| 精品久久久久久久毛片微露脸| 12—13女人毛片做爰片一| 日韩欧美国产在线观看| 麻豆av在线久日| 国产av不卡久久| av福利片在线观看| 国产成人精品久久二区二区91| 精品国产乱子伦一区二区三区| 日韩欧美三级三区| 成人精品一区二区免费| 久久精品影院6| 麻豆一二三区av精品| 国产一区二区三区在线臀色熟女| 国产午夜精品论理片| 国产精品久久久久久人妻精品电影| 熟妇人妻久久中文字幕3abv| 亚洲国产欧美一区二区综合| 一区二区三区激情视频| 成人性生交大片免费视频hd| 日本 av在线| 亚洲中文字幕一区二区三区有码在线看 | 啦啦啦观看免费观看视频高清| 日本黄色片子视频| 给我免费播放毛片高清在线观看| 欧美性猛交╳xxx乱大交人| 日本 av在线| 日韩免费av在线播放| 老熟妇乱子伦视频在线观看| 成人欧美大片| 国产一区二区三区在线臀色熟女| 男女下面进入的视频免费午夜| 国产一级毛片七仙女欲春2| 三级男女做爰猛烈吃奶摸视频| 国产1区2区3区精品| 日韩欧美在线二视频| 美女 人体艺术 gogo| 青草久久国产| 日本a在线网址| 天堂√8在线中文| 97人妻精品一区二区三区麻豆| 国产一区二区三区视频了| 日韩精品中文字幕看吧| 欧美最黄视频在线播放免费| 亚洲欧洲精品一区二区精品久久久| 99精品久久久久人妻精品| 欧美在线一区亚洲| 757午夜福利合集在线观看| 嫩草影院精品99| 日韩精品中文字幕看吧| 国产成人av教育| 国内揄拍国产精品人妻在线| 波多野结衣巨乳人妻| 亚洲av第一区精品v没综合| 亚洲无线观看免费| 香蕉国产在线看| 日韩欧美三级三区| 老汉色∧v一级毛片| 91久久精品国产一区二区成人 | 日本与韩国留学比较| 观看免费一级毛片| 亚洲精华国产精华精| 岛国在线免费视频观看| 舔av片在线| 99在线人妻在线中文字幕| 国产真人三级小视频在线观看| 悠悠久久av| 在线观看免费视频日本深夜| 国产精品乱码一区二三区的特点| 精品国内亚洲2022精品成人| 丰满人妻熟妇乱又伦精品不卡| 国产v大片淫在线免费观看| 中文字幕精品亚洲无线码一区| 黄片大片在线免费观看| 国产成人精品久久二区二区91| 日韩三级视频一区二区三区| 亚洲激情在线av| 人人妻人人澡欧美一区二区| 国产精品综合久久久久久久免费| 色老头精品视频在线观看| 99国产精品99久久久久| 精品久久久久久久久久免费视频| 18禁裸乳无遮挡免费网站照片| 狠狠狠狠99中文字幕| 亚洲av电影在线进入| 美女高潮喷水抽搐中文字幕| 精品国产三级普通话版| 人人妻,人人澡人人爽秒播| 老汉色∧v一级毛片| 成年版毛片免费区| 老熟妇乱子伦视频在线观看| 搡老熟女国产l中国老女人| 91老司机精品| 成人av一区二区三区在线看| 亚洲va日本ⅴa欧美va伊人久久| 香蕉久久夜色| 国产亚洲精品综合一区在线观看| 怎么达到女性高潮| 两个人看的免费小视频| 在线播放国产精品三级| 国产成人啪精品午夜网站| 97超级碰碰碰精品色视频在线观看| 一本一本综合久久| 亚洲黑人精品在线| 国产高清有码在线观看视频| 久久天堂一区二区三区四区| 欧美精品啪啪一区二区三区| 精品一区二区三区av网在线观看| 国产成人精品久久二区二区91| 久久久久久九九精品二区国产| www.精华液| 国产精品av久久久久免费| 精品国产亚洲在线| 黑人欧美特级aaaaaa片| 国产精品久久视频播放| 午夜福利欧美成人| 韩国av一区二区三区四区| 国产淫片久久久久久久久 | 国产激情久久老熟女| 偷拍熟女少妇极品色| 亚洲av成人av| 丰满人妻一区二区三区视频av | 黑人欧美特级aaaaaa片| 老司机福利观看| 国产单亲对白刺激| 在线免费观看不下载黄p国产 | 国产精品一及| 99re在线观看精品视频| 日韩欧美精品v在线| 国产乱人视频| av在线蜜桃| 国产午夜福利久久久久久| 中文亚洲av片在线观看爽| 精品不卡国产一区二区三区| 女同久久另类99精品国产91| 在线a可以看的网站| 999精品在线视频| 老汉色∧v一级毛片| 黑人巨大精品欧美一区二区mp4| 欧美不卡视频在线免费观看| 老熟妇仑乱视频hdxx| 久久国产精品影院| 亚洲自偷自拍图片 自拍| 亚洲一区二区三区不卡视频| 中文字幕高清在线视频| 一区二区三区国产精品乱码| 国产在线精品亚洲第一网站| 一本精品99久久精品77| 久久热在线av| 99精品在免费线老司机午夜| 天堂av国产一区二区熟女人妻| 99久久精品一区二区三区| 国产伦精品一区二区三区四那| 亚洲 欧美一区二区三区| 99热精品在线国产| 99热这里只有是精品50| 亚洲人成网站高清观看| 一级毛片女人18水好多| 免费观看精品视频网站| 亚洲欧美日韩高清在线视频| 国产成人av激情在线播放| 中文字幕最新亚洲高清| 伦理电影免费视频| 两个人看的免费小视频| 亚洲在线观看片| 亚洲在线观看片| 亚洲专区字幕在线| 免费无遮挡裸体视频| 不卡一级毛片| 久久久久亚洲av毛片大全| 最新中文字幕久久久久 | 日韩欧美免费精品| 嫩草影院入口| 男人的好看免费观看在线视频| 午夜激情福利司机影院| 精品久久久久久成人av| 九色国产91popny在线| 91在线观看av| 网址你懂的国产日韩在线| 国产高潮美女av| 97人妻精品一区二区三区麻豆| 日韩欧美一区二区三区在线观看| 亚洲精品美女久久av网站| 欧美性猛交╳xxx乱大交人| 老鸭窝网址在线观看| 两性夫妻黄色片| 很黄的视频免费| 99热这里只有精品一区 | 免费在线观看视频国产中文字幕亚洲| 婷婷精品国产亚洲av| 99久久久亚洲精品蜜臀av| 色在线成人网| 亚洲激情在线av| 舔av片在线| 男插女下体视频免费在线播放| 成人国产一区最新在线观看| 亚洲成人精品中文字幕电影| 色在线成人网| 他把我摸到了高潮在线观看| 亚洲精品国产精品久久久不卡| 九色成人免费人妻av| 国产午夜精品久久久久久| a级毛片a级免费在线| 天天躁狠狠躁夜夜躁狠狠躁| 免费看a级黄色片| 久久精品91蜜桃| 夜夜看夜夜爽夜夜摸| 欧美av亚洲av综合av国产av| 18禁国产床啪视频网站| 国产亚洲欧美在线一区二区| 欧美日韩乱码在线| 18禁裸乳无遮挡免费网站照片| 久99久视频精品免费| 国产单亲对白刺激| 久久国产精品人妻蜜桃| 看黄色毛片网站| 国产黄片美女视频| 精品一区二区三区视频在线观看免费| 大型黄色视频在线免费观看| 男女下面进入的视频免费午夜| 一进一出好大好爽视频| 国产一区二区激情短视频| 久久精品影院6| 亚洲欧美日韩卡通动漫| 久久久国产成人免费| 狂野欧美激情性xxxx| 久久午夜综合久久蜜桃| 亚洲午夜理论影院| 亚洲专区字幕在线| 成年女人看的毛片在线观看| 日韩欧美国产一区二区入口| 无遮挡黄片免费观看| 哪里可以看免费的av片| 啪啪无遮挡十八禁网站| 狂野欧美白嫩少妇大欣赏| 国产精品久久电影中文字幕| 中文字幕av在线有码专区| 日韩中文字幕欧美一区二区| 91av网一区二区| 国产亚洲欧美在线一区二区| 亚洲av第一区精品v没综合| 在线观看免费午夜福利视频| 国产精品 欧美亚洲| 亚洲无线观看免费| 美女扒开内裤让男人捅视频| 精品一区二区三区av网在线观看| 在线十欧美十亚洲十日本专区| 女同久久另类99精品国产91| 免费看a级黄色片| 日本精品一区二区三区蜜桃| 国产一区二区三区在线臀色熟女| 久久天躁狠狠躁夜夜2o2o| 成年免费大片在线观看| 国产探花在线观看一区二区| 黄色日韩在线| 亚洲精品一区av在线观看| 男女床上黄色一级片免费看| 日日干狠狠操夜夜爽| 高清毛片免费观看视频网站| 婷婷丁香在线五月| netflix在线观看网站| 久久久久久九九精品二区国产| 色综合亚洲欧美另类图片| 一个人免费在线观看电影 | 两性午夜刺激爽爽歪歪视频在线观看| 亚洲人成网站在线播放欧美日韩| 搡老岳熟女国产| 五月玫瑰六月丁香| 色播亚洲综合网| 国产97色在线日韩免费| 日本一本二区三区精品| 又黄又爽又免费观看的视频| 好男人电影高清在线观看| 久久精品影院6| 午夜福利视频1000在线观看| 欧美中文日本在线观看视频| 精品久久久久久久人妻蜜臀av| 变态另类成人亚洲欧美熟女| 搡老岳熟女国产| 琪琪午夜伦伦电影理论片6080| 国产亚洲精品久久久com| 蜜桃久久精品国产亚洲av| 国产午夜精品论理片| 男人舔女人的私密视频| 色综合亚洲欧美另类图片| 观看美女的网站| 欧美日韩福利视频一区二区| 成年免费大片在线观看| 精品国内亚洲2022精品成人| 久久国产精品影院| 国产成人系列免费观看| 制服丝袜大香蕉在线| 一进一出抽搐动态| 亚洲欧洲精品一区二区精品久久久| 国产激情偷乱视频一区二区| 麻豆一二三区av精品| 国内久久婷婷六月综合欲色啪| 成人无遮挡网站| 国产一区在线观看成人免费| 99视频精品全部免费 在线 | 亚洲精品粉嫩美女一区| 亚洲精品粉嫩美女一区| 麻豆av在线久日| 国产99白浆流出| 村上凉子中文字幕在线| 黄色视频,在线免费观看| 美女 人体艺术 gogo| 亚洲av中文字字幕乱码综合| 日韩人妻高清精品专区| 免费高清视频大片| 好男人在线观看高清免费视频| 真实男女啪啪啪动态图| 亚洲电影在线观看av| 国产蜜桃级精品一区二区三区| 母亲3免费完整高清在线观看| 国产麻豆成人av免费视频| 美女高潮的动态| 亚洲中文日韩欧美视频| 国产成+人综合+亚洲专区| 91字幕亚洲| 国产高清视频在线观看网站| av视频在线观看入口| 老司机午夜十八禁免费视频| 少妇的逼水好多| 免费看a级黄色片| 极品教师在线免费播放| 午夜激情欧美在线| 99久久99久久久精品蜜桃| 99久久久亚洲精品蜜臀av| 日韩三级视频一区二区三区| 麻豆成人午夜福利视频| 亚洲天堂国产精品一区在线| 亚洲成人免费电影在线观看| 国产精品免费一区二区三区在线| 麻豆av在线久日| 亚洲欧洲精品一区二区精品久久久| 亚洲最大成人中文| 最新中文字幕久久久久 | 国产视频内射| 黑人操中国人逼视频| 亚洲美女黄片视频| 亚洲九九香蕉| 精品午夜福利视频在线观看一区| 亚洲色图av天堂| 九九久久精品国产亚洲av麻豆 | 在线看三级毛片| 好看av亚洲va欧美ⅴa在| 欧美av亚洲av综合av国产av| 性色av乱码一区二区三区2| 狂野欧美白嫩少妇大欣赏| 国产伦在线观看视频一区| 丁香欧美五月| 欧美激情在线99| 99久久成人亚洲精品观看| 母亲3免费完整高清在线观看| 黄色视频,在线免费观看| 老司机午夜十八禁免费视频| 99国产精品99久久久久| 亚洲av成人不卡在线观看播放网| 淫秽高清视频在线观看| 国产日本99.免费观看| 99久久99久久久精品蜜桃| 午夜福利在线在线| 俄罗斯特黄特色一大片| 亚洲国产中文字幕在线视频| 成人永久免费在线观看视频| 久久久国产欧美日韩av| 91在线精品国自产拍蜜月 | 九九热线精品视视频播放| 热99re8久久精品国产| 亚洲中文日韩欧美视频| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 精品午夜福利视频在线观看一区| 亚洲色图av天堂| 日韩欧美一区二区三区在线观看| 久久久久九九精品影院| 丁香六月欧美| av中文乱码字幕在线| 中文字幕熟女人妻在线| 中出人妻视频一区二区| 黑人巨大精品欧美一区二区mp4| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 成人永久免费在线观看视频| 久久亚洲精品不卡| 久久久久九九精品影院| 精品国内亚洲2022精品成人| 欧美一区二区精品小视频在线| 一个人免费在线观看的高清视频| 久久草成人影院| 女同久久另类99精品国产91| 国产av不卡久久| 美女 人体艺术 gogo| 757午夜福利合集在线观看| 午夜免费观看网址| 国产麻豆成人av免费视频| 好看av亚洲va欧美ⅴa在| 中文字幕人成人乱码亚洲影| а√天堂www在线а√下载| av黄色大香蕉| 色哟哟哟哟哟哟| 噜噜噜噜噜久久久久久91| 亚洲男人的天堂狠狠| 日本黄大片高清| 99国产综合亚洲精品| 日韩欧美国产在线观看| 国内少妇人妻偷人精品xxx网站 | 亚洲 欧美一区二区三区| 亚洲美女黄片视频| 在线观看午夜福利视频| 中文字幕最新亚洲高清| 身体一侧抽搐| 成在线人永久免费视频| 成人亚洲精品av一区二区| 在线a可以看的网站| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器 | 国产三级中文精品| 免费观看人在逋| 日本三级黄在线观看| 国产一区二区三区视频了| 亚洲av日韩精品久久久久久密| 亚洲第一欧美日韩一区二区三区| 国产精品一区二区精品视频观看| 18禁国产床啪视频网站| 亚洲欧美激情综合另类| 国内精品一区二区在线观看| 亚洲,欧美精品.| 三级毛片av免费| 色老头精品视频在线观看| 欧美黑人欧美精品刺激| 欧美色欧美亚洲另类二区| 欧美日韩瑟瑟在线播放| 国产精品一区二区三区四区久久| 久久这里只有精品19| 国产精品1区2区在线观看.| www.999成人在线观看| 91麻豆av在线| 国产一级毛片七仙女欲春2| 国产精品自产拍在线观看55亚洲| 国产激情久久老熟女| 麻豆久久精品国产亚洲av| 观看美女的网站| 精品国产乱子伦一区二区三区| 人人妻人人澡欧美一区二区| 中国美女看黄片| 亚洲欧美日韩东京热| 久久精品影院6| 国产精品 欧美亚洲| 亚洲精品一区av在线观看| 美女cb高潮喷水在线观看 | 久久精品国产清高在天天线| 很黄的视频免费| 日本黄大片高清| 亚洲国产精品成人综合色| 国产亚洲精品综合一区在线观看| 久久久成人免费电影| 久久精品国产99精品国产亚洲性色| 精品不卡国产一区二区三区| 免费看美女性在线毛片视频| 亚洲人成网站高清观看| 亚洲国产精品999在线| 久久久国产成人免费| 亚洲专区国产一区二区| 国产精品av久久久久免费| 免费在线观看日本一区| 草草在线视频免费看| 成人鲁丝片一二三区免费| 亚洲成人久久性| 这个男人来自地球电影免费观看| 91九色精品人成在线观看| 免费高清视频大片| 丰满人妻熟妇乱又伦精品不卡| 亚洲自拍偷在线| 可以在线观看的亚洲视频| 精品不卡国产一区二区三区| 性色av乱码一区二区三区2| 国产欧美日韩精品一区二区| 亚洲激情在线av| 99热这里只有精品一区 | 99热只有精品国产| 99热这里只有精品一区 | 国产黄a三级三级三级人| 国产aⅴ精品一区二区三区波| av国产免费在线观看| 成年免费大片在线观看| 国产精品影院久久| 亚洲欧美日韩高清专用| 老熟妇乱子伦视频在线观看| 久久久久久大精品| 亚洲 欧美一区二区三区| 黄片小视频在线播放| 国产av一区在线观看免费| 搡老岳熟女国产| 夜夜夜夜夜久久久久| 国内少妇人妻偷人精品xxx网站 | 亚洲av成人一区二区三| 变态另类丝袜制服| 中文在线观看免费www的网站| 国产一区在线观看成人免费| 精品无人区乱码1区二区| 老熟妇仑乱视频hdxx| 日日干狠狠操夜夜爽| 国产精品99久久久久久久久| 久久久水蜜桃国产精品网| 国产人伦9x9x在线观看| 最新中文字幕久久久久 | 亚洲 欧美一区二区三区| 亚洲av免费在线观看| 国产在线精品亚洲第一网站| avwww免费| 欧美一级毛片孕妇| 村上凉子中文字幕在线| 99久久国产精品久久久| 国产高清视频在线播放一区| 亚洲成av人片免费观看| 国产一级毛片七仙女欲春2| a在线观看视频网站| 国内精品一区二区在线观看| 免费观看人在逋| 人人妻人人看人人澡| 女生性感内裤真人,穿戴方法视频| 国产在线精品亚洲第一网站| 黄片小视频在线播放| 久久久久久久久免费视频了| 丰满的人妻完整版| 国产精品永久免费网站| 国产主播在线观看一区二区| 夜夜夜夜夜久久久久| 久久国产乱子伦精品免费另类| 国产日本99.免费观看| 久久欧美精品欧美久久欧美| 亚洲无线在线观看| 亚洲精品中文字幕一二三四区| 在线观看午夜福利视频|