999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

金昌市≥10 ℃積溫變化特征及對農作物的影響

2025-03-11 00:00:00劉騰嬌劉小瑛張軍鵬陳文博黃建文
甘肅農業科技 2025年2期
關鍵詞:趨勢

摘要:通過探究金昌市氣候變化背景下活動積溫(≥10 ℃)對農作物生長的影響,因地制宜地安排作物布局、合理部署農業生產。選用金昌市金川區、永昌縣1991 — 2020年30 a的觀測數據,運用趨勢分析法、M-K突變檢驗法、五日滑動平均法等氣候統計診斷方法對金昌市≥10 ℃積溫時空變化特征及玉米、春小麥的種植面積、產量進行分析。結果表明,氣候變暖背景下,金昌市熱量條件有較好的改善,對農作物產生了一定影響,喜溫作物玉米的種植比例和產量顯著增加,喜涼作物春小麥的種植比例減少。金昌市≥10 ℃積溫、積溫持續日數總體呈波動增加的趨勢,但金川區≥10 ℃積溫持續日數的增加主要由初日提前造成,永昌縣由終日延后造成。金昌市≥10 ℃積溫異常概率較小,偏高年均在積溫突變之后,金川區積溫突變發生在1996年,永昌縣發生在1997年。金昌市≥10 ℃積溫與海拔密切相關,海拔每升高100 m,≥10 ℃積溫減少187.5 ℃·d、終日提早2.4 d、初日推遲2.4 d、持續日數延長5.1 d。

關鍵詞:趨勢分析;M-K突變檢驗;五日滑動平均;均方差判別;≥10 ℃積溫

中圖分類號:P468" " " " " " " " 文獻標志碼:A" " " " " " " " 文章編號:2097-2172(2025)02-0159-07

doi:10.3969/j.issn.2097-2172.2025.02.011

Variation Characteristics of ≥10 ℃ Accumulated Temperature and

Its Effect on Crops in Jinchang

LIU Tengjiao, LIU Xiaoying, ZHANG Junpeng, CHEN Wenbo, HUANG Jianwen

(Jinchang Meteorological Bureau, Jinchang Gansu 737100, China)

Abstract: This study explored the impact of the ≥10 ℃ accumulated temperature values on crop growth under the background of climate change in Jinchang City, aiming to arrange crop layouts and plan agricultural production appropriately. Using 30 years of observational data(1991 to 2020) from Jinchuan District and Yongchang County, Jinchang City's, the study analyzed the spatial and temporal variation characteristics of ≥10 ℃ accumulated temperature and the planting area and yields of maize and spring wheat using climate statistical methods such as trend analysis, M-K mutation test, and five-day moving average. Results showed that under the background of cilmate warming, Jinchang's thermal conditions were improved significantly, affecting crops. The planting proportion and yield of warm-loving crops like maize were significantly increased, while the planting proportion of cool-loving crops like spring wheat were decreased. The ≥10 ℃ accumulated temperature and the number of days with accumulated temperature in Jinchang City showed an overall fluctuating increase. However, the increase in the number of days with accumulated temperature in Jinchuan District was mainly caused by an earlier start, while in Yongchang County, it was caused by a later end. The probability of abnormal ≥10 ℃ accumulated temperature in Jinchang City was relatively low, with higher-than-usual years occurring after a temperature mutation. The temperature mutation in Jinchuan District occurred in 1996, and in Yongchang County, it occurred in 1997. The ≥10 ℃ accumulated temperature in Jinchang City was closely related to altitude. For every 100 meters increase in altitude, the ≥10 ℃ accumulated temperature decreased by 187.5 ℃·d, the end day advanced by 2.4 days, the start day was delayed by 2.4 days, and the duration of the accumulated temperature increased by 5.1 days.

Key words: Trend analysis; M-K mutation test; Five-day moving average; Mean square error discrimination; ≥10 ℃ accumulated temperature

