• <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网站的网址| 蜜桃在线观看..| 涩涩av久久男人的天堂| 久久午夜福利片| 国产1区2区3区精品| 狠狠婷婷综合久久久久久88av| 男女无遮挡免费网站观看| 嫩草影院入口| 美女国产视频在线观看| 大话2 男鬼变身卡| 国产在视频线精品| 中国国产av一级| 精品久久蜜臀av无| 18在线观看网站| 人妻一区二区av| 亚洲精品国产av成人精品| 1024视频免费在线观看| 成人亚洲欧美一区二区av| 日韩在线高清观看一区二区三区| 乱人伦中国视频| 精品久久久精品久久久| 久久影院123| 草草在线视频免费看| 永久免费av网站大全| 五月开心婷婷网| 久久久精品免费免费高清| 精品亚洲成国产av| 日韩中字成人| 久久精品国产亚洲av涩爱| av在线老鸭窝| 黄色视频在线播放观看不卡| 久久久久久久精品精品| 精品午夜福利在线看| av免费观看日本| 久久精品人人爽人人爽视色| 成人午夜精彩视频在线观看| 国产高清三级在线| 精品少妇内射三级| 免费观看无遮挡的男女| 欧美精品av麻豆av| 国产成人免费观看mmmm| 国产精品国产三级专区第一集| 亚洲人成77777在线视频| 亚洲成av片中文字幕在线观看 | 一个人免费看片子| 亚洲成人手机| √禁漫天堂资源中文www| 七月丁香在线播放| 波野结衣二区三区在线| 国产高清三级在线| 亚洲国产成人一精品久久久| 精品亚洲成国产av| 秋霞在线观看毛片| 成人无遮挡网站| 色哟哟·www| 亚洲av福利一区| 一区二区三区精品91| 亚洲精品美女久久av网站| 欧美97在线视频| 国产欧美日韩一区二区三区在线| 两个人看的免费小视频| 美女国产视频在线观看| 国产精品久久久久久精品电影小说| 美女国产高潮福利片在线看| 一边摸一边做爽爽视频免费| 毛片一级片免费看久久久久| 人人澡人人妻人| 母亲3免费完整高清在线观看 | 有码 亚洲区| 深夜精品福利| 精品亚洲成a人片在线观看| 欧美成人精品欧美一级黄| 亚洲五月色婷婷综合| 一级毛片 在线播放| 亚洲av中文av极速乱| 日韩av在线免费看完整版不卡| 免费观看av网站的网址| 国产精品国产av在线观看| 亚洲欧美清纯卡通| 9色porny在线观看| www.熟女人妻精品国产 | 亚洲欧美日韩卡通动漫| 亚洲av男天堂| 岛国毛片在线播放| 麻豆精品久久久久久蜜桃| 搡老乐熟女国产| 各种免费的搞黄视频| 中文字幕人妻熟女乱码| 欧美亚洲 丝袜 人妻 在线| 一级毛片黄色毛片免费观看视频| a级片在线免费高清观看视频| 黑人巨大精品欧美一区二区蜜桃 | 9191精品国产免费久久| 少妇人妻久久综合中文| 美女主播在线视频| 18在线观看网站| 男女边摸边吃奶| 90打野战视频偷拍视频| 夫妻性生交免费视频一级片| 国产日韩欧美在线精品| 亚洲少妇的诱惑av| 超色免费av| 日本-黄色视频高清免费观看| 久久久国产欧美日韩av| 婷婷色av中文字幕| 三级国产精品片| 精品一区在线观看国产| 中文天堂在线官网| 高清毛片免费看| 建设人人有责人人尽责人人享有的| a 毛片基地| 日韩三级伦理在线观看| 国产成人精品久久久久久| 免费高清在线观看日韩| 国产精品一区二区在线观看99| 97超碰精品成人国产| 免费看不卡的av| 国产精品国产av在线观看| 极品少妇高潮喷水抽搐| 女人精品久久久久毛片| 国产精品偷伦视频观看了| 26uuu在线亚洲综合色| 欧美精品人与动牲交sv欧美| 日本免费在线观看一区| 纵有疾风起免费观看全集完整版| 制服丝袜香蕉在线| 十八禁网站网址无遮挡| 国产精品一区二区在线观看99| 欧美老熟妇乱子伦牲交| 老女人水多毛片| 成人综合一区亚洲| 丁香六月天网| 国产精品国产三级国产专区5o| 欧美日韩亚洲高清精品| 最近中文字幕高清免费大全6| 久久久久久久大尺度免费视频| 亚洲av男天堂| 国产午夜精品一二区理论片| 在线观看美女被高潮喷水网站| av免费观看日本| h视频一区二区三区| 精品视频人人做人人爽| 乱码一卡2卡4卡精品| av免费观看日本| 亚洲精品久久成人aⅴ小说| 成人亚洲精品一区在线观看| 1024视频免费在线观看| 国产成人精品无人区| 亚洲av.av天堂| 中文字幕av电影在线播放| 国产精品国产三级国产专区5o| 日韩视频在线欧美| 久久久久精品人妻al黑| 免费观看av网站的网址| 中国三级夫妇交换| 多毛熟女@视频| 亚洲成人手机| 女人被躁到高潮嗷嗷叫费观| 亚洲av电影在线进入| 亚洲精品乱码久久久久久按摩| 成人免费观看视频高清| 一级黄片播放器| 99re6热这里在线精品视频| 伊人久久国产一区二区| 亚洲美女搞黄在线观看| 国产精品蜜桃在线观看| 内地一区二区视频在线| 黑人巨大精品欧美一区二区蜜桃 | 夫妻性生交免费视频一级片| a级毛色黄片| 欧美3d第一页| 国产亚洲最大av| 国产精品人妻久久久影院| 亚洲综合色网址| 精品久久国产蜜桃| 侵犯人妻中文字幕一二三四区| 男人操女人黄网站| 天堂俺去俺来也www色官网| 激情视频va一区二区三区| 美女内射精品一级片tv| 成人综合一区亚洲| 亚洲天堂av无毛| 中文字幕另类日韩欧美亚洲嫩草| 久久av网站| 一个人免费看片子| 丝袜在线中文字幕| 五月伊人婷婷丁香| 日韩av在线免费看完整版不卡| 视频在线观看一区二区三区| 国产在线视频一区二区| av有码第一页| 成人午夜精彩视频在线观看| 亚洲少妇的诱惑av| 免费看av在线观看网站| 在线观看www视频免费| 亚洲色图综合在线观看| 亚洲av成人精品一二三区| 中文天堂在线官网| 国产视频首页在线观看| 久久久久久人人人人人| 日韩人妻精品一区2区三区| 日日啪夜夜爽| 国产 精品1| 永久网站在线| 欧美日本中文国产一区发布| 国产精品国产av在线观看| 97人妻天天添夜夜摸| 97人妻天天添夜夜摸| av免费在线看不卡| 国产午夜精品一二区理论片| 人人妻人人澡人人爽人人夜夜| 欧美日韩一区二区视频在线观看视频在线| 十分钟在线观看高清视频www| 欧美日韩成人在线一区二区| 亚洲成人一二三区av| 久久99热6这里只有精品| a 毛片基地| 国产精品99久久99久久久不卡 | 国产免费视频播放在线视频| 91久久精品国产一区二区三区| 香蕉国产在线看| 成人亚洲欧美一区二区av| 国产成人精品福利久久| 国产探花极品一区二区| 美女国产高潮福利片在线看| 国产黄色视频一区二区在线观看| 国产精品久久久久久av不卡| 国产视频首页在线观看| 人成视频在线观看免费观看| 99热全是精品| 欧美国产精品va在线观看不卡| 日本爱情动作片www.在线观看| 激情五月婷婷亚洲| 中国国产av一级| 国产一区二区在线观看av| 精品午夜福利在线看| 成人影院久久| 成人无遮挡网站| 曰老女人黄片| 精品国产国语对白av| freevideosex欧美| 女人被躁到高潮嗷嗷叫费观| 丝袜美足系列| 最近中文字幕高清免费大全6| 免费黄色在线免费观看| 日韩视频在线欧美| 九草在线视频观看| 在线天堂中文资源库| 午夜福利,免费看| 国产免费福利视频在线观看| 亚洲精品视频女| 亚洲欧美成人综合另类久久久| 九九在线视频观看精品| 捣出白浆h1v1| 亚洲中文av在线| 制服诱惑二区| 亚洲国产日韩一区二区| 在线观看免费日韩欧美大片| 亚洲内射少妇av| 国产高清国产精品国产三级| 亚洲伊人久久精品综合| xxx大片免费视频| 校园人妻丝袜中文字幕| 久久韩国三级中文字幕| 99久国产av精品国产电影| 激情视频va一区二区三区| 国产成人a∨麻豆精品| 亚洲国产最新在线播放| 22中文网久久字幕| 22中文网久久字幕| 午夜福利网站1000一区二区三区| 男人舔女人的私密视频| 乱人伦中国视频| 国产精品久久久久久精品古装| 高清欧美精品videossex| av不卡在线播放| 日本av手机在线免费观看| 色视频在线一区二区三区| 午夜老司机福利剧场| 啦啦啦啦在线视频资源| 黑丝袜美女国产一区| 国产成人精品无人区| 国产在线视频一区二区| 亚洲精品一区蜜桃| 色吧在线观看| 亚洲美女视频黄频| 少妇人妻 视频| 久久精品国产自在天天线| 少妇的逼水好多| 一区二区三区乱码不卡18| 高清黄色对白视频在线免费看| 天天躁夜夜躁狠狠躁躁| 国产乱人偷精品视频| a级片在线免费高清观看视频| 免费高清在线观看日韩| 久久久久久伊人网av| 日韩成人av中文字幕在线观看| 国产乱来视频区| 日本欧美视频一区| 亚洲第一区二区三区不卡| 亚洲激情五月婷婷啪啪| 亚洲精品乱久久久久久| 亚洲精品久久午夜乱码| 亚洲高清免费不卡视频| 少妇人妻精品综合一区二区| 久久久久久久精品精品| 亚洲内射少妇av| 1024视频免费在线观看| 成人综合一区亚洲| 精品人妻熟女毛片av久久网站| 久久99热这里只频精品6学生| 亚洲高清免费不卡视频| 日韩一本色道免费dvd| 在线精品无人区一区二区三| 亚洲在久久综合| 国产精品一区二区在线观看99| 岛国毛片在线播放| 久久久久久久国产电影| 国产成人免费观看mmmm| 亚洲中文av在线| 久久免费观看电影| 一本色道久久久久久精品综合| 久久精品久久精品一区二区三区| 欧美丝袜亚洲另类| 日本午夜av视频| 久久久久久久久久久久大奶| 国产熟女欧美一区二区| 亚洲一区二区三区欧美精品| a级片在线免费高清观看视频| 熟女av电影| 男女啪啪激烈高潮av片| 免费不卡的大黄色大毛片视频在线观看| 在线观看一区二区三区激情| 黑丝袜美女国产一区| 人成视频在线观看免费观看| 菩萨蛮人人尽说江南好唐韦庄| 久久久久久伊人网av| 国产成人a∨麻豆精品| 中文字幕人妻熟女乱码| 国产精品国产av在线观看| 满18在线观看网站| 日韩大片免费观看网站| 久久久久视频综合| 亚洲精品久久午夜乱码| 成人影院久久| 亚洲欧美日韩卡通动漫| 黄色毛片三级朝国网站| 99久久综合免费| 欧美性感艳星| 国产 一区精品| 午夜视频国产福利| 久久久国产欧美日韩av| 欧美成人午夜精品| 亚洲国产成人一精品久久久| 青春草视频在线免费观看| 老司机影院毛片| 成人亚洲精品一区在线观看| 午夜日本视频在线| 午夜91福利影院| 亚洲综合色惰| 午夜久久久在线观看| 欧美精品国产亚洲| 久久99热这里只频精品6学生| 18禁动态无遮挡网站| 大陆偷拍与自拍| 亚洲美女视频黄频| av在线播放精品| 男人舔女人的私密视频| 国产高清三级在线| 国产精品国产av在线观看| 欧美激情 高清一区二区三区| 欧美精品一区二区大全| 亚洲av男天堂| 黑人欧美特级aaaaaa片| 99久久精品国产国产毛片| 9热在线视频观看99| 男人添女人高潮全过程视频| 亚洲国产精品999| 黄片无遮挡物在线观看| av片东京热男人的天堂| 中文字幕最新亚洲高清| 在现免费观看毛片| 亚洲美女搞黄在线观看| 亚洲精品久久成人aⅴ小说| 老司机影院毛片| 丝袜人妻中文字幕| 久久久国产一区二区| av国产久精品久网站免费入址| 水蜜桃什么品种好| 精品人妻一区二区三区麻豆| 国产伦理片在线播放av一区| 中国美白少妇内射xxxbb| 亚洲精品美女久久久久99蜜臀 | 91国产中文字幕| 如何舔出高潮| 亚洲人成网站在线观看播放| 日本免费在线观看一区| 婷婷色综合大香蕉| 人人澡人人妻人| 纯流量卡能插随身wifi吗| 成人午夜精彩视频在线观看| 如日韩欧美国产精品一区二区三区| 国产男女内射视频| 成年人午夜在线观看视频| 男女午夜视频在线观看 | 亚洲精品aⅴ在线观看| 日韩中文字幕视频在线看片| 大片免费播放器 马上看| 丝袜在线中文字幕| 成人午夜精彩视频在线观看| 久久久久国产网址| 精品国产一区二区久久| 少妇猛男粗大的猛烈进出视频| 国产精品久久久久久精品古装| 在线观看免费视频网站a站| 亚洲图色成人| 日本黄色日本黄色录像| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 这个男人来自地球电影免费观看 | 人体艺术视频欧美日本| 免费在线观看完整版高清| 宅男免费午夜| 国产高清三级在线| 男女午夜视频在线观看 | 少妇高潮的动态图| 免费观看av网站的网址| 在线观看一区二区三区激情| 午夜视频国产福利| 一级,二级,三级黄色视频| 桃花免费在线播放| 99久国产av精品国产电影| 午夜精品国产一区二区电影| 亚洲av电影在线进入| 美女主播在线视频| 中文字幕精品免费在线观看视频 | 久久99蜜桃精品久久| 日本欧美视频一区| 亚洲色图综合在线观看| 视频中文字幕在线观看| 五月伊人婷婷丁香| 搡老乐熟女国产| 亚洲成国产人片在线观看| 又黄又爽又刺激的免费视频.| 2021少妇久久久久久久久久久| av国产久精品久网站免费入址| 最近手机中文字幕大全| av线在线观看网站| 国产成人精品在线电影| 久久精品国产自在天天线| 亚洲五月色婷婷综合| 你懂的网址亚洲精品在线观看| 天天操日日干夜夜撸| 国产熟女欧美一区二区| www.色视频.com| 狂野欧美激情性bbbbbb| 亚洲欧美成人精品一区二区| 成人黄色视频免费在线看| 最近手机中文字幕大全| 黄色一级大片看看| 久久久久人妻精品一区果冻| 午夜免费观看性视频| 国产欧美日韩综合在线一区二区| av国产精品久久久久影院| 日韩欧美一区视频在线观看| 免费观看a级毛片全部| 在线观看免费日韩欧美大片| 国产永久视频网站| 最近手机中文字幕大全| 精品酒店卫生间| 成人影院久久| 女人被躁到高潮嗷嗷叫费观| 只有这里有精品99| 免费看av在线观看网站| 中文欧美无线码| 少妇 在线观看| 国产成人免费无遮挡视频| 中国三级夫妇交换| xxxhd国产人妻xxx| 亚洲美女搞黄在线观看| 少妇被粗大猛烈的视频| 如何舔出高潮| 国产精品成人在线| 欧美国产精品va在线观看不卡| 免费高清在线观看日韩| 午夜福利视频精品| av有码第一页| 亚洲色图综合在线观看| 热99久久久久精品小说推荐| 亚洲国产av影院在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美激情极品国产一区二区三区 | a 毛片基地| 母亲3免费完整高清在线观看 | www.熟女人妻精品国产 | 色婷婷av一区二区三区视频| 日韩一本色道免费dvd| 免费观看性生交大片5| 人妻一区二区av| 久久久精品94久久精品| 考比视频在线观看| 国产麻豆69| 亚洲精品国产av蜜桃| 亚洲激情五月婷婷啪啪| av卡一久久| 久久亚洲国产成人精品v| 欧美另类一区| 日韩视频在线欧美| 欧美精品国产亚洲| av在线老鸭窝| 国产黄色免费在线视频| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区 | 日本欧美视频一区| 亚洲美女黄色视频免费看| 久久久久网色| 色网站视频免费| 成人手机av| 久久久久精品性色| 精品酒店卫生间| 久久热在线av| 国产伦理片在线播放av一区| 国产亚洲午夜精品一区二区久久| 伊人久久国产一区二区| 又大又黄又爽视频免费| 少妇人妻精品综合一区二区| 中文字幕人妻丝袜制服| 久久99热这里只频精品6学生| 少妇被粗大的猛进出69影院 | 丰满迷人的少妇在线观看| 日本91视频免费播放| 晚上一个人看的免费电影| a级片在线免费高清观看视频| 一本色道久久久久久精品综合| 免费看av在线观看网站| 精品酒店卫生间| 中文精品一卡2卡3卡4更新| 国产亚洲一区二区精品| 亚洲av国产av综合av卡| 老司机影院成人| 伦理电影大哥的女人| 亚洲色图综合在线观看| 欧美老熟妇乱子伦牲交| 一本色道久久久久久精品综合| 在线观看人妻少妇| √禁漫天堂资源中文www| 免费在线观看完整版高清| 人人澡人人妻人| 不卡视频在线观看欧美| 久久精品久久精品一区二区三区| 国产伦理片在线播放av一区| 国产亚洲最大av| 捣出白浆h1v1| 日韩欧美精品免费久久| 久久精品人人爽人人爽视色| 卡戴珊不雅视频在线播放| 久久久久久久久久人人人人人人| 少妇猛男粗大的猛烈进出视频| 不卡视频在线观看欧美| 精品一区在线观看国产| 深夜精品福利| 搡老乐熟女国产| 精品人妻熟女毛片av久久网站| 午夜福利视频精品| 伦理电影大哥的女人| 美国免费a级毛片| 欧美日韩一区二区视频在线观看视频在线| 午夜老司机福利剧场| 国产成人精品久久久久久| 高清欧美精品videossex| 国产精品嫩草影院av在线观看| 男人舔女人的私密视频| 69精品国产乱码久久久| 在线亚洲精品国产二区图片欧美| 免费观看a级毛片全部| 中文字幕亚洲精品专区| 26uuu在线亚洲综合色| 欧美亚洲日本最大视频资源| 最新中文字幕久久久久| 丁香六月天网| 久久人人爽人人爽人人片va| 成人综合一区亚洲| av免费在线看不卡| 国产免费又黄又爽又色| 国产成人免费观看mmmm| 国产在视频线精品| 免费少妇av软件| 精品午夜福利在线看| 91成人精品电影| 在线观看三级黄色| 在线观看免费视频网站a站| 精品一区在线观看国产| 欧美+日韩+精品| 色视频在线一区二区三区| 久久精品国产亚洲av涩爱| 在线观看免费高清a一片| 国产精品一区二区在线观看99| 亚洲丝袜综合中文字幕| 2021少妇久久久久久久久久久| 久久精品国产亚洲av天美| 制服诱惑二区| 国产精品蜜桃在线观看| 国产高清三级在线| 精品一区在线观看国产| 国产成人精品久久久久久| 免费看不卡的av| 国产熟女午夜一区二区三区| 亚洲国产av新网站| 成人二区视频| 久久人妻熟女aⅴ| 男女国产视频网站| 日韩欧美一区视频在线观看| 色吧在线观看|