• <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)境影響分析
    99热精品在线国产| 黄色女人牲交| 乱人视频在线观看| 亚洲成人久久爱视频| 99热这里只有精品一区| xxx96com| 精品不卡国产一区二区三区| www.色视频.com| 99国产综合亚洲精品| 精品一区二区三区视频在线 | ponron亚洲| 色综合亚洲欧美另类图片| 18禁美女被吸乳视频| 可以在线观看的亚洲视频| 两个人看的免费小视频| 精品国产超薄肉色丝袜足j| 精品国产三级普通话版| 成人午夜高清在线视频| 特大巨黑吊av在线直播| 亚洲最大成人中文| 少妇人妻一区二区三区视频| 日本a在线网址| 亚洲成a人片在线一区二区| 国产精品野战在线观看| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久午夜电影| av女优亚洲男人天堂| 国产在视频线在精品| 欧美激情久久久久久爽电影| avwww免费| 国产黄色小视频在线观看| 午夜激情欧美在线| 久久精品国产综合久久久| 我的老师免费观看完整版| 日韩欧美一区二区三区在线观看| 日本 欧美在线| 特大巨黑吊av在线直播| 老司机在亚洲福利影院| 亚洲国产精品sss在线观看| 九九热线精品视视频播放| 亚洲熟妇熟女久久| 国产亚洲精品久久久com| 亚洲人成网站在线播| 中文字幕av在线有码专区| 91字幕亚洲| 床上黄色一级片| 给我免费播放毛片高清在线观看| 精品久久久久久久久久久久久| 欧美不卡视频在线免费观看| 九九久久精品国产亚洲av麻豆| 岛国在线免费视频观看| 国产精品 欧美亚洲| 久久久久国产精品人妻aⅴ院| 少妇人妻一区二区三区视频| 精品久久久久久久末码| 日本在线视频免费播放| 三级男女做爰猛烈吃奶摸视频| av欧美777| 69av精品久久久久久| 99热这里只有精品一区| 亚洲国产高清在线一区二区三| 免费av毛片视频| 男女视频在线观看网站免费| 国内久久婷婷六月综合欲色啪| 网址你懂的国产日韩在线| 午夜福利欧美成人| 成人无遮挡网站| 久久伊人香网站| 白带黄色成豆腐渣| 国产精品,欧美在线| 在线观看午夜福利视频| 日韩大尺度精品在线看网址| 老汉色∧v一级毛片| 男女床上黄色一级片免费看| bbb黄色大片| 亚洲av中文字字幕乱码综合| 日本与韩国留学比较| 国产精品自产拍在线观看55亚洲| 熟女少妇亚洲综合色aaa.| 亚洲最大成人中文| 97超视频在线观看视频| 日韩有码中文字幕| 波多野结衣巨乳人妻| 嫩草影视91久久| 久久久国产精品麻豆| 国产日本99.免费观看| 99久久99久久久精品蜜桃| 国产精品影院久久| 舔av片在线| 亚洲成a人片在线一区二区| 国产免费男女视频| 国产蜜桃级精品一区二区三区| 久久亚洲真实| 丰满人妻熟妇乱又伦精品不卡| 一级黄色大片毛片| 亚洲熟妇熟女久久| 亚洲第一欧美日韩一区二区三区| 午夜两性在线视频| 99精品欧美一区二区三区四区| 麻豆成人午夜福利视频| 麻豆一二三区av精品| 久久精品人妻少妇| а√天堂www在线а√下载| 亚洲欧美精品综合久久99| 色老头精品视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久亚洲真实| 天堂网av新在线| 不卡一级毛片| 日韩欧美免费精品| 精品99又大又爽又粗少妇毛片 | 一个人免费在线观看电影| 国产亚洲精品久久久com| 两个人看的免费小视频| 不卡一级毛片| 噜噜噜噜噜久久久久久91| 国产精品美女特级片免费视频播放器| 亚洲国产日韩欧美精品在线观看 | 午夜激情欧美在线| 国产免费男女视频| 成人三级黄色视频| 18美女黄网站色大片免费观看| 亚洲国产中文字幕在线视频| 757午夜福利合集在线观看| av国产免费在线观看| 女人十人毛片免费观看3o分钟| 成年女人毛片免费观看观看9| 亚洲一区二区三区不卡视频| 欧美bdsm另类| 性色avwww在线观看| 两个人看的免费小视频| 大型黄色视频在线免费观看| 亚洲av一区综合| 97超级碰碰碰精品色视频在线观看| 韩国av一区二区三区四区| 婷婷精品国产亚洲av在线| 露出奶头的视频| 51国产日韩欧美| 日韩欧美在线乱码| 3wmmmm亚洲av在线观看| 在线观看午夜福利视频| 成人欧美大片| 人人妻,人人澡人人爽秒播| 91在线精品国自产拍蜜月 | 成年女人永久免费观看视频| 特大巨黑吊av在线直播| 国产高潮美女av| 亚洲狠狠婷婷综合久久图片| 亚洲成a人片在线一区二区| 精品一区二区三区av网在线观看| 午夜免费成人在线视频| 欧美日韩中文字幕国产精品一区二区三区| 在线观看日韩欧美| 亚洲国产日韩欧美精品在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 久久久久久九九精品二区国产| 国产精品自产拍在线观看55亚洲| 九色国产91popny在线| 一区二区三区高清视频在线| 黄色视频,在线免费观看| 老司机福利观看| 欧美区成人在线视频| 国产三级黄色录像| 欧美一级毛片孕妇| 每晚都被弄得嗷嗷叫到高潮| 母亲3免费完整高清在线观看| 国产高清三级在线| 国产午夜精品久久久久久一区二区三区 | www.www免费av| 看免费av毛片| 国产精品一及| 波野结衣二区三区在线 | 色av中文字幕| 91久久精品电影网| 一个人观看的视频www高清免费观看| 国产成人av激情在线播放| 人人妻人人澡欧美一区二区| 欧美成人一区二区免费高清观看| 欧美高清成人免费视频www| 国内揄拍国产精品人妻在线| 久久午夜亚洲精品久久| 成人高潮视频无遮挡免费网站| 免费av毛片视频| 51国产日韩欧美| 白带黄色成豆腐渣| 舔av片在线| 黄色丝袜av网址大全| 91字幕亚洲| 91久久精品国产一区二区成人 | 午夜福利欧美成人| 少妇人妻精品综合一区二区 | 国产精品综合久久久久久久免费| 99久久无色码亚洲精品果冻| 看片在线看免费视频| 91av网一区二区| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 丁香六月欧美| 日本 av在线| 18禁黄网站禁片午夜丰满| 欧美日韩精品网址| 国产一区二区三区在线臀色熟女| 黑人欧美特级aaaaaa片| 一级a爱片免费观看的视频| 午夜免费激情av| 午夜免费男女啪啪视频观看 | 色综合欧美亚洲国产小说| 亚洲 欧美 日韩 在线 免费| 91在线精品国自产拍蜜月 | av在线天堂中文字幕| 欧美国产日韩亚洲一区| 午夜福利18| 久久九九热精品免费| 一本精品99久久精品77| 精品人妻1区二区| 露出奶头的视频| 亚洲avbb在线观看| 亚洲精品在线观看二区| 亚洲国产色片| 在线a可以看的网站| 欧美日韩中文字幕国产精品一区二区三区| 国产色爽女视频免费观看| 亚洲精品亚洲一区二区| 亚洲精品影视一区二区三区av| 成人欧美大片| 欧美黄色片欧美黄色片| 久久伊人香网站| 国产乱人视频| 免费看十八禁软件| 欧美乱码精品一区二区三区| 免费av毛片视频| 老司机午夜十八禁免费视频| 12—13女人毛片做爰片一| 日韩av在线大香蕉| 老司机深夜福利视频在线观看| 男女床上黄色一级片免费看| 精品久久久久久久久久久久久| 男女午夜视频在线观看| 国产欧美日韩一区二区精品| 熟女电影av网| 亚洲av熟女| 亚洲欧美日韩高清在线视频| 黄色女人牲交| 午夜福利成人在线免费观看| 女人高潮潮喷娇喘18禁视频| 亚洲国产色片| 亚洲精品在线美女| 丰满人妻一区二区三区视频av | 亚洲一区二区三区色噜噜| 精品免费久久久久久久清纯| 人人妻人人看人人澡| 色吧在线观看| 亚洲男人的天堂狠狠| 午夜福利在线在线| 午夜福利视频1000在线观看| 国产精品99久久99久久久不卡| 久久久精品欧美日韩精品| 日本一本二区三区精品| 欧美午夜高清在线| 蜜桃久久精品国产亚洲av| av天堂在线播放| 一本久久中文字幕| 免费无遮挡裸体视频| 亚洲精品国产精品久久久不卡| 级片在线观看| 久久99热这里只有精品18| aaaaa片日本免费| 熟女电影av网| 成人无遮挡网站| 特级一级黄色大片| 欧美性猛交╳xxx乱大交人| 黑人欧美特级aaaaaa片| 国产欧美日韩一区二区三| 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 性色avwww在线观看| 少妇熟女aⅴ在线视频| 一个人免费在线观看的高清视频| 成人一区二区视频在线观看| 日韩欧美三级三区| 在线播放国产精品三级| 婷婷六月久久综合丁香| 人妻久久中文字幕网| 波多野结衣高清无吗| 男插女下体视频免费在线播放| 欧美激情久久久久久爽电影| 亚洲av第一区精品v没综合| 人人妻人人看人人澡| 五月玫瑰六月丁香| 中文在线观看免费www的网站| 国产一区在线观看成人免费| 亚洲第一欧美日韩一区二区三区| 日韩成人在线观看一区二区三区| xxxwww97欧美| 又粗又爽又猛毛片免费看| 成人无遮挡网站| 国产精品一区二区三区四区免费观看 | 最新在线观看一区二区三区| 成人国产一区最新在线观看| 午夜精品久久久久久毛片777| 又紧又爽又黄一区二区| 性色av乱码一区二区三区2| 99久久精品一区二区三区| 亚洲av美国av| 欧美大码av| 国产不卡一卡二| 51国产日韩欧美| 我的老师免费观看完整版| 乱人视频在线观看| 中文在线观看免费www的网站| 搡老岳熟女国产| 欧美成狂野欧美在线观看| 国产精品野战在线观看| 99精品欧美一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 3wmmmm亚洲av在线观看| 亚洲av电影不卡..在线观看| 一级a爱片免费观看的视频| aaaaa片日本免费| 国内少妇人妻偷人精品xxx网站| www.www免费av| 一个人看的www免费观看视频| 少妇裸体淫交视频免费看高清| 国产伦精品一区二区三区四那| 国产高清视频在线观看网站| 亚洲人成网站高清观看| 中文字幕人妻丝袜一区二区| 中文在线观看免费www的网站| 日本 av在线| 色综合欧美亚洲国产小说| 动漫黄色视频在线观看| 一本久久中文字幕| 久久香蕉精品热| 一区二区三区免费毛片| 国产老妇女一区| 国产精品 国内视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品国产精品久久久不卡| 淫妇啪啪啪对白视频| 宅男免费午夜| 最近最新中文字幕大全电影3| 色综合亚洲欧美另类图片| 亚洲不卡免费看| 午夜a级毛片| 成年女人毛片免费观看观看9| 村上凉子中文字幕在线| 免费看日本二区| 人人妻人人澡欧美一区二区| 久久久国产成人精品二区| 国产精品影院久久| 日本黄色片子视频| 搡老岳熟女国产| www国产在线视频色| 婷婷精品国产亚洲av在线| 午夜精品久久久久久毛片777| 精品一区二区三区人妻视频| 亚洲人成电影免费在线| 国产三级中文精品| 在线免费观看不下载黄p国产 | 亚洲最大成人中文| 免费人成视频x8x8入口观看| 又紧又爽又黄一区二区| 久久性视频一级片| 国产一区二区三区视频了| 桃红色精品国产亚洲av| 国产精品亚洲一级av第二区| 亚洲成av人片免费观看| 一进一出抽搐gif免费好疼| 成年人黄色毛片网站| 不卡一级毛片| 国产美女午夜福利| 亚洲av成人不卡在线观看播放网| 日本黄色片子视频| 国产一区二区在线观看日韩 | 亚洲美女黄片视频| 欧美成人免费av一区二区三区| 日韩成人在线观看一区二区三区| 人妻久久中文字幕网| 一a级毛片在线观看| 婷婷丁香在线五月| av欧美777| 久久伊人香网站| www.色视频.com| svipshipincom国产片| 国产主播在线观看一区二区| 国产aⅴ精品一区二区三区波| 丰满人妻一区二区三区视频av | 国产精品av视频在线免费观看| 国产精品 欧美亚洲| 国产免费av片在线观看野外av| 亚洲不卡免费看| 国产免费av片在线观看野外av| 天堂影院成人在线观看| 亚洲av成人av| 亚洲国产欧美网| 俄罗斯特黄特色一大片| 色老头精品视频在线观看| 久久亚洲真实| 精品国产亚洲在线| 可以在线观看的亚洲视频| 国产三级中文精品| www国产在线视频色| 国产一区二区在线观看日韩 | 精品久久久久久,| 首页视频小说图片口味搜索| av国产免费在线观看| 精品99又大又爽又粗少妇毛片 | 美女免费视频网站| 丰满人妻熟妇乱又伦精品不卡| 18禁国产床啪视频网站| 人妻丰满熟妇av一区二区三区| 有码 亚洲区| 99久久99久久久精品蜜桃| 偷拍熟女少妇极品色| 免费看光身美女| 免费搜索国产男女视频| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 大型黄色视频在线免费观看| 国产成人欧美在线观看| 久久婷婷人人爽人人干人人爱| 制服人妻中文乱码| 国产三级黄色录像| 特级一级黄色大片| 久久天躁狠狠躁夜夜2o2o| av天堂中文字幕网| 最近在线观看免费完整版| 精品一区二区三区人妻视频| 老司机午夜十八禁免费视频| 日韩 欧美 亚洲 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 岛国视频午夜一区免费看| 国产色婷婷99| 综合色av麻豆| 欧美绝顶高潮抽搐喷水| 欧美xxxx黑人xx丫x性爽| 熟女电影av网| 亚洲人与动物交配视频| 91av网一区二区| 色尼玛亚洲综合影院| 久久久久久久久大av| 国产精品一区二区免费欧美| 国产三级黄色录像| 两性午夜刺激爽爽歪歪视频在线观看| 国模一区二区三区四区视频| 国产高潮美女av| 日本撒尿小便嘘嘘汇集6| www日本黄色视频网| 亚洲不卡免费看| 深夜精品福利| 美女大奶头视频| 欧美极品一区二区三区四区| 他把我摸到了高潮在线观看| 美女被艹到高潮喷水动态| 国产一区二区三区在线臀色熟女| 国产毛片a区久久久久| 午夜福利在线观看吧| 动漫黄色视频在线观看| x7x7x7水蜜桃| www日本黄色视频网| 亚洲天堂国产精品一区在线| 欧美日韩中文字幕国产精品一区二区三区| 国产精品乱码一区二三区的特点| 日韩欧美精品免费久久 | 免费av毛片视频| 日韩大尺度精品在线看网址| 免费大片18禁| 黄色日韩在线| 最近在线观看免费完整版| 免费看a级黄色片| 亚洲男人的天堂狠狠| 亚洲精品色激情综合| 色综合婷婷激情| 亚洲欧美日韩卡通动漫| 亚洲中文字幕一区二区三区有码在线看| 少妇人妻精品综合一区二区 | 国产国拍精品亚洲av在线观看 | 久久久久久久午夜电影| 日韩精品中文字幕看吧| 国产精华一区二区三区| 婷婷精品国产亚洲av在线| 亚洲久久久久久中文字幕| 99精品欧美一区二区三区四区| 免费av观看视频| 免费搜索国产男女视频| 成人av一区二区三区在线看| 首页视频小说图片口味搜索| 宅男免费午夜| 91在线观看av| 久久久久国内视频| 波多野结衣高清无吗| 国产精品久久电影中文字幕| 中文字幕人妻丝袜一区二区| 欧美乱色亚洲激情| 一区二区三区免费毛片| 亚洲 国产 在线| 成人亚洲精品av一区二区| 欧美黑人欧美精品刺激| 韩国av一区二区三区四区| 亚洲精品影视一区二区三区av| 日韩欧美 国产精品| 蜜桃亚洲精品一区二区三区| 日本精品一区二区三区蜜桃| 91九色精品人成在线观看| 欧美xxxx黑人xx丫x性爽| 大型黄色视频在线免费观看| 国产精品久久电影中文字幕| 国产精品三级大全| 少妇熟女aⅴ在线视频| 亚洲最大成人手机在线| 色在线成人网| 国产高清三级在线| 我要搜黄色片| 在线天堂最新版资源| 国产日本99.免费观看| 国产蜜桃级精品一区二区三区| a级一级毛片免费在线观看| 性色av乱码一区二区三区2| 久久精品91蜜桃| 国产亚洲欧美98| 18禁国产床啪视频网站| 亚洲无线观看免费| 男人的好看免费观看在线视频| 97碰自拍视频| 在线播放国产精品三级| 精品国内亚洲2022精品成人| 久久精品亚洲精品国产色婷小说| 中国美女看黄片| 大型黄色视频在线免费观看| 精品久久久久久成人av| 国产一区二区激情短视频| 色老头精品视频在线观看| 99在线视频只有这里精品首页| 久久久久亚洲av毛片大全| 99久久成人亚洲精品观看| 一个人免费在线观看的高清视频| 精品不卡国产一区二区三区| 两人在一起打扑克的视频| 99久久精品热视频| 亚洲成av人片在线播放无| www.www免费av| 国产野战对白在线观看| 精品一区二区三区av网在线观看| 亚洲七黄色美女视频| 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久成人av| 午夜福利免费观看在线| 欧美区成人在线视频| 夜夜夜夜夜久久久久| 精品国产美女av久久久久小说| 中文字幕av在线有码专区| 欧美成人一区二区免费高清观看| 最好的美女福利视频网| 国产又黄又爽又无遮挡在线| 观看美女的网站| 三级国产精品欧美在线观看| 国产v大片淫在线免费观看| 中文字幕人妻丝袜一区二区| 又黄又爽又免费观看的视频| 一进一出抽搐动态| 精品国产美女av久久久久小说| 好看av亚洲va欧美ⅴa在| 亚洲成av人片免费观看| 欧美bdsm另类| 成人精品一区二区免费| a在线观看视频网站| 少妇人妻一区二区三区视频| 国产激情欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 精品人妻1区二区| 99热这里只有是精品50| 男女之事视频高清在线观看| 国产高清视频在线播放一区| 国内久久婷婷六月综合欲色啪| 国产精品精品国产色婷婷| 午夜福利在线观看免费完整高清在 | 国产欧美日韩一区二区精品| 精品国产三级普通话版| 亚洲成av人片在线播放无| 国产亚洲精品久久久com| 日韩 欧美 亚洲 中文字幕| 少妇的逼水好多| 一本久久中文字幕| 五月伊人婷婷丁香| 一区二区三区国产精品乱码| 黄片大片在线免费观看| av国产免费在线观看| 一级毛片高清免费大全| 亚洲五月婷婷丁香| 欧美bdsm另类| 99久久99久久久精品蜜桃| 欧美乱妇无乱码| 亚洲欧美日韩无卡精品| 国产精品野战在线观看| 午夜免费激情av| 国产一区二区在线观看日韩 | 亚洲精品美女久久久久99蜜臀| 欧美日韩精品网址| 桃色一区二区三区在线观看| 日本精品一区二区三区蜜桃| h日本视频在线播放| 美女高潮的动态| 亚洲精品国产精品久久久不卡| 国产老妇女一区| 婷婷丁香在线五月| 国产精品爽爽va在线观看网站| 三级男女做爰猛烈吃奶摸视频| 少妇的丰满在线观看| 变态另类丝袜制服| 国产视频一区二区在线看| 亚洲av不卡在线观看| 嫩草影院精品99| 悠悠久久av| 色尼玛亚洲综合影院| 精品免费久久久久久久清纯|