氣候變化是全球關注的熱點問題,聯合國氣候變化政府間專家委員會指出,1880 — 2012年全球地表平均溫度大約增加了0.85 ℃,氣候增暖會引起熱量資源的變化,對農業生產將造成一定的影響[1 ]。就甘肅省而言,氣候變暖帶來的積溫增加及積溫帶北移東擴使得進入20世紀90年代后各地作物種植結構均發生很大的變化,具體表現在河西地區玉米和棉花種植比例明顯增加,而小麥種植比例明顯減少,從主要以小麥為主的糧食作物種植結構演化為以玉米和棉花為主的格局[2 - 3 ]。基于平均氣溫和活動積溫(≥10 ℃)分析農作物熱量資源變化,探究氣候變化背景下活動積溫的變化值對農作物生長的影響,因地制宜地安排作物布局、合理部署農業生產非常重要[4 - 5 ]。

金昌市別稱鎳都,地處中國西北地區、甘肅省河西走廊中段、祁連山北麓,地勢南高北低,總面積96萬hm2,主要種植農作物為玉米和小麥。通過對金昌市≥10 ℃積溫時空變化特征分析,揭示了氣候變化背景下熱量資源新格局及其對主要農作物的影響,將對本地現代農業結構規劃、農作物品種調整以及農業的定量化評估具有重要意義。

1" "資料來源與方法

1.1" "資料來源

資料來源于甘肅省金昌氣象觀測站、永昌縣觀測站1991 — 2020年的氣象觀測資料,氣候基準值采用1991 — 2020年的平均值。觀測數據均來自地面氣象記錄年報表,相關天氣日數的定義和統計方法參照《地面氣象觀測規范》[6 ],農作物的統計來自《永昌縣統計年鑒》《金川區統計年鑒》。

1.2" "研究方法

1.2.1" " Mann-Kendall突變檢驗法" " M-K檢驗法是一種非參數統計檢驗方法,與參數統計檢驗方法相比,該方法不需要樣本遵從一定的分布規律,也不受少數異常值的干擾,且計算也相對比較簡單,是目前比較常用的趨勢診斷方法[7 ]。用M-K檢驗法對金昌熱量要素突變特征進行分析。

1.2.2" " 線性趨勢分析法" " 用一元線性方程描述線性趨勢變化,建立氣候變量與其所對應時間的一元線性回歸方程。

xi=a+bti(i=1,2,…,n),

式中,xi為≥10 ℃年積溫,a為回歸常數,b為回歸系數,a和b可以用最小二乘進行估計。b的符號表示氣候變量x的趨勢傾向,b值的大小反映上升或下降的速率,ti為時間[8 ]。

1.2.3" " 均方差判別法" " 運用≥10 ℃積溫多年平均值(x)和均方差(δ)判斷指標對xi的異常性進行判斷,若xi在(x ±δ)之內為正常年;xi在(x±δ)與(x±2δ)之間為偏高或偏低年;xi在(x±2δ)之外為特高或特低年[9 ]。

1.2.4" " 五日滑動平均法" " 采用五日滑動平均法計算日均氣溫,通過≥10 ℃界限溫度起止日期,并通過初日和終日計算≥10 ℃積溫和持續日數[10 ],五日滑動平均法具體過程為:在時長為1 a的日平均溫度序列中,選取任意連續5天日平均溫度大于或等于某一界限溫度最長的一段時期,在第1個5天中,選取第1個日平均溫度大于或等于某一界限溫度的日期,作為穩定通過某一界限溫度的初日,在最后1個5天中,選取最后1個日平均溫度大于或等于某一界限溫度的日期,作為穩定通過某一界限溫度的終日,初日和終日中間的溫度之和即為穩定通過某一界限溫度的積溫[11 - 12 ]。

2" "結果與分析

2.1" "金昌市≥10 ℃積溫年際變化特征

2.1.1" " 初日" " 金川區≥10 ℃積溫初日多年平均在4月9日,最早為3月21日(2013年),最晚為4月24日(2008年);永昌縣初日平均在4月23日,最早為4月4日(2008年),最晚為5月12日(2014年)。兩地的初日總體呈提前的趨勢,但金川初日呈顯著提前趨勢,變化速率為4.3 d/10 a,永昌縣僅為0.3 d/10 a(圖1 a)。

