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

    梅雨期長江中下游地區(qū)濕位渦的特征及其與青藏高原的聯(lián)系

    2022-12-03 02:40:30由麗蓉姜繼蘭陳笑晨馬婷婷劉屹岷吳國雄
    大氣科學(xué) 2022年6期
    關(guān)鍵詞:長江中下游地區(qū)雨帶負(fù)值

    由麗蓉 2, 姜繼蘭 2 陳笑晨 馬婷婷 劉屹岷 2 吳國雄 2

    1 中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室(LASG), 北京 100029

    2 中國科學(xué)院大學(xué),北京 100049

    3 University of Exeter, United Kingdom

    4 福建省氣候中心,福州 350008

    1 引言

    長江中下游地區(qū)地處東亞季風(fēng)區(qū),面積數(shù)十萬平方公里,是中國最主要的工農(nóng)業(yè)中心,也是中國人口最稠密的地區(qū)之一。每年6月中旬左右,東亞夏季風(fēng)向北推進(jìn),冷暖空氣在此交匯,形成梅雨(Ding, 1992; Si et al., 2009)。作為東亞夏季風(fēng)進(jìn)程中最重要的一個雨季,長江中下游梅雨(以下簡稱“梅雨”)一般持續(xù)一個月左右,降水通常較為集中,且多伴隨災(zāi)害性暴雨(陶詩言, 1980)。梅雨異常所帶來的降水多寡往往是引起長江中下游地區(qū)旱澇異常的主要原因,對這里的人民生活和社會經(jīng)濟(jì)發(fā)展產(chǎn)生重要影響。如第二階段的梅雨,即所謂的“二度梅”,對1998年長江流域特大洪澇災(zāi)害的發(fā)生有重要影響(貝耐芳和趙思雄, 2002),而2020年夏季持續(xù)時間長達(dá)62天的超長梅雨帶來了自1961年以來最大的梅雨降水量(Ding et al.,2021),造成長江流域自1998年以來最為嚴(yán)重的流域性洪澇災(zāi)害。

    前人研究指出,梅雨受到來自熱帶、副熱帶和中高緯多個不同影響因子直接或間接的影響,表現(xiàn)出多時間尺度的特征(Weng et al., 1999; 丁一匯等,2007; Ding et al., 2020)。在 年 代 際 時 間 尺 度 上,梅雨主要受太平洋年代際振蕩(Zhu et al., 2011,2015; Sun et al., 2019)和大西洋多年代際振蕩(Lu et al., 2006)的調(diào)制。另外,ENSO的年代際變化對梅雨也有影響(Wang et al., 2012)。在年際尺度上,梅 雨 主 要 受ENSO(Huang and Wu, 1989;Tanaka, 1997; Wang et al., 2008; Ding et al., 2020)、熱帶印度洋海溫異常(Kosaka et al., 2011; Wang et al., 2013; Yamaura and Tomita, 2014)、北大西洋濤動(Wu et al., 2009; Wang et al., 2018)、歐亞大陸中高緯波列(Liu et al., 2019)、急流(李崇銀等,2004; Sampe and Xie, 2010; Li and Zhang, 2014)等的影響。歷史上,長江流域的兩次特大洪澇災(zāi)害(1998年和2016年)均伴隨著前冬El Ni?o事件的發(fā)展(鮑媛媛, 2021),而2020年長江流域超強(qiáng)梅雨的發(fā)生與印度洋海溫異常增暖有關(guān)(Takaya et al., 2020; Zhou et al., 2021)。在次季節(jié)尺度上,熱帶季節(jié)內(nèi)振蕩(Zhu et al., 2003; 毛江玉和吳國雄,2005; 齊 艷 軍 等, 2016; Hsu et al., 2016; Shao et al.,2018; Zhang W J et al., 2021)和 準(zhǔn) 雙 周 振 蕩(Ouyang and Liu, 2020; Ding et al., 2021; Zhang G S et al., 2021)是影響梅雨的兩個重要因子。另有研究指出,北大西洋濤動位相的次季節(jié)調(diào)整是導(dǎo)致2020年梅雨期我國東部雨帶南北擺動的重要原因(Liu et al., 2020)。

    青藏高原是世界上海拔最高的高原,且地形復(fù)雜,其表面巨大的熱力作用對東亞夏季風(fēng)的影響已經(jīng)達(dá)成 廣 泛共識(Wu and Zhang, 1998; Yasunari and Miwa, 2006; Wang et al., 2018; Duan et al.,2020)。研究表明,春季和夏季高原熱源偏強(qiáng)時,都對應(yīng)著長江中下游地區(qū)夏季降水偏多(羅會邦和陳蓉, 1995; Zhao and Chen, 2001; Duan et al., 2005;Hsu and Liu, 2003; Wang et al., 2009; He C et al.,2019)。Dong et al.(2019)的研究指出青藏高原東南部地表潛熱通量的異常是長江流域梅雨降水異常的關(guān)鍵信號源。另外,高原冬春季積雪對夏季長江流域的降水也有影響,吳統(tǒng)文和錢正安(2000)指出當(dāng)冬春季高原積雪偏多時,夏季長江流域和江南北部降水偏多;當(dāng)高原積雪偏少時,夏季江淮地區(qū)降水偏少。另一方面,高原近地表和邊界層內(nèi)存在強(qiáng)烈的地氣相互作用和復(fù)雜的地形效應(yīng),形成了獨(dú)特的天氣系統(tǒng)(如高原低渦和青藏高原切變線等)。這些天氣系統(tǒng)的發(fā)展和東移往往造成梅雨區(qū)災(zāi)害性天氣過程(Tao and Ding, 1981; Wang and Orlanski, 1987; Zheng et al., 2013; Ni et al., 2017;馬婷等, 2020)。如1998年夏季長江流域出現(xiàn)了罕見的持續(xù)性暴雨和洪澇災(zāi)害過程,而大部分暴雨的發(fā)生是由青藏高原發(fā)展東移的低渦或切變線造成的(師春香等, 2000; 張順利和陶詩言, 2002)。另外,青藏高原上空對流云團(tuán)和中尺度對流系統(tǒng)十分活躍(江吉喜等, 1996; Li et al., 2008),一些對流云團(tuán)在合適的環(huán)流條件下發(fā)展東移也會導(dǎo)致長江中下游地區(qū)大暴雨的發(fā)生(卓嘎等, 2002; 傅慎明等,2011)。最新的研究表明青藏高原可以通過影響位渦的制造在高原近地層形成低壓系統(tǒng),這些系統(tǒng)在東移的過程中形成自高原向長江中下游地區(qū)伸展的高位渦(或正渦度)帶,對2020年梅雨降水異常有重要貢獻(xiàn)(Li et al., 2021; Ma et al., 2022)。然而,目前關(guān)于高原通過影響位渦進(jìn)而影響長江中下游地區(qū)梅雨的研究還多以個例分析為主。

    位渦理論是大氣動力學(xué)的核心內(nèi)容之一,被廣泛應(yīng)用于天氣分析和研究工作中。一般而言,位渦對許多中高緯度天氣系統(tǒng)的發(fā)展和移動具有指示意義,但在對流層中低層,尤其是在低緯度地區(qū),位渦變得很弱。同時,位渦不包含水汽的特征,在考慮降水特別是梅雨降水等濕空氣的運(yùn)動時具有局限性,由此,濕位渦(MPV)應(yīng)運(yùn)而生。吳國雄等(1995)從原始方程出發(fā),在導(dǎo)出精確形式濕位渦方程的基礎(chǔ)上,證明了絕熱、無摩擦的飽和濕空氣中濕位渦守恒。在這種守恒性質(zhì)的制約下,渦旋垂直渦度的發(fā)展由大氣穩(wěn)定度、斜壓性和風(fēng)的垂直切變等因素決定,由此提出了傾斜渦度發(fā)展(SVD)理論。濕位渦守恒性的提出和SVD理論的建立對暴雨過程的診斷提供了理論依據(jù),大量研究結(jié)果表明,無論是濕位渦(高守亭等, 2002; 劉峰等, 2011;肖云清等, 2020)和還是SVD理論(吳國雄和蔡雅萍, 1997; 吳國雄和劉還珠, 1999; 王子謙等, 2010;李靜楠等, 2016)都是暴雨診斷中的有力工具,對強(qiáng)降水的落區(qū)具有較好的指示意義。如高守亭等(2002)從理論分析和資料診斷兩個方面證明了濕位渦異常區(qū)與強(qiáng)降水區(qū)有很好的對應(yīng)關(guān)系。目前,濕位渦和SVD理論多用于個例和局地暴雨的分析,較少被用于氣候態(tài)梅雨的分析中。根據(jù)SVD理論,由于梅雨鋒附近濕等熵面傾斜很大,十分有利于垂直渦度增長從而導(dǎo)致暴雨的發(fā)生。因此,從大尺度環(huán)流的角度去討論梅雨期濕位渦的分布特征是十分有意義的。

    綜上所述,前人對于引起梅雨異常的因子已有很多研究,但對于梅雨期長江中下游地區(qū)濕位渦的分布特征以及青藏高原在其中所起的作用還少有研究。因此,本文將圍繞如下問題開展研究:(1)長江中下游地區(qū)梅雨期濕位渦的分布特征如何?與梅雨降水的關(guān)系如何?(2)青藏高原對長江中下游地區(qū)濕位渦的分布有何影響?本文其余部分安排如下:第二節(jié)介紹本文所用的資料、方法和模式;第三節(jié)主要討論氣候態(tài)梅雨期濕位渦的分布特征,發(fā)現(xiàn)濕位渦斜壓項(xiàng)MPV2與梅雨雨帶有較好的對應(yīng)關(guān)系,并簡要分析其中的原因;第四節(jié)利用FGOALSf3-L海氣耦合模式設(shè)計數(shù)值試驗(yàn),分析青藏高原對長江中下游地區(qū)MPV2分布的影響;第五節(jié)對全文進(jìn)行總結(jié)。

    2 資料、方法和模式介紹

    2.1 資料

    本文所使用的資料包括:(1)歐洲中期天氣預(yù)報中心提供的ERA-Interim 1980~2017年逐日再分析數(shù)據(jù)集(Dee et al., 2011),變量包括溫度、比濕、經(jīng)向風(fēng)分量、緯向風(fēng)分量、垂直速度、位勢高度,水平空間分辨率為1.5°×1.5°,垂直方向有37層(https://apps.ecmwf.int/datasets/data/interimfull-daily/levtype=pl/ [2022-06-02]);(2)中 國 氣象局提供的中國地面氣候資料日值數(shù)據(jù)集的逐日降水資料;(3)日本氣象廳提供的JRA-55再分析數(shù)據(jù)集(Kobayashi et al., 2015; Harada et al., 2016)中的全球降水資料(https://rda.ucar.edu/datasets/ds628.0 [2022-03-21])。

    2.2 梅雨期的定義及分析

    本文采用丁一匯等(2007)使用的長江流域梅雨季定義標(biāo)準(zhǔn):(1)選擇長江中下游地區(qū)5個具有代表性的氣象臺站(武漢、九江、蕪湖、南京和上海)來代表梅雨期的降水量。首先,如果該日滿足至少有2個臺站發(fā)生降水,且5個臺站總的日降水量大于或等于10 mm,則被定義為雨日。(2)在梅雨季,西太平洋副熱帶高壓脊線的位置在20°~25°N。(3)入梅日需滿足以下條件:從該日起之后的10個連續(xù)的雨日中,若條件(1)中的雨日等于或大于一半,則開始降水的第一天定義為入梅日。(4)在梅雨季期間,有一次或多次降水事件發(fā)生,每次事件必須持續(xù)6天以上的雨日,且5個測站的日平均降水量大于25 mm。(5)最后一次降雨事件中的第2天定義為出梅日。由此得到的1980~2017年入、出梅日期見表1。

    根據(jù)表1可以計算出1980~2017年間平均入梅日期為6月18日,平均出梅日期為7月12日。本文將每年入梅日期作為第0天,在其之前(后)的天數(shù)記為負(fù)(正)天數(shù)進(jìn)行合成分析,由此得到各變量在入梅日前后的特征。

    表1 1980~2017年入梅、出梅日期Table 1 Dates of the onset and retreat of the Meiyu season from 1980 to 2017

    2.3 濕位渦和SVD理論

    濕位渦是一個包含大氣動力、熱力以及水汽特性的綜合變量。本文采用吳國雄等(1995)提出的濕位渦和SVD理論進(jìn)行診斷分析。在靜力近似條件下,且假定垂直速度的水平變化比水平速度的垂直切變小得多時,濕位渦的計算公式如下:

    其中,MPV是濕位渦,單位為PVU(1 PVU=10?6m2K s?1kg?1);g是重力加速度,g=9.8 m s?2;θe為相當(dāng)位溫。MPV又可分解為兩項(xiàng)之和:MPV1和MPV2,其中MPV1是濕位渦的垂直分量積;MPV2是水平分量積,公式如下:

    (2)式中 ζap=f+(?v/?x??u/?y)p=f+ζp,為絕對渦度的垂直分量??梢姡琈PV1為絕對渦度的垂直分量與相當(dāng)位溫垂直梯度的乘積,稱作濕正壓項(xiàng)。在北半球,由于絕對渦度的垂直分量f+ζp在大多數(shù)情況下為正值,所以 ?θe/?p的符號決定了MPV1的正負(fù):當(dāng) (?θe/?p)<0時,大氣是對流穩(wěn)定狀態(tài),此時MPV1>0;當(dāng) (?θe/?p)>0時,大氣是對流不穩(wěn)定狀態(tài),此時MPV1<0。MPV2是水平風(fēng)的垂直切變與相當(dāng)位溫水平梯度的乘積,稱作濕斜壓項(xiàng)。根據(jù)SVD理論,在濕位渦守恒的前提下,不管大氣的對流穩(wěn)定性或是濕對稱穩(wěn)定性如何,由于濕等熵面的傾斜,水平風(fēng)的垂直切變增大或者濕斜壓性的增加能夠?qū)е麓怪睖u度的顯著發(fā)展,且濕等熵面的傾斜越大,垂直渦度的發(fā)展越激烈。傾斜渦度發(fā)展的必要條件是M PV2/(?θe/?p)>0。

    2.4 數(shù)值試驗(yàn)介紹

    本文采用由中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點(diǎn)實(shí)驗(yàn)室發(fā)展的新一代氣候系統(tǒng)模式FGOALS-f3-L(The lowresolution version of Flexible Global Ocean-Atmosphere-Land System model, finite-volume version 3)進(jìn)行數(shù)值模擬研究,該版本的模式參加了最新的第六次國際耦合模式比較計劃(CMIP6)。評估顯示,該版本的模式在模擬能力上較之前版本有了明顯提升,可以更好地模擬出歷史氣候中各種要素的特征,包括長期趨勢和氣候平均態(tài)的分布,有助于研究氣候變化機(jī)理、提高預(yù)測氣候的能力(He B et al., 2019; Guo et al., 2020; Tang et al., 2020)。

    本文一共設(shè)計了兩組試驗(yàn),控制試驗(yàn)是CMIP6的piControl(pre-industrial control simulation)試驗(yàn),即將外強(qiáng)迫(溫室氣體、太陽輻射、氣溶膠、土地利用等)維持在1850年的水平下驅(qū)動全球耦合模式進(jìn)行500年以上的長期積分;敏感性試驗(yàn)與控制試驗(yàn)相同,但將青藏—伊朗高原(TIP)500 m以上的地形削除,以考察高原大地形的影響。

    3 入梅前后濕位渦與降水的聯(lián)系

    3.1 濕位渦在入梅前后的特征

    濕位渦正壓項(xiàng)MPV1和斜壓項(xiàng)MPV2有不同的物理意義,二者的配置對降水的發(fā)生發(fā)展有重要影響。濕正壓項(xiàng)MPV1表示慣性穩(wěn)定性和對流穩(wěn)定性的作用(劉峰等, 2011),在對流層中上層MPV1一般為正,在對流層低層如果存在對流不穩(wěn)定((?θe/?p)>0),則MPV1為負(fù)。從850 hPa MPV1的水平分布(圖1a,b和c)可以看出,入梅前后MPV1在中國東部和西太平洋地區(qū)基本為負(fù)值分布。入梅前第3天MPV1負(fù)值中心出現(xiàn)在梅雨區(qū)北部的江淮地區(qū),量值小于?0.3 PVU(圖1a)。根據(jù)趙亮和丁一匯(2008),表明有弱冷空氣在此聚集,此時有較強(qiáng)的對流不穩(wěn)定,為入梅提供了條件。高原北部和西北太平洋地區(qū)有較強(qiáng)的正MPV1分布。此時東亞降水中心在日本南部,沿著降水帶也有較弱的正MPV1分布,中國的降水大值中心位于華南。入梅當(dāng)天隨著中國東部降水中心北移至長江中下游地區(qū),之前位于江淮地區(qū)的MPV1負(fù)值中心減弱北移,長江中下游地區(qū)為較弱的MPV1負(fù)值區(qū)(圖1b)。高原北部和西北太平洋地區(qū)仍然維持較強(qiáng)的MPV1正值中心。入梅以后,歐亞大陸和西北太平洋地區(qū)的正負(fù)MPV1的分布與入梅當(dāng)天基本一致,負(fù)值中心位于梅雨雨帶北側(cè)(圖1c)。

    圖1 1980~2017年850 hPa (a–c)濕位渦正壓項(xiàng)(MPV1,填色,單位:PVU)和(d–f)濕位渦斜壓項(xiàng)(MPV2,填色,單位:PVU)及降水量(等值線,單位:mm d?1)在(a、d)入梅前第三天、(b、e)入梅當(dāng)天和(c、f)入梅后第3天的分布。紅色方框表示長江中下游梅雨區(qū)(28°~34°N,112°~123°E),下同F(xiàn)ig. 1 Distributions of (a–c) MPV1 (shaded, units: PVU) and (d–f) MPV2 (shaded, units: PVU) at 850 hPa and precipitation (contours, units: mm d?1)on (a, d) the 3rd day before Meiyu onset, (b, e) the Meiyu onset day, and (c, f) the 3rd day after Meiyu onset during the period of 1980–2017. The red box indicates the Meiyu area in the middle and lower Yangtze River basin (28°–34°N, 112°–123°E), the same below

    850 hPa MPV2(圖1d,e和f)的分布表明,入梅前后MPV2負(fù)值帶均與東亞主要降水中心重合:入梅前3天,雨帶和MPV2負(fù)值中心主要分布在日本南部和我國華南地區(qū)(圖1d);入梅當(dāng)天,雨帶和MPV2負(fù)值帶北跳到長江中下游地區(qū)(圖1e),且在整個梅雨期內(nèi)一直維持在這里。這一特征從中國東部(112°~123°E)MPV2和降水隨時間的演變(圖2)中可以更加直觀地看出,在入梅前5天左右,中國東部主要的降水中心和MPV2負(fù)值中心位于北緯25°~30°N之間的長江以南地區(qū),此時,MPV2負(fù)值中心較降水中心略微偏北。入梅后,隨著長江中下游地區(qū)(28°~34°N)降水中心的出現(xiàn)和增強(qiáng),該地區(qū)負(fù)MPV2的強(qiáng)度也增強(qiáng)。此時,MPV2負(fù)值中心與降水中心更為接近。之后二者的變化基本保持一致,在28°~34°N之間南北擺動,直至出梅以后,二者均在長江中下游地區(qū)顯著減弱。根據(jù)SVD理論(吳國雄等, 1995),入梅之后,長江中下游地區(qū)MPV2負(fù)值帶的形成暗示著這里斜壓性增強(qiáng),有利于垂直渦度發(fā)展增強(qiáng)從而造成梅雨降水。另外值得注意的是,在入梅前后,青藏高原南部南亞一帶也有較強(qiáng)的降水,但MPV2在該區(qū)域很小,且沒有明顯的負(fù)值中心,表明高原南部降水和梅雨降水的性質(zhì)可能存在不同。

    圖2 1980~2017年梅雨期前后112°~123°E平均的850 hPa MPV2(填色,單位:PVU)和降水量(等值線,單位:mm d?1)的時間—緯度剖面。黑色虛線表示長江中下游梅雨區(qū)域(28°~34°N),紅色虛線表示出梅日Fig. 2 Time–latitude cross-section of MPV2 at 850 hPa (shaded, units: PVU) and precipitation (contours, units: mm d?1), averaged over 112°–123°E around the Meiyu season during the period of 1980–2017. The black dotted lines indicate the Meiyu area in the middle and lower Yangtze River basin(28°–34°N), and the red line indicates the date of the Meiyu retreat

    500 hPa MPV1的分布在一定程度上可以表示天氣形勢。入梅前后,除青藏高原南部外,東亞和西北太平洋地區(qū)基本為正MPV1分布(圖3a,b和c),表示對流層中層主要是以對流穩(wěn)定的層結(jié)為主。MPV1大值中心在入梅前后均位于東亞降水中心的北側(cè)。入梅前3天,江淮地區(qū)向東至日本一帶為MPV1大值帶,強(qiáng)度維持在0.3 PVU以上,表明這些地區(qū)有弱冷空氣入侵。入梅當(dāng)天,MPV1的分布與入梅前3天基本一致,只是MPV1大值帶在中國東部略有加強(qiáng)且向北移動,而長江中下游梅雨區(qū)正MPV1略有減弱,說明入梅當(dāng)天梅雨區(qū)冷空氣活動減弱、而梅雨區(qū)以北冷空氣活動增強(qiáng)(圖3b)。入梅后第3天,梅雨區(qū)正MPV1的分布與入梅當(dāng)天基本一致,只是梅雨區(qū)以北MPV1正值中心強(qiáng)度有所減弱(圖3c)。這與入梅之后,梅雨區(qū)盛行的南風(fēng)上升氣流造成負(fù)的MPV1經(jīng)向和垂直方向平流、以及來自中高緯的冷空氣入侵減弱有關(guān)(后文將詳細(xì)分析)。

    圖3 同圖1,但為500 hPa(a–c)MPV1(填色,單位:PVU)和(d–f)MPV2(填色,單位:PVU)及降水量(等值線,單位:mm d?1)的分布Fig. 3 As in Fig. 1, but for distributions of (a–c) MPV1 (shaded, units: PVU) and (d–f) MPV2 (shaded, units: PVU) at 500 hPa and precipitation(contours, units: mm d?1)

    與MPV1相反,除青藏高原南部外,MPV2在500 hPa基本為負(fù)值分布(圖3d,e和f)。負(fù)MPV2的大值中心呈帶狀分布在高原北部和東亞降水中心北側(cè)。高原北部的MPV2負(fù)值帶在入梅前后沒有明顯的位置和強(qiáng)度的變化,而東亞降水中心北側(cè)的MPV2負(fù)值帶隨雨帶的北抬而北抬。與MPV1相似,梅雨區(qū)MPV2負(fù)值帶的強(qiáng)度在入梅之后也有所減弱。

    以上分析可知,入梅前后,長江中下游地區(qū)MPV1與MPV2大值帶均與梅雨雨帶有較好的對應(yīng),且隨雨帶的北抬而北抬。在對流層中層二者均位于雨帶北側(cè),而對流層低層,MPV2與雨帶的對應(yīng)更好,幾乎與雨帶重合。

    濕位渦的垂直結(jié)構(gòu)對降水的發(fā)生發(fā)展也有較好的指示作用。圖4給出了入梅前后MPV1和相當(dāng)位溫在中國東部(112°~123°E)的垂直剖面??梢钥闯鋈朊非昂螅瑵竦褥孛嬖陂L江中下游地區(qū)傾斜最大,而在其南北兩側(cè)基本為水平分布。高層濕等熵面的分布相對較為密集,對應(yīng)大的?θe/?p,導(dǎo)致高層MPV1較大。這主要是由于高層水汽含量較少,?θe/?p的分布主要受靜力穩(wěn)定度(?θ/?p)的影響,氣壓的變化在高層相對較小,導(dǎo)致高層靜力穩(wěn)定度較大,因此,?θe/?p在高層較大。在30°N以南的對流層高層,MPV1大于2 PVU的層次較薄,而30°N以北,該層次顯著變厚,長江中下游地區(qū)高層為劇烈傾斜分布的MPV1。

    圖4 1980~2017年入梅前第6天至入梅后第3天中國東部(112°~123°E)MPV1(填色,單位:PVU)、相當(dāng)位溫(灰色等值線,單位:K)、經(jīng)向和垂直方向的MPV1平流(綠色等值線,單位:10?5 PVU s?1)和風(fēng)場 [矢量,v單位:m s?1,ω單位:(?1/50)Pa s?1] 的垂直剖面(紅色豎虛線表示梅雨區(qū)緯度范圍28°~34°N,下同)Fig. 4 Vertical cross section of MPV1 (shaded, units: PVU), equivalent potential temperature (gray contours, units: K), meridional and vertical MPV1 advection (green contours, units: 10 ?5 PVU s?1), and wind [vector, v units: m s?1, ω units: (?1/50) Pa s?1] over eastern China (112°–123°E) from the 6th day before the Meiyu onset to the 3rd day after the Meiyu onset during the period of 1980–2017. The red dotted lines indicate the Meiyu area in the middle and lower Yangtze River basin (28°–34°N), the same below

    入梅前,長江中下游梅雨區(qū)以南盛行偏南風(fēng)上升氣流,而梅雨區(qū)主要以偏北風(fēng)下沉氣流為主,冷暖空氣交匯在梅雨區(qū)南部。在偏北風(fēng)下沉氣流較強(qiáng)的梅雨區(qū)北側(cè)MPV1經(jīng)向和垂直方向平流為正,有利于來自高緯度地區(qū)對流層中上層的高M(jìn)PV1沿濕等熵面向梅雨區(qū)滲透(圖4a和b),造成冷空氣在此聚集以及MPV1在梅雨區(qū)北部的大值中心,與前文分析相對應(yīng)。入梅前第6天,梅雨區(qū)0.2 PVU等值線下傳至700 hPa,而700 hPa以下是微弱的負(fù)MPV1分布(圖4a)。入梅前第3天,梅雨區(qū)北部對流層中層的偏北風(fēng)下沉氣流略有減弱,其所對應(yīng)的經(jīng)向和垂直方向的正MPV1平流也相應(yīng)減弱(圖4b),但正的MPV1經(jīng)向和垂直方向平流仍能導(dǎo)致梅雨區(qū)北側(cè)MPV1的增強(qiáng),這也解釋了圖3所示的梅雨區(qū)北側(cè)MPV1在入梅之前有所增強(qiáng)的原因。入梅當(dāng)天,偏南風(fēng)上升氣流占據(jù)整個梅雨區(qū),冷暖空氣交匯于梅雨區(qū)北側(cè)。由于MPV1在35°N以南呈現(xiàn)出由低緯向高緯、由低層向高層逐漸增大的特征,入梅后梅雨區(qū)偏南風(fēng)上升氣流的發(fā)展導(dǎo)致整個梅雨區(qū)以負(fù)MPV1經(jīng)向和垂直方向平流為主,使得梅雨區(qū)MPV1較入梅前有所減弱,對流層中層MPV1大值區(qū)出現(xiàn)在梅雨區(qū)北側(cè)(圖4c),與前文對應(yīng)。另外,入梅后梅雨區(qū)以北上升氣流也有所發(fā)展,冷空氣活動減弱,導(dǎo)致該處MPV1經(jīng)向和垂直方向的平流由入梅前的正平流轉(zhuǎn)為入梅后的負(fù)平流,使得入梅后梅雨區(qū)以北MPV1的強(qiáng)度也有所減弱。之后,偏南風(fēng)上升氣流一直維持在梅雨區(qū)內(nèi),直至梅雨結(jié)束。

    入梅前后,濕斜壓項(xiàng)MPV2在整個對流層基本為負(fù)值分布(圖5),其分布結(jié)構(gòu)與MPV1存在顯著差別。梅雨區(qū)高層200 hPa以上有一明顯的負(fù)值中心,該中心在入梅之后略有北移。200 hPa以下,MPV2負(fù)值帶沿梅雨區(qū)及其北側(cè)傾斜的濕等熵面密集帶分布,隨高度向北傾斜。入梅前第6天,該MPV2負(fù)值帶的極大值中心位于300 hPa梅雨區(qū)北界,中心值低于?0.24 PVU,而低層850 hPa附近MPV2負(fù)極大值位于梅雨區(qū)南界,中心值低于?0.08 PVU(圖5a)。入梅前第3天,MPV2負(fù)值帶在200 hPa以下均有所減弱,且整體略微向北移動,梅 雨 區(qū)700 hPa以 下 為MPV2界 于?0.04~?0.08 PVU的分布(圖5b)。入梅當(dāng)天,MPV2負(fù)值帶進(jìn)一步北移,300 hPa的負(fù)極大值中心北移至40°N附近,與對流層中層MPV1大值帶在入梅當(dāng)天北移至梅雨區(qū)以北相對應(yīng)。對流層中層MPV2負(fù)值帶強(qiáng)度較之前變化不大,但梅雨區(qū)700 hPa以下負(fù)MPV2顯著增強(qiáng)至低于?0.08 PVU(圖5c),表明入梅后梅雨區(qū)斜壓性增強(qiáng)。從濕等熵面的分布可以推斷,梅雨區(qū)800 hPa以上存在(?θe/?p)<0,因此,有MPV2/(?θe/?p)>0,滿足SVD的必要條件。根據(jù)SVD理論,入梅當(dāng)天梅雨區(qū)低層斜壓性增強(qiáng),有利于暖濕氣流沿濕等熵面上滑,從而導(dǎo)致暖濕氣流的垂直渦度顯著發(fā)展(Cui et al., 2003),造成梅雨降水。入梅后第3天,MPV2負(fù)值帶在200 hPa以下基本維持不再北移,強(qiáng)度較入梅當(dāng)天有所減弱(圖5d)。MPV2在梅雨區(qū)南部的減弱也與偏南風(fēng)上升氣流引起的MPV2經(jīng)向和垂直方向的平流有關(guān)。

    圖5 同圖4,但為MPV2(填色,單位:PVU)和相當(dāng)位溫(等值線,單位:K)的分布Fig. 5 As in Fig. 4, but for the distributions of MPV2 (shaded, units: PVU) and equivalent potential temperature (contours, units: K)

    3.2 濕位渦斜壓項(xiàng)MPV2與梅雨雨帶對應(yīng)的原因分析

    前文分析可知,濕位渦斜壓項(xiàng)MPV2負(fù)值帶與梅雨雨帶有較好的對應(yīng),在850 hPa基本與梅雨雨帶重合,在500 hPa位于梅雨雨帶北側(cè)。由公式(3)可知,MPV2主要由水平風(fēng)的垂直切變和相當(dāng)位溫的水平梯度構(gòu)成,下面我們將分別討論這兩者的相對貢獻(xiàn)。由于入梅前后相當(dāng)位溫在梅雨區(qū)大致為平行于緯圈的帶狀分布(圖6),尤其是入梅后,這種帶狀分布更為平直,使得梅雨區(qū)相當(dāng)位溫的水平梯度主要以經(jīng)向梯度(圖6a,b和c)為主,緯向梯度(圖6d,e和f)較經(jīng)向梯度小一個量級左右。所以公式(3)等號右邊第二項(xiàng)(g(?v/?p)(?θe/?x)) 量級相對于第一項(xiàng)( ?g(?u/?p)(?θe/?y))較小,可以忽略不計,即梅雨期MPV2主要是由相當(dāng)位溫經(jīng)向梯度( ?θe/?y)和緯向風(fēng)垂直切變(?u/?p) 的乘積決定(M PV2≈?g(?u/?p)(?θe/?y))。下面分別就二者對MPV2負(fù)值帶的相對貢獻(xiàn)進(jìn)行定性分析。

    圖6 1980~2017年500 hPa相當(dāng)位溫(等值線,單位:K)及其(a–c)經(jīng)向和(d–f)緯向梯度(填色,單位:10?5 K m?1)在(a,d)入梅前第3天、(b,e)入梅當(dāng)天和(c,f)入梅后第3天的分布Fig. 6 Distributions of equivalent potential temperature (contours, units: K) and its (a–c) meridional and (d–f) zonal gradient (shaded, units: 10?5 K m?1) at 500 hPa on (a, d) the 3rd day before the Meiyu onset, (b, e) the Meiyu onset day, and (c, f) the 3rd day after the Meiyu onset during the period of 1980–2017

    圖7分別給出了入梅前后500 hPa和850 hPa相當(dāng)位溫經(jīng)向梯度? θe/?y的分布。從圖中可以看出,除低緯地區(qū)外,對流層中低層 ?θe/?y的分布均與MPV2相似。具體而言,高原上空500 hPa MPV2北負(fù)南正的偶極型分布主要由? θe/?y所決定(圖7a,b和c),而 ?u/?p在高原上均為負(fù)值分布,且負(fù)值中心更偏東(圖8a,b和c)。高原上空 ?θe/?y的分布主要是由于夏季高原是熱源,其南北兩側(cè)的溫度梯度剛好相反;同時由于高原上水汽含量相對較少,相當(dāng)位溫的分布與位溫的分布較為相似,因此?θe/?y在高原上呈現(xiàn)出北負(fù)南正的偶極型分布。對于長江中下游地區(qū),對流層中層? θe/?y和 ?u/?p的分布共同決定了MPV2在雨帶北側(cè)的大值帶。而對流層低層,相對于 ?u/?p位于雨帶偏南側(cè)的負(fù)中心(圖8d,e和f),? θe/?y在 雨帶北側(cè)的負(fù)中心(圖7d,e和f)更好地對應(yīng)了MPV2負(fù)值帶。由此可見,無論是高原地區(qū)還是長江中下游地區(qū), ?θe/?y的分布與MPV2更加類似,其對MPV2負(fù)值帶的貢獻(xiàn)要更大一些。

    圖7 1980~2017年(a–c)500 hPa和(d–f)850 hPa ? θe/?y(填色,單位:10?5 K m?1)和降水量(等值線,單位:mm d?1)在(a,d)入梅前第3天、(b,e)入梅當(dāng)天和(c,f)入梅后第3天的分布Fig. 7 Distributions of ? θe/?y (shaded, units: 10?5 K m?1 at (a–c) 500 hPa and (d–f) 850 hPa and precipitation (contours, units: mm d?1) on (a, d) the 3rd day before the Meiyu onset, (b, e) the Meiyu onset day, and (c, f) the 3rd day after the Meiyu onset during the period of 1980–2017

    圖8 同圖7,但為? u/?p(填色,單位:10?5 m s?1 Pa?1)和降水量(等值線,單位:mm d?1)的分布Fig. 8 As in Fig. 7, but for ? u/?p (shaded, units: 10?5 m s?1 Pa?1) and precipitation (contours, units: mm d?1)

    下面用SVD理論來定性地解釋MPV2負(fù)值帶與雨帶的對應(yīng)關(guān)系。由SVD的必要條件(MPV2/(?θe/?p)>0) 可知,對流不穩(wěn)定( (?θe/?p)>0)的情況發(fā)生時,若MPV2>0,傾斜渦度可能增長;對流穩(wěn)定( (?θe/?p)<0)的情況下,若MPV2<0,傾斜渦度也可能增長。

    由前文分析可知,梅雨發(fā)生前后MPV2負(fù)值帶沿傾斜的濕等熵面密集帶分布,且在整個梅雨期基本維持。梅雨區(qū)850 hPa以上MPV2負(fù)值帶滿足(?θe/?p)<0,且濕等熵面在此處傾斜最為劇烈(圖5),所以該區(qū)域滿足SVD的必要條件。入梅后,梅雨區(qū)低空斜壓性增強(qiáng),根據(jù)SVD理論,當(dāng)南風(fēng)氣流沿850 hPa以上傾斜的濕等熵面上滑(圖4c和d)時垂直渦度可以得到顯著發(fā)展加強(qiáng),有利于梅雨降水的產(chǎn)生。

    4 濕位渦的分布與青藏高原的可能聯(lián)系

    前文指出,在500 hPa,濕位渦斜壓項(xiàng)MPV2在梅雨期的分布特征為一條西起青藏高原北部、向東延伸至西北太平洋的狹長負(fù)值帶,大值中心位于青藏高原北部。在入梅之前(圖3d),該負(fù)值帶在中國東部很強(qiáng),位置偏南,位于30°N附近,中心強(qiáng)于?0.14 PVU,達(dá)到與MPV1相同的量級。入梅以后,該負(fù)值帶在中國東部有所減弱并北移,但高原北部的負(fù)值帶仍穩(wěn)定維持(圖3e和f)。那么,青藏高原是否影響了MPV2負(fù)值帶的分布呢?

    為了回答這一問題,本節(jié)利用FGOALS-f3-L海氣耦合模式設(shè)計了去掉青藏—伊朗高原地形的模擬試驗(yàn)??刂圃囼?yàn)的結(jié)果表明FGOALS-f3-L模式基本模擬出了夏季MPV2負(fù)值帶自高原北部向東延伸至西北太平洋地區(qū)的分布,且負(fù)值帶的極大值中心分布在高原北部,只是強(qiáng)度略弱于觀測(圖9a)。在無高原試驗(yàn)中,該負(fù)值帶在高原北部至西北太平洋地區(qū)顯著減弱甚至消失,同時中國東部至日本一帶的降水也隨之減弱(圖9b)。兩個試驗(yàn)的差值場進(jìn)一步突出了高原的存在對MPV2負(fù)值帶分布的影響(圖9c),說明青藏高原大地形的存在使高原北部向東至西北太平洋地區(qū)形成了一條強(qiáng)MPV2負(fù)值帶,影響中國東部的降水。圖中還可以看出青藏高原大地形的存在為6月中南半島和南海一帶帶來更多降水。

    圖9 FGOALS-f3-L模式輸出的6月500 hPa MPV2(填色,單位:PVU)和降水量(等值線,單位:mm d?1)在(a)控制試驗(yàn)和(b)無高原試驗(yàn)中的分布,(c)為二者的差異[(a)-(b)]Fig. 9 Distributions of the June MPV2 (shaded, units: PVU) at 500 hPa and precipitation (contours, units: mm d?1) in (a) the control run and (b) the no Tibetan–Iranian Plateau run obtained from FGOALS-f3-L. (c) is the difference between (a) and (b) [(a)-(b)]

    對比兩組試驗(yàn)中國東部MPV2的垂直分布(圖10)可以看出,在高原的作用下,6月梅雨區(qū)附近出現(xiàn)和觀測(圖5)一致的由低層向高層向北傾斜的MPV2負(fù)值帶(圖10a和c),中心位于梅雨區(qū)北側(cè)300~400 hPa附近,略微低于觀測的300 hPa(圖5)。由前文分析可知,這種垂直結(jié)構(gòu)對于傾斜渦度的發(fā)展十分重要。而在沒有高原的情況下(圖10b),MPV2的傾斜結(jié)構(gòu)更為偏北并且強(qiáng)度很弱,對長江中下游地區(qū)的影響較小。此外,青藏高原的存在也對200 hPa以上MPV2的分布有很大影響,但其中的聯(lián)系需要更深入的研究與討論。

    圖10 同圖9,但為112°~123°E平均的MPV2的垂直分布Fig. 10 As in Fig. 9, but for the vertical cross section of MPV2 averaged over 112°–123°E

    前文分析可知,相當(dāng)位溫的經(jīng)向梯度(? θe/?y)和緯向風(fēng)的垂直切變( ?u/?p)共同決定了青藏高原和中國東部MPV2負(fù)值帶的形成,那么,青藏高原是通過影響何種物理過程進(jìn)而影響MPV2的分布呢?為了回答這一問題,我們對比了兩個試驗(yàn)中?θe/?y和 ?u/?p的分布(圖11)。從圖中可以看出,控制試驗(yàn)也較好地模擬出了 ?θe/?y(圖11a)和?u/?p(圖11d)的大值帶分布,與觀測十分相似。而無高原試驗(yàn)中,由于高原熱源作用的消失,?θe/?y在高原上北負(fù)南正的偶極型分布消失,取而代之的是高原和中國東部較弱的 ?θe/?y分布(圖11b)。在水汽含量較少的對流層中上層,θe的 分布與 θ的分布較為相似。根據(jù)熱成風(fēng)關(guān)系,溫度的經(jīng)向梯度與緯向風(fēng)的垂直切變成正比,因此,無高原試驗(yàn)中?u/?p在高原和中國東部也顯著減弱(圖11e)。由此可見,高原的存在能夠同時影響到? θe/?y和 ?u/?p的分布進(jìn)而影響到MPV2的分布。

    圖11 FGOALS-f3-L模式輸出的500 hPa ?θe/?y (填色,單位:10?5 K m?1)在(a)控制試驗(yàn)和(b)無高原試驗(yàn)中的分布,(c)為二者的差異 [(a)-(b)]。(d)、(e)和(f)同(a)、(b)和(c),但為? u/?p(填色,單位:10?5 m s?1 Pa?1)的分布Fig. 11 Distributions of ? θe/?y (shaded, units: 10?5 K m?1) at 500 hPa in (a) the control run and (b) the no Tibetan–Iranian Plateau run obtained from FGOALS-f3-L. (c) is the difference between (a) and (b) [(a)-(b)]. (d)–(f) are identical to (a)–(c) but for ? u/?p (shaded, units: 10?5 m s?1 Pa?1)

    5 結(jié)論

    梅雨異常所帶來的旱澇災(zāi)害對東亞地區(qū)經(jīng)濟(jì)有很大影響,因此,研究梅雨所伴隨的大尺度環(huán)流及其異常機(jī)制具有重大意義。濕位渦分析是暴雨研究中常用的方法,將其運(yùn)用到氣候態(tài)梅雨的研究可以從不同的角度對梅雨發(fā)生發(fā)展的機(jī)理做出解釋。本文分析了梅雨期濕位渦在長江中下游地區(qū)的分布特征以及青藏高原對氣候態(tài)梅雨期濕位渦分布的影響,主要得到以下結(jié)論:

    (1)入梅前后,MPV1與MPV2大值帶均沿傾斜的濕等熵面分布,且與梅雨雨帶有較好的對應(yīng),隨雨帶的北抬而北抬。在對流層中層二者均位于雨帶北側(cè),而對流層低層,MPV2與梅雨雨帶近乎重合。

    (2)MPV2負(fù)值帶和梅雨雨帶對應(yīng)較好的原因是:梅雨發(fā)生前后MPV2負(fù)值帶沿傾斜的濕等熵面密集帶分布,滿足SVD的必要條件,當(dāng)暖濕空氣沿傾斜的濕等熵面上滑時,其垂直渦度可以得到顯著發(fā)展增強(qiáng),有利于降水的產(chǎn)生。

    (3)梅雨期MPV2負(fù)值帶西起青藏高原向東延伸至西北太平洋。FGOALS-f3-L無高原試驗(yàn)的結(jié)果表明,青藏高原大地形條件對MPV2負(fù)值帶的形成有重要影響,在無青藏—伊朗高原的試驗(yàn)中,長江中下游地區(qū)MPV2負(fù)值帶顯著減弱甚至消失。

    本研究發(fā)現(xiàn)除了暴雨分析之外,濕位渦運(yùn)用于氣候態(tài)梅雨的分析也有顯著的意義,揭示了青藏高原可以通過影響濕位渦斜壓項(xiàng)MPV2的分布進(jìn)而影響到梅雨。但是部分研究內(nèi)容還需要進(jìn)一步完善,如本文分析了梅雨期濕位渦分布的特征發(fā)現(xiàn)MPV2負(fù)值帶的形成可能受到青藏高原的影響,但沒有深入討論該影響的機(jī)理,在下一步研究中可以就該問題進(jìn)行探究。

    猜你喜歡
    長江中下游地區(qū)雨帶負(fù)值
    王金成:秋風(fēng)
    石油過剩:一桶油如何突然跌至負(fù)值
    英語文摘(2020年7期)2020-09-21 03:40:56
    2017年長江中下游地區(qū)一次持續(xù)性異常降水過程分析
    回味暑假生活,看看動物小伙伴們的表現(xiàn)
    長江中下游地區(qū)梅雨期異常年降水及大氣熱源分布特征
    青藏高原春季土壤濕度對長江中下游地區(qū)初夏短期氣候影響的數(shù)值模擬
    論青藏高原溫度對長江中下游地區(qū)降水的影響
    各種各樣的雨
    最大熵算法在氣象雨量預(yù)測中的應(yīng)用分析
    小型蒸發(fā)器蒸發(fā)量測定為負(fù)值的原因及改進(jìn)措施
    18禁黄网站禁片午夜丰满| 一级片免费观看大全| 19禁男女啪啪无遮挡网站| 成人免费观看视频高清| 久久久久九九精品影院| 女人被狂操c到高潮| 怎么达到女性高潮| 国产精品1区2区在线观看.| 91成人精品电影| 深夜精品福利| 日韩高清综合在线| 伦理电影免费视频| 亚洲专区国产一区二区| 美女午夜性视频免费| 一个人免费在线观看的高清视频| 女警被强在线播放| 91精品国产国语对白视频| 亚洲成人免费电影在线观看| 免费av毛片视频| 亚洲av第一区精品v没综合| 一本久久中文字幕| 亚洲精品中文字幕一二三四区| 国产99白浆流出| 国产精品野战在线观看| 夜夜夜夜夜久久久久| 麻豆国产av国片精品| 亚洲精品久久国产高清桃花| 欧美在线一区亚洲| 国产成人精品久久二区二区91| 精品国产乱码久久久久久男人| 91精品三级在线观看| 亚洲男人天堂网一区| 国产欧美日韩一区二区三| 亚洲激情在线av| 乱人伦中国视频| 一二三四社区在线视频社区8| 一个人观看的视频www高清免费观看 | videosex国产| 可以在线观看的亚洲视频| xxx96com| 亚洲一区二区三区色噜噜| 日韩欧美一区二区三区在线观看| 大型av网站在线播放| 女人精品久久久久毛片| 久久精品亚洲精品国产色婷小说| 亚洲美女黄片视频| 狠狠狠狠99中文字幕| 首页视频小说图片口味搜索| 亚洲自偷自拍图片 自拍| 人人妻人人澡欧美一区二区 | 国产成人免费无遮挡视频| 999久久久精品免费观看国产| 国产精品久久久人人做人人爽| 男女下面插进去视频免费观看| 国产区一区二久久| 亚洲一区二区三区不卡视频| 国产xxxxx性猛交| 亚洲九九香蕉| 91麻豆av在线| 手机成人av网站| 波多野结衣av一区二区av| 久久久久久大精品| 日本vs欧美在线观看视频| 成人欧美大片| 亚洲av成人不卡在线观看播放网| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 午夜日韩欧美国产| 亚洲自拍偷在线| 国产精品自产拍在线观看55亚洲| 很黄的视频免费| 看片在线看免费视频| 少妇粗大呻吟视频| 国产精品日韩av在线免费观看 | 精品久久久久久,| 成人亚洲精品一区在线观看| 热99re8久久精品国产| 啦啦啦韩国在线观看视频| 激情视频va一区二区三区| 久久天堂一区二区三区四区| 成人手机av| xxx96com| 成人国产一区最新在线观看| 黄频高清免费视频| 法律面前人人平等表现在哪些方面| 亚洲国产精品合色在线| 国产精品久久久久久精品电影 | 国产精品一区二区在线不卡| 欧美亚洲日本最大视频资源| 国产精品日韩av在线免费观看 | 国产又色又爽无遮挡免费看| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 999久久久精品免费观看国产| 日本精品一区二区三区蜜桃| 可以在线观看的亚洲视频| 久久 成人 亚洲| 纯流量卡能插随身wifi吗| 免费高清视频大片| 老司机午夜十八禁免费视频| 欧美精品啪啪一区二区三区| 两个人免费观看高清视频| 最近最新免费中文字幕在线| 欧美久久黑人一区二区| 亚洲国产高清在线一区二区三 | 国产精品精品国产色婷婷| 精品一区二区三区四区五区乱码| 免费在线观看视频国产中文字幕亚洲| 9191精品国产免费久久| 嫩草影视91久久| 亚洲狠狠婷婷综合久久图片| 中出人妻视频一区二区| 国产亚洲精品久久久久久毛片| 成人三级黄色视频| 欧美激情极品国产一区二区三区| 97超级碰碰碰精品色视频在线观看| 亚洲国产中文字幕在线视频| 一级黄色大片毛片| 国产一区二区在线av高清观看| 女警被强在线播放| 女性被躁到高潮视频| 亚洲国产看品久久| 桃红色精品国产亚洲av| 久久青草综合色| 国产精品综合久久久久久久免费 | 成人18禁高潮啪啪吃奶动态图| 亚洲熟女毛片儿| 国产成年人精品一区二区| 国产亚洲av嫩草精品影院| 亚洲 欧美一区二区三区| 久久久久久人人人人人| 精品免费久久久久久久清纯| 熟女少妇亚洲综合色aaa.| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久精品电影 | 精品一区二区三区av网在线观看| 99国产极品粉嫩在线观看| 俄罗斯特黄特色一大片| 91麻豆av在线| 欧美日韩黄片免| 精品国产美女av久久久久小说| 国产成人精品在线电影| 在线观看免费视频网站a站| 国产一区在线观看成人免费| 亚洲精品中文字幕在线视频| 亚洲 欧美 日韩 在线 免费| 亚洲专区国产一区二区| 97碰自拍视频| 精品久久久久久,| 宅男免费午夜| 国产一区二区三区在线臀色熟女| 久久中文字幕一级| 欧美成人一区二区免费高清观看 | 无人区码免费观看不卡| 老司机午夜福利在线观看视频| 久久久久精品国产欧美久久久| 亚洲av五月六月丁香网| 国产一卡二卡三卡精品| e午夜精品久久久久久久| 搞女人的毛片| 欧美午夜高清在线| 国产国语露脸激情在线看| 咕卡用的链子| 久久精品国产亚洲av香蕉五月| 一边摸一边抽搐一进一小说| 精品国产乱子伦一区二区三区| 欧美亚洲日本最大视频资源| 日本免费a在线| 欧美中文综合在线视频| 亚洲中文字幕日韩| 亚洲人成77777在线视频| 狂野欧美激情性xxxx| 国产精品国产高清国产av| 国产蜜桃级精品一区二区三区| 青草久久国产| 免费在线观看视频国产中文字幕亚洲| 欧美黄色片欧美黄色片| 欧美丝袜亚洲另类 | 色精品久久人妻99蜜桃| 在线国产一区二区在线| 色综合亚洲欧美另类图片| 欧美乱妇无乱码| 热re99久久国产66热| 亚洲国产看品久久| 国产国语露脸激情在线看| 99久久99久久久精品蜜桃| avwww免费| 1024视频免费在线观看| 精品久久久久久久人妻蜜臀av | 午夜精品久久久久久毛片777| 精品卡一卡二卡四卡免费| 在线观看一区二区三区| 给我免费播放毛片高清在线观看| 搞女人的毛片| 黄色视频不卡| 多毛熟女@视频| 色哟哟哟哟哟哟| 日韩视频一区二区在线观看| 欧美大码av| 男人舔女人的私密视频| 精品免费久久久久久久清纯| 免费观看精品视频网站| 午夜福利欧美成人| 一进一出抽搐动态| 岛国在线观看网站| 99在线人妻在线中文字幕| 欧美丝袜亚洲另类 | 午夜两性在线视频| 国产精品亚洲美女久久久| 美女国产高潮福利片在线看| 国产精品亚洲美女久久久| 国产精品亚洲美女久久久| 一级a爱片免费观看的视频| 免费看十八禁软件| 亚洲国产中文字幕在线视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲伊人色综图| 一进一出好大好爽视频| 久久久久国内视频| 欧美大码av| 夜夜看夜夜爽夜夜摸| 日本五十路高清| 亚洲精品久久成人aⅴ小说| 午夜免费观看网址| 日本免费a在线| aaaaa片日本免费| 亚洲第一电影网av| 亚洲最大成人中文| 欧美丝袜亚洲另类 | 欧美老熟妇乱子伦牲交| АⅤ资源中文在线天堂| 免费高清在线观看日韩| 亚洲电影在线观看av| 国产精品综合久久久久久久免费 | 男女下面进入的视频免费午夜 | 亚洲狠狠婷婷综合久久图片| 性欧美人与动物交配| 丁香六月欧美| 亚洲av成人不卡在线观看播放网| 久久中文字幕一级| 欧美中文综合在线视频| 欧美国产精品va在线观看不卡| 国产精品98久久久久久宅男小说| 美女免费视频网站| 90打野战视频偷拍视频| 成人国语在线视频| av电影中文网址| 婷婷精品国产亚洲av在线| 天天添夜夜摸| 成年人黄色毛片网站| 成人18禁高潮啪啪吃奶动态图| 女同久久另类99精品国产91| 午夜免费观看网址| 18禁裸乳无遮挡免费网站照片 | 国产成人一区二区三区免费视频网站| 啦啦啦韩国在线观看视频| 久久人人97超碰香蕉20202| 视频区欧美日本亚洲| 中亚洲国语对白在线视频| 久久香蕉激情| 久久久久精品国产欧美久久久| 久久久久国内视频| 国产麻豆成人av免费视频| 真人做人爱边吃奶动态| 久久久久久人人人人人| 精品国产超薄肉色丝袜足j| 久久精品影院6| 久久精品91无色码中文字幕| 两性夫妻黄色片| 夜夜躁狠狠躁天天躁| 99国产精品99久久久久| 咕卡用的链子| 色哟哟哟哟哟哟| 啦啦啦免费观看视频1| 黄片小视频在线播放| 国产av精品麻豆| 首页视频小说图片口味搜索| 9色porny在线观看| 国产精品乱码一区二三区的特点 | 久久久久亚洲av毛片大全| 久久婷婷成人综合色麻豆| 黄色女人牲交| 大陆偷拍与自拍| 国产色视频综合| 在线观看免费日韩欧美大片| 99在线视频只有这里精品首页| 正在播放国产对白刺激| 看免费av毛片| 又黄又粗又硬又大视频| 女人高潮潮喷娇喘18禁视频| 欧美日本亚洲视频在线播放| 亚洲国产欧美日韩在线播放| 欧美av亚洲av综合av国产av| 国产一区二区三区视频了| 黑人巨大精品欧美一区二区蜜桃| 久久狼人影院| 91精品三级在线观看| 国产av又大| 大码成人一级视频| 国产成人免费无遮挡视频| 手机成人av网站| 国产一级毛片七仙女欲春2 | 久久青草综合色| 黑人欧美特级aaaaaa片| 韩国av一区二区三区四区| 国产成人免费无遮挡视频| 久久精品国产综合久久久| 在线观看www视频免费| 免费少妇av软件| 宅男免费午夜| 一二三四社区在线视频社区8| 久久精品成人免费网站| 母亲3免费完整高清在线观看| 精品一区二区三区四区五区乱码| 很黄的视频免费| 久久影院123| 老司机深夜福利视频在线观看| 黄色 视频免费看| 国产一卡二卡三卡精品| 90打野战视频偷拍视频| 中国美女看黄片| 亚洲成av人片免费观看| 国产1区2区3区精品| 国产熟女xx| 欧美精品亚洲一区二区| 丝袜在线中文字幕| 99国产精品免费福利视频| av网站免费在线观看视频| 又大又爽又粗| 免费在线观看亚洲国产| 极品人妻少妇av视频| 看黄色毛片网站| 国语自产精品视频在线第100页| 91字幕亚洲| 精品久久久久久久毛片微露脸| 高清黄色对白视频在线免费看| 国产成人一区二区三区免费视频网站| 大码成人一级视频| 国产aⅴ精品一区二区三区波| 男男h啪啪无遮挡| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品亚洲av一区麻豆| 色av中文字幕| 亚洲国产中文字幕在线视频| 国产熟女xx| 精品一品国产午夜福利视频| 国产精品久久久久久人妻精品电影| 欧美一级a爱片免费观看看 | 亚洲熟妇熟女久久| 亚洲七黄色美女视频| 精品国产乱子伦一区二区三区| 纯流量卡能插随身wifi吗| АⅤ资源中文在线天堂| 国产亚洲欧美精品永久| 黄网站色视频无遮挡免费观看| 久久国产乱子伦精品免费另类| 日韩视频一区二区在线观看| 午夜久久久在线观看| 国产精品影院久久| 久久精品国产亚洲av高清一级| 亚洲一区高清亚洲精品| 在线观看免费视频日本深夜| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全免费视频| 日本撒尿小便嘘嘘汇集6| 最新美女视频免费是黄的| 九色亚洲精品在线播放| 91麻豆精品激情在线观看国产| av超薄肉色丝袜交足视频| 午夜福利成人在线免费观看| 国产精品电影一区二区三区| 亚洲成人精品中文字幕电影| 国产精品av久久久久免费| 97超级碰碰碰精品色视频在线观看| 免费无遮挡裸体视频| 搡老岳熟女国产| 亚洲精品中文字幕一二三四区| 久久草成人影院| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 黄片播放在线免费| www.自偷自拍.com| 国产精品秋霞免费鲁丝片| 一边摸一边做爽爽视频免费| 国产成人精品久久二区二区91| netflix在线观看网站| 美女午夜性视频免费| 日本vs欧美在线观看视频| 女同久久另类99精品国产91| 一级片免费观看大全| 变态另类成人亚洲欧美熟女 | 麻豆一二三区av精品| 日本 av在线| 桃红色精品国产亚洲av| 女同久久另类99精品国产91| 男人的好看免费观看在线视频 | 国产精品一区二区在线不卡| 欧美黑人精品巨大| 十分钟在线观看高清视频www| 黄色成人免费大全| 男人的好看免费观看在线视频 | 999久久久精品免费观看国产| 国产亚洲av嫩草精品影院| 国产成+人综合+亚洲专区| 国产精品av久久久久免费| 国产成人av激情在线播放| 日韩三级视频一区二区三区| 色尼玛亚洲综合影院| 国产精品久久久人人做人人爽| 男人操女人黄网站| 九色亚洲精品在线播放| 国产熟女xx| 久久精品亚洲精品国产色婷小说| 九色亚洲精品在线播放| 亚洲av美国av| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女 | 日日摸夜夜添夜夜添小说| 少妇熟女aⅴ在线视频| 精品人妻1区二区| АⅤ资源中文在线天堂| 久久久久国内视频| 亚洲av片天天在线观看| 亚洲欧美日韩无卡精品| 一区二区三区精品91| 免费在线观看视频国产中文字幕亚洲| 亚洲国产精品999在线| 99国产极品粉嫩在线观看| 亚洲第一av免费看| 黑丝袜美女国产一区| 一进一出抽搐动态| 精品电影一区二区在线| 亚洲aⅴ乱码一区二区在线播放 | 男女下面进入的视频免费午夜 | 每晚都被弄得嗷嗷叫到高潮| 亚洲九九香蕉| 精品国产乱子伦一区二区三区| 色婷婷久久久亚洲欧美| 丰满的人妻完整版| 精品国产一区二区久久| 成人免费观看视频高清| 久久久久国内视频| 脱女人内裤的视频| 亚洲第一av免费看| 亚洲人成伊人成综合网2020| netflix在线观看网站| 美女大奶头视频| 一级,二级,三级黄色视频| 午夜免费鲁丝| 黄色女人牲交| 69av精品久久久久久| 很黄的视频免费| 黑人巨大精品欧美一区二区mp4| 午夜a级毛片| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| 一夜夜www| 好看av亚洲va欧美ⅴa在| 亚洲欧美一区二区三区黑人| 精品国产国语对白av| 丝袜美腿诱惑在线| 亚洲avbb在线观看| 黄色成人免费大全| 久久久精品欧美日韩精品| 午夜精品久久久久久毛片777| 国产成+人综合+亚洲专区| 国产成人精品在线电影| 成熟少妇高潮喷水视频| 久久久久久免费高清国产稀缺| 国产精品九九99| 狂野欧美激情性xxxx| 神马国产精品三级电影在线观看 | 久久久久九九精品影院| 日韩国内少妇激情av| 91麻豆精品激情在线观看国产| 成人三级做爰电影| 国产一区二区三区视频了| 性色av乱码一区二区三区2| 午夜日韩欧美国产| 久久精品国产清高在天天线| 色综合欧美亚洲国产小说| 午夜两性在线视频| 少妇裸体淫交视频免费看高清 | 久久久国产精品麻豆| 99久久99久久久精品蜜桃| 99热只有精品国产| 色播在线永久视频| 色播亚洲综合网| 国产精品99久久99久久久不卡| 91国产中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 久久香蕉激情| 午夜精品国产一区二区电影| 欧美成人一区二区免费高清观看 | 日韩欧美免费精品| 国产成人精品久久二区二区免费| 国产成人精品无人区| 99国产精品一区二区三区| 午夜免费激情av| 夜夜看夜夜爽夜夜摸| 很黄的视频免费| av天堂在线播放| 欧美另类亚洲清纯唯美| 亚洲美女黄片视频| 亚洲欧洲精品一区二区精品久久久| 国产一区二区三区综合在线观看| 国产精品野战在线观看| cao死你这个sao货| 亚洲最大成人中文| 国产熟女午夜一区二区三区| 亚洲三区欧美一区| 最新在线观看一区二区三区| 麻豆久久精品国产亚洲av| 亚洲一卡2卡3卡4卡5卡精品中文| 搡老岳熟女国产| 人妻久久中文字幕网| 免费高清视频大片| 91老司机精品| 国内毛片毛片毛片毛片毛片| 欧美激情久久久久久爽电影 | 欧美日韩一级在线毛片| 999久久久精品免费观看国产| 免费搜索国产男女视频| 亚洲最大成人中文| 亚洲无线在线观看| 麻豆一二三区av精品| 日本精品一区二区三区蜜桃| 国产人伦9x9x在线观看| 波多野结衣av一区二区av| 99精品久久久久人妻精品| 日本免费一区二区三区高清不卡 | 免费无遮挡裸体视频| 亚洲人成77777在线视频| 91在线观看av| 亚洲熟妇熟女久久| 99精品在免费线老司机午夜| 欧美黑人精品巨大| 色尼玛亚洲综合影院| 亚洲美女黄片视频| 天天添夜夜摸| 欧美 亚洲 国产 日韩一| 亚洲av电影在线进入| 在线国产一区二区在线| 国产欧美日韩精品亚洲av| 亚洲 欧美一区二区三区| 欧美黑人精品巨大| 国产精品永久免费网站| 亚洲成人精品中文字幕电影| 亚洲免费av在线视频| 国产午夜福利久久久久久| 免费高清在线观看日韩| 中文字幕最新亚洲高清| 国产午夜福利久久久久久| 99香蕉大伊视频| 亚洲av五月六月丁香网| 麻豆av在线久日| 久久精品亚洲熟妇少妇任你| 欧美激情高清一区二区三区| 精品卡一卡二卡四卡免费| 黄色视频,在线免费观看| 看片在线看免费视频| 久久久国产精品麻豆| 国产激情欧美一区二区| 天天添夜夜摸| 9色porny在线观看| 国产xxxxx性猛交| 国内精品久久久久久久电影| 制服人妻中文乱码| 欧美在线黄色| 正在播放国产对白刺激| 亚洲久久久国产精品| 国产成人av激情在线播放| 亚洲精品粉嫩美女一区| 男男h啪啪无遮挡| 十八禁网站免费在线| 看片在线看免费视频| 欧美老熟妇乱子伦牲交| 国内精品久久久久精免费| av在线播放免费不卡| 亚洲国产中文字幕在线视频| 午夜福利欧美成人| svipshipincom国产片| 国产在线观看jvid| 多毛熟女@视频| 久热这里只有精品99| 国产人伦9x9x在线观看| 男人舔女人下体高潮全视频| xxx96com| 99热只有精品国产| 可以在线观看毛片的网站| 亚洲久久久国产精品| 一级毛片高清免费大全| 长腿黑丝高跟| 精品卡一卡二卡四卡免费| 国产黄a三级三级三级人| 亚洲国产欧美网| videosex国产| 欧美绝顶高潮抽搐喷水| 日韩欧美免费精品| 午夜激情av网站| 欧美成人性av电影在线观看| av免费在线观看网站| 母亲3免费完整高清在线观看| 亚洲熟女毛片儿| a在线观看视频网站| 一区二区三区高清视频在线| 亚洲熟妇中文字幕五十中出| 脱女人内裤的视频| 国产欧美日韩精品亚洲av| 看黄色毛片网站| 国产成人精品在线电影| 久久精品亚洲熟妇少妇任你| 两性夫妻黄色片| 色婷婷久久久亚洲欧美|