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

    丹江鸚鵡溝小流域土壤有效錳質(zhì)量分?jǐn)?shù)分布特征

    2015-02-02 02:41:16唐輝李占斌宋曉強(qiáng)同新奇王添楊媛媛姚京威
    中國水土保持科學(xué) 2015年6期
    關(guān)鍵詞:維數(shù)分形海拔

    唐輝,李占斌,2?,宋曉強(qiáng),同新奇,王添,楊媛媛,姚京威

    (1. 西安理工大學(xué) 西北水資源與環(huán)境生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,710048,西安;2.中國科學(xué)院 水利部 水土保持

    研究所 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點(diǎn)實(shí)驗(yàn)室,712100,楊西楊凌;3.陜西省水土保持局,710004,西安)

    ?

    丹江鸚鵡溝小流域土壤有效錳質(zhì)量分?jǐn)?shù)分布特征

    唐輝1,李占斌1,2?,宋曉強(qiáng)3,同新奇3,王添1,楊媛媛1,姚京威1

    (1. 西安理工大學(xué) 西北水資源與環(huán)境生態(tài)教育部重點(diǎn)實(shí)驗(yàn)室,710048,西安;2.中國科學(xué)院 水利部 水土保持

    研究所 黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點(diǎn)實(shí)驗(yàn)室,712100,楊西楊凌;3.陜西省水土保持局,710004,西安)

    摘要:錳是土壤中的一種重要元素,其遷移過程可在一定程度上反映水土流失過程。為進(jìn)一步闡明小流域內(nèi)水土流失過程和小流域土壤有效錳的分布特征,采用DTPA-原子吸收火焰分光光度法對鸚鵡溝小流域不同地形因子、土地利用方式下土壤有效錳質(zhì)量分?jǐn)?shù)的測定,作出流域內(nèi)土壤有效錳質(zhì)量分?jǐn)?shù)的各個(gè)方向分形維數(shù),得出該小流域有效態(tài)錳的分布特征。結(jié)果表明:鸚鵡溝小流域土壤有效錳質(zhì)量分?jǐn)?shù)范圍為7.43~57.67 mg/kg,均值為39.52 mg/kg;海拔480~540 m地區(qū)質(zhì)量分?jǐn)?shù)最高,為41.61 mg/kg,變異系數(shù)從大到小為低海拔、中海拔、高海拔;8~14°坡度質(zhì)量分?jǐn)?shù)分布最高,為43.32 mg/kg,在不同坡度上變異系數(shù)從大到小為:較緩坡、微坡、急陡坡、較陡坡、緩坡;上坡位質(zhì)量分?jǐn)?shù)最高42.16 mg/kg,變異系數(shù)從大到小上坡位、中坡位、下坡位;林地質(zhì)量分?jǐn)?shù)最高42.07 mg/kg,變異系數(shù)從大到小為農(nóng)地、林地、草地;在45°方向上具有較顯著的分形特征,且變異程度最高。

    關(guān)鍵詞:土壤有效錳; 地形; 土地利用; 分形維數(shù)

    項(xiàng)目名稱: 國家科技支撐計(jì)劃課題“農(nóng)田水土保持工程與耕作關(guān)鍵技術(shù)研究” (2011BAD31B01);國家自然科學(xué)基金“土石山區(qū)坡面水沙調(diào)控對土壤氮素遷移的作用機(jī)制研究” (41401316)

    錳(Mn)是植物生長所必需的微量元素之一[1]。錳在土壤中的存在形態(tài)包括殘留態(tài)、有機(jī)態(tài)、氧化鐵結(jié)合態(tài)、易還原態(tài)、交換態(tài)和水溶態(tài),而對植物有直接效應(yīng)的主要是后三者,統(tǒng)稱為有效態(tài)錳[2],其質(zhì)量分?jǐn)?shù)能較好的反映土壤的供錳強(qiáng)度。目前對土壤有效態(tài)錳的研究主要在質(zhì)量分?jǐn)?shù)的空間分布[3]、空間變異特征[4]、影響因素[5-7]、變化規(guī)律[8],形態(tài)分布[9],以及生物氣候和人類活動對土壤有效錳質(zhì)量分?jǐn)?shù)分布的影響[10-12]等;對土壤有效錳質(zhì)量分?jǐn)?shù)的研究集中在面積較大的大尺度上,或是在田間地塊上,對于小流域尺度上不同海拔、地形、土地利用對土壤有效錳的質(zhì)量分?jǐn)?shù)分布影響研究較少。丹江中游是南水北調(diào)工程的水源區(qū)之一,該區(qū)的養(yǎng)分和水土流失備受關(guān)注,水土流失過程伴隨著土壤中養(yǎng)分和金屬元素的遷移,在一定程度上可用養(yǎng)分流失過程來表示。筆者研究丹江中游鸚鵡小流域的土壤有效錳質(zhì)量分?jǐn)?shù)分布,以闡明其分布局部的特征,為進(jìn)一步研究小流域尺度上水土流失過程奠定基礎(chǔ),也為小流域尺度上田間管理和灌水施肥提供一定參考。

    1材料與方法

    1.1研究區(qū)概況

    鸚鵡溝小流域位于陜西省商洛市商南縣城關(guān)鎮(zhèn)五里鋪村(E 110°52′16″~110°55′30″,N 33°29′55″~33°33′50″),流域總面積1.86 km2,流域內(nèi)主溝長3 232.90 m,最大主溝道比降0.01 m/m,流域坡面比降0.33 m/m,屬于多邊形水系;流域最高海拔600 m,最低海拔464 m,大部分為低山丘陵地貌,溝谷開闊;流域內(nèi)土壤以黃棕壤、風(fēng)化砂壤土為主,坡面有效土層厚度20~70 cm,溝道土質(zhì)以黃棕壤為主,土壤潮濕,土層較厚,土質(zhì)肥沃,相對耐旱,是該流域的主要產(chǎn)糧區(qū);黃棕壤分布在河道兩岸,砂壤土分布于坡面上;土壤多呈微酸性,硼、錳、鋅等微量元素較缺乏;土地覆蓋類型以農(nóng)地、林地和草地為主。喬木以櫟樹、松樹為主,灌木樹種較多且雜,草地以禾本科的草為主,農(nóng)地以小麥、玉米和花生為主;年均氣溫14 ℃,最高氣溫40.9 ℃,最低-12.2 ℃,多年平均日照時(shí)間為1 974 h,無霜期216 d,年均降水量803.2 mm,其中7—9月份的降水量占全年降水量50%左右[13]。

    1.2樣品采集及測定

    2012年對丹江鸚鵡溝小流域進(jìn)行網(wǎng)格布點(diǎn),網(wǎng)格距離100 m,采用GPS對每個(gè)采樣點(diǎn)定位(精度為±3 m)。采樣時(shí)以網(wǎng)格結(jié)點(diǎn)為圓心,清除地表植物和枯枝落葉后在5 m范圍內(nèi)采集10鉆(每隔40°方向上距圓心點(diǎn)2~5 m之間取樣,視具體環(huán)境而定)耕層土壤(0~20 cm),均勻混合后作為該點(diǎn)的土樣,共采集195個(gè)樣品。在實(shí)驗(yàn)室中,分揀去除石塊及植物根系等雜物,室溫下自然風(fēng)干后碾磨,過2 mm篩后,密封于聚乙烯塑料袋中保存待測。土壤有效錳采用DTPA浸提,原子火焰分光光度法測定,測定儀器為國產(chǎn)TAS-990原子吸收火焰分光光度計(jì)。

    1.3數(shù)據(jù)處理

    采用SPSS18.0、Excel 2010等軟件對土壤中有效猛質(zhì)量分?jǐn)?shù)的數(shù)據(jù)進(jìn)行描述統(tǒng)計(jì)、相關(guān)分析及回歸分析,GS+7.0作分形維數(shù)和半方差分析,ArcGIS 9.3地統(tǒng)計(jì)分析模塊制作空間分布圖。

    半方差函數(shù)(Semivariogram)是描述空間變量的關(guān)鍵函數(shù),用于描述變量的空間變異結(jié)構(gòu),反映不同距離觀測值之間的變化:

    (1)

    式中:r(h)為是間距為h的半方差函數(shù)值,在一定范圍內(nèi)隨h的增加而增大;N(h)為間隔距離為h的樣點(diǎn);Z(Xi)和Z(Xi+h)分別是區(qū)域化變量Z(X)在空間位置Xi和Xi+h的實(shí)測值。

    分形維數(shù)[14]

    D=2-H/2。

    (2)

    式中H為步長與半方差取對數(shù)值擬合直線關(guān)系中的斜率。

    圖1  研究區(qū)土地利用及采樣點(diǎn)分布圖Fig.1  Distribution of land use and sampling sites in  the study area

    2結(jié)果與分析

    2.1地形對土壤有效錳質(zhì)量分?jǐn)?shù)的影響

    2.1.1不同海拔土壤有效錳質(zhì)量分?jǐn)?shù)的分布選擇了已知海拔值的165個(gè)樣本進(jìn)行分析,該流域內(nèi)海拔為420~600 m,按照海拔跨度以每60 m為間隔把流域海拔值分為3個(gè)等值段,420~480 m為低海拔,480~540 m為中海拔,540~600 m為高海拔,分別分析每個(gè)海拔段上的樣本特性。如表1所示,海拔在480~540 m之間的有效錳質(zhì)量分?jǐn)?shù)均值要大于海拔420~480 m與540~600 m,即土壤有效錳在流域內(nèi)中海拔區(qū)域質(zhì)量分?jǐn)?shù)分布較高。土壤有效錳質(zhì)量分?jǐn)?shù)在每個(gè)海拔段的變異程度不斷增大,分別為12%、23%、32%,都在中等強(qiáng)度變異范圍內(nèi)。說明高海拔區(qū)土壤有效錳質(zhì)量分?jǐn)?shù)離散程度較大,中海拔區(qū)次之,低海拔區(qū)最小。顯著性檢驗(yàn)表明420~540 m與540~600 m海拔上的有效錳質(zhì)量分?jǐn)?shù)差異不顯著,480~540 m海拔上的有效錳質(zhì)量分?jǐn)?shù)較其他2個(gè)海拔段上的有效錳質(zhì)量分?jǐn)?shù)差異顯著,3個(gè)海拔段的有效錳質(zhì)量分?jǐn)?shù)概率分布都呈尖瘦型,420~480 m海拔段分布右偏,480~540 m和540~600 m段分布左偏。高海拔區(qū)的有效錳被淋溶隨徑流向下遷移,有效錳質(zhì)量分?jǐn)?shù)在中下部累積,隨海拔的降低,土壤有效錳質(zhì)量分?jǐn)?shù)有一定程度的增加趨勢[15]。中海拔區(qū)的地勢比高海拔區(qū)平緩,而低海拔區(qū)有河流分布,有效錳隨水流遷移,導(dǎo)致了流域內(nèi)有效錳質(zhì)量分?jǐn)?shù)在海拔上分呈現(xiàn)中海拔區(qū)高,低海拔和高海拔區(qū)較低的趨勢。

    2.1.2不同坡度土壤有效錳質(zhì)量分?jǐn)?shù)分布對所取得的樣本按水土保持綜合治理通則(GBT 15772—2008)中坡度等級分類,統(tǒng)計(jì)結(jié)果如表2所示,較緩坡(5~7°)的變異系數(shù)最大為41%,微坡(0~4°)次之,較陡坡度上土壤有效錳的分布比小坡度上的分布均勻。緩坡(8~14°)上有效錳質(zhì)量分?jǐn)?shù)平均值最大,且分布最均勻,急陡坡(≥35°)上質(zhì)量分?jǐn)?shù)次之;微坡、較緩坡、陡坡(25~34°)上土壤有效錳質(zhì)量分?jǐn)?shù)均值差異不大。由此可見,土壤有效錳在緩坡、較陡坡、急坡上易發(fā)生富集。

    表1 不同海拔土壤有效錳質(zhì)量分?jǐn)?shù)顯著性檢驗(yàn)及基本統(tǒng)計(jì)

    注:顯著性差異采用Ducan’s multiple range test方法在α=0.05水平下分析,a和b指顯著性檢驗(yàn)差異順序,不同順序表示差異性顯著,相同順序表示差異性不顯著,下同。Note:The significance difference is analyzed by Ducan’s multiple range test. Values within each column with the same letter are not significantly different (α=0.05). The same below.

    表2 不同坡度土壤有效錳質(zhì)量分?jǐn)?shù)顯著性及基本統(tǒng)計(jì)

    2.1.3不同坡位土壤有效錳質(zhì)量分?jǐn)?shù)分布把所取得的樣本按照上、中和下坡位[16]分類,它們的土壤有效錳質(zhì)量分?jǐn)?shù)均值分別為44.16 、39.24、40.87 mg/kg,呈現(xiàn)出中坡位質(zhì)量分?jǐn)?shù)最低的現(xiàn)象;中坡位質(zhì)量分?jǐn)?shù)分布呈平緩型,上坡位和下坡位呈尖瘦型,均值概率分布均是偏左。上、中、下坡位的變異系數(shù)分別為31%、26%、6%,表明土壤有效錳質(zhì)量分?jǐn)?shù)的離散程度逐漸減低,即土壤有效錳的質(zhì)量分?jǐn)?shù)分布在坡面上從上到下均勻程度越來越高。

    由于地表徑流和地下徑流對土壤的淋溶作用,上坡位降雨產(chǎn)生的地表徑流和地下徑流量較小,對土壤的淋溶作用較弱,來自上坡位的地表徑流和地下徑流使得中坡位徑流量增大,導(dǎo)致徑流對土壤的淋溶作用增強(qiáng),而且對坡面的沖刷作用增強(qiáng),土壤中的有效錳隨徑流泥沙的遷移量增大。下坡面是徑流和泥沙的匯集,土壤有效錳質(zhì)量分?jǐn)?shù)理應(yīng)最高[17];但由于坡腳處徑流量達(dá)到了坡面的最大值,徑流泥沙不會全部停留在坡腳處,會一直沿平地前進(jìn),直到?jīng)]有流速和流量為止。這解釋了下坡位土壤有效錳質(zhì)量分?jǐn)?shù)比較均勻的特征。

    表3 不同坡位土壤有效錳質(zhì)量分?jǐn)?shù)顯著性及基本統(tǒng)計(jì)

    2.2不同土地利用方式下土壤有效錳質(zhì)量分?jǐn)?shù)分布

    鸚鵡溝小流域的土地利用方式大致分為草地 、林地、農(nóng)地3種,土地利用圖(來源于Quickbird遙感影像)顯示農(nóng)地面積最大,占總面積的61.9%,其次是林地,占總面積的28.7%,面積最小的是草地,占總面積的9.4%。林地的土壤錳質(zhì)量分?jǐn)?shù)均值最大,為42.07 mg/kg,農(nóng)地次之,為40.04 mg/kg,最小的是草地,為31.36 mg/kg。

    由此可得出,草地土壤錳質(zhì)量分?jǐn)?shù)明顯低于農(nóng)地和林地,這主要是因?yàn)椴荼局参锏母鶎τ行B(tài)錳有很強(qiáng)的富集能力[18],引起草地上土壤中有效錳的質(zhì)量分?jǐn)?shù)降低。草地、林地、農(nóng)地土壤有效錳質(zhì)量分?jǐn)?shù)的變異系數(shù)分別為43%、20%、2%,三者的差異非常大。這是由于流域內(nèi)草地疏密程度不均,造成了單位面積上草本植物的根對土壤有效錳富集能力不一,最終導(dǎo)致土壤中的有效錳質(zhì)量分?jǐn)?shù)不均,離散程度大。對于林地來說,農(nóng)地土壤有效錳質(zhì)量分?jǐn)?shù)的變異程度非常小,農(nóng)地相對于草地和林地來說土壤結(jié)構(gòu)單一,說明耕作改變了土壤的空間變異特征,常年的耕作使得土壤的性質(zhì)差異不大,土壤有效錳質(zhì)量分?jǐn)?shù)分布均勻[19]。

    2.3土壤有效錳質(zhì)量分?jǐn)?shù)空間結(jié)構(gòu)分析

    2.3.1分形維數(shù)采用GS+7.0地統(tǒng)計(jì)軟件中的分析維數(shù)計(jì)算模塊(Fractal)對該小流域的土壤有效錳質(zhì)量分?jǐn)?shù)進(jìn)行了分形維數(shù)的分析(圖2),目前應(yīng)用分形理論研究土壤化學(xué)性質(zhì)的空間分布特征[20]時(shí),多數(shù)情況下用單一分形維數(shù)表征空間分布特征[21]。土壤有效錳質(zhì)量分?jǐn)?shù)分形維數(shù)為1.965(R2=0.357)。

    表4 不同土地利用方式下土壤有效錳質(zhì)量分?jǐn)?shù)顯著性及基本統(tǒng)計(jì)

    D0—分形維數(shù);SE—標(biāo)準(zhǔn)誤差;R2—確定系數(shù);n—參與計(jì)算點(diǎn)個(gè)數(shù)。D0: fractal dimension; SE: standar error; R2: coefficient of determination; n: the points number which participate calculation. 圖2 土壤有效錳質(zhì)量分?jǐn)?shù)的各向同性分形維數(shù)Fig.1  Isotropic fractal dimension of the content of  soil available Mn

    圖3 0°、45°、90°、135°方向上土壤有效錳質(zhì)量分?jǐn)?shù)各向異性分形維數(shù)Fig.3 Anisotropic fractal dimension of the content of soil available Mn at 0°, 45°, 90°, 135°

    地統(tǒng)計(jì)學(xué)中的各向同性是指在計(jì)算變異函數(shù)值時(shí)不區(qū)分方向,只要距離為h的點(diǎn)均可參與計(jì)算;而各向異性是相對各向同性而言的,樣本在不同方向的變異性有所差異。土壤的分形維數(shù)也有此特征,在不同的方向上表現(xiàn)出不同的分形特征。

    以平面直角坐標(biāo)系橫縱軸正向分別代表正東和正北方向,對土壤有效錳質(zhì)量分?jǐn)?shù)的分形維數(shù)在0°、45°、90°、135°方向進(jìn)行了分析,其中45°方向上分形特征最明顯(R2=0.758,為4個(gè)方向上的最大值,且更為接近1),其他幾個(gè)方向的分形特征不太明顯。4個(gè)方向上的分形維數(shù)值分別為1.989、1.840、1.964、1.983,從數(shù)值大小上可以看出,土壤錳質(zhì)量分?jǐn)?shù)45°方向上均勻程度最低,0°方向上均勻程度最高。如圖3所示,流域內(nèi)河流的流向應(yīng)是造成這一現(xiàn)象的原因之一,河流的流向大致是45°方向,結(jié)合各向異性的分形維數(shù)特征,表明了徑流的淋溶作用是造成土壤有效錳質(zhì)量分?jǐn)?shù)各向異性的原因之一。

    2.3.2空間格局用GS+7.0土壤有效錳質(zhì)量分?jǐn)?shù)做半方差分析,得出有效錳質(zhì)量分?jǐn)?shù)分布的最優(yōu)模型為指數(shù)模型,基于ArcGIS 9.3地統(tǒng)計(jì)模塊進(jìn)行插值(圖4)。質(zhì)量分?jǐn)?shù)在28.82~31.32 mg/kg之間的面積的占流域總面積的4.02%,31.32~33.82 mg/kg之間的面積占5.69%,33.82~36.32 mg/kg之間的面積占6.71%,36.32~38.82 mg/kg之間的面積占15.34%,32.82~41.32 mg/kg之間的占17.37%,41.32~43.82 mg/kg之間的占29.17%,43.82~47.28 mg/kg之間的占21.72%。由土壤有效錳質(zhì)量分?jǐn)?shù)Kriging簡單預(yù)測插值圖可看出,鸚鵡溝小流域有效錳質(zhì)量分?jǐn)?shù)較大的區(qū)域主要分布在流域的中西部,質(zhì)量分?jǐn)?shù)分布最大區(qū)域(41.32~43.82 mg/kg)主要分布在流域的正南、正北以及正西方向,呈現(xiàn)出從東部的逐漸上升到中部然后中部到西部逐漸下降的趨勢。該流域水系的流向是自東向西,徑流的匯聚過程伴隨著重金屬元素的遷移,成為導(dǎo)致鸚鵡溝小流域土壤有效錳質(zhì)量分?jǐn)?shù)分布不均勻的因素之一,東部集水區(qū)質(zhì)量分?jǐn)?shù)較少到中部匯聚到高峰;因西部處于流域內(nèi)水系的下游,地勢平坦且不是主要匯水區(qū),土壤中有效錳的質(zhì)量分?jǐn)?shù)反而沒有中部高。

    圖4 研究區(qū)土壤有效錳質(zhì)量分?jǐn)?shù)分布Fig.4  Distribution of the content of soil available Mn  in study area

    3討論

    該研究區(qū)內(nèi)本研究結(jié)果是在小流域的范圍內(nèi)取得的,對相近地區(qū)或者是土壤類型相似地區(qū)土壤的有效錳質(zhì)量分?jǐn)?shù)值有一定的參考。土壤有效錳質(zhì)量分?jǐn)?shù)與土壤母質(zhì)、土地利用方式、地形、土壤pH值顯著相關(guān)[22],這與本研究中結(jié)果一致。王小兵等[23]得出有效錳質(zhì)量分?jǐn)?shù)與海拔顯著正相關(guān);本文中海拔高程差小,可能在較大高程上統(tǒng)計(jì)分類不顯著,所以筆者研究中海拔上的分布結(jié)果僅適合于小尺度上的應(yīng)用。文中土地利用上的分布趨勢與文獻(xiàn)[6]中的一致;但本文中有效錳質(zhì)量分?jǐn)?shù)值高,文獻(xiàn)[6]研究區(qū)域在西藏,原因可能在于氣候方面,氣候條件影響土壤性質(zhì)。除了以上提到的影響因素外,測試樣品的干濕程度[24]也會影響到最終的測定值。盡管如此,本研究從地形和土地利用2個(gè)方面探討了土壤有效錳在小流域上的分布特征,并得出了其質(zhì)量分?jǐn)?shù)的空間分布圖,對于指導(dǎo)田間施肥管理和研究水土流失過程有著重要的意義,施肥與產(chǎn)量的關(guān)系、有效錳質(zhì)量分?jǐn)?shù)與土壤侵蝕量間的關(guān)系有待進(jìn)一步研究。

    4結(jié)論

    該研究區(qū)內(nèi)土壤有效錳總體平均質(zhì)量分?jǐn)?shù)是39.52 mg/kg,變異系數(shù)為28%,統(tǒng)計(jì)分布對照劉錚等[1]給出的土壤易還原態(tài)錳質(zhì)量分?jǐn)?shù)臨界值100 mg/kg,鸚鵡溝小流域的土壤有效態(tài)錳缺乏,局部區(qū)域更為缺乏。以60m為間隔的海拔區(qū)間上土壤有效錳的分布仍有統(tǒng)計(jì)學(xué)上的差異,這對于研究海拔因素對有效錳的質(zhì)量分?jǐn)?shù)有一定的參考價(jià)值,以后可在更大尺度上進(jìn)行相關(guān)研究。

    1)在此小流域研究得出,中海拔(480~540 m)土壤有效錳質(zhì)量分?jǐn)?shù)最高,為41.61 mg/kg;有效錳在8~15°的坡上發(fā)生富集,不同的坡位上也有顯著的差異,上坡位其質(zhì)量分?jǐn)?shù)最高,并表現(xiàn)為從上到下均勻程度越來越高。

    2)不同土地利用上有效錳的分布最為明顯,林地的土壤錳質(zhì)量分?jǐn)?shù)均值最大,為42.07 mg/kg,農(nóng)地次之,為40.04 mg/kg,最小的是草地,31.36 mg/kg,草地土壤錳質(zhì)量分?jǐn)?shù)明顯低于農(nóng)地和林地。

    3)空間分析表明,流域土壤有效錳質(zhì)量分?jǐn)?shù)的分形維數(shù)為1.965,且表現(xiàn)出各向異性,流域內(nèi)土壤有效錳隨徑流遷移,表現(xiàn)在東部分布較低,中西部較高。

    參考文獻(xiàn)5

    [1]劉錚.土壤與植物中錳的研究進(jìn)展[J].土壤科學(xué)進(jìn)展,1991,19(6):1-10

    [2]丁維新.土壤中錳賦存形態(tài)的地域性差異研究[J].長江流域資源與環(huán)境,1995,4(1):38-44

    [3]王新中,劉國順,楊程,等.小尺度下煙田土壤有效態(tài)微量元素空間分布格局研究[J].中國煙草科學(xué),2011,32(1):12-16

    [4]王雪婧,常棟.植煙前后土壤有效錳含量的空間變異特征[J]河南農(nóng)業(yè)大學(xué)學(xué)報(bào),2013,47(3):334-339

    [5]李飛,張亞歡,周聰,等.不同區(qū)域土壤中錳鉛和鉻的化學(xué)形態(tài)分析研究[J].微量元素與健康研究,2013,30(5):42-43

    [6]蒲玉林,龍高飛,茍文平,等.西藏土壤有效錳含量及其影響因子分析[J].西南師范大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,35(6):163-168

    [7]Heal K V. Manganese and land-use in upland catchments in Scotland[J]. Science of the Total Environment, 2001, 265(1/2/3):169-179

    [8]田茂成,黎娟,田峰,等.湘西植煙土壤有效錳含量及變化規(guī)律研究[J].湖北農(nóng)業(yè)科學(xué),2013,52(17):4013-4016

    [9]魏孝榮,邵明安.黃土溝壑區(qū)小流域不同地形條件下土壤錳的形態(tài)分布及其有效性[J].植物營養(yǎng)與肥料學(xué)報(bào),2008,14 (3):439-444

    [10] 丁維新.生物氣候帶對土壤中錳的賦存形態(tài)及其有效性的影響[J].地理科學(xué),1995,15(4):365-373

    [11] 楊野,王麗,郭蘭萍,等.三七不同間隔年限種植土壤中微量元素動態(tài)變化規(guī)律研究[J].中國中藥雜志,2014,39(4):580-587

    [12] 張楊珠,吳名宇,李順義,等.稻作制、有機(jī)肥和地下水位對紅壤性水稻土全錳及不同形態(tài)有效錳含量的影響[J].農(nóng)業(yè)現(xiàn)代化研究,2008,29(3):357-360

    [13] 徐國策,李占斌,李鵬,等. 丹江鸚鵡溝小流域土壤侵蝕和養(yǎng)分損失定量分析[J]. 農(nóng)業(yè)工程學(xué)報(bào),2013,29(10):160-167

    [14] Burrough P A. Multiscale sources of spatial variation in soil—The application of fractal concepts to nested levels of soil avaliation [J].European Journal of Soil science,1983,34 (3):577-597

    [15] 郜文軍,王印傳,霍習(xí)良,等.山區(qū)小流域不同海拔和坡位土壤養(yǎng)分分布及相關(guān)性分析[J].土壤肥料科學(xué),2008(3):248-252

    [16] 高雪松,鄧良基,張世熔.不同利用方式與坡位土壤物理性質(zhì)及養(yǎng)分特征分析[J].水土保持學(xué)報(bào),2005,19(2):53-56

    [17] Pierson F B, Mulla D J. Aggregate stability in the Palouse region of Washington: effect of landscape position [J]. Soil Science Society of America Journal, 1990, 54(5):1407-1412

    [18] 周曉梅,郭繼勛,李玉.松嫩平原羊草草地土壤-植物間微量元素錳的時(shí)空動態(tài)[J].東北林業(yè)大學(xué)學(xué)報(bào),2008,36(5):36-38

    [19] Brubaker S C, Jones A J, Lewis D T, et al. Soil properties associated with landscape position [J].Soil Science Society of America Journal,1993,57(1):235-239

    [20] 龔元石,廖超子,李保國. 土壤含水量和容重的空間變異及其分形特征[J].土壤學(xué)報(bào),1998,35(1):10-15

    [21] 方萍,呂成文,朱艾莉.分形方法在土壤特性空間變異研究中的應(yīng)用[J].土壤,2011,43(5):710-713

    [22] Xu Xiaoxun, Zhang Shirong, Ye Wen, et al. Spatial variability characteristics of soil available Mn and Zn in the middle reaches of Tuojiang River Basin[C]∥Advances in Environmental Science and Engineering. Part 3, 2012: 2204-2212

    [23] 王小兵,周翼衡,李強(qiáng),等. 曲靖煙區(qū)海拔高度分布特點(diǎn)及其與土壤養(yǎng)分的關(guān)系[J].煙草農(nóng)學(xué),2013,11:86-90

    [24] Anjos C, Magalh?es M C F, Abreu M M. Metal (Al, Mn, Pb and Zn) soils extractable reagents for available fraction assessment: comparison using plants, and dry and moist soils from the Bra?al abandoned lead mine area, Portugal[J]. Journal of Geochemical Exploration, 2012, 113: 45-55

    (責(zé)任編輯:程云郭雪芳)

    Distribution characteristics of the content of available Mn in Yingwugou

    watershed, Shaanxi Province

    Tang Hui1, Li Zhanbin1,2, Song Xiaoqiang3, Tong Xinqi3, Wang Tian1, Yang yuanyuan1, Yao Jingwei1

    (1.Key Laboratory of Northwest Water Resources and Environment Ecology of Ministry of Education, Xi′an University of Technology,

    710048, Xi′an, China; 2. Key Laboratory of Soil Erosion and Dry-land Farming on the Loess Plateau, Institute of Soil and Water

    Conservation, Chinese Academy of Sciences and Ministry of Water Resources, 712100, Yangling, Shaanxi, China;3. Shaanxi Bureau of

    Soil and Water Conservation, 710004, Xi′an, China)

    Abstract:Manganese (Mn) is one of important elements in soil and its transportation process is an reflection of soil and water loss. In order to understand more deeply about the process of soil and water loss, the content of available Mn in Yingwugou watershed located in Shaanxi Province was investigated. This study employed DTPA-atomic absorption fire spectrophotometry to measure the content of available Mn in different topography and land use types of the watershed. Then fractal dimension of the content of available Mn in a few directions was calculated and finally the distribution of the content of available Mn was clarified. The results showed that the content of soil available Mn at the Yingwugou watershed ranged from 7.43-57.67 mg/kg, 39.52 mg/kg on average. 1) Middle elevation had the highest content (41.61 mg/kg) and the value of the coefficient variation (CV) from high to low was ordered distinctively as low, middle, and high, suggesting that the available Mn in soil is moved from high to low elevations by runoff in the watershed after leaching, gathers at middle or low elevation areas, and finally the available Mn at low elevation areas is carried away by flowing river. The highest content of available Mn (43.32 mg/kg) could be found on slopes with degree ranging from 8-14° and the value of CV from high to low was ranked as very gentle slope, micro-slope, very steep slope, steep slope and gentle slope; the available Mn was easier to aggregate at very gentle slope, micro-slope, and steep slope. The upslope position had the highest content (42.16 mg/kg) and the value of CV from high to low was ordered as down slope, middle slope, upslope, indicating that runoff is the main force to move the available Mn. 2) The woodland had the highest content of 42.07 mg/kg and the value of CV from high to low was ranked as farmland, woodland, grassland. The woodland and the grassland had more complicated soil structures because of more biocenoses and less disturbance. The grass roots have a high ability to utilize the available Mn and the soil texture of farmland is highly uniform. All of these reasons lead to the current distribution characteristics of content of available Mn in the different landuse types. 3) The fractal characteristics were obvious at the 45° direction, and extent of variation was the highest. The high content of Mn was distributed in the west, south, and north of the watershed. The total spatial pattern of the watershed was that the content of available Mn increased from east to west and then decreased, but still higher than the east. These results would provide reference to describe the process of soil and water loss.

    Keywords:soil available Mn; topography; land use; fractal dimension

    通信作者?簡介: 李占斌(1962—),男,博士,研究員。主要研究方向:水文水資源、土壤侵蝕與水土保持等。E-mail:zhanbinli@126.com

    作者簡介:第一 唐輝(1989—),男,碩士。主要研究方向:土壤侵蝕與水土保持。E-mail:tanghuicq@126.com

    收稿日期:2015-05-15修回日期: 2015-11-03

    中圖分類號:S157.1; S158.2

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1672-3007(2015)06-0076-07

    猜你喜歡
    維數(shù)分形海拔
    β-變換中一致丟番圖逼近問題的維數(shù)理論
    高海拔的甘南
    散文詩(2021年22期)2022-01-12 06:14:04
    感受分形
    一類齊次Moran集的上盒維數(shù)
    南美三大高海拔球場
    分形之美
    海拔4000米的父愛
    分形空間上廣義凸函數(shù)的新Simpson型不等式及應(yīng)用
    關(guān)于齊次Moran集的packing維數(shù)結(jié)果
    涉及相變問題Julia集的Hausdorff維數(shù)
    色综合色国产| 日韩欧美在线乱码| 成人漫画全彩无遮挡| 男女啪啪激烈高潮av片| 成年版毛片免费区| 国产单亲对白刺激| 精品欧美国产一区二区三| 国产国拍精品亚洲av在线观看| 精品人妻视频免费看| 欧美日韩一区二区视频在线观看视频在线 | 成人亚洲精品av一区二区| 性色avwww在线观看| 久久久精品欧美日韩精品| 亚洲激情五月婷婷啪啪| 免费av不卡在线播放| 青春草视频在线免费观看| 国产精品伦人一区二区| 亚洲五月天丁香| 成人特级av手机在线观看| 成人三级黄色视频| 亚洲成人av在线免费| 十八禁国产超污无遮挡网站| 51国产日韩欧美| 国产精品三级大全| avwww免费| 亚洲成人精品中文字幕电影| 国产精品久久久久久久电影| 伦理电影大哥的女人| 天堂网av新在线| 国产精品嫩草影院av在线观看| 欧美性感艳星| 亚洲一级一片aⅴ在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产爱豆传媒在线观看| 国产白丝娇喘喷水9色精品| 日韩欧美在线乱码| 国产综合懂色| 成人鲁丝片一二三区免费| 午夜激情欧美在线| 老司机福利观看| 观看免费一级毛片| 国产成人一区二区在线| 成人美女网站在线观看视频| 麻豆乱淫一区二区| 老司机影院成人| 18禁裸乳无遮挡免费网站照片| 国产精品日韩av在线免费观看| 在线免费观看的www视频| 亚洲精品亚洲一区二区| 日韩av不卡免费在线播放| 精品熟女少妇av免费看| 小说图片视频综合网站| 日韩在线高清观看一区二区三区| 久久久久久久久久久丰满| 色视频www国产| 美女 人体艺术 gogo| 欧美色视频一区免费| 午夜福利在线观看吧| 综合色丁香网| 精品久久久久久久久av| 亚洲一级一片aⅴ在线观看| 深夜a级毛片| 我的老师免费观看完整版| 久久久久免费精品人妻一区二区| 十八禁国产超污无遮挡网站| 特大巨黑吊av在线直播| 老司机福利观看| 亚洲欧美精品综合久久99| av在线播放精品| 成年女人永久免费观看视频| 国产成年人精品一区二区| 日日啪夜夜撸| 亚洲激情五月婷婷啪啪| 波野结衣二区三区在线| 亚洲经典国产精华液单| 国产精品1区2区在线观看.| 日日摸夜夜添夜夜添小说| 国产 一区 欧美 日韩| 成年女人毛片免费观看观看9| av女优亚洲男人天堂| 一个人观看的视频www高清免费观看| 日本免费a在线| 欧美成人a在线观看| 91久久精品国产一区二区三区| 成人三级黄色视频| 搡老岳熟女国产| 日韩精品有码人妻一区| 男人舔奶头视频| 三级国产精品欧美在线观看| 色哟哟哟哟哟哟| 在线天堂最新版资源| 国产白丝娇喘喷水9色精品| 久久鲁丝午夜福利片| 国产伦精品一区二区三区四那| 亚洲图色成人| 晚上一个人看的免费电影| 免费在线观看成人毛片| 夜夜夜夜夜久久久久| 久久人人爽人人片av| 少妇熟女欧美另类| 99热6这里只有精品| 久久久国产成人免费| 简卡轻食公司| 99久国产av精品国产电影| 老司机影院成人| 国产精品久久视频播放| 18禁在线播放成人免费| 成人三级黄色视频| 成人高潮视频无遮挡免费网站| 亚洲国产日韩欧美精品在线观看| 少妇的逼水好多| 人妻丰满熟妇av一区二区三区| 午夜久久久久精精品| 久久精品国产99精品国产亚洲性色| 亚洲国产精品成人久久小说 | 在线免费观看的www视频| 欧美bdsm另类| 一级a爱片免费观看的视频| 久久人人精品亚洲av| 亚洲中文字幕一区二区三区有码在线看| 高清午夜精品一区二区三区 | 亚洲精华国产精华液的使用体验 | 欧美高清性xxxxhd video| 极品教师在线视频| 久久精品国产99精品国产亚洲性色| 国产单亲对白刺激| 久久草成人影院| 免费在线观看成人毛片| 秋霞在线观看毛片| 亚洲av不卡在线观看| 欧洲精品卡2卡3卡4卡5卡区| 最新在线观看一区二区三区| 麻豆精品久久久久久蜜桃| 夜夜夜夜夜久久久久| 精品久久久久久成人av| 久久久久久久久久久丰满| 中文字幕免费在线视频6| 国产成人aa在线观看| 国产精品美女特级片免费视频播放器| 午夜老司机福利剧场| 99在线视频只有这里精品首页| 1024手机看黄色片| 一区二区三区免费毛片| 18禁在线无遮挡免费观看视频 | 亚洲中文日韩欧美视频| 大又大粗又爽又黄少妇毛片口| 亚洲丝袜综合中文字幕| 国产激情偷乱视频一区二区| 久久精品国产鲁丝片午夜精品| 成年女人毛片免费观看观看9| 男人的好看免费观看在线视频| 精品欧美国产一区二区三| 久久精品国产亚洲av天美| 国产久久久一区二区三区| 午夜久久久久精精品| 久久久久久久久久久丰满| 国产v大片淫在线免费观看| 直男gayav资源| 中文字幕久久专区| 国产av麻豆久久久久久久| 成人亚洲欧美一区二区av| av中文乱码字幕在线| 久久这里只有精品中国| 国产精品三级大全| 国产蜜桃级精品一区二区三区| 少妇人妻一区二区三区视频| 国产真实伦视频高清在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品一区av在线观看| 永久网站在线| 人人妻,人人澡人人爽秒播| 午夜福利高清视频| 亚洲真实伦在线观看| 午夜a级毛片| 十八禁网站免费在线| 噜噜噜噜噜久久久久久91| or卡值多少钱| 99国产精品一区二区蜜桃av| 免费大片18禁| 亚洲七黄色美女视频| 一区二区三区免费毛片| 免费在线观看成人毛片| 欧美高清性xxxxhd video| 久久久久久久亚洲中文字幕| 97超碰精品成人国产| 免费av不卡在线播放| 看免费成人av毛片| 亚洲av五月六月丁香网| av在线蜜桃| av国产免费在线观看| 亚洲av免费在线观看| 淫秽高清视频在线观看| 最近中文字幕高清免费大全6| 99久久精品国产国产毛片| 性插视频无遮挡在线免费观看| 精品人妻熟女av久视频| 精品日产1卡2卡| 联通29元200g的流量卡| 久久久久国产网址| 国产精品久久久久久亚洲av鲁大| 男人舔女人下体高潮全视频| 国产蜜桃级精品一区二区三区| 国产伦精品一区二区三区视频9| 床上黄色一级片| 国产一区二区三区在线臀色熟女| 亚洲精品影视一区二区三区av| 久久久久免费精品人妻一区二区| 日日摸夜夜添夜夜添av毛片| 成年女人看的毛片在线观看| 日本与韩国留学比较| 亚洲美女视频黄频| 成年女人永久免费观看视频| 日本爱情动作片www.在线观看 | 欧美日韩国产亚洲二区| 中文亚洲av片在线观看爽| 久久人人精品亚洲av| 黑人高潮一二区| 偷拍熟女少妇极品色| 97热精品久久久久久| 少妇人妻精品综合一区二区 | 精品免费久久久久久久清纯| 亚洲va在线va天堂va国产| 最新在线观看一区二区三区| 亚洲在线观看片| 国产精品不卡视频一区二区| 不卡一级毛片| 黄片wwwwww| 听说在线观看完整版免费高清| 狠狠狠狠99中文字幕| 韩国av在线不卡| 淫妇啪啪啪对白视频| 亚洲最大成人中文| 国产精品1区2区在线观看.| 神马国产精品三级电影在线观看| 蜜臀久久99精品久久宅男| 国产乱人偷精品视频| 亚洲乱码一区二区免费版| 成年女人永久免费观看视频| 卡戴珊不雅视频在线播放| 免费在线观看影片大全网站| 日韩一本色道免费dvd| 色综合亚洲欧美另类图片| 露出奶头的视频| 免费人成视频x8x8入口观看| 校园春色视频在线观看| 国产av麻豆久久久久久久| 夜夜看夜夜爽夜夜摸| 97人妻精品一区二区三区麻豆| 舔av片在线| 国产淫片久久久久久久久| 内地一区二区视频在线| 97超碰精品成人国产| 精品一区二区三区人妻视频| 美女免费视频网站| 亚洲综合色惰| 天美传媒精品一区二区| 亚洲成人久久爱视频| 精品久久久久久久久久免费视频| 99视频精品全部免费 在线| 国产不卡一卡二| 有码 亚洲区| 春色校园在线视频观看| 女的被弄到高潮叫床怎么办| 尾随美女入室| 久久欧美精品欧美久久欧美| 国产一区二区在线av高清观看| 精品人妻熟女av久视频| 十八禁国产超污无遮挡网站| 国产欧美日韩精品一区二区| 午夜福利在线观看免费完整高清在 | 国产不卡一卡二| 少妇的逼好多水| 久久天躁狠狠躁夜夜2o2o| 精品99又大又爽又粗少妇毛片| 国产精品不卡视频一区二区| 国产精品综合久久久久久久免费| 黄片wwwwww| av在线天堂中文字幕| 免费看日本二区| 亚洲性久久影院| 久久欧美精品欧美久久欧美| or卡值多少钱| 午夜激情福利司机影院| av专区在线播放| 国产老妇女一区| 男女下面进入的视频免费午夜| 18禁在线播放成人免费| 最近的中文字幕免费完整| 啦啦啦观看免费观看视频高清| 看十八女毛片水多多多| 久久精品国产自在天天线| 69av精品久久久久久| 最新在线观看一区二区三区| 国产精品人妻久久久久久| 男女啪啪激烈高潮av片| 国产精品国产三级国产av玫瑰| 久久韩国三级中文字幕| 中出人妻视频一区二区| 97超碰精品成人国产| 白带黄色成豆腐渣| 人人妻,人人澡人人爽秒播| 亚洲性夜色夜夜综合| 99热6这里只有精品| 国产探花在线观看一区二区| 国内少妇人妻偷人精品xxx网站| 五月玫瑰六月丁香| 最后的刺客免费高清国语| 在线播放国产精品三级| 一级黄色大片毛片| 91av网一区二区| 久久韩国三级中文字幕| 国产69精品久久久久777片| 欧美高清成人免费视频www| 99国产极品粉嫩在线观看| 午夜亚洲福利在线播放| 亚洲精品456在线播放app| 亚洲,欧美,日韩| 欧美日韩精品成人综合77777| 最新中文字幕久久久久| 丰满乱子伦码专区| 麻豆成人午夜福利视频| 国产在视频线在精品| 在线观看免费视频日本深夜| 搡女人真爽免费视频火全软件 | 少妇熟女欧美另类| 丝袜喷水一区| 国产精品一区二区性色av| 久久久久久久午夜电影| 简卡轻食公司| 亚洲成av人片在线播放无| 亚洲av中文av极速乱| 日韩成人av中文字幕在线观看 | 婷婷精品国产亚洲av在线| 日本爱情动作片www.在线观看 | 免费观看在线日韩| 国产黄a三级三级三级人| 国产成人aa在线观看| 久久亚洲国产成人精品v| 国产三级在线视频| 天美传媒精品一区二区| 免费看av在线观看网站| 大香蕉久久网| 国内揄拍国产精品人妻在线| av在线播放精品| 亚洲精品456在线播放app| 欧美人与善性xxx| av天堂在线播放| а√天堂www在线а√下载| 直男gayav资源| 一进一出抽搐gif免费好疼| 国产精品爽爽va在线观看网站| 中出人妻视频一区二区| 国产黄片美女视频| 国产精品国产三级国产av玫瑰| 久久久色成人| 欧美最黄视频在线播放免费| 日本五十路高清| 日日干狠狠操夜夜爽| 性欧美人与动物交配| 伦理电影大哥的女人| av免费在线看不卡| 久久久国产成人免费| 国产伦精品一区二区三区四那| 国产av在哪里看| 国产成人a∨麻豆精品| 亚洲av成人精品一区久久| ponron亚洲| 天美传媒精品一区二区| 亚洲不卡免费看| 国产免费一级a男人的天堂| 噜噜噜噜噜久久久久久91| 精品无人区乱码1区二区| 免费无遮挡裸体视频| 少妇的逼好多水| 免费在线观看成人毛片| 三级经典国产精品| 天堂√8在线中文| 99热全是精品| 一个人看视频在线观看www免费| 久久久久久久久大av| 久99久视频精品免费| 婷婷亚洲欧美| 国产欧美日韩精品亚洲av| 日本撒尿小便嘘嘘汇集6| 国内精品美女久久久久久| 亚洲av第一区精品v没综合| av黄色大香蕉| 悠悠久久av| 久久综合国产亚洲精品| 国产单亲对白刺激| 欧美日韩一区二区视频在线观看视频在线 | 美女xxoo啪啪120秒动态图| 人妻久久中文字幕网| 欧美日韩综合久久久久久| 国产精品,欧美在线| 最近手机中文字幕大全| 天堂av国产一区二区熟女人妻| 人人妻,人人澡人人爽秒播| 性插视频无遮挡在线免费观看| 久久人妻av系列| 亚洲av中文av极速乱| 日本爱情动作片www.在线观看 | 国产真实伦视频高清在线观看| 日韩成人av中文字幕在线观看 | 日韩欧美精品v在线| 搡老妇女老女人老熟妇| 简卡轻食公司| 亚洲国产日韩欧美精品在线观看| 婷婷色综合大香蕉| 亚洲欧美清纯卡通| 亚洲欧美精品综合久久99| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区三区av在线 | 国产精品三级大全| 亚洲精品456在线播放app| 亚洲av五月六月丁香网| 国产高清视频在线观看网站| 久久久久性生活片| 国产精品一区二区免费欧美| 日韩av不卡免费在线播放| 日日摸夜夜添夜夜爱| 久久国内精品自在自线图片| 高清毛片免费看| 一区二区三区四区激情视频 | 日韩亚洲欧美综合| 国产美女午夜福利| or卡值多少钱| 99久国产av精品国产电影| 午夜激情福利司机影院| 午夜免费激情av| 成年免费大片在线观看| 99热精品在线国产| av天堂中文字幕网| 中文亚洲av片在线观看爽| 亚洲人成网站高清观看| 网址你懂的国产日韩在线| 婷婷精品国产亚洲av在线| 国产高清视频在线观看网站| 天美传媒精品一区二区| 久久精品夜色国产| av在线亚洲专区| 免费看a级黄色片| 久久精品影院6| 国产黄色视频一区二区在线观看 | 亚洲精品一卡2卡三卡4卡5卡| 蜜臀久久99精品久久宅男| 乱人视频在线观看| 国产亚洲欧美98| 十八禁国产超污无遮挡网站| 久久精品国产亚洲av涩爱 | 亚洲人成网站在线播| 午夜久久久久精精品| 国产精品不卡视频一区二区| 久久热精品热| 亚洲熟妇熟女久久| 亚洲国产色片| 亚洲丝袜综合中文字幕| ponron亚洲| 长腿黑丝高跟| 美女 人体艺术 gogo| 99热精品在线国产| 国产精品综合久久久久久久免费| 国产精品久久久久久精品电影| 嫩草影院入口| 国产又黄又爽又无遮挡在线| 欧美激情久久久久久爽电影| 99国产极品粉嫩在线观看| 亚洲最大成人av| 亚洲欧美日韩高清专用| 丝袜喷水一区| 1000部很黄的大片| 亚洲熟妇熟女久久| 成人性生交大片免费视频hd| 中文字幕人妻熟人妻熟丝袜美| 午夜日韩欧美国产| 蜜臀久久99精品久久宅男| 一级毛片久久久久久久久女| 免费看日本二区| 少妇高潮的动态图| 午夜a级毛片| 国产精品久久久久久av不卡| 如何舔出高潮| 三级经典国产精品| 久久这里只有精品中国| 日日摸夜夜添夜夜爱| 老司机午夜福利在线观看视频| 欧美日韩国产亚洲二区| 久久国产乱子免费精品| 一本一本综合久久| 国产精品99久久久久久久久| 国产三级在线视频| 乱人视频在线观看| 免费看日本二区| 国语自产精品视频在线第100页| 男女做爰动态图高潮gif福利片| 欧美成人a在线观看| 日韩欧美国产在线观看| 两个人的视频大全免费| 国产成年人精品一区二区| 1024手机看黄色片| 欧美日本视频| www日本黄色视频网| 日韩高清综合在线| 黄色欧美视频在线观看| av天堂在线播放| 偷拍熟女少妇极品色| 亚洲一区高清亚洲精品| 内射极品少妇av片p| 日韩欧美一区二区三区在线观看| 免费在线观看成人毛片| 九九在线视频观看精品| 少妇裸体淫交视频免费看高清| 亚洲在线自拍视频| 日韩中字成人| 韩国av在线不卡| 日本-黄色视频高清免费观看| 久久韩国三级中文字幕| 天堂av国产一区二区熟女人妻| 一区福利在线观看| 桃色一区二区三区在线观看| 日本黄大片高清| 六月丁香七月| 日韩国内少妇激情av| 久久久精品94久久精品| 男女那种视频在线观看| 国产一区二区在线观看日韩| 午夜激情欧美在线| 成人特级黄色片久久久久久久| 一区二区三区免费毛片| 日韩国内少妇激情av| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲网站| 国产乱人偷精品视频| 青春草视频在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 国产精品久久久久久久久免| 性欧美人与动物交配| 男人狂女人下面高潮的视频| 中文字幕免费在线视频6| 性欧美人与动物交配| 成人av在线播放网站| 能在线免费观看的黄片| av女优亚洲男人天堂| 亚洲色图av天堂| 亚洲性夜色夜夜综合| 啦啦啦观看免费观看视频高清| 日日撸夜夜添| 欧洲精品卡2卡3卡4卡5卡区| 午夜福利在线在线| 欧美在线一区亚洲| 日韩欧美国产在线观看| 久久精品夜夜夜夜夜久久蜜豆| 国产精品伦人一区二区| 国产精品福利在线免费观看| 老司机影院成人| 精品一区二区三区视频在线观看免费| 91久久精品国产一区二区成人| 中文字幕av在线有码专区| 亚洲18禁久久av| 十八禁国产超污无遮挡网站| 国产成人91sexporn| 99热这里只有精品一区| 亚洲成人中文字幕在线播放| 麻豆精品久久久久久蜜桃| 此物有八面人人有两片| 精品无人区乱码1区二区| 91午夜精品亚洲一区二区三区| 非洲黑人性xxxx精品又粗又长| 欧美一区二区精品小视频在线| 别揉我奶头~嗯~啊~动态视频| 国产一级毛片七仙女欲春2| 国产精华一区二区三区| 成人亚洲欧美一区二区av| 欧美zozozo另类| 伊人久久精品亚洲午夜| 亚洲av第一区精品v没综合| 亚洲精华国产精华液的使用体验 | 热99在线观看视频| 亚洲av不卡在线观看| 免费在线观看影片大全网站| 欧美色视频一区免费| 成人特级黄色片久久久久久久| 美女cb高潮喷水在线观看| 在线观看66精品国产| 97热精品久久久久久| av天堂中文字幕网| 欧美xxxx黑人xx丫x性爽| 一级黄色大片毛片| 精品久久久久久久末码| 亚洲av二区三区四区| videossex国产| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 九九在线视频观看精品| 日韩欧美精品免费久久| 全区人妻精品视频| 男人的好看免费观看在线视频| 国产熟女欧美一区二区| 亚洲婷婷狠狠爱综合网| 老熟妇乱子伦视频在线观看| 久久6这里有精品| 网址你懂的国产日韩在线| 亚洲精品影视一区二区三区av| 三级经典国产精品| 亚洲欧美清纯卡通| 欧美三级亚洲精品| 亚洲国产精品成人久久小说 | 精品人妻偷拍中文字幕| 在线免费观看不下载黄p国产| 日韩精品青青久久久久久| 99热这里只有是精品在线观看| 国产午夜精品久久久久久一区二区三区 | www.色视频.com| 老女人水多毛片| АⅤ资源中文在线天堂|