2.1.2" " 終日" " 金川區≥10 ℃積溫終日多年平均在10月8日,最早為9月16日(2011年),最晚為10月22日(2009、2010年);永昌縣終日平均在9月25日,最早為9月13日(2002、2018年),最遲10月21日(2008年)。金川區≥10 ℃積溫終日變化幅度不大,呈微弱延后的趨勢,變化速率為0.9 d/10 a,近30 a提早2.7 d;永昌縣≥10 ℃積溫終日變化速率為3.3 d/10 a,近30 a延后9.9 d,呈顯著延后的趨勢(圖1 b)。

2.1.3" " 持續日數" " 金川區≥10 ℃積溫持續日數平均為176 d,最少為157 d(2011年),最多為193 d(2009年),相差36 d;永昌縣平均147 d,最少為120 d(2002年),最多為164 d(2009年)。金川區、永昌縣≥10 ℃的持續日數均表現出增加趨勢,變化速率分別為0.37、0.44 d/10 a (圖1 c)。

2.1.4" " ≥10 ℃積溫" " 金川區≥10 ℃積溫平均為" "3 483.7 ℃·d,最小值為3 217.0 ℃·d(1995年),最大值為3 768.5 ℃·d(2013年),兩者相差493.6 ℃·d;永昌縣≥10 ℃積溫平均為2 386.4 ℃·d,最小值為2 132.2 ℃·d(2002年),最大值為2 632.1 ℃·d(2016年),兩者相差499.9 ℃·d;金昌市≥10 ℃積溫總體呈波動上升的趨勢,金川區增加速率為42.2 ℃·d/10 a,近30 a增加1 266 ℃·d,永昌縣增加速率為116.7 ℃·d/10 a,近30 a增加3 501 ℃·d(圖1 d)。

2.2" "金昌市≥10 ℃積溫年代際變化特征

2.2.1" " 初日" " 20世紀90年代至21世紀00年代金川區和永昌縣的變化并不明顯;21世紀00年代至10年代金川區平均提前了10 d,永昌縣略有延遲(圖2 a)。

2.2.2" " 終日" " 20世紀90年代至21世紀00年代金川區基本無變化,永昌縣延后6 d;21世紀00年代至10年代金川區、永昌縣均無變化(圖2 b)。

2.2.3" " 持續日數" " 各年代金昌市均表現出增加的趨勢,20世紀90年代至21世紀00年代變化更為明顯,21世紀00年代至10年代進入平穩期(圖2 c)。

2.2.4" " ≥10 ℃積溫" " 從圖2 d可知,金川區20世紀90年代≥10 ℃積溫平均值為3 439.5 ℃·d;21世紀00年代≥10 ℃平均積溫為3 507.4 ℃·d,較20世紀90年代上升了67.9 ℃·d;21世紀10年代≥10 ℃平均積溫為3 504.4 ℃·d,較21世紀00年代持平,較20世紀90年代增加64.9 ℃·d。永昌縣20世紀90年代≥10 ℃平均積溫為2 266.0 ℃·d,21世紀00年代≥10 ℃平均積溫為2 413.4 ℃·d,較20世紀90年代上升了147.4 ℃·d;21世紀10年代≥10 ℃平均積溫為2 479.7 ℃·d,較21世紀00年代上升66.3 ℃·d,較20世紀90年代增加213.7 ℃·d,增幅明顯。

2.3" "金昌市≥10 ℃積溫突變分析

