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

    太湖不同湖區(qū)風(fēng)浪的季節(jié)變化特征*

    2016-02-23 10:21:07王震,吳挺峰,鄒華
    湖泊科學(xué) 2016年1期
    關(guān)鍵詞:水生植物空間分布風(fēng)浪

    ?

    太湖不同湖區(qū)風(fēng)浪的季節(jié)變化特征*

    (1:江南大學(xué)環(huán)境與土木工程學(xué)院,無(wú)錫 214122)

    (2:中國(guó)科學(xué)院南京地理與湖泊研究所湖泊與環(huán)境國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京 210008)

    (3:江蘇省水文水資源勘測(cè)局無(wú)錫分局,無(wú)錫 214031)

    (4:河海大學(xué)水文水資源學(xué)院,南京 210098)

    摘要:為明晰太湖風(fēng)浪的空間分布及季節(jié)變化,在湖心區(qū)設(shè)立波浪觀測(cè)站,利用其記錄的波浪數(shù)據(jù)證明SWAN模型能夠較好地模擬太湖風(fēng)浪.基于所建模型,對(duì)2013年自然風(fēng)場(chǎng)條件下太湖不同湖區(qū)風(fēng)浪季節(jié)動(dòng)態(tài)進(jìn)行模擬分析,結(jié)果表明:受岸線、地形和島嶼等地理因素影響,大太湖的風(fēng)浪總是最強(qiáng),其有效波高均值為0.523m;而東太湖風(fēng)浪最小,有效波高均值為0.305m.受盛行風(fēng)場(chǎng)季節(jié)變化影響,太湖春、夏季有效波高均值明顯大于秋、冬季.太湖波浪的能量主要來(lái)源于風(fēng)場(chǎng),其有效波高隨風(fēng)速增大而增大,兩者呈極顯著正相關(guān).而風(fēng)向則可以通過(guò)改變風(fēng)區(qū)長(zhǎng)度來(lái)影響風(fēng)浪生消.在偏東風(fēng)作用下,太湖湖西區(qū)的風(fēng)浪大于東部湖區(qū);而受盛行于冬季的偏北風(fēng)影響,太湖南部水域風(fēng)浪要大于北部.同時(shí),太湖風(fēng)浪的時(shí)空分布特征是造成太湖水質(zhì)參數(shù)、沉積物和水生植物空間分布差異的重要原因之一。

    關(guān)鍵詞:太湖分區(qū);風(fēng)浪;季節(jié)動(dòng)態(tài);空間分布;SWAN模型;水生植物

    Changes in seasonal characteristics of wind and wave in different regions of Lake Taihu

    在淺水湖泊中,水動(dòng)力擾動(dòng)引起的底泥懸浮不僅導(dǎo)致水體透明度下降,還可造成賦存在沉積物中的營(yíng)養(yǎng)鹽向上覆水層釋放,增加淺水湖泊富營(yíng)養(yǎng)化的風(fēng)險(xiǎn).平均水深 1.9m 的太湖是典型的淺水富營(yíng)養(yǎng)化湖泊[1].該湖泊水面開(kāi)闊,水底地形平坦,多數(shù)岸線暴露在長(zhǎng)風(fēng)區(qū)條件下,有利于風(fēng)浪的生成和發(fā)展.研究已經(jīng)證明,太湖底泥懸浮的主要能量來(lái)自波浪[2-3],湖流的作用幾乎可以忽略[4-5].因此,解析太湖波浪特征對(duì)研究其水動(dòng)力過(guò)程及其生態(tài)環(huán)境效應(yīng)具有重要意義。

    由于風(fēng)浪對(duì)太湖生態(tài)環(huán)境具有重要意義,近年來(lái)也有不少學(xué)者對(duì)風(fēng)浪開(kāi)展了研究.李一平等[6]在太湖實(shí)際波浪觀測(cè)的基礎(chǔ)上,對(duì)太湖波浪進(jìn)行數(shù)值模擬,研究了太湖波浪的主要影響因素.劉興平[7]研究了太湖波浪的生消動(dòng)態(tài)過(guò)程并結(jié)合湖流三維模型,尋求太湖湖流對(duì)波浪的影響.同時(shí),學(xué)者們也積極探討太湖波浪對(duì)生態(tài)環(huán)境的影響,如波浪和湖流共同作用下太湖水體懸浮物輸運(yùn)過(guò)程[8],太湖波浪與湖流對(duì)沉積物再懸浮的影響[9-11],太湖白帽覆蓋率的變化與特征[11].這些研究對(duì)探求太湖波浪數(shù)值模擬成果有一定意義,但都是在假想的定常風(fēng)或較短時(shí)間內(nèi)的實(shí)際風(fēng)場(chǎng)作用下對(duì)太湖波浪的數(shù)值模擬.而利用基于動(dòng)譜平衡方程開(kāi)發(fā)的SWAN模型,就實(shí)際風(fēng)場(chǎng)周年變化條件下太湖波浪時(shí)空分布的特征研究還鮮有報(bào)道。

    本文基于太湖實(shí)測(cè)風(fēng)場(chǎng)和波浪數(shù)據(jù),通過(guò)第3代淺水波浪模型SWAN,分析 2013 年自然風(fēng)場(chǎng)條件下太湖不同湖區(qū)風(fēng)浪時(shí)空分布特征,并從較大的時(shí)空尺度上討論風(fēng)浪的生態(tài)環(huán)境效應(yīng),本研究可為太湖生態(tài)環(huán)境問(wèn)題的動(dòng)力學(xué)成因研究提供幫助。

    1 數(shù)據(jù)與方法

    1.1 太湖概況

    太湖(30°55′40″~31°32′58″ N,119°52′32″~120°36′10″ E)是長(zhǎng)江三角洲的一個(gè)大型淺水湖泊(圖1).太湖水域面積 2338km2,南北長(zhǎng) 68.5km,東西平均寬 34km,最寬處 56km.太湖平均水深 1.9m,最大水深不超過(guò) 3m,其最深處位于湖心偏西、平臺(tái)山以北的區(qū)域.太湖湖面開(kāi)闊,水下地形平坦,平均坡度只有19.7″.按照岸線、水質(zhì)和水下地形等因素,太湖可分為大太湖、梅梁灣、貢湖、竺山湖、胥湖、西部沿岸、南部沿岸、東太湖和箭湖東茭咀 9 個(gè)湖區(qū)[12](圖1).太湖流域夏季受海陸氣溫差異影響,盛行東南風(fēng);冬季受來(lái)自西伯利亞的寒流影響,盛行西北風(fēng).由于全球氣候變暖、城鎮(zhèn)化等原因,太湖的年平均風(fēng)速和最大風(fēng)速均呈逐年下降的趨勢(shì)。

    圖1 觀測(cè)站及太湖分區(qū)示意Fig.1 Sketch of observation stations and nine regions in Lake Taihu

    1.2 第三代淺水波浪模型

    SWAN模型建立在波浪能量平衡方程基礎(chǔ)上,考慮了波浪作用過(guò)程中的波浪折射、反射、淺化、破碎、白帽、底摩擦及波-波相互作用,其中未考慮湖流的作用和波浪繞射的影響.SWAN模型能夠很好地運(yùn)用于海岸、湖泊、河口等淺水區(qū)域,為波浪的傳播以及帶來(lái)的影響提供充分的依據(jù)。

    1.2.1 動(dòng)譜平衡方程SWAN 模型不是以二維能譜密度而是以二維動(dòng)譜密度表示隨機(jī)波,因在流場(chǎng)中,動(dòng)譜密度守恒,而能譜密度不守恒,動(dòng)譜密度與能譜密度和相對(duì)頻率的關(guān)系為:N(σ,θ)=E(σ,θ)/σ,N(σ,θ)隨時(shí)間、空間而變化.在笛卡爾直角坐標(biāo)系下,動(dòng)譜平衡方程可表示為[13]:

    (1)

    式中,方程左邊第1項(xiàng)為N隨時(shí)間的變化率;第2、3項(xiàng)分別表示N在地理空間坐標(biāo)x、y方向上的傳播(Cx和Cy分別為x和y方向的傳播速率);第4項(xiàng)表示流場(chǎng)和水深所引起的N在相對(duì)頻率σ空間的變化;第5項(xiàng)為N在譜分布方向θ空間(譜方向分布范圍)的傳播;S為以譜密度表示的源匯項(xiàng),包括風(fēng)能輸入、波與波之間非線性相互作用和由于底摩擦、白浪、破碎等引起的能量損耗;Cx、Cy、Cσ和Cθ分別代表在x、y、σ和θ空間的波浪傳播速度。

    1.2.2 波浪能量輸入和損耗項(xiàng)風(fēng)能輸入使用共振機(jī)制[14]和反饋機(jī)制[15]來(lái)描述,相應(yīng)的源函數(shù)可表示為線性增長(zhǎng)和指數(shù)增長(zhǎng)兩部分之和.波浪能量的耗散項(xiàng)主要是白浪、底部摩檫和由水深引起的波破碎.白浪主要由波陡度控制,根據(jù)Hasselmann等提出的脈動(dòng)平均模型進(jìn)行計(jì)算[13];底摩擦采用Collins拖曳理論模型;水深引起的波破碎基于Miche準(zhǔn)則,Battjes等[16]根據(jù)大量試驗(yàn)數(shù)據(jù)及現(xiàn)場(chǎng)資料,認(rèn)為在淺水域,對(duì)于不同類型的地貌,隨機(jī)波的最大破碎波高Hm與水深d的關(guān)系可表示為:Hm=γd。

    1.2.3 波與波之間非線性相互作用在深水域,四相波與波非線性相互作用占主導(dǎo)地位,譜能由譜峰處向低頻轉(zhuǎn)移(使得峰頻變小)和高頻轉(zhuǎn)移(高頻處能量由于白浪而耗散掉).在淺水域,三相波與波之間非線性相互作用是主要影響因素,能量由低頻向高頻處轉(zhuǎn)移.在SWAN模型中,四相波-波相互作用采用Hasselmann等[13]提出的離散相互作用近似法(DIA)計(jì)算,在三相波相互作用的計(jì)算中,每個(gè)譜方向上均采用Eldeberky的集合三相近似模型(LTA)[17],它由Eldeberky等的離散三相近似模型[18]改進(jìn)而得。

    1.2.4 SWAN參數(shù)設(shè)置SWAN模型物理過(guò)程考慮底摩擦效應(yīng)、白浪損耗和非線性波-波之間相互作用等耗散機(jī)制.底摩擦引起的消耗采用Collins模型,底摩擦系數(shù)取0.025;破碎系數(shù)取均值0.73.沿逆時(shí)針?lè)较虻淖V分布方向θ空間(即譜方向分布范圍)的譜方向步長(zhǎng)取為10°,譜頻率σ計(jì)算范圍為0.04~1.00Hz,網(wǎng)格數(shù)為15.以JONSWAP譜和4階余弦函數(shù)表示二維譜.因入湖河道對(duì)太湖波浪的影響較小,可忽略不計(jì),作陸地邊界處理.以X10、Y10風(fēng)場(chǎng)作為模型風(fēng)輸入.計(jì)算網(wǎng)格為矩形網(wǎng)格,計(jì)算步長(zhǎng)1km,由此可將太湖劃分為69×69個(gè)網(wǎng)格,x、y方向上計(jì)算域長(zhǎng)度均為69km.同時(shí),模型驗(yàn)證的計(jì)算時(shí)間域?yàn)?014年8月14日6:10至8月17日5:40,計(jì)算時(shí)間步長(zhǎng)設(shè)置為10min,模型計(jì)算結(jié)果輸出的時(shí)間間隔為30min.模型模擬的計(jì)算時(shí)間域?yàn)?013年1月1日0:00至12月31日23:30,計(jì)算時(shí)間步長(zhǎng)設(shè)置為30min,模型計(jì)算結(jié)果輸出的時(shí)間間隔為 6h。

    1.3 數(shù)據(jù)采集與處理

    風(fēng)速、風(fēng)向觀測(cè)站位于太湖湖心(31°16′40″N,120°9′16″E),于2013年1月1日0:00至2013年12月31日23:30及2014年8月14日6:10至8月17日5:40逐時(shí)記錄水面以上10m高度處的10min平均風(fēng)速,采樣時(shí)間間隔30min,全湖采用同一風(fēng)場(chǎng)數(shù)據(jù)作為輸入,將此風(fēng)場(chǎng)分解至笛卡爾坐標(biāo)系的x軸和y軸,表現(xiàn)形式為X10、Y10:

    (2)

    (3)

    式中,WD代表站點(diǎn)處的風(fēng)向,WS代表站點(diǎn)處的風(fēng)速。

    波浪觀測(cè)使用英國(guó)Valeport有限公司生產(chǎn)的采樣頻率為8Hz的高精度Valeport,觀測(cè)站位于太湖湖心(31°18′33″N,120°12′22″E),于2014年8月14日6:10至8月17日05:40進(jìn)行觀測(cè),采樣時(shí)間間隔30min,得到有效波高(Hs)及譜峰周期(Tp)等。

    數(shù)據(jù)統(tǒng)計(jì)分析使用Excel和SPSS 20.0軟件,作圖使用Origin 8.1和Surfer 8軟件.采用Pearson相關(guān)分析法表示相關(guān)性,P<0.01表示兩者呈極顯著相關(guān),P<0.05表示兩者呈顯著相關(guān).使用確定性系數(shù)表示模擬值與實(shí)測(cè)值之間的擬合程度:

    (4)

    2 結(jié)果與討論

    2.1 模型驗(yàn)證

    圖2 模擬與實(shí)測(cè)有效波高的線性關(guān)系Fig.2 Linear relationship between calculated and measured values of significant wave height

    圖3 2014 年 8 月 14日 6:10 至 17 日 5:40 湖心觀測(cè)站有效波高實(shí)測(cè)值與計(jì)算值的比較Fig.3 Comparison of measured and calculated values of significant wave height in the center of Lake Taihu from August 14th 6:10 to 17th 5:40 of 2014

    通過(guò)建立2014年8月14日6:10至8月17日5:40間模擬與實(shí)測(cè)有效波高值之間的線性關(guān)系(圖2),表明模擬與實(shí)測(cè)有效波高值能夠緊密分布在直線y=0.946x兩側(cè),判定系數(shù)R2=0.960(樣本數(shù)n=288);相對(duì)誤差范圍為0.03%~28.7%,均值為10.7%;實(shí)測(cè)值與模擬值呈極顯著正相關(guān)(r=0.83,P<0.01),確定性系數(shù)為0.84.可見(jiàn)有效波高的模擬值能夠較好地反映太湖波浪的實(shí)際情況,SWAN模型能夠較好地表現(xiàn)太湖的風(fēng)浪特征。

    驗(yàn)證結(jié)果表明,模擬的有效波高與實(shí)測(cè)波高變化趨勢(shì)基本一致(圖3),但兩者之間存在一定的誤差.在模擬初期,實(shí)測(cè)值與模擬值誤差較大,隨著模擬時(shí)間的推移,風(fēng)浪場(chǎng)逐漸形成,模擬值與實(shí)測(cè)值越來(lái)越接近.誤差來(lái)源[19]可能是:(1) 由于全湖采用的是空間均勻的逐時(shí)風(fēng)場(chǎng),與實(shí)際風(fēng)場(chǎng)有一定的誤差;(2) 在阻擋物較多的狹長(zhǎng)水域以及湖灣內(nèi),以此空間均勻風(fēng)場(chǎng)帶來(lái)的誤差影響就會(huì)大些;(3) 未考慮波流相互作用;(4) 盡管風(fēng)浪是太湖波浪的主要能量來(lái)源,但是實(shí)測(cè)波浪還包括船行波和潮汐等。

    2.2 太湖不同湖區(qū)風(fēng)浪的季節(jié)分布特征

    圖4 太湖湖心氣象觀測(cè)站記錄的2013年各月月平均風(fēng)速變化過(guò)程Fig.4 Changes of monthly average wind speed in 2013 recorded by meteorological station in center of Lake Taihu

    2.2.1 2013年太湖風(fēng)場(chǎng)特征分析對(duì)太湖湖心觀測(cè)站(圖1)記錄的2013年風(fēng)速、風(fēng)向開(kāi)展統(tǒng)計(jì)分析.2013年太湖月均風(fēng)速(圖4)在3.0~6.0m/s之間,其中春季(3-5月)、夏季(6-8月)和秋季(9-11月)風(fēng)速相對(duì)較高,平均約4.9m/s;冬季(2013年1、2和12月)風(fēng)速低,為3.8m/s. 2013年太湖四季風(fēng)向分布情況見(jiàn)表1,春季和夏季以東南風(fēng)為主,秋季以東北風(fēng)為主,冬季以西北風(fēng)為主.依據(jù)蒲福風(fēng)級(jí)分類標(biāo)準(zhǔn),2013年太湖風(fēng)場(chǎng)序列統(tǒng)計(jì)情況見(jiàn)表2. 2013年太湖流域主要以輕風(fēng)、微風(fēng)、和風(fēng)為主,風(fēng)速介于1.6~7.9m/s.其中尤以微風(fēng)發(fā)生頻率最高,占全年的31.87%.全年觀測(cè)到的最大風(fēng)速為15.0m/s,達(dá)到疾風(fēng)標(biāo)準(zhǔn)。

    2.2.2 2013年太湖風(fēng)浪時(shí)空分布特征基于2013年實(shí)測(cè)風(fēng)場(chǎng),利用SWAN模型對(duì)太湖2013年1月1日0:00至12月31日23:30的風(fēng)浪進(jìn)行模擬,統(tǒng)計(jì)分析太湖風(fēng)浪的時(shí)空分布特征。

    2013年9個(gè)湖區(qū)12個(gè)月的有效波高的最大值見(jiàn)表3.2013年太湖的有效波高最大值在0.217~0.652m之間.大太湖的風(fēng)浪總是最強(qiáng),其有效波高均值為0.523m;而東太湖風(fēng)浪最小,有效波高均值為0.305m.有文獻(xiàn)報(bào)道[6,20],風(fēng)浪的大小與風(fēng)速、風(fēng)時(shí)、風(fēng)區(qū)長(zhǎng)度有著密切的關(guān)系;另外,風(fēng)浪的影響因子還包括水深、地形、岸線形狀、沉水植物覆蓋度等.大太湖平均水深2.3m,水面開(kāi)闊;東太湖平均水深僅為1.05m,水域呈狹長(zhǎng)型;這兩個(gè)湖區(qū)的特征有著本質(zhì)區(qū)別,東太湖的水力條件為沉水植物的生長(zhǎng)創(chuàng)造了有利的環(huán)境.事實(shí)上,通過(guò)多年的野外觀測(cè)發(fā)現(xiàn)大太湖水底幾乎無(wú)水草,而東太湖水草生長(zhǎng)茂盛.Jackson等[21]和Hamilton等[22]研究認(rèn)為水生植物能夠阻礙波浪的形成,并降低水流速度.因此,水生植物能阻礙波浪的發(fā)展,起到消浪的作用;進(jìn)而可以降低底泥再懸浮[23]和內(nèi)源釋放[24]的可能性.這又恰恰說(shuō)明草型湖泊(胥湖和東太湖)的透明度和水質(zhì)優(yōu)于藻型湖泊(梅梁灣、貢湖和竺山湖)[25-26].太湖北部3個(gè)湖灣有效波高值大小順序?yàn)椋好妨簽?貢湖>竺山湖.竺山湖水域面積最小,風(fēng)吹程短,不利于風(fēng)浪的生長(zhǎng),這是該湖灣有效波高小的主要原因之一.在盛行東南風(fēng)的夏季,西部沿岸區(qū)的有效波高值均比南部沿岸區(qū)高,這與西部沿岸區(qū)風(fēng)區(qū)長(zhǎng)度、岸線形狀規(guī)則有著密切的關(guān)系。

    表1 2013年太湖湖心氣象站記錄的風(fēng)向統(tǒng)計(jì)結(jié)果

    表2 2013 年太湖湖心氣象站記錄的風(fēng)速等級(jí)統(tǒng)計(jì)結(jié)果

    受控于風(fēng)速、風(fēng)向的季節(jié)波動(dòng),太湖春、夏季有效波高均值大于秋、冬季.SWAN模擬結(jié)果表明:2013年全太湖春、夏、秋和冬季有效波高平均值分別為0.433、0.444、0.418和0.343m.各湖區(qū)有效波高季節(jié)分布與此相似,例如:大太湖春、夏、秋、冬季有效波高分別為0.565、0.543、0.546和0.435m;東太湖春、夏、秋、冬季有效波高分別為:0.424、0.419、0.416和0.329m.風(fēng)速是導(dǎo)致太湖風(fēng)浪這種季節(jié)變化的主要因素.統(tǒng)計(jì)分析結(jié)果顯示,太湖有效波高與風(fēng)速呈極顯著正相關(guān)(r=0.47,P<0.01),即有效波高與風(fēng)速的季節(jié)變化一致,風(fēng)速越大,有效波高越大.此外,湖流[27]、水位[28]、行船[29]及水生植物等季節(jié)變化也可能對(duì)風(fēng)浪季節(jié)動(dòng)態(tài)產(chǎn)生影響。

    表3 2013年不同湖區(qū)各月有效波高最大值(m)

    2.3 典型風(fēng)向下微風(fēng)引起的波浪季節(jié)分布特征

    由于2013年太湖盛行微風(fēng),因此借助所建SWAN模型分析微風(fēng)風(fēng)速條件下西南、東南、東北和西北風(fēng)在太湖形成的風(fēng)浪空間分布.在西南風(fēng)作用下,波高由西南沿岸向湖心逐漸變大,繞過(guò)西山島,在梅梁灣和貢湖口達(dá)到最大(圖5A).在夏季盛行的東南風(fēng)作用下,波高由箭湖東茭咀-南部沿岸-湖心逐漸變大,在梅梁湖口以西的湖心區(qū)域達(dá)到最大(圖5B).在東北風(fēng)作用下,波高最大值在西部沿岸附近的湖心區(qū)域,其次為大太湖和西部沿岸,各湖灣的有效波高基本在0.15~0.18m(圖5C).而西北風(fēng)引起的風(fēng)浪的有效波高值則由西北沿岸區(qū)向湖心逐漸增大,西山島附近達(dá)到最大;在西山島的阻礙作用下,胥湖和東太湖的波高明顯減小(圖5D)。

    風(fēng)向可以通過(guò)改變風(fēng)區(qū)長(zhǎng)度來(lái)影響風(fēng)浪生消.不過(guò),不論何種風(fēng)向下,大太湖的有效波高始終最大,在0.16~0.22m之間.而該湖區(qū)也是太湖水動(dòng)力擾動(dòng)最為頻繁和劇烈的區(qū)域.沈吉等[30]觀測(cè)到湖心大部分湖底露出黃褐色硬質(zhì)黃土層,局部地區(qū)有小于3cm的上覆浮泥,這可能與該區(qū)域始終頻繁遭受強(qiáng)風(fēng)浪擾動(dòng)有關(guān).東太湖和胥湖的有效波高都是最小的,在0.10~0.16m之間(圖5).這兩塊水域相對(duì)封閉,風(fēng)區(qū)長(zhǎng)度較小,不利于風(fēng)浪的成長(zhǎng)和發(fā)展.這樣的水動(dòng)力條件大大降低了沉水植物的機(jī)械損傷,為沉水植物的生長(zhǎng)提供了良好的環(huán)境;同時(shí),而沉水植物繁盛又可以起到消浪的作用,兩者間存在著正反饋機(jī)制。

    圖5 2013年典型風(fēng)向下微風(fēng)引起的太湖有效波高的季節(jié)分布Fig.5 Seasonal distribution of significant wave height in Lake Taihu under different wind directions in 2013

    3 結(jié)論

    1) 通過(guò)2014年8月14日6:10至8月17日5:40實(shí)測(cè)有效波高與模擬有效波高比較發(fā)現(xiàn),其變化趨勢(shì)基本一致,能夠緊密分布在直線y=0.946x兩側(cè),判定系數(shù)R2=0.960(樣本數(shù)n=288);相對(duì)誤差范圍為0.03%~28.7%,均值為10.7%;實(shí)測(cè)值與模擬值呈極顯著正相關(guān)(r=0.83,P<0.01),確定性系數(shù)為0.84.可見(jiàn)SWAN模型能夠較好地體現(xiàn)太湖風(fēng)浪特征。

    2) 受岸線、地形和島嶼等地理因素影響,2013年大太湖的風(fēng)浪總是最強(qiáng)的,其有效波高均值為0.523m;而東太湖風(fēng)浪最小,有效波高均值為0.305m.受盛行風(fēng)場(chǎng)季節(jié)變化影響,太湖春、夏季有效波高均值明顯大于秋、冬季.太湖波浪的能量主要來(lái)源于風(fēng)場(chǎng),其有效波高隨風(fēng)速增大,兩者呈極顯著正相關(guān)(r=0.47,P<0.01)。

    3) 風(fēng)向可以通過(guò)改變風(fēng)區(qū)長(zhǎng)度來(lái)影響風(fēng)浪生消.在偏東風(fēng)作用下,太湖湖西區(qū)的風(fēng)浪大于東部湖區(qū);而受盛行于冬季的偏北風(fēng)影響,太湖南部水域風(fēng)浪要大于北部.同時(shí),太湖風(fēng)浪的時(shí)空分布特征是造成太湖水質(zhì)參數(shù)、沉積物和水生植物空間分布差異的重要原因之一。

    參考文獻(xiàn)4

    [1]孫順才, 黃漪平. 太湖. 北京:海洋出版社, 1993。

    [2]秦伯強(qiáng), 胡維平, 高光等. 太湖沉積物懸浮的動(dòng)力機(jī)制及內(nèi)源釋放的概念性模式. 科學(xué)通報(bào), 2003, 48(17): 1822-1831。

    [3]Wu TF, Qin BQ, Zhu GWetal. Modeling of turbidity dynamics caused by wind-induced waves and current in the Taihu Lake.InternationalJournalofSedimentResearch, 2013, 28(2): 139-148。

    [4]Jin KR, Ji ZG. Calibration and verification of a spectral wind-wave model for Lake Okeechobee.OceanEngineering, 2001, 28(5): 571-584。

    [5]Bailey MC, Hamilton DP. Wind induced sediment resuspension: a lake-wide model.EcologicalModeling, 1997, 99(2): 217-228。

    [6]李一平, 逄勇, 劉興平等. 太湖波浪數(shù)值模擬. 湖泊科學(xué), 2008, 20(1): 117-122. DOI 10. 18307/2008. 0118。

    [7]劉興平. 太湖波浪過(guò)程的數(shù)值模擬. 江蘇大學(xué)學(xué)報(bào): 自然科學(xué)版, 2009, 30(1): 80-85。

    [8]李一平, 逄勇, 羅瀲蔥. 波流作用下太湖水體懸浮物輸運(yùn)實(shí)驗(yàn)及模擬. 水科學(xué)進(jìn)展, 2009, 20(5): 701-706。

    [9]羅瀲蔥, 秦伯強(qiáng). 太湖波浪與湖流對(duì)沉積物再懸浮不同影響的研究. 水文, 2003, 23(3): 1-4。

    [10]湯露露, 王鵬, 姚琪. 太湖湖流、波浪、沉積物的三維數(shù)值模擬. 水資源保護(hù), 2011, 27(2): 1-5。

    [11]陶蓉茵. 不同風(fēng)場(chǎng)條件下太湖波浪數(shù)值模擬的研究及其應(yīng)用[學(xué)位論文]. 南京: 南京信息工程大學(xué), 2012。

    [12]劉聚濤, 高俊峰, 姜加虎等. 基于突變理論的太湖藍(lán)藻水華危險(xiǎn)性分區(qū)評(píng)價(jià). 湖泊科學(xué), 2010, 22(4): 488-494. DOI 10. 18307/2010. 0403。

    [13]Hasselmann S, Hasselmann K, Allender JHetal. Computation and parameterizations of the nonlinear energy transfer in a gravity wave spectrum. Part Ⅱ: Parameterizations of the nonlinear energy transfer for application in wave models.JournalofPhysicalOceangraphy, 1985, 15(11): 1378-1391。

    [14]Phillips OM. On the generation of waves by turbulent wind.JournalofFluidMechanics, 1957, 2(5): 417-445。

    [15]Miles JW. Hamiltonian formulations for surface waves.AppliedScientificResearch, 1981, 37(1/2): 103-110。

    [16]Battjes JA, Stive MJF. Calibration and verification of a dissipation model for random breaking waves.JournalofGeophysicalResearch, 1985, 90(C5): 9159-9167。

    [17]Eldeberky Y. Nonlinear transformation of wave spectra in the nearshore zone[Dissertation]. Netherlands: Delft University of Technology, 1996。

    [18]Eldeberky Y, Battjes JA. Parameterization of triad interactions in wave energy models. Gdansk, Poland: Proceeding Coastal Dynamics Conference, 1995: 140-148。

    [19]許遐禎, 陶蓉茵, 趙巧華等. 大型淺水湖泊太湖波浪特征及其對(duì)風(fēng)場(chǎng)的敏感性分析. 湖泊科學(xué), 2013, 25(1): 55-64. DOI 10. 18307/2013. 0108。

    [20]曹勇, 孫從軍, 高陽(yáng)俊等. 淀山湖水域風(fēng)浪要素計(jì)算研究. 長(zhǎng)江流域資源與環(huán)境, 2011, 20(Z1): 39-43。

    [21]Jackson HO, Starrett WC. Turbidity and sedimentation at Lake Chautauqua, lllinois.JournalofWildlifeManagement, 1959, 23: 157-168。

    [22]Hamilton DP, Mitchell SF. An empirical model for sediment resuspension in shallow lakes.Hydrobiologia, 1996, 317(3): 209-220。

    [23]李一平, 逄勇, 李勇. 水動(dòng)力作用下太湖底泥的再懸浮通量. 水利學(xué)報(bào), 2007, 38(5): 558-564。

    [24]逄勇, 李一平, 顏潤(rùn)潤(rùn)等. 內(nèi)外源共同作用對(duì)太湖營(yíng)養(yǎng)鹽貢獻(xiàn)量研究. 水利學(xué)報(bào), 2008, 39(9): 1051-1059。

    [25]王震, 鄒華, 楊桂軍等. 太湖葉綠素a的時(shí)空分布特征及其與環(huán)境因子的相關(guān)關(guān)系. 湖泊科學(xué), 2014, 26(4): 567-575. DOI 10. 18307/2014. 0411。

    [26]李英杰, 胡小貞, 胡社榮等. 草、藻型湖泊水體生態(tài)及理化特性的實(shí)驗(yàn)對(duì)比. 生態(tài)環(huán)境學(xué)報(bào), 2009, 18(5): 1649-1654。

    [27]李一平, 逄勇, 劉興平. 太湖湖流對(duì)波浪的影響機(jī)制研究. 水利學(xué)報(bào), 2007, 38(S1): 303-308。

    [28]Lind OT, Davaloslind LO, Chrzanowski THetal. Inorganic turbidity and the failure of fishery models.InternationaleRevuederGesamtenHydrobiologie, 1994, 79(1): 7-16。

    [29]Hofmann H, Lorke A, Peeters F. The relative importance of wind and ship waves in the littoral zone of a large lake.LimnologyandOceanography, 2008, 53(1): 368-380。

    [30]沈吉, 袁和忠, 劉恩峰等. 太湖表層沉積物的空間分布與層序特征分析. 科學(xué)通報(bào), 2010, 55(36): 3516-3524。

    J.LakeSci.(湖泊科學(xué)), 2016, 28(1): 217-224

    ?2016 byJournalofLakeSciences

    WANG Zhen1, WU Tingfeng2**, ZOU Hua1***, JIA Xiaowang3, HUANG Lie4, LIANG Chaorong4& ZHANG Zhihao4

    (1:SchoolofEnvironmentandCivilEngineering,JiangnanUniversity,Wuxi214122,P.R.China)

    (2:StateKeyLaboratoryofLakeScienceandEnvironment,NanjingInstituteofGeographyandLimnology,ChineseAcademyofSciences,Nanjing210008,P.R.China)

    (3:WuxiBranchofJiangsuProvinceHydrologyandWaterResourceInvestigationBureau,Wuxi214031,P.R.China)

    (4:CollegeofHydrologyandWaterResources,HohaiUniversity,Nanjing210098,P.R.China)

    Abstract:To understand the spatial and temporal distribution of wind and wave in Lake Taihu, wave observatory was set up in the center of lake. It was found that SWAN model could simulate the waves well in Lake Taihu by using the wave data recorded. Based on the calibrated model, waves of different seasons in Lake Taihu were simulated under the natural wind conditions in 2013. The results showed that due to shoreline, topography and islands, waves in the large Taihu water always showed the strongest with the mean significant wave height of 0.523m, while waves in east Lake Taihu were the smallest with the mean significant wave height of 0.305m. The values of mean significant wave height were higher in spring and summer than in autumn and winter because of seasonal changes. Wave energy in Lake Taihu was driven mainly by wind, leading to the significant wave height increased with wind speed. A significant positive correlation was found between significant wave height and wind speed. The wind direction affected growth and extinction of waves with changing the length of wind fetch. Under the easterly wind, waves in western lake was greater than that in eastern lake; while under northerly wind in winter, waves in southern lake was greater than that in northern lake. Meanwhile, the spatial and temporal distribution of waves was one of the major reasons to form a different spatial distribution of water quality, sediment and aquatic macrophytes in Lake Taihu。

    Keywords:Partitions of Lake Taihu; wind and wave; seasonal dynamics; spatial distribution; SWAN model; aquatic macrophytes

    通信作者王震1,吳挺峰2*;E-mail:tfwu@niglas.ac.cn.,鄒華1*共同通信作者;E-mail: zouhua@jiangnan.edu.cn.,賈小網(wǎng)3,黃列4,梁朝榮4,張志浩4

    基金項(xiàng)目收修改稿.王震(1990~),男,碩士研究生;E-mail:wzh.wz8@163.com.

    DOI10.18307/2016.0125

    猜你喜歡
    水生植物空間分布風(fēng)浪
    坐在同一條船上
    風(fēng)浪干擾條件下艦船航向保持非線性控制系統(tǒng)
    不怕風(fēng)浪的小船
    East Side Story
    鐵嶺市凡河流域水生植物調(diào)查研究
    白龜山濕地重金屬元素分布特征及其來(lái)源分析
    綠色科技(2016年20期)2016-12-27 17:34:13
    基于GIS技術(shù)的福建省柳葉白前資源適宜性空間分布研究
    江蘇省臭氧污染變化特征
    科技視界(2016年18期)2016-11-03 23:51:58
    鐵路客流時(shí)空分布研究綜述
    城市生活污水現(xiàn)狀及植物凈化污水研究
    亚洲无线在线观看| 黄色片一级片一级黄色片| av天堂在线播放| 香蕉国产在线看| 久久中文字幕人妻熟女| 黄片小视频在线播放| 国产精品综合久久久久久久免费 | 波多野结衣高清无吗| 国产在线观看jvid| 亚洲精品粉嫩美女一区| 国产精华一区二区三区| 两人在一起打扑克的视频| 人妻久久中文字幕网| 亚洲人成电影免费在线| 日本在线视频免费播放| 91大片在线观看| 级片在线观看| 久久热在线av| 9热在线视频观看99| 午夜精品国产一区二区电影| 久热这里只有精品99| av有码第一页| 国产精品av久久久久免费| 69av精品久久久久久| 欧美日韩福利视频一区二区| 黄片播放在线免费| 女人爽到高潮嗷嗷叫在线视频| 国产精品亚洲美女久久久| 欧美久久黑人一区二区| 国产三级在线视频| 精品国产一区二区三区四区第35| 人人妻,人人澡人人爽秒播| 精品人妻1区二区| 成人国产一区最新在线观看| 日韩成人在线观看一区二区三区| 国产高清视频在线播放一区| 国产精品99久久99久久久不卡| 久久精品亚洲熟妇少妇任你| 亚洲国产欧美网| www.熟女人妻精品国产| 最新在线观看一区二区三区| 国产国语露脸激情在线看| 中国美女看黄片| 亚洲国产欧美一区二区综合| 国产高清有码在线观看视频 | 国产精品一区二区三区四区久久 | 亚洲第一欧美日韩一区二区三区| 国产亚洲欧美98| 久久这里只有精品19| 成年人黄色毛片网站| 免费在线观看亚洲国产| 黄色a级毛片大全视频| 色婷婷久久久亚洲欧美| 欧美av亚洲av综合av国产av| 黄色视频,在线免费观看| 51午夜福利影视在线观看| 久久久国产成人免费| 国产精品久久电影中文字幕| 免费无遮挡裸体视频| 乱人伦中国视频| 18禁国产床啪视频网站| 麻豆av在线久日| 欧美黄色淫秽网站| 亚洲色图综合在线观看| 亚洲精品国产区一区二| 91麻豆精品激情在线观看国产| 亚洲九九香蕉| av福利片在线| 亚洲欧美日韩高清在线视频| 亚洲国产精品合色在线| 精品久久久久久,| videosex国产| 真人一进一出gif抽搐免费| 欧美日韩亚洲国产一区二区在线观看| 9191精品国产免费久久| 日韩中文字幕欧美一区二区| 国产高清激情床上av| 大陆偷拍与自拍| 久久九九热精品免费| 成熟少妇高潮喷水视频| 男男h啪啪无遮挡| 激情视频va一区二区三区| 国产主播在线观看一区二区| 国产成人精品无人区| 国产麻豆69| 丰满人妻熟妇乱又伦精品不卡| 搞女人的毛片| 人人妻,人人澡人人爽秒播| АⅤ资源中文在线天堂| 欧美色视频一区免费| 制服人妻中文乱码| 国产精品免费视频内射| 精品一品国产午夜福利视频| 久久久久久久久久久久大奶| 成人18禁在线播放| 在线观看66精品国产| 国产又爽黄色视频| 日本免费a在线| 欧美日韩福利视频一区二区| 高潮久久久久久久久久久不卡| 看免费av毛片| 欧美成人午夜精品| 在线观看免费视频网站a站| 国产精品爽爽va在线观看网站 | 18禁国产床啪视频网站| 丝袜美腿诱惑在线| av在线天堂中文字幕| 丝袜人妻中文字幕| 男人操女人黄网站| 日韩有码中文字幕| tocl精华| 一区二区三区精品91| 国产熟女午夜一区二区三区| 操美女的视频在线观看| 午夜久久久久精精品| 日韩欧美国产在线观看| 久久久久久久久免费视频了| 啦啦啦观看免费观看视频高清 | x7x7x7水蜜桃| 亚洲第一欧美日韩一区二区三区| 窝窝影院91人妻| 黄色成人免费大全| 999精品在线视频| 大香蕉久久成人网| 一卡2卡三卡四卡精品乱码亚洲| 搡老熟女国产l中国老女人| 午夜福利高清视频| 亚洲精品久久成人aⅴ小说| 亚洲狠狠婷婷综合久久图片| 亚洲精品久久国产高清桃花| 国产区一区二久久| 日本vs欧美在线观看视频| 免费久久久久久久精品成人欧美视频| 中文字幕最新亚洲高清| 一级毛片女人18水好多| 久久精品国产亚洲av高清一级| 国产三级在线视频| 精品人妻1区二区| 一级a爱片免费观看的视频| 成人特级黄色片久久久久久久| 久久久久久免费高清国产稀缺| 神马国产精品三级电影在线观看 | 天天添夜夜摸| 高清毛片免费观看视频网站| 精品无人区乱码1区二区| 熟女少妇亚洲综合色aaa.| 怎么达到女性高潮| 欧美成人午夜精品| 国产精品二区激情视频| 日本 欧美在线| 亚洲avbb在线观看| 亚洲色图av天堂| 男人舔女人下体高潮全视频| 大型av网站在线播放| 91麻豆av在线| 成人国产综合亚洲| 国产主播在线观看一区二区| 男女下面插进去视频免费观看| 亚洲成国产人片在线观看| 中文字幕av电影在线播放| 午夜视频精品福利| 国产精品日韩av在线免费观看 | 国产亚洲精品久久久久久毛片| 亚洲国产精品sss在线观看| 亚洲狠狠婷婷综合久久图片| 日本撒尿小便嘘嘘汇集6| 91精品国产国语对白视频| www.自偷自拍.com| 国产一区二区三区在线臀色熟女| 成人亚洲精品av一区二区| 看免费av毛片| 亚洲专区国产一区二区| 他把我摸到了高潮在线观看| 精品熟女少妇八av免费久了| 亚洲午夜理论影院| 欧美av亚洲av综合av国产av| 精品免费久久久久久久清纯| 欧美大码av| 国产精品亚洲一级av第二区| 国产欧美日韩精品亚洲av| 国产欧美日韩一区二区精品| 国产亚洲精品一区二区www| 国产野战对白在线观看| 亚洲欧美一区二区三区黑人| 国内久久婷婷六月综合欲色啪| 久久婷婷成人综合色麻豆| 国产一区二区三区综合在线观看| 两人在一起打扑克的视频| 亚洲专区字幕在线| 免费少妇av软件| 国产成人精品无人区| www.精华液| 九色亚洲精品在线播放| 国产精品久久久久久精品电影 | 国产区一区二久久| 亚洲一区高清亚洲精品| 女人被狂操c到高潮| 国内精品久久久久精免费| 国产精品久久久人人做人人爽| 午夜福利,免费看| 老鸭窝网址在线观看| 亚洲第一欧美日韩一区二区三区| 91九色精品人成在线观看| 999久久久国产精品视频| 国产精品免费一区二区三区在线| 免费看美女性在线毛片视频| 国产成人影院久久av| 日本撒尿小便嘘嘘汇集6| 亚洲精品国产色婷婷电影| 50天的宝宝边吃奶边哭怎么回事| 老司机深夜福利视频在线观看| 国产亚洲精品久久久久久毛片| 亚洲av成人av| 一区在线观看完整版| 国产激情欧美一区二区| 丁香欧美五月| 国产亚洲av嫩草精品影院| 黄色a级毛片大全视频| 大型黄色视频在线免费观看| 日韩视频一区二区在线观看| 日韩欧美在线二视频| 午夜福利影视在线免费观看| 婷婷六月久久综合丁香| 99在线视频只有这里精品首页| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 高清黄色对白视频在线免费看| 欧美国产日韩亚洲一区| 一区二区三区激情视频| 国产一区二区三区综合在线观看| 免费在线观看亚洲国产| 精品少妇一区二区三区视频日本电影| 日韩三级视频一区二区三区| 中出人妻视频一区二区| 国产99白浆流出| 99在线人妻在线中文字幕| 久久久精品国产亚洲av高清涩受| 少妇熟女aⅴ在线视频| 日韩免费av在线播放| 男女之事视频高清在线观看| 久久草成人影院| 变态另类成人亚洲欧美熟女 | 久久精品成人免费网站| 亚洲色图综合在线观看| 国产在线观看jvid| 日韩大码丰满熟妇| 久久伊人香网站| 欧美日韩黄片免| 久久精品人人爽人人爽视色| 精品久久蜜臀av无| 黄频高清免费视频| 精品卡一卡二卡四卡免费| 精品国产乱子伦一区二区三区| 人妻丰满熟妇av一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看日本一区| 欧美最黄视频在线播放免费| 丁香六月欧美| 一本大道久久a久久精品| 搡老熟女国产l中国老女人| 韩国av一区二区三区四区| 91在线观看av| 女人爽到高潮嗷嗷叫在线视频| 黄色成人免费大全| 大型黄色视频在线免费观看| 男女下面进入的视频免费午夜 | 久久青草综合色| 成年版毛片免费区| 天天躁夜夜躁狠狠躁躁| 久久婷婷人人爽人人干人人爱 | 亚洲 欧美一区二区三区| 国产精品乱码一区二三区的特点 | 91国产中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 国产成人免费无遮挡视频| 国产乱人伦免费视频| 国产野战对白在线观看| 欧美激情极品国产一区二区三区| 久久亚洲真实| 亚洲成国产人片在线观看| 看片在线看免费视频| 在线观看免费视频日本深夜| 国产xxxxx性猛交| 淫秽高清视频在线观看| 欧美乱妇无乱码| √禁漫天堂资源中文www| 亚洲自偷自拍图片 自拍| 天天躁夜夜躁狠狠躁躁| 男女床上黄色一级片免费看| av超薄肉色丝袜交足视频| 欧美一级毛片孕妇| 色播在线永久视频| 久久婷婷成人综合色麻豆| 精品国产美女av久久久久小说| 51午夜福利影视在线观看| 非洲黑人性xxxx精品又粗又长| 久久热在线av| 亚洲专区国产一区二区| 精品国产一区二区三区四区第35| 黄色视频,在线免费观看| 久99久视频精品免费| 人人澡人人妻人| 少妇 在线观看| 午夜a级毛片| 日日干狠狠操夜夜爽| 欧美色欧美亚洲另类二区 | 18禁观看日本| 人人妻人人爽人人添夜夜欢视频| 久久天躁狠狠躁夜夜2o2o| 可以免费在线观看a视频的电影网站| 亚洲av片天天在线观看| 亚洲成人久久性| 午夜福利,免费看| 18禁裸乳无遮挡免费网站照片 | 色精品久久人妻99蜜桃| 成人亚洲精品av一区二区| 中文字幕精品免费在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 中文字幕人妻丝袜一区二区| 久久久久国内视频| 法律面前人人平等表现在哪些方面| 国产成+人综合+亚洲专区| 神马国产精品三级电影在线观看 | 久久 成人 亚洲| 亚洲精品国产区一区二| 深夜精品福利| 1024视频免费在线观看| svipshipincom国产片| 久久精品国产综合久久久| 亚洲国产欧美网| 久久婷婷成人综合色麻豆| 嫩草影视91久久| 欧美激情极品国产一区二区三区| 好男人在线观看高清免费视频 | 久久热在线av| 夜夜看夜夜爽夜夜摸| 免费av毛片视频| 免费看十八禁软件| 亚洲精品国产色婷婷电影| 麻豆成人av在线观看| 国产一卡二卡三卡精品| a在线观看视频网站| 制服人妻中文乱码| 女同久久另类99精品国产91| 日本 av在线| 1024视频免费在线观看| 日韩欧美国产在线观看| 亚洲成人免费电影在线观看| 亚洲成人国产一区在线观看| 嫁个100分男人电影在线观看| 国产av一区在线观看免费| 麻豆成人av在线观看| 欧美老熟妇乱子伦牲交| 国产精华一区二区三区| 国产精品一区二区免费欧美| 色哟哟哟哟哟哟| 黄色视频不卡| 非洲黑人性xxxx精品又粗又长| 久久精品国产综合久久久| 国产亚洲精品第一综合不卡| 91av网站免费观看| 国产精品亚洲av一区麻豆| 无限看片的www在线观看| 国产xxxxx性猛交| 国产精品久久久久久人妻精品电影| 亚洲精品国产一区二区精华液| 国产精品一区二区精品视频观看| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 午夜免费鲁丝| 伊人久久大香线蕉亚洲五| 99国产综合亚洲精品| 亚洲欧美精品综合久久99| 欧美日韩精品网址| 一进一出好大好爽视频| 成年人黄色毛片网站| 可以在线观看毛片的网站| 日韩欧美国产一区二区入口| 久久人妻福利社区极品人妻图片| 国产区一区二久久| 亚洲午夜精品一区,二区,三区| 男女下面进入的视频免费午夜 | 国产欧美日韩一区二区三| 国产亚洲av嫩草精品影院| 亚洲精品av麻豆狂野| 午夜免费成人在线视频| 性色av乱码一区二区三区2| 亚洲国产中文字幕在线视频| 亚洲欧美激情在线| 国产午夜精品久久久久久| 97超级碰碰碰精品色视频在线观看| 又黄又粗又硬又大视频| 一区福利在线观看| 91国产中文字幕| 免费在线观看日本一区| 日韩精品青青久久久久久| 亚洲七黄色美女视频| 正在播放国产对白刺激| 国产一区二区三区视频了| 成人免费观看视频高清| www国产在线视频色| 黑人巨大精品欧美一区二区mp4| 天堂动漫精品| 日本一区二区免费在线视频| 国产免费男女视频| 欧美午夜高清在线| 69精品国产乱码久久久| 一本综合久久免费| 精品国产美女av久久久久小说| 99精品久久久久人妻精品| 久久精品亚洲精品国产色婷小说| 久久婷婷人人爽人人干人人爱 | 成人特级黄色片久久久久久久| 久久久久久久午夜电影| 精品午夜福利视频在线观看一区| 一二三四在线观看免费中文在| 变态另类丝袜制服| 可以免费在线观看a视频的电影网站| 亚洲欧美日韩另类电影网站| 亚洲精品美女久久av网站| 性少妇av在线| 一级a爱视频在线免费观看| 久久久水蜜桃国产精品网| 两人在一起打扑克的视频| 久久国产精品影院| av有码第一页| 欧美另类亚洲清纯唯美| 在线观看www视频免费| 国产一区二区三区在线臀色熟女| 免费久久久久久久精品成人欧美视频| 妹子高潮喷水视频| 免费观看人在逋| 两个人看的免费小视频| 涩涩av久久男人的天堂| 欧美日韩精品网址| 久久欧美精品欧美久久欧美| 亚洲av熟女| 午夜免费鲁丝| 国产精品久久久av美女十八| 大香蕉久久成人网| 好看av亚洲va欧美ⅴa在| 一进一出好大好爽视频| 精品一区二区三区四区五区乱码| 一边摸一边抽搐一进一出视频| 99精品欧美一区二区三区四区| 夜夜夜夜夜久久久久| 午夜两性在线视频| 国产一区二区三区在线臀色熟女| 成熟少妇高潮喷水视频| 女人被躁到高潮嗷嗷叫费观| 久久精品影院6| 男人的好看免费观看在线视频 | 亚洲精品一区av在线观看| 亚洲精华国产精华精| 免费在线观看日本一区| 成人三级黄色视频| 欧美成狂野欧美在线观看| 亚洲专区字幕在线| 色综合亚洲欧美另类图片| 日韩一卡2卡3卡4卡2021年| 国产精品1区2区在线观看.| 亚洲男人天堂网一区| 精品欧美一区二区三区在线| 丝袜在线中文字幕| 久久久久久久精品吃奶| 国产1区2区3区精品| 日本 av在线| 很黄的视频免费| 日韩中文字幕欧美一区二区| 欧美日韩乱码在线| 国产精品电影一区二区三区| 久久国产精品男人的天堂亚洲| av福利片在线| 久热爱精品视频在线9| 日韩视频一区二区在线观看| xxx96com| 精品国产一区二区三区四区第35| 国产欧美日韩一区二区精品| 欧美日韩瑟瑟在线播放| 国产欧美日韩精品亚洲av| 欧美另类亚洲清纯唯美| 日韩免费av在线播放| 亚洲精品久久国产高清桃花| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| 色综合欧美亚洲国产小说| 91字幕亚洲| 又大又爽又粗| 99精品在免费线老司机午夜| 久久精品aⅴ一区二区三区四区| 亚洲欧美激情综合另类| 亚洲中文字幕日韩| 99国产精品一区二区蜜桃av| 久久国产亚洲av麻豆专区| 宅男免费午夜| 大码成人一级视频| 黑人巨大精品欧美一区二区蜜桃| 一区二区三区精品91| 老熟妇乱子伦视频在线观看| 亚洲av电影在线进入| 日韩大尺度精品在线看网址 | 91麻豆av在线| 此物有八面人人有两片| 久久午夜亚洲精品久久| 老司机在亚洲福利影院| 国产av精品麻豆| 99精品在免费线老司机午夜| 免费av毛片视频| 亚洲av电影不卡..在线观看| av中文乱码字幕在线| 99精品久久久久人妻精品| 夜夜看夜夜爽夜夜摸| 久久亚洲真实| 久久精品亚洲熟妇少妇任你| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 午夜免费观看网址| 成熟少妇高潮喷水视频| 国产三级在线视频| 三级毛片av免费| 12—13女人毛片做爰片一| 桃红色精品国产亚洲av| 国产精品亚洲美女久久久| 亚洲无线在线观看| 在线国产一区二区在线| 少妇裸体淫交视频免费看高清 | 国产午夜精品久久久久久| 久久热在线av| 免费av毛片视频| 亚洲国产中文字幕在线视频| 亚洲欧美日韩无卡精品| 性少妇av在线| www.自偷自拍.com| 久99久视频精品免费| 国产一区二区三区视频了| www.www免费av| 久热爱精品视频在线9| av网站免费在线观看视频| www日本在线高清视频| 久久中文字幕人妻熟女| 精品卡一卡二卡四卡免费| 女性生殖器流出的白浆| 99精品在免费线老司机午夜| 欧美黑人欧美精品刺激| 国产人伦9x9x在线观看| 国产精品久久久人人做人人爽| 精品日产1卡2卡| 啦啦啦免费观看视频1| 中文字幕色久视频| 日韩成人在线观看一区二区三区| 操美女的视频在线观看| av中文乱码字幕在线| 日本vs欧美在线观看视频| 国产精品久久久久久亚洲av鲁大| 女警被强在线播放| 大型黄色视频在线免费观看| 免费在线观看黄色视频的| 午夜视频精品福利| 在线观看日韩欧美| 老熟妇仑乱视频hdxx| 在线观看日韩欧美| 国产精品综合久久久久久久免费 | 精品一区二区三区av网在线观看| 丰满的人妻完整版| 国产精品一区二区在线不卡| 少妇被粗大的猛进出69影院| 高清在线国产一区| 1024香蕉在线观看| 精品第一国产精品| 国产成年人精品一区二区| 久久久久久久久免费视频了| 日韩欧美三级三区| 精品久久久久久久久久免费视频| 男女床上黄色一级片免费看| 午夜免费鲁丝| 国产午夜精品久久久久久| 真人做人爱边吃奶动态| 国产一区二区激情短视频| 中文字幕高清在线视频| 亚洲无线在线观看| 咕卡用的链子| 一夜夜www| 曰老女人黄片| 亚洲五月天丁香| 午夜两性在线视频| 制服诱惑二区| 18禁观看日本| av有码第一页| 精品久久久久久久人妻蜜臀av | 亚洲人成电影免费在线| 嫩草影视91久久| 亚洲人成伊人成综合网2020| 午夜福利18| 欧美日韩亚洲国产一区二区在线观看| 在线观看免费视频网站a站| 国产片内射在线| 99热只有精品国产| 可以在线观看毛片的网站| 国产精品爽爽va在线观看网站 | 99精品久久久久人妻精品| 好男人电影高清在线观看| 老司机午夜十八禁免费视频| 在线视频色国产色| 成年女人毛片免费观看观看9| avwww免费| 两个人看的免费小视频| 19禁男女啪啪无遮挡网站| 真人做人爱边吃奶动态| 久久精品人人爽人人爽视色| 在线国产一区二区在线| 搡老妇女老女人老熟妇| 国产乱人伦免费视频| 久久久久九九精品影院| 少妇粗大呻吟视频| 亚洲成a人片在线一区二区|