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

    風(fēng)暴潮、大潮對廣西潿洲島西南沙灘侵蝕的影響分析

    2015-01-05 02:32:38李明杰吳少華劉秋興董劍希
    海洋學(xué)報(bào) 2015年9期
    關(guān)鍵詞:風(fēng)暴潮大潮潮位

    李明杰,吳少華,劉秋興,董劍希

    (1.國家海洋環(huán)境預(yù)報(bào)中心,北京 100081)

    風(fēng)暴潮、大潮對廣西潿洲島西南沙灘侵蝕的影響分析

    李明杰1,吳少華1,劉秋興1,董劍希1

    (1.國家海洋環(huán)境預(yù)報(bào)中心,北京 100081)

    本文統(tǒng)計(jì)分析了廣西潿洲島沿海氣候、潮汐和風(fēng)暴潮等歷史資料,利用耿貝爾方法推算了潿洲島多年一遇年極值高潮位,并估算了其漫灘范圍分布,指出近幾年高潮位出現(xiàn)的頻次和極值均越來越高是潿洲島西南部沙灘侵蝕愈加嚴(yán)重的重要原因,最后結(jié)合風(fēng)暴潮-海浪耦合數(shù)值模擬了研究區(qū)域內(nèi)“0312”號臺風(fēng)風(fēng)暴潮漫灘的情況,分析了風(fēng)暴潮和大潮對潿洲島西南部沙灘侵蝕的影響,對當(dāng)?shù)匕稙┬迯?fù)和防護(hù)具有一定的指導(dǎo)意義。

    風(fēng)暴潮;大潮;沙灘侵蝕;耿貝爾方法;數(shù)值模擬

    1 引言

    潿洲島是廣西眾多島嶼中最大的海島,也是中國最大、地質(zhì)年齡最年輕的火山島。該島位于北部灣中部,北臨廣西北海市,東望雷州半島,東南與斜陽島毗鄰,南與海南島隔海相望,西面面向越南,距北海市區(qū)大約36海里(見圖1)。潿洲島南北方向的長度為6.5 km,東西方向?qū)? km,總面積24.74 km2,島的最高海拔79 m。潿洲島附近水深較淺,水深普遍小于10 m。在島嶼的東側(cè)及西北側(cè)分別有兩條深約30~40 m的海槽。

    潿洲島及其周邊海域的旅游資源、水產(chǎn)資源、油汽資源等非常豐富,自然環(huán)境優(yōu)越。近幾年來,潿洲島岸線整體遭受侵蝕,尤其西南部岸灘作為潿洲島重要的濱海浴場,侵蝕更為嚴(yán)重,2006-2013年間年均下蝕可達(dá)0.18 m[1]。潿洲島西南部潮間淺灘屬沙質(zhì)海灘,灘面坡度為1°~6°,總長度約2.5 km,寬50~150 m。沙灘等高線大體呈南北走向,地形呈東高西低。高潮線以上為木麻黃防風(fēng)林和山體,地形起伏較大。海岸侵蝕使得西南部海灘岸線后退,低灘下蝕致坡度變陡,沙灘下層礫石裸露,并在離岸區(qū)堆積成沿岸沙壩,臨海部分樹木枯死,對當(dāng)?shù)氐穆糜尉坝^造成了嚴(yán)重的影響。

    沿海岸灘侵蝕的因素多種多樣[2—5],既有自然因素如風(fēng)暴潮和近岸浪[6—10]、大潮、海平面上升[11—12]、入海泥沙減少[13]等,又有人為因素如采砂挖沙、破壞海灘植被或砍伐岸灘防護(hù)林、不合理的海岸工程建設(shè)等[14]。

    在諸多侵蝕因素中,風(fēng)暴潮和近岸浪對岸灘地形的破壞性最為嚴(yán)重?!?216”號臺風(fēng)風(fēng)暴潮期間,長時(shí)間的增水和大浪的作用給山東省沿海造成了嚴(yán)重災(zāi)害,沙質(zhì)海岸后退3~12 m,最大達(dá)30 m,損失土地4.66 km2[6—7]。1962年7月下旬至8月中旬,長江口連續(xù)遭受兩次強(qiáng)臺風(fēng)襲擊,某些岸段灘面平均刷低20~30 cm,岸線后退10~20 m[15]。不同于萊州灣等陸地海岸帶的侵蝕[13],潿洲島離陸地較遠(yuǎn)且沒有入海河流,沒有入海泥沙的影響,風(fēng)暴潮、潮流和近岸浪等是造成潿洲島海岸侵蝕的主要海洋動力[1]。風(fēng)暴潮過程以后,泥沙會從離岸區(qū)向岸區(qū)推移,但沙灘循環(huán)作用的不對稱性,沙灘沉積物總體會處于虧損狀態(tài),從而導(dǎo)致沙灘的侵蝕后退[9,16—17]。此外,冬季大風(fēng)引起的近岸浪配合天文潮大潮也會顯著改變當(dāng)?shù)匕稙┑默F(xiàn)狀。

    近幾年來,風(fēng)暴潮致海岸侵蝕的報(bào)道及其侵蝕過程的研究工作較多[8,10,16—19],利用詳實(shí)的歷史資料結(jié)合數(shù)值模式來分析風(fēng)暴潮和大潮的影響強(qiáng)度以及時(shí)空變化的研究尚不多。本文從風(fēng)暴潮和大潮的角度去探討其對潿洲島西南部岸灘侵蝕的影響。

    本文搜集了潿洲島氣候、潮汐和風(fēng)暴潮等歷史資料以及潿洲島附近精細(xì)化的水深地形和陸地高程數(shù)據(jù),并進(jìn)行了統(tǒng)計(jì)與計(jì)算分析,最后利用風(fēng)暴潮-海浪耦合模式模擬了潿洲島西南部沙灘風(fēng)暴潮漫灘的過程。

    圖1 模式模擬區(qū)域地形及區(qū)域內(nèi)典型潮位站分布Fig.1 Bathymetry of our model domain and the distribution of typical tidal stations therein

    2 廣西潿洲島氣候以及潮汐特征分析

    2.1 廣西沿海氣候特征及相關(guān)分析

    廣西海岸帶地處低緯度,屬我國南亞熱帶季風(fēng)型海洋性氣候。海岸帶風(fēng)向具有明顯的季節(jié)性變化。每年9、10月至翌年3、4月受北方大陸干冷氣團(tuán)的控制,盛行偏北氣流,廣西各地影響的風(fēng)向均是北-東北風(fēng);從4、5月至8、9月,在海洋暖溫氣團(tuán)的主宰下,盛行偏南氣流,多吹西南-東東南風(fēng)[20]。

    夏秋兩季每年平均受2~4次臺風(fēng)影響。臺風(fēng)由南海進(jìn)入北部灣時(shí),因受海南島和雷州半島的阻擋,一般只有6~10級,風(fēng)速大于8級的大風(fēng)平均每年出現(xiàn)11.1 d。近幾年進(jìn)入北部灣的臺風(fēng)強(qiáng)度普遍偏強(qiáng),“1223”強(qiáng)臺風(fēng)“山神”、“1330”超強(qiáng)臺風(fēng)“海燕”和“1415”臺風(fēng)“海鷗”進(jìn)入北部灣風(fēng)速均達(dá)到了13級,“1409”超強(qiáng)臺風(fēng)“威馬遜”進(jìn)入北部灣更是達(dá)到16級風(fēng)速。

    王慧和隋偉輝[21]分析了北部灣多年平均大風(fēng)(6級以上)天數(shù)的季節(jié)性變化,結(jié)果顯示12月份為全年中大風(fēng)天數(shù)最多的月份,月均大風(fēng)天數(shù)大于6天的月份主要集中在10月份至翌年的1月,夏秋兩季(5-9月)雖然受臺風(fēng)的影響,最大風(fēng)速為全年中最大,但大風(fēng)天數(shù)卻為全年中最少。從這個(gè)角度來說,夏秋季臺風(fēng)影響期間,因風(fēng)速較大,可能會對研究區(qū)域產(chǎn)生突變式的影響;10月至翌年1月雖然很少受到臺風(fēng)影響,但因大風(fēng)天數(shù)較多,其影響也不容忽視,可能會對研究區(qū)域產(chǎn)生持續(xù)性的影響。

    從物質(zhì)輸運(yùn)的角度分析,10月至翌年1月主要盛行偏北風(fēng),此時(shí)淺海的風(fēng)海流體積輸運(yùn)為偏西南方向,即研究區(qū)域的物質(zhì)輸運(yùn)的趨勢是離岸、向南;5-9月若無臺風(fēng)影響,因盛行偏南風(fēng),此時(shí)淺海的風(fēng)海流體積輸運(yùn)為偏東北方向,即研究區(qū)域的物質(zhì)輸運(yùn)的趨勢是向岸、向北。這個(gè)結(jié)論恰好解釋了潿洲島西南部形成沙壩以后泥沙存在季節(jié)性南北方向搬運(yùn)的現(xiàn)象。

    2.2 廣西潿洲島潮汐特征及相關(guān)分析

    潿洲站基本分潮振幅由潮汐調(diào)和分析得出:HM2太陰半日分潮振幅為36.73 cm,HK1太陰太陽合成日分潮振幅84.78 cm,HO1太陰日分潮振幅92.71 cm。根據(jù)GB/T 17839-1999《警戒潮位核定方法》的潮汐類型判別標(biāo)準(zhǔn),潿洲站A=(HK1+HO1)/HM2值為4.83,A>4.0,因此潿洲站潮汐性質(zhì)是正規(guī)全日潮性質(zhì)。

    圖3 潿洲站1973-2014年高潮位年極值序列Fig.3 Annual extreme high tide of Weizhou Island during 1973 to 2014

    圖2給出了2012年潿洲島天文潮的情況。從圖中可以看出,不同于我國北方沿海夏季潮位高,冬季潮位低的情況,潿洲站5-7月以及10-12月為全年中的大潮期。這兩個(gè)時(shí)間段恰好分別是臺風(fēng)影響時(shí)段和大風(fēng)天數(shù)最多時(shí)段,因此高潮位配合大風(fēng)天氣引起的風(fēng)暴增水和大浪過程使得這兩段時(shí)間成為全年中沙灘侵蝕的高危時(shí)段。

    2.3 廣西潿洲島年極值高潮位推算分析

    本文統(tǒng)計(jì)了潿洲站1973-2014年共42年潮位年極值,其中有18年超過當(dāng)?shù)鼐涑蔽唬?2年超過當(dāng)?shù)乩碚摯蟪本€(見圖3)。超過當(dāng)?shù)乩碚摯蟪本€,即意味著潮水漫過潮間帶,此時(shí)極易打破沙灘泥沙的動態(tài)平衡。從圖3可以看出,2000年以前水位年極值超過當(dāng)?shù)卮蟪本€僅3年,近十幾年竟出現(xiàn)了9次。結(jié)合藍(lán)色趨勢線可看出近幾年來高潮位出現(xiàn)的頻次和極值均越來越高。在相同的波浪參數(shù)下,水層的厚度與海底摩擦消能成反比。高潮位使得近海區(qū)水層增厚,向岸傳遞波浪的能量就更多。這應(yīng)該是造成近幾年研究區(qū)域沙灘侵蝕愈加嚴(yán)重的重要原因之一。

    圖4展示了1973-2014年中,高潮位年極值出現(xiàn)的月份分布。從圖中可看出,年極值高潮位主要集中出現(xiàn)在11月、12月、1月和7月,占全年的80.0%,其中12月份為全年中最易出現(xiàn)年極值高潮位的月份(占26.7%)。從某種意義來說,雖然冬季大風(fēng)造成的風(fēng)暴增水和波浪遠(yuǎn)不如臺風(fēng)引起的劇烈,但是其對研究區(qū)域內(nèi)沙灘的侵蝕影響跟臺風(fēng)相比有過之而無不及。

    圖4 潿洲站1973-2014年高潮位年極值出現(xiàn)的月份分布Fig.4 Monthly probability of occurrence of the annual extreme high tide of Weizhou Island during 1973 to 2014

    3 廣西潿洲島多年一遇潮位及相應(yīng)淹沒范圍推算

    本文采用耿貝爾(Gumbel)曲線法[22—23]推算了潿洲站多年一遇最高潮位。Gumbe1分布曲線是以極值分布理論為基礎(chǔ),針對極值統(tǒng)計(jì)需要而提出的。對水文氣象要素多年一遇的極值推算有較大適用性,只要隨機(jī)變量的原始分布屬于正態(tài)分布或部分分布等指數(shù)分布族,其極值分布就漸近于Gumbel曲線。

    將潿洲站年最高潮位的經(jīng)驗(yàn)頻率點(diǎn)繪到圖上與理論曲線進(jìn)行對照,經(jīng)驗(yàn)頻率與理論曲線擬合較好(圖5)。

    圖5 潿洲站耿貝爾曲線與實(shí)測年極值經(jīng)驗(yàn)頻率分布Fig.5 Gumbel curve and empirical frequency in Weizhou tidal station

    表1給出了用耿貝爾方法計(jì)算所得潿洲站多年一遇高潮位。潿洲站建站至今最高潮位為517 cm,從表2中可知該高潮位為30年一遇的標(biāo)準(zhǔn)。50年一遇的高潮位為524.2 cm,100年一遇的高潮位為534.4 cm。需要說明的是,一定重現(xiàn)期(如50年)所對應(yīng)的高潮位(524.2 cm)并非是指今后這段時(shí)間(50年)一定會出現(xiàn),當(dāng)然也可能近幾年就會出現(xiàn)。從海洋防災(zāi)減災(zāi)的方面考慮,海洋工程可以依據(jù)推算結(jié)果做相應(yīng)重現(xiàn)期的防護(hù)。

    表1 潿洲站不同重現(xiàn)期年最高潮位Tab.1 Annual highest tidal level of Weizhou tidal station at different return periods

    圖6展示了50年和100年一遇高潮位在研究區(qū)域內(nèi)漫灘的情況??紤]到潿洲島西南部研究區(qū)域縱深較小,基本在50~150 m之間,潿洲島西南部漫灘范圍可以由所推算的高潮位與測點(diǎn)高程對比得出(測點(diǎn)最高分辨率達(dá)到0.5 m)。從圖中可看出,50年(100年)一遇高潮位在北部越過大潮線50~60(100) m,中部越過大潮線20~50(30~50) m,南部越過大潮線50~60(60~70) m。圖中標(biāo)注了歷史臺風(fēng)最大淹沒范圍以供參照。需要說明的是,歷史淹沒點(diǎn)為近期通過走訪當(dāng)?shù)鼐用窈筒榭催z留痕跡等方式獲取的數(shù)據(jù),并非某一次過程的淹沒范圍。另外,歷史淹沒點(diǎn)是高潮位和大浪聯(lián)合作用的結(jié)果,即如果多年一遇高潮疊加大浪,淹沒范圍肯定會比圖6所給范圍更加大。

    圖6 潿洲島西南部50年(a)和100年(b)一遇高潮位淹沒范圍示意圖Fig.6 Submerged areas of the southwestern Weizhou Island at 50-year (a) and 100-year (b) frequency tidal levels

    4 廣西潿洲島歷史風(fēng)暴潮的統(tǒng)計(jì)分析

    據(jù)統(tǒng)計(jì),1956-2014年中,廣西潿洲島沿海大于50 cm的臺風(fēng)風(fēng)暴潮過程共發(fā)生83次,其中4次超過當(dāng)?shù)鼐涑蔽?480 cm),分別是“8609”、“9204”、“0312”和“0814”號臺風(fēng)所引起。從圖7中可以看出,潿洲島的風(fēng)暴潮主要發(fā)生在7月份(占32.9%),其次是8月份(占27.6%),然后分別是10月(17.1%)和9月(15.8%)。

    圖7 1956-2014年潿洲島風(fēng)暴潮季節(jié)分布Fig.7 Seasonal distributions of the storm surge of Weizhou Island during 1956 to 2014

    臺風(fēng)風(fēng)暴潮過程中,增水大于80 cm,且最高潮位大于455 cm的典型風(fēng)暴潮過程共6次(表2)。從表2中也可以看出,影響潿洲島較強(qiáng)的臺風(fēng)風(fēng)暴潮過程主要集中在近20年內(nèi)。

    2003年第12號臺風(fēng)“科羅旺”(Krovanh)8月25日(農(nóng)歷七月廿八)4時(shí)前后登陸海南省文昌縣,登陸時(shí)臺風(fēng)最大風(fēng)速35 m/s,近中心最低氣壓970 hPa。受其影響廣東、廣西、海南沿海有12個(gè)站的最大增水超過1.0 m,有9個(gè)站的最高潮位超過當(dāng)?shù)鼐涑蔽?。廣東省南渡站增水最大,達(dá)3.59 m,廣西潿洲站最大增水達(dá)1.78 m。

    “科羅旺”影響期間,廣西沿海天文潮并不特別高,距離當(dāng)?shù)鼐涑蔽?0 cm左右。由于臺風(fēng)路徑和強(qiáng)度均特別有利,潿洲站于8月26日0時(shí)出現(xiàn)了達(dá)178 cm的最大風(fēng)暴增水,為該站自建站以來歷史上最大風(fēng)暴增水。同時(shí),大于1 m以上的風(fēng)暴增水持續(xù)21 h左右,疊加到高潮位上,最高潮位為517 cm,超過當(dāng)?shù)鼐涑蔽?7 cm,達(dá)到了橙色預(yù)警級別。該高潮位亦為歷史上最高潮位,經(jīng)過計(jì)算該高潮位達(dá)到了30年一遇的標(biāo)準(zhǔn)。

    表2 歷史上潿洲島顯著臺風(fēng)風(fēng)暴潮典型個(gè)例統(tǒng)計(jì)Tab.2 Statistics of the typical storm surge of Weizhou Island

    5 廣西潿洲島風(fēng)暴潮漫灘數(shù)值模擬

    針對于廣西沿海復(fù)雜的水深地形現(xiàn)狀,本文采用了目前被國際上廣泛接受的水動力模型——ADCIRC模式并耦合海浪模式(SWAN)對“0312”臺風(fēng)風(fēng)暴潮過程進(jìn)行模擬。ADCIRC模型是可以應(yīng)用于海洋、海岸、河口區(qū)域的多尺度水動力計(jì)算模型。模型基于非結(jié)構(gòu)網(wǎng)格,在空間上采用有限元方法離散,時(shí)間上采用有限差分法。該模型被美國工程兵部隊(duì)(ACE)和美國海軍研究實(shí)驗(yàn)室(NRL)廣泛應(yīng)用于各個(gè)軍港的潮汐、海流和風(fēng)暴潮預(yù)報(bào)中[24]。

    5.1 模式設(shè)置及數(shù)據(jù)來源

    模型的水深數(shù)據(jù)采用國家海洋環(huán)境預(yù)報(bào)中心業(yè)務(wù)化運(yùn)行的分辨率為2′×2′的南海水深數(shù)據(jù)插值得到。廣西沿海及其附近海域,尤其潿洲島沿海的水深數(shù)據(jù)采用最新加密測量的水深地形資料插值得到。

    為了減少或避免計(jì)算區(qū)域開邊界對模擬結(jié)果的影響,更好地模擬風(fēng)暴潮在北部灣海域的傳播過程,本文中所做網(wǎng)格的計(jì)算區(qū)域選在涵蓋了南海大部分范圍,只采用一條開邊界。邊界上添加NAO.99 8個(gè)主要分潮(M2,S2,K2,N2,K1,O1,P1,Q1)作為潮驅(qū)動。模式時(shí)間步長為10 s,模擬時(shí)段設(shè)置為5 d,涵蓋整個(gè)風(fēng)暴潮過程。模式進(jìn)行數(shù)值模擬時(shí),采用球坐標(biāo)、二維模式類型、混合底摩擦形式等,風(fēng)場驅(qū)動采用Holland風(fēng)場。

    SWAN模型中,采用二維動譜能密度表示隨機(jī)波浪場。模型考慮了風(fēng)能輸入,波浪的折射和反射,波浪破碎以及波-波相互作用等。海浪模式為風(fēng)暴潮模式提供輻射應(yīng)力,并作為外力驅(qū)動加入到風(fēng)應(yīng)力中,波長和周期信息也將用來參與風(fēng)暴潮模式的風(fēng)應(yīng)力計(jì)算。波浪譜方向的分辨率為10°,模式時(shí)間步長為1 200 s,每隔1 h與風(fēng)暴潮模式交換一次數(shù)據(jù)。

    圖8和圖9分別展示了潿洲島附近以及其西南岸灘的高分辨率非結(jié)構(gòu)網(wǎng)格分布情況。整個(gè)計(jì)算區(qū)域包括了949 728個(gè)三角形單元,共計(jì)483 043個(gè)節(jié)點(diǎn)。對廣西沿海區(qū)域,尤其是潿洲島海域進(jìn)行了重點(diǎn)加密,廣西沿海分辨率達(dá)到200 m左右,潿洲島域附近海域分辨率達(dá)到100 m左右,能夠較好地刻畫重點(diǎn)地物。

    圖8 潿洲島附近精細(xì)化網(wǎng)格剖分(自北向南分別是潿洲島和斜陽島)Fig.8 Fine-resolution model gird around Weizhou Island (from north to south: Weizhou Island and Xieyang Island)

    圖9 潿洲島西南部沙灘位置及網(wǎng)格剖分(方框內(nèi)為西南部沙灘)Fig.9 The location of the beach in the southwestern Weizhou Island (the domain in the rectangular box) and the surrounding model gird

    圖10 模式模擬天文潮與調(diào)和分析天文潮對比圖Fig.10 Comparison of astronomical tides from the model result and harmonic analysis

    5.2 模式驗(yàn)證

    圖10為2003年第12號臺風(fēng)“科羅旺”影響期間,風(fēng)暴潮-海浪耦合模式中廣東、廣西和海南各代表潮位站天文潮模擬結(jié)果與調(diào)和分析結(jié)果比較驗(yàn)證圖。從圖中可看出,各潮位站天文潮的振幅和遲角與調(diào)和分析結(jié)果符合較好。北部灣沿海因?yàn)槌毕匦员容^復(fù)雜,天文潮模擬效果較廣東、海南沿海稍差一些,但主要的潮汐特征也較好地模擬出來。

    圖11 模式模擬風(fēng)暴增水與實(shí)況增水對比圖Fig.11 Comparison of storm surge from the model result and the measured data

    圖11為“科羅旺”影響期間,廣西和海南各代表潮位站模擬增水與實(shí)況增水比較驗(yàn)證圖。從圖中可看出,各站的增水值和最大增水相位與實(shí)況符合較好。需要說明的是,模式模擬的風(fēng)暴增水由耦合模式加風(fēng)場驅(qū)動的結(jié)果減去不加風(fēng)場(即天文潮)結(jié)果所得。

    5.3 風(fēng)暴潮漫灘結(jié)果

    風(fēng)暴潮漫灘結(jié)果同樣是通過密集的陸地高程測點(diǎn)(最小間距0.5 m)是否淹沒來展示的。測點(diǎn)是否被淹沒則依據(jù)調(diào)和分析天文潮疊加模擬所得風(fēng)暴增水與陸地測點(diǎn)高程對比所得。研究區(qū)域內(nèi)沙灘縱深較小(50~150 m之間),在動態(tài)網(wǎng)格較難刻畫的情況下,這種判別方法是比較經(jīng)濟(jì)且合理可信的。

    從前人的研究可以看出,較嚴(yán)重的岸灘侵蝕事件均由高潮位配合大浪造成。因?yàn)轱L(fēng)暴增水或天文大潮引起的高潮位,使得近海區(qū)水層增厚,波浪能增加,同時(shí)海底摩擦消能作用也相應(yīng)減小,“浪借潮勢”,海浪向岸傳遞的能量就更多,對海岸侵蝕作用也會越嚴(yán)重。除了水層增厚帶來的影響,高潮位持續(xù)時(shí)間長短也決定著侵蝕的嚴(yán)重程度[10]。蔡鋒等[9]分析比較了臺風(fēng)“艾利”對福建沿海9個(gè)岸灘剖面的侵蝕作用,指出沙灘侵蝕區(qū)域基本上與大波高所處位置一致,而大波高位置與岸灘坡度和高潮位置緊密相關(guān)。根據(jù)測量計(jì)算,潿洲島西南部沙灘灘面坡度約為1∶30,與崇武青山剖面比較類似,大波高作用的位置也是灘面沙體主要損失位置應(yīng)在高潮帶。

    從圖12中看出,2003年第12號臺風(fēng)“科羅旺”影響期間,長時(shí)間的風(fēng)暴增水疊加到高潮上,使得潮水的大范圍地突破潮間帶,最遠(yuǎn)漫過大潮線近100 m;沙灘北部、中部及南部均有出現(xiàn)漫過大潮線甚至漫過海邊小路的情況。這次風(fēng)暴潮過程淹沒范圍和歷史淹沒點(diǎn)大致相近,中南部沙灘甚至超過歷史淹沒點(diǎn)。考慮海浪和風(fēng)暴潮或大潮的聯(lián)合作用,此時(shí)勢必打破沿海沙灘泥沙的動態(tài)平衡,造成沙灘北部、中部及南部漫灘地區(qū)泥沙損失,進(jìn)而形成對海岸的侵蝕。

    參照圖6和圖12中潮位漫灘的位置與范圍,可以有針對性地栽種防風(fēng)林樹種和固沙植被,并適當(dāng)采取防護(hù)性工程建設(shè),如人工補(bǔ)沙、修建岸邊道路護(hù)堤等,以阻止沙灘繼續(xù)后退。

    圖12 0312號臺風(fēng)影響期間最大淹沒區(qū)域分布Fig.12 Maximum submerged area during the process of the 0312 typhoon

    6 結(jié)論

    本文通過對廣西潿洲島氣候、潮汐和風(fēng)暴潮等歷史資料統(tǒng)計(jì)分析并結(jié)合風(fēng)暴潮漫灘數(shù)值模擬來分析風(fēng)暴潮和大潮對潿洲島西南部沙灘侵蝕的影響,得出如下結(jié)論:

    (1)根據(jù)風(fēng)暴潮歷史資料統(tǒng)計(jì),廣西潿洲島的臺風(fēng)風(fēng)暴潮主要發(fā)生在7-8月(占60.5%),其中7月份發(fā)生頻次最多(占32.9%)。同時(shí),每年5-7月和10-12月為潿洲島的天文大潮期,這兩個(gè)時(shí)間段恰好分別是臺風(fēng)影響時(shí)段和大風(fēng)天數(shù)最多時(shí)段,因此大風(fēng)天氣引起的風(fēng)暴增水和大浪配合天文高潮位使得這兩段時(shí)間成為全年中沙灘侵蝕的高危時(shí)段。

    (2)近幾年來高潮位出現(xiàn)的頻次和極值均越來越高。年極值高潮位主要集中出現(xiàn)在11月、12月、1月和7月,其中12月份為全年中最易出現(xiàn)年極值高潮位的月份。從某種意義來說,雖然冬季大風(fēng)造成的風(fēng)暴增水和波浪遠(yuǎn)不如臺風(fēng)引起的劇烈,但是其對潿洲島西南部沙灘的侵蝕影響跟臺風(fēng)相比有過之而無不及。

    (3)經(jīng)耿貝爾方法推算,廣西潿洲島50(100)年一遇的水位為524.2(534.4) cm。50年(100年)一遇高水位在研究區(qū)域北部越過大潮線50~60(100) m,中部越過大潮線20~50(30~50) m,南部越過大潮線50~60(60~70) m。海洋工程可以依據(jù)推算結(jié)果做相應(yīng)重現(xiàn)期的防護(hù)。

    (4)風(fēng)暴潮和大潮對沙質(zhì)海岸侵蝕的影響不容忽視。通過歷史資料統(tǒng)計(jì)以及風(fēng)暴潮-海浪耦合模式結(jié)果的對比發(fā)現(xiàn),較強(qiáng)的風(fēng)暴潮過程疊加在天文潮高潮以后,潮水就會漫過潿洲島西南部大潮線,在大浪的配合作用下加劇對海灘的侵蝕作用??梢栽谝浊治g區(qū)域有針對性地栽種防風(fēng)林樹種和固沙植被,并適當(dāng)采取防護(hù)性工程建設(shè)以阻止沙灘繼續(xù)后退。

    致謝:本文中使用的廣西潿洲島附近水深地形數(shù)據(jù)、陸地高程數(shù)據(jù)以及部分潮位資料由國家海洋局北海海洋環(huán)境監(jiān)測中心站測量并提供,在此表示感謝!

    [1] 姚子恒,高偉,高珊,等. 廣西北海潿洲島海岸侵蝕特征[J] . 海岸工程,2013,32(4):31-40.

    Yao Ziheng,Gao Wei,Gao Shan,et al. Coastal erosion of the Weizhou Island in Beihai,Guangxi Province[J]. Journal of Coastal Research,2013,32(4): 31-40.

    [2] 張忍順,陸麗云,王艷紅. 江蘇海岸侵蝕過程及其趨勢[J]. 地理研究,2002,21(4):469-478.

    Zhang Renshun,Lu Liyun,Wang Yanhong. The mechanism and trend of coastal erosion of Jiangsu Province in China[J]. Geographical Research,2002,21(4): 469-478.

    [3] 蔡鋒,蘇賢澤,劉建輝,等. 全球氣候變化背景下我國海岸侵蝕問題及防范對策[J]. 自然科學(xué)進(jìn)展,2008,18(10):1093-1103.

    Cai Feng,Su Xianze,Liu Jianhui,et al. China coastal erosion in global climate change and its countermeasures[J]. Progress in Natural Science,2008,18(10): 1093-1103.

    [4] 崔猛,牛茜茹,張緒良. 青島市海岸侵蝕的原因與防治[J]. 中國農(nóng)學(xué)通報(bào),2012,28(05):283-288.

    Cui Meng,Niu Xiru,Zhang Xuliang. The origin and countermeasures of coastal erosion of Qingdao City[J]. Chinese Agricultural Science Bulletin,2012,28(05): 283-288.

    [5] Liu J,Cai F,Qi H,et al. Coastal erosion along the west coast of the Taiwan Strait and its influencing factors[J]. Journal of Ocean University of China,2011,10(1): 23-34.

    [6] 王文海,吳桑云,陳雪英. 山東省9216號強(qiáng)熱帶氣旋風(fēng)暴期間的海岸侵蝕災(zāi)害[J]. 海洋地質(zhì)與第四紀(jì)地質(zhì),1996,14(4):71-78.

    Wang Wenhai,Wu Sangyun,Chen Xueying. Erosional disasters caused by storm surge during No. 9216 strong tropical cyclone along Shandong coast[J]. Marine Geology & Quaternary Geology,1996,14(4): 71-78.

    [7] 吳桑云,王文海,武桂秋. 山東沿海的9216號臺風(fēng)暴潮災(zāi)害[J]. 災(zāi)害學(xué),1996,9(1):44-47.

    Wu Sangyun,Wang Wenhai,Wu Guiqiu. Storm surge disasters caused by No. 9216 typhoon along Shandong coast[J]. Journal of Catastrophology,1996,9(1): 44-47.

    [8] 蔡鋒,蘇賢澤,楊順良,等. 廈門島海灘剖面對9914號臺風(fēng)大浪波動力的快速響應(yīng)[J]. 海洋工程,2002,20(2):85-90.

    Cai Feng,Su Xianze,Yang Shunliang,et al. A rapid response to 9914 typhoon-induced storm wave force made by the beach profiles of Xiamen Island[J]. The Ocean Engineering,2002,20(2): 85-90.

    [9] 蔡鋒,雷剛,蘇賢澤,等. 臺風(fēng)“艾利”對福建沙質(zhì)海灘影響過程研究[J]. 海洋工程,2006,24(1):98-109.

    Cai Feng,Lei Gang,Su Xianze,et al. Study on process response of Fujian beach geomorphology to typhoon Aere[J]. The Ocean Engineering,2006,24(1): 98-109.

    [10] Morton R A,Sallenger Jr A H. Morphological impacts of extreme storms on sandy beaches and barriers[J]. Journal of Coastal Research,2003,19(3): 560-573.

    [11] Cazenave A,Cozannet G L. Sea level rise and its coastal impacts[J]. Earth’s Future,2014,2(2): 15-34.

    [12] Nicholls R J,Cazenave A. Sea-level rise and its impact on coastal zones[J]. Science,2010,328(5985): 1517-1520.

    [13] 豐愛平,夏東興,谷東起,等. 萊州灣南岸海岸侵蝕過程與原因研究[J]. 海洋科學(xué)進(jìn)展,2006,24(1):83-90.

    Feng Aiping,Xia Dongxing,Gu Dongqi,et al. Study on process and cause of the coastal erosion along the south coast of the Laizhou Bay[J]. Advances in Marine Science,2006,24(1): 83-90.

    [14] Hsu T W,Lin T Y,Tseng I F. Human impact on coastal erosion in Taiwan[J]. Journal of Coastal Research,2007,23(4): 961-973.

    [15] 季子修. 中國海岸侵蝕特點(diǎn)及侵蝕加劇原因分析[J]. 自然災(zāi)害學(xué)報(bào),1996,5(2):65.

    Ji Zixiu. The characteristics of coastal erosion and cause of erosion[J]. Journal of Natural Disasters,1996,5(2): 65.

    [16] 石謙. 風(fēng)暴潮疊加天文大潮下的沙灘循環(huán)與海岸侵蝕[J]. 廈門理工學(xué)院學(xué)報(bào),2010,18(4):51-55.

    Shi Qian. Beach cycle and coastal erosion under the conditions of a storm surge combined with a spring tide[J]. Journal of Xiamen University of Technology,2010,18(4): 51-55.

    [17] Otvos E G. Beach aggradation following hurricane landfall: impact comparisons from two contrasting hurricanes,northern Gulf of Mexico[J]. Journal of Coastal Research,2004,20(1): 326-339.

    [18] 蔡鋒,蘇賢澤,夏東興. 熱帶氣旋前進(jìn)方向兩側(cè)海灘風(fēng)暴效應(yīng)差異研究——以海灘對0307號臺風(fēng)“伊布都”的響應(yīng)為例[J]. 海洋科學(xué)進(jìn)展,2004,22(4): 436-445.

    Cai Feng,Su Xianze,Xia Dongxing. Study on the difference between storm effects of beaches on two sides of the tropical cyclone track: taking the responses of Beaches to NO. 0307 typhoon Imbudo as an example[J]. Advances in Marine Science,2004,22(4): 436-445.

    [19] Morio O,Garlan T,Suanez S,et al. Modelling impacts of “Johanna storm” on an open-beach with ECORS simulators[C]//Coastal Dynamics 2013-7th International Conference on Coastal Dynamics. Arcachon,F(xiàn)rance,2013:1251-1262.

    [20] 劉仲桂. 廣西沿海風(fēng)能概況[J]. 武漢水利電力學(xué)院學(xué)報(bào),1989,22(6):130-134.

    Liu Zhongui. General situation of the wind energy along the coast of Guangxi[J]. Engineering Journal of Wuhan University,1989,22(6): 130-134.

    [21] 王慧,隋偉輝. 中國近海18個(gè)海區(qū)的海面大風(fēng)季節(jié)變化特征[J]. 天氣預(yù)報(bào)技術(shù)總結(jié)???,2012,4(3):61-65.

    Wang Hui,Sui Weihui. The seasonal variation of sea surface wind in 18 sea areas of China offshore[J]. Weather Forecast Review,2012,4(3): 61-65.

    [22] 候瑞科. 利用耿貝爾極值分布計(jì)算年最高水位[J]. 海洋通報(bào),1993,12(1): 126-129.

    Hou Ruike. Computation of annual highest water levels using Gumbel curve method[J]. Marine Science Bulletin,1993,12(1): 126-129.

    [23] 李明杰,齊鵬,侯一筠. 山東沿岸多年一遇最高水位計(jì)算[J]. 海洋科學(xué),2009,33(11): 78-81.

    Li Mingjie,Qi Peng,Hou Yijun. Computation of highest water levels of multiyear return periods along the coast of Shandong[J]. Marine Sciences,2009,33(11): 78-81.

    [24] Luettich R A Jr,Westerink J J,Scheffner N W. ADCIRC: an advanced three-dimensional circulation model for shelves coasts and estuaries,report 1: theory and methodology of ADCIRC-2DDI and ADCIRC-3DL[R]. 1992.

    Impacts of storm surge and spring tide on the beach erosion ofthe southwestern Weizhou Island,Guangxi Province

    Li Mingjie1,Wu Shaohua1,Liu Qiuxing1,Dong Jianxi1

    (1.NationalMarineEnvironmentalForecastingCenter,Beijing100081,China)

    In this paper,various history datum,such as meteorological data,tidal data,and storm surge data,etc,are collected and analyzed. Highest tidal levels of multiyear return periods are calculated using the Gumbel curve method,and the related inundation areas are also showed. It is revealed that the worsening beach erosion of the southwestern Weizhou Island is mainly due to the frequently occurrence of the extreme high tide in recent years. In addition,a coupled storm surge-wave numerical model is used to present the inundation process by storm surge during the 0312 typhoon. Impacts of storm surge and spring tide on the beach erosion is analyzed and elaborated,which is of certain directive significance to the local beach restoration and protection.

    storm surge; spring tide; beach erosion; Gumbel curve method; numerical simulation

    2014-10-17;

    2015-01-20。

    海洋公益性行業(yè)科研專項(xiàng)經(jīng)費(fèi)項(xiàng)目——沿海重點(diǎn)保障區(qū)域精細(xì)化綜合預(yù)報(bào)系統(tǒng)研制與應(yīng)用(201305031);廣西北海潿洲島整治修復(fù)項(xiàng)目——北海市潿洲島西南海域岸灘侵蝕研究。

    李明杰(1984—),男,江蘇省泰興市人,博士,助理研究員,主要從事風(fēng)暴潮預(yù)報(bào)及相關(guān)研究。E-mail:limj@nmefc.gov.cn

    10.3969/j.issn.0253-4193.2015.09.013

    P731.23

    A

    0253-4193(2015)09-0126-12

    李明杰,吳少華,劉秋興,等. 風(fēng)暴潮、大潮對廣西潿洲島西南沙灘侵蝕的影響分析[J]. 海洋學(xué)報(bào),2015,37(9):126-137,

    Li Mingjie,Wu Shaohua,Liu Qiuxing,et al. Impacts of storm surge and spring tide on the beach erosion of the southwestern Weizhou Island,Guangxi Province[J]. Haiyang Xuebao,2015,37(9):126-137,doi:10.3969/j.issn.0253-4193.2015.09.013

    猜你喜歡
    風(fēng)暴潮大潮潮位
    第五課:中國青年投身改革開放的大潮
    基于距離倒數(shù)加權(quán)的多站潮位改正方法可行性分析
    2012年“蘇拉”和“達(dá)維”雙臺風(fēng)影響的近海風(fēng)暴潮過程
    唐山市警戒潮位標(biāo)志物維護(hù)研究
    防范未來風(fēng)暴潮災(zāi)害的綠色海堤藍(lán)圖
    科學(xué)(2020年4期)2020-11-26 08:27:00
    基于多變量LSTM神經(jīng)網(wǎng)絡(luò)模型的風(fēng)暴潮臨近預(yù)報(bào)
    5G手機(jī)大潮即來
    海峽姐妹(2019年1期)2019-03-23 02:42:54
    多潮位站海道地形測量潮位控制方法研究
    基于改進(jìn)的OLS-RBF模型的感潮河段潮位預(yù)測研究
    《神威啟示錄》系列報(bào)道四 大潮退去看“裸泳”
    天堂影院成人在线观看| 在线观看免费视频日本深夜| 国产一区二区三区视频了| 色av中文字幕| 亚洲av电影不卡..在线观看| 免费观看的影片在线观看| 一区二区三区激情视频| 国产淫片久久久久久久久 | 国产主播在线观看一区二区| 3wmmmm亚洲av在线观看| 男插女下体视频免费在线播放| 男插女下体视频免费在线播放| 亚洲五月婷婷丁香| 亚洲精品久久国产高清桃花| 最近最新免费中文字幕在线| 久久久成人免费电影| 亚洲电影在线观看av| 午夜日韩欧美国产| 日韩大尺度精品在线看网址| 两人在一起打扑克的视频| 色播亚洲综合网| 男人舔奶头视频| 亚洲片人在线观看| 精品久久久久久久久亚洲 | 婷婷亚洲欧美| 精品国内亚洲2022精品成人| 性插视频无遮挡在线免费观看| 五月玫瑰六月丁香| 一区福利在线观看| 美女被艹到高潮喷水动态| 欧美中文日本在线观看视频| 一进一出抽搐动态| 国产不卡一卡二| 成人特级黄色片久久久久久久| 在线观看舔阴道视频| 又爽又黄无遮挡网站| 久久性视频一级片| www.色视频.com| 变态另类成人亚洲欧美熟女| 国产成人啪精品午夜网站| 禁无遮挡网站| 免费高清视频大片| 日韩欧美国产一区二区入口| 真人做人爱边吃奶动态| 免费在线观看亚洲国产| av视频在线观看入口| 精品久久久久久久久亚洲 | 国产探花在线观看一区二区| 男人和女人高潮做爰伦理| 国产精品爽爽va在线观看网站| 国产主播在线观看一区二区| 国内揄拍国产精品人妻在线| 三级男女做爰猛烈吃奶摸视频| 天堂影院成人在线观看| 最后的刺客免费高清国语| 国产91精品成人一区二区三区| 午夜两性在线视频| 国内精品一区二区在线观看| 午夜日韩欧美国产| 小蜜桃在线观看免费完整版高清| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲第一电影网av| 午夜激情福利司机影院| 亚洲精华国产精华精| 女人十人毛片免费观看3o分钟| 久久中文看片网| av中文乱码字幕在线| 久久精品人妻少妇| 一级黄色大片毛片| 亚洲av中文字字幕乱码综合| 真实男女啪啪啪动态图| 午夜久久久久精精品| 波多野结衣巨乳人妻| 黄色丝袜av网址大全| 亚洲第一欧美日韩一区二区三区| 性欧美人与动物交配| 一级作爱视频免费观看| 色精品久久人妻99蜜桃| 俄罗斯特黄特色一大片| 成人特级av手机在线观看| 深夜精品福利| 欧美在线黄色| 亚洲 国产 在线| 成人精品一区二区免费| 国内精品久久久久精免费| 日日摸夜夜添夜夜添小说| 亚洲av中文字字幕乱码综合| 精品日产1卡2卡| 两人在一起打扑克的视频| 色av中文字幕| 久久草成人影院| 别揉我奶头~嗯~啊~动态视频| 97热精品久久久久久| 90打野战视频偷拍视频| 91麻豆精品激情在线观看国产| 亚洲av免费在线观看| 欧美+亚洲+日韩+国产| 免费看光身美女| 久久久久精品国产欧美久久久| 精品日产1卡2卡| 首页视频小说图片口味搜索| 99热6这里只有精品| 日本黄色片子视频| 日本三级黄在线观看| 精品久久久久久久人妻蜜臀av| 国产探花在线观看一区二区| 人人妻人人看人人澡| 国产精品国产高清国产av| 天天一区二区日本电影三级| 国产免费av片在线观看野外av| 一进一出好大好爽视频| 美女黄网站色视频| 99热只有精品国产| 深爱激情五月婷婷| 亚洲成人久久爱视频| 久久伊人香网站| 淫妇啪啪啪对白视频| 特级一级黄色大片| 国产白丝娇喘喷水9色精品| 亚洲av五月六月丁香网| 成人高潮视频无遮挡免费网站| 亚洲 欧美 日韩 在线 免费| 中文亚洲av片在线观看爽| 国产高清视频在线播放一区| 永久网站在线| av天堂在线播放| 国产三级中文精品| 男人舔奶头视频| 国产伦精品一区二区三区视频9| 噜噜噜噜噜久久久久久91| 亚洲欧美日韩卡通动漫| 日韩亚洲欧美综合| 特大巨黑吊av在线直播| a级一级毛片免费在线观看| 在线观看66精品国产| 亚洲专区中文字幕在线| 国产野战对白在线观看| 在线免费观看的www视频| 最近视频中文字幕2019在线8| 真人一进一出gif抽搐免费| 91麻豆av在线| 久久精品国产亚洲av涩爱 | 国产精品99久久久久久久久| 国产一区二区三区视频了| 草草在线视频免费看| 舔av片在线| 国产精品嫩草影院av在线观看 | 日韩av在线大香蕉| 午夜激情欧美在线| 国产精品一区二区免费欧美| .国产精品久久| 欧美+亚洲+日韩+国产| avwww免费| 日本精品一区二区三区蜜桃| 亚洲七黄色美女视频| 国产精品美女特级片免费视频播放器| 黄色丝袜av网址大全| 日日摸夜夜添夜夜添av毛片 | 亚洲第一电影网av| netflix在线观看网站| 欧美不卡视频在线免费观看| 白带黄色成豆腐渣| 免费无遮挡裸体视频| 网址你懂的国产日韩在线| 日本成人三级电影网站| 小说图片视频综合网站| 99视频精品全部免费 在线| 精品国产亚洲在线| 成人av一区二区三区在线看| 日韩成人在线观看一区二区三区| 九色成人免费人妻av| 成人特级av手机在线观看| 99国产综合亚洲精品| 亚洲aⅴ乱码一区二区在线播放| 亚洲精品久久国产高清桃花| 可以在线观看的亚洲视频| 在线免费观看不下载黄p国产 | 一个人观看的视频www高清免费观看| 最后的刺客免费高清国语| 国产麻豆成人av免费视频| 麻豆av噜噜一区二区三区| 国产视频内射| 国产成人欧美在线观看| 如何舔出高潮| 尤物成人国产欧美一区二区三区| 高清在线国产一区| 久久人妻av系列| 国内精品久久久久精免费| 韩国av一区二区三区四区| 在线十欧美十亚洲十日本专区| 免费av不卡在线播放| 精品一区二区三区人妻视频| 久久久久亚洲av毛片大全| 久久婷婷人人爽人人干人人爱| 国产伦精品一区二区三区视频9| 亚洲18禁久久av| 久久久久久国产a免费观看| 国产精品久久久久久人妻精品电影| 男人的好看免费观看在线视频| 99精品久久久久人妻精品| 别揉我奶头~嗯~啊~动态视频| 美女高潮的动态| 精品久久久久久成人av| 真实男女啪啪啪动态图| 99久久99久久久精品蜜桃| 国产精品1区2区在线观看.| 欧美三级亚洲精品| 波多野结衣高清无吗| 欧美bdsm另类| 久久精品综合一区二区三区| 俄罗斯特黄特色一大片| 久久6这里有精品| 午夜精品在线福利| 又粗又爽又猛毛片免费看| 自拍偷自拍亚洲精品老妇| 久久国产乱子免费精品| 久久99热这里只有精品18| 亚洲av成人av| 色综合站精品国产| 99riav亚洲国产免费| 毛片女人毛片| 亚洲久久久久久中文字幕| 搡老熟女国产l中国老女人| 欧美激情国产日韩精品一区| 久久精品国产亚洲av涩爱 | 欧美不卡视频在线免费观看| 动漫黄色视频在线观看| 亚洲精品一区av在线观看| 97超视频在线观看视频| 午夜福利免费观看在线| 真人一进一出gif抽搐免费| 校园春色视频在线观看| 免费人成在线观看视频色| 久久国产精品影院| 久久热精品热| 国产主播在线观看一区二区| or卡值多少钱| 欧美成人性av电影在线观看| 可以在线观看的亚洲视频| 99热这里只有是精品在线观看 | 国产真实伦视频高清在线观看 | 草草在线视频免费看| 国产一区二区激情短视频| 国产黄色小视频在线观看| 禁无遮挡网站| 日本黄色片子视频| 一本精品99久久精品77| 国产av麻豆久久久久久久| 国语自产精品视频在线第100页| 亚洲一区二区三区色噜噜| 自拍偷自拍亚洲精品老妇| 草草在线视频免费看| 哪里可以看免费的av片| 国内毛片毛片毛片毛片毛片| 午夜精品在线福利| 黄色丝袜av网址大全| 成人欧美大片| 在线十欧美十亚洲十日本专区| 身体一侧抽搐| 免费看a级黄色片| 美女黄网站色视频| 欧美日韩乱码在线| 欧美色欧美亚洲另类二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 在线十欧美十亚洲十日本专区| 不卡一级毛片| 无遮挡黄片免费观看| 变态另类丝袜制服| 午夜亚洲福利在线播放| 免费观看的影片在线观看| 蜜桃亚洲精品一区二区三区| 真实男女啪啪啪动态图| 国产午夜精品论理片| 色哟哟哟哟哟哟| 少妇丰满av| 亚洲性夜色夜夜综合| 青草久久国产| 精品午夜福利视频在线观看一区| 亚洲欧美日韩高清在线视频| 伊人久久精品亚洲午夜| 亚洲成人免费电影在线观看| 精品人妻视频免费看| 精品午夜福利在线看| 欧美成人性av电影在线观看| 可以在线观看的亚洲视频| 日本在线视频免费播放| 日日摸夜夜添夜夜添av毛片 | 757午夜福利合集在线观看| 99久久久亚洲精品蜜臀av| 一区福利在线观看| 舔av片在线| 久久久国产成人免费| www.熟女人妻精品国产| 久久精品国产亚洲av涩爱 | 神马国产精品三级电影在线观看| 久久精品国产亚洲av香蕉五月| 免费观看人在逋| 在线看三级毛片| 久久国产精品人妻蜜桃| av福利片在线观看| 精品久久久久久,| 最近中文字幕高清免费大全6 | 男女视频在线观看网站免费| 精品一区二区三区视频在线| www.熟女人妻精品国产| 内射极品少妇av片p| 神马国产精品三级电影在线观看| 亚洲18禁久久av| 亚洲五月天丁香| 欧美黄色淫秽网站| 精品国内亚洲2022精品成人| 久久国产精品影院| 亚洲美女视频黄频| 一个人免费在线观看电影| 三级毛片av免费| 午夜久久久久精精品| 老司机深夜福利视频在线观看| 日本免费一区二区三区高清不卡| 少妇的逼好多水| 日本 欧美在线| av在线蜜桃| 丰满的人妻完整版| 天天躁日日操中文字幕| 亚洲第一欧美日韩一区二区三区| 国产精品不卡视频一区二区 | 久久久色成人| 国产 一区 欧美 日韩| 久久国产精品影院| 久久精品久久久久久噜噜老黄 | 日日摸夜夜添夜夜添小说| 免费无遮挡裸体视频| 精品人妻视频免费看| 欧美性感艳星| 99在线视频只有这里精品首页| 国产精品永久免费网站| 女同久久另类99精品国产91| 日韩欧美精品v在线| 国产一区二区三区在线臀色熟女| 色视频www国产| 脱女人内裤的视频| 国产精品久久久久久久电影| 国产高清激情床上av| 一级毛片久久久久久久久女| 99在线视频只有这里精品首页| 精品人妻一区二区三区麻豆 | 97人妻精品一区二区三区麻豆| 国产一区二区三区在线臀色熟女| 国产成人影院久久av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲一区二区三区不卡视频| 不卡一级毛片| 黄片小视频在线播放| 嫩草影院入口| 色综合欧美亚洲国产小说| 成人精品一区二区免费| 人妻夜夜爽99麻豆av| 欧美黑人欧美精品刺激| www.www免费av| 精品久久久久久成人av| 日日摸夜夜添夜夜添av毛片 | 999久久久精品免费观看国产| 我的女老师完整版在线观看| 久久久久久久精品吃奶| 国产精品永久免费网站| 老熟妇乱子伦视频在线观看| 亚洲成人中文字幕在线播放| 亚洲狠狠婷婷综合久久图片| 在线免费观看不下载黄p国产 | 国产久久久一区二区三区| 亚洲国产精品合色在线| 亚洲精品日韩av片在线观看| 国内毛片毛片毛片毛片毛片| 亚洲人成伊人成综合网2020| 国产亚洲精品综合一区在线观看| 国产爱豆传媒在线观看| 色综合站精品国产| 久久久久性生活片| 国产精品亚洲美女久久久| 亚洲成人精品中文字幕电影| 观看美女的网站| 国产单亲对白刺激| 久久久久免费精品人妻一区二区| 99久久精品热视频| 欧美xxxx性猛交bbbb| 精品久久久久久久人妻蜜臀av| 成年女人看的毛片在线观看| 国产亚洲精品久久久com| 欧美高清性xxxxhd video| 国产伦一二天堂av在线观看| 亚洲自拍偷在线| 亚洲av一区综合| 亚洲第一欧美日韩一区二区三区| 免费无遮挡裸体视频| 久久亚洲精品不卡| 最近最新免费中文字幕在线| 成人鲁丝片一二三区免费| 亚洲,欧美,日韩| 欧美在线黄色| 狠狠狠狠99中文字幕| 亚洲avbb在线观看| 少妇高潮的动态图| 欧美午夜高清在线| 男插女下体视频免费在线播放| 日本一本二区三区精品| 人妻夜夜爽99麻豆av| 深夜精品福利| 男人和女人高潮做爰伦理| 日韩欧美在线乱码| 人人妻人人澡欧美一区二区| 最近中文字幕高清免费大全6 | 淫秽高清视频在线观看| 桃色一区二区三区在线观看| 国产精品亚洲美女久久久| 亚洲成人中文字幕在线播放| 性插视频无遮挡在线免费观看| 尤物成人国产欧美一区二区三区| 免费av不卡在线播放| 午夜福利欧美成人| 日韩av在线大香蕉| 尤物成人国产欧美一区二区三区| 又紧又爽又黄一区二区| 免费av毛片视频| 青草久久国产| 成人一区二区视频在线观看| 成人国产一区最新在线观看| 18禁在线播放成人免费| 中文字幕免费在线视频6| 国产黄片美女视频| 欧美黑人欧美精品刺激| 精品熟女少妇八av免费久了| 亚洲成av人片在线播放无| 一级a爱片免费观看的视频| 亚洲五月天丁香| 中文字幕久久专区| 亚洲专区中文字幕在线| 欧美性猛交╳xxx乱大交人| 精品人妻一区二区三区麻豆 | 美女大奶头视频| 天天躁日日操中文字幕| 午夜激情欧美在线| 亚洲午夜理论影院| av欧美777| 亚洲精品久久国产高清桃花| 久久久久国内视频| 欧美色视频一区免费| 99久久99久久久精品蜜桃| 亚洲乱码一区二区免费版| 国产黄色小视频在线观看| 久久久久精品国产欧美久久久| 一区二区三区高清视频在线| 亚洲无线在线观看| 免费黄网站久久成人精品 | 毛片一级片免费看久久久久 | 欧美xxxx黑人xx丫x性爽| 高潮久久久久久久久久久不卡| 亚洲 国产 在线| 免费搜索国产男女视频| 亚洲av成人精品一区久久| 人人妻,人人澡人人爽秒播| 亚洲国产精品成人综合色| 夜夜躁狠狠躁天天躁| 精品人妻视频免费看| 欧美成人一区二区免费高清观看| 欧美一级a爱片免费观看看| 亚洲精品一区av在线观看| 夜夜爽天天搞| 啦啦啦观看免费观看视频高清| av天堂在线播放| 亚洲国产精品成人综合色| 国产 一区 欧美 日韩| 嫩草影视91久久| 亚洲七黄色美女视频| 国内久久婷婷六月综合欲色啪| 国产高清视频在线播放一区| 97热精品久久久久久| 亚洲欧美日韩高清在线视频| 国产三级在线视频| 性插视频无遮挡在线免费观看| 日本一本二区三区精品| 最好的美女福利视频网| 色精品久久人妻99蜜桃| 高潮久久久久久久久久久不卡| 欧美又色又爽又黄视频| 一个人观看的视频www高清免费观看| 脱女人内裤的视频| 欧美bdsm另类| xxxwww97欧美| 一个人看的www免费观看视频| 国产成年人精品一区二区| 色视频www国产| 国产视频一区二区在线看| 亚洲aⅴ乱码一区二区在线播放| 欧美潮喷喷水| 免费av毛片视频| 国产成人欧美在线观看| 色哟哟哟哟哟哟| 亚洲精品影视一区二区三区av| 国产激情偷乱视频一区二区| 男女做爰动态图高潮gif福利片| 51国产日韩欧美| 国产老妇女一区| 日本免费a在线| 免费观看精品视频网站| 男女下面进入的视频免费午夜| 97超级碰碰碰精品色视频在线观看| 黄色日韩在线| 天美传媒精品一区二区| a级一级毛片免费在线观看| 一进一出好大好爽视频| 首页视频小说图片口味搜索| 精品一区二区三区人妻视频| 91久久精品国产一区二区成人| 少妇的逼水好多| 少妇丰满av| 少妇人妻精品综合一区二区 | 国产探花在线观看一区二区| 亚洲精品成人久久久久久| 午夜福利视频1000在线观看| 好看av亚洲va欧美ⅴa在| 女生性感内裤真人,穿戴方法视频| 特大巨黑吊av在线直播| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区成人| 看片在线看免费视频| av国产免费在线观看| 精品人妻熟女av久视频| 国产一区二区三区视频了| 美女xxoo啪啪120秒动态图 | 国产高清激情床上av| 国产成人aa在线观看| 欧美三级亚洲精品| 毛片一级片免费看久久久久 | 欧美xxxx性猛交bbbb| 一区二区三区激情视频| 啪啪无遮挡十八禁网站| 亚洲成av人片在线播放无| 欧美性猛交黑人性爽| 亚洲av成人精品一区久久| 一进一出好大好爽视频| 成年女人毛片免费观看观看9| 午夜免费成人在线视频| 好男人在线观看高清免费视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩乱码在线| av专区在线播放| 床上黄色一级片| 日本一本二区三区精品| 免费高清视频大片| 国产伦精品一区二区三区四那| 久久中文看片网| 久久久久精品国产欧美久久久| 一级a爱片免费观看的视频| 免费在线观看成人毛片| 一区二区三区高清视频在线| 久久久久九九精品影院| 久久精品国产亚洲av涩爱 | 久久久精品欧美日韩精品| 99精品在免费线老司机午夜| 天堂√8在线中文| 真实男女啪啪啪动态图| 欧美bdsm另类| 免费看日本二区| 日本与韩国留学比较| 在线免费观看不下载黄p国产 | 神马国产精品三级电影在线观看| 青草久久国产| 国内精品美女久久久久久| 在线看三级毛片| 1000部很黄的大片| 久9热在线精品视频| 日韩中文字幕欧美一区二区| 午夜免费激情av| 国产精品亚洲一级av第二区| 69av精品久久久久久| 欧美+日韩+精品| 51午夜福利影视在线观看| 色5月婷婷丁香| 免费看美女性在线毛片视频| 国产伦精品一区二区三区四那| 欧美一区二区国产精品久久精品| 看免费av毛片| 色哟哟·www| 欧美高清成人免费视频www| 色哟哟·www| 国产激情偷乱视频一区二区| 床上黄色一级片| 成人一区二区视频在线观看| 亚洲av成人av| 国产探花在线观看一区二区| 亚洲熟妇中文字幕五十中出| 天堂网av新在线| 亚洲熟妇中文字幕五十中出| 精品一区二区三区av网在线观看| 在线免费观看的www视频| 国产aⅴ精品一区二区三区波| 日韩欧美一区二区三区在线观看| 国产精品永久免费网站| 99久久精品热视频| 别揉我奶头 嗯啊视频| 欧美午夜高清在线| 免费一级毛片在线播放高清视频| 欧美一区二区国产精品久久精品| 国产伦在线观看视频一区| 国产伦人伦偷精品视频| 久久久精品大字幕| 97人妻精品一区二区三区麻豆| 亚洲精品456在线播放app | 真人做人爱边吃奶动态| 色综合站精品国产| 小蜜桃在线观看免费完整版高清| 露出奶头的视频| 91狼人影院| 中亚洲国语对白在线视频| 欧美不卡视频在线免费观看| 嫩草影院入口|