2.3.1" " 金川區≥10℃積溫突變分析" " 從≥10 ℃積溫累積距平曲線(圖3 a)可以看出,金川區累積距平近30 a呈波動上升的趨勢,1991 — 1996年積溫累積距平均為負距平,且呈下降的趨勢,1996年以后以正距平為主,且呈上升的趨勢,表明1996年為≥10 ℃積溫的主要轉折年。對近30 a積溫進行M-K突變檢驗(圖3 b),正序列特征曲線(UF)與反序列特征曲線(UB)在1996年相交于信度之間,通過0.05顯著性檢驗,與累積距平變化轉折年份一致。≥10 ℃積溫突變后(1997 — 2020年)較突變前(1991 — 1996年)增加215.9 ℃·d。從≥10 ℃積溫持續日數累積距平變化曲線(圖4 a)可以看出,金川區1991 — 1996年持續日數以負距平為主,呈下降趨勢,1996年之后累積距平總體呈上升趨勢,且以正距平為主,表明1996年為≥10 ℃積溫持續日數的轉折年。突變后(1997 — 2020年)較突變前(1991 — 1996年)增加10.9 d。

2.3.2" " 永昌縣≥10 ℃積溫突變分析" " 運用同樣的方法對永昌縣≥10 ℃積溫進行M-K突變分析結果可知,永昌縣≥10 ℃積溫突變出現在1997年,突變后(1998 — 2020年)較突變前(1991 — 1997年)增加225.6 ℃·d;積溫持續日數突變后較突變前增加8.3 d。

2.4" "金昌市≥10 ℃積溫異常分析

對金川區、永昌縣近30 a≥10 ℃積溫進行均方差判別,并統計積溫異常年結果(表1)表明,金川區≥10 ℃積溫正常年為18 a,概率為60.0%;偏低年5 a,概率為16.7%;偏高年7 a,概率為23.3%。對永昌縣近30 a≥10 ℃積溫進行均方差判別,并統計積溫異常年,結果表明,≥10 ℃積溫正常年為22 a,概率為73.3%;偏低年2 a,概率為0.7%;偏高年5 a,概率為16.7%;特少年僅有1 a(2002年)。積溫異常年并沒有明顯的年際及年代際的變化特征,偏高年均發生在積溫突變之后。

2.5" "積溫與海拔的關系

由表2可看出,積溫、持續日數隨海拔的升高遞減,初日隨海拔的升高推遲,終日則相反,隨海拔升高提早。金昌市海拔每升高100 m,≥10 ℃積溫減少187.6 ℃·d;終日提早2.4 d;初日推遲2.4 d;持續日數延長5.1 d。殷雪蓮等[13 ]的研究表明,≥10 ℃積溫與海拔密切相關。

3" "積溫變化對主要作物的影響

3.1" "積溫變化對主要作物種植面積的影響

氣候變暖使作物種植格局發生較大演變[14 ],調查發現,1990 — 2005年金昌市的主要種植作物結構發生了改變,金川由小麥變為小麥、玉米;永昌由小麥變成了小麥、油菜[12 ]。金昌市種植的主要農作物為春小麥和玉米,且80%的作物在永昌縣種植。用春小麥、玉米作為喜涼作物、喜溫作物的代表進行分析,1991 — 2020年金昌市玉米的種植面積呈顯著性增加趨勢,金川區、永昌縣的增加速率分別為0.36萬、0.62萬hm2/10 a,但金川區春小麥種植面積略有減少。永昌縣主要農作物的種植面積均呈增加的趨勢,玉米種植比例顯著增大(圖5)。說明在氣候變暖背景下,農作物的種植結構發生了調整,喜溫的玉米種植比例迅速擴大。

3.2" "積溫變化對主要農作物產量的影響

從圖6可知,主要農作物的產量均呈增加趨勢,玉米增產更為明顯,永昌縣速率為4.95 t/10 a,金川區為2.93 t/10 a;金川區小麥的產量基本穩定。說明了氣候變暖,積溫的增加會對農作物產量產生一定影響,對發展喜溫的高產作物非常有利,這與曹玲等[15 ]的研究結果一致。金昌市3月氣溫近10 a異常偏高,春季氣溫回升較快(金昌市氣象局提供),使得玉米適播期提前,因此,應適時早播,金川區的積溫初日顯著提前,永昌終日延后,金川區的玉米應考慮更早播種,永昌縣可根據長短期氣候預測,適當提早播種,以提高農作物產量。近30 a金昌市春小麥的播種面積和產量的增加速率均不明顯,應適當減少種植比例,或者先在金川區嘗試種植冬小麥。

