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

    南海夏季風(fēng)建立日期的預(yù)報(bào)研究*

    2010-06-05 08:03:30藍(lán)光東溫之平程炳巖李永華
    關(guān)鍵詞:季風(fēng)南海風(fēng)速

    周 浩,藍(lán)光東,溫之平,程炳巖,李永華

    (1.重慶市氣候中心,重慶 401147;2.中山大學(xué)季風(fēng)與環(huán)境研究中心∥大氣科學(xué)系,廣東 廣州 510275)

    對東亞季風(fēng)系統(tǒng)的研究[1]表明,南海夏季風(fēng)是亞洲季風(fēng)中爆發(fā)最早的一支季風(fēng)。南海夏季風(fēng)的爆發(fā)往往征兆亞洲季風(fēng)區(qū)大氣環(huán)流由冬季形勢向夏季形勢的轉(zhuǎn)變,東亞夏季風(fēng)的來臨和中國東部雨季的開始。南海夏季風(fēng)爆發(fā)的早晚及其強(qiáng)弱聯(lián)系著大氣環(huán)流的異常并關(guān)系到中國東部地區(qū)的天氣氣候變化。因此,關(guān)于南海夏季風(fēng)爆發(fā)的研究得到廣大專家學(xué)者的高度重視。許多工作從動(dòng)力學(xué)和熱力學(xué)出發(fā)研究南海夏季風(fēng)爆發(fā)的物理機(jī)理,取得了很多重要的成果,但針對季風(fēng)建立日期的預(yù)報(bào)研究則偏少。孫淑清等[2]分析發(fā)現(xiàn)冬季風(fēng)異常與南海夏季風(fēng)存在很好的相關(guān)性,但沒有指出預(yù)報(bào)季風(fēng)爆發(fā)的具體方法。藍(lán)關(guān)東等[3]雖提出了有預(yù)報(bào)作用的前期因子,卻存在所用資料時(shí)間序列較短的問題。周文等[4]建立了預(yù)報(bào)季風(fēng)爆發(fā)的統(tǒng)計(jì)方法,并選取1968-1998年的資料對該方法進(jìn)行了驗(yàn)證,但其預(yù)報(bào)因子缺乏明確的物理含義,研究所用資料也較短。因此,本文利用ECMWF和NCEP 兩套42 a長時(shí)間序列資料,試圖找出與南海夏季風(fēng)爆發(fā)遲早具有較好關(guān)系的預(yù)報(bào)因子,通過使用SAS統(tǒng)計(jì)分析系統(tǒng)建立南海夏季風(fēng)建立日期預(yù)報(bào)方程,從而為準(zhǔn)確的掌握南海夏季風(fēng)建立日期提供科學(xué)的方法。

    1 資料和方法

    1.1 資料

    1)歐洲中期天氣預(yù)報(bào)中心(ECMWF)的全球格點(diǎn)再分析資料,以下簡稱ERA資料。ERA資料分為18層(分別為1 000,925,850,775,700,600,500,400,300,250,200,150,100,70,50,30,20和10 hPa)、水平分辨率為2.5°×2.5°經(jīng)緯度,資料總長度為528個(gè)月(1958年1月-2001年12月)的全球格點(diǎn)再分析資料。ERA分析系統(tǒng)采用ECMWF業(yè)務(wù)資料同化系統(tǒng)的特別版,它包括一個(gè)T106譜預(yù)報(bào)模式,垂直為31層,采用三維半拉格朗日平流格式[5]。

    2)美國環(huán)境預(yù)報(bào)中心(NCEP)/美國國家大氣研究中心(NCAR)第1版全球格點(diǎn)再分析資料(以下簡稱NCEP 資料)。NCEP資料垂直方向分為17層(分別為1 000,925,850,700,600,500,400,300,250,200,150,100,70,50,30,20和10 hPa),水平分辨率為2.5°×2.5°經(jīng)緯度,總長度為698個(gè)月(1948年1月-2006年2月)的全球格點(diǎn)再分析資料。

    3)美國航天航空局(NOAA)極軌氣象衛(wèi)星觀測的全球月平均向外長波輻射(OLR)資料。資料長度為252個(gè)月(1978年1月-1999年12月),水平分辨率為2.5°×2.5°經(jīng)緯度。

    4)南海夏季風(fēng)建立日期選用文獻(xiàn)[6-13]等確定的日期。

    1.2 方法

    本文除采用常規(guī)氣象統(tǒng)計(jì)分析方法外,主要采用SAS統(tǒng)計(jì)分析系統(tǒng)建立南海夏季風(fēng)建立日期預(yù)報(bào)方程。SAS作為全球使用最為廣泛的三大著名統(tǒng)計(jì)分析軟件(SAS,SPSS和SYSTAT)之一,是目前國際上主流的一種大型統(tǒng)計(jì)分析系統(tǒng),被譽(yù)為統(tǒng)計(jì)分析的標(biāo)準(zhǔn)軟件。SAS提供了多種先進(jìn)的回歸分析方案,結(jié)合本文實(shí)際情況采用了SAS提供的線性回歸全模型(Linear Regression-Full mode)建立南海夏季風(fēng)建立日期回歸預(yù)報(bào)方程。該方法能最大程度的運(yùn)用到預(yù)報(bào)因子的有效信息,使回歸預(yù)報(bào)方程能更加準(zhǔn)確的反映出南海夏季風(fēng)建立日期的變化特征,為預(yù)報(bào)奠定良好的基礎(chǔ)。

    2 南海夏季風(fēng)預(yù)報(bào)因子的建立

    利用ERA 和NCEP兩套資料分析南海夏季風(fēng)建立日期(選用文獻(xiàn)[13]確定的建立日期)與前期各氣象要素場(u、v、ω風(fēng)場、高度場、氣溫場)之間的相關(guān)性。結(jié)果發(fā)現(xiàn),前期2月份各要素場和后期季風(fēng)的建立日期都有顯著相關(guān)區(qū)域存在。各氣象要素場和建立日期之間存在密切的聯(lián)系,反映出要素場的前期變化對后期季風(fēng)爆發(fā)有一定的指示能力和預(yù)報(bào)作用。因此,定義相關(guān)系數(shù)超過0.3(r>rα=0.304 4,α=0.05)的區(qū)域?yàn)轭A(yù)報(bào)關(guān)鍵區(qū),并把關(guān)鍵區(qū)內(nèi)氣象要素的平均值做為預(yù)報(bào)因子,即預(yù)報(bào)因子必須是預(yù)報(bào)關(guān)鍵區(qū)內(nèi)氣象要素的平均值。如果存在多個(gè)預(yù)報(bào)關(guān)鍵區(qū),則采用線性組合方式將各個(gè)關(guān)鍵區(qū)內(nèi)要素平均值組合成一個(gè)預(yù)報(bào)因子。預(yù)報(bào)關(guān)鍵區(qū)在水平方向和垂直方向都必須足夠廣闊和深厚,具體定義原則如下:①預(yù)報(bào)關(guān)鍵區(qū)是相關(guān)系數(shù)全部超過0.3(r>rα=0.304 4,α=0.05)的區(qū)域;②同一個(gè)預(yù)報(bào)關(guān)鍵區(qū)要同時(shí)出現(xiàn)在兩個(gè)以上的氣壓層上;③同一個(gè)預(yù)報(bào)關(guān)鍵區(qū)的經(jīng)度和緯度跨度都必須大于10°;④同一個(gè)預(yù)報(bào)關(guān)鍵區(qū)最好具有在時(shí)間尺度上的延續(xù)性,即在連續(xù)幾個(gè)月中都有出現(xiàn)。

    根據(jù)上述方式選出10個(gè)預(yù)報(bào)因子,具體定義方法如下:

    1) 印度洋地區(qū)(50°-90°E,10°S-10°N)和太平洋地區(qū)(160°E-160°W,10°S-10°N)的顯著相關(guān)區(qū)域內(nèi)的低層(1 000~700 hPa)和高層(200~100 hPa)平均緯向風(fēng)速的差值定義為預(yù)報(bào)因子F1:

    F1=u1-u2

    其中,ua1為2月份50-90°E,10°S-10°N區(qū)域內(nèi)1000~700 hPa平均的緯向風(fēng)速,ua2為2月份160°E-160°W,10°S-10°N區(qū)域內(nèi)1000~700 hPa平均的緯向風(fēng)速,ub1為2月份50°E-90°E,10°S-10°N區(qū)域內(nèi)200~100 hPa平均的緯向風(fēng)速,ub2為2月份160°E-160°W,10°S-10°N區(qū)域內(nèi)200~100 hPa平均的緯向風(fēng)速。u1為ua1與ua2之差,u2為ub1與ub2之差。

    2) 南太平洋地區(qū)的顯著關(guān)鍵1區(qū)(140 -170°E,50°-35°S)和2區(qū) (145°-165°E,27.5-20°S)的整個(gè)對流層(1 000~100 hPa)平均緯向風(fēng)之差定義為預(yù)報(bào)因子F2:

    F2=u1-u2

    其中,u1為2月份140°-170°E,50°-35°S區(qū)域內(nèi)1 000~100 hPa平均的緯向風(fēng)速,u2為2月份145°-165°E,27.5°-20°S區(qū)域內(nèi)1 000~100 hPa平均的緯向風(fēng)速。

    3) 130°-145°E,40°-20°S區(qū)域和155°-175°E,40°-20°S區(qū)域?qū)α鲗?1 000~100 hPa)的平均經(jīng)向風(fēng)速之差定義為預(yù)報(bào)因子F3:

    F3=v1-v2

    其中,v1為2月份130°-145°E,40°-20°S區(qū)域1 000~100 hPa平均的經(jīng)向風(fēng)速,v2為2月份155°-180°E,40°-20°S區(qū)域1 000~100 hPa平均的經(jīng)向風(fēng)速。

    4) 40°-60°E,5°-20°N,1 000~700 hPa區(qū)域的平均經(jīng)向風(fēng)速定義為預(yù)報(bào)因子F4:

    F4=vave

    5) 140°E-170°W,5°-25°N,250~100 hPa區(qū)域的平均經(jīng)向風(fēng)速定義為預(yù)報(bào)因子F5:

    F5=vave

    6) 印度洋區(qū)域(50°-100°E,5°S-25°N)和太平洋區(qū)域(180°-80°W,15°S-10°N)整個(gè)對流層(1 000~100 hPa)的平均垂直風(fēng)速ω之差定義為預(yù)報(bào)因子F6:

    F6= ω1- ω2

    其中,ω1為2月份50°-100°E,5°-25°S區(qū)域1 000~100 hPa平均的垂直風(fēng)速,ω2為2月份180°-80°W,15°S-10°N區(qū)域1 000~100 hPa平均的垂直風(fēng)速。

    7)50°-120°E,15°S-25°N區(qū)域的對流層低層(1 000~700 hPa)的平均位勢高度定義為預(yù)報(bào)因子F7:

    F7=zave

    8)南半球區(qū)域1(145°-165°E,35°-25°S)和區(qū)域2(155°-180°E,60°-50°S)對流層的平均高度之差定義為預(yù)報(bào)因子F8:

    F8=z1-z2

    其中,z1為2月份145°-165°E,35°-25°S區(qū)域1 000~200 hPa平均的位勢高度值,z2為2月份155°-180°E,60°-50°S區(qū)域1 000~200 hPa平均的位勢高度值。

    9)120°-150°E,2.5°-15°N區(qū)域的平均氣溫定義為預(yù)報(bào)因子F9:

    F9=Tave

    其中,Tave是120°-150°E,2.5°-15°N區(qū)域1 000~850 hPa的平均氣溫。

    10)亞歐大陸地區(qū)(25°-120°E,10°-45°N)和太平洋地區(qū)(180°-40°W,30°S-40°N)對流層中上層(500~200 hPa)的平均氣溫之差定義為預(yù)報(bào)因子F10:

    F10=T1-T2

    其中,T1為2月份25°-120°E,10°-45°N區(qū)域500~200 hPa多層平均的氣溫,T2為2月份180°-40°W,30°S-40°N區(qū)域500~200 hPa多層平均的氣溫。

    根據(jù)以上方法定義的10個(gè)預(yù)報(bào)因子主要來自于5個(gè)氣象要素場,即u風(fēng)場(F1、F2)、v風(fēng)場(F3、F4、F5)、ω風(fēng)場(F6)、位勢高度場(F7,F(xiàn)8)、氣溫場(F9、F10)。10個(gè)預(yù)報(bào)因子能較全面的反映出對南海夏季風(fēng)建立日期有預(yù)報(bào)作用的關(guān)鍵區(qū),能較好的代表與南海夏季風(fēng)建立日期有密切關(guān)系的2月份各要素場的異常變化,為研究前期氣象要素場的異常對南海夏季風(fēng)建立日期的影響作用提供定量分析的途徑。分析預(yù)報(bào)因子與南海夏季風(fēng)爆發(fā)遲早的相關(guān)性,結(jié)果如表1,可見用ERA資料計(jì)算的10個(gè)預(yù)報(bào)因子與不同學(xué)者定義的南海夏季風(fēng)建立日期之間的相關(guān)系數(shù)幾乎都超過0.3(r>rα=0.304 4,α=0.05),達(dá)到95%的置信度,表明這10個(gè)預(yù)報(bào)因子不僅與文獻(xiàn)[13]定義的建立日期存在顯著相關(guān)關(guān)系,而且與其他學(xué)者的建立日期也存在較好的相關(guān)性,說明這10個(gè)預(yù)報(bào)因子對南海夏季風(fēng)建立日期有較好的指示作用,且具有作為預(yù)報(bào)因子的可能。

    表1 10個(gè)預(yù)報(bào)因子分別與各學(xué)者確定的南海夏季風(fēng)建立日期之間的相關(guān)系數(shù)(1958-1998年)

    *表示相關(guān)系數(shù)未達(dá)到95%的置信度

    3 南海夏季風(fēng)建立日期的預(yù)報(bào)

    3.1 預(yù)報(bào)方程的建立

    上述分析所用的ERA資料時(shí)間長度只到2001年12月,無法計(jì)算2001年以后的預(yù)報(bào)因子,所以下面將用NCEP資料(1948年1月-2006年2月)按照前面的方法計(jì)算10個(gè)預(yù)報(bào)因子,得到的預(yù)報(bào)因子時(shí)間序列為1948-2006年,樣本數(shù)為59個(gè)。使用相關(guān)分析對兩套再分析資料計(jì)算得到的預(yù)報(bào)因子進(jìn)行對比分析,發(fā)現(xiàn)2套預(yù)報(bào)因子之間的相關(guān)系數(shù)超過99%的置信度,2套預(yù)報(bào)因子的變化特征基本相同,演變趨勢相一致。

    用10個(gè)預(yù)報(bào)因子作為回歸因子,采用SAS提供的線性回歸全模型建立南海夏季風(fēng)建立日期回歸預(yù)報(bào)方程。經(jīng)過計(jì)算得出結(jié)果,具體分析如下:

    (1) 方差分析:回歸平方和=3 389.265 50,誤差平方和=1 396.344 25,總平方和=4 785.609 76;三者的自由度分別為10、30、40?;貧w平方和的均值=338.626 55,誤差平方和的均值=46.544 81。實(shí)際的F值就=7.28,大于Fα(0.05,10,30)=2.16,回歸方程顯著。

    (2) 模型檢驗(yàn)P值=0.000 1,小于顯著水平α=0.05,說明模型擬合合理,不必修改回歸模型。

    (3) 判定系數(shù)R2=0.708 2。說明建立日期的變化有71%是來自前期因子的變化。因此,用這10個(gè)因子能很好的反映南海夏季風(fēng)建立日期。

    (4) 回歸方程:

    Y=203.642 84-0.302 01X1-0.547 38X2-1.514 77X3+2.462 09X4-0.312 57X5+528.119 86X6-0.172 86X7-0.093 54X8-0.410 77X9-1.828 95X10

    3.2 南海夏季風(fēng)建立日期的擬合預(yù)報(bào)

    用前述的回歸預(yù)報(bào)方程回報(bào)1948-2006年南海夏季風(fēng)建立日期(見表2)。擬合的季風(fēng)建立日期平均值為5月第4候,其中最早為4月第5候(1985年),最晚為6月第1候(1968年)。若定義早于5月第2候?yàn)榧撅L(fēng)建立早年,共有8年為早年,分別是:1964、1966、1984、1985、1986、1994、1996、1999;晚于5月6候?yàn)榧撅L(fēng)建立晚年,共有10年,分別是:1953、1959、1968、1973、1982、1983、1987、1991、1992、1993 (有下劃線的年份為公認(rèn)的季風(fēng)爆發(fā)異常年)。分析發(fā)現(xiàn),原始建立日期和擬合建立日期之間關(guān)于季風(fēng)爆發(fā)異常年的確定存在一定的差異,如1964年在擬合日期中為偏早年,而在原始建立日期中則為正常年份。這種差異可能與兩序列對異常的定義標(biāo)準(zhǔn)不同有關(guān)系。原始建立日期以5月第3候和第5候?yàn)楸l(fā)早、晚年的閾值,而擬合日期則以5月第2候和第6候?yàn)榉纸缇€。圖1為原始建立日期和擬合日期的年際變化曲線,兩序列的相關(guān)系數(shù)高達(dá)0.73,說明擬合日期可以較好的反映出南海夏季風(fēng)建立的變化特征。圖2是原始日期和擬合日期差值的絕對值的時(shí)間序列變化曲線。從圖2可看,兩日期之間最大相差7候,相差0~1候的情況占64.9%,反映出擬合日期不僅能很好的刻畫原始日期的變化規(guī)律,還能較好較準(zhǔn)確的預(yù)報(bào)季風(fēng)建立的具體日期,且精確到候。分析擬合日期與各種南海夏季風(fēng)建立日期的相關(guān)性發(fā)現(xiàn),所有的相關(guān)系數(shù)均超過0.4,即達(dá)到99%的置信度,而且有6/8的相關(guān)系數(shù)達(dá)到0.6以上,說明回歸預(yù)報(bào)方程能較好的擬合季風(fēng)建立日期,同時(shí)也具有較好的預(yù)測效果。

    表2 1948-2006年南海夏季風(fēng)建立日期1)

    1)5.5表示5月第5候,其他依此類推

    圖1 1950-1999年標(biāo)準(zhǔn)化的南海夏季風(fēng)建立日期和擬合日期

    圖2 南海夏季風(fēng)建立日期與擬合日期的差值時(shí)間序列(絕對值)

    從圖2可見,原始日期和擬合日期之間相差3候以上的年份共有6 a(1951、1953、1957、1970、2001、2002),其中有4年(1951、1953、1957、1970)的建立日期本身存在較大的爭議(不同學(xué)者確定的建立日期存在較大差異)。以1970年為例,文獻(xiàn)[6,7,9-11]確定的南海夏季風(fēng)建立日期為5月第6候以后,均認(rèn)為1970年是季風(fēng)建立晚年,文獻(xiàn)[8,12-13]卻認(rèn)為1970年為季風(fēng)建立早年,說明各文獻(xiàn)對1970年南海夏季風(fēng)建立日期沒有統(tǒng)一的定論且存在較大的爭議,根據(jù)何金海等[14]的定義1970年季風(fēng)建立日期不存在優(yōu)勢候。除此之外,文獻(xiàn)[15-16]認(rèn)為2002年南海夏季風(fēng)在5月15日爆發(fā),這一結(jié)果與梁建茵等[13]的結(jié)果也存在較大的差異。如果以梁建茵等[13]的建立日期作為南海夏季風(fēng)建立日期的標(biāo)準(zhǔn)值,前面建立的回歸預(yù)報(bào)方程對有爭議的季風(fēng)建立年份的擬合效果則不理想,這可能是有兩方面的原因造成的。一方面,回歸方案本身不夠完善,存在系統(tǒng)性的誤差;另一方面,現(xiàn)有的季風(fēng)建立日期還不能完全反映南海夏季風(fēng)的爆發(fā)早晚,原始建立日期本身存在的不足直接造成擬合結(jié)果的不準(zhǔn)確。

    圖2中反映出2001年原始建立日期和擬合日期相差7候,為擬合效果最差年份。下面將具體分析2001年南海夏季風(fēng)爆發(fā)的氣候特征,從而進(jìn)一步訂正季風(fēng)建立的具體日期。

    從2001年5月第3候850 hPa環(huán)流場圖(圖3a)上見,西太平洋副熱帶高壓已經(jīng)退出南海海域,但該地區(qū)仍然由氣旋性環(huán)流控制,盛行東北風(fēng)。但在5月第4候(圖3b),整個(gè)南海區(qū)域由西南風(fēng)控制。西南氣流主要來源于孟加拉灣南部的低緯熱帶地區(qū),西太平洋副熱帶高壓脊此時(shí)已經(jīng)退到菲律賓以東洋面,索馬里急流顯著加強(qiáng),整個(gè)北印度洋為強(qiáng)西風(fēng)帶,南海夏季風(fēng)爆發(fā)形勢已經(jīng)存在。從5月第6候(圖3c)850hPa流場圖可見,夏季風(fēng)的環(huán)流形勢在第5候的基礎(chǔ)上更加明顯。

    圖3 2001年5月第3(a)、4(b)、5(c)候850 hPa風(fēng)場

    分析2001年南海地區(qū)(10°-20°N,110°-120°E)850 hPa風(fēng)場u、v分量以及OLR資料的候平均變化曲線(圖略)可知,南海地區(qū)第一次由偏東風(fēng)轉(zhuǎn)為偏西風(fēng)的日期為5月第1候。雖然南海地區(qū)早在4月份就出現(xiàn)由偏北風(fēng)轉(zhuǎn)為偏南風(fēng)的現(xiàn)象,但一直到5月第3候偏南風(fēng)才穩(wěn)定的出現(xiàn),這時(shí)整個(gè)南海地區(qū)才以偏南風(fēng)為主。南海地區(qū)的平均OLR值在5月第4候才超過235 W/m2,5月第5候后就一直穩(wěn)定在高值區(qū),到9月份才開始下降。從以上3種描述南海季風(fēng)建立的主要指標(biāo)來看,垂直對流指標(biāo)(OLR小于235 W/m2)和動(dòng)力指標(biāo)(u、v值大于0,以偏南風(fēng)和偏西風(fēng)為主)同時(shí)達(dá)到季風(fēng)建立條件的日期是5月第5候,說明2001年南海夏季風(fēng)建立日期應(yīng)該確定為5月第5候,而擬合日期為5月第6候,兩者雖然相差1候,但都指示出2001年南海夏季風(fēng)建立日期偏晚,屬于偏晚年。

    綜上所述,回歸預(yù)報(bào)方程能準(zhǔn)確預(yù)報(bào)出南海夏季風(fēng)平均建立日期為5月第4候,這與現(xiàn)有大多數(shù)研究成果相吻合。除此之外,通過分析2001年季風(fēng)爆發(fā)的特征可以知道,預(yù)報(bào)方程對季風(fēng)爆發(fā)遲早存在爭議的年份的預(yù)報(bào)也較好,南海夏季風(fēng)建立日期回歸預(yù)報(bào)方程對季風(fēng)建立日期的擬合和預(yù)測較好,能刻畫出建立日期的年際變化特征,也有較好的預(yù)測能力。根據(jù)NCEP資料,采用回歸預(yù)報(bào)方程得到的2006年南海夏季風(fēng)建立日期為5月第4候,這一結(jié)果與客觀觀測事實(shí)基本相符。

    4 結(jié) 論

    1)分析發(fā)現(xiàn),前期2月份風(fēng)場(u、v、ω分量)、位勢高度場、氣溫場等氣象要素場與后期南海夏季風(fēng)建立日期存在密切的聯(lián)系。選取要素場與建立日期的典型顯著相關(guān)區(qū)為建立日期預(yù)報(bào)關(guān)鍵區(qū),定義關(guān)鍵區(qū)中的要素平均值為預(yù)報(bào)因子,并從中選取10個(gè)對季風(fēng)建立日期指示能力和預(yù)測作用最好的預(yù)報(bào)因子建立預(yù)測方程。選出來的10個(gè)預(yù)報(bào)因子與大多數(shù)學(xué)者確定的季風(fēng)建立日期都存在顯著相關(guān)關(guān)系。

    2)采用SAS統(tǒng)計(jì)分析系統(tǒng)分析和計(jì)算了10個(gè)預(yù)報(bào)因子和建立日期的統(tǒng)計(jì)關(guān)系,并在分析結(jié)果的基礎(chǔ)上選取最好的回歸擬合方案建立南海夏季風(fēng)建立日期的回歸預(yù)報(bào)方程。對預(yù)報(bào)方程進(jìn)行回報(bào)檢驗(yàn)發(fā)現(xiàn),回歸預(yù)報(bào)方程得到的擬合日期與各位學(xué)者確定的各種季風(fēng)建立日期的相關(guān)系數(shù)都超過0.4,達(dá)到99%以上的置信度,其中大多數(shù)相關(guān)系數(shù)高達(dá)0.6,說明該回歸預(yù)報(bào)方程能較好的刻畫出季風(fēng)建立日期的變化特征。預(yù)報(bào)方程得到的南海夏季風(fēng)平均建立日期為5月第4候,這與現(xiàn)有的大多研究工作的結(jié)果相同。對比分析擬合日期和原始季風(fēng)建立日期發(fā)現(xiàn),回歸預(yù)報(bào)方程能較準(zhǔn)確的預(yù)報(bào)出季風(fēng)建立的具體日期,特別是對公認(rèn)的南海夏季風(fēng)爆發(fā)異常年份的預(yù)報(bào)十分準(zhǔn)確,但對季風(fēng)爆發(fā)遲早還存在較大爭議的年份的預(yù)報(bào)效果較差。這可能是由于各學(xué)者對季風(fēng)建立日期本身存在明顯的爭議而無法準(zhǔn)確確定所造成。通過對建立日期存在爭議的季風(fēng)爆發(fā)年份的各種季風(fēng)建立指標(biāo)以及環(huán)流特征的分析,發(fā)現(xiàn)回歸預(yù)報(bào)方程對這些年份的預(yù)報(bào)誤差較小,如2001年季風(fēng)建立日期為5月第5候,預(yù)報(bào)方程的預(yù)測結(jié)果為5月第6候,兩者雖然相差1候,卻都反映出當(dāng)年南海夏季風(fēng)爆發(fā)屬于偏晚年。

    回歸預(yù)報(bào)方程作為南海夏季風(fēng)建立日期預(yù)報(bào)方法的一種嘗試,雖然能較好的預(yù)報(bào)出季風(fēng)建立日期,但也存在一些問題。如無法清晰的描述出各個(gè)預(yù)報(bào)因子與季風(fēng)爆發(fā)遲早的物理聯(lián)系,這些都需要在日后的工作進(jìn)一步研究分析。

    參考文獻(xiàn):

    [1] 丁一匯,馬鶴年. 東亞季風(fēng)研究的現(xiàn)狀[M]∥何金海,丁一匯,陳隆勛. 亞洲季風(fēng)研究的新進(jìn)展.北京:氣象出版社,1996:1-14.

    [2] 孫淑清,陳雋.異常東亞冬季風(fēng)對夏季南海地區(qū)風(fēng)場及熱力場的影響[J]. 氣候與環(huán)境研究,2000,5(4):400-416.

    [3] 藍(lán)光東,溫之平,賀海晏. 南海夏季風(fēng)爆發(fā)的大氣熱源特征及其爆發(fā)遲早原因的探討[J]. 熱帶氣象學(xué)報(bào),2004,20(3): 271-280.

    [4] 周文,溫之平,陳創(chuàng)買. 南海西南季風(fēng)爆發(fā)的預(yù)測研究[J]. 中山大學(xué)學(xué)報(bào)(自然科學(xué)版),2002,41(3):95-98.

    [5] GIBSON J K, KALLBERG P,UPPALA S, et al. ECMWF Re-Analysis project report series 1[R]. ERA-15 description, version 2 R,1999:73.

    [6] 戴念軍,謝安,張勇. 南海夏季風(fēng)活動(dòng)的年際和年代際特征[J].氣候與環(huán)境研究,2000,5(4):363-374.

    [7] 賀海晏,溫之平,簡茂球. 近50年南海熱帶季風(fēng)建立遲早的氣候特征[M]//何金海,丁一匯,高輝. 南海夏季風(fēng)建立日期的確定與季風(fēng)指數(shù). 北京:氣象出版社,2001:49-54.

    [8] 馮瑞權(quán),王安宇,吳池勝. 南海夏季風(fēng)建立的氣候特征I.40年平均[J]. 熱帶氣象學(xué)報(bào),2001,17(4):345-354.

    [9] 李崇銀,張利平. 南海夏季風(fēng)特征及其指數(shù)[J]. 自然科學(xué)進(jìn)展,1999,9(6):536-541.

    [10] 姚永紅,錢永甫.用濕位渦定義的南海西南季風(fēng)指數(shù)及其與我國區(qū)域降水的關(guān)系研究[J]. 南京大學(xué)學(xué)報(bào):自然科學(xué),2001,37(6):781-788.

    [11] 高輝,何金海,譚言科,等. 40a南海夏季風(fēng)建立日期的確定[J]. 南京氣象學(xué)院學(xué)報(bào),2001,24(3):379-383.

    [12] 張秀芝,李江龍,閆俊岳,等. 南海夏季風(fēng)爆發(fā)的環(huán)流特征及指標(biāo)研究[J].氣候與環(huán)境研究,2002,7(3):321-331.

    [13] 梁建茵,吳尚森. 南海西南季風(fēng)爆發(fā)日期及其影響因子[J]. 大氣科學(xué),2002,26(6):829-844.

    [14] 何金海,丁一匯,高輝,等. 南海夏季風(fēng)建立日期的確定與季風(fēng)指數(shù)[M].北京:氣象出版社,2001:118-119.

    [15] 閆俊岳,唐志毅,姚華棟,等.2002年南海季風(fēng)建立及其雨帶變化的天氣學(xué)研究[J].氣象學(xué)報(bào),2003,61(5): 569-579.

    [16] 閆俊岳,唐志毅,姚華棟,等.2002年南海季風(fēng)爆發(fā)前后西沙海區(qū)海-氣通量交換及其變化[J].地球物理學(xué)報(bào),2005,48(5): 1000-1010.

    猜你喜歡
    季風(fēng)南海風(fēng)速
    南海明珠
    北海北、南海南
    黃河之聲(2021年10期)2021-09-18 03:07:18
    基于Kmeans-VMD-LSTM的短期風(fēng)速預(yù)測
    基于最優(yōu)TS評分和頻率匹配的江蘇近海風(fēng)速訂正
    綠水青山圖——海洋季風(fēng)的贊歌
    戶撒刀
    南海的虎斑貝
    基于GARCH的短時(shí)風(fēng)速預(yù)測方法
    南海隨筆
    草原(2016年1期)2016-01-31 21:21:51
    萬馬奔騰
    中華奇石(2015年9期)2015-07-09 18:32:58
    亚洲成av人片免费观看| 久久草成人影院| 夜夜夜夜夜久久久久| 色尼玛亚洲综合影院| 美女被艹到高潮喷水动态| 免费av不卡在线播放| 此物有八面人人有两片| 国产高清视频在线观看网站| 十八禁人妻一区二区| 国内精品久久久久久久电影| 99久久精品一区二区三区| 亚洲欧美精品综合一区二区三区| 亚洲国产看品久久| 免费大片18禁| 亚洲人成伊人成综合网2020| 99久久99久久久精品蜜桃| 亚洲一区二区三区色噜噜| 午夜福利欧美成人| 国产成人精品久久二区二区91| 成人高潮视频无遮挡免费网站| 好看av亚洲va欧美ⅴa在| 每晚都被弄得嗷嗷叫到高潮| 日韩三级视频一区二区三区| 亚洲人成电影免费在线| 国产精品野战在线观看| 日韩精品青青久久久久久| 国产成年人精品一区二区| 欧美成人一区二区免费高清观看 | 国产精品98久久久久久宅男小说| 男女下面进入的视频免费午夜| 国产精品精品国产色婷婷| av女优亚洲男人天堂 | 色综合站精品国产| 午夜激情欧美在线| 欧美乱码精品一区二区三区| 亚洲乱码一区二区免费版| 亚洲精品美女久久久久99蜜臀| 中国美女看黄片| 色综合欧美亚洲国产小说| 特级一级黄色大片| av欧美777| 18禁美女被吸乳视频| 亚洲欧美激情综合另类| 香蕉久久夜色| 丁香欧美五月| cao死你这个sao货| 黄色丝袜av网址大全| 欧美日韩瑟瑟在线播放| 女人高潮潮喷娇喘18禁视频| 一个人免费在线观看电影 | 日本 av在线| 在线观看美女被高潮喷水网站 | 久久久精品欧美日韩精品| 国产精品美女特级片免费视频播放器 | 脱女人内裤的视频| 999久久久国产精品视频| 美女高潮喷水抽搐中文字幕| 真人一进一出gif抽搐免费| 国产 一区 欧美 日韩| 免费av不卡在线播放| avwww免费| 精品免费久久久久久久清纯| 蜜桃久久精品国产亚洲av| 亚洲最大成人中文| 99久久精品热视频| 免费观看的影片在线观看| 午夜久久久久精精品| 亚洲精品一卡2卡三卡4卡5卡| 色噜噜av男人的天堂激情| 国产成人精品久久二区二区91| 欧美中文日本在线观看视频| 少妇熟女aⅴ在线视频| 亚洲自拍偷在线| 香蕉丝袜av| 亚洲精华国产精华精| 色综合欧美亚洲国产小说| 无限看片的www在线观看| 国内少妇人妻偷人精品xxx网站 | 少妇熟女aⅴ在线视频| 此物有八面人人有两片| 欧美在线黄色| 舔av片在线| 18美女黄网站色大片免费观看| 免费人成视频x8x8入口观看| bbb黄色大片| 成人三级做爰电影| 亚洲自偷自拍图片 自拍| 国产一级毛片七仙女欲春2| 在线视频色国产色| 99国产综合亚洲精品| 日本 欧美在线| 在线观看免费午夜福利视频| 91老司机精品| 又紧又爽又黄一区二区| 99国产极品粉嫩在线观看| 国产主播在线观看一区二区| 听说在线观看完整版免费高清| 国产成人精品无人区| 亚洲成人久久爱视频| 蜜桃久久精品国产亚洲av| 99热精品在线国产| 中文字幕av在线有码专区| 中文字幕最新亚洲高清| 久久久久精品国产欧美久久久| 久久精品国产99精品国产亚洲性色| 99在线视频只有这里精品首页| 又大又爽又粗| 国产精品 国内视频| 国产黄色小视频在线观看| 亚洲无线在线观看| 色综合亚洲欧美另类图片| 国产私拍福利视频在线观看| 国产精品综合久久久久久久免费| 久久精品91蜜桃| 香蕉久久夜色| 在线看三级毛片| 毛片女人毛片| 免费高清视频大片| 成在线人永久免费视频| 欧美性猛交╳xxx乱大交人| 99国产精品一区二区蜜桃av| 在线观看午夜福利视频| 成熟少妇高潮喷水视频| 看免费av毛片| 一级毛片女人18水好多| 免费看光身美女| 在线观看66精品国产| 好男人电影高清在线观看| 亚洲美女黄片视频| 欧美日韩中文字幕国产精品一区二区三区| 国产精品美女特级片免费视频播放器 | 精品99又大又爽又粗少妇毛片 | 两个人的视频大全免费| 色综合亚洲欧美另类图片| 麻豆av在线久日| 色视频www国产| 国产伦精品一区二区三区四那| 熟妇人妻久久中文字幕3abv| 欧美又色又爽又黄视频| 长腿黑丝高跟| 热99re8久久精品国产| 中文在线观看免费www的网站| 亚洲无线观看免费| 91字幕亚洲| 亚洲,欧美精品.| 亚洲五月婷婷丁香| 欧美国产日韩亚洲一区| 免费无遮挡裸体视频| 国产高清视频在线播放一区| 俄罗斯特黄特色一大片| 狂野欧美白嫩少妇大欣赏| 中文在线观看免费www的网站| 在线十欧美十亚洲十日本专区| 亚洲国产精品sss在线观看| 琪琪午夜伦伦电影理论片6080| 久久久久国内视频| 神马国产精品三级电影在线观看| 嫩草影院入口| 国产精品久久电影中文字幕| 99国产精品一区二区三区| 欧美在线黄色| 精品免费久久久久久久清纯| 女人被狂操c到高潮| 国产av麻豆久久久久久久| 99视频精品全部免费 在线 | 在线看三级毛片| 香蕉久久夜色| 人人妻人人澡欧美一区二区| 亚洲熟妇熟女久久| 女警被强在线播放| 国产av麻豆久久久久久久| 精品一区二区三区视频在线观看免费| 国产高清激情床上av| 淫秽高清视频在线观看| 熟女人妻精品中文字幕| 可以在线观看毛片的网站| 女同久久另类99精品国产91| 久久99热这里只有精品18| 在线观看美女被高潮喷水网站 | 亚洲专区中文字幕在线| 九色成人免费人妻av| 美女cb高潮喷水在线观看 | 曰老女人黄片| 一个人看视频在线观看www免费 | 丁香欧美五月| 亚洲一区二区三区不卡视频| 好男人电影高清在线观看| 老鸭窝网址在线观看| 日韩大尺度精品在线看网址| 男女床上黄色一级片免费看| av欧美777| 亚洲美女黄片视频| 美女午夜性视频免费| 国产欧美日韩精品亚洲av| 亚洲av电影在线进入| 黄频高清免费视频| 变态另类丝袜制服| 国产精品一及| 男女之事视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美成人一区二区免费高清观看 | 日韩国内少妇激情av| 午夜两性在线视频| 精品国产亚洲在线| 大型黄色视频在线免费观看| 99久久精品热视频| 国产97色在线日韩免费| 悠悠久久av| 国产成年人精品一区二区| 69av精品久久久久久| 99热精品在线国产| 1000部很黄的大片| 观看美女的网站| 午夜福利视频1000在线观看| 日本三级黄在线观看| 欧美一区二区精品小视频在线| 亚洲中文av在线| x7x7x7水蜜桃| 老司机在亚洲福利影院| 宅男免费午夜| 毛片女人毛片| 国产精品一区二区精品视频观看| 欧美黑人欧美精品刺激| 精品日产1卡2卡| 又黄又爽又免费观看的视频| 九九久久精品国产亚洲av麻豆 | 欧美3d第一页| 国产精品98久久久久久宅男小说| 九九热线精品视视频播放| 日本撒尿小便嘘嘘汇集6| 国产探花在线观看一区二区| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站高清观看| 国内揄拍国产精品人妻在线| 国产亚洲av嫩草精品影院| 夜夜躁狠狠躁天天躁| 日本黄色视频三级网站网址| 热99在线观看视频| 男人和女人高潮做爰伦理| 黑人欧美特级aaaaaa片| 国产精品久久久久久亚洲av鲁大| 国产高清有码在线观看视频| 亚洲一区高清亚洲精品| 老熟妇仑乱视频hdxx| 97超视频在线观看视频| 黄色片一级片一级黄色片| 欧美乱码精品一区二区三区| 国产不卡一卡二| 老鸭窝网址在线观看| 亚洲欧美精品综合久久99| 国产精品一区二区精品视频观看| 女生性感内裤真人,穿戴方法视频| 久久久久性生活片| 国产激情欧美一区二区| 高潮久久久久久久久久久不卡| 欧美日韩国产亚洲二区| 欧美国产日韩亚洲一区| 级片在线观看| 在线观看日韩欧美| 男插女下体视频免费在线播放| 午夜福利在线观看吧| 白带黄色成豆腐渣| 亚洲人成电影免费在线| 欧美中文综合在线视频| 国产激情久久老熟女| 成年人黄色毛片网站| 夜夜夜夜夜久久久久| 啦啦啦观看免费观看视频高清| 女警被强在线播放| 国产精品女同一区二区软件 | 不卡一级毛片| www.自偷自拍.com| 变态另类丝袜制服| 日韩欧美 国产精品| av在线蜜桃| 此物有八面人人有两片| 亚洲成人免费电影在线观看| 老汉色∧v一级毛片| 老汉色∧v一级毛片| 黄色 视频免费看| 精品国产超薄肉色丝袜足j| 国产av麻豆久久久久久久| 国产伦在线观看视频一区| 一区二区三区激情视频| 久久国产乱子伦精品免费另类| www日本黄色视频网| 久久久水蜜桃国产精品网| 狠狠狠狠99中文字幕| 青草久久国产| 色视频www国产| av黄色大香蕉| 真人做人爱边吃奶动态| 小蜜桃在线观看免费完整版高清| 激情在线观看视频在线高清| 两个人的视频大全免费| 国产精品99久久99久久久不卡| 一级毛片女人18水好多| 欧美三级亚洲精品| 最好的美女福利视频网| av欧美777| 九九热线精品视视频播放| 成在线人永久免费视频| 国产精品 国内视频| 少妇的丰满在线观看| 美女被艹到高潮喷水动态| 国产精品1区2区在线观看.| 午夜福利18| 最近视频中文字幕2019在线8| 1024手机看黄色片| 天堂影院成人在线观看| 国产精品1区2区在线观看.| av国产免费在线观看| 少妇裸体淫交视频免费看高清| 久99久视频精品免费| 91老司机精品| 国产毛片a区久久久久| 日韩欧美精品v在线| 叶爱在线成人免费视频播放| 校园春色视频在线观看| 亚洲av美国av| 久久午夜亚洲精品久久| 窝窝影院91人妻| 日韩欧美国产一区二区入口| 两个人视频免费观看高清| 精品久久久久久成人av| av天堂中文字幕网| 国产一区二区激情短视频| 精品午夜福利视频在线观看一区| 午夜福利高清视频| 成人无遮挡网站| 午夜日韩欧美国产| 免费av毛片视频| 亚洲人成伊人成综合网2020| 精品福利观看| 国产精品爽爽va在线观看网站| 午夜福利成人在线免费观看| 亚洲中文字幕日韩| 亚洲国产精品成人综合色| 中文字幕高清在线视频| 国产野战对白在线观看| 母亲3免费完整高清在线观看| 久久久久精品国产欧美久久久| 亚洲真实伦在线观看| 精品久久久久久久末码| 黄色日韩在线| 午夜免费观看网址| 又爽又黄无遮挡网站| 高清在线国产一区| 免费在线观看日本一区| av欧美777| 给我免费播放毛片高清在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 日韩三级视频一区二区三区| 午夜福利成人在线免费观看| 伊人久久大香线蕉亚洲五| 久久久国产精品麻豆| 国产久久久一区二区三区| 老熟妇仑乱视频hdxx| 在线观看午夜福利视频| 激情在线观看视频在线高清| 亚洲av电影在线进入| 色视频www国产| 99久久久亚洲精品蜜臀av| 69av精品久久久久久| 国产成人精品久久二区二区免费| 亚洲熟妇中文字幕五十中出| 丰满人妻熟妇乱又伦精品不卡| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| 真人做人爱边吃奶动态| www.999成人在线观看| 老司机午夜十八禁免费视频| 亚洲专区国产一区二区| 香蕉丝袜av| 日韩大尺度精品在线看网址| av在线蜜桃| 国产精品久久久人人做人人爽| 一本久久中文字幕| 高清在线国产一区| 免费看光身美女| www日本黄色视频网| 91久久精品国产一区二区成人 | 亚洲黑人精品在线| 精品久久久久久成人av| 小说图片视频综合网站| 欧美黑人欧美精品刺激| 黄色片一级片一级黄色片| 欧美日韩国产亚洲二区| av欧美777| 这个男人来自地球电影免费观看| 精品久久久久久久久久久久久| 亚洲一区二区三区色噜噜| 99久久精品热视频| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器 | 日韩人妻高清精品专区| 极品教师在线免费播放| 午夜成年电影在线免费观看| 成年女人看的毛片在线观看| 女生性感内裤真人,穿戴方法视频| 又爽又黄无遮挡网站| 亚洲午夜精品一区,二区,三区| 日本精品一区二区三区蜜桃| 欧美乱妇无乱码| 欧美日韩亚洲国产一区二区在线观看| 全区人妻精品视频| 舔av片在线| 香蕉国产在线看| 精品不卡国产一区二区三区| 听说在线观看完整版免费高清| 精品久久久久久久久久久久久| 国产午夜福利久久久久久| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 夜夜看夜夜爽夜夜摸| 免费大片18禁| 国产激情欧美一区二区| 久久精品亚洲精品国产色婷小说| 免费av不卡在线播放| 久久精品国产亚洲av香蕉五月| 男女下面进入的视频免费午夜| 久久久水蜜桃国产精品网| 欧美性猛交╳xxx乱大交人| 宅男免费午夜| 天天躁日日操中文字幕| 此物有八面人人有两片| 视频区欧美日本亚洲| 在线永久观看黄色视频| 久久精品aⅴ一区二区三区四区| 五月玫瑰六月丁香| 又大又爽又粗| 亚洲一区二区三区不卡视频| 亚洲人与动物交配视频| 很黄的视频免费| 国产精品av久久久久免费| 中文字幕精品亚洲无线码一区| 两性午夜刺激爽爽歪歪视频在线观看| 五月伊人婷婷丁香| 欧美成人性av电影在线观看| 女人被狂操c到高潮| 日本三级黄在线观看| 女人被狂操c到高潮| 亚洲无线观看免费| 久久欧美精品欧美久久欧美| 精品国产三级普通话版| 免费高清视频大片| 精品无人区乱码1区二区| 一级毛片精品| 99久久99久久久精品蜜桃| 国产成人一区二区三区免费视频网站| 日本撒尿小便嘘嘘汇集6| 岛国在线观看网站| 亚洲av成人av| 性欧美人与动物交配| 夜夜看夜夜爽夜夜摸| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 人人妻,人人澡人人爽秒播| av女优亚洲男人天堂 | 日本 av在线| 欧美中文日本在线观看视频| 亚洲美女黄片视频| 男女视频在线观看网站免费| 999久久久精品免费观看国产| 久久久久国内视频| 激情在线观看视频在线高清| 青草久久国产| 精品一区二区三区视频在线观看免费| 亚洲成av人片免费观看| 最好的美女福利视频网| 女人被狂操c到高潮| 国产三级中文精品| 成人亚洲精品av一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品亚洲av一区麻豆| 99久久精品热视频| 蜜桃久久精品国产亚洲av| 好男人电影高清在线观看| 国产成人精品久久二区二区91| 久久精品亚洲精品国产色婷小说| 九九久久精品国产亚洲av麻豆 | 狂野欧美激情性xxxx| 97碰自拍视频| 欧美三级亚洲精品| 国产亚洲欧美在线一区二区| 51午夜福利影视在线观看| 精品久久久久久久久久久久久| 亚洲 国产 在线| 亚洲性夜色夜夜综合| 97超级碰碰碰精品色视频在线观看| 黄色日韩在线| 19禁男女啪啪无遮挡网站| 国产成人一区二区三区免费视频网站| 精品久久久久久久人妻蜜臀av| 男女做爰动态图高潮gif福利片| 免费在线观看成人毛片| 日韩中文字幕欧美一区二区| 久久精品综合一区二区三区| 国产69精品久久久久777片 | 色综合婷婷激情| 精品电影一区二区在线| 亚洲熟妇熟女久久| 国产69精品久久久久777片 | 色哟哟哟哟哟哟| 国产一区二区在线av高清观看| 国产久久久一区二区三区| 亚洲 国产 在线| 岛国在线免费视频观看| 国产高清三级在线| 老熟妇仑乱视频hdxx| 日本 欧美在线| 99久久精品一区二区三区| 国产精品香港三级国产av潘金莲| 成人av在线播放网站| 少妇熟女aⅴ在线视频| 19禁男女啪啪无遮挡网站| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 日韩欧美精品v在线| 黄色视频,在线免费观看| 日本一二三区视频观看| 久久久精品大字幕| 午夜福利欧美成人| 成人无遮挡网站| 99热这里只有是精品50| 亚洲精品美女久久久久99蜜臀| 真人一进一出gif抽搐免费| 悠悠久久av| 国内精品久久久久久久电影| 婷婷精品国产亚洲av在线| 午夜成年电影在线免费观看| 欧美性猛交黑人性爽| 亚洲专区国产一区二区| 美女 人体艺术 gogo| 久久亚洲精品不卡| 亚洲av中文字字幕乱码综合| 99国产精品一区二区蜜桃av| 日本免费一区二区三区高清不卡| 国产久久久一区二区三区| 久久久久久大精品| 女同久久另类99精品国产91| 色视频www国产| 久久久久久久久久黄片| 亚洲va日本ⅴa欧美va伊人久久| 青草久久国产| 美女免费视频网站| 精品国内亚洲2022精品成人| 国产成人精品久久二区二区91| 国产高清视频在线播放一区| 国产久久久一区二区三区| 亚洲av日韩精品久久久久久密| 国产乱人伦免费视频| 中文字幕精品亚洲无线码一区| 男人舔女人下体高潮全视频| 久久草成人影院| avwww免费| 成人av在线播放网站| 午夜免费成人在线视频| 97超级碰碰碰精品色视频在线观看| 91字幕亚洲| 好男人在线观看高清免费视频| 欧美精品啪啪一区二区三区| 国产高清激情床上av| 九色成人免费人妻av| 国产aⅴ精品一区二区三区波| 国内少妇人妻偷人精品xxx网站 | 男人和女人高潮做爰伦理| 两性午夜刺激爽爽歪歪视频在线观看| 午夜福利免费观看在线| 国产综合懂色| 久久国产乱子伦精品免费另类| 男插女下体视频免费在线播放| 欧美不卡视频在线免费观看| 国产三级黄色录像| 一本综合久久免费| 男女视频在线观看网站免费| 国产亚洲av高清不卡| 嫁个100分男人电影在线观看| 精品国产超薄肉色丝袜足j| 在线观看免费午夜福利视频| 久久久国产欧美日韩av| 蜜桃久久精品国产亚洲av| 国产高清激情床上av| 国产成人av激情在线播放| 极品教师在线免费播放| 一二三四在线观看免费中文在| 国产精品美女特级片免费视频播放器 | 村上凉子中文字幕在线| 日韩欧美国产在线观看| 好看av亚洲va欧美ⅴa在| 亚洲av成人av| 亚洲男人的天堂狠狠| 精品国产三级普通话版| 中文字幕精品亚洲无线码一区| 国产乱人伦免费视频| 亚洲国产欧美一区二区综合| 两个人的视频大全免费| 午夜免费激情av| 免费av不卡在线播放| 男女之事视频高清在线观看| 国产亚洲精品综合一区在线观看| 国产成人精品无人区| 熟妇人妻久久中文字幕3abv| 免费观看人在逋| 久久精品国产99精品国产亚洲性色| 亚洲精品一区av在线观看| x7x7x7水蜜桃| 久久亚洲精品不卡| 亚洲av电影不卡..在线观看| 两个人的视频大全免费| 深夜精品福利| 国产亚洲欧美98| 日本与韩国留学比较| 1024手机看黄色片|