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

    土地利用結(jié)構(gòu)與景觀格局對開都河水質(zhì)的影響

    2025-01-28 00:00:00王凱迪麗努爾·阿吉李平平胡叢巧康佳佳彭清清
    水生態(tài)學(xué)雜志 2025年1期
    關(guān)鍵詞:開都河景觀格局土地利用

    摘要:探究土地利用結(jié)構(gòu)與景觀格局對干旱區(qū)河流水質(zhì)的影響機(jī)制,對維護(hù)干旱區(qū)流域生態(tài)安全具有參考意義。基于開都河2021年豐水期和枯水期15個采樣點(diǎn)的水質(zhì)數(shù)據(jù),以流域景觀類型占比表征土地利用結(jié)構(gòu),景觀指數(shù)表征景觀格局,采用bioenv分析、冗余分析(RDA)等方法探討土地利用結(jié)構(gòu)與景觀格局對開都河水質(zhì)的影響。結(jié)果表明:(1)開都河屬清潔、較清潔水質(zhì),枯水期水質(zhì)略優(yōu)于豐水期;(2)豐、枯水期的2 000、5 000 m半徑緩沖區(qū)土地利用對水質(zhì)解釋率最佳,分別為63.81%、64.91%;景觀格局對水質(zhì)解釋率的最佳緩沖區(qū)分別為500、200 m,解釋率分別為73.88%、63.33%;(3)豐水期耕地對水質(zhì)的貢獻(xiàn)率最高為57.9%,景觀指數(shù)LPI對水質(zhì)的貢獻(xiàn)率最高為48.8%;枯水期對水質(zhì)貢獻(xiàn)率最高的為草地68%,景觀指數(shù)LSI對水質(zhì)貢獻(xiàn)率最高為21.6%;(4)豐水期,耕地、林地、水域的組合與水質(zhì)的相關(guān)系數(shù)最高為0.424 5,枯水期單一草地的組合與水質(zhì)的相關(guān)系數(shù)最高為0.512 9;豐水期NP、LPI、ENN_MN與水質(zhì)的相關(guān)系數(shù)最高為0.076 5,枯水期NP、PD、LPI與水質(zhì)的相關(guān)系數(shù)最高為0.216 3;(5)土地利用結(jié)構(gòu)與景觀格局對水質(zhì)的交互作用在豐水期2 000 m半徑緩沖區(qū)內(nèi)為18%,高于枯水期200 m半徑緩沖區(qū)的11%。研究結(jié)果表明土地利用結(jié)構(gòu)與景觀格局都對開都河水質(zhì)有較大影響,二者的交互作用在該影響中占有重要地位。

    關(guān)鍵詞:土地利用;景觀格局;水質(zhì);開都河

    中圖分類號:X522" " " " 文獻(xiàn)標(biāo)志碼:A" " " nbsp; 文章編號:1674-3075(2025)01-0045-13

    河流水質(zhì)是公共衛(wèi)生和社會經(jīng)濟(jì)可持續(xù)發(fā)展的基礎(chǔ)(Castrillo amp; García,2020),而河流水環(huán)境受周圍土地利用與景觀格局的影響顯著(Allan et al,2004;Sarah amp; Hobbi et al,2017)。人類活動使土地利用方式和自然景觀格局發(fā)生了改變。Tu(2011)使用地理加權(quán)回歸(GWR)分析不同城市流域中6種土地利用與14種水質(zhì)指標(biāo)之間的空間變化關(guān)系,發(fā)現(xiàn)土地利用結(jié)構(gòu)和景觀格局是影響河流水質(zhì)的主要原因。土地利用對水質(zhì)的影響主要是非點(diǎn)源污染,景觀格局受非點(diǎn)源污染的影響同樣十分顯著(Basnyat et al,1999)。目前,關(guān)于土地利用與景觀格局對水質(zhì)的影響研究越來越多,但是對于哪種景觀模式和土地利用情況更有利于水質(zhì)保護(hù)有很大不確定性。因此,深入研究土地利用與景觀格局對水質(zhì)的影響,不僅為土地利用和景觀格局規(guī)劃提供參考依據(jù),也有利于地表水環(huán)境的可持續(xù)發(fā)展(崔丹等,2019)。

    土地利用和景觀格局對水環(huán)境的影響一般建立在一定空間尺度上(祁蘭蘭等,2021),空間尺度是指河流和景觀斑塊之間的距離,常以河岸帶和子流域表現(xiàn)。而從景觀格局角度研究土地利用對水質(zhì)影響,主要分析景觀格局指數(shù)與水質(zhì)指標(biāo)的關(guān)系。郭羽羽等(2021)在黃河流域的研究中發(fā)現(xiàn)耕地、草地、林地與水質(zhì)的相關(guān)性最高。李明濤等(2013)、Bu等(2014)分別對密云水庫和太子江流域水質(zhì)展開研究,結(jié)果表明景觀格局指標(biāo)中,聚集度、破碎度、多樣性特征以及物理連接度與水質(zhì)具有顯著相關(guān)性。此外由于空間尺度效應(yīng)存在地域差異性,土地利用與景觀格局對水質(zhì)影響最強(qiáng)的空間尺度存在差異(張宇碩等,2020)。胡叢巧等(2023)對博斯騰湖土地利用對水質(zhì)的影響研究發(fā)現(xiàn)豐、平、枯水期的1 000、3 000、500 m半徑緩沖區(qū)對水質(zhì)解釋率最佳,而王小平等(2017)的研究則得出相反結(jié)論,其認(rèn)為4 000 m緩沖區(qū)尺度下兩者關(guān)系較全流域尺度強(qiáng)。Meneses等(2015)認(rèn)為,相較于流域尺度,景觀格局對水環(huán)境質(zhì)量的依賴性在河岸帶緩沖區(qū)尺度更小,且Xu等(2019)在研究貴州烏江流域后發(fā)現(xiàn)影響河流水質(zhì)的關(guān)鍵區(qū)域?yàn)楹影稁Ь彌_區(qū)尺度內(nèi)的景觀格局,而曹燦等(2018)對干旱區(qū)艾比湖的研究,發(fā)現(xiàn)景觀格局對水質(zhì)影響的最佳緩沖區(qū)為300 m。這些研究雖然取得了一定的成果,但對于干旱區(qū)土地利用與景觀格局對水質(zhì)的影響研究較少。同時,前人研究主要聚焦于單一土地利用類型或景觀格局指數(shù)對水質(zhì)的影響,而對于景觀指數(shù)組合與土地利用類型組合對水質(zhì)影響的研究較少。

    因此,本文利用2021年開都河豐水期和枯水期15個樣點(diǎn)的水質(zhì)數(shù)據(jù),結(jié)合2020年土地利用數(shù)據(jù),分析土地利用與景觀格局對開都河水質(zhì)的影響,探討不同時空尺度下的開都河水質(zhì)對土地利用與景觀格局的相應(yīng)關(guān)系,找出影響水質(zhì)的最佳緩沖區(qū)和關(guān)鍵變量,以期為開都河水質(zhì)管理提供理論依據(jù)。

    1" "材料與方法

    1.1" "研究區(qū)概況與數(shù)據(jù)采集

    開都河流域位于中國西北的天山南坡中部(82°58′~86°05′ E,42°14′~43°21′ N),發(fā)源于海拔5 000 m的薩爾明山,流域面積22 000 km2,全長560 km(Ba et al,2018;Zhao et al,2021),海拔871~5 000 m,南部與塔里木盆地接壤,中部為低山丘陵、沖積扇戈壁區(qū),地形復(fù)雜,交叉分布。該流域氣候類型為大陸性溫帶干旱氣候,年平均氣溫為-1.82 ℃,年總降水量為328.8 mm,年平均降水量61.2 mm(李茹霞等,2023),冬、夏兩季降雨量季節(jié)差異顯著,水位在夏季較高,為豐水期,冬季為枯水期。開都河流經(jīng)尤爾都斯和焉耆盆地,最終匯入博斯騰湖。

    本研究在開都河取15個采樣點(diǎn)(圖1),下游K1~K8、中游K9~K10、上游K11~K15,2021年6月(豐水期)與10月(枯水期)進(jìn)行樣品采集,水樣選取50 cm深度的流動水,每個點(diǎn)位采集3個平行樣本進(jìn)行分析。用500 mL潔凈聚乙烯塑料瓶采集水樣,在0~4 ℃冷藏密封帶回實(shí)驗(yàn)室。使用重鉻酸鹽法測定化學(xué)需氧量(COD),鉬酸氨分光光度法測定總磷(TP),納氏試劑分光光度法測定氨氮(NH3-N),以堿性過硫酸鉀消解紫外分光光度法測定總氮(TN),多水質(zhì)參數(shù)檢測儀測定溶解氧(DO)與濁度。以上操作均依據(jù)《國家水質(zhì)采樣技術(shù)指導(dǎo)》(HJ494—2009)。

    開都河2020年30 m×30 m精度土地利用遙感監(jiān)測數(shù)據(jù)源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http://www.resdc.cn)。按照國家土地利用現(xiàn)狀分類標(biāo)準(zhǔn)(GB/T21010—2017)并結(jié)合開都河土地利用實(shí)況,將研究區(qū)土地利用類型劃分為耕地、林地、草地、水域、建設(shè)用地和未利用地6類。以采樣點(diǎn)為圓心的圓形緩沖區(qū)作為分析單元,根據(jù)已有水質(zhì)與土地利用尺度效應(yīng)的相關(guān)研究為依據(jù)確定緩沖區(qū)半徑(張志敏等,2022;溫嘉偉等,2023),作200、500、1 000、2 000和5 000 m緩沖區(qū),并與土地利用類型數(shù)據(jù)疊加,獲取不同空間尺度范圍內(nèi)的土地利用類型面積占比。

    1.2" "景觀指數(shù)選取

    本文從景觀破碎度、聚集度、優(yōu)勢度、多樣性和物理連接度等景觀指數(shù)中,選取了7個通用性高、代表性較好的指數(shù)(表1),其中NP、PD、LSI分別表示景觀中斑塊的數(shù)量、密度、形狀,通常用來評估景觀破碎化的程度(徐啟渝等,2020);COHESION表示景觀中同一類型斑塊的聚集程度和給定閾值的斑塊連接度;LPI表示景觀中面積最大斑塊及其面積占比程度,用來形容流域內(nèi)優(yōu)勢景觀類型;SHDI、ENN_MN分別表示流域內(nèi)景觀類型的豐富度以及同類型斑塊之間的平均距離。

    1.3" "水質(zhì)評價

    內(nèi)梅羅污染指數(shù)法是評價河流湖泊水質(zhì)的方法,本文依照《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838—2002)中Ⅲ類水質(zhì)標(biāo)準(zhǔn)進(jìn)行計算,綜合分析各評價指標(biāo),具體方法見參考文獻(xiàn)(李梁婷等,2022),評價標(biāo)準(zhǔn)見表2。

    1.4" "數(shù)據(jù)分析

    采用Arcgis10.8和Fragstats4.2軟件分別計算各土地利用類型面積和景觀格局指數(shù),數(shù)據(jù)通過Kolmogorov-Smirnov檢驗(yàn)正態(tài)性,Pgt;0.05表明數(shù)據(jù)符合正態(tài)分布特征。通過Canoco5軟件進(jìn)行冗余分析(RDA),判定開都河流域土地利用與景觀格局對水質(zhì)的影響。借助R語言“vegan”包中的“bioenv”函數(shù),分別得到土地利用類型和景觀指數(shù)的最佳組合,并以此來表征土地利用結(jié)構(gòu)和景觀格局;通過Mantel檢驗(yàn)分析土地利用和景觀格局與水質(zhì)指標(biāo)相關(guān)性,方差分解(variation partitioning)是將響應(yīng)變量組成數(shù)據(jù)表的總方差無偏分解成由各個解釋變量所決定的子方差,使用originPro2021進(jìn)行繪圖。

    2" "結(jié)果與分析

    2.1" "水質(zhì)指標(biāo)特征

    單因子污染評價結(jié)果(表3)顯示,開都河上、中、下游皆未超出Ⅲ類水標(biāo)準(zhǔn)。內(nèi)梅羅污染指數(shù)評價結(jié)果顯示:豐水期開都河上、中、下游均屬較清潔水質(zhì),枯水期上游、下游屬較清潔水質(zhì),而中游屬清潔水質(zhì)。豐水期各采樣點(diǎn)水質(zhì)指標(biāo)濃度分別為DO(7.45~10.01 mg/L)、TP(0.01~0.10 mg/L)、TN(0.41~0.78 mg/L)、NH3-N(0.16~0.98 mg/L)、COD(11.23~29.22 mg/L)、濁度(0.31~4.13 UTN)、pH(7.02~7.86);枯水期水質(zhì)指標(biāo)濃度分別為DO(6.04~9.93 mg/L)、TP(0.02~0.08 mg/L)、TN(0.219~0.45 mg/L)、NH3-N(0.11~0.76 mg/L)、COD(5.36~22.43 mg/L)、濁度(0.03~8.28 UTN)、pH(7.02~8.10)。豐水期濁度和COD超標(biāo)率分別為33.33%、26.67%,但其在枯水期的超標(biāo)率均為13.33%。TP除豐水期K4采樣點(diǎn),其余采樣點(diǎn)均未超出Ⅲ類水標(biāo)準(zhǔn)。NH3-N、TN在各采樣點(diǎn)(枯水期、豐水期)均未超標(biāo)。

    2.2" "土地利用類型對水質(zhì)影響的緩沖區(qū)半徑

    不同空間尺度下土地利用類型面積占比如圖2所示??傮w而言,隨著緩沖區(qū)半徑的擴(kuò)大,水域面積占比逐漸減小,未利用地與耕地面積占比逐漸增加;建設(shè)用地與林地面積在所有緩沖區(qū)內(nèi)都維持較小占比;草地面積在所有緩沖區(qū)內(nèi)都維持較高占比。各尺度主要景觀類型均由草地與未利用地組成。在自然因素與人類活動的影響下,土地利用比例在不同采樣點(diǎn)具有空間分異性。

    利用RDA分析探討影響開都河流域的最佳空間尺度,發(fā)現(xiàn)豐、枯水期半徑緩沖區(qū)分別為2 000、5 000 m內(nèi)的土地利用方式對水質(zhì)的解釋率最高(表4),分別為63.81%和64.91%,說明在該空間尺度內(nèi)的土地利用類型對開都河水質(zhì)的影響最為顯著。

    2.3" "不同空間尺度景觀異質(zhì)性

    隨著空間尺度的增加(圖3),NP、LSI、ENN_NM、COHESION、SHDI逐漸增大,PD、LPI逐漸減小。NP增大,景觀斑塊數(shù)量變多,土地利用破碎度增大,空間異質(zhì)性增強(qiáng);LSI表示景觀斑塊形狀的規(guī)則化程度,本研究中隨空間尺度的增加LSI呈上升趨勢,說明景觀形狀復(fù)雜性隨空間尺度增加而增強(qiáng),形狀多樣性增大。ENN_NM代表平均最近鄰體距離,其值越大,代表同類景觀之間的距離越遠(yuǎn)。COHESION代表斑塊結(jié)合度,其值越大,景觀的空間連通性越好。SHDI屬于敏感性指標(biāo),其值越高表示土地利用類型越豐富,景觀類型趨向均衡化;本研究中隨著空間尺度的增加,SHDI呈逐漸上升的趨勢,表明尺度空間越大,其景觀多樣性越高,空間異質(zhì)性越強(qiáng)。PD代表斑塊個數(shù),其值隨空間尺度的增加而逐漸變小,表明景觀破碎化程度隨緩沖區(qū)的增大而減小,緩沖區(qū)半徑越大,景觀異質(zhì)性越小。LPI最大斑塊,其值越小,說明隨著緩沖區(qū)的增大,人類活動對景觀的干擾逐漸減小。利用RDA分析,探討影響開都河流域的最佳空間尺度,發(fā)現(xiàn)豐、枯水期緩沖區(qū)半徑分別為500、200 m內(nèi)的景觀格局對水質(zhì)的解釋率最高(表5),分別為73.88%和63.33%,說明在該空間尺度內(nèi)的景觀格局對開都河水質(zhì)的影響最顯著。

    2.4" "土地利用與景觀格局對水質(zhì)的影響

    2.4.1" "土地利用類型對水質(zhì)的影響" "圖4為不同土地利用類型對水質(zhì)各參數(shù)影響的二位排序結(jié)果,其與水質(zhì)指標(biāo)之間形成的夾角表示相關(guān)性,角度越小,相關(guān)性越強(qiáng),當(dāng)夾角小于90°時,兩者為正相關(guān),接近或等于90°時,兩者相關(guān)性較小或者不相關(guān)(周俊菊等,2019)。豐水期2 000 m半徑緩沖區(qū)內(nèi),耕地與NH3-N、COD、TN、TP、濁度呈正相關(guān);未利用地與耕地結(jié)果相反;草地與DO呈正相關(guān),與COD、濁度、NH3-N、TN呈負(fù)相關(guān),與TP相關(guān)性不強(qiáng);水域與TN、NH3-N呈正相關(guān),與濁度、DO、TP、COD呈負(fù)相關(guān);建設(shè)用地與COD的相關(guān)性最強(qiáng),與NH3-N、TN、濁度、TP呈正相關(guān),與DO呈負(fù)相關(guān);林地與COD相關(guān)性不強(qiáng),與TN、NH3-N呈正相關(guān),與濁度、TP、DO呈負(fù)相關(guān)(圖4a)??菟? 000 m半徑緩沖區(qū)內(nèi),草地和未利用地面積比例與TP、濁度均呈正相關(guān),與DO、TN、COD、NH3-N呈負(fù)相關(guān);耕地和建設(shè)用地與COD、NH3-N、TN、DO呈正相關(guān),與濁度、TP呈負(fù)相關(guān);水域與濁度、DO、NH3-N、TN呈正相關(guān),與TP、COD呈負(fù)相關(guān);林地面積比例與未利用地及草地的結(jié)果相反(圖4b)。由上述分析可知,在6種土地利用類型中,無論是豐水期還是枯水期,建設(shè)用地均與各水質(zhì)指標(biāo)之間存在一定程度的正相關(guān)性,表明建設(shè)用地對水質(zhì)的影響最大,是主要的污染輸出源;草地與未利用地在不同水期與大部分水質(zhì)指標(biāo)之間呈現(xiàn)負(fù)相關(guān)性。

    豐水期各土地利用對水質(zhì)的貢獻(xiàn)率見圖5,耕地對水質(zhì)的貢獻(xiàn)率最高為57.9%,其次為林地,對水質(zhì)的貢獻(xiàn)率為17.6%;枯水期對水質(zhì)貢獻(xiàn)率最高的為草地,其對水質(zhì)的貢獻(xiàn)率為68%,其次是水域,對水質(zhì)的貢獻(xiàn)率為11%。豐水期土地利用對水質(zhì)的貢獻(xiàn)率大小為:耕地gt;林地gt;水域gt;未利用地gt;草地gt;建設(shè)用地;枯水期土地利用對水質(zhì)的貢獻(xiàn)率大小為:草地gt;水域gt;建設(shè)用地gt;林地gt;未利用地gt;耕地。

    2.4.2" "景觀格局指數(shù)與水質(zhì)指標(biāo)的相關(guān)性" "不同景觀格局指數(shù)對水質(zhì)各參數(shù)影響的二維排序結(jié)果(圖6)顯示,豐水期500 m半徑緩沖區(qū)內(nèi),LSI、NP、PD與NH3-N、COD、TP呈正相關(guān),與濁度、TN相關(guān)性不強(qiáng),與DO呈負(fù)相關(guān);COHESION與LSI、NP、PD對水質(zhì)的影響結(jié)果相反;LPI與DO呈正相關(guān),與TP相關(guān)性不強(qiáng),與COD、NH3-N、TN、濁度呈負(fù)相關(guān);ENN_M(jìn)N與LPI對水質(zhì)影響的結(jié)果相反,其與濁度的相關(guān)性最強(qiáng);SHDI與TN、濁度呈負(fù)相關(guān),與NH3-N、COD、TP呈正相關(guān)??菟?00 m半徑緩沖區(qū)內(nèi),LPI、SHDI與濁度、TP呈正相關(guān),與DO相關(guān)性不強(qiáng),與COD、NH3-N、TN呈負(fù)相關(guān);NP、PD、LSI與TP、DO、COD、NH3-N、TN呈正相關(guān),與濁度呈負(fù)相關(guān);COHESION與NP對水質(zhì)的影響結(jié)果相反。由上述分析可知,在7種景觀格局指數(shù)中,無論是豐水期還是枯水期,NP、PD、LSI均與各水質(zhì)指標(biāo)之間存在一定程度的正相關(guān)性,COHESION與各水質(zhì)指標(biāo)之間均存在一定程度的負(fù)相關(guān)性;景觀的空間連通性越好,對研究區(qū)水質(zhì)的改善則會越好。

    豐水期各土地利用對水質(zhì)的貢獻(xiàn)率見圖7,LPI對水質(zhì)的貢獻(xiàn)率最高為48.8%,其次為COHESION,對水質(zhì)的貢獻(xiàn)率為18.0%;枯水期對水質(zhì)貢獻(xiàn)率最高的為LSI,其對水質(zhì)的貢獻(xiàn)率為21.6%,其次是COHESION,對水質(zhì)的貢獻(xiàn)率為15.5%。豐水期景觀格局對水質(zhì)的貢獻(xiàn)率大小均為:LPIgt;COHESIONgt;PDgt;SHDIgt;NPgt;LSI>ENN_M(jìn)N??菟诰坝^格局對水質(zhì)的貢獻(xiàn)率大小均為:LSIgt;COHESIONgt;LPIgt;NPgt;ENN_M(jìn)Ngt;SHDI>PD。

    2.5" "水質(zhì)與土地利用結(jié)構(gòu)和景觀指數(shù)組合的相關(guān)關(guān)系

    根據(jù)“bioenv”分析結(jié)果(表6),豐水期2 000 m半徑緩沖區(qū)尺度上,耕地、林地、水域的組合與水質(zhì)的相關(guān)系數(shù)最高為0.424 5;枯水期5 000 m半徑緩沖區(qū)尺度上,單一草地與水質(zhì)的相關(guān)系數(shù)最高為0.512 9。豐水期500 m半徑緩沖區(qū)尺度上,NP、LPI、ENN_MN與水質(zhì)的相關(guān)系數(shù)最高為0.076 5;枯水期200 m半徑緩沖區(qū)尺度上,NP、PD、LPI與水質(zhì)的相關(guān)系數(shù)最高為0.216 3。

    2.6" "土地利用結(jié)構(gòu)和景觀格局對水質(zhì)的交互影響

    Mantel檢驗(yàn)分析結(jié)果表明(表7),豐水期500、2 000 m半徑緩沖區(qū)范圍內(nèi)土地利用和景觀格局與水質(zhì)的相關(guān)性均達(dá)到顯著水平(Plt;0.05),其中2 000 m半徑緩沖區(qū)優(yōu)于500 m半徑緩沖區(qū)??菟?00、5 000 m半徑緩沖區(qū)范圍內(nèi)土地利用與水質(zhì)的相關(guān)性顯著,且200 m半徑緩沖區(qū)與水質(zhì)的相關(guān)性優(yōu)于5 000 m半徑緩沖區(qū),而景觀格局與水質(zhì)的相關(guān)性均不顯著??傮w而言,不同緩沖區(qū)、不同水期土地利用與水質(zhì)的相關(guān)性均優(yōu)于景觀格局。

    對豐水期2 000 m半徑緩沖區(qū)和枯水期200 m半徑緩沖區(qū)內(nèi)土地利用和景觀格局與水質(zhì)進(jìn)行方差分解分析,結(jié)果見圖8。土地利用結(jié)構(gòu)與景觀格局對水質(zhì)的交互作用在豐水期2 000 m半徑緩沖區(qū)內(nèi)為18%,其中土地利用單獨(dú)解釋了16%的水質(zhì)變化,景觀格局單獨(dú)解釋了9%的水質(zhì)變化。交互作用部分分別占土地利用結(jié)構(gòu)總貢獻(xiàn)率(總貢獻(xiàn)率=單獨(dú)解釋+交互作用)的52.94%、占景觀格局總貢獻(xiàn)率的66.67%??菟?00 m半徑緩沖區(qū)范圍內(nèi)土地利用單獨(dú)解釋了10%的水質(zhì)變化,景觀格局單獨(dú)解釋了6%的水質(zhì)變化,交互作用部分分別占土地利用總貢獻(xiàn)率的52.38%、占景觀格局總貢獻(xiàn)率的65.71%。可見土地利用結(jié)構(gòu)與景觀格局的交互作用是水質(zhì)變化的重要因素。

    3" "討論

    3.1" "土地利用類型對水質(zhì)的影響

    不同的土地利用類型通過影響流域特征,如水文循環(huán)、生物多樣性、土壤侵蝕、污染物遷移轉(zhuǎn)化及人類活動等最終影響水質(zhì)變化,但二者之間的最佳緩沖區(qū)尚未形成統(tǒng)一結(jié)論(項(xiàng)頌等,2016)。研究區(qū)的差異、土地利用類型的劃分及數(shù)據(jù)集的分辨率不同皆會導(dǎo)致研究結(jié)果不同。本研究通過RDA分析發(fā)現(xiàn)豐水期土地利用對水質(zhì)的解釋率最高的緩沖區(qū)為2 000 m,這與易帆等(2023)在大渡河流域的研究結(jié)果一致;枯水期土地利用對水質(zhì)最佳解釋率為5 000 m半徑緩沖區(qū),而夏明珠等(2024)發(fā)現(xiàn)北京白河流域土地利用對水質(zhì)的最佳緩沖區(qū)為1 000 m,代孟均等(2024)以旱雨季的平均解釋率分析得出,3 000 m緩沖區(qū)是海河流域天津段土地利用對水質(zhì)指標(biāo)的最佳影響尺度。結(jié)論的不一致,可能與研究區(qū)概況相關(guān),開都河位于干旱區(qū),且流域內(nèi)海拔差異顯著,景觀類型復(fù)雜,同時緩沖區(qū)的尺度選擇也是造成差異的一個原因。已有研究大多建立在5 000 m半徑緩沖區(qū)內(nèi),而緩沖區(qū)與緩沖區(qū)之間的間隔選擇存在多樣性。

    開都河2 000 m與5 000 m半徑緩沖內(nèi)以草地、未利用地與耕地為主(圖2)。草地主要集中在開都河的中游和上游(K9~K15采樣點(diǎn)),其與TN、COD、NH3-N呈負(fù)相關(guān)(圖4),這主要是由于草地屬于自然景觀類型,徑流匯入河網(wǎng)的過程中自然景觀類型起過濾作用,這使開都河水質(zhì)處于較清潔水平。耕地主要集中開都河的下游(K1~K8),與NH3-N、COD、TN呈正相關(guān),這是由于耕作和施肥等農(nóng)業(yè)活動,導(dǎo)致土壤中氮和農(nóng)藥等殘留含量高(Mishra et al,2021),而未被吸收的養(yǎng)分和有機(jī)質(zhì)等經(jīng)過降雨徑流沖刷后形成污染源進(jìn)入附近的河流(李紀(jì)周,2011),從而導(dǎo)致水體水質(zhì)下降(Mello et al,2020;Wei et al,2020;Zhang et al,2020)。未利用地主要集中在開都河上、中游,在不同水期與水質(zhì)指標(biāo)濁度呈正效應(yīng),與其他水質(zhì)指標(biāo)表現(xiàn)為負(fù)效應(yīng),這主要由于開都河流域未利用地主要為沙漠、戈壁,風(fēng)和徑流將泥沙帶入河流使得河流變渾濁,同時沙礫等類型的未利用地由于受人類干擾弱,對水質(zhì)有凈化效果。

    3.2" "景觀格局對水質(zhì)的影響

    景觀格局在河流生態(tài)水文過程、物質(zhì)循環(huán)和能量轉(zhuǎn)換等過程中發(fā)揮了重要作用(徐明珠等,2023)。從空間上來看,景觀格局對水質(zhì)的影響在空間尺度上有較大依賴性,每條河流均存在一個臨界緩沖區(qū),即流域內(nèi)景觀格局和水質(zhì)指標(biāo)之間影響關(guān)系最大的緩沖區(qū),本研究表明在200、500 m半徑緩沖區(qū)內(nèi)景觀格局對水質(zhì)的解釋能力最強(qiáng),這與米秋菊等(2024)、劉可暄等(2022)結(jié)果相似,而朱穎等(2021)認(rèn)為在陽澄湖小流域景觀格局與水質(zhì)指標(biāo)最相關(guān)的空間尺度為800 m半徑緩沖區(qū),這主要是由于景觀格局指數(shù)的有效性不僅受尺度效應(yīng)、生態(tài)意義可解釋性,還受土地分類的不確定性影響,導(dǎo)致不同的研究者會有不同的結(jié)論。同時,土地利用對水質(zhì)的影響側(cè)重于土地單元的多樣性,而景觀格局側(cè)重于空間構(gòu)型特征(劉金寶等,2023),其對空間變化更敏感,能在較小的空間尺度上反映對水質(zhì)的影響(Wu et al,2002),這也造成了本研究中景觀格局和土地利用對河流水質(zhì)影響最優(yōu)緩沖區(qū)不一致的原因。景觀格局對于水質(zhì)的影響仍具有很大的不確定性,對于最優(yōu)空間尺度目前仍無統(tǒng)一定論,有待進(jìn)一步探討。

    本研究表明豐水期500 m半徑緩沖區(qū)內(nèi),與水質(zhì)表現(xiàn)為負(fù)效應(yīng)的LPI為主要影響因子,其反映最大斑塊占整個景觀面積的比例,傳達(dá)出優(yōu)勢類的信息,已有研究表明當(dāng)優(yōu)勢景觀類型為“源”景觀,則LPI可能與水質(zhì)污染程度呈正相關(guān),若為“匯”景觀,則有利于水質(zhì)的提高(楊強(qiáng)強(qiáng)等,2020),開都河流域500 m緩沖區(qū)內(nèi)草地占比大,作為“匯”景觀,對開都河水質(zhì)凈化效果顯著;枯水期200 m半徑緩沖區(qū)內(nèi),LSI為主要的影響因子,與水質(zhì)表現(xiàn)為正效應(yīng),其反映了地表景觀的破碎度和復(fù)雜度,景觀破碎度越高,受人類干擾活動越大,斑塊破碎也有利于污染物的遷移輸出,會加劇流域水質(zhì)惡化。

    3.3" "土地利用和景觀格局與水質(zhì)的交互作用

    本研究發(fā)現(xiàn)在豐水期2 000 m、枯水期200 m半徑緩沖區(qū)內(nèi),土地利用和景觀格局對水質(zhì)交互作用顯著,但未解釋率占比均較大。大部分污染物通過徑流匯入河流,流域內(nèi)的徑流對河流水質(zhì)產(chǎn)生重要影響(鄭一和王學(xué)軍,2002;Zhang et al,2011)。這一產(chǎn)匯流過程受到土地利用結(jié)構(gòu)和景觀格局之間的交互作用影響(Ouyang et al,2010;Jiang et al,2017)。而影響產(chǎn)流過程的關(guān)鍵因素包括降雨強(qiáng)度、前期土壤含水量、土壤物理性質(zhì)以及植被覆蓋比例(Liu et al,2014)。土壤的特性,如飽和土壤含水量、粒度和孔隙度,對水體下滲和壤中流生成機(jī)制有直接影響。不同的土地利用類型在污染物的“源”和“匯”方面發(fā)揮不同作用,導(dǎo)致了不同土地利用類型產(chǎn)生不同水質(zhì)的徑流(Germer et al,2010)。景觀格局也在匯流過程中與水質(zhì)相關(guān)性顯著(潘釗等,2018)。自然景觀類型如林地和草地可截留泥沙、重金屬和有機(jī)物等污染物,凈化水質(zhì)。而人工景觀類型如耕地和建設(shè)用地可能富集污染物(Shigaki et al,2007)。土壤是地表污染物的主要承載體,影響水體下滲和壤中流的產(chǎn)生機(jī)制。壤中流在流動過程中受到自然景觀類型的截留作用,有助于吸附和固定污染物。但在流經(jīng)人工景觀時,由于下墊面的不同,污染物可能改變流向。這些因素共同導(dǎo)致了水質(zhì)的變化,因此開都河土地利用結(jié)構(gòu)和景觀格局的交互作用是水質(zhì)變化的重要原因。圖8中水質(zhì)變化的未解釋部分可能是未考慮土壤性質(zhì)差異造成的,已有研究表明,即使是相同土地利用類型,不同的土壤特征也會對徑流形成產(chǎn)生較大影響,尤其對壤中流的產(chǎn)流機(jī)制影響巨大(Li et al,2018),且土地利用類型的相鄰區(qū)域受到周邊多種土地利用共同作用,對污染物的源匯機(jī)制難以分析可能是未解釋占比高的原因之一。

    參考文獻(xiàn)

    曹燦, 張飛, 阿依尼格爾·亞力坤, 等, 2018. 艾比湖區(qū)域景觀格局與河流水質(zhì)關(guān)系探討[J]. 環(huán)境科學(xué), 39(4):1568-1577.

    CAO C, ZHANG F, AYINIGEER Y, et al, 2018. Relationship between landscape pattern and water quality in the Ebinur lake region[J]. Environmental Science, 39(4):1568-1577.

    崔丹, 陳巖, 馬冰然, 等, 2019. 土地利用/景觀格局對水環(huán)境質(zhì)量的影響[J]. 水科學(xué)進(jìn)展, 30(3):423-433.

    CUI D, CHEN Y, MA B R, et al, 2019. Effects of land use/landscape patterns on the water quality[J]. Advances in Water Science, 30(3):423-433.

    代孟均, 張兵, 杜倩倩, 等, 2024. 不同緩沖區(qū)的土地利用方式對地表水水質(zhì)的影響:以海河流域天津段為例[J]. 環(huán)境科學(xué), 45(3):1512-1524.

    DAI M J, ZHANG B, DU Q Q, et al, 2024. Effects of land use types on water quality at different buffer scales: Tianjin section of the Haihe river basin as an example[J]. Environmental Science, 45(3):1512-1524.

    郭羽羽, 李思悅, 劉睿, 等, 2021. 黃河流域多時空尺度土地利用與水質(zhì)的關(guān)系[J]. 湖泊科學(xué), 33(3):737-748.

    GUO Y Y, LI S Y, LIU R, et al, 2021. Relationship between landscape pattern and water quality of the multi-scale effects in the Yellow River basin[J]. Journal of Lake Sciences, 33(3):737-748.

    胡叢巧, 迪麗努爾·阿吉, 李茹霞, 等, 2023. 博斯騰湖不同時空尺度下土地利用景觀格局對水質(zhì)的影響[J/OL]. 水生態(tài)學(xué)雜志:1-11[2023-11-09].https://doi.org/10.15928/j.1674-3075.202211020436.

    HU C Q, DILINUER·A, LI R X, et al, 2023. Effects of land use landscape pattern on water quality in Bosten lake at different temporal and spatial scales[J/OL]. Journal of Hydroecology:1-11[2023-11-09]. https://doi. org/10.15928/j.1674-3075.202211020436.

    李紀(jì)周, 2011. 天津市規(guī)?;笄蒺B(yǎng)殖場糞污治理及資源化利用調(diào)查研究[D]. 北京:中國農(nóng)業(yè)科學(xué)院.

    李梁婷, 邵蓓, 卓瑪曲西, 等, 2022. 西藏拉魯濕地水質(zhì)綜合狀況時空變化特征及污染源分析[J]. 安徽農(nóng)業(yè)科學(xué), 50(12):67-74,78.

    LI L T, SHAO B, ZHUOMAQUXI, et al, 2022. Analysis of spatial-temporal variation characteristics and pollution source of water quality comprehensive conditions in Lalu wetland, Tibet[J]. Journal of Anhui Agricultural Sciences, 50(12):67-74,78.

    李明濤, 王曉燕, 劉文竹, 2013. 潮河流域景觀格局與非點(diǎn)源污染負(fù)荷關(guān)系研究[J]. 環(huán)境科學(xué)學(xué)報, 33(8):2296-2306.

    LI M T, WANG X Y, LIU W Z, et al, 2013. Relationship between landscape pattern and non-point source pollution loads in the Chaohe river watershed[J]. Acta Scientiae Circumstantiae, 33(8):2296-2306.

    李茹霞, 迪麗努爾·阿吉, 賽米熱·托合提, 等, 2023. 不同空間尺度土地利用方式對開都河水質(zhì)的影響[J]. 海洋湖沼通報, 45(1):98-106.

    LI R X, DILINUER·A, SAIMIRE·T, et al, 2023. Effects of land use patterns at different spatial scales on water quality of Kaidu river[J]. Transactions of Oceanology and Limnology, 45(1):98-106.

    劉金寶, 劉軒, 孫增慧, 等, 2023. 丹江口水庫土地利用與水質(zhì)變化響應(yīng)關(guān)系研究[J]. 遙感技術(shù)與應(yīng)用, 38(2):285-296.

    LIU J B, LIU X, SUN Z H, et al, 2023. Study on the relationship between land use and water quality change in Danjiangkou reservoir[J]. Remote Sensing Technology and Application, 38(2):285-296.

    劉可暄, 王冬梅, 常國梁, 等, 2022. 多空間尺度景觀格局與地表水質(zhì)響應(yīng)關(guān)系研究[J]. 環(huán)境科學(xué)學(xué)報, 42(2):23-31.

    LIU K X, WANG D M, CHANG G L, et al, 2022. Study on the relationship between landscape pattern and surface water quality at multiple spatial scales[J]. Acta Scientiae Circumstantiae, 42(2):23-31.

    米秋菊, 邊玉明, 甄立鵬, 等, 2024. 土地利用結(jié)構(gòu)與景觀格局對黎河水質(zhì)的影響[J]. 水生態(tài)學(xué)雜志, 45(3):78-85.

    MI Q J, BIAN Y M, ZHEN L P, et al, 2024. Impacts of land use structure and landscape pattern on water quality in Lihe river[J]. Journal of Hydroecology, 45(3):78-85.

    潘釗, 孫自永, 馬瑞, 等, 2018. 黑河上游高寒山區(qū)降雨—徑流形成過程的同位素示蹤[J]. 地球科學(xué), 43(11):4226-4236.

    PAN Z, SUN Z Y, MA R, et al, 2018. Isotopic investigation of rainfall-runoff generation in an alpine catchment in headwater regions of Heihe river, northeast Qinghai-Tibet plateau[J]. Earth Science, 43(11):4226-4236.

    祁蘭蘭, 王金亮, 葉輝, 等, 2021. 滇中“三湖流域”土地利用景觀格局與水質(zhì)變化關(guān)系研究[J]. 水土保持研究, 28(6):199-208.

    QI L L, WANG J L, YE H, et al, 2021. Study on the land use landscape pattern and water quality change in the ‘three lakes’ basin in central Yunnan[J]. Research of Soil and Water Conservation, 28(6):199-208.

    王小平, 張飛, 李曉航, 等, 2017. 艾比湖區(qū)域景觀格局空間特征與地表水質(zhì)的關(guān)聯(lián)分析[J]. 生態(tài)學(xué)報, 37(22):7438-7452.

    WANG X P, ZHANG F, LI X H, et al, 2017. Correlation analysis between the spatial characteristics of land use/cover-landscape pattern and surface-water quality in the Ebinur lake area[J]. Acta Ecologica Sinica, 37(22):7438-7452.

    溫嘉偉, 王鵬, 黃高翔, 等, 2023. 土地利用結(jié)構(gòu)與空間格局對鄱陽湖流域中小河流水質(zhì)的影響[J]. 環(huán)境科學(xué), 44(12):6728-6743.

    WEN J W, WANG P, HUANG G X, et al, 2023. Influence of land use structure and spatial pattern on water quality of small and medium-sized rivers in Poyang lake basin[J]. Environmental Science, 44(12):6728-6743.

    夏明珠, 薛萬來, 黃俊雄, 等, 2024. 北京白河流域不同空間尺度下水質(zhì)對土地利用的響應(yīng)關(guān)系[J]. 環(huán)境科學(xué)研究, 37(4):764-775.

    XIA M Z, XUE W L, HUANG J X, et al, 2024. Relationship between water quality and land use responses at different spatial scales in white river basin in Beijing[J]. Research of Environmental Sciences, 37(4):764-775.

    項(xiàng)頌, 龐燕, 儲昭升, 等, 2016. 入湖河流水質(zhì)對土地利用時空格局的響應(yīng)研究:以洱海北部流域?yàn)槔齕J]. 環(huán)境科學(xué), 37(8):2947-2956.

    XIANG S, PANG Y, CHU Z S, et al, 2016. Response of inflow water quality to land use pattern in northern watershed of lake Erhai[J]. Environmental Science, 37(8):2947-2956.

    項(xiàng)頌, 龐燕, 竇嘉順, 等, 2018. 不同時空尺度下土地利用對洱海入湖河流水質(zhì)的影響[J]. 生態(tài)學(xué)報, 38(3):876-885.

    XIANG S, PANG Y, DOU J S, et al, 2018. Impact of land use on the water quality of inflow river to Erhai lake at different temporal and spatial scales[J]. Acta Ecologica Sinica, 38(3):876-885.

    徐明珠, 徐國策, 喬海亮, 等, 2023. 秦嶺南麓小流域不同空間尺度景觀格局對水質(zhì)的影響分析[J]. 環(huán)境科學(xué)學(xué)報, 43(10):396-406.

    XU M Z, XU G C, QIAO H L, et al, 2023. Influence of different spatial scale landscape patterns on water quality in a small watershed at the southern foot of the Qinling mountains[J]. Acta Scientiae Circumstantiae, 43(10):396-406.

    徐啟渝, 王鵬, 王濤, 等, 2020. 土地利用結(jié)構(gòu)與景觀格局對鄱陽湖流域贛江水質(zhì)的影響[J]. 湖泊科學(xué), 32(4):1008-1019.

    XU Q Y, WANG P, WANG T, et al, 2020. Investigation of the impacts of land use structure and landscape pattern on water quality in the Ganjiang river, lake Poyang basin[J]. Journal of Lake Sciences, 32(4):1008-1019.

    楊強(qiáng)強(qiáng), 徐光來, 楊先成, 等, 2020. 青弋江流域土地利用/景觀格局對水質(zhì)的影響[J]. 生態(tài)學(xué)報, 40(24):9048-9058.

    YANG Q Q, XU G L, YANG X C, et al, 2020. Responses of water quality to land use amp; landscape pattern in the Qingyijiang river watershed[J]. Acta Ecologica Sinica, 40(24):9048-9058.

    易帆, 陳旻, 何曉楓, 等, 2023. 流域土地利用分析中空間尺度對水質(zhì)的影響[J]. 中國環(huán)境科學(xué),43(8):4280-4291.

    YI F, CHEN M, HE X F, et al, 2023. Effect of spatial scale on water quality in watershed land use analysis[J]. China Environmental Science, 43(8):4280-4291.

    張宇碩, 吳殿廷, 呂曉, 2020. 土地利用/覆蓋變化對生態(tài)系統(tǒng)服務(wù)的影響:空間尺度視角的研究綜述[J]. 自然資源學(xué)報, 35(5):1172-1189.

    ZHANG Y S, WU D T, LYU X, et al, 2020. A review on the impact of land use/land cover change on ecosystem services from a spatial scale perspective[J]. Journal of Natural Resources, 35(5):1172-1189.

    張志敏, 杜景龍, 陳德超, 等, 2022. 典型網(wǎng)狀河網(wǎng)區(qū)域土地利用和景觀格局特征對地表季節(jié)水質(zhì)的影響:以江蘇省溧陽市為例[J]. 湖泊科學(xué), 34(5):1524-1540.

    ZHANG Z M, DU J L, CHEN D C, et al, 2022. Effects of land use and landscape pattern characteristics on seasonal surface water quality in a typical reticulated river network area: a case study of Liyang city, Jiangsu Province[J]. Journal of Lake Sciences, 34(5):1524-1540.

    鄭一, 王學(xué)軍, 2002. 非點(diǎn)源污染研究的進(jìn)展與展望[J]. 水科學(xué)進(jìn)展, (1):105-110.

    ZHENG Y, WANG X J, 2002. Advances and prospects for nonpoint source pollution studies[J]. Advances in Water Science, (1):105-110.

    周俊菊, 向鵑, 王蘭英, 等, 2019. 祁連山東部冰溝河流域景觀格局與河流水化學(xué)特征關(guān)系[J]. 生態(tài)學(xué)雜志, 38(12):3779-3788.

    ZHOU J J, XIANG J, WANG L Y, et al, 2019. Relationship between landscape pattern and hydrochemical characteristics of Binggou river basin in eastern Qilian mountains[J]. Chinese Journal of Ecology, 38(12):3779-3788.

    朱穎, 王怡, 賀風(fēng)春, 等, 2021. 景觀格局對陽澄湖小流域水質(zhì)變化影響的空間尺度效應(yīng)[J]. 水土保持通報, 41(2):105-113, 353.

    ZHU Y, WANG Y, HE F C, et al, 2021. Spatial scale effect of landscape pattern on water quality change in Yangcheng lake watershed[J]. Bulletin of Soil and Water Conservation, 41(2):105-113,353.

    ALLAN J D, 2004. Landscapes and riverscapes: the Influence of land use on stream ecosystems[J]. Annual Review of Ecology Evolution and Systematics, 35(1):257-284.

    BA W, DU P, LIU T, et al, 2018. Simulating hydrological responses to climate change using dynamic and statistical downscaling methods: a case study in the Kaidu River Basin, Xinjiang, China[J]. Journal of Arid Land, 10(6):905-920.

    BASNYAT P, TEETER L D, FLYNN K M, et al, 1999. Relationships between landscape characteristics and nonpoint source pollution inputs to coastal estuaries[J]. Environmental Management, 23(4):539-549.

    BU H, MENG W, ZHANG Y, et al, 2014. Relationships between land use patterns and water quality in the Taizi River basin China[J]. Ecological Indicators, 41:187-197.

    CASTRILLO M, GARCíA á L, 2020. Estimation of high frequency nutrient concentrations from water quality surrogates using machine learning methods[J]. Water Research, 172:115490.

    GERMER S, NEILL C, KRUSCHE A V, et al, 2010. Influence of land-use change on near-surface hydrological processes: undisturbed forest to pasture[J]. Journal of Hydrology, 380(3/4):473-480.

    HOBBIE S E, FINLAY J C, JANKE B D, et al, 2017. Contrasting nitrogen and phosphorus budgets in urban watersheds and implications for managing urban water pollution[J]. Proceedings of the National Academy of Sciences, 114(16):4177-4182.

    JIANG Y, XIE Z, ZHANG H, et al, 2017. Effects of land use types on dissolved trace metal concentrations in the Le’an River Basin, China[J]. Environmental monitoring and assessment, 189: 1-19.

    LI C, ZHANG Y, KHAREL G, et al, 2018. Impact of climate variability and landscape patterns on water budget and nutrient loads in a peri-urban watershed: a coupled analysis using process-based hydrological model and landscape indices[J]. Environmental Management, 61(6):954-967.

    LIU R, WANG J, SHI J, et al, 2014. Runoff characteristics and nutrient loss mechanism from plain farmland under simulated rainfall conditions[J]. Science of the Total Environment, 468/469:1069-1077.

    MELLO K, TANIWAKI R H, PAULA F R, et al, 2020. Multiscale land use impacts on water quality: assessment, planning, and future perspectives in Brazil[J]. Journal of Environmental Management, 270:110879.

    MENESES B M, REIS R, VALE M J, et al, 2015. Land use and land cover changes in Zêzere watershed (Portugal)-water quality implications[J]. Science of the Total Environment, 527:439-447.

    MISHRA A, ALNAHIT A, CAMPBELL B, 2021. Impact of land uses, drought, flood, wildfire, and cascading events on water quality and microbial communities: a review and analysis[J]. Journal of Hydrology, 596:125707.

    OUYANG W, SKIDMORE A K, TOXOPEUS A G, et al, 2010. Long-term vegetation landscape pattern with non-point source nutrient pollution in upper stream of Yellow River basin[J]. Journal of Hydrology,389(3/4):373-380.

    SHIGAKI F, SHARPLEY A, PROCHNOW L I, 2007. Rainfall intensity and phosphorus source effects on phosphorus transport in surface runoff from soil trays[J]. Science of the Total Environment, 373(1):334-343.

    TU J, 2011. Spatially varying relationships between land use and water quality across an urbanization gradient explored by geographically weighted regression[J]. Applied Geography, 31(1):376-392.

    WEI W, GAO Y N, HUANG J C, et al, 2020. Exploring the effect of basin land degradation on lake and reservoir water quality in China[J]. Journal of Cleaner Production, 268:122249.

    WU J, SHEN W, SUN W, et al, 2002. Empirical patterns of the effects of changing scale on landscape metrics[J]. Landsc. Ecol, 17(8): 761-782.

    XU G, REN X, YANG Z, et al, 2019. Influence of landscape structures on water quality at multiple temporal and spatial scales: a case study of Wujiang River Watershed in Guizhou[J]. Water, 11(1):159.

    ZHANG J, LI S Y, JIANG C S, 2020. Effects of land use on water quality in a River Basin (Daning) of the Three Gorges Reservoir Area, China: Watershed versus riparian zone[J]. Ecological Indicators, 113:106226.

    ZHANG X, LUO Y, GOH K S, 2018. Modelingspray drift and runoff-related inputs of pesticides to receiving water[J]. Environmental Pollution, 234:48-58.

    ZHAO B, SUN H, YAN D, et al, 2021. Quantifying changes and drivers of runoff in the Kaidu River Basin associated with plausible climate scenarios[J]. Journal of Hydrology: Regional Studies, 38:100968.

    Impact of Land Use and Landscape Patterns on Water Quality in Kaidu River

    Abstract:Understanding the influence of land use and landscape on river water quality in arid regions is crucial for maintaining ecological security. In this study, Kaidu River was selected for a case study, and we investigated the impact of land use and landscape on water quality in the river using BIOENV and redundancy analysis (RDA), aiming to provide a theoretical reference for water quality management in rivers in arid regions. The study was based on a water quality investigation at 15 sampling sites along the upper, middle and lower reaches of Kaidu River in 2021 during the wet and dry seasons. Based on remote sensing images of the Kaidu River basin in 2020, land use data of circular buffer regions around the sampling points at five scales (200 m, 500 m, 1 000 m, 2 000 m, 5 000 m) were generated as water quality response units, and the land use structure was characterized by the proportion of different land use types in the watershed. The landscape pattern indices at each of the five spatial scales were extracted to characterize the landscape pattern. Results show that: (1) Water quality of Kaidu River was clean or relatively clean at each sampling site, and water quality in the dry season was slightly better than that in the wet season. (2) Land use in the buffer zones at 2 000 m and 5 000 m scales explained water quality well during both the wet and seasons, with explanation rates of 63.81% and 64.91%, respectively. The landscape within the 500 m and 200 m circular buffer zones had the strongest explanatory power for water quality, with explanation rates of 73.88% and 63.33%, respectively. (3) During the wet season, the contribution rates of cropland (57.9%) and the largest patch index (LPI: 48.8%) for water quality was highest, while the land use type and landscape index with the highest contribution rates to water quality in dry season was grassland (68%) and landscape shape index (LSI: 21.6%). (4) During the wet season, the combination of cropland, woodland and water had the highest correlation coefficients (0.424 5) with water quality, and the highest correlation coefficient in dry season was 0.512 9 between grassland and water quality. The combination of landscape indices (NP, LPI, and ENN_MN) had the highest correlation coefficient (0.076 5) with water quality in wet season, and the highest correlation coefficient in dry season was 0.216 3 between the combination of landscape indices (NP, PD, and LPI) and water quality. (5) The explanatory power of the interaction between land use structure and landscape pattern on water quality was 18% within the 2 000 m circular buffer zone during the wet season, higher than the 11% within the 200 m circular buffer zone during the dry season. These results show that both land use structure and landscape pattern have large impacts on water quality in Kaidu River, and their interaction plays an even larger role.

    Key words:land use; landscape pattern; water quality; Kaidu River

    猜你喜歡
    開都河景觀格局土地利用
    開都河流域新疆裸重唇魚胚胎發(fā)育觀察
    基于逐步聚類分析法的開都河徑流預(yù)測模型
    人民黃河(2019年2期)2019-09-10 07:22:44
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    資源型城市轉(zhuǎn)型中景規(guī)格局的優(yōu)化
    開都河諸小河流水資源管理存在的問題及建議
    濱??h土地利用挖潛方向在哪里
    大都市邊緣區(qū)土地利用規(guī)劃對景觀格局的影響
    氣候變化和人類活動對開都河徑流的影響
    近十年我國城市景觀生態(tài)學(xué)研究新進(jìn)展
    綠色科技(2014年4期)2015-01-06 20:15:01
    鴨綠江口濱海濕地保護(hù)區(qū)景觀變化分析
    99热6这里只有精品| 老熟妇仑乱视频hdxx| 在线视频色国产色| 丝袜美腿诱惑在线| 国内精品久久久久久久电影| 91九色精品人成在线观看| 免费在线观看影片大全网站| 波多野结衣高清作品| 最近最新免费中文字幕在线| 99热这里只有精品一区 | 久久久久九九精品影院| 正在播放国产对白刺激| 又紧又爽又黄一区二区| 18禁国产床啪视频网站| a在线观看视频网站| 国产亚洲欧美精品永久| 精品福利观看| 久久精品aⅴ一区二区三区四区| 欧美日韩乱码在线| 国产亚洲欧美98| 中文字幕高清在线视频| 此物有八面人人有两片| 哪里可以看免费的av片| 三级毛片av免费| 脱女人内裤的视频| 日本免费一区二区三区高清不卡| 天堂动漫精品| 最近在线观看免费完整版| 国产av一区在线观看免费| 又黄又爽又免费观看的视频| 午夜福利一区二区在线看| 国产精品久久久久久人妻精品电影| 亚洲成国产人片在线观看| 50天的宝宝边吃奶边哭怎么回事| 不卡一级毛片| 久久久久久大精品| 国产精品亚洲av一区麻豆| 中文资源天堂在线| 99re在线观看精品视频| 亚洲国产精品成人综合色| 欧美一级毛片孕妇| 日本熟妇午夜| 国产免费av片在线观看野外av| 一级片免费观看大全| 老司机福利观看| 免费在线观看视频国产中文字幕亚洲| 男女下面进入的视频免费午夜 | 麻豆一二三区av精品| 日韩欧美 国产精品| 波多野结衣高清作品| 制服诱惑二区| netflix在线观看网站| 欧美+亚洲+日韩+国产| 国产亚洲欧美在线一区二区| 亚洲免费av在线视频| 亚洲av成人一区二区三| 国产又爽黄色视频| av超薄肉色丝袜交足视频| 亚洲成a人片在线一区二区| 精品高清国产在线一区| 12—13女人毛片做爰片一| 十八禁网站免费在线| tocl精华| 9191精品国产免费久久| 国产不卡一卡二| 亚洲aⅴ乱码一区二区在线播放 | 麻豆久久精品国产亚洲av| 国产精品二区激情视频| 人人妻人人澡人人看| 99国产精品99久久久久| 精华霜和精华液先用哪个| 亚洲欧美一区二区三区黑人| 国产精品永久免费网站| 色哟哟哟哟哟哟| 久久精品影院6| 国产真实乱freesex| 级片在线观看| 精品高清国产在线一区| 亚洲欧美精品综合一区二区三区| 深夜精品福利| 久99久视频精品免费| 成人精品一区二区免费| 一级a爱片免费观看的视频| 两个人视频免费观看高清| 精品日产1卡2卡| 美女国产高潮福利片在线看| 国产精品免费一区二区三区在线| 可以在线观看的亚洲视频| 免费一级毛片在线播放高清视频| 色在线成人网| 欧美成人一区二区免费高清观看 | 色综合欧美亚洲国产小说| 亚洲国产精品合色在线| 国产91精品成人一区二区三区| 给我免费播放毛片高清在线观看| 一级毛片高清免费大全| www日本黄色视频网| 大型av网站在线播放| 色尼玛亚洲综合影院| 一级黄色大片毛片| 亚洲精品粉嫩美女一区| 欧美色欧美亚洲另类二区| 国产伦在线观看视频一区| 波多野结衣av一区二区av| av视频在线观看入口| 国产精品久久电影中文字幕| 每晚都被弄得嗷嗷叫到高潮| 久久久久国内视频| 黄色毛片三级朝国网站| 亚洲成人精品中文字幕电影| 又大又爽又粗| 99国产精品99久久久久| 亚洲国产看品久久| 男人操女人黄网站| 成人永久免费在线观看视频| 丰满的人妻完整版| 国产精品美女特级片免费视频播放器 | 女性生殖器流出的白浆| 精品午夜福利视频在线观看一区| 一区二区三区激情视频| 国产精品日韩av在线免费观看| 日本撒尿小便嘘嘘汇集6| 无限看片的www在线观看| 午夜a级毛片| 高清毛片免费观看视频网站| 99久久99久久久精品蜜桃| 成人三级黄色视频| 一本久久中文字幕| 人人妻人人看人人澡| 久久久国产欧美日韩av| 亚洲五月婷婷丁香| 久久 成人 亚洲| 最好的美女福利视频网| 视频在线观看一区二区三区| 淫妇啪啪啪对白视频| 亚洲专区中文字幕在线| 国产伦人伦偷精品视频| 最近最新中文字幕大全电影3 | 午夜福利视频1000在线观看| 美国免费a级毛片| 一区福利在线观看| 欧美日韩福利视频一区二区| 激情在线观看视频在线高清| 十八禁网站免费在线| 精品一区二区三区四区五区乱码| 亚洲专区国产一区二区| 久久久久久大精品| 级片在线观看| 91成人精品电影| 人成视频在线观看免费观看| 亚洲av成人一区二区三| 婷婷六月久久综合丁香| 一进一出抽搐动态| 激情在线观看视频在线高清| 欧美性猛交黑人性爽| 国产v大片淫在线免费观看| 一二三四在线观看免费中文在| 免费高清视频大片| 成年女人毛片免费观看观看9| 又紧又爽又黄一区二区| 亚洲九九香蕉| 国产精品一区二区免费欧美| 国产在线观看jvid| 久9热在线精品视频| 搡老熟女国产l中国老女人| 久久亚洲精品不卡| 美女 人体艺术 gogo| 99国产精品一区二区三区| 久久天堂一区二区三区四区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品啪啪一区二区三区| 久久久水蜜桃国产精品网| 国产亚洲欧美98| 麻豆成人av在线观看| 女性生殖器流出的白浆| 久久午夜综合久久蜜桃| 嫩草影视91久久| 国产成人精品无人区| 我的亚洲天堂| 成人三级黄色视频| 18禁黄网站禁片午夜丰满| 不卡av一区二区三区| 亚洲,欧美精品.| 日本撒尿小便嘘嘘汇集6| 国产黄片美女视频| 极品教师在线免费播放| 亚洲欧美日韩高清在线视频| 国产精品久久久av美女十八| 久久久久九九精品影院| 超碰成人久久| videosex国产| 中文字幕最新亚洲高清| 亚洲国产毛片av蜜桃av| 国产国语露脸激情在线看| 亚洲av五月六月丁香网| 日韩三级视频一区二区三区| 欧美成人免费av一区二区三区| 亚洲精品国产区一区二| а√天堂www在线а√下载| 精品久久久久久久久久免费视频| 18禁黄网站禁片免费观看直播| 91国产中文字幕| 国产激情偷乱视频一区二区| 午夜免费观看网址| 午夜福利成人在线免费观看| 亚洲精品国产区一区二| 最近最新免费中文字幕在线| 国产免费男女视频| 欧美黑人欧美精品刺激| 国产激情久久老熟女| 一进一出抽搐gif免费好疼| 99国产精品一区二区蜜桃av| 精品卡一卡二卡四卡免费| 国产免费av片在线观看野外av| 国产麻豆成人av免费视频| 亚洲午夜精品一区,二区,三区| 亚洲精品国产精品久久久不卡| 色综合欧美亚洲国产小说| 12—13女人毛片做爰片一| 大型av网站在线播放| 免费一级毛片在线播放高清视频| 日韩欧美三级三区| 法律面前人人平等表现在哪些方面| 久久久久久九九精品二区国产 | 日本在线视频免费播放| 亚洲欧美日韩无卡精品| 国产精品美女特级片免费视频播放器 | 成人亚洲精品av一区二区| 日本免费一区二区三区高清不卡| 午夜福利视频1000在线观看| av免费在线观看网站| 露出奶头的视频| 国产精品久久电影中文字幕| 成人三级黄色视频| 日韩一卡2卡3卡4卡2021年| 夜夜看夜夜爽夜夜摸| 日本免费a在线| 淫妇啪啪啪对白视频| 午夜福利欧美成人| 国产极品粉嫩免费观看在线| 精品国产超薄肉色丝袜足j| 欧美日韩乱码在线| 欧美最黄视频在线播放免费| 国产精品免费一区二区三区在线| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利18| 久久九九热精品免费| 亚洲自偷自拍图片 自拍| 国产精品一区二区三区四区久久 | 女人高潮潮喷娇喘18禁视频| 成年人黄色毛片网站| 欧美激情久久久久久爽电影| 日韩欧美免费精品| 一级毛片高清免费大全| 在线天堂中文资源库| 欧洲精品卡2卡3卡4卡5卡区| 99re在线观看精品视频| 久久久精品国产亚洲av高清涩受| 一进一出抽搐gif免费好疼| 午夜福利在线观看吧| 国产亚洲精品一区二区www| 欧美黑人欧美精品刺激| 日韩欧美在线二视频| av天堂在线播放| 亚洲色图 男人天堂 中文字幕| 白带黄色成豆腐渣| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 亚洲熟妇熟女久久| 级片在线观看| 美女午夜性视频免费| 久久久精品国产亚洲av高清涩受| 在线观看www视频免费| 国产高清有码在线观看视频 | 亚洲国产精品合色在线| 国产黄色小视频在线观看| 成人欧美大片| 亚洲精品美女久久久久99蜜臀| 国产精品 欧美亚洲| 无遮挡黄片免费观看| 在线观看免费午夜福利视频| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 熟妇人妻久久中文字幕3abv| 嫩草影院精品99| 欧美一级a爱片免费观看看 | 午夜精品在线福利| 亚洲国产高清在线一区二区三 | 韩国av一区二区三区四区| 国产欧美日韩精品亚洲av| 国产高清videossex| 天天躁狠狠躁夜夜躁狠狠躁| 最新在线观看一区二区三区| 波多野结衣高清无吗| 亚洲色图 男人天堂 中文字幕| 亚洲专区国产一区二区| 亚洲精华国产精华精| 国产精品亚洲美女久久久| 美女高潮到喷水免费观看| 黄片播放在线免费| 黄片小视频在线播放| 亚洲国产中文字幕在线视频| av在线播放免费不卡| 日韩一卡2卡3卡4卡2021年| 成人av一区二区三区在线看| 国产熟女午夜一区二区三区| 色哟哟哟哟哟哟| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 久久久久国内视频| 中文亚洲av片在线观看爽| 欧美在线一区亚洲| 白带黄色成豆腐渣| 老鸭窝网址在线观看| 99热6这里只有精品| 久久久久久免费高清国产稀缺| 久久天躁狠狠躁夜夜2o2o| 欧美国产精品va在线观看不卡| 欧美激情极品国产一区二区三区| 亚洲精品av麻豆狂野| 99热只有精品国产| xxx96com| 日韩欧美三级三区| 精品一区二区三区av网在线观看| 一进一出好大好爽视频| www.熟女人妻精品国产| 丰满的人妻完整版| 欧美大码av| 久久99热这里只有精品18| 色综合欧美亚洲国产小说| 色在线成人网| 精品久久蜜臀av无| 美女大奶头视频| 97超级碰碰碰精品色视频在线观看| 久久久久久九九精品二区国产 | 成人18禁高潮啪啪吃奶动态图| 久久精品aⅴ一区二区三区四区| 老司机福利观看| 成人三级黄色视频| 成熟少妇高潮喷水视频| 亚洲精品久久成人aⅴ小说| 老司机福利观看| 国产又爽黄色视频| 精品久久蜜臀av无| 在线观看一区二区三区| 法律面前人人平等表现在哪些方面| 97人妻精品一区二区三区麻豆 | 日韩免费av在线播放| 国产熟女午夜一区二区三区| 日本免费一区二区三区高清不卡| 无人区码免费观看不卡| 韩国精品一区二区三区| e午夜精品久久久久久久| av天堂在线播放| 99热6这里只有精品| 欧美+亚洲+日韩+国产| 日本一区二区免费在线视频| 97碰自拍视频| 操出白浆在线播放| 神马国产精品三级电影在线观看 | 搡老熟女国产l中国老女人| 亚洲人成77777在线视频| 看片在线看免费视频| 久久久水蜜桃国产精品网| 免费观看精品视频网站| 亚洲中文av在线| 亚洲人成网站高清观看| 久久久水蜜桃国产精品网| 日韩大尺度精品在线看网址| 国产一区二区激情短视频| 久久 成人 亚洲| 久久久久久久午夜电影| 亚洲av片天天在线观看| 久久中文字幕人妻熟女| 757午夜福利合集在线观看| 满18在线观看网站| 一级a爱视频在线免费观看| 免费人成视频x8x8入口观看| 女同久久另类99精品国产91| 青草久久国产| 丝袜在线中文字幕| 午夜激情福利司机影院| 波多野结衣高清无吗| 在线看三级毛片| 久久人人精品亚洲av| 少妇的丰满在线观看| 搡老岳熟女国产| 成熟少妇高潮喷水视频| 国产精品久久久久久精品电影 | 久久久国产精品麻豆| 亚洲国产欧美网| 亚洲性夜色夜夜综合| 一本综合久久免费| 亚洲黑人精品在线| 一区二区三区激情视频| av福利片在线| 欧美+亚洲+日韩+国产| 老司机午夜十八禁免费视频| 波多野结衣av一区二区av| 国产v大片淫在线免费观看| 国产又色又爽无遮挡免费看| 天堂影院成人在线观看| 国产精品一区二区三区四区久久 | 宅男免费午夜| 啦啦啦韩国在线观看视频| 女生性感内裤真人,穿戴方法视频| 一级毛片高清免费大全| 午夜亚洲福利在线播放| 国产精品1区2区在线观看.| 18禁裸乳无遮挡免费网站照片 | 啪啪无遮挡十八禁网站| 亚洲三区欧美一区| 99在线视频只有这里精品首页| 亚洲全国av大片| 女警被强在线播放| 久久人人精品亚洲av| 91麻豆av在线| 国产私拍福利视频在线观看| 啦啦啦观看免费观看视频高清| 久久香蕉国产精品| av电影中文网址| 久久精品国产亚洲av高清一级| 欧美日韩精品网址| 两性夫妻黄色片| 999久久久国产精品视频| 中文字幕精品亚洲无线码一区 | 三级毛片av免费| 欧美日本亚洲视频在线播放| 特大巨黑吊av在线直播 | 可以在线观看毛片的网站| 精品久久久久久久毛片微露脸| 一本大道久久a久久精品| 狠狠狠狠99中文字幕| av免费在线观看网站| 给我免费播放毛片高清在线观看| 嫩草影视91久久| 国产成人欧美在线观看| 精品国产美女av久久久久小说| 欧美久久黑人一区二区| 国产亚洲精品一区二区www| 香蕉丝袜av| 女性被躁到高潮视频| 欧美成人一区二区免费高清观看 | 一级a爱视频在线免费观看| 国产亚洲av嫩草精品影院| 日本免费一区二区三区高清不卡| 久久久国产成人精品二区| 午夜a级毛片| 精品第一国产精品| 成人免费观看视频高清| 欧美精品啪啪一区二区三区| 在线观看免费日韩欧美大片| 亚洲全国av大片| 亚洲成国产人片在线观看| 国产高清激情床上av| 91大片在线观看| 欧美激情极品国产一区二区三区| 欧美绝顶高潮抽搐喷水| 国产精品久久视频播放| 日韩欧美 国产精品| 最新在线观看一区二区三区| 国产免费av片在线观看野外av| 丰满人妻熟妇乱又伦精品不卡| 国产成年人精品一区二区| 中文字幕高清在线视频| 欧美在线一区亚洲| 999精品在线视频| aaaaa片日本免费| 啦啦啦免费观看视频1| 亚洲天堂国产精品一区在线| 欧美性猛交黑人性爽| 天天一区二区日本电影三级| 久久性视频一级片| av中文乱码字幕在线| 熟女少妇亚洲综合色aaa.| 亚洲精品色激情综合| 久久久久久久久中文| 欧美zozozo另类| 亚洲中文字幕日韩| 日韩一卡2卡3卡4卡2021年| 男女下面进入的视频免费午夜 | 国产免费av片在线观看野外av| 看免费av毛片| 18禁美女被吸乳视频| 精品不卡国产一区二区三区| 18禁黄网站禁片午夜丰满| 国产又爽黄色视频| 精品一区二区三区av网在线观看| 丝袜美腿诱惑在线| 国产亚洲精品久久久久久毛片| 中文字幕人妻熟女乱码| www日本黄色视频网| 91大片在线观看| 露出奶头的视频| 欧美日韩精品网址| 欧美日韩福利视频一区二区| 中文字幕人成人乱码亚洲影| 精品国产乱子伦一区二区三区| 日韩视频一区二区在线观看| 精品福利观看| 久久九九热精品免费| 免费看a级黄色片| 午夜成年电影在线免费观看| 久久久水蜜桃国产精品网| 两性午夜刺激爽爽歪歪视频在线观看 | 精品熟女少妇八av免费久了| 精华霜和精华液先用哪个| 亚洲人成网站高清观看| 日韩欧美免费精品| 欧美三级亚洲精品| 国产亚洲欧美98| 午夜免费观看网址| 亚洲国产日韩欧美精品在线观看 | 欧美日韩黄片免| 欧美黄色淫秽网站| 久久中文看片网| 欧美乱码精品一区二区三区| 国产成人影院久久av| 亚洲av成人不卡在线观看播放网| 国产精品久久电影中文字幕| 国产亚洲欧美在线一区二区| 午夜激情av网站| 欧美日韩福利视频一区二区| 一区福利在线观看| 国产1区2区3区精品| 男女之事视频高清在线观看| 国产午夜福利久久久久久| 成年免费大片在线观看| 在线免费观看的www视频| 美国免费a级毛片| 嫁个100分男人电影在线观看| 女人高潮潮喷娇喘18禁视频| 最近最新中文字幕大全电影3 | 久久 成人 亚洲| 免费人成视频x8x8入口观看| 91字幕亚洲| 日本在线视频免费播放| 国产一区二区三区在线臀色熟女| 一二三四在线观看免费中文在| 久久这里只有精品19| 丰满的人妻完整版| 韩国av一区二区三区四区| 18禁裸乳无遮挡免费网站照片 | 在线观看日韩欧美| 亚洲片人在线观看| 久久国产亚洲av麻豆专区| 母亲3免费完整高清在线观看| 国产精品综合久久久久久久免费| 精品久久久久久成人av| 欧美日韩中文字幕国产精品一区二区三区| 国产精品1区2区在线观看.| 久久这里只有精品19| 久久精品国产清高在天天线| 好看av亚洲va欧美ⅴa在| 久久天躁狠狠躁夜夜2o2o| 欧美又色又爽又黄视频| 亚洲天堂国产精品一区在线| 香蕉丝袜av| 身体一侧抽搐| 99热这里只有精品一区 | 村上凉子中文字幕在线| 50天的宝宝边吃奶边哭怎么回事| 国产成人精品久久二区二区91| 亚洲国产欧美日韩在线播放| www.精华液| 一个人免费在线观看的高清视频| 亚洲三区欧美一区| 国产v大片淫在线免费观看| 变态另类丝袜制服| 激情在线观看视频在线高清| 久久这里只有精品19| 亚洲成人国产一区在线观看| 超碰成人久久| 久99久视频精品免费| 日韩欧美一区二区三区在线观看| netflix在线观看网站| 十八禁人妻一区二区| 亚洲自拍偷在线| 啦啦啦观看免费观看视频高清| 99久久无色码亚洲精品果冻| 国产精品影院久久| 亚洲人成网站高清观看| 99国产精品一区二区蜜桃av| 国产精品久久久久久亚洲av鲁大| 欧美久久黑人一区二区| 欧美zozozo另类| 国产精品影院久久| 91老司机精品| 99国产精品一区二区蜜桃av| 国产精品影院久久| 国产精品一区二区免费欧美| 99国产精品一区二区蜜桃av| 亚洲精品在线观看二区| 黄片播放在线免费| 亚洲黑人精品在线| 最好的美女福利视频网| 欧美最黄视频在线播放免费| 18美女黄网站色大片免费观看| 亚洲自拍偷在线| 成人欧美大片| 18禁裸乳无遮挡免费网站照片 | 欧美色欧美亚洲另类二区| 日韩欧美三级三区| 曰老女人黄片| 两人在一起打扑克的视频| 一二三四社区在线视频社区8| 国产精品一区二区精品视频观看| 亚洲成a人片在线一区二区| 十八禁网站免费在线| 亚洲欧美精品综合一区二区三区| 女人被狂操c到高潮| 欧美激情极品国产一区二区三区| 一边摸一边做爽爽视频免费| 老司机午夜福利在线观看视频| 国产亚洲欧美在线一区二区|