4" "討論與結論

受全球氣候變化影響,近年來極端天氣比常年更頻發,給農業生產造成了極大的威脅和損" " 失[4, 16 - 17 ]。通常以日平均氣溫穩定通過≥10 ℃的持續天數作為劃分氣候帶的主要指標,反映某區喜溫作物氣候生長期的長短,是決定喜溫作物生長和分布的重要指標之一。甘肅省從南到北分為北亞熱帶和溫帶,溫帶再細分為暖溫帶和中溫帶。在氣候變暖背景下,氣候帶的分布呈現出整體向高海拔擴張和高緯度北移的趨勢[3 ]。按照中國氣候區劃新方案對不同氣候帶的劃分指標,結合≥10 ℃積溫研究發現金川區屬于暖溫帶干旱區,永昌屬于中溫帶半干旱區[18 ],若是氣候持續增暖,金昌市的暖溫帶地區將擴大。張強等[19 ]研究表明,氣候變暖對越冬作物的冬前生長發育及喜溫作物和喜熱作物的全生育期生長發育均比較有利。根據我國氣象專家預測,21世紀我國氣候將明顯繼續變暖[20 - 22 ],氣候變暖將使農作物生長活動10 ℃以上的積溫增加,對發展喜溫的高產作物非常有利。因此,可以進一步擴大玉米種植面積和種植比例,提高糧食總產水平。雖然未來氣候將呈持續變暖趨勢,但仍會有相對低溫年份。由于玉米對熱量資源的要求較高[23 ],應在低溫年型降低玉米種植比例,增暖年型在原有基礎上加大種植面積。

氣候變暖使得金昌市熱量資源改善,對農作物產生了一定影響,主要影響了玉米種植比例和產量。永昌縣農作物種植面積和產量均增加,金川區小麥的種植比例和產量增幅很小,說明永昌縣的氣候對農作物的種植更為適宜。從提升經濟效益和優化作物種植的角度考慮,金昌應減少春小麥的種植面積,加大玉米的種植面積,或者嘗試冬小麥的種植;應適時提早播種,以提高玉米產量。氣候變暖雖是不爭的事實,但還是會出現相對低溫時段,應根據不同的氣候年適當調整玉米的種植結構和種植比例。

金昌市1991 — 2020年≥10 ℃積溫總體呈波動上升的趨勢,永昌縣的增加速率是金川區的2.8倍;≥10 ℃積溫持續日數呈增加的態勢,但金川區≥10 ℃積溫持續日數的增加主要由初日提前引起,永昌縣則由終日延后造成。金川區≥10 ℃積溫平均為3 483.7 ℃·d,永昌縣為2 386.4 ℃·d。≥10 ℃積溫突變金川發生在1996年,永昌發生在1997年,表現出局地特征。金川區和永昌縣≥10 ℃積溫年代際變化并不完全一致,20世紀90年代均表現出顯著增加的趨勢,21世紀10年代下降,21世紀00年代金川區增加緩慢,永昌仍顯著增加。金川區積溫增加20世紀90年代貢獻最大,永昌縣20世紀90年代至21世紀00年代均有顯著貢獻。金川區與永昌縣10 ℃積溫正常年最多,概率為60%以上,偏低年和偏高年概率小,金川區近30 a并無特高和特低年,永昌縣2002年為特低年,偏高年均發生在積溫突變之后。≥10 ℃積溫與海拔密切相關。積溫大小、持續日數隨海拔的升高遞減,初日隨海拔的升高推遲,終日則相反,隨海拔升高提早。金昌市≥10 ℃積溫隨海拔每升高100 m,≥10 ℃積溫減少187.5 ℃·d;終日提早2.4 d;初日推遲2.4 d;持續日數延長5.1 d。隨著氣候持續變暖,金昌市的氣候區可能會發生變化,具體表現為暖溫帶區域擴大。

