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

    快速城市化地區(qū)昆明市“城—郊—鄉(xiāng)”梯度土壤表層碳氮磷化學(xué)計量特征

    2020-06-08 10:56:56何云玲張淑潔鄧福英榮麗張林艷
    生態(tài)科學(xué) 2020年3期
    關(guān)鍵詞:研究

    何云玲, 張淑潔, 鄧福英, 榮麗, 張林艷

    快速城市化地區(qū)昆明市“城—郊—鄉(xiāng)”梯度土壤表層碳氮磷化學(xué)計量特征

    何云玲1,*, 張淑潔1, 鄧福英1, 榮麗2, 張林艷3

    1. 云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院, 昆明 650091 2. 云南大學(xué)國際河流與生態(tài)安全研究院, 昆明 650091 3. 云南大學(xué)建筑與規(guī)劃學(xué)院, 昆明 650091

    基于航空遙感影像, 利用地理信息系統(tǒng)(ArcGIS 10.2)和景觀格局分析(Fragstats 4.2)軟件, 通過主成分分析(PCA)和多元線性回歸分析對研究區(qū)城市化水平進行量化, 選取覆蓋研究區(qū)典型的“城—郊—鄉(xiāng)”樣帶, 進行表層土壤取樣, 測定總有機碳(C)、全氮(N)、全磷(P)的含量并分析其空間變化特征。結(jié)果表明: 快速的城市化進程導(dǎo)致了土壤C、N、P化學(xué)計量特征的空間分異, “城—郊—鄉(xiāng)”梯度上變異系數(shù)都較大, 尤其城區(qū)各元素的變異系數(shù)分別為84.08%、76.95%和81.34%, 屬于強變異, 表明土壤表層C、N、P含量空間分布極不均衡。研究區(qū)土壤C、N、P含量均表現(xiàn)為城區(qū)顯著高于郊區(qū)的特征, 與郊區(qū)土壤相比, 城區(qū)土壤C含量較高, 說明人為活動造成城區(qū)土壤明顯的有機質(zhì)富集現(xiàn)象。城區(qū)土壤P含量高于郊區(qū)和鄉(xiāng)村, 城市化過程加劇了城區(qū)土壤P的富集, 也使該地區(qū)水體(主要是滇池)面臨富營養(yǎng)化風(fēng)險。分析城市化過程中強烈的人類活動干擾對土壤C、N、P含量及其化學(xué)計量特征, 對于提高城市土壤管理水平、改善生態(tài)環(huán)境、促進城市生態(tài)系統(tǒng)健康持續(xù)發(fā)展具有重要意義。今后昆明市在進行城市綠化時可有針對性地選擇部分豆科固氮植物加強生物固氮, 目的是提高土壤N含量的同時促進城區(qū)土壤豐富的P被吸收利用, 在一定程度上預(yù)防以及減輕P對水體的危害。

    快速城市化; “城—郊—鄉(xiāng)”; C:N:P; 空間變異; 昆明中部

    0 前言

    土壤有機碳庫的細微變化都將會影響全球氣候變化以及全球碳平衡[1], 而且土壤碳(C)、氮(N)、磷(P)之比是反映土壤環(huán)境質(zhì)量的靈敏指標(biāo)之一, 常作為診斷或預(yù)測養(yǎng)分限制或飽和的指標(biāo)[2]。森林、農(nóng)田、草原等土地利用方式改變對土壤碳庫及通量影響, 已得到較廣泛的關(guān)注, 而由人類活動主導(dǎo)的劇烈環(huán)境變化所驅(qū)動的土壤環(huán)境演化及其生態(tài)響應(yīng)仍是一個尚未系統(tǒng)回答的科學(xué)問題[3]。伴隨著城市化進程的加快和城市規(guī)模的不斷擴張, 越來越多的自然用地和農(nóng)業(yè)用地轉(zhuǎn)變?yōu)槌鞘薪ㄔO(shè)用地, 改變了城市土壤的理化性質(zhì), 必將影響土壤C、N、P含量的變化, 這對全球氣候變化及陸地生態(tài)系統(tǒng)將產(chǎn)生深刻的影響[4]。

    近年來國內(nèi)外學(xué)者對城市區(qū)域土壤C、N、P的研究熱度不斷提高, Rabbi等分析了澳大利亞東部不同類型土壤有機碳含量的差異及成因[5], Lidia等分析了波蘭城市化對土壤有機碳的影響[6]。在人為和自然因素的交互影響下, 城市綠地C、N、P時空異質(zhì)性顯著[7–8]。從目前的研究結(jié)果來看, 不同時間、不同城市、城市內(nèi)部不同功能區(qū), 綠地土壤C、N、P的分布均存在變異, 例如許乃政等對上海市土壤有機碳含量研究發(fā)現(xiàn)隨著城市用地年限延長, 土壤有機碳密度逐漸增大[9]; 呂麗平等研究表明海南的土壤有機碳含量隨城市化進程總體呈下降趨勢[10]; Chen等對南昌市“城—郊—鄉(xiāng)”梯度上土壤P素研究結(jié)果發(fā)現(xiàn)全P含量城區(qū)明顯高于郊區(qū)[11]。導(dǎo)致研究結(jié)果不一致的其中一個原因是對城、郊、鄉(xiāng)的地域選擇多是基于行政區(qū)劃進行, 但是“城—郊—鄉(xiāng)”梯度的形成應(yīng)該是多種因素相互作用的結(jié)果, 對其準(zhǔn)確描述是開展城市化強度差異對生態(tài)系統(tǒng)產(chǎn)生影響的重要前提[12]。另外, 土壤C、N、P化學(xué)計量特征伴隨城市化進程在“城—郊—鄉(xiāng)”梯度下的研究還少有報道[13]。

    昆明市作為云南省省會, 中國面向東南亞、南亞開放的門戶城市, 自2000年以來隨著舊城改造不斷推進以及新區(qū)建設(shè), 城市快速發(fā)展, 建設(shè)用地擴張速度逐年增加, 城市土地利用和景觀格局以前所未有的速度發(fā)生著劇烈的變化, 城市環(huán)境發(fā)生了很大的改變。本研究在遙感影像資料和GIS技術(shù)支持下, 選擇昆明市包含城—郊—鄉(xiāng)梯度差異的典型樣帶為研究區(qū)域, 由中心城區(qū)到郊區(qū)到鄉(xiāng)村, 分析城市化過程中強烈的人類活動干擾對樣帶土壤C、N、P含量及其化學(xué)計量特征, 對于提高城市土壤管理水平、改善生態(tài)環(huán)境、促進城市生態(tài)系統(tǒng)健康持續(xù)發(fā)展具有重要意義。

    1 材料與方法

    昆明市位于中國西南云貴高原中部, 東經(jīng)102°10′—103°41′, 北緯24°24′—26°33′之間, 屬北亞熱帶低緯高原山地季風(fēng)氣候, 氣候溫和, 四季如春, 年平均氣溫16.5 ℃, 年降雨量1450 mm, 5—10月為雨季, 降水量占全年的85%左右; 11月至次年4月為干季, 降水量僅占全年的15%左右; 市中心海拔1891 m, 土壤主要為紅壤和磚紅壤。2017年末, 全市常住人口678.3萬人, 建成區(qū)面積達到420.5 km2。由于其特殊的地理位置, 是中國面向東南亞、南亞等的前沿和門戶, 伴隨著昆明市人口不斷增多與城市建設(shè)用地面積的急劇增長, 城市化進程進一步加劇??紤]到其山間盆地的地貌類型, 本文以昆明市中部主要行政區(qū)和建成區(qū)所在地的西山、五華、盤龍和官渡這四個區(qū)內(nèi)具有連續(xù)的城鄉(xiāng)交錯帶作為主體研究區(qū), 這個區(qū)域也是昆明市快速城市化擴張的代表區(qū)域。

    2 研究方法

    2.1 城-郊-鄉(xiāng)梯度的劃分

    “城—郊—鄉(xiāng)”梯度的劃分是研究人類活動對生態(tài)環(huán)境影響的基礎(chǔ), 以往研究表明城市化物理變量與景觀變量的結(jié)合可在一定程度上較好地完成“城—郊—鄉(xiāng)”梯度的量化[14–15]。本文利用精度2.5 m的衛(wèi)星影像數(shù)據(jù)與Landsat8遙感影像數(shù)據(jù), 使用ENVI軟件獲取研究區(qū)土地利用信息, 并結(jié)合ArcGIS10.2與Fragstats4.2軟件以1 km×1 km為網(wǎng)格單位提取不透水下墊面比例、道路網(wǎng)密度、景觀形狀指數(shù)與香農(nóng)多樣性指數(shù)4個城市化指標(biāo), 并在SPSS22.0軟件里采用主成分分析(PCA)與多元線性回歸分析構(gòu)建“城—郊—鄉(xiāng)”梯度模型對研究區(qū)的城市化水平指數(shù)(UI)進行量化。其中, 不透水下墊面比例和道路網(wǎng)密度是城市化進程中最顯著和直觀的特征, 可反映城市化發(fā)展方向的變化趨勢; 景觀形狀指數(shù)和香農(nóng)多樣性指數(shù)分別是景觀形狀的復(fù)雜程度和異質(zhì)性, 數(shù)值越大, 表明斑塊越復(fù)雜, 土地利用的破碎化程度越高。根據(jù)研究區(qū)城市化指數(shù)UI值, 將其導(dǎo)入ArcGIS網(wǎng)格面與矢量化后的土地分類圖的疊加圖屬性表中, 并與昆明市衛(wèi)星影像數(shù)據(jù)進行對比分析, 最后將-5.804

    根據(jù)城—郊—鄉(xiāng)劃分的網(wǎng)格單位, 在研究區(qū)內(nèi)選取面積相等的連續(xù)格網(wǎng)(東西方向1 km, 南北方向1 km)作為梯度帶, 自西向東, 共設(shè)置44個連續(xù)梯度網(wǎng)格(以涵蓋城鄉(xiāng)交錯土地景觀空間格局連續(xù)區(qū)域為原則), 分別記做序號1, 2, ...... , 44, 如圖1所示。將樣帶上各城市化指標(biāo)進行分析(圖2), 可見隨著從鄉(xiāng)村到郊區(qū)到城區(qū)梯度的變化, 不透水下墊面比例、道路網(wǎng)密度、景觀形狀指數(shù)和香農(nóng)多樣性指數(shù)各指標(biāo)數(shù)值逐漸增大; 鄉(xiāng)村的不透水下墊面比例與道路網(wǎng)密度均明顯低于城區(qū)與郊區(qū), 且城區(qū)最高(其中不透水下墊面比例高達90%以上); 城區(qū)的景觀形狀指數(shù)與香農(nóng)多樣性指數(shù)最大, 表明城區(qū)的景觀異質(zhì)性最大, 且景觀破碎化程度明顯高于郊區(qū)與鄉(xiāng)村。經(jīng)過實地驗證和遙感影像比對, 該連續(xù)梯度網(wǎng)格帶符合城—郊—鄉(xiāng)梯度的基本特征。

    2.2 土壤采樣及分析方法

    土壤樣品的采集于2018年5月初開始進行, 針對上述研究區(qū)劃分出來的自西向東樣帶內(nèi)的44個1 km×1 km的網(wǎng)格, 在每個網(wǎng)格內(nèi)采用對角線采樣法采集表層(0—30 cm)各5個土壤樣品, 共采集220個土樣, 土壤采樣點綠地類型見表1所示。采樣時使用GPS定位采樣點, 記錄經(jīng)緯度坐標(biāo)。使用工兵鏟在地表進行深度超過30 cm的采樣坑樣品采集, 沿坑壁垂直向下均勻采集0—30 cm的土壤樣品約1 kg, 去除草根、石塊。將每個網(wǎng)格內(nèi)的5個土壤樣品混合, 按照四分法取混合土樣1—2 kg, 密封在聚乙烯塑料樣品袋內(nèi)并標(biāo)記。同時在每個采樣點現(xiàn)場使用100 cm3的環(huán)刀采集環(huán)刀樣品用以測定土壤容重, 每個采樣點進行3次重復(fù)采樣, 共采集樣品660個; 并在每個采樣點按照上層(0—10 cm)、中層(0—20 cm)、下層(20—30 cm)分別采集10—20 g左右土壤, 放入鋁盒中蓋緊, 帶回實驗室立即稱重并進行烘干用以計算土壤含水量。將采集的土壤樣品風(fēng)干后, 經(jīng)木棒研磨去除動植物殘體和碎石等異物, 過2 mm尼龍篩混勻, 裝入聚乙烯塑料樣品袋內(nèi)標(biāo)記、裝箱, 運送至資質(zhì)檢測機構(gòu)對總有機碳含量(C)、全氮(N)、全磷(P)等進行測定。土壤總有機碳采用重鉻酸鉀—外加熱法測定, 全氮采用半微量開氏法測定, 全磷采用氫氟酸—高氯酸酸溶—鉬銻抗比色法測定。樣品分析結(jié)果采用SPSS22.0進行相關(guān)分析、回歸分析、以及方差分析與多重比較。

    圖1 昆明市中部樣帶分布圖

    Figure 1 Distribution of soil sampling transect along urban- suburb-rural gradient in central Kunming

    圖2 沿城—郊—鄉(xiāng)樣帶城市化指標(biāo)值變化特征

    Figure 2 Variations of urbanization index along the urban-suburb-rural gradient belt in central Kunming

    表1 土壤采樣點植被覆蓋類型

    3 結(jié)果與分析

    3.1 碳氮磷含量特征

    昆明市中部調(diào)查樣方內(nèi)土壤表層總有機碳(C)、全氮(N)、全磷(P)含量的頻數(shù)分布表明(圖3), C含量的范圍為1.85—48.58 g·kg–1, 平均值為16.18 g·kg–1, 變異系數(shù)(CV)達到72.25%; N含量的范圍為0.32—3.88 g·kg–1, 平均值為1.37 g·kg–1, 變異系數(shù)69.78%; P含量的范圍為0.29—3.15 g·kg–1, 平均值為0.80 g·kg–1,變異系數(shù)58.70%??梢? 昆明市中部調(diào)查樣帶內(nèi)土壤表層C、N、P含量變異系數(shù)都較大, 屬于強變異, 空間分布極不均衡。

    昆明市中部城—郊—鄉(xiāng)梯度上土壤表層總有機碳(C)、全氮(N)、全磷(P)含量變化(圖4)可見, 不同城市化水平的區(qū)域土壤C、N、P含量變化趨勢大體一致, 整體呈現(xiàn)出鄉(xiāng)村的含量數(shù)值較大, 伴隨城市化水平的增大, C、N、P含量向城區(qū)遞減。研究區(qū)土壤表層C含量在“城—郊—鄉(xiāng)”梯度上呈現(xiàn)的變化特征與李家熙等[17]對蘇州的研究結(jié)果一致: 城市綠地土壤有機碳含量低于農(nóng)業(yè)土壤與山地土壤; 但與Luo等對北京[18]、Liu等[19]對重慶的城市土壤研究結(jié)果不一致: 土壤有機碳含量隨著距離城市中心距離的增加而逐漸減小。

    3.2 碳氮磷化學(xué)計量關(guān)系

    土壤C與N含量相關(guān)系數(shù)高達0.931(<0.01)呈極顯著正相關(guān)關(guān)系(圖5), 土壤C與P含量相關(guān)系數(shù)0.728(<0.01)呈顯著正相關(guān)關(guān)系; 表明土壤N、P含量是影響C含量變化的重要因素之一, 土壤中N、P含量的升高能夠加快植物生長, 致使更多的凋落物在地面聚集, 增加土壤有機質(zhì)的輸入, 這一結(jié)果與馮錦等對新疆[20]、王天樂等對內(nèi)蒙[21]的土壤C、N、P計量特征研究結(jié)果相一致。

    圖3 調(diào)查樣方土壤總有機碳(C)、全氮(N)、全磷(P)的含量頻數(shù)分布圖

    Figure 3 Frequency distribution histograms of topsoil total organic carbon, total nitrogen and total phosphorus in survey plots

    圖4 城—郊—鄉(xiāng)梯度上土壤表層總有機碳(C)、全氮(N)、全磷(P)含量分布圖

    Figure 4 Contents of topsoil total organic carbon, total nitrogen and total phosphorus along the urban-suburb -rural gradient belt

    圖5 土壤表層總有機碳(C)和全氮(N)、全磷(P)含量的相關(guān)分析(直線為線性擬合, Tau為Mann-Kendall秩相關(guān)系數(shù), 帶**表示達到了99%的置信水平)

    Figure 5 The correlation analysis of total organic carbon and total nitrogen, total phosphorus

    昆明市中部城—郊—鄉(xiāng)不同區(qū)域土壤表層總有機碳(C)、全氮(N)、全磷(P)化學(xué)計量的統(tǒng)計特征如圖6所示。研究區(qū)內(nèi)土壤表層C含量大小表現(xiàn)為鄉(xiāng)村>城區(qū)>郊區(qū), 鄉(xiāng)村土壤表層C含量分別為郊區(qū)、城區(qū)的3.06倍、1.38倍; N含量大小也表現(xiàn)為鄉(xiāng)村>城區(qū)>郊區(qū), 鄉(xiāng)村土壤表層N含量分別是郊區(qū)、城區(qū)的2.48倍、1.31倍; P含量大小則表現(xiàn)為城區(qū)>鄉(xiāng)村>郊區(qū), 城區(qū)土壤表層P含量分別是郊區(qū)、鄉(xiāng)村的1.69倍、1.15倍。另外, 變異系數(shù)計算結(jié)果顯示土壤表層C、N和P含量在城區(qū)的變異系數(shù)均為最大, 分別為84.08%、76.95%和81.34%, 屬于強變異, 空間分布極不均衡。再從不同區(qū)域間的差異顯著性檢驗來看, 土壤表層C、N、P含量均表現(xiàn)為城區(qū)顯著高于郊區(qū)的變化特征。

    圖6 城—郊—鄉(xiāng)不同區(qū)域土壤表層總有機碳(C)、全氮(N)、全磷(P)化學(xué)計量(圖中不同大寫字母、小寫字母和數(shù)字分別表示各要素不同區(qū)域間在p<0.05水平顯著)

    Figure 6 Stoichiometry of topsoil total organic carbon, total nitrogen and total phosphorus in different areas

    研究區(qū)內(nèi)土壤表層的C:N最大值(16.62)和最小值(2.94), C:P最大值(44.80)和最小值(2.03), N:P最大值(3.17)和最小值(0.38), 數(shù)值上均低于國內(nèi)平均水平, 空間上也均出現(xiàn)在鄉(xiāng)村區(qū)域; 化學(xué)計量特征平均值大小順序均表現(xiàn)為鄉(xiāng)村>城區(qū)>郊區(qū)。從不同區(qū)域間的差異顯著性檢驗來看, 土壤表層C:N和N:P均表現(xiàn)為郊區(qū)顯著低于其他區(qū)域, 城區(qū)和鄉(xiāng)村之間則無明顯差異; 土壤表層C:P則在城區(qū)、郊區(qū)和鄉(xiāng)村兩兩之間無顯著差異。

    4 討論與結(jié)論

    4.1 討論

    (1) 土壤碳氮磷含量特征

    城市土壤常常遭受各種自然和人為活動導(dǎo)致的物理干擾、化學(xué)污染, 通常具有高度空間異質(zhì)性[22]。不少研究表明, 隨著城市—郊區(qū)—鄉(xiāng)村的空間梯度演替, 表層土壤有機碳含量漸趨減少[23], 被普遍認為是碳源的城市區(qū)域, 其土壤可能具有較大的碳吸存潛力[24]。國內(nèi)外現(xiàn)有研究所用土壤有機碳數(shù)據(jù)采樣深度多集中于表層20 cm(10—40 cm), 昆明市中部城區(qū)土壤表層C含量13.84 g·kg–1, 高于柳云龍等[25]和段保正等[26]分別對上海市(11.91 g·kg–1)和西安市(5.59 g·kg–1)的研究結(jié)果; 低于張小萌等[27]對烏魯木齊(33.05 g·kg–1)的研究結(jié)果。昆明市中部“城—郊—鄉(xiāng)”土壤表層N含量略高于同處于中國西南地區(qū)鄰近的成都東部“城—郊—鄉(xiāng)”土壤調(diào)查結(jié)果, 但是P含量均<1 g·kg–1, 低于成都的調(diào)查結(jié)果[13]。從整體上看, 昆明市中部土壤C、N、P含量均表現(xiàn)為城區(qū)顯著高于郊區(qū)的特征, 與王秋兵等[28]、曾宏達等[29]分別對沈陽市、福州市的研究結(jié)果相同。

    與郊區(qū)土壤相比, 城區(qū)土壤有機碳含量較高, 說明城市土壤表現(xiàn)出明顯的有機質(zhì)富集現(xiàn)象。城區(qū)土壤有機質(zhì)除少量來自本身植物枯枝落葉和雜草, 主要由工業(yè)生產(chǎn)、生活垃圾排放等人類活動產(chǎn)生的有機廢棄物帶入; 而且城區(qū)土壤不以生產(chǎn)為目的, 大部分綠地受到較好的人為養(yǎng)護與管理, 較為堅實, 土壤有機碳分解較慢, 有利于有機碳的積累[9]。土壤N、P含量在城區(qū)顯著高于郊區(qū), 一方面可能與大部分城區(qū)以前多為菜地等農(nóng)田, 長期施肥導(dǎo)致氮素積累, 另外一方面與城市化進程中人為富P垃圾堆填過多, 且P的遷移能力差有關(guān)[30]。

    與鄉(xiāng)村土壤相比, 城區(qū)土壤C、N、P含量的區(qū)域差異性雖然未達到顯著性檢驗, 其中C、N含量總體表現(xiàn)為鄉(xiāng)村大于城區(qū)的特征, 主要原因在于本研究中鄉(xiāng)村土壤采樣點多為自然狀態(tài)下的林地土壤, 人為干擾較小, 郁閉度較高, 表層土壤C、N受植被枯落物歸還等影響而富集[31]。研究結(jié)果與劉偉玲等[32]對深圳市不同土地利用類型下土壤有機碳含量的研究結(jié)果相一致: 林地土壤有機碳含量高于綠地與農(nóng)用地的土壤有機碳含量。

    (2) 土壤碳氮磷計量關(guān)系

    土壤具有高度的結(jié)構(gòu)復(fù)雜性、空間異質(zhì)性和生物多樣性, 但是卻有研究表明, 在大的區(qū)域尺度上土壤中C:N:P具有較為穩(wěn)定的特征, 通常在一定范圍之內(nèi)波動[33]。例如Cleveland和Liptzin研究發(fā)現(xiàn)全世界范圍內(nèi)土壤C、N、P三者含量之間具有顯著的正相關(guān)關(guān)系, 比值約為186:13:1[34]。但是土地利用和管理措施的變化等因素能顯著改變元素的生物地球化學(xué)循環(huán), 導(dǎo)致生態(tài)系統(tǒng)C、N、P的化學(xué)計量特征發(fā)生改變[35]。昆明市中部“城—郊—鄉(xiāng)”不同區(qū)域土壤養(yǎng)分含量變化很大, 使得土壤C、N、P之比存在較大的空間變異性, 而且土壤C:N:P指標(biāo)均遠低于全球以及國內(nèi)平均水平。

    土壤C:N可用于判斷土壤中有機質(zhì)的分解程度及其對土壤肥力的潛在貢獻, 適當(dāng)?shù)腃:N比例有助于微生物發(fā)酵分解; 高的C:N(>25)說明有機質(zhì)正在累積, 其累積速率大于分解速率, C:N在12—16之間說明有機質(zhì)已經(jīng)被微生物很好地分解[36]。Cusack在波多黎各亞熱帶土壤研究中發(fā)現(xiàn)城區(qū)森林土壤C:N低于鄉(xiāng)村, 城市化過程強烈影響土壤N素的初級轉(zhuǎn)化速率[37]。昆明市中部“城—郊—鄉(xiāng)”土壤C:N 在10—12之間, 和一般耕作土壤C:N數(shù)值表現(xiàn)一致[36]; 而且研究區(qū)結(jié)果表明土壤C:N有郊區(qū)顯著低于城區(qū)的趨勢, 可能與城區(qū)的人為活動造成綠地土壤C富集, 而N相對穩(wěn)定有關(guān)。

    土壤C:P是衡量微生物礦化、土壤有機物質(zhì)釋放P或從環(huán)境中吸收固持P素潛力的一個指標(biāo)[38]; 土壤N:P則指示了植物生長過程中土壤營養(yǎng)成分的供應(yīng)情況, 可用來判斷養(yǎng)分限制狀況及哪種養(yǎng)分限制了有機質(zhì)的分解[36]。昆明市“城—郊—鄉(xiāng)”土壤C:P和N:P都相對較小, 尤其是N:P遠低于我國土壤5.2的平均值[39]; 而且總體表現(xiàn)為城區(qū)土壤P含量高于郊區(qū)和鄉(xiāng)村, 城市化過程加劇了城區(qū)土壤P的富集, 較低的N:P與P顯著積累有關(guān), 較低的C:P也是P有效性高的一個指標(biāo), P較多, 也使該地區(qū)水體(主要是滇池)面臨富營養(yǎng)化風(fēng)險。因此, 今后昆明市在進行城市綠化時可有針對性地選擇部分豆科固氮植物加強生物固氮, 目的是提高土壤N含量的同時促進城區(qū)土壤豐富的P被吸收利用, 在一定程度上預(yù)防以及減輕P對水體的危害[13]。

    4.2 結(jié)論

    通過對研究區(qū)土壤表層總有機碳(C)、總氮(N)、總磷(P)含量與計量關(guān)系在“城—郊—鄉(xiāng)”梯度上空間分布特征的研究, 得出以下結(jié)論:

    (1)研究區(qū)“城—郊—鄉(xiāng)”土壤C含量1.85—48.58 g·kg–1, 平均值16.18 g·kg–1, 變異系數(shù)72.25%; N含量0.32—3.88 g·kg–1, 平均值1.37 g·kg–1, 變異系數(shù)69.78%; P含量0.29—3.15 g·kg–1, 平均值0.80 g·kg–1,變異系數(shù)58.70%。土壤C、N、P含量屬于強變異, 尤其是城區(qū), 表明空間分布極不均衡。

    (2) 研究區(qū)“城—郊—鄉(xiāng)”土壤C與N含量呈極顯著正相關(guān)關(guān)系, C與P呈顯著正相關(guān)關(guān)系。表明土壤N、P含量是影響C含量變化的重要因素之一, N、P含量升高能夠加快植物生長, 致使更多的凋落物在地面聚集, 增加土壤有機質(zhì)的輸入。

    (3) 研究區(qū)“城—郊—鄉(xiāng)”土壤C、N含量均表現(xiàn)為城區(qū)顯著高于郊區(qū)的特征, 與郊區(qū)土壤相比, 城區(qū)土壤C含量較高, 說明人為活動造成城區(qū)土壤明顯的有機質(zhì)富集現(xiàn)象。城區(qū)土壤P含量高于郊區(qū)和鄉(xiāng)村, 城市化過程加劇了城區(qū)土壤P的富集, 也使該地區(qū)水體(主要是滇池)面臨富營養(yǎng)化風(fēng)險。

    致謝:云南大學(xué)資源環(huán)境與地球科學(xué)學(xué)院自然地理學(xué)專業(yè)碩士研究生葉鑫、王振、袁會珍、李同艷同學(xué)在野外土壤采樣過程中給予了大力幫助, 作者在此表示衷心的感謝!

    [1] 周莉, 李保國, 周廣勝. 土壤有機碳的主導(dǎo)影響因子及其研究進展[J]. 地球科學(xué)進展, 2005, 20(1): 99–105.

    [2] WU Gang, WEI Jing, DENG Hongbing, et al. Nutrient cycling in an Alpine tundra ecosystem on Changbai Mountain, Northeast China[J]. Applied Soil Ecology, 2006, 32(2): 199–209.

    [3] PICKETT S T A, CADENASSO M L, GROVE M J, et al. Urban ecological systems: scientific foundations and a decade of progress[J]. Journal of Environmental Mana-gement, 2011, 92(11): 331–362.

    [4] 張甘霖, 朱永官, 傅伯杰. 城市土壤質(zhì)量演變及其生態(tài)環(huán)境效應(yīng)[J]. 生態(tài)學(xué)報, 2003, 23(3): 539–546.

    [5] RABBI S M F, TIGHE M, DELGADO-BAQUERIZO M, et al. Climate and soil properties limit the positive effects of land use reversion on carbon storage in Eastern Au-stra-lia[J]. Scientific Reports, 2015, 5: 17866.

    [6] LIDIA O, KAMIL P, WOJCIECH K, et al. The effect of urbanization on soil properties and soil organic carbon accumulation in tosoil of Pruszkow-a medium-sized city in the Warsaw Metropolitan Area, Poland[J]. Soil Science Annual, 2014, 65(1): 10–17.

    [7] 羅上華, 毛齊正, 馬克明, 等. 北京城市綠地表層土壤碳氮分布特征[J]. 生態(tài)學(xué)報, 2014, 34(2): 6011–6019.

    [8] 沈非, 任雅茹, 黃艷萍, 等. 蕪湖城市綠地表層土壤有機碳密度分布特征, 2018, 49(5): 1123–1129.

    [9] 許乃政, 張?zhí)伊? 王興祥, 等. 城市化進程中的土壤有機碳庫演變趨勢分析[J]. 土壤通報, 2011, 42(3): 659–663.

    [10] 呂麗平, 王登峰, 魏志遠, 等. 近20年海南島土壤有機碳時空變異[J]. 熱帶作物學(xué)報, 2016, 37(2): 338–344.

    [11] CHEN Fusheng, LI Xi, NAGLE Greg, et al. Topsoil phosphorus signature in five forest types along an urban– suburban–rural gradient in Nanchang, southern China[J]. Journal of Forestry Research, 2010, 21(1): 39–44.

    [12] YE Huichun, HUANG Yuanfang, CHEN Pengfei, et al. Effects of land use change on the spatiotemporal variability of soil organic carbon in an urban-rural ecotone of Beijing[J]. Journal of Integrative Agriculture, 2016, 15(4): 918–928.

    [13] 王啟, 李艷, 王連維, 等. 成都東部“城–郊–鄉(xiāng)”梯度綠地土壤碳氮磷化學(xué)計量特征[J]. 土壤, 2017, 49(2): 358– 363.

    [14] TOIT M J D, CILLIERS S S. Aspects influencing the selection of representative urbanization measures to quantify urban–rural gradients[J]. Landscape Ecology, 2011, 26(2): 169–181.

    [15] GIANOTTI A G S, GETSON J M, HUTYRA L R, et al. Defining urban, suburban, and rural: a method to link perceptual definitions with geospatial measures of urbani-zation in central and eastern Massachusetts[J]. Urban Ecosystems, 2016, 19(2): 823–833.

    [16] 何云玲, 張淑潔, 李同艷. 昆明市中部城–郊–鄉(xiāng)土地梯度帶的劃分及景觀分異[J]. 長江流域資源與環(huán)境. 2018, 27(9): 2022–2030.

    [17] 李家熙, 葛曉立. 城市土壤環(huán)境地球化學(xué)研究–以蘇州市為例[J]. 地質(zhì)通報, 2005, 24(8): 710–714.

    [18] LUO Shanghua, MAO Qizheng, MA Keming. Comparison on soil carbon stocks between urban and suburban topsoil in Beijing, China[J]. Chinese Geographical Science, 2014, 24(5): 551–561.

    [19] LIU Yong, WANG Cheng, YUE Wenze, et al. Storage and density of soil organic carbon in urban topsoil of hilly cities: a case study of Chongqing municipality of China. Chinese Geographical Science, 2013, 23(1): 26–34.

    [20] 馮錦, 崔東, 孫國軍, 等. 新疆土壤有機碳與土壤理化性質(zhì)的相關(guān)性[J]. 草業(yè)科學(xué), 2017, 34(4): 692–697.

    [21] 王天樂, 衛(wèi)智軍, 閆瑞瑞, 等. 氮、磷、鉀不同施肥配比對土壤養(yǎng)分的影響[J]. 北方園藝, 2017, 23: 147–154.

    [22] HU Peilei, LIU Shujuan, Ye Yingying, et al. Effects of environmental factors on soil organic carbon under natural or managed vegetation restoration[J]. Land Degradation and Development, 2018, 1: 1–11.

    [23] ARUNRAT N, PUMIJUMNONG N, HATANO R. Predicting local-scale impact of climate change on rice yield and soil organic carbon sequestration: A case study in Roi Et Province, Northeast Thailand[J]. Agricultural Systems, 2018, 164: 58–70.

    [24] ASSEFA D, REWALD B, SANDéN H, et al. Deforestation and land use strongly effect soil organic carbon and nitrogen stock in Northwest Ethiopia[J]. Catena, 2017, 153: 89–99.

    [25] 柳云龍, 章立佳, 施振香, 等. 上海城市樣帶土壤有機碳空間變異性研究[J]. 長江流域資源與環(huán)境, 2011, 20(12): 1488–1494.

    [26] 段保正, 石輝, 魏小芳, 等. 西安市城區(qū)表層土壤碳儲量與分布特征[J]. 水土保持通報, 2016, 36(6): 293–297.

    [27] 張小萌, 李艷紅, 王盼盼, 等. 烏魯木齊城市土壤有機碳空間變異研究[J]. 干旱區(qū)資源與環(huán)境, 2016, 30(2): 117–121.

    [28] 王秋兵, 段迎秋, 魏忠義, 等. 沈陽市城市土壤有機碳空間變異特征研究[J]. 土壤通報, 2009, 40(2): 252–257.

    [29] 曾宏達, 徐涵秋, 劉智才, 等. 快速城市化地區(qū)住宅用地表層土壤有機碳的變異性及其影響因素—以福州南臺島為例[J]. 生態(tài)學(xué)報, 2018, 38(4): 1427–1435.

    [30] YUAN Dagang, ZHANG Ganlin, GONG Zitong, et a1. Variations of soil phosphorus accumulation in Nanjing, China as affected by urban development[J]. Journal of Plant Nutrition& Soil Science, 2007, 170(2): 244-249.

    [31] 楊佳佳, 張向茹, 馬露莎, 等. 黃土高原刺槐林不同組分生態(tài)化學(xué)計量關(guān)系研究[J]. 土壤學(xué)報, 2014, 5l(1): 133– 142.

    [32] 劉偉玲, 張林波, 葉有華, 等. 深圳市不同土地利用類型土壤有機碳與密度特征[J]. 生態(tài)科學(xué), 2011, 30(5): 486– 492.

    [33] OKTABA L, PAZIEWSKI K, KWASOWSKI W, et al. The effect of urbanization on soil properties and soil organic carbon accumulation in topsoil of Pruszków – a medium- sized city in the Warsaw Metropolitan Area, Poland[J]. Soil Science Annual, 2014, 65(1): 10–17.

    [34] CLEVELAND C C, LIPTZIN D. C: N: P stoichiometry in soil: Is there a“Redfield ratio”for the microbial biomass?[J]. Biogeochemistry, 2007, 85(3): 235–252.

    [35] GAO Yang, HE Nianpeng, YU Guirui, et a1. Long-term effects of different land use types on C, N, and P stoi-chiometry an d storage in subtropical ecosystems: A case study in China[J]. Ecological Engineering, 2014, 67(6): 171–181.

    [36] PAUL E. Soil Microbiology, Ecology, and Bioche-mi-stry[M]. Boston: Academic, Amsterdam, 2007.

    [37] CUSACK D F. Soil nitrogen levels are linked to deco-mposition enzyme activities along an urban-remote tropical forest gradient[J]. Soil Biology & Biochemistry, 2013, 57(57): 192–203.

    [38] 曹娟, 閏文德, 項文化, 等. 湖南會同3個林齡杉木人工林土壤碳、氮、磷化學(xué)計量特征[J]. 林業(yè)科學(xué), 2015, 51(7): 1–8.

    [39] TIAN Hanqin, CHEN Guangsheng, ZHANG Chi, et a1. Pattern and variation of C: N: P ratios in China’s soils: A synthesis of observational data[J]. Biogeochemistry, 2010, 98(1/2/3): 139–151.

    Stoichiometric characteristics of topsoil C, N and P along urban-suburb-rural gradient in a rapid urbanizing area of Kunming

    HE Yunling1,*, ZHANG Shujie1, DENG Fuying1, RONG Li2, ZHANG Linyan3

    1. College of Resources Environment & Earth Science, Yunnan University, Kunming 650091, China 2. Institute of International Rivers and Eco-security, Yunnan University, Kunming 650091, China 3. College of Architecture and Urban Planning, Yunnan University, Kunming 650091, China

    Rapid urbanization is associated with physical changes in the landscape, such as increasingly impervious surface area and the replacement of natural vegetation with lawns, which can positively or negatively alter soil organic carbon (SOC) stocks. Based on remote sensing image, geographic information system (ArcGIS 10.2) and landscape pattern analysis software (Fragstats 4.2), regional urbanization levels were analyzed using principal component analysis (PCA) and multiple linear regression in the central Kunming. Then soil sampling plots were selected from typical sites based on the urban-suburb-rural gradient covering the central Kunming, and stoichiometric variability of topsoil C, N and P was analyzed. The results indicated that the spatial differentiation of soil C: N: P characteristic is caused by the rapid urbanization process. The coefficient of variation (CV) in the urban area is the highest. The CV for each element in the urban area is 84.08%, 76.95% and 81.34%, respectively, indicating that the spatial distribution of C, N and P contents in the soil surface is extremely uneven.The content of C, N and P in the topsoil is significantly higher in the urban area than in the suburbs. Compared with the soil in the suburbs, the content of C in the urban area is higher, indicating that human activities have caused obvious organic matter enrichment in the urban area. The content of topsoil P in urban areas is higher than that in the suburbs and rural areas. The urbanization process aggravates the enrichment of topsoil P in urban areas, and also exposes the water bodies (mainly Dianchi Lake) to the risk of eutrophication. The analysis of soil C, N and P contents and stoichiometric characteristics caused by human activities in the process of urbanization is of great significance for improving the level of urban soil management, improving the ecological environment and promoting the healthy and sustainable development of urban ecosystem. In the future, some leguminous nitrogen fixation plants can be selected in Kunming for urban greening to strengthen biological nitrogen fixation. The purpose is to improve soil N content and promote the absorption and utilization of P in urban soil, so as to prevent and reduce the harm of P to water.

    rapid urbanization; urban-suburb-rural gradient; C: N: P; spatial variability; central Kunming

    10.14108/j.cnki.1008-8873.2020.03.003

    S154.1

    A

    1008-8873(2020)03-017-08

    2019-04-19;

    2019-05-15

    國家自然科學(xué)基金項目(41961044); 云南省科技廳-云南大學(xué)“雙一流”建設(shè)聯(lián)合基金項目(2018FY001-017)

    何云玲(1978—), 女, 博士, 副教授, 主要從事區(qū)域生態(tài)環(huán)境變化及其影響研究, E-mail: hyl610@126.com

    何云玲

    何云玲, 張淑潔, 鄧福英, 等. 快速城市化地區(qū)昆明市“城—郊—鄉(xiāng)”梯度土壤表層碳氮磷化學(xué)計量特征[J]. 生態(tài)科學(xué), 2020, 39(3): 17–24.

    HE Yunling, ZHANG Shujie, DENG Fuying, et al. Stoichiometric characteristics of topsoil C, N and P along urban-suburb-rural gradient in a rapid urbanizing area of Kunming[J]. Ecological Science, 2020, 39(3): 17–24.

    猜你喜歡
    研究
    FMS與YBT相關(guān)性的實證研究
    2020年國內(nèi)翻譯研究述評
    遼代千人邑研究述論
    視錯覺在平面設(shè)計中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    關(guān)于遼朝“一國兩制”研究的回顧與思考
    EMA伺服控制系統(tǒng)研究
    基于聲、光、磁、觸摸多功能控制的研究
    電子制作(2018年11期)2018-08-04 03:26:04
    新版C-NCAP側(cè)面碰撞假人損傷研究
    關(guān)于反傾銷會計研究的思考
    焊接膜層脫落的攻關(guān)研究
    電子制作(2017年23期)2017-02-02 07:17:19
    久热久热在线精品观看| 啦啦啦在线观看免费高清www| 成年女人在线观看亚洲视频 | 国产精品一区www在线观看| 777米奇影视久久| 美女脱内裤让男人舔精品视频| 国产成人a区在线观看| 2021少妇久久久久久久久久久| 国产视频首页在线观看| 亚洲美女视频黄频| 国产熟女欧美一区二区| 欧美成人午夜免费资源| 久久影院123| 日韩中字成人| 国产一区二区三区综合在线观看 | 亚洲国产精品999| 久热久热在线精品观看| 国产免费一区二区三区四区乱码| 国产成人91sexporn| 在线看a的网站| 国产一区二区三区综合在线观看 | 精品人妻视频免费看| 在线播放无遮挡| 新久久久久国产一级毛片| 精品久久久久久久末码| av又黄又爽大尺度在线免费看| 日韩欧美 国产精品| 免费看光身美女| 国产色爽女视频免费观看| 97在线视频观看| 亚洲国产av新网站| 在线看a的网站| 最近最新中文字幕大全电影3| 我的女老师完整版在线观看| 日韩av不卡免费在线播放| 一个人看的www免费观看视频| 国产老妇女一区| 免费人成在线观看视频色| 国产色婷婷99| 国产成人精品一,二区| 亚洲天堂国产精品一区在线| 国产淫片久久久久久久久| 日本爱情动作片www.在线观看| 国产男女内射视频| 日本午夜av视频| 亚洲图色成人| 在线观看美女被高潮喷水网站| 黑人高潮一二区| 国产精品偷伦视频观看了| 国产av码专区亚洲av| 国产精品人妻久久久久久| 国产大屁股一区二区在线视频| 最近中文字幕高清免费大全6| 日本与韩国留学比较| 美女脱内裤让男人舔精品视频| 韩国av在线不卡| 美女xxoo啪啪120秒动态图| 中文天堂在线官网| 欧美激情久久久久久爽电影| 欧美潮喷喷水| 丝瓜视频免费看黄片| 丝袜脚勾引网站| 国产高清有码在线观看视频| 大码成人一级视频| 国产亚洲5aaaaa淫片| 最近中文字幕2019免费版| 亚洲欧美成人精品一区二区| tube8黄色片| 国产老妇伦熟女老妇高清| 成人漫画全彩无遮挡| 丰满人妻一区二区三区视频av| 亚洲av中文av极速乱| 美女脱内裤让男人舔精品视频| 亚洲最大成人中文| 中文资源天堂在线| 国产午夜精品一二区理论片| 可以在线观看毛片的网站| 亚洲国产精品国产精品| 在线免费观看不下载黄p国产| 2021少妇久久久久久久久久久| av在线播放精品| 青春草亚洲视频在线观看| 亚洲精品色激情综合| av播播在线观看一区| 你懂的网址亚洲精品在线观看| 不卡视频在线观看欧美| av在线app专区| 国产精品久久久久久精品古装| 简卡轻食公司| 男插女下体视频免费在线播放| 丰满乱子伦码专区| eeuss影院久久| 国产乱人视频| 大陆偷拍与自拍| 成年女人在线观看亚洲视频 | 亚洲,欧美,日韩| 国产色婷婷99| 丝袜喷水一区| 在线免费观看不下载黄p国产| 2021少妇久久久久久久久久久| 免费大片18禁| 精品一区在线观看国产| 麻豆国产97在线/欧美| 97在线视频观看| 七月丁香在线播放| 免费观看a级毛片全部| 制服丝袜香蕉在线| 岛国毛片在线播放| 国产伦理片在线播放av一区| 国产成人精品一,二区| 高清毛片免费看| 最近2019中文字幕mv第一页| 欧美日韩国产mv在线观看视频 | 欧美日韩视频精品一区| 国产永久视频网站| 天堂中文最新版在线下载 | freevideosex欧美| 亚洲精品成人久久久久久| 国产精品久久久久久精品电影| 搡女人真爽免费视频火全软件| 亚洲成人精品中文字幕电影| 永久免费av网站大全| 亚洲经典国产精华液单| 亚洲在久久综合| 亚洲精品一二三| 亚洲人成网站在线观看播放| 亚洲人成网站高清观看| 欧美日韩亚洲高清精品| 超碰97精品在线观看| 午夜爱爱视频在线播放| 精品久久久精品久久久| 亚洲,一卡二卡三卡| 校园人妻丝袜中文字幕| 国产视频首页在线观看| 80岁老熟妇乱子伦牲交| 日韩一区二区三区影片| 麻豆成人午夜福利视频| 禁无遮挡网站| 天天躁日日操中文字幕| 久久综合国产亚洲精品| 色哟哟·www| 亚洲av.av天堂| 全区人妻精品视频| 国产黄片美女视频| 久久久精品免费免费高清| 男人和女人高潮做爰伦理| 一级a做视频免费观看| 日本色播在线视频| 免费高清在线观看视频在线观看| 国产精品一区二区性色av| 美女脱内裤让男人舔精品视频| 国产白丝娇喘喷水9色精品| 欧美高清成人免费视频www| 夜夜爽夜夜爽视频| 国产日韩欧美在线精品| 国产成人精品一,二区| 亚洲欧美日韩东京热| 国产伦理片在线播放av一区| 国产伦精品一区二区三区四那| 欧美97在线视频| 一级爰片在线观看| 久久久久久久亚洲中文字幕| 精品亚洲乱码少妇综合久久| 久久久久久久国产电影| 中文字幕亚洲精品专区| 国产精品蜜桃在线观看| 嘟嘟电影网在线观看| 国产精品久久久久久精品古装| 热re99久久精品国产66热6| 日韩一区二区三区影片| 国产视频内射| 少妇熟女欧美另类| 亚洲国产精品成人久久小说| 18禁裸乳无遮挡免费网站照片| 欧美zozozo另类| av国产精品久久久久影院| 国产综合精华液| 超碰av人人做人人爽久久| 狂野欧美激情性xxxx在线观看| 有码 亚洲区| 久久久久久国产a免费观看| 久久久久久久午夜电影| 国产综合精华液| 最后的刺客免费高清国语| 韩国av在线不卡| 晚上一个人看的免费电影| 成人亚洲精品av一区二区| 欧美成人精品欧美一级黄| 一级av片app| 嘟嘟电影网在线观看| 国产乱人偷精品视频| 国产一区有黄有色的免费视频| 日韩电影二区| 秋霞伦理黄片| 精品视频人人做人人爽| 免费看不卡的av| 亚洲在久久综合| 色综合色国产| 你懂的网址亚洲精品在线观看| 天堂中文最新版在线下载 | 午夜日本视频在线| 国产精品偷伦视频观看了| 少妇人妻一区二区三区视频| 日韩中字成人| 国产精品av视频在线免费观看| 国内揄拍国产精品人妻在线| av免费在线看不卡| 一本色道久久久久久精品综合| 91久久精品国产一区二区三区| 麻豆国产97在线/欧美| 噜噜噜噜噜久久久久久91| 亚洲最大成人av| 亚洲精品久久午夜乱码| 超碰av人人做人人爽久久| 亚洲,一卡二卡三卡| 男男h啪啪无遮挡| 一个人观看的视频www高清免费观看| 色播亚洲综合网| 色婷婷久久久亚洲欧美| 婷婷色av中文字幕| 内地一区二区视频在线| 又黄又爽又刺激的免费视频.| 成人欧美大片| 国产毛片a区久久久久| 人妻夜夜爽99麻豆av| 亚洲精品久久久久久婷婷小说| 免费高清在线观看视频在线观看| 成人国产麻豆网| 久久久国产一区二区| 青青草视频在线视频观看| 久久久久九九精品影院| 中文乱码字字幕精品一区二区三区| 狂野欧美白嫩少妇大欣赏| 18禁动态无遮挡网站| 欧美3d第一页| 18禁动态无遮挡网站| 国产综合懂色| 人人妻人人爽人人添夜夜欢视频 | 日本黄大片高清| 成人毛片a级毛片在线播放| 亚洲综合色惰| 少妇的逼好多水| 日日摸夜夜添夜夜爱| 成年人午夜在线观看视频| 三级国产精品片| 涩涩av久久男人的天堂| 人人妻人人看人人澡| 日韩,欧美,国产一区二区三区| 精品国产一区二区三区久久久樱花 | a级一级毛片免费在线观看| 麻豆成人av视频| 九九在线视频观看精品| 日本-黄色视频高清免费观看| 国产永久视频网站| 精品久久久久久久久av| 校园人妻丝袜中文字幕| 亚洲精品日本国产第一区| 亚洲av在线观看美女高潮| 成人午夜精彩视频在线观看| 久久久久久久久久成人| 午夜福利网站1000一区二区三区| av播播在线观看一区| 一区二区av电影网| 亚洲精品国产成人久久av| 国产成人a∨麻豆精品| 少妇 在线观看| 亚洲av不卡在线观看| 午夜视频国产福利| h日本视频在线播放| 亚洲精品视频女| 免费观看a级毛片全部| 国产亚洲午夜精品一区二区久久 | 久久久久久久国产电影| 久久这里有精品视频免费| 国产成人精品福利久久| 亚洲自拍偷在线| 中文精品一卡2卡3卡4更新| 久久综合国产亚洲精品| 王馨瑶露胸无遮挡在线观看| 国产老妇伦熟女老妇高清| 成人亚洲精品av一区二区| 成人美女网站在线观看视频| 亚洲最大成人中文| 22中文网久久字幕| 久久久成人免费电影| av一本久久久久| videossex国产| 狂野欧美激情性bbbbbb| av线在线观看网站| 可以在线观看毛片的网站| 色综合色国产| 久久久亚洲精品成人影院| 久久久久久伊人网av| 狂野欧美激情性xxxx在线观看| 天美传媒精品一区二区| 中文乱码字字幕精品一区二区三区| 亚洲国产成人一精品久久久| 欧美极品一区二区三区四区| 免费看av在线观看网站| 亚洲不卡免费看| 日本三级黄在线观看| 亚洲精华国产精华液的使用体验| 国产欧美亚洲国产| 在线观看av片永久免费下载| 啦啦啦中文免费视频观看日本| 中文乱码字字幕精品一区二区三区| 晚上一个人看的免费电影| 精品久久久久久久人妻蜜臀av| 免费观看av网站的网址| 视频区图区小说| 深爱激情五月婷婷| 亚洲欧美清纯卡通| 欧美 日韩 精品 国产| 日韩一区二区视频免费看| 精品酒店卫生间| 免费看a级黄色片| 日日撸夜夜添| 国产在视频线精品| 新久久久久国产一级毛片| 国产欧美亚洲国产| www.色视频.com| 简卡轻食公司| 99re6热这里在线精品视频| 我的老师免费观看完整版| 国产伦精品一区二区三区视频9| 久久影院123| 好男人在线观看高清免费视频| 嘟嘟电影网在线观看| 久久久久久久亚洲中文字幕| 亚洲av二区三区四区| 国产 一区精品| 激情 狠狠 欧美| 久久久久网色| 久久久a久久爽久久v久久| 在线播放无遮挡| 日韩人妻高清精品专区| 国产精品99久久久久久久久| 99视频精品全部免费 在线| 国内揄拍国产精品人妻在线| 免费大片18禁| 亚洲欧美日韩另类电影网站 | 少妇人妻久久综合中文| 一区二区三区四区激情视频| 成人鲁丝片一二三区免费| 亚洲三级黄色毛片| 亚洲经典国产精华液单| 久久精品国产自在天天线| 嫩草影院精品99| 亚洲人与动物交配视频| 免费人成在线观看视频色| 国产有黄有色有爽视频| 国产精品av视频在线免费观看| 最近中文字幕高清免费大全6| 久久精品夜色国产| 免费人成在线观看视频色| 成人午夜精彩视频在线观看| 99re6热这里在线精品视频| 日韩av在线免费看完整版不卡| 日日撸夜夜添| xxx大片免费视频| 六月丁香七月| 99久久精品国产国产毛片| 国产精品麻豆人妻色哟哟久久| 人体艺术视频欧美日本| 国产精品成人在线| 久久久精品94久久精品| 成年人午夜在线观看视频| 在线观看一区二区三区激情| 久久久久久久久久人人人人人人| 久久精品久久久久久噜噜老黄| 日韩强制内射视频| 精品人妻熟女av久视频| av在线老鸭窝| 国产精品久久久久久久电影| 国产成人91sexporn| 国产毛片在线视频| 国产毛片a区久久久久| 免费看a级黄色片| 欧美性感艳星| 色视频www国产| 99re6热这里在线精品视频| 国产一区二区三区av在线| av在线天堂中文字幕| 欧美丝袜亚洲另类| 看黄色毛片网站| 精品人妻视频免费看| 欧美日韩视频高清一区二区三区二| 成年人午夜在线观看视频| av专区在线播放| 亚洲精品自拍成人| 我要看日韩黄色一级片| 超碰97精品在线观看| 亚洲国产日韩一区二区| 26uuu在线亚洲综合色| 亚洲国产精品成人综合色| 国产91av在线免费观看| 搞女人的毛片| 国产精品成人在线| 亚洲,一卡二卡三卡| 欧美日韩亚洲高清精品| 亚洲成色77777| 成人国产av品久久久| 纵有疾风起免费观看全集完整版| 男女国产视频网站| 日本黄色片子视频| 内射极品少妇av片p| 九九久久精品国产亚洲av麻豆| 亚洲va在线va天堂va国产| 精品一区二区三卡| 精品人妻偷拍中文字幕| 搡老乐熟女国产| 久热久热在线精品观看| 久久久久久久久久久免费av| 国产成人免费观看mmmm| 亚洲最大成人手机在线| 日韩成人伦理影院| 我的女老师完整版在线观看| 色婷婷久久久亚洲欧美| 国产伦精品一区二区三区四那| 免费看a级黄色片| 欧美成人a在线观看| 亚洲精品一区蜜桃| 国产精品一区二区在线观看99| 99热网站在线观看| 简卡轻食公司| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| 毛片一级片免费看久久久久| 亚洲丝袜综合中文字幕| 亚洲成人精品中文字幕电影| 深爱激情五月婷婷| 久久久久久久久久久丰满| 国产精品久久久久久久电影| 欧美另类一区| 又爽又黄无遮挡网站| 五月开心婷婷网| 综合色av麻豆| 欧美极品一区二区三区四区| 亚洲国产高清在线一区二区三| 舔av片在线| 九九爱精品视频在线观看| 亚洲av日韩在线播放| 超碰av人人做人人爽久久| 欧美极品一区二区三区四区| 全区人妻精品视频| 日韩强制内射视频| 丝瓜视频免费看黄片| 中文天堂在线官网| 久久人人爽av亚洲精品天堂 | 性插视频无遮挡在线免费观看| 一个人观看的视频www高清免费观看| av女优亚洲男人天堂| 日韩在线高清观看一区二区三区| 2021天堂中文幕一二区在线观| 国产av国产精品国产| 又黄又爽又刺激的免费视频.| 日韩一本色道免费dvd| 免费观看无遮挡的男女| 在线观看一区二区三区激情| 免费看a级黄色片| 性插视频无遮挡在线免费观看| 亚洲在线观看片| 亚洲自拍偷在线| 日韩欧美精品免费久久| 人妻 亚洲 视频| 久久精品国产自在天天线| 国产亚洲精品久久久com| 中国美白少妇内射xxxbb| 搡女人真爽免费视频火全软件| 在线 av 中文字幕| 99re6热这里在线精品视频| 七月丁香在线播放| 麻豆乱淫一区二区| av福利片在线观看| 久久久精品94久久精品| 亚洲在线观看片| 99热国产这里只有精品6| 欧美xxxx黑人xx丫x性爽| 国产成人午夜福利电影在线观看| 日韩精品有码人妻一区| 精品一区二区三区视频在线| 日韩国内少妇激情av| 97超视频在线观看视频| 人妻系列 视频| 人体艺术视频欧美日本| 我要看日韩黄色一级片| 高清在线视频一区二区三区| 亚洲激情五月婷婷啪啪| 亚洲精品影视一区二区三区av| 久久久欧美国产精品| 国产成人91sexporn| 99热全是精品| 欧美xxⅹ黑人| 久久6这里有精品| 国产成人午夜福利电影在线观看| 亚洲国产av新网站| 国产视频内射| 大话2 男鬼变身卡| 色播亚洲综合网| 国产又色又爽无遮挡免| 丰满人妻一区二区三区视频av| 欧美丝袜亚洲另类| 你懂的网址亚洲精品在线观看| 人人妻人人添人人爽欧美一区卜| 国产xxxxx性猛交| 亚洲一卡2卡3卡4卡5卡精品中文| 在线 av 中文字幕| 成人免费观看视频高清| 精品一区在线观看国产| 极品少妇高潮喷水抽搐| 久久精品亚洲av国产电影网| 欧美 亚洲 国产 日韩一| 亚洲精品美女久久av网站| 日本猛色少妇xxxxx猛交久久| 伊人久久大香线蕉亚洲五| 亚洲精品自拍成人| 亚洲综合精品二区| 精品人妻熟女毛片av久久网站| 国产亚洲欧美精品永久| 熟妇人妻不卡中文字幕| 国产精品国产三级专区第一集| 亚洲欧美清纯卡通| 国产一区二区三区综合在线观看| √禁漫天堂资源中文www| 十八禁高潮呻吟视频| 黄片无遮挡物在线观看| 久久久久网色| 大陆偷拍与自拍| 国产99久久九九免费精品| 久久久精品免费免费高清| 日日撸夜夜添| 欧美成人精品欧美一级黄| 久久天躁狠狠躁夜夜2o2o | www.熟女人妻精品国产| 欧美av亚洲av综合av国产av | 精品少妇一区二区三区视频日本电影 | 丰满饥渴人妻一区二区三| 狠狠婷婷综合久久久久久88av| 日本爱情动作片www.在线观看| 1024香蕉在线观看| 夫妻午夜视频| 精品一区二区三卡| 中国三级夫妇交换| 国产乱人偷精品视频| 免费在线观看黄色视频的| 另类精品久久| 国产极品天堂在线| 久久久久精品性色| 免费女性裸体啪啪无遮挡网站| 丝袜美足系列| 岛国毛片在线播放| 亚洲 欧美一区二区三区| 亚洲久久久国产精品| 国产片内射在线| 男女无遮挡免费网站观看| 久久久久国产一级毛片高清牌| 精品人妻一区二区三区麻豆| 日韩精品有码人妻一区| 中文字幕亚洲精品专区| 制服诱惑二区| 久久久久久免费高清国产稀缺| 日韩大码丰满熟妇| 水蜜桃什么品种好| 亚洲av日韩在线播放| 精品久久久精品久久久| 久久久久久久久久久免费av| 久久99热这里只频精品6学生| 深夜精品福利| 天堂8中文在线网| 免费看不卡的av| 精品一区在线观看国产| 精品一品国产午夜福利视频| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品国产av成人精品| 成人亚洲精品一区在线观看| 91精品国产国语对白视频| 90打野战视频偷拍视频| 中文字幕高清在线视频| 国产一区二区激情短视频 | 1024香蕉在线观看| 大片电影免费在线观看免费| 波多野结衣av一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 自线自在国产av| 成人国产麻豆网| 国产欧美日韩一区二区三区在线| 成人手机av| 亚洲成人一二三区av| 波多野结衣一区麻豆| 青春草国产在线视频| 日韩大码丰满熟妇| 韩国av在线不卡| 只有这里有精品99| 国产精品久久久久久久久免| 看非洲黑人一级黄片| 国产精品久久久久久精品古装| 欧美精品一区二区免费开放| 精品少妇黑人巨大在线播放| 国产精品 欧美亚洲| 人妻 亚洲 视频| 男女边吃奶边做爰视频| 成年动漫av网址| 久久久国产一区二区| 欧美黑人精品巨大| 久久ye,这里只有精品| 18禁国产床啪视频网站| 成人漫画全彩无遮挡| 亚洲在久久综合| 丁香六月欧美| 90打野战视频偷拍视频| 亚洲婷婷狠狠爱综合网| 色94色欧美一区二区| 夜夜骑夜夜射夜夜干| 国产福利在线免费观看视频| 九色亚洲精品在线播放|