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

    2021年漢江秋汛的水雨情及環(huán)流異常特征研究*

    2023-08-09 06:23:34吳翠紅劉佳明
    氣象 2023年7期
    關(guān)鍵詞:丹江口水庫洪峰漢江

    高 琦 吳翠紅 徐 明 劉佳明 彭 濤,4

    1 武漢中心氣象臺(tái),武漢 430074 2 中國氣象局水文氣象重點(diǎn)開放實(shí)驗(yàn)室,北京 100081 3 中國氣象局武漢暴雨研究所暴雨監(jiān)測(cè)預(yù)警湖北省重點(diǎn)實(shí)驗(yàn)室,武漢 430205 4 三峽國家氣候觀象臺(tái),湖北宜昌 443099 5 長江勘測(cè)規(guī)劃設(shè)計(jì)研究有限責(zé)任公司,武漢 430010

    提 要:利用NCEP/NCAR再分析資料及常規(guī)氣象水文實(shí)況觀測(cè)資料,研究2021年漢江秋汛的水雨情及環(huán)流異常特征,并初步分析致洪成因,結(jié)果表明:2021年漢江流域降水呈“西多東少”分布,漢江上游累計(jì)降水量為400~985 mm,累計(jì)面雨量居1961年以來第一位。秋汛期間丹江口水庫出現(xiàn)10 000 m3·s-1以上的洪峰過程7次,漢江中下游主要水文控制站點(diǎn)累計(jì)超警時(shí)間為151~580 h,最大超警幅度為0.29~1.56 m,呈現(xiàn)越到下游超警時(shí)間越長、超警幅度越大的特征。與2011年和2017年水情對(duì)比表明,2021年漢江上游洪峰出現(xiàn)時(shí)間最早,洪峰次數(shù)最多,上游高水位流量持續(xù)時(shí)間最長,漢江中下游累計(jì)超警時(shí)間也最長,水情異常程度明顯更強(qiáng)。在大氣環(huán)流形勢(shì)方面,巴爾喀什湖—貝加爾湖地區(qū)為負(fù)異常,我國大陸東部至西北太平洋為正異常,日本島以東地區(qū)為負(fù)異常,副熱帶高壓、南亞高壓、東亞副熱帶西風(fēng)急流較歷史同期偏強(qiáng),異常的低層流場(chǎng)及水汽輻合中心的形成,以及對(duì)流層高層的輻散抽吸,是漢江上游降水得以持續(xù)并階段性增強(qiáng)的有利環(huán)流配置。當(dāng)?shù)そ谒畮烊霂炝髁拷咏?000 m3·s-1后,漢江上游超過60 mm的5 d累計(jì)面雨量與10 000 m3·s-1以上洪峰的出現(xiàn)密切相關(guān)。上游來水增大后,漢江中下游河道逐級(jí)收窄是導(dǎo)致此地降水并不明顯卻出現(xiàn)超警洪水的主要原因。

    引 言

    漢江為長江中游最大支流,其地勢(shì)西高東低(圖1),以湖北省丹江口水庫以上為上游,河谷狹窄,長約925 km;丹江口至鐘祥為中游,河谷較寬,沙灘多,長約270 km;鐘祥至漢口為下游,長約382 km,流經(jīng)江漢平原。漢江中下游河道呈“漏斗型”,河道逐級(jí)收窄,其中皇莊至澤口段河道寬度為600~2500 m,澤口至仙桃段河道寬度為300~600 m,至仙桃以下河道寬度僅為200~300 m,因而過洪能力有限。漢江上游大部地區(qū)處于華西秋雨區(qū)內(nèi),某些年份由于降水持續(xù)時(shí)間長、累計(jì)雨量大,易造成秋汛。漢江上游來水增大后,由于河道斷面越到下游越窄,又受長江水位頂托的影響,往往導(dǎo)致泄洪能力不足,造成漢江中下游也泛濫成災(zāi),歷史上為了確保漢江中下游城市安全,曾21次啟用分蓄洪區(qū)。

    圖1 漢江流域分區(qū)及氣象、水文站點(diǎn)分布Fig.1 Distribution of sub-basin and meteorological and hydrological stations in the Hanjiang River Basin

    漢江秋汛一般發(fā)生于華西秋雨偏強(qiáng)的年份(殷淑燕和黃春長,2012;肖鶯等,2013)。一系列研究表明,華西秋雨的強(qiáng)弱多寡與東亞大氣環(huán)流異常密切相關(guān),西太平洋副熱帶高壓(以下簡稱副高)、南亞高壓和東亞副熱帶西風(fēng)急流等季風(fēng)環(huán)流成員對(duì)華西秋雨的維持及增強(qiáng)具有重要作用。鮑媛媛等(2003)、李傳浩等(2015)的研究均表明,巴爾喀什湖地區(qū)500 hPa準(zhǔn)穩(wěn)定的低壓槽及距平場(chǎng)上的負(fù)異常,以及偏西偏北的副高,是導(dǎo)致華西地區(qū)持續(xù)陰雨的典型環(huán)流特征。白虎志和董文杰(2004)指出,副高、印緬槽、貝加爾湖低槽是華西秋雨的主要影響系統(tǒng),當(dāng)貝加爾湖、印緬槽深且副高強(qiáng)時(shí),有利于華西多秋雨。俞亞勛等(2013)、吳瑤和杜良敏(2016)的研究表明,副高的西伸滯留及其脊線的準(zhǔn)雙周振蕩對(duì)華西秋雨具有重要作用,當(dāng)脊線準(zhǔn)雙周周期較強(qiáng)時(shí),漢江上游9—10月容易降水偏多。羅霄等(2013)總結(jié)了近50 年華西秋雨的典型環(huán)流后發(fā)現(xiàn),東亞副熱帶西風(fēng)急流對(duì)華西地區(qū)的垂直運(yùn)動(dòng)有較大影響,當(dāng)東亞副熱帶西風(fēng)急流中心西移,中心強(qiáng)風(fēng)速帶偏窄時(shí),華西北部地區(qū)上升運(yùn)動(dòng)加強(qiáng),華西秋雨出現(xiàn)北多南少的異常型。進(jìn)一步研究表明,青藏高原地表熱狀況(陳忠明等,2001;鮑媛媛等,2020)、海溫(韓晉平等,2013;鄭然等,2018;彭韻萌和徐涂明,2021)、印度洋偶極子(劉佳等,2015)、熱帶大氣季節(jié)內(nèi)振蕩(趙佳玉等,2016)等是季風(fēng)系統(tǒng)及大氣環(huán)流異常形成的深層機(jī)制。此外,華西秋雨年際及年代際差異明顯,存在4~8 a(王春學(xué)等,2015),5~8 a(蔣竹將等,2014),3、13、17 a(馮麗文和郭其蘊(yùn),1983)以及10~15 a(梁健洪,1989)的周期。由于這種周期性的存在,使得漢江秋汛發(fā)生頻次相對(duì)較少,因此對(duì)于秋汛期降水特征、洪水演變及可能造成的影響,預(yù)報(bào)員往往認(rèn)識(shí)相對(duì)有限。隨著氣候變暖背景加強(qiáng),進(jìn)入21世紀(jì)后,華西秋雨強(qiáng)度明顯增強(qiáng)(羅霄等,2013;俞亞勛等,2013),這使得進(jìn)一步深入認(rèn)識(shí)和了解漢江秋汛的特點(diǎn)和成因迫在眉睫。

    2021年秋季漢江上游出現(xiàn)了超過20年一遇的大洪水,丹江口水庫發(fā)生7次超過10 000 m3·s-1的入庫洪水過程,8月下旬至10月上旬,丹江口水庫秋汛累計(jì)來水量約為340億m3,較常年同期偏多約4倍,為1969年建庫以來歷史同期第一位。通過研究2021年漢江秋汛期間降水的時(shí)空分布及歷史極端性特征,漢江流域水位、流量的演變特征,并與典型秋汛年(2011年和2017年)的水情進(jìn)行對(duì)比,分析致洪降水的大尺度環(huán)流及水汽輸送異常,從氣象、水文等多角度揭示了2021年漢江秋汛的致洪成因,有助于加深對(duì)于漢江秋汛發(fā)生發(fā)展規(guī)律的認(rèn)識(shí),促進(jìn)流域上下游、左右岸聯(lián)動(dòng)防洪工作的開展。

    1 資料和方法

    根據(jù)長江水利委員會(huì)水文局提供的漢江流域水系劃分,將漢江流域細(xì)分為石泉以上、石泉安康、安康丹江口、丹皇區(qū)間和皇莊以下等 5 個(gè)子流域(圖1),該區(qū)域內(nèi)包含氣象站63個(gè)。文中所指的漢江上游為石泉以上、石泉安康、安康丹江口3個(gè)子流域組成的區(qū)域。文中涉及的漢江流域水文站有6個(gè),其中丹江口水庫代表漢江上游來水情況,皇莊、沙洋、岳口、仙桃、漢川為漢江中下游干流主要水文控制站。

    采用的資料包括:NCEP/NCAR再分析資料,1日4次(02、08、14、20時(shí),北京時(shí),下同),分辨率為2.5°×2.5°;漢江流域1961—2021年國家氣象站逐日降水資料;漢江流域主要水文控制站的水位、流量資料。文中所指的歷史平均為1991—2020年的30年平均。面雨量(畢寶貴等,2003;徐晶和姚學(xué)祥,2007;高琦等,2014)的計(jì)算采用站點(diǎn)平均法。

    2 致洪降水的時(shí)空分布特征

    根據(jù)國家氣候中心監(jiān)測(cè),2021年華西秋雨北區(qū)雨期自8月23日開始,至10月16日結(jié)束。漢江流域上游所有站點(diǎn)均處于華西秋雨監(jiān)測(cè)區(qū)內(nèi),其中大部分站點(diǎn)屬華西秋雨監(jiān)測(cè)北區(qū),少數(shù)站點(diǎn)屬監(jiān)測(cè)南區(qū)。若以整個(gè)漢江上游為研究對(duì)象,參照華西秋雨行業(yè)標(biāo)準(zhǔn)(中國氣象局,2019),則自8月21日起,即已滿足50%以上站點(diǎn)降水量≥0.1 mm的秋雨日判定標(biāo)準(zhǔn),且10月8日后不再滿足,考慮該時(shí)間與漢江流域汛情的發(fā)展有較好的對(duì)應(yīng)關(guān)系,因此本文選取8月21日至10月7月為降水研究時(shí)段。

    2.1 致洪降水的空間分布特征

    從2021年8月21日至10月7日漢江流域累計(jì)降水量空間分布可見(圖2a),漢江流域降水呈現(xiàn)西多東少分布特征,漢江下游累計(jì)降水量僅為50~150 mm,中游為150~400 mm,上游則為400~800 mm,局部站點(diǎn)可達(dá)900 mm以上,最強(qiáng)降水中心位于石泉以上子流域的鎮(zhèn)巴站,達(dá)985 mm。分析累計(jì)降水量距平分布(圖2b),除漢江中下游部分地區(qū)接近歷史同期或略偏少外,漢江中游北部及漢江上游的大部地區(qū)較歷史同期偏多1~2倍,局地偏多2倍以上,由此可見2021年秋汛期的降水主要來自于漢江上游地區(qū),因此下文將重點(diǎn)針對(duì)漢江上游,研究其降水的時(shí)間演變及極端性特征。

    圖2 2021年8月21日至10月7日漢江流域(a)總降水量和(b)降水距平百分率分布Fig.2 (a) Accumulated precipitation and (b) precipitation anomaly percentage in the Hanjiang River Basin from 21 August to 7 October 2021

    2.2 致洪降水的時(shí)間演變特征

    從漢江上游出現(xiàn)0.1 mm以上降水站數(shù)占比的逐日演變可見(圖3,虛線),2021年8月21日至10月7日,漢江上游降水過程頻繁,雨日較多,定義單日50%以上站點(diǎn)出現(xiàn)0.1 mm以上降水為一個(gè)雨日,則期間雨日約占總?cè)諗?shù)的60.4%。從逐日面雨量演變可以看出(圖3,柱狀),大致可分為4個(gè)集中期,分別在8月21日至9月6日、9月15—19日、9月22—28日、10月3—6日,其中8月21日至9月6日為降水持續(xù)時(shí)間最長、過程最為頻繁和集中的時(shí)段。從5 d累計(jì)面雨量的演變(圖3,實(shí)線)可以看出,漢江上游5 d累計(jì)面雨量≥60 mm的峰值有7次,其中≥90 mm峰值有3次,分別為8月28日至9月1日(98.7 mm)、9月15—19日(98.6 mm)、9月23—27日(94.6 mm),對(duì)比丹江口水庫的入庫流量發(fā)現(xiàn),在入庫流量接近2000 m3·s-1后,60 mm以上的5 d累計(jì)面雨量與10 000 m3·s-1以上的洪峰有較好的對(duì)應(yīng)關(guān)系,而90 mm以上的5 d累計(jì)面雨量與20 000 m3·s-1以上的洪峰有較好的對(duì)應(yīng)關(guān)系。

    圖3 2021年8月21日至10月7日漢江上游日面雨量、5 d累計(jì)面雨量、≥0.1 mm站點(diǎn)百分率的逐日演變Fig.3 Daily variation of areal rainfall, 5 d accumulated areal rainfall and percentage of site rainfall ≥0.1 mm from 21 August to 7 October 2021

    2.3 致洪降水的極端性特征

    高琦等(2020)基于近60 年降水?dāng)?shù)據(jù),利用相對(duì)閾值法(百分位法)統(tǒng)計(jì)了漢江子流域極端面雨量閾值及5 d累計(jì)面雨量極值,2021年8月21日至10月7日漢江上游3個(gè)子流域(石泉以上、石泉安康、安康丹江口)的逐日面雨量超過極端面雨量95%分位閾值的日數(shù)分別占總?cè)諗?shù)的29.1%、22.9%、16.7%,表明秋汛期間漢江上游子流域約有1/5~1/3的日數(shù)有強(qiáng)降水發(fā)生。9月1—5日石泉以上的5 d累計(jì)面雨量達(dá)117.2 mm, 8月28日至9月1日安康丹江口的5 d累計(jì)面雨量達(dá)115.5 mm,均突破了兩子流域自1961年以來5 d累計(jì)面雨量的極值。

    從1961年以來漢江5個(gè)子流域8月21日至10月7日累計(jì)面雨量的時(shí)間演變(圖4)可以看出,2021年漢江上游的石泉以上、石泉安康、安康丹江口累計(jì)面雨量分別為641.8、641.5、438.4 mm,分別位居1961年以來的第一、第一、第二位,且由石泉以上、石泉安康、安康丹江口3個(gè)子流域組成的漢江上游流域累計(jì)面雨量位居1961年以來第一位。

    綜上所述,2021年漢江秋汛致洪降水主要有以下特征:一是降水開始時(shí)間早、持續(xù)時(shí)間長,從8月21日開始發(fā)展,至10月上旬減弱;二是降水的累計(jì)強(qiáng)度大,漢江上游累計(jì)降水量普遍在400~800 mm,其中局部站點(diǎn)達(dá)900 mm以上;三是降水歷史極端性強(qiáng),漢江上游流域累計(jì)面雨量居1961年以來第一位。

    3 漢江秋汛的水情演變特征

    3.1 2021年漢江流域水情特征分析

    丹江口水庫的入庫流量可反映漢江上游來水的總體情況,從其8月初至10月上旬的演變(圖5a)可以看出,丹江口的入庫流量大致可分為4個(gè)階段,8月21日起,隨著降水進(jìn)入第一個(gè)集中期,丹江口的入庫流量也出現(xiàn)顯著上漲,并在短短的17 d內(nèi)先后在8月23日(14 400 m3·s-1)、8月30日(23 400 m3·s-1)、9月2日(16 400 m3·s-1)、9月6日(18 800 m3·s-1)出現(xiàn)4次10 000 m3·s-1以上的洪峰過程;9月15日起,隨著第2~4個(gè)降水集中期發(fā)展,漢江上游來水再度出現(xiàn)3輪明顯上漲,丹江口水庫分別在9月19日(22 800 m3·s-1)、9月29日(24 900 m3·s-1)、10月7日(10 500 m3·s-1)出現(xiàn)3次洪峰過程。結(jié)合丹江口水庫的出庫流量可以看出,丹江口水庫對(duì)7次洪峰的削峰率分別為46.5%、67.0%、47.0%、46.3%、71.0%、55.8%、2.0%,有效減輕了上游洪水對(duì)漢江中下游的影響,并在10月10日14時(shí)首次實(shí)現(xiàn)170 m的蓄水目標(biāo)。

    注:藍(lán)色實(shí)線:水位,紫色實(shí)線:入庫流量及流量,黃色實(shí)線:出庫流量,橙色實(shí)線:警戒水位。圖5 2021年8月初至10月上旬洪峰前后漢江主要控制站水位流量演變(a)丹江口水庫,(b)皇莊,(c)沙洋,(d)岳口,(e)仙桃,(f)漢川Fig.5 The temporal variation of water level (blue line), inflow and flow (purple line), outflow (yellow line) with its alert water level (orange line) before and after the flood peak in the Hanjing River Basin in 2021(a) Danjiangkou Reservoir, (b) Huangzhuang, (c) Shayang, (d) Yuekou, (e) Xiantao, (f) Hanchuan

    從中下游干流主要水文控制站的水位流量演變(圖5b~5f)可以看出,受漢江上游來水影響,皇莊以下主要水文控制站點(diǎn)在8月下旬至9月上中旬、9月下旬至10月上旬均存在2個(gè)超警戒階段,其中皇莊、沙洋、岳口、仙桃、漢川累計(jì)超警時(shí)間分別為151、255、228、232、580 h,最大超警幅度分別為0.29、0.40、0.50、0.53、1.56 m,呈現(xiàn)越到下游超警時(shí)間越長、超警幅度越大的特征,這與中下游河道落差小,斷面愈下愈小,導(dǎo)致的泄洪能力不足有關(guān)。此外,由于長江上游秋汛也較為明顯,長江上游來水造成漢口水位升高,也對(duì)漢江來水形成頂托作用。

    3.2 2021年漢江流域水情與2011年、2017年的對(duì)比

    2010年以來,漢江秋汛比較嚴(yán)重的年份還出現(xiàn)在2011年(蔡薌寧等,2012;柳艷菊等,2012;蔡新玲等,2013)、2017年(梁萍等,2019;李玉榮等,2017),其中2011年為減輕漢江上游來水對(duì)下游城市尤其是對(duì)武漢的影響,曾采用杜家臺(tái)分蓄洪區(qū)進(jìn)行分洪(劉志文,2012)。對(duì)比此3年丹江口水庫超過10 000 m3·s-1洪峰出現(xiàn)次數(shù)可見(表1),2011年次數(shù)最少,僅為2次;2017年為3次;2021年最多,達(dá)7次。從首次洪峰出現(xiàn)時(shí)間來看,2021年最早,在8月末就有洪峰出現(xiàn),而2011年、2017年則分別出現(xiàn)在9月中旬和9月下旬。在3個(gè)典型秋汛年中,最大洪峰流量出現(xiàn)在2011年9月19日,達(dá)26 600 m3·s-1,且此次過程流量漲幅也最大,達(dá)22 980 m3·s-1,2021年9月29日洪峰次之,為24 900 m3·s-1,流量過程漲幅居第二位,為21 850 m3·s-1,2017年的洪峰流量均在20 000 m3·s-1以下。從丹江口水庫洪峰水位來看,2021年普遍高于2011年和2017年。

    表1 2021、2017、2011年丹江口水庫水文特征值Table 1 Key hydrological values of the Danjiangkou Reservoir in 2021, 2017 and 2011

    對(duì)2011、2017、2021年漢江中下游主要水位控制站的水文特征值(表2)進(jìn)行對(duì)比,從超警戒站點(diǎn)數(shù)量來看,2021年、2017年皇莊以下江段均為全線超警,但在2011年超警戒站點(diǎn)略少,警戒主要發(fā)生在沙洋以下江段。從最高水位的情況來看,皇莊站以2017年為最高,為48.62 m;沙洋、岳口、仙桃站均是以2011年為最高;而漢川站,則以2021年為最高,為30.56 m,2011年次之,為30.55 m。從累計(jì)超警戒時(shí)長來看,均是以2021年為最長,但從超保證時(shí)長來看,在2011年仙桃站還出現(xiàn)了超保證的情況。

    表2 2021、2017、2011年漢江中下游主要水文控制站特征值Table 2 Key hydrological values of the control stations in the middle and lower reaches of Hanjiang River in 2021, 2017 and 2011

    對(duì)比此3年的水情特征后發(fā)現(xiàn),2011年漢江上游洪峰集中發(fā)生于9月中旬,洪峰次數(shù)最少(2次),但流量漲幅迅猛,洪峰流量及過程流量漲幅最大均出現(xiàn)在這一年。此外,2011年漢江中下游水位高、流量大,其超警戒站數(shù)及持續(xù)時(shí)間雖不及2021年,但沙洋、岳口、仙桃等站水位均為3年中最高,且出現(xiàn)了1個(gè)站超保證;2017年漢江上游洪峰發(fā)生在9月末至10月中旬前期,次數(shù)為3次,略高于2011年,但在峰值流量、中下游主要站點(diǎn)洪峰水位及超警戒時(shí)長上來看,普遍不及2011年和2021年;2021年漢江上游洪峰出現(xiàn)在8月下旬至10月上旬,洪峰出現(xiàn)時(shí)間最早,洪峰次數(shù)最多(7次),上游高水位流量持續(xù)時(shí)間最長,中下游累計(jì)超警時(shí)間也最長,其下游的漢川站水位為3年中最高,其超警幅度也最大。

    4 2021年漢江秋汛的異常環(huán)流特征

    4.1 500 hPa環(huán)流異常特征

    秋汛的形成與大尺度環(huán)流的異常穩(wěn)定維持密切相關(guān)。2021年8月21日至10月7日(圖6,等值線),極渦呈偏心型,其中心位于78°N、180°E附近,亞歐中高緯呈現(xiàn)“兩槽一脊”環(huán)流型,我國東北及其以北地區(qū)有高壓脊維持,巴爾喀什湖—貝加爾湖地區(qū)、日本以東—白令海地區(qū)各有一低槽發(fā)展。中低緯地區(qū),副高較歷史同期明顯偏強(qiáng)、偏西,其平均西伸脊點(diǎn)位于27°N、100°E附近。對(duì)應(yīng)的距平場(chǎng)顯示(圖6,填色),歐亞地區(qū)中高緯從西南至東北呈現(xiàn)經(jīng)向“-+-”波列分布,即巴爾喀什湖—貝加爾湖地區(qū)為明顯的負(fù)異常,歐洲大陸北部為顯著的正異常,白令海峽以北的大部地區(qū)為顯著的負(fù)異常,其負(fù)異常中心較歷史同期差值在-7 dagpm以上。此外,緯向上同樣存在明顯的“-+-”波列分布,即巴爾喀什湖—貝加爾湖地區(qū)為負(fù)異常,我國大陸東部至西北太平洋為正異常,日本島以東地區(qū)為負(fù)異常。中高緯地區(qū)的環(huán)流異常有利于冷空氣在巴爾喀什湖—貝加爾湖地區(qū)堆積并不斷分裂南下,與異常偏西的副高外圍暖濕氣流交匯于漢江上游地區(qū)。2021年在巴爾喀什湖—貝加爾湖地區(qū)高度場(chǎng)的負(fù)異常、歐洲北部的正異常、我國大陸東部至西北太平洋的正異常,與白虎志和董文杰(2004)、羅霄等(2013)、王春學(xué)等(2015)研究的華西秋雨偏強(qiáng)年的典型大尺度異常環(huán)流相比,基本一致,主要差異表現(xiàn)在低緯度孟加拉灣至南海地區(qū)高度場(chǎng)仍以正距平為主,而非典型環(huán)流中的明顯負(fù)距平。

    注:紅色實(shí)線:漢江流域邊界。圖6 2021年8月21日至10月7日500 hPa平均位勢(shì)高度(等值線,單位:dagpm)及距平場(chǎng)(填色)Fig.6 The 500 hPa mean geopotential height (isoline, unit: dagpm) and anomaly field (colored) from 21 August to 7 October 2021

    4.2 西太平洋副高異常特征

    500 hPa副高的變化是影響秋雨多少的主要因子之一,從2021年8月21日至10月7日期間逐日副高強(qiáng)度及面積指數(shù)演變(圖7a)可見,副高異常偏強(qiáng)、面積偏大,秋汛期間平均的強(qiáng)度指數(shù)和面積指數(shù)分別為歷史同期的1.5倍和1.7倍。從秋汛期間副高脊線、西伸脊點(diǎn)演變可見(圖7b),副高呈現(xiàn)出階段性的南北及東西擺動(dòng),對(duì)比逐日降水演變可見,副高西伸脊線在26°~28°N、西伸脊點(diǎn)在90°~110°E的階段,也是漢江上游降水顯著增強(qiáng)的階段。

    圖7 2021年8月21日至10月7日逐日500 hPa(a)副高面積指數(shù)和強(qiáng)度指數(shù)及(b)西伸脊點(diǎn)和脊線演變Fig.7 Daily variation of (a) subtropical high intensity and area index and (b) westward extension ridge point and ridge line at 500 hPa from 21 August to 7 October 2021

    4.3 東亞副熱帶西風(fēng)急流和南亞高壓

    南亞高壓及其北側(cè)的東亞副熱帶西風(fēng)急流是中高緯地區(qū)對(duì)流層中上層大氣環(huán)流的重要組成部分,其對(duì)大范圍雨帶形成及維持具有重要作用(張瓊和吳國雄,2001;況雪源和張耀存,2006;陸曉娟等,2022)。從歷史平均的8月21日至10月7日200 hPa 南亞高壓及平均緯向風(fēng)的歷史平均分布特征(圖8a)可見,南亞高壓1252 dagpm所控制的范圍很小,僅在高原西側(cè)90°E附近區(qū)域東亞副熱帶西風(fēng)急流軸(u≥30 m·s-1)呈準(zhǔn)東西向帶狀分布,在90°E、150°E存在兩個(gè)35 m·s-1以上的急流中心。從2021年8月21日至10月7日秋汛期間200 hPa南亞高壓及平均緯向風(fēng)(圖8b)分布可見,南亞高壓較歷史同期異常偏強(qiáng)、范圍偏大,1252 dagpm覆蓋了20°~130°E的廣大地區(qū),且其中心值超過1256 dagpm。東亞副熱帶西風(fēng)急流較歷史同期范圍更廣,覆蓋了37°~42°N、 35°~160°E大部區(qū)域,并且強(qiáng)度更強(qiáng),在60°~140°E的大范圍區(qū)域內(nèi)其中心強(qiáng)度均超過了35 m·s-1。從2021年秋汛期105°~112°E平均的200 hPa位勢(shì)高度、緯向風(fēng)及水平散度隨時(shí)間演變來看(圖8c),南亞高壓及東亞副熱帶西風(fēng)急流經(jīng)歷了4次加強(qiáng)北抬隨后又減弱南撤的過程,當(dāng)東亞副熱帶西風(fēng)急流軸與南亞高壓脊線分別穩(wěn)定在35°~40°N、25°~30°N時(shí),在30°~35°N(即漢江上游所在區(qū)域)有明顯的高空輻散發(fā)展,對(duì)應(yīng)的降水也進(jìn)入顯著增強(qiáng)階段。

    注:圖a,b中紅色實(shí)線:漢江流域邊界;圖c中紅色實(shí)線:南亞高壓脊線。圖8 (a)1991—2021年歷史平均和(b)2021年8月21日至10月7日200 hPa平均位勢(shì)高度(黑色實(shí)線,單位:dagpm,綠色實(shí)線:1252 dagpm)及緯向風(fēng)(填色:≥30 m·s-1),(c)2021年105°~112°E的200 hPa平均位勢(shì)高度(黑色虛線,單位:dagpm)、緯向風(fēng)(黑色實(shí)線,單位:m·s-1)和水平散度(填色)的時(shí)間-緯度剖面Fig.8 (a, b) The 200 hPa average geopotential height (black solid line, unit: dagpm: green solid line: 1252 dagpm) and u-wind (colored: ≥30 m·s-1) from 21 August to 7 October (a) of climate mean (1991-2021) and (b) in 2021, (c) the time-latitude cross-section of geopotential height (black dashed line, unit: dagpm), u-wind (black solid line, unit: m·s-1) and horizontal divergence (colored)

    4.4 低層流場(chǎng)及水汽輸送異常

    由2021年8月21日至10月7日700 hPa距平流場(chǎng)(圖9)可見,在中高緯度巴爾喀什湖—貝加爾湖地區(qū)有異常的氣旋性環(huán)流存在,而我國東部沿海地區(qū)處于與副高相關(guān)的異常反氣旋環(huán)流控制下,在其南側(cè)的0°~20°N、 100°~120°E區(qū)域內(nèi),可以看到存在明顯的氣旋性環(huán)流,其與越赤道氣流、高壓底部的偏東氣流匯合成一支偏東氣流,在進(jìn)入陸地后轉(zhuǎn)向形成一支異常的西南氣流,最終與中高緯氣旋性環(huán)流底部的偏北氣流在漢江上游地區(qū)交匯。

    注:D:氣旋性環(huán)流中心,G:反氣旋性環(huán)流中心,紅色實(shí)線:漢江流域邊界。圖9 2021年8月21日至10月7日700 hPa距平流場(chǎng)Fig.9 The 700 hPa streamline field of anomalies from 21 August to 7 October 2021

    從氣候平均的8月21日至10月7日700 hPa水汽通量及水汽通量散度(圖10a)分布特征可以看出,水汽輻合區(qū)主要位于28°~40°N、 80°~100°E附近,即漢江流域以西的四川盆地附近地區(qū),此外在漢江流域以東地區(qū)也存在一個(gè)弱水汽輻合中心,漢江上游流域無明顯的水汽輻合。但在2021年秋汛期(圖10b),水汽輻合中心范圍明顯增大,東擴(kuò)至118°E附近,漢江上游地區(qū)處于明顯的水汽輻合區(qū)中,這與李多和顧薇(2022)的研究結(jié)論較為一致。水汽的輸送路徑主要有4支,來自南海、阿拉伯海的水汽首先在0°~20°N、 100°~120°E形成一個(gè)水汽輻合中心,而后經(jīng)由副高南側(cè)的偏東南氣流與孟加拉灣水汽匯合成一條西南水汽輸送帶,最終與來自西風(fēng)帶的水汽在漢江上游形成強(qiáng)烈的輻合,為秋汛期持續(xù)降水提供了充分的水汽供應(yīng)。

    注:紅色實(shí)線:漢江流域邊界。圖10 (a)1991—2021年歷史平均和(b)2021年8月21日至10月7日700 hPa平均水汽通量(箭矢,單位:kg·m-1·s-1)及水汽通量散度(填色,單位:10-7 g·cm-2·hPa-1·s-1)Fig.10 The 700 hPa water vapor flux (arrow, unit: kg·m-1·s-1) and water vapor flux divergence (colored, unit: 10-7 g·cm-2·hPa-1·s-1) from 21 August to 7 October (a) of climate mean (1991-2021) and (b) in 2021

    5 2021年漢江秋汛成因概念模型

    綜合上文對(duì)2021年漢江秋汛水文氣象演變特征和大氣環(huán)流異常分析結(jié)果,得到2021年漢江秋汛成因概念模型(圖11)。

    圖11 2021年漢江秋汛成因示意圖Fig.11 Schematic diagram of causes for autumn floods in Hanjiang River in 2021

    當(dāng)200 hPa南亞高壓、東亞副熱帶西風(fēng)急流較同期異常偏強(qiáng);500 hPa巴爾喀什湖—貝加爾湖地區(qū)有低槽發(fā)展,其距平場(chǎng)呈明顯負(fù)距平,我國大陸東部有高壓脊發(fā)展、副高異常偏西偏北,且中東部大部地區(qū)在距平場(chǎng)上呈明顯正距平;對(duì)流層中低層,漢江上游地區(qū)有異常的水汽輻合中心發(fā)展,并且上述異常大尺度環(huán)流配置穩(wěn)定維持時(shí),有利于漢江上游出現(xiàn)大范圍持續(xù)性降水。當(dāng)東亞副熱帶西風(fēng)急流軸位于35°~40°N,南亞高壓脊線位于25°~30°N時(shí),副高脊線在26°~28°N、西伸脊點(diǎn)在90°~110°E時(shí),有利于漢江上游降水顯著增強(qiáng)。在丹江口水庫入庫流量接近2000 m3·s-1后,當(dāng)漢江上游5 d累計(jì)面雨量達(dá)60 mm以上時(shí),往往造成漢江上游流域10 000 m3·s-1以上的洪峰過程。漢江上游洪峰形成后,在其向中下游的移動(dòng)過程中,由于河道落差小,斷面逐級(jí)收窄,導(dǎo)致降水并不明顯的漢江中下游地區(qū)出現(xiàn)超警戒洪水。

    6 結(jié) 論

    本文分析了2021年漢江秋汛降水的時(shí)空分布及歷史極端性特征,漢江流域主要水文控制站水情演變及大尺度環(huán)流異常特征,并初步探討了致洪的環(huán)流特征,得出以下結(jié)論:

    (1)2021年漢江秋汛致洪降水主要分布在漢江上游地區(qū),具有開始時(shí)間早、持續(xù)時(shí)間長、降水的累計(jì)強(qiáng)度大及歷史極端性強(qiáng)等特征,降水自8月21日開始發(fā)展,至10月上旬方減弱,持續(xù)近50 d,期間經(jīng)歷了4個(gè)降水的集中時(shí)段。漢江上游3個(gè)子流域石泉以上、石泉安康、安康丹江口累計(jì)面雨量分別為1961年以來的第一、第一、第二位,且整個(gè)漢江上游流域累計(jì)面雨量居1961年以來第一位。

    (2)受漢江上游降水影響,丹江口水庫出現(xiàn)10 000 m3·s-1以上的洪峰過程7次,最大洪峰流量出現(xiàn)在9月29日,為24 900 m3·s-1。漢江中下游主要水文控制站點(diǎn)累計(jì)超警時(shí)間在151~580 h,最大超警幅度分別為0.29~1.56 m,呈現(xiàn)越到下游超警時(shí)間越長、超警幅度越大的特征。與典型秋汛年水情的對(duì)比分析表明,2021年水情的異常程度明顯超過2011年、2017年。2021年漢江上游洪峰出現(xiàn)在8月下旬至10月上旬,洪峰出現(xiàn)時(shí)間最早,洪峰次數(shù)最多(7次),上游高水位流量持續(xù)時(shí)間最長,中下游累計(jì)超警時(shí)間也最長,其下游的漢川站水位為3年中最高,超警幅度也最大。

    (3)副高異常偏西偏北,是導(dǎo)致秋汛期間漢江流域降水“西多東少”的主要原因。在同期異常偏強(qiáng)的南亞高壓、東亞副熱帶西風(fēng)急流及對(duì)流層中低層異常的西南氣流的共同作用下,漢江上游地區(qū)形成了持續(xù)的高空輻散及低層水汽輻合中心,這為降水的長期維持起到重要作用。當(dāng)副高脊線在26°~28°N、西伸脊點(diǎn)在90°~110°E,且西風(fēng)急流軸與南亞高壓脊線分別穩(wěn)定在35°~40°N、25°~30°N時(shí),對(duì)應(yīng)的漢江上游降水也進(jìn)入顯著增強(qiáng)階段。當(dāng)?shù)そ谒畮烊霂炝髁拷咏?000 m3·s-1后,漢江上游60 mm 以上的5 d累計(jì)面雨量往往能造成10 000 m3·s-1以上的洪峰過程。上游來水增大后,漢江中下游河道逐級(jí)收窄,是導(dǎo)致降水并不明顯的漢江中下游地區(qū)出現(xiàn)超警洪水的主要原因。

    猜你喜歡
    丹江口水庫洪峰漢江
    漢江春曉
    南風(fēng)(2021年32期)2021-12-31 05:57:16
    漢江,為你梳妝
    守望漢江
    漢江之歌
    丹江口水庫的水資源管理與保護(hù)概述
    水庫移民農(nóng)地流轉(zhuǎn)與生計(jì)恢復(fù)問題及對(duì)策研究
    淡定!
    基于ETM+數(shù)據(jù)的水體信息提取
    綠色科技(2017年10期)2017-07-05 12:18:01
    解禁洪峰
    丹江口水庫河南外遷農(nóng)村移民安置效果研究
    科技資訊(2017年7期)2017-05-06 11:42:44
    色5月婷婷丁香| 久久久精品区二区三区| 国产亚洲av片在线观看秒播厂| 国产精品国产av在线观看| 伦理电影大哥的女人| 久久99精品国语久久久| 多毛熟女@视频| 欧美日韩国产mv在线观看视频| 三级国产精品片| 久久久久久久亚洲中文字幕| 寂寞人妻少妇视频99o| 菩萨蛮人人尽说江南好唐韦庄| 九色成人免费人妻av| 国语对白做爰xxxⅹ性视频网站| 在线观看www视频免费| 成人午夜精彩视频在线观看| 一级,二级,三级黄色视频| 亚洲精品色激情综合| 亚洲高清免费不卡视频| 亚洲不卡免费看| 91精品三级在线观看| 亚洲精品aⅴ在线观看| 欧美人与性动交α欧美精品济南到 | 好男人视频免费观看在线| 性色av一级| 国产午夜精品一二区理论片| 亚洲精品日韩在线中文字幕| 青春草亚洲视频在线观看| 欧美精品亚洲一区二区| 少妇被粗大的猛进出69影院 | 免费人妻精品一区二区三区视频| 99九九线精品视频在线观看视频| 丰满饥渴人妻一区二区三| 亚洲人成网站在线观看播放| 两个人的视频大全免费| 国国产精品蜜臀av免费| 一边摸一边做爽爽视频免费| 两个人免费观看高清视频| 高清黄色对白视频在线免费看| av免费在线看不卡| av不卡在线播放| av在线app专区| 国产白丝娇喘喷水9色精品| av播播在线观看一区| 制服人妻中文乱码| 欧美亚洲日本最大视频资源| 热99久久久久精品小说推荐| 日韩av在线免费看完整版不卡| 精品久久蜜臀av无| 美女cb高潮喷水在线观看| 久久精品国产鲁丝片午夜精品| 九草在线视频观看| 极品人妻少妇av视频| 精品卡一卡二卡四卡免费| 国产一级毛片在线| 欧美日本中文国产一区发布| 午夜av观看不卡| 亚洲欧美成人综合另类久久久| 一区在线观看完整版| 婷婷色麻豆天堂久久| 一区在线观看完整版| av免费在线看不卡| 中文天堂在线官网| 18在线观看网站| 亚洲欧美精品自产自拍| 最新的欧美精品一区二区| 成人国产麻豆网| 午夜激情久久久久久久| 久久青草综合色| 色视频在线一区二区三区| 亚洲精品av麻豆狂野| 校园人妻丝袜中文字幕| 我要看黄色一级片免费的| 亚洲精品成人av观看孕妇| a级毛色黄片| 午夜精品国产一区二区电影| av在线app专区| 亚洲精华国产精华液的使用体验| 亚洲久久久国产精品| 精品人妻熟女毛片av久久网站| 国产免费现黄频在线看| 国产精品一区二区三区四区免费观看| 男女啪啪激烈高潮av片| 日韩强制内射视频| 亚洲av中文av极速乱| 亚洲精品乱久久久久久| 狂野欧美白嫩少妇大欣赏| 久久久精品免费免费高清| 免费观看在线日韩| a级毛片免费高清观看在线播放| 老司机亚洲免费影院| 日本欧美视频一区| freevideosex欧美| 黄色怎么调成土黄色| 美女视频免费永久观看网站| videosex国产| 国产欧美另类精品又又久久亚洲欧美| 麻豆乱淫一区二区| 国产精品一区二区三区四区免费观看| 精品久久蜜臀av无| 秋霞在线观看毛片| 人人妻人人爽人人添夜夜欢视频| 极品人妻少妇av视频| 一本大道久久a久久精品| 伊人久久国产一区二区| 亚洲中文av在线| 五月伊人婷婷丁香| 亚洲欧美日韩卡通动漫| 久久精品国产亚洲av天美| 全区人妻精品视频| 亚洲欧洲国产日韩| 一本久久精品| 欧美97在线视频| 99国产综合亚洲精品| 又大又黄又爽视频免费| 免费观看的影片在线观看| 国产免费又黄又爽又色| 男女边摸边吃奶| 久久影院123| 欧美3d第一页| 日本黄色日本黄色录像| 国产精品国产三级国产av玫瑰| 中文字幕精品免费在线观看视频 | 全区人妻精品视频| 午夜久久久在线观看| 99热全是精品| 国产白丝娇喘喷水9色精品| 一级毛片 在线播放| 日本欧美国产在线视频| 尾随美女入室| 色吧在线观看| 久久久久精品性色| 在线播放无遮挡| 一区在线观看完整版| 国产免费一区二区三区四区乱码| 考比视频在线观看| 最新中文字幕久久久久| 春色校园在线视频观看| 日韩成人伦理影院| av播播在线观看一区| 国产免费视频播放在线视频| 高清毛片免费看| 考比视频在线观看| 亚洲怡红院男人天堂| 中文字幕最新亚洲高清| 欧美丝袜亚洲另类| 精品熟女少妇av免费看| 久久久久视频综合| 亚洲欧美日韩另类电影网站| 国产在线免费精品| 有码 亚洲区| 99热这里只有是精品在线观看| 欧美日本中文国产一区发布| 一级片'在线观看视频| 秋霞在线观看毛片| 全区人妻精品视频| 亚洲欧洲精品一区二区精品久久久 | 91精品一卡2卡3卡4卡| 欧美日韩视频精品一区| 在线观看免费高清a一片| 乱码一卡2卡4卡精品| 全区人妻精品视频| 精品少妇久久久久久888优播| av黄色大香蕉| 欧美精品一区二区大全| 色哟哟·www| 亚洲精品456在线播放app| .国产精品久久| 九色成人免费人妻av| 亚洲成人一二三区av| 成年人免费黄色播放视频| 亚洲怡红院男人天堂| 下体分泌物呈黄色| 看十八女毛片水多多多| 啦啦啦在线观看免费高清www| 少妇熟女欧美另类| 色网站视频免费| 久久久精品94久久精品| 极品少妇高潮喷水抽搐| 91久久精品电影网| 高清不卡的av网站| 日本91视频免费播放| 午夜日本视频在线| 人妻人人澡人人爽人人| 91精品一卡2卡3卡4卡| 精品久久久久久电影网| 最近2019中文字幕mv第一页| 日本黄色日本黄色录像| 免费av中文字幕在线| 国产无遮挡羞羞视频在线观看| 99视频精品全部免费 在线| √禁漫天堂资源中文www| 欧美亚洲日本最大视频资源| 天天躁夜夜躁狠狠久久av| 日本爱情动作片www.在线观看| 久久久国产欧美日韩av| 欧美日韩视频精品一区| 最近2019中文字幕mv第一页| 日本免费在线观看一区| 免费观看性生交大片5| av在线老鸭窝| 国产一区二区在线观看日韩| 久久久久国产精品人妻一区二区| 波野结衣二区三区在线| 一级二级三级毛片免费看| 亚洲国产精品国产精品| 亚洲色图 男人天堂 中文字幕 | 免费黄色在线免费观看| 国产熟女欧美一区二区| 91久久精品国产一区二区三区| 老司机影院成人| 在线观看一区二区三区激情| 欧美日韩综合久久久久久| 亚州av有码| 中文字幕精品免费在线观看视频 | 大香蕉久久成人网| 欧美 亚洲 国产 日韩一| 在线播放无遮挡| 我的女老师完整版在线观看| 在线亚洲精品国产二区图片欧美 | 欧美日韩亚洲高清精品| 国产免费福利视频在线观看| 好男人视频免费观看在线| 蜜桃在线观看..| 亚洲国产av新网站| 国产无遮挡羞羞视频在线观看| 免费观看在线日韩| 黑人欧美特级aaaaaa片| 性色av一级| 国产深夜福利视频在线观看| av有码第一页| 国产伦精品一区二区三区视频9| av在线播放精品| 99热6这里只有精品| av专区在线播放| 免费高清在线观看日韩| 精品久久久久久久久亚洲| 国产高清国产精品国产三级| 亚洲精品第二区| 免费看不卡的av| 三级国产精品欧美在线观看| 亚洲av男天堂| 久久韩国三级中文字幕| 18禁动态无遮挡网站| 亚洲无线观看免费| 满18在线观看网站| 国产视频首页在线观看| 午夜影院在线不卡| 亚洲伊人久久精品综合| 久久久久久久国产电影| 久久免费观看电影| 2021少妇久久久久久久久久久| 不卡视频在线观看欧美| 国产亚洲最大av| 国产日韩欧美亚洲二区| 26uuu在线亚洲综合色| 人妻少妇偷人精品九色| 亚洲欧美精品自产自拍| 性色avwww在线观看| 午夜福利,免费看| 午夜激情久久久久久久| 少妇精品久久久久久久| 亚洲欧美成人综合另类久久久| 人妻系列 视频| 大话2 男鬼变身卡| 国产免费视频播放在线视频| kizo精华| 精品久久久久久久久av| 黄色欧美视频在线观看| 七月丁香在线播放| 丝瓜视频免费看黄片| 2022亚洲国产成人精品| 亚洲欧美中文字幕日韩二区| 国产精品一区二区在线观看99| 国产精品女同一区二区软件| 日韩在线高清观看一区二区三区| 午夜久久久在线观看| 精品亚洲成国产av| 亚洲精品日韩av片在线观看| 最近2019中文字幕mv第一页| 欧美xxⅹ黑人| 一本久久精品| 国产av精品麻豆| 人人妻人人添人人爽欧美一区卜| av视频免费观看在线观看| 七月丁香在线播放| av在线app专区| 婷婷色综合大香蕉| 久久精品国产亚洲网站| 久久久久久久久久久久大奶| av女优亚洲男人天堂| 亚洲人成77777在线视频| 丰满迷人的少妇在线观看| 亚洲第一av免费看| 亚洲精品乱久久久久久| 精品亚洲成国产av| 国产成人免费观看mmmm| 亚洲精品乱码久久久v下载方式| 亚洲美女黄色视频免费看| 国产片特级美女逼逼视频| 免费高清在线观看日韩| 久久午夜福利片| 亚洲在久久综合| 国国产精品蜜臀av免费| 精品人妻在线不人妻| 成人国产av品久久久| 欧美日韩av久久| 夜夜看夜夜爽夜夜摸| 亚洲色图综合在线观看| 国产日韩欧美在线精品| 亚洲美女黄色视频免费看| 亚洲高清免费不卡视频| 考比视频在线观看| 亚洲怡红院男人天堂| 丝袜在线中文字幕| 超色免费av| 亚洲国产精品成人久久小说| 三级国产精品欧美在线观看| 在线 av 中文字幕| 天天操日日干夜夜撸| 一区二区三区乱码不卡18| 亚洲国产日韩一区二区| 亚洲精品av麻豆狂野| 国产精品99久久久久久久久| 欧美性感艳星| 国产伦理片在线播放av一区| av天堂久久9| 青春草视频在线免费观看| 26uuu在线亚洲综合色| 青春草国产在线视频| 22中文网久久字幕| 下体分泌物呈黄色| 日韩人妻高清精品专区| 欧美精品亚洲一区二区| 日韩av在线免费看完整版不卡| av女优亚洲男人天堂| 亚洲国产精品一区二区三区在线| av不卡在线播放| 高清午夜精品一区二区三区| 最近手机中文字幕大全| 亚洲成人av在线免费| 美女内射精品一级片tv| 美女脱内裤让男人舔精品视频| 色视频在线一区二区三区| 少妇精品久久久久久久| 在线观看三级黄色| 国产免费又黄又爽又色| 国产男女内射视频| 在线观看免费视频网站a站| 我的女老师完整版在线观看| 少妇熟女欧美另类| 麻豆成人av视频| 亚洲欧美成人综合另类久久久| 韩国av在线不卡| 天堂中文最新版在线下载| 18禁在线播放成人免费| 国产高清三级在线| 这个男人来自地球电影免费观看 | 国精品久久久久久国模美| 国产亚洲精品第一综合不卡 | 亚洲av欧美aⅴ国产| 18禁在线播放成人免费| 啦啦啦啦在线视频资源| 日本黄色日本黄色录像| 久久久久久伊人网av| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲av涩爱| 最近最新中文字幕免费大全7| 欧美3d第一页| 国产有黄有色有爽视频| 精品久久国产蜜桃| 国产 精品1| 精品少妇黑人巨大在线播放| 亚洲怡红院男人天堂| 黄色视频在线播放观看不卡| 精品少妇黑人巨大在线播放| 大码成人一级视频| 麻豆精品久久久久久蜜桃| 成人毛片60女人毛片免费| 成人漫画全彩无遮挡| 国产毛片在线视频| 免费观看性生交大片5| 日韩,欧美,国产一区二区三区| 国产日韩欧美亚洲二区| 久久久久视频综合| 亚洲精品久久成人aⅴ小说 | 在线观看一区二区三区激情| 国产精品国产三级国产av玫瑰| 热re99久久精品国产66热6| 国产免费一级a男人的天堂| 国产精品嫩草影院av在线观看| 伦精品一区二区三区| 久久久久久久精品精品| 国产在线免费精品| 一区二区三区乱码不卡18| 中文字幕人妻丝袜制服| 春色校园在线视频观看| 中文天堂在线官网| 建设人人有责人人尽责人人享有的| 国产午夜精品久久久久久一区二区三区| 亚洲成人一二三区av| 国产一区二区在线观看日韩| 日日啪夜夜爽| 黑人高潮一二区| 天天影视国产精品| 高清av免费在线| 赤兔流量卡办理| 制服人妻中文乱码| 国产成人精品福利久久| 日韩欧美精品免费久久| 午夜福利在线观看免费完整高清在| 亚洲四区av| av网站免费在线观看视频| 少妇 在线观看| 2022亚洲国产成人精品| 如何舔出高潮| 久久久久精品久久久久真实原创| 高清在线视频一区二区三区| 黑人高潮一二区| 在线观看免费视频网站a站| 男女国产视频网站| 日韩,欧美,国产一区二区三区| 国产精品久久久久久av不卡| 亚洲欧美精品自产自拍| 亚洲五月色婷婷综合| 街头女战士在线观看网站| 女的被弄到高潮叫床怎么办| 美女国产视频在线观看| 国产日韩欧美在线精品| 搡老乐熟女国产| 国产精品99久久久久久久久| 精品卡一卡二卡四卡免费| 精品熟女少妇av免费看| 亚洲五月色婷婷综合| 欧美日韩在线观看h| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 午夜福利视频在线观看免费| 欧美另类一区| 日韩,欧美,国产一区二区三区| 18禁在线无遮挡免费观看视频| 免费观看av网站的网址| 亚洲精品第二区| 一区在线观看完整版| 国产在线一区二区三区精| av又黄又爽大尺度在线免费看| 国产免费又黄又爽又色| 国产日韩一区二区三区精品不卡 | 国产日韩一区二区三区精品不卡 | 18禁在线播放成人免费| 色94色欧美一区二区| 国产亚洲精品久久久com| 久久免费观看电影| 嘟嘟电影网在线观看| 成年女人在线观看亚洲视频| 亚洲精品乱久久久久久| 国产成人午夜福利电影在线观看| 亚洲,一卡二卡三卡| 人体艺术视频欧美日本| 狂野欧美激情性bbbbbb| 日韩中字成人| 免费少妇av软件| 97在线视频观看| 色5月婷婷丁香| 国产精品久久久久久av不卡| 两个人的视频大全免费| 亚洲精品中文字幕在线视频| 久久久久久伊人网av| 精品视频人人做人人爽| 日韩强制内射视频| 高清视频免费观看一区二区| 精品国产国语对白av| 天堂中文最新版在线下载| 99热6这里只有精品| 国产一区二区在线观看日韩| 亚洲第一av免费看| 久久人人爽av亚洲精品天堂| 成人亚洲精品一区在线观看| 免费日韩欧美在线观看| 街头女战士在线观看网站| 亚洲图色成人| √禁漫天堂资源中文www| 99国产精品免费福利视频| 国模一区二区三区四区视频| 永久免费av网站大全| 麻豆精品久久久久久蜜桃| 成人国产麻豆网| 欧美97在线视频| 黄片播放在线免费| 久久久久久久久久人人人人人人| 亚洲人与动物交配视频| 国产国语露脸激情在线看| av在线老鸭窝| 黑人欧美特级aaaaaa片| 人体艺术视频欧美日本| 韩国av在线不卡| 国产精品秋霞免费鲁丝片| 美女中出高潮动态图| 色视频在线一区二区三区| 九九在线视频观看精品| 午夜老司机福利剧场| 九九在线视频观看精品| 日韩免费高清中文字幕av| 波野结衣二区三区在线| 久久97久久精品| 久久鲁丝午夜福利片| 亚洲色图综合在线观看| 久久久久久久久久成人| 国产高清国产精品国产三级| 热99久久久久精品小说推荐| a级毛片免费高清观看在线播放| 久久精品熟女亚洲av麻豆精品| 啦啦啦在线观看免费高清www| 日本欧美视频一区| 爱豆传媒免费全集在线观看| 欧美精品高潮呻吟av久久| 国产一区亚洲一区在线观看| 亚洲欧美色中文字幕在线| 国产一区亚洲一区在线观看| 欧美一级a爱片免费观看看| 自拍欧美九色日韩亚洲蝌蚪91| 一个人免费看片子| 国产欧美另类精品又又久久亚洲欧美| 欧美精品高潮呻吟av久久| 自拍欧美九色日韩亚洲蝌蚪91| 成年女人在线观看亚洲视频| 午夜久久久在线观看| 国产av一区二区精品久久| 搡老乐熟女国产| 人妻一区二区av| 亚洲精品久久午夜乱码| 十八禁网站网址无遮挡| 插阴视频在线观看视频| 在线观看免费日韩欧美大片 | 国产片内射在线| 国产成人精品在线电影| 一级黄片播放器| 美女国产视频在线观看| 日日摸夜夜添夜夜添av毛片| 一本久久精品| 日本av免费视频播放| 国产亚洲av片在线观看秒播厂| 国产日韩一区二区三区精品不卡 | 免费观看无遮挡的男女| 22中文网久久字幕| 久久精品国产a三级三级三级| 91久久精品国产一区二区三区| 免费高清在线观看日韩| 青春草国产在线视频| 一本色道久久久久久精品综合| 高清视频免费观看一区二区| 热99国产精品久久久久久7| 亚洲天堂av无毛| 免费人成在线观看视频色| 啦啦啦中文免费视频观看日本| 午夜福利网站1000一区二区三区| 亚洲精品国产av蜜桃| 欧美国产精品一级二级三级| 国产高清不卡午夜福利| 99热全是精品| 18禁裸乳无遮挡动漫免费视频| 亚洲久久久国产精品| 日韩亚洲欧美综合| 五月伊人婷婷丁香| 精品人妻在线不人妻| 精品国产国语对白av| 久久99精品国语久久久| 国产一区有黄有色的免费视频| 久久99热这里只频精品6学生| 亚洲国产精品一区二区三区在线| 在线亚洲精品国产二区图片欧美 | 伦理电影大哥的女人| 久久久久国产精品人妻一区二区| 日日撸夜夜添| 色5月婷婷丁香| 国产精品嫩草影院av在线观看| 日韩伦理黄色片| 亚洲精品久久成人aⅴ小说 | 看免费成人av毛片| 国产精品久久久久久精品电影小说| 少妇丰满av| 成人午夜精彩视频在线观看| 性色avwww在线观看| 欧美日韩视频精品一区| 在线观看免费视频网站a站| 亚洲婷婷狠狠爱综合网| 街头女战士在线观看网站| 美女中出高潮动态图| 丝袜美足系列| av播播在线观看一区| 精品视频人人做人人爽| 好男人视频免费观看在线| 精品一区二区三区视频在线| 精品国产一区二区久久| 久久毛片免费看一区二区三区| 婷婷色综合www| 国产一区二区在线观看日韩| 热99久久久久精品小说推荐| 国产精品久久久久成人av| 哪个播放器可以免费观看大片| 精品视频人人做人人爽| av国产久精品久网站免费入址| 国产精品一国产av| 精品人妻在线不人妻| 日韩免费高清中文字幕av| 最黄视频免费看| 欧美变态另类bdsm刘玥| 99久久精品国产国产毛片| 国产极品粉嫩免费观看在线 | 欧美精品国产亚洲| 亚洲成色77777| 建设人人有责人人尽责人人享有的| 亚洲av福利一区| 国产白丝娇喘喷水9色精品| 热re99久久精品国产66热6|