本文根據近30 a積溫變化及農作物種植面積、產量,對規劃不同地域種植結構、調整玉米種植比例給出了針對性建議,但針對不同區域定量的種植規劃還需做進一步的研究。農業生產活動受多方因素的影響,氣候變化只是影響因素之一,此外還與農業技術、勞動力、市場和政策等密切相關。以現有的收集資料,無法單純研究氣候要素變化對農業的影響,只能初步分析氣候變化對金昌農作物產量的影響,并不能科學分析其對農業實際產量的客觀影響規律。

參考文獻:

[1] IPCC. Climate change 2013: The physical science basis [M]." Cambridge: Cambridge University Press, 2013.

[2] 王鶴齡,王潤元,張" "強,等." 甘肅省作物布局演變及其對區域氣候變暖的響應[J]." 自然資源學報,2012,27(3):413-421.

[3] 王國強,張" "勃,張耀宗,等." 甘肅省近55年來積溫變化趨勢特征分析[J]." 水土保持研究,2016,23(5):193-198.

[4] 李春艷,張" "慧,于永奎,等." 2022年平羅縣氣象因子對主要農作物生長發育的影響[J]." 寒旱農業科學,2024,3(3):245-250.

[5] 張美蘭,郭世乾,賈蕊鴻,等." 甘肅省馬鈴薯種植適宜性評價及影響因素分析[J]." 寒旱農業科學,2023,2(8):731-735.

[6] 中國氣象局." 地面氣象觀測規范[M]." 北京:氣象出版社,2003.

[7] 魏風英." 現代氣候診斷與預測技術[M]." 北京:氣象出版社,1999.

[8] 閆" "俊,王海功,李紅梅." 氣候變暖背景下安徽氣溫變化趨勢分析[J]." 安徽農業科學,2013,41(36):13953-13955;13971.

[9] 胡啟瑞,谷" "然,李新建,等." 1987 — 2020年石河子棉區≥10 ℃積溫變化特征分析[J]." 沙漠綠洲與氣象,2023,17(1):171-175.

[10] 張" "茹,杜金玲,馬秀枝,等. 大興安嶺林區1970 — 2019年積溫的時空變化特征[J]." 中國農學通報,2023,39(19):93-101.

[11] 孟艷靈,殷淑燕,楊" "鋒,等." 晉陜蒙地區≥10 ℃積溫的時空變化特征[J]." 中國農業氣象,2016,37(6):615-622.

[12] 孟艷靈." 1960 — 2015 年晉陜蒙地區農業水熱氣候條件時空變化研究[D]." 西安:陜西師范大學,2017.

[13] 殷雪蓮,何金梅,郭萍萍." 河西走廊≥10 ℃界限溫度演變特征及其對玉米生產的影響[J]." 干旱地區農業研究,2014,32(6):236-243.

[14] 鄧振鏞,張" "強,劉德祥,等." 氣候變暖對甘肅種植業結構和農作物生長的影響[J]." 中國沙漠,2007,27(4):627-632.

[15] 曹" "玲,鄧振鏞,竇永祥,等." 氣候變暖對河西走廊綠洲灌區玉米產量影響及對策研究[J]." 西北植物學報,2008(5):1043-1048.

[16] 張" "新,謝翠婷,周慧僚." 淺析天峨縣2021年氣象條件對農業生產的影響[J]." 農業災害研究,2022,12(8):94-96.

[17] 馬" "銳,江" "敏,石春林." 氣候變化對農業影響評價研究進展[J]." 江蘇農業科學,2015,43(10):16-19.

[18] 鄭景云,尹云鶴,李炳元." 中國氣候區劃新方案[J]." 地理學報,2010,65(1):3-12.

[19] 張" "強,鄧振鏞,趙映東,等." 全球氣候變化對我國西北地區農業的影響[J]." 生態學報,2008,28(3):1210-1218.

[20] 韓林君,白愛娟,蒲學敏." 基于 CMIP6 的祁連山氣候變化特征預估[J]." 高原氣象,2022,41(4):864-875.

