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

    長江三角洲地區(qū)極端氣溫事件變化特征及其與ENSO的關(guān)系

    2017-11-02 02:01:18賈艷青張耀宗王國強(qiáng)
    生態(tài)學(xué)報 2017年19期
    關(guān)鍵詞:長江三角洲地區(qū)拉尼厄爾尼諾

    賈艷青, 張 勃,張耀宗,唐 敏,馬 彬,王國強(qiáng)

    1 西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,蘭州 730070 2 忻州師范學(xué)院地理系,忻州 034000

    長江三角洲地區(qū)極端氣溫事件變化特征及其與ENSO的關(guān)系

    賈艷青1,2, 張 勃1,*,張耀宗1,唐 敏1,馬 彬1,王國強(qiáng)1

    1 西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院,蘭州 730070 2 忻州師范學(xué)院地理系,忻州 034000

    基于1960—2014年65個氣象站點(diǎn)逐日最高、最低氣溫和平均氣溫資料,分析了長江三角洲地區(qū)極端氣溫事件的變化規(guī)律和ENSO事件強(qiáng)度對極端氣溫指數(shù)變化趨勢的影響。結(jié)果表明:近55年長江三角洲地區(qū)夏季日數(shù)(SU)、熱夜日數(shù)(TR)、暖晝?nèi)諗?shù)(TX90)、暖夜日數(shù)(TN90)、異常暖晝持續(xù)指數(shù)(WSDI)、生長期(GSL)均呈增加趨勢,其中暖夜日數(shù)(TN90)增加幅度最大,增幅為8.55d/10a;極值指數(shù)也呈上升趨勢,其中月最低氣溫極小值(TNn)上升幅度最大為(0.53℃/10a);冰凍日數(shù)(ID)、霜凍日數(shù)(FD)、冷晝?nèi)諗?shù)(TX10)、冷夜日數(shù)(TN10)、異常冷晝持續(xù)指數(shù)(CSDI)均呈減少趨勢,其中冷夜日數(shù)(TN10)減少幅度最大(-6.06d/10a);月平均日較差(DTR)以0.11℃/10a的速率呈下降趨勢??臻g上,所有站點(diǎn)SU、TXn、TNx呈增加趨勢;TR、TX90、TN90、TNn、TXx、WSDI、GSL分別有97%、85%、98%、95%、78%、92%、94%的站點(diǎn)呈增加趨勢;所有站點(diǎn)ID、FD、TX10、TN10呈減少趨勢;CSDI、DTR分別有87%、77%的站點(diǎn)呈減少趨勢。多數(shù)極端氣溫指數(shù)與緯度、經(jīng)度、海拔顯著相關(guān)。氣候變暖突變后,極端暖指數(shù)明顯增加,極端冷指數(shù)明顯減少??傮w上,厄爾尼諾對極端氣溫指數(shù)的影響大于拉尼娜的影響。

    極端氣溫指數(shù);空間分布;氣候變化;ENSO

    近100年,全球地表溫度升高了0.85℃[1],近50年中國年平均地面溫度以0.22℃/10a的速率上升了1.1℃,增幅比同期全球或北半球高很多[2]。全球氣溫升高不僅導(dǎo)致極端氣溫發(fā)生變化,而且引起極端氣候事件頻率和強(qiáng)度的變化。近年來,中國群發(fā)性或區(qū)域性的極端氣候事件頻次和范圍不斷增大[3]。長江三角洲地區(qū)地處中國東部,由于地勢低洼,是中國自然災(zāi)害高風(fēng)險地區(qū),尤其臺風(fēng)、洪澇、熱浪等極端氣候災(zāi)害頻繁,對當(dāng)?shù)亟?jīng)濟(jì)發(fā)展和人民生活造成不利影響。近年來不同學(xué)者從極端高溫、低溫事件發(fā)生頻率、強(qiáng)度[4],極端高、低溫災(zāi)害特征[5],極端氣溫事件發(fā)生原因[6]等不同角度對長江三角洲地區(qū)極端氣溫事件開展了研究,但這些研究僅關(guān)注氣溫極值態(tài)下的某一類極端氣溫事件,綜合應(yīng)用多指標(biāo)對長江三角洲地區(qū)極端氣溫事件的研究較少且研究時間序列截止年份也較早。本研究采用更長時間序列資料和最新的極端氣溫指數(shù)研究極端氣溫變化特征,以期進(jìn)一步認(rèn)識長江三角洲地區(qū)極端天氣氣候事件變化規(guī)律,為構(gòu)建自然災(zāi)害防御應(yīng)急機(jī)制和氣候安全體系提供基礎(chǔ)資料。

    圖1 長江三角洲地區(qū)氣象站點(diǎn)分布圖 Fig.1 Distribution of meteorological stations in the Yangtze River Delta

    研究表明,厄爾尼諾/南方濤動(ENSO)是海洋與大氣環(huán)流耦合過程中的強(qiáng)烈信號,常引起大氣環(huán)流的改變并對全球范圍的氣候異常都有明顯的影響[7- 9],同時許多觀測也證實了赤道中東太平洋海表溫度異常與長江中下游地區(qū)氣溫升高有很好的對應(yīng)關(guān)系。近年來,在ENSO對我國冬季氣溫變化的影響方面取得了大量的研究成果[10-14],但長江三角洲地區(qū)極端氣溫與ENSO之間關(guān)系的研究未見報道。因此,本文分析ENSO對長江三角洲地區(qū)極端氣溫變化趨勢的影響,為深入認(rèn)識長江三角洲地區(qū)極端天氣事件的變化特征和成因提供科學(xué)依據(jù)。

    1 數(shù)據(jù)與方法

    1.1 數(shù)據(jù)來源

    氣象數(shù)據(jù)來源于中國氣象局國家氣象信息中心,本文所指的長江三角洲地區(qū)在行政區(qū)劃上包括安徽省、江蘇省、浙江省和上海市。選取了該區(qū)70個氣象站1960—2014年逐日平均氣溫、最高氣溫、最低氣溫資料,依據(jù)每月累計缺測不超過10%,每年累計缺測不超過5%的標(biāo)準(zhǔn)剔除了不符合要求的數(shù)據(jù),共選取了65個氣象站(圖1)。記錄顯示,自1960年以來,碭山、阜寧、南通、義烏、銅陵、射陽、如皋、金華8站均發(fā)生過遷站。由于臺站遷移和觀測方式發(fā)生改變等因素,使得這些站點(diǎn)的氣候序列不可避免地存在不均一性現(xiàn)象,因此研究中使用RHtest方法對上述8個臺站的數(shù)據(jù)進(jìn)行了均一化檢驗。ENSO數(shù)據(jù)來源于美國海洋大氣局(NOAA)氣候預(yù)測中心(CPC),時間段為1960—2014年,選用赤道中、東太平洋El Nio 3.4區(qū)(5°N—5°S,120°—160°W)的海洋表面溫度距平SSTA數(shù)據(jù)來表征ENSO事件。本文根據(jù)研究區(qū)氣候特點(diǎn)將長江三角洲地區(qū)劃分為3個子區(qū)域,Ⅰ區(qū)為南溫帶亞濕潤區(qū),包括安徽和江蘇北部,共分布12個氣象站點(diǎn);Ⅱ區(qū)為北亞熱帶濕潤區(qū),包括上海、安徽和江蘇中南部及浙江北中部,共分布40個氣象站點(diǎn);Ⅲ區(qū)為中亞熱帶濕潤區(qū),包括浙江中南部,共分布13個氣象站點(diǎn)。

    1.2 研究方法

    極端氣溫指數(shù)采用的是WMO氣候委員會與CLIVAR計劃聯(lián)合設(shè)立的氣候變化檢測、監(jiān)測和指標(biāo)專家組(ETCCDMI)確定的極端氣溫變化指標(biāo)(表1),使用Rclimdex軟件計算了長江三角洲地區(qū)65個氣象站點(diǎn)的16個極端氣溫指數(shù)。采用一元線性趨勢法分析極端氣溫指數(shù)的時間變化趨勢[15],并使用 Mann-Kendall(M-K)方法進(jìn)行變化趨勢的顯著性檢驗[16],空間插值采用反距離加權(quán)插值法,突變檢驗采用Mann-Kendall (M-K)法。運(yùn)用Morlet 連續(xù)小波分析方法[17],分析研究區(qū)近55a極端氣溫指數(shù)的周期變化特征。選用El Nio 3.4區(qū)海表溫度距平SSAT來表征ENSO事件,當(dāng)SSAT連續(xù)6個月高于0.5℃時為1次El Nio事件或ENSO暖事件,當(dāng)連續(xù)6個月低于-0.5℃時為1次La Nina事件或ENSO冷事件。根據(jù)李曉燕等[18]對ENSO暖(冷)事件劃分標(biāo)準(zhǔn)的研究,對ENSO事件強(qiáng)度等級進(jìn)行量化,將El Nio事件年強(qiáng)弱程度按照極強(qiáng)(5)、強(qiáng)(4)、中等(3)、弱(2)、極弱(1),La Nina事件年強(qiáng)弱程度按照極強(qiáng)(-5)、強(qiáng)(-4)、中等(-3)、弱(-2)、極弱(-1)的標(biāo)準(zhǔn)分別劃分為5類強(qiáng)度年,未受ENSO事件影響的年份定義為0。

    表1 極端氣溫指標(biāo)的定義

    2 結(jié)果分析

    2.1 長江三角洲地區(qū)極端氣溫指數(shù)時空變化

    2.1.1 年際變化趨勢

    長江三角洲地區(qū)16種極端氣溫指數(shù)的年際變化趨勢表明(表2),絕對指數(shù)ID與FD以-0.40d/10a和-4.03d/10a的速率顯著下降,SU和TR以3.29d/10a和3.14d/10a的速率顯著上升。相對指數(shù)TX10和TN10呈顯著減少趨勢,年際傾向率分別為-2.56d/10a和-6.06d/10a;TX90和TN90呈顯著增加趨勢,年際傾向率分別為4.27d/10a和8.55d/10a。TXn、TNn、TXx和TXn均表現(xiàn)為上升趨勢,其年際傾向率分別為0.35、0.53、0.17、0.27℃/10a。其它指數(shù)中CSDI和DTR呈減少趨勢,年際傾向率分別為-0.33 d/10a和-0.11℃/10a;WSDI和GSL呈增加趨勢,年際傾向率分別為1.34d/10a和3.11d/10a。以上極端氣溫指數(shù)的變化趨勢中,除TXx和CSDI僅通過0.05的顯著性檢驗外,其它指數(shù)均通過0.01或0.001的顯著性水平檢驗。

    3個氣候區(qū)的極端氣溫指數(shù)與全區(qū)變化趨勢一致,冷指日數(shù)呈上升趨勢,暖指日數(shù)呈下降趨勢,極值指數(shù)微弱上升。絕對指數(shù)ID和FD在Ⅰ區(qū)減少趨勢最明顯,年際傾向率分別為-0.84d/10a和-6.20 d/10a;SU在Ⅲ區(qū)增加趨勢最明顯(4.30 d/10a),TR在Ⅰ區(qū)增加趨勢最明顯(3.90 d/10a)。相對指數(shù)TX10在Ⅲ區(qū)減少趨勢最大(-3.52d/10a),在Ⅰ區(qū)減少趨勢最小。

    表2長江三角洲地區(qū)極端氣溫指數(shù)年際變化趨勢、顯著性檢驗及與其它區(qū)域?qū)Ρ?/p>

    Table2Comparisonbetweentheinter-annualvariationtrends,significancetestoftemperatureextremesoverYangtzeRiverDeltaandotherdifferentclimaticregions

    指數(shù)Indices長江三角洲地區(qū)YangtzeRiverDelta(1960—2014)全區(qū)WholeregionⅠ區(qū)ⅠregionⅡ區(qū)ⅡregionⅢ區(qū)ⅢregionKPKPKPKP中國China(1961—2008)K長江流域YangtzeRiverBasin(1962—2011)KID-0.400.010-0.840.016-0.360.013-0.090.033-2.32-0.48FD-4.03<0.0001-6.20<0.0001-3.85<0.0001-2.58<0.0001-3.48-3.29SU3.29<0.00011.550.051*3.48<0.00014.30<0.00011.182.93TR3.14<0.00013.90<0.00013.02<0.00012.82<0.00012.001.80TX10-2.560.003-1.870.055*-2.450.008-3.52<0.0001-0.47-0.84TN10-6.06<0.0001-7.31<0.0001-5.86<0.0001-5.50<0.0001-2.06-2.78TX904.27<0.001-0.220.832*4.910.00016.43<0.00010.622.24TN908.55<0.000111.15<0.00017.81<0.00018.41<0.00011.752.86TXn0.350.0050.340.0190.59<0.00010.360.0030.350.33TNn0.53<0.00010.750.00020.520.00010.360.00030.630.47TXx0.170.0320.000.968*0.200.0240.260.0030.070.16TNx0.27<0.00010.28<0.00010.26<0.00010.27<0.00010.210.19WSDI1.340.0040.400.333*1.520.0041.660.0003/0.83CSDI-0.330.015-0.300.072*-0.300.038-0.470.009/-0.67GSL3.110.0054.850.00022.880.0302.210.017/0.23DTR-0.110.0006-0.33<0.0001-0.080.0160.000.914*-0.18-0.07

    K表示年際傾向率,P表示顯著性水平,*未通過0.05顯著性檢驗

    TN10在Ⅰ區(qū)減少趨勢最大(-7.31d/10a);TX90在Ⅱ區(qū)和Ⅲ區(qū)均為增加趨勢,其中Ⅲ區(qū)增加趨勢最明顯(6.43 d/10a);TN90在Ⅰ區(qū)增加趨勢最大(11.15 d/10a)。極值指數(shù)除TXx在Ⅰ區(qū)沒有明顯變化外,其它指數(shù)在3個分區(qū)均為微弱增加趨勢。TXn在Ⅱ區(qū)上升趨勢最顯著(0.59℃/10a),在Ⅲ區(qū)次之;TNn在Ⅰ區(qū)上升趨勢最大(0.75℃/10a);TXx在Ⅲ區(qū)上升趨勢最大(0.26℃/10a);TNx在3個分區(qū)上升趨勢差異不大。其它指數(shù)中除DTR在Ⅲ區(qū)變化趨勢不明顯且未通過顯著性檢驗外,在各氣候區(qū)中,CSDI和DTR呈減少趨勢,WSDI和GSL呈增加趨勢。CSDI在Ⅲ區(qū)減少趨勢最顯著(-0.47 d/10a);DTR在Ⅰ區(qū)減少趨勢最顯著(-0.33℃/10a);WSDI和GSL分別在Ⅲ區(qū)和Ⅰ區(qū)增加趨勢顯著,年際傾向率分別為1.66和4.85 d/10a。

    對比分析長江三角洲地區(qū)與全國、長江流域大致相同時段極端氣溫指數(shù)的變化趨勢值發(fā)現(xiàn)(表2),長江三角洲地區(qū)暖指日數(shù)(SU、TR、TX90、TN90、WSDI、GSL)增加趨勢、冷指日數(shù)(ID、FD、TX10、TN10、CSDI)減少趨勢和多數(shù)極值指數(shù)的上升趨勢均比全國和長江流域明顯,表明長江三角洲地區(qū)氣溫上升趨勢要比全國和長江流域顯著,這與崔林麗、謝志清[4,19]等的研究結(jié)論一致,原因可能是長江三角洲地區(qū)人為熱排放量隨著城市化的快速發(fā)展不斷增加造成的。

    2.1.2 空間變化

    圖2是長江三角洲地區(qū)極端氣溫指數(shù)變化趨勢的空間分布。可以看出,除瑞安站不存在冰凍日數(shù)外,所有站點(diǎn)ID均處于減少趨勢,但變幅較小,多數(shù)站點(diǎn)的傾向率在-1.2—0.1d/10a之間,其中Ⅰ區(qū)減少趨勢比較明顯,平均為-0.82d/10a,贛榆減少速率達(dá)-1.2d/10a。所有站點(diǎn)FD均呈減少趨勢,Ⅰ區(qū)減少幅度較大,平均為-5.94 d/10a。所有站點(diǎn)SU均為增加趨勢,Ⅱ區(qū)和Ⅲ區(qū)站點(diǎn)增加趨勢顯著,尤其義烏、仙居、鄞縣站,平均增加幅度為6.22 d/10a。97%的站點(diǎn)TR為增加趨勢,長江三角洲北部和中部多數(shù)站點(diǎn)增加趨勢較明顯,平均為5.65 d/10a。

    圖2 長江三角洲地區(qū)各項極端氣溫指數(shù)變化趨勢空間分布Fig.2 Spatial distribution of inter-annual variation of temperature extremes in the Yangtze River Delta

    所有站點(diǎn)TX10均呈減少趨勢,長江三角洲東部部分站點(diǎn)減少趨勢比較明顯,平均為-3.94 d/10a。所有站點(diǎn)TN10均處于減少趨勢,減幅為-11.8—-0.66 d/10a,Ⅰ區(qū)和Ⅱ區(qū)中東部減少趨勢明顯,平均為7.62 d/10a。TX90有85%的站點(diǎn)呈增加趨勢,增加幅度為0.3—11.95 d/10a;Ⅱ區(qū)和Ⅲ區(qū)多數(shù)站點(diǎn)增加趨勢較明顯,平均為6.79 d/10a,其余15%的站點(diǎn)處于減少趨勢。TN90除義烏站(1.88 d/10a)趨于減少外,其余站點(diǎn)以0.20—23.62 d/10a的速率呈增加趨勢;Ⅰ區(qū)和Ⅱ區(qū)大部分站點(diǎn)增加趨勢顯著,平均為10.15 d/10a。

    所有站點(diǎn)TXn均呈上升趨勢,平均為0.38℃/10a。Ⅰ區(qū)西部和Ⅱ區(qū)東部上升趨勢微弱,平均為0.22℃/10a。TNn除瑞安、定海和如皋3個站點(diǎn)外,均表現(xiàn)為上升趨勢,平均為0.59℃/10a;TXx有78%的站點(diǎn)呈上升趨勢,主要集中在Ⅱ區(qū)和Ⅲ區(qū)東部,其余地區(qū)上升趨勢比較微弱,Ⅰ區(qū)西北部部分站點(diǎn)還表現(xiàn)為下降趨勢。全區(qū)TNx呈上升趨勢,其中Ⅱ區(qū)東南部部分地區(qū)上升速率較大,上升速率較小的地區(qū)零散分布在各區(qū)??傮w來看,TNx變化趨勢的空間差異表現(xiàn)不明顯。

    92%的站點(diǎn)WSDI處于增加趨勢,主要分布在Ⅲ區(qū)和Ⅱ區(qū)除西南部之外的大部分地區(qū),平均增加趨勢2.61 d/10a。87%的站點(diǎn)CSDI呈減少趨勢。GSL總體分布為增加趨勢,大約占到總站點(diǎn)數(shù)的94%,增加趨勢較明顯的區(qū)域是Ⅰ區(qū)和Ⅱ區(qū)東北部,平均增幅為4.90 d/10a。DTR的空間分布下降趨勢的站點(diǎn)數(shù)多于上升的站點(diǎn)數(shù),占總站點(diǎn)數(shù)的77%,Ⅰ區(qū)北部下降趨勢比較明顯,碭山站下降趨勢達(dá)-0.56℃/10a。

    2.2 極端氣溫指數(shù)的突變分析

    M-K突變檢驗結(jié)果表明(表3),近55a長江三角洲地區(qū)16個極端氣溫指數(shù)均有突變發(fā)生。極端冷指數(shù)(ID、FD、TX10、TN10、CSDI)在突變時間點(diǎn)前后均呈現(xiàn)由增加到減少的轉(zhuǎn)變;極端暖指數(shù)除SU在突變點(diǎn)后微弱減少外, TR、TX90、TN90、WSDI、GSL在突變點(diǎn)前后均表現(xiàn)出由減少到增加的轉(zhuǎn)變;極值指數(shù)中,除TXx在突變點(diǎn)后呈下降趨勢外,其它極值指數(shù)在突變點(diǎn)前后都經(jīng)歷了由低到高的轉(zhuǎn)變;DTR在突變點(diǎn)前后呈由上升到下降的轉(zhuǎn)變。多數(shù)極端氣溫指數(shù)在突變點(diǎn)前后的變化特征表明長江三角洲地區(qū)對全球變暖有較好的響應(yīng)。對比長江三角洲地區(qū)與長江流域、珠江流域及全國極端氣溫的突變情況[21-22],發(fā)現(xiàn)多數(shù)突變均發(fā)生在20世紀(jì)80年代至2000年初,并且只發(fā)生過一次突變。

    表3 長江三角洲地區(qū)極端氣溫指數(shù)突變分析

    2.3 極端氣溫指數(shù)周期變化分析

    圖3為長江三角洲地區(qū)ID和FD小波方差圖。ID在1968—1978年以4a尺度為中心存在明顯的高低震蕩。1960—1983年在5—10a的時間尺度上存在明顯的高低震蕩,1960—1990年以17a為中心存在明顯高低震蕩。由小波方差圖可知,ID存在4a、6a、11a、17a的震蕩周期,其中17a為主周期。FD在1960—2014年有4—6a的高低震蕩變化,1960—1990年以21a為中心存在兩次明顯的高低震蕩變化。從小波方差圖得知,FD存在2a、5a、11a、21a的震蕩周期,其中21a為主周期。其它極端氣溫指數(shù)小波分析特征見表4。

    圖3 長江三角洲地區(qū)ID和FD小波方差圖Fig.3 The wavelet variances of ID and FD in the Yangtze River Delta

    2.4 極端氣溫指數(shù)與海拔、經(jīng)緯度的關(guān)系

    表5是16項極端氣溫指數(shù)與地理位置的相關(guān)系數(shù)。多數(shù)極端氣溫指數(shù)與緯度、經(jīng)度和海拔顯著相關(guān)。ID、FD 、SU、TX90、DTR、TN10、TXx、WSDI與緯度顯著負(fù)相關(guān),與經(jīng)度顯著正相關(guān)。ID 與緯度的相關(guān)系數(shù)為-0.86,TX90與經(jīng)度的相關(guān)系數(shù)為0.61。TX10、TNn與緯度顯著正相關(guān),與經(jīng)度顯著正相關(guān)。TR、TN90、TNx、GSL與緯度顯著正相關(guān),TNx與經(jīng)度顯著正相關(guān)。TR、TN90、TXx、TNx、WSDI、GSL與海拔負(fù)相關(guān),ID、FD、TN10與海拔正相關(guān)。

    表4 長江三角洲地區(qū)極端氣溫指數(shù)周期分析

    表5長江三角洲地區(qū)極端氣溫指數(shù)與地理位置的相關(guān)系數(shù)

    Table5ComparisonofcorrelationcoefficientbetweenthelineartrendsinextremetemperatureindicesandgeographicalparametersintheYangtzeRiverDelta

    IDFDSUTRTX10TN10TX90TN90緯度Latitude-0.86***-0.72***-0.66***0.29*0.33**-0.33**-0.59***0.29*經(jīng)度Longitude0.55***0.28*0.50***0.06-0.71***0.050.61***0.12海拔Altitude0.30*0.38**0.01-0.30*0.150.25*-0.15-0.34**TXnTNnTXxTNxWSDICSDIGSLDTR緯度Latitude-0.180.54***-0.40**0.04**-0.27*0.190.65***-0.67***經(jīng)度Longitude-0.19-0.53***0.56***0.13**0.42***-0.090.090.40***海拔Altitude0.15-0.17-0.32**-0.28**-0.41***-0.19-0.43***0.22

    ***、**、*表示通過0.001、0.01和0.05水平的顯著性檢驗

    從不同緯度帶極端氣溫指數(shù)的變化趨勢來看(表6),大多數(shù)極端氣溫指數(shù)在高緯度帶變幅大,在低緯度帶變幅小。ID、FD、TN10、TX90隨緯度的增加而減少,即緯度每增加1°,ID、FD、TN10、TX90的變化幅度分別減少0.15d/10a、0.68d/10a、0.47 d/10a、-0.99 d/10a。SU和TR的增加幅度在29°N分別是最大值和最小值。TX10、TN90的變幅在29°N最大,在29°N—35°N之間隨著緯度的增加而增加。TNx和GSL的變幅隨著緯度的增加而增加,TNx在28°N—32°N變幅呈增加趨勢,GSL在31°N—35°N變幅呈增加趨勢。WSDI和DTR的變幅隨著緯度的增加而減少,其變幅分別在31°N—35°N、30°N—35°N隨緯度增加呈減少趨勢,TXn變幅在28°N—32°N之間變化平穩(wěn),30°N以南變幅增加。TXx和CSDI隨緯度的變化并未表現(xiàn)出明顯的增減趨勢。

    2.5 長江三角洲地區(qū)極端氣溫指數(shù)與區(qū)域氣候變化的相關(guān)分析

    為了解長江三角洲地區(qū)極端氣溫指數(shù)變化與區(qū)域氣溫增暖的關(guān)系,計算了1960—2014年年平均氣溫與各極端氣溫指數(shù)的相關(guān)系數(shù)(表7),可以看出12項極端氣溫指數(shù)與平均氣溫的相關(guān)性較好,其中ID、FD、SU、GSL、DTR與平均氣溫的相關(guān)系數(shù)在0.6以上,表明長江三角洲地區(qū)極端氣溫指數(shù)變化對全球氣候變暖有較好的響應(yīng)。

    表6 長江三角洲地區(qū)不同緯度帶極端氣溫指數(shù)的平均變化趨勢

    表7 長江三角洲地區(qū)極端氣溫指數(shù)變化趨勢與平均氣溫的相關(guān)系數(shù)

    ***、**、*分別表示通過0.001、0.01和0.05的顯著性檢驗

    為進(jìn)一步探討區(qū)域氣候變暖對極端氣溫指數(shù)的影響,通過M-K方法對長江三角洲地區(qū)年平均氣溫做突變檢驗發(fā)現(xiàn)1996年是明顯突變點(diǎn),依此為節(jié)點(diǎn)分別統(tǒng)計突變點(diǎn)后平均氣溫和各極端氣溫指數(shù)的變化趨勢(表8)。可看出,在1996年突變后,平均氣溫、各項暖指數(shù)SU、TR、TX90、TN90、WSDI、GSL和極值指數(shù)均表現(xiàn)為上升趨勢,冷指數(shù)ID、FD、TX10、TN10、CSDI和DTR均呈下降趨勢??梢娙蜃兣尘跋?長江三角洲地區(qū)變暖趨勢明顯。

    表8 長江三角洲地區(qū)平均氣溫突變后平均氣溫和極端氣溫指數(shù)變化

    2.6 極端氣溫指數(shù)與ENSO的相關(guān)分析

    圖4 1960—2014年厄爾尼諾/拉尼娜事件的發(fā)生強(qiáng)度 Fig.4 The occurrence intensity of El Nio /La Nina events from 1960 to 2014

    統(tǒng)計1960—2014年的ENSO事件,共發(fā)生35次ENSO事件,其中13個年份發(fā)生暖事件,22個年份發(fā)生冷事件。ENSO事件的發(fā)生具有周期性特點(diǎn)(圖4),El Nio 事件有3個高峰期,分別為1982—1983年、1991—1992年、1997—1998年。La Nina事件有3個高峰期,分別是1973—1976年、1988—1989年、1999—2000年。

    圖5為長江三角洲地區(qū)ENSO事件與極端氣溫指數(shù)ID和WSDI距平關(guān)系圖。由圖可知,多數(shù)厄爾尼諾年ID在減少,負(fù)距平幾率是85%,69%的厄爾尼諾年下一年ID距平小于0,54%的厄爾尼諾年下一年ID距平大于厄爾尼諾年。68%的拉尼娜年ID距平小于0,32%的拉尼娜年下一年ID距平大于0,50%的拉尼娜年下一年ID距平大于拉尼娜年。半數(shù)厄爾尼諾年WSDI在增加,正距平的幾率是54%,54%的厄爾尼諾年下一年WSDI距平大于0,46%的厄爾尼諾年下一年WSDI距平大于厄爾尼諾年。55%的拉尼娜年WSDI距平大于0,59%的拉尼娜年下一年WSDI距平小于0,59%的拉尼娜年下一年WSDI距平小于拉尼娜年。

    圖5 ENSO事件與長江三角洲地區(qū)冰凍日數(shù)、異常暖晝持續(xù)指數(shù)距平關(guān)系圖Fig.5 The relationship between ENSO events and ice days, warm spell duration days anomaly

    表9是厄爾尼諾年、拉尼娜年及各自滯后一年16個極端氣溫指數(shù)的統(tǒng)計值。分析可知,暖指日數(shù)SU、TR、TN90、WSDI、DSL和極值指數(shù)TXn、TXx、TNx在厄爾尼諾年或厄爾尼諾滯后一年大于正常年。在厄爾尼諾年或厄爾尼諾滯后一年,多數(shù)暖指日數(shù)和極值指數(shù)的值在5類年份中最大,表明厄爾尼諾現(xiàn)象加劇了長江三角洲地區(qū)暖指日數(shù)增加和極值指數(shù)上升。在拉尼娜滯后一年,冷指數(shù)ID、FD、TX10、TN10多年平均值大于厄爾尼諾年(厄爾尼諾滯后一年)和正常年,表明拉尼娜現(xiàn)象使得冷指日數(shù)增加。冷持續(xù)日數(shù)(CSDI)和月平均日較差(DTR)的值在厄爾尼諾年大于正常年和拉尼娜年,表明厄爾尼諾現(xiàn)象對這兩個指數(shù)的影響較大。總體上,厄爾尼諾對極端暖指數(shù)和極值指數(shù)的變化影響較大,對極端冷指數(shù)的影響較小。厄爾尼諾對極端氣溫指數(shù)的影響發(fā)生在當(dāng)年或次年,拉尼娜對極端氣溫指數(shù)的影響發(fā)生在拉尼娜滯后一年,厄爾尼諾對極端氣溫的影響大于拉尼娜對極端氣溫的影響。

    表9長江三角洲地區(qū)ENSO事件發(fā)生前后極端氣溫指數(shù)統(tǒng)計

    Table9ComparisonsofaverageextremetemperatureindicesinENSOyears(and lag 1 years)andnormalyearsoverYangtzeRiverDelta

    指數(shù)類型Indicestype指數(shù)Indices暖事件年Warmeventsyears暖事件年+1Warmeventsoneyearlag正常年Normalyears冷事件年Coldeventsyears冷事件年+1Coldeventsoneyearlag絕對指數(shù)ID(d)1.281.632.061.662.08AbsoluteindicesFD(d)40.5043.2443.1246.0046.48SU(d)142.52136.07141.67135.77138.96TR(d)96.1897.4998.0098.1495.95相對指數(shù)TX10(d)33.8335.0932.3534.6734.83RelativeindicesTN10(d)31.5227.9529.1329.8433.49TX90(d)44.9345.4646.4742.5539.87TN90(d)53.7350.6751.0646.2944.48極值指數(shù)TXn(℃)0.190.11-0.50-0.02-0.42RelativeindicesTNn(℃)-7.30-7.45-7.47-7.23-7.65TXx(℃)36.9137.3837.3537.0136.78TNx(℃)27.7228.7228.1327.6827.64其它指數(shù)WSDI(d)6.438.056.265.473.39OtherindicesCSDI(d)1.921.561.261.311.48GSL(d)324.22319.76321.23319.36319.72DTR(℃)8.418.308.368.318.35

    3 結(jié)論與討論

    3.1 結(jié)論

    (1)近55年長江三角洲地區(qū)極端暖指日數(shù)(SU、TR、TX90、TN90、WSDI、GSL)和極值指數(shù)(TXn、TNn、TXx、TNx)呈顯著上升趨勢,冷指日數(shù)(ID、FD、TX10、TN10、CSDI)和月平均日較差(DTR)呈顯著下降趨勢??傮w上,暖指數(shù)的上升趨勢大于冷指數(shù)的下降趨勢,全區(qū)范圍所有指數(shù)的變化趨勢均通過了0.05 的顯著性檢驗。

    (2)空間分布上,3個氣候區(qū)極端氣溫指數(shù)變化趨勢有所差異,多數(shù)極端氣溫指數(shù)在長江三角洲地區(qū)北部和南部變化較明顯。92%以上的站點(diǎn)暖指日數(shù)SU、TR、TX90、TN90、WSDI、GSL呈增加趨勢,85%的站點(diǎn)TX90呈增加趨勢;所有站點(diǎn)冷指日數(shù)ID、FD、TX10、TN10呈減少趨勢,87%的站點(diǎn)CSDI呈減少趨勢;除TXx外,95%以上的站點(diǎn)其它極值指數(shù)為上升趨勢。

    (3)多數(shù)極端氣溫指數(shù)突變發(fā)生在20世紀(jì)80年代至21世紀(jì)初。冷指日數(shù)突變前后經(jīng)歷了由增加到減少的轉(zhuǎn)變,多數(shù)暖指日數(shù)在突變點(diǎn)前后呈現(xiàn)由減少到增加的轉(zhuǎn)變。周期分析表明,多數(shù)極端氣溫指數(shù)呈現(xiàn)3—5個周期震蕩,主周期均在10a以上。

    (4)多數(shù)極端氣溫指數(shù)與海拔、緯度、經(jīng)度顯著相關(guān),不同指數(shù)與地理位置的相關(guān)系數(shù)大小有差異。極端氣溫指數(shù)變化與區(qū)域氣候變暖密切相關(guān),平均氣溫突變后,ID、FD、TX10、TN10等冷指日數(shù)明顯減少,SU、TR、TX90、TN90等暖指日數(shù)明顯增加,表明極端氣溫指數(shù)對氣候變暖有較好的響應(yīng)。

    (5)厄爾尼諾對極端氣溫的影響發(fā)生在當(dāng)年或次年,拉尼娜對極端氣溫的影響發(fā)生在拉尼娜滯后1年,總體上,厄爾尼諾對極端氣溫的影響大于拉尼娜的影響。

    3.2 討論

    采用16個極端氣溫指數(shù),分析了長江三角洲地區(qū)1960—2014年極端氣溫指數(shù)的年際變化趨勢和空間分布差異,結(jié)果表明該區(qū)最低氣溫與最高氣溫相關(guān)的極端氣溫指數(shù)變化趨勢與總體變暖一致,這與全球升溫背景下很多區(qū)域的研究結(jié)論基本一致[20- 25]。聶安祺[26]研究中國三大城市帶城市化氣候效應(yīng),得出長江三角洲地區(qū)的氣候效應(yīng)最強(qiáng)。本文得出的結(jié)論也印證了這一點(diǎn),表明快速城市化背景下長江三角洲地區(qū)氣溫上升明顯。有研究[27-29]發(fā)現(xiàn),地面氣溫的上升是不對稱的,即最低氣溫上升明顯,而最高氣溫上升趨勢較弱,氣溫日較差明顯變小。本研究得出長江三角洲地區(qū)多數(shù)基于最低氣溫的極端氣溫指數(shù)變化幅度大于基于最高氣溫的極端氣溫指數(shù)變化幅度,氣溫日較差明顯下降,進(jìn)一步驗證了地面氣溫上升中的“非對稱性”現(xiàn)象。周雅清和張雷等人[20,29]的研究表明,國內(nèi)外大量研究證實的“氣溫非對稱性變化”現(xiàn)象,主要與城市化對地面最低氣溫的明顯增加效應(yīng)有關(guān)。由于臺站變遷、觀測儀器變更及城市化的影響,導(dǎo)致不同規(guī)模城市站地面氣溫記錄序列中存在偏差。有研究[19-20,26,29]表明,城市熱島效應(yīng)對區(qū)域平均氣溫有明顯影響,今后計算區(qū)域極端氣溫時,需要對臺站氣溫資料進(jìn)行訂正以獲取消除了城市化效應(yīng)的氣溫數(shù)據(jù),這樣才能更好地分析區(qū)域極端氣溫的變化趨勢。

    長江三角洲地處氣候系統(tǒng)復(fù)雜的亞熱帶季風(fēng)氣候區(qū),ENSO通過影響東亞季風(fēng)的強(qiáng)度和太平洋副熱帶高壓的位置,進(jìn)一步影響季風(fēng)的強(qiáng)度和長江三角洲地區(qū)的氣候。在強(qiáng)El Nio年,長江三角洲極端氣溫暖事件增加,意味著熱浪的發(fā)生頻率增加。在強(qiáng)La Nina年,極端氣溫冷事件增加,寒潮的發(fā)生頻率增加。李崇銀等[30]研究表明,ENSO主要是由東亞季風(fēng)異常造成的赤道西太平洋異常緯向風(fēng)所驅(qū)動的熱帶太平洋次表層海溫距平的循環(huán)。在厄爾尼諾年冬季,長江三角洲為異常偏南風(fēng),東亞冬季風(fēng)偏弱,導(dǎo)致降水偏多氣溫冷濕;拉尼娜年冬季,長江三角洲為異常偏北風(fēng),冬季風(fēng)偏強(qiáng),導(dǎo)致降水偏少氣溫暖干。此外,長江三角洲地區(qū)的極端氣溫事件與ENSO的持續(xù)時間有關(guān),ENSO持續(xù)時間的長短又與大氣環(huán)流的持續(xù)時間密切相關(guān)。ENSO對長江三角洲地區(qū)氣溫影響機(jī)理復(fù)雜,未來要全面了解ENSO對極端氣溫的影響,有必要研究ENSO與低緯和中、高緯度環(huán)流的關(guān)系及其變化的物理過程。雖然ENSO不同位相長江三角洲地區(qū)極端氣溫差異明顯,但顯然ENSO不是影響區(qū)域極端氣溫變化的唯一原因,除了全球變暖與地理位置因素外,氣候系統(tǒng)內(nèi)部變率、城市溫室氣體和大氣成分變化在內(nèi)的外部強(qiáng)迫都與極端氣溫事件變化密切相關(guān),這些要素對極端氣溫的影響有待于做進(jìn)一步研究。

    [1] 秦大河, Stocker T, 259名作者和TSU (駐伯爾尼和北京). IPCC第五次評估報告第一工作組報告的亮點(diǎn)結(jié)論. 氣候變化研究進(jìn)展, 2014, 10(1): 1- 6.

    [2] 鄭祚芳. 北京極端氣溫變化特征及其對城市化的響應(yīng). 地理科學(xué), 2011, 31(4): 459- 463.

    [3] 秦大河. 中國極端天氣氣候事件和災(zāi)害風(fēng)險管理與適應(yīng)國家評估報告. 北京: 科學(xué)出版社, 2015:5-11.

    [4] 崔林麗, 史軍, 楊引明, 李貴才, 范文義. 長江三角洲氣溫變化特征及城市化影響. 地理研究, 2008, 27(4): 775- 786.

    [5] 吳榮軍, 鄭有飛, 劉建軍, 談建國, 許遐禎, 余永江. 長江三角洲主要城市高溫災(zāi)害的趨勢分析. 自然災(zāi)害學(xué)報, 2010, 19(5): 56- 63.

    [6] 閆俊, 江俊杰. 長江三角洲近51年極端氣候事件變化及其原因分析. 廣東氣象, 2012, 34(5): 29- 34.

    [7] Timilsena J, Piechota T, Tootle G, Singh A. Associations of interdecadal/interannual climate variability and long-term Colorado River Basin streamflow. Journal of Hydrology, 2009, 365(3/4): 289- 301.

    [8] 孫建奇, 袁薇, 高玉中. 阿拉伯半島-北太平洋型遙相關(guān)及其與亞洲夏季風(fēng)的關(guān)系. 中國科學(xué) D輯: 地球科學(xué), 2008, 38(6): 750- 762.

    [9] 陳文, 蘭曉青, 王林, 馬音. ENSO和北極濤動對東亞冬季氣候異常的綜合影響. 科學(xué)通報, 2013, 58(8): 634- 641.

    [10] 韓文韜, 衛(wèi)捷, 沈新勇. 近50年中國冬季氣溫對ENSO響應(yīng)的時空穩(wěn)定性分析研究. 氣候與環(huán)境研究, 2014, 19(1): 97- 106.

    [11] 伍紅雨, 潘蔚娟, 王婷. 華南冬季氣溫異常與ENSO的關(guān)系. 氣象, 2014, 40(10): 1230- 1239.

    [12] 徐同, 郭品文. 中國冬季氣溫對ENSO的非線性響應(yīng)特征. 氣象與減災(zāi)研究, 2007, 30(4): 9- 13.

    [13] 解明恩, 張萬誠. ENSO事件與云南冬季氣溫異常. 應(yīng)用氣象學(xué)報, 2000, 11(1): 115- 122.

    [14] 梁蘇潔, 丁一匯, 趙南, 孫穎. 近50年中國大陸冬季氣溫和區(qū)域環(huán)流的年代際變化研究. 大氣科學(xué), 2014, 38(5): 974- 992.

    [15] 魏鳳英. 現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù). 北京: 氣象出版社, 1999:37-40.

    [16] Hamed K H. Trend detection in hydrologic data: The Mann-Kendall trend test under the scaling hypothesis. Journal of Hydrology, 2008, 349(3/4): 350- 363.

    [17] 劉賢趙, 張安定, 李嘉竹. 地理學(xué)數(shù)學(xué)方法. 北京: 科學(xué)出版社, 2009:190-202.

    [18] 李曉燕, 翟盤茂. ENSO事件指數(shù)與指標(biāo)研究. 氣象學(xué)報, 2000, 58(1): 102- 109.

    [19] 謝志清, 杜銀, 曾燕, 施雅風(fēng), 武金崗. 長江三角洲城市帶擴(kuò)展對區(qū)域溫度變化的影響. 地理學(xué)報, 2007, 62(7): 717- 727.

    [20] 周雅清, 任國玉. 城市化對華北地區(qū)極端氣溫事件頻率的影響. 高原氣象, 2014, 33(6): 1589- 1598.

    [21] 王瓊, 張明軍, 王圣杰, 駱書飛, 汪寶龍, 朱小凡. 1962—2011年長江流域極端氣溫事件分析. 地理學(xué)報, 2013, 68(5): 611- 625.

    [22] 劉青娥, 吳孝情, 陳曉宏, 楊冰. 珠江流域1960—2012年極端氣溫的時空變化特征. 自然資源學(xué)報, 2015, 30(8): 1356- 1366.

    [23] 趙安周, 劉憲鋒, 朱秀芳, 潘耀忠, 趙玉玲, 王冬利. 1965—2013年黃土高原地區(qū)極端氣溫趨勢變化及空間差異. 地理研究, 2016, 35(4): 639- 652.

    [24] 閆慧敏, 陳偉娜, 楊方興, 劉紀(jì)遠(yuǎn), 胡云鋒, 冀詠贊. 過去50年內(nèi)蒙古極端氣候事件時空格局特征. 地理研究, 2014, 33(1): 13- 22.

    [25] 汪寶龍, 張明軍, 魏軍林, 王圣杰, 馬潛, 李小飛. 西北地區(qū)近50a氣溫和降水極端事件的變化特征. 自然資源學(xué)報, 2012, 27(10): 1720- 1733.

    [26] 聶安琪, 陳星, 馮志剛. 中國三大城市帶城市化效應(yīng)的檢測與對比. 氣象科學(xué), 2011, 31(4): 372- 383.

    [27] Qian W H, Lin X. Regional trends in recent temperature and indices in China. Climate Research, 2004, 27(2): 119- 134.

    [28] 華麗娟, 馬柱國, 羅德海. 1961—2000年中國區(qū)域氣溫較差分析. 地理學(xué)報, 2004, 59(5): 680- 688.

    [29] 張雷, 任國玉, 劉江, 周雅清, 任玉玉, 張愛英, 馮亞文. 城市化對北京氣象站極端氣溫指數(shù)趨勢變化的影響. 地球物理學(xué)報, 2011, 54(5): 1150- 1159.

    [30] 李崇銀, 穆穆, 周廣慶, 楊輝. ENSO機(jī)理及其預(yù)測研究. 大氣科學(xué), 2008, 32(4): 761- 781.

    CorrelationanalysisofvariationofextremetemperatureeventsandElNioSouthernOscillation(ENSO)inYangtzeRiverDeltaregionduring1960—2014

    JIA Yanqing1,2, ZHANG Bo1,*, ZHANG Yaozong1, TANG Min1, MA Bin1, WANG Guoqiang1

    1CollegeofGeographyandEnvironmentalScience,NorthwestNormalUniversity,Lanzhou730070,China2DepartmentofGeography,XinzhouNormalUniversity,Xinzhou034000,China

    Based on daily temperature (maximum, minimum, and average) data of 65 stations, spatial and temporal changes of extreme temperature events and the El Nio Southern Oscillation (ENSO) influence on the extreme temperature threshold over Yangtze River Delta were analyzed for the period 1960—2014. The results showed that in the last 55 years, the extreme temperature indices, i.e. summer days (SU), tropical nights (TR), warm days (TX90), warm nights (TN90), warm spell duration days (WSDI), and growing season length (GSL) have all increased significantly, and TN90 showed an increasing rate of 8.55d per decade. The extremal indices (TXn, TNn, TXx, and TNx) had increasing trends in the Yangtze River Delta region, and TNn showed an increasing rate of 0.53℃ per decade. The ice days (ID), frost days (FD), cold days (TX10), cold nights (TN10), and cold spell duration days (CSDI) have decreased significantly by 0.40 d, 4.03 d, 2.56 d, 6.06 d and 0.33 d per decade, respectively, and TN10 showed a decreasing rate of 6.06 d per decade. DTR has decreased by 0.11℃ per decade. For the spatial distribution of SU, TXn and TNx, all meteorological stations showed increasing trends during 1960—2014. TR, TX90, TN90, TNn, TXx, WSDI, and GSL showed increasing trends at 97%, 85%, 98%, 95%, 78%, 92%, and 94% of meteorological stations, respectively. ID, FD, TX10, and TN10 showed decreasing trends for all stations. For CSDI and DTR, 87% and 77% stations showed decreasing trend. Most extreme temperature indices show good correlation with latitude, longitude, and altitude. After the detected abrupt climate changes, extreme warm indices increased, and extreme cold indices decreased significantly. Overall, the effect of El Nio on extreme temperature indices was greater than the effect of La Nina.

    extreme temperature indices; spatial distribution; climate change; ENSO

    國家自然科學(xué)基金項目(41561024);高校博士學(xué)科點(diǎn)專項科研基金項目(20136203110002);山西省高等學(xué)校哲學(xué)社會科學(xué)研究項目(2015265)

    2016- 07- 09; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 05- 27

    *通訊作者Corresponding author.E-mail: zhangbo@nwnu.edu.cn

    10.5846/stxb201607091409

    賈艷青, 張勃,張耀宗,唐敏,馬彬,王國強(qiáng).長江三角洲地區(qū)極端氣溫事件變化特征及其與ENSO的關(guān)系.生態(tài)學(xué)報,2017,37(19):6402- 6414.

    Jia Y Q, Zhang B, Zhang Y Z, Tang M, Ma B, Wang G Q.Correlation analysis of variation of extreme temperature events and El Nio Southern Oscillation(ENSO) in Yangtze River Delta region during 1960—2014.Acta Ecologica Sinica,2017,37(19):6402- 6414.

    猜你喜歡
    長江三角洲地區(qū)拉尼厄爾尼諾
    長江三角洲地區(qū)大氣污染物時空分布特征研究
    拉尼娜來了,這個“小女孩”不簡單
    “一帶一路”倡議下長江三角洲地區(qū)對非貿(mào)易的進(jìn)展與展望
    厄爾尼諾現(xiàn)象橫行中國
    梅拉尼婭:世界上最受欺凌的人是我
    最小的一個
    厄爾尼諾現(xiàn)象及其影響
    長江三角洲地區(qū)工業(yè)設(shè)計咨詢企業(yè)設(shè)計流程與方法初探
    三步圖解厄爾尼諾與拉尼娜
    基于引力模型的區(qū)域農(nóng)產(chǎn)品物流網(wǎng)絡(luò)構(gòu)建研究——以長江三角洲地區(qū)為例
    404 Not Found

    404 Not Found


    nginx
    国内揄拍国产精品人妻在线| 特大巨黑吊av在线直播| 亚洲精品国产精品久久久不卡| 成人特级黄色片久久久久久久| 脱女人内裤的视频| 日本黄色片子视频| 桃红色精品国产亚洲av| 国模一区二区三区四区视频 | 中出人妻视频一区二区| 禁无遮挡网站| 人妻丰满熟妇av一区二区三区| 18美女黄网站色大片免费观看| 精品久久久久久久毛片微露脸| 久久天堂一区二区三区四区| 国产又色又爽无遮挡免费看| 国产亚洲精品久久久久久毛片| 啦啦啦韩国在线观看视频| 亚洲天堂国产精品一区在线| 久久久国产成人精品二区| 一边摸一边抽搐一进一小说| 宅男免费午夜| 一级毛片高清免费大全| 欧美日韩一级在线毛片| 成人无遮挡网站| 欧美zozozo另类| 久久中文字幕人妻熟女| 天堂动漫精品| 男女做爰动态图高潮gif福利片| 亚洲成人中文字幕在线播放| 亚洲熟女毛片儿| 欧美乱色亚洲激情| 人妻丰满熟妇av一区二区三区| 亚洲欧美精品综合久久99| 欧美日韩福利视频一区二区| 少妇人妻一区二区三区视频| 757午夜福利合集在线观看| 国产av麻豆久久久久久久| 看黄色毛片网站| 啦啦啦韩国在线观看视频| 免费观看精品视频网站| 日韩欧美在线乱码| 亚洲无线观看免费| 国内精品久久久久精免费| 亚洲第一电影网av| 久9热在线精品视频| 久久精品91无色码中文字幕| 欧美日韩乱码在线| 午夜亚洲福利在线播放| 中文字幕熟女人妻在线| 在线观看美女被高潮喷水网站 | 人人妻,人人澡人人爽秒播| 国产成+人综合+亚洲专区| 精品国产乱子伦一区二区三区| 人妻久久中文字幕网| 日本黄色片子视频| 在线观看免费午夜福利视频| 亚洲国产色片| 午夜精品在线福利| 久9热在线精品视频| 欧美日韩福利视频一区二区| 日韩高清综合在线| 国产精品 国内视频| 亚洲无线在线观看| 午夜免费观看网址| 精华霜和精华液先用哪个| 免费在线观看影片大全网站| 激情在线观看视频在线高清| 国产欧美日韩精品一区二区| 法律面前人人平等表现在哪些方面| 美女 人体艺术 gogo| 午夜福利在线观看吧| 一夜夜www| 嫩草影院精品99| 99精品欧美一区二区三区四区| 欧美在线黄色| 国产成人系列免费观看| 日本免费一区二区三区高清不卡| 一个人看的www免费观看视频| 亚洲专区字幕在线| 最近最新中文字幕大全电影3| 久久久久国产一级毛片高清牌| 国产美女午夜福利| 一进一出抽搐动态| 亚洲五月天丁香| 日本a在线网址| 亚洲欧美日韩卡通动漫| 91久久精品国产一区二区成人 | 在线a可以看的网站| 欧美日本视频| 免费在线观看日本一区| 日日夜夜操网爽| 日日摸夜夜添夜夜添小说| 精品99又大又爽又粗少妇毛片 | 黄片大片在线免费观看| 久久久久久久精品吃奶| 亚洲av第一区精品v没综合| 亚洲狠狠婷婷综合久久图片| 亚洲国产欧洲综合997久久,| 国产精品一及| 一进一出好大好爽视频| 草草在线视频免费看| 婷婷六月久久综合丁香| 国产三级在线视频| 亚洲乱码一区二区免费版| 国产成人系列免费观看| 国产成人系列免费观看| 欧美+亚洲+日韩+国产| 2021天堂中文幕一二区在线观| 午夜福利18| 大型黄色视频在线免费观看| 亚洲色图av天堂| 国产欧美日韩精品一区二区| 伦理电影免费视频| 99久久综合精品五月天人人| 99久国产av精品| 精品久久久久久成人av| 级片在线观看| 国内精品久久久久精免费| 高潮久久久久久久久久久不卡| 国产精品电影一区二区三区| 天天一区二区日本电影三级| 色老头精品视频在线观看| 久久精品亚洲精品国产色婷小说| 国产免费男女视频| 窝窝影院91人妻| 日本一本二区三区精品| 制服人妻中文乱码| 日韩欧美国产在线观看| svipshipincom国产片| 亚洲欧美精品综合久久99| 亚洲欧美日韩高清在线视频| 日韩国内少妇激情av| 日本成人三级电影网站| 一个人看的www免费观看视频| 亚洲色图av天堂| h日本视频在线播放| 日本成人三级电影网站| 全区人妻精品视频| 亚洲中文字幕日韩| 午夜两性在线视频| 日本成人三级电影网站| 全区人妻精品视频| 国内精品美女久久久久久| 国产午夜福利久久久久久| 精品国产亚洲在线| 两个人的视频大全免费| 国产精品综合久久久久久久免费| 国产爱豆传媒在线观看| 伊人久久大香线蕉亚洲五| 欧美+亚洲+日韩+国产| 神马国产精品三级电影在线观看| 啦啦啦免费观看视频1| 亚洲精品色激情综合| 天天添夜夜摸| 国产成人精品无人区| 国产亚洲精品久久久com| 好男人在线观看高清免费视频| 国产成人av教育| 无人区码免费观看不卡| 欧美丝袜亚洲另类 | 亚洲中文av在线| 91在线精品国自产拍蜜月 | 亚洲精品粉嫩美女一区| 欧美性猛交黑人性爽| 综合色av麻豆| 欧美成人性av电影在线观看| 又大又爽又粗| 亚洲国产色片| e午夜精品久久久久久久| 丰满人妻熟妇乱又伦精品不卡| 国产高清videossex| 日本三级黄在线观看| 日韩大尺度精品在线看网址| 欧美最黄视频在线播放免费| 欧美日韩瑟瑟在线播放| 不卡av一区二区三区| 三级男女做爰猛烈吃奶摸视频| 男人舔奶头视频| 久久久国产欧美日韩av| 看片在线看免费视频| 欧美一级毛片孕妇| 亚洲在线观看片| 欧美日本视频| 欧美日韩综合久久久久久 | 丰满人妻熟妇乱又伦精品不卡| 国产精品99久久久久久久久| 精品福利观看| 精品电影一区二区在线| 亚洲美女黄片视频| 国产免费av片在线观看野外av| 精品乱码久久久久久99久播| 九色成人免费人妻av| 午夜福利免费观看在线| 99国产综合亚洲精品| 国产伦人伦偷精品视频| 最近最新免费中文字幕在线| 别揉我奶头~嗯~啊~动态视频| 老司机午夜十八禁免费视频| 首页视频小说图片口味搜索| 精品99又大又爽又粗少妇毛片 | 成人三级做爰电影| 午夜久久久久精精品| 丁香六月欧美| 欧美绝顶高潮抽搐喷水| 免费av不卡在线播放| 1000部很黄的大片| 久久人人精品亚洲av| av欧美777| 长腿黑丝高跟| 在线观看美女被高潮喷水网站 | 欧美一级毛片孕妇| 日韩大尺度精品在线看网址| 757午夜福利合集在线观看| 99国产极品粉嫩在线观看| 他把我摸到了高潮在线观看| 日日摸夜夜添夜夜添小说| 怎么达到女性高潮| 国产成人精品久久二区二区免费| 亚洲色图av天堂| 亚洲精品一区av在线观看| 久久精品aⅴ一区二区三区四区| 90打野战视频偷拍视频| 国产精品久久久久久人妻精品电影| 精品无人区乱码1区二区| 免费在线观看日本一区| 日韩精品中文字幕看吧| 三级国产精品欧美在线观看 | a级毛片在线看网站| 日本在线视频免费播放| 欧美日韩国产亚洲二区| 国产精品98久久久久久宅男小说| 伊人久久大香线蕉亚洲五| 婷婷六月久久综合丁香| 亚洲精品美女久久av网站| 国产伦人伦偷精品视频| 亚洲五月婷婷丁香| 性欧美人与动物交配| xxx96com| 欧美日韩国产亚洲二区| 国产亚洲精品久久久久久毛片| 好男人电影高清在线观看| 一a级毛片在线观看| 亚洲精品美女久久av网站| 这个男人来自地球电影免费观看| 免费在线观看亚洲国产| 黑人欧美特级aaaaaa片| 国产精品一区二区精品视频观看| 三级国产精品欧美在线观看 | 色综合欧美亚洲国产小说| 色综合婷婷激情| 亚洲 欧美 日韩 在线 免费| 国产精品久久久av美女十八| 亚洲激情在线av| 亚洲精品粉嫩美女一区| 91av网站免费观看| 1024手机看黄色片| 级片在线观看| 日本熟妇午夜| 999精品在线视频| 久久精品国产清高在天天线| 亚洲无线在线观看| 久久午夜亚洲精品久久| 法律面前人人平等表现在哪些方面| 在线国产一区二区在线| 男人舔女人的私密视频| 一二三四社区在线视频社区8| 日韩av在线大香蕉| 最新在线观看一区二区三区| 午夜精品一区二区三区免费看| 成人国产综合亚洲| 老熟妇仑乱视频hdxx| 国产主播在线观看一区二区| 露出奶头的视频| 一本久久中文字幕| 99久久综合精品五月天人人| 波多野结衣巨乳人妻| 黄色成人免费大全| 精品久久久久久,| 亚洲18禁久久av| 十八禁网站免费在线| 亚洲熟妇熟女久久| 国产精品久久视频播放| 国产69精品久久久久777片 | 女人高潮潮喷娇喘18禁视频| 男女床上黄色一级片免费看| 高清毛片免费观看视频网站| 九色国产91popny在线| 亚洲国产高清在线一区二区三| 五月伊人婷婷丁香| 嫩草影院精品99| 麻豆国产av国片精品| 高清在线国产一区| 欧美日韩瑟瑟在线播放| 人人妻,人人澡人人爽秒播| 男女之事视频高清在线观看| 99久久国产精品久久久| 99热这里只有是精品50| 黄色成人免费大全| 国产一区二区在线观看日韩 | 国产精品一区二区三区四区免费观看 | 国产伦在线观看视频一区| 午夜两性在线视频| 久久99热这里只有精品18| 99国产综合亚洲精品| 久久这里只有精品19| 午夜日韩欧美国产| 亚洲欧美精品综合一区二区三区| 精华霜和精华液先用哪个| 俄罗斯特黄特色一大片| 久久久国产成人免费| 脱女人内裤的视频| av天堂中文字幕网| 久久人人精品亚洲av| 中文字幕人妻丝袜一区二区| 99久久无色码亚洲精品果冻| 国产成人欧美在线观看| 久久久久久久久中文| ponron亚洲| 国产欧美日韩精品亚洲av| 欧美色视频一区免费| 亚洲男人的天堂狠狠| 成人国产综合亚洲| 黄色 视频免费看| 欧美午夜高清在线| 色播亚洲综合网| 色视频www国产| 国产熟女xx| 波多野结衣高清无吗| 国产精品一区二区三区四区久久| 久久久久国内视频| 亚洲男人的天堂狠狠| 不卡av一区二区三区| 岛国视频午夜一区免费看| 狂野欧美白嫩少妇大欣赏| 一个人看视频在线观看www免费 | 精品福利观看| 国产精品久久久av美女十八| 叶爱在线成人免费视频播放| 亚洲人成电影免费在线| 国产探花在线观看一区二区| 一区二区三区激情视频| 国产精品,欧美在线| 国产精品久久视频播放| 午夜成年电影在线免费观看| 日韩国内少妇激情av| 久久国产乱子伦精品免费另类| 久久久久久大精品| 国产成人av教育| 中文字幕av在线有码专区| 日本免费一区二区三区高清不卡| 19禁男女啪啪无遮挡网站| 久久久色成人| 欧美最黄视频在线播放免费| 三级国产精品欧美在线观看 | 成熟少妇高潮喷水视频| 一个人看视频在线观看www免费 | 国内少妇人妻偷人精品xxx网站 | 日韩大尺度精品在线看网址| 午夜福利在线在线| 制服丝袜大香蕉在线| 色视频www国产| 老司机深夜福利视频在线观看| 黄色 视频免费看| 麻豆国产97在线/欧美| 老司机午夜十八禁免费视频| 日本熟妇午夜| 男女做爰动态图高潮gif福利片| 亚洲五月婷婷丁香| 中文字幕久久专区| 国产精品综合久久久久久久免费| 九色成人免费人妻av| 国产成人一区二区三区免费视频网站| 免费一级毛片在线播放高清视频| 黄色女人牲交| 波多野结衣巨乳人妻| 美女午夜性视频免费| 亚洲精华国产精华精| 亚洲精品乱码久久久v下载方式 | 看免费av毛片| 三级男女做爰猛烈吃奶摸视频| 久久人人精品亚洲av| 国产精品自产拍在线观看55亚洲| 老司机午夜十八禁免费视频| 免费看光身美女| 少妇丰满av| www.www免费av| 一级毛片精品| 久久久国产欧美日韩av| 国产精品自产拍在线观看55亚洲| 中出人妻视频一区二区| 亚洲人成网站在线播放欧美日韩| 国产高清视频在线观看网站| x7x7x7水蜜桃| 老熟妇乱子伦视频在线观看| 午夜免费观看网址| 后天国语完整版免费观看| 中文字幕高清在线视频| 日本五十路高清| 免费观看人在逋| 精品久久久久久成人av| 久久国产精品影院| 亚洲国产欧洲综合997久久,| av女优亚洲男人天堂 | 天天躁狠狠躁夜夜躁狠狠躁| 国产成人精品无人区| www.熟女人妻精品国产| 夜夜看夜夜爽夜夜摸| 亚洲熟妇中文字幕五十中出| 国产日本99.免费观看| 亚洲最大成人中文| 久久婷婷人人爽人人干人人爱| 毛片女人毛片| 小蜜桃在线观看免费完整版高清| 桃色一区二区三区在线观看| 欧美乱色亚洲激情| 波多野结衣高清作品| 国产精品久久电影中文字幕| 一区二区三区激情视频| 国产一区二区激情短视频| 欧美3d第一页| 亚洲成人免费电影在线观看| 一个人免费在线观看的高清视频| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美 国产精品| 18禁美女被吸乳视频| 亚洲在线自拍视频| 日韩成人在线观看一区二区三区| 成年人黄色毛片网站| 亚洲av成人一区二区三| 男女下面进入的视频免费午夜| 国产亚洲av嫩草精品影院| 久久久久久久久中文| 男女下面进入的视频免费午夜| tocl精华| 亚洲精品美女久久av网站| 最近在线观看免费完整版| 老司机深夜福利视频在线观看| 精品人妻1区二区| 久久香蕉国产精品| 欧美成人免费av一区二区三区| 成人永久免费在线观看视频| 日韩欧美精品v在线| 久久性视频一级片| 欧美日本视频| 长腿黑丝高跟| 欧美一级毛片孕妇| 欧美一区二区精品小视频在线| 精品国产乱码久久久久久男人| 国内久久婷婷六月综合欲色啪| 我要搜黄色片| 最近视频中文字幕2019在线8| 亚洲欧美精品综合久久99| 999精品在线视频| 免费大片18禁| 国产精品一区二区免费欧美| 男女午夜视频在线观看| 国产精品久久视频播放| 在线观看舔阴道视频| 日韩欧美一区二区三区在线观看| 午夜激情欧美在线| 又紧又爽又黄一区二区| 男女之事视频高清在线观看| 丁香六月欧美| 丰满的人妻完整版| 无人区码免费观看不卡| 精品福利观看| 久久久久久国产a免费观看| 99热这里只有精品一区 | ponron亚洲| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清| 无遮挡黄片免费观看| 午夜免费激情av| 9191精品国产免费久久| 亚洲片人在线观看| 国产麻豆成人av免费视频| 老熟妇乱子伦视频在线观看| 女生性感内裤真人,穿戴方法视频| 熟女少妇亚洲综合色aaa.| 18禁观看日本| 亚洲黑人精品在线| 久久精品国产综合久久久| 国产私拍福利视频在线观看| 亚洲欧美激情综合另类| 一进一出抽搐gif免费好疼| 亚洲aⅴ乱码一区二区在线播放| 欧美最黄视频在线播放免费| 日韩 欧美 亚洲 中文字幕| tocl精华| h日本视频在线播放| 久久久久久大精品| 国产高清三级在线| 国产不卡一卡二| 黄色片一级片一级黄色片| 亚洲人成网站在线播放欧美日韩| 制服丝袜大香蕉在线| 精品一区二区三区av网在线观看| 午夜影院日韩av| 麻豆国产av国片精品| 亚洲精品一卡2卡三卡4卡5卡| 欧美性猛交黑人性爽| 制服人妻中文乱码| 搡老岳熟女国产| 成人性生交大片免费视频hd| 久久国产精品影院| 久久久精品欧美日韩精品| 午夜精品一区二区三区免费看| 成人三级黄色视频| 99精品久久久久人妻精品| 久久国产精品影院| 欧美激情在线99| 欧美日本视频| 深夜精品福利| 国产成人啪精品午夜网站| 午夜福利18| 精品熟女少妇八av免费久了| 国产精品久久久久久久电影 | 88av欧美| 久久精品国产99精品国产亚洲性色| 成人av一区二区三区在线看| 欧美黑人巨大hd| 亚洲va日本ⅴa欧美va伊人久久| 天天添夜夜摸| 脱女人内裤的视频| 日本精品一区二区三区蜜桃| 欧美另类亚洲清纯唯美| 村上凉子中文字幕在线| 最新在线观看一区二区三区| 黄色 视频免费看| 十八禁人妻一区二区| 成人特级黄色片久久久久久久| 青草久久国产| 不卡一级毛片| 免费观看的影片在线观看| 曰老女人黄片| 免费看美女性在线毛片视频| 精品国产乱子伦一区二区三区| 国产欧美日韩一区二区三| 国产又黄又爽又无遮挡在线| 香蕉av资源在线| 校园春色视频在线观看| 亚洲真实伦在线观看| 国产淫片久久久久久久久 | 少妇人妻一区二区三区视频| 亚洲欧美日韩高清专用| 一区二区三区激情视频| a在线观看视频网站| 免费无遮挡裸体视频| 日本免费a在线| 在线国产一区二区在线| 国产午夜精品久久久久久| 脱女人内裤的视频| 99精品欧美一区二区三区四区| 午夜成年电影在线免费观看| 男人舔女人下体高潮全视频| 欧美中文日本在线观看视频| 日本与韩国留学比较| 久久这里只有精品19| 男人舔女人的私密视频| 在线免费观看不下载黄p国产 | 999精品在线视频| 99riav亚洲国产免费| 国产精品一区二区精品视频观看| 国产成人精品无人区| 婷婷亚洲欧美| 欧美xxxx黑人xx丫x性爽| 全区人妻精品视频| 桃红色精品国产亚洲av| 亚洲中文字幕日韩| 午夜福利视频1000在线观看| 日本黄色视频三级网站网址| 午夜免费激情av| 国产极品精品免费视频能看的| 久久午夜亚洲精品久久| 午夜视频精品福利| 性欧美人与动物交配| 亚洲av成人不卡在线观看播放网| 操出白浆在线播放| 国内精品美女久久久久久| 十八禁人妻一区二区| 变态另类成人亚洲欧美熟女| 最近最新中文字幕大全免费视频| 成年版毛片免费区| 久久天躁狠狠躁夜夜2o2o| 在线观看66精品国产| 这个男人来自地球电影免费观看| 日本成人三级电影网站| a级毛片在线看网站| 黑人巨大精品欧美一区二区mp4| 久久天躁狠狠躁夜夜2o2o| 精品国产亚洲在线| 国产成人aa在线观看| 欧美+亚洲+日韩+国产| 欧美丝袜亚洲另类 | 999久久久精品免费观看国产| 亚洲成av人片免费观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲色图 男人天堂 中文字幕| 亚洲精品在线美女| tocl精华| 国产精品综合久久久久久久免费| 国产成人aa在线观看| 久久欧美精品欧美久久欧美| 日本撒尿小便嘘嘘汇集6| 午夜影院日韩av| 少妇的丰满在线观看| 18禁黄网站禁片免费观看直播| 人妻久久中文字幕网| av欧美777| 国产高清videossex| 99热只有精品国产| 精品国产乱子伦一区二区三区| 长腿黑丝高跟| 一级毛片高清免费大全| 国产精品av久久久久免费|