[21] 何雨芩,范立張,楊曉鵬,等." 云南省甘蔗凍害風險分布及未來趨勢預估[J]." 氣象與環境科學,2023,46(2):58-66.

[22] 蘇正娥,劉志娟,楊婉蓉,等." 未來氣候變化情景下東北地區中晚熟春玉米機械粒收氣候適宜區分析[J]." 中國農業氣象,2023,44(8):649-663.

[23] 韓耀杰,彭記永,張溪荷,等." 氣候變化背景下河北夏玉米生長季熱量資源利用率[J]." 干旱氣象,2023,41(4):639-647.

猜你喜歡
趨勢
趨勢
第一財經(2025年5期)2025-05-16 00:00:00
退休的未來趨勢
英語世界(2023年12期)2023-12-28 03:36:16
趨勢
第一財經(2021年6期)2021-06-10 13:19:08
趨勢
初秋唇妝趨勢
Coco薇(2017年9期)2017-09-07 21:23:49
未來直銷的七大趨勢
趨勢
流行色(2016年10期)2016-12-05 02:27:24
SPINEXPO?2017春夏流行趨勢
關注醫改新趨勢
中國衛生(2015年2期)2015-11-12 13:14:02
“去編”大趨勢
中國衛生(2015年7期)2015-11-08 11:09:38
主站蜘蛛池模板: 波多野结衣一区二区三视频| 久久久精品国产SM调教网站| 乱人伦视频中文字幕在线| 国产婬乱a一级毛片多女| 日本高清有码人妻| 国产精品网曝门免费视频| 超清无码熟妇人妻AV在线绿巨人| 国产资源免费观看| 久久公开视频| 国产精品九九视频| 秋霞一区二区三区| 国产精品自在线天天看片| 久久青草视频| 免费国产福利| 国产乱码精品一区二区三区中文| 欧美精品综合视频一区二区| 欧美午夜网| 四虎精品国产AV二区| 男人天堂亚洲天堂| 亚洲av无码片一区二区三区| 亚洲色大成网站www国产| 97精品伊人久久大香线蕉| 日韩国产一区二区三区无码| 美女潮喷出白浆在线观看视频| 精品国产污污免费网站| 国产亚卅精品无码| 欧美性久久久久| 亚洲日韩Av中文字幕无码| 国产一区二区三区夜色| 国产美女精品一区二区| 久久中文电影| 日本黄网在线观看| AV无码无在线观看免费| 国产黄在线免费观看| 欧美午夜性视频| 天天色天天操综合网| 国产精品99在线观看| 久热re国产手机在线观看| 国产免费观看av大片的网站| 国产波多野结衣中文在线播放| 波多野结衣视频一区二区| 欧美在线观看不卡| 国产女人在线| 久久久久无码精品| 国产精品福利导航| 国产免费羞羞视频| 99ri国产在线| 毛片免费视频| 国产激情在线视频| 熟女成人国产精品视频| 五月天丁香婷婷综合久久| 国产成人AV综合久久| 亚洲福利一区二区三区| 伊人福利视频| 亚洲三级a| 99无码熟妇丰满人妻啪啪| 538国产在线| 亚洲天堂网视频| 国产一级在线观看www色 | 国产无码高清视频不卡| 免费人成网站在线高清| 国产色婷婷视频在线观看| 国产精品视屏| 国产在线视频自拍| 成人午夜久久| 青青草原偷拍视频| 国产成人精品在线| 亚洲乱码在线播放| 亚洲日韩在线满18点击进入| 亚洲精品第一页不卡| 国产精品永久在线| 四虎AV麻豆| 国产av一码二码三码无码| 91无码人妻精品一区二区蜜桃| 乱人伦视频中文字幕在线| 午夜欧美理论2019理论| 久久国产精品电影| 亚洲不卡网| 亚洲精品你懂的| 高清欧美性猛交XXXX黑人猛交| 日韩欧美国产精品| 欧美特级AAAAAA视频免费观看|