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

    海南島南渡江河口枯季大小潮的鹽度變化特征

    2012-08-14 01:21:56龔文平王道儒趙軍鵬莫文淵
    海洋通報(bào) 2012年6期
    關(guān)鍵詞:南渡江小潮落潮

    龔文平,王道儒,趙軍鵬,莫文淵

    (1.中山大學(xué) 海洋學(xué)院近岸海洋科技研究中心,廣東 廣州 510275;2.海南省海洋開發(fā)規(guī)劃設(shè)計(jì)研究院,海南 ???570125)

    南渡江是海南省最大的河流,其河口長(zhǎng)度較短(25 km左右),口外是流向復(fù)雜的瓊州海峽,口內(nèi)發(fā)育駐波??菁緯r(shí),南渡江河口呈現(xiàn)出鹽水楔河口的特征(趙軍鵬,2011),相應(yīng)地其鹽分輸運(yùn)與鹽水入侵有其自身的特點(diǎn)。近幾十年來(lái)因河口段的人類活動(dòng)(主要為河道采砂)而產(chǎn)生鹽水入侵的加劇,對(duì)當(dāng)?shù)氐纳a(chǎn)和生活用水產(chǎn)生不利影響,有必要加強(qiáng)對(duì)南渡江鹽水入侵的研究。

    國(guó)際上對(duì)各種類型河口的水動(dòng)力與鹽分輸運(yùn)的過(guò)程與機(jī)理進(jìn)行了大量的研究(MacCready,2007;Lerczak etal,2006;Ralston etal,2010),國(guó)內(nèi)對(duì)長(zhǎng)江、珠江等大河口的鹽水入侵的研究也相當(dāng)深入(Wu etal,2010;Gong etal,2011)。但對(duì)小河口鹽水入侵的研究相對(duì)較少。南渡江河口為典型的小河口,潮汐作用弱,徑流作用相對(duì)較強(qiáng),其鹽分向陸輸運(yùn)的機(jī)理如何?鹽水入侵在枯季時(shí)大小潮期的變化怎樣?人類活動(dòng)對(duì)鹽水入侵的影響又如何?這些問(wèn)題的解決將豐富國(guó)內(nèi)外關(guān)于河口動(dòng)力學(xué)的研究,同時(shí)對(duì)南渡江鹽水入侵的防控具有重要意義。

    南渡江發(fā)源于海南島中部山區(qū),向北注入瓊州海峽。全長(zhǎng)334 km,流域面積達(dá)7 033 km2。龍?zhí)翂沃芬韵聻楹涌诙危ê佣伍L(zhǎng)26.5 km)。南渡江在麻余村附近分為三汊,即干流(長(zhǎng)4.2 km),橫溝河(長(zhǎng)5 km)和海甸溪(長(zhǎng)6 km)。干流與橫溝河之間為新埠島,橫溝河與海甸溪之間為海甸島(圖1)。干流是南渡江主要的水沙排泄通道,橫溝河內(nèi)的漲落潮勢(shì)均力敵,而海甸溪?jiǎng)t主要受漲潮流控制。

    圖1 南渡江區(qū)位圖及觀測(cè)站位分布與縱斷面位置

    河口段徑流主要來(lái)自南渡江上游地區(qū),且洪枯季變化明顯,南渡江約80%的徑流量來(lái)自洪季的6-10月。據(jù)龍?zhí)琳径嗄晁馁Y料統(tǒng)計(jì),實(shí)測(cè)最大年平均流量為296m3/s(1973年),最小年平均流量為74.6m3/s(2004年)。2004-2007年期間,年平均流量為108.1m3/s,徑流量為34.7×108m3。其中干流河口輸水量占總徑流量的70%以上,而橫溝河占30%左右,海甸溪水沙則表現(xiàn)為凈進(jìn)。

    本區(qū)多年平均潮差為1.21m,最大潮差2.5m左右,為弱潮河口。潮汐系數(shù)為3.86,屬于不正規(guī)全日潮。瓊州海峽潮流為往復(fù)流性質(zhì),有漲潮東流、西流,落潮東流、西流4種流動(dòng)形式,以漲潮東流和落潮西流為主,轉(zhuǎn)流一般發(fā)生在平均水位附近,平均水位以上以東流為主,平均水位以下以西流為主。

    圖2 沿河口的縱斷面的水深變化及站3、4所在位置

    1 材料和方法

    本研究主要根據(jù)實(shí)測(cè)資料進(jìn)行分析。趙軍鵬(2011)采用FVCOM模型對(duì)南渡江河口的鹽水入侵過(guò)程及機(jī)理進(jìn)行了細(xì)致的研究,但模型由于時(shí)空尺度的限制以及對(duì)物理過(guò)程的簡(jiǎn)化,對(duì)一些細(xì)微的物理過(guò)程(如水躍(jump))尚不能完全反映,因此有必要對(duì)實(shí)測(cè)資料進(jìn)行更細(xì)致的分析。

    海南省海洋開發(fā)規(guī)劃設(shè)計(jì)研究院于2009年枯水期在南渡江河口進(jìn)行了兩航次的現(xiàn)場(chǎng)水文同步觀測(cè)(小潮期為2009年02月07日09:00至2009年02月08日10∶00;大潮期為2009年02月11日10∶00至2009年02月12日11∶00),觀測(cè)內(nèi)容有水位、水流、鹽度、懸沙含量等。

    觀測(cè)期間設(shè)定6個(gè)海流觀測(cè)站,多船同步作業(yè),大小潮期分別連續(xù)觀測(cè)26 h,采樣的時(shí)間間隔為1 h,測(cè)量?jī)x器采用日本ALEC公司EM自容式海流計(jì)。根據(jù)水深的不同,垂向上采用不同的分層方法測(cè)量水流,當(dāng)水深≤5m采用三點(diǎn)法(表、0.6H、底層);1號(hào)站采用六點(diǎn)法(表層、0.2 H、0.4H、0.6 H、0.8 H、底層),其中表層為距表面0.5m、底層為距底0.5 m,H為瞬時(shí)水深。同時(shí)在2、3、4站采用RBR CTD進(jìn)行每小時(shí)1次的溫鹽剖面觀測(cè)。

    在近岸設(shè)立4個(gè)臨時(shí)水尺進(jìn)行26 h的短期同步潮位觀測(cè),同時(shí)2號(hào)站隨船放置自動(dòng)潮位記錄儀進(jìn)行水位測(cè)量。觀測(cè)期間基本為靜風(fēng)條件,風(fēng)對(duì)水動(dòng)力的影響很小。觀測(cè)期間上游龍?zhí)了畨翁幍膹搅髁繛?0m3/s左右。

    此外,海南省海洋開發(fā)規(guī)劃設(shè)計(jì)研究院還于2009年11月20日至2009年12月28日在站位3附近利用漁民搭建的捕魚平臺(tái),采用自容式溫鹽儀進(jìn)行了連續(xù)38天的表層溫鹽觀測(cè)。

    通過(guò)分析上述實(shí)測(cè)數(shù)據(jù),本文主要研究南渡江干流(包括2至4號(hào)站)在枯季大小潮期間的流速、鹽度、水體分層等的潮內(nèi)變化。對(duì)于流速數(shù)據(jù),根據(jù)一個(gè)潮周期內(nèi)垂向平均流速的變化確定主流向(Emery etal,2001)。由于河道較窄,水流基本為往復(fù)流,本研究中的流速分析以主流向上的流速變化為主。水體分層采用分層系數(shù)(Sp=δs/s0,其中δs為表底層的鹽度差,s0為垂向平均鹽度)加以量化。結(jié)合38天的連續(xù)溫鹽觀測(cè),采用EFDC(Environmental Fluid Dynamics Code)(Hamrick et al,1992)模型模擬了觀測(cè)期間觀測(cè)點(diǎn)的流速過(guò)程。該模型在南渡江河口已進(jìn)行過(guò)良好的水位與流速驗(yàn)證。

    2 結(jié)果

    2.1 小潮期的流速、鹽度、水體分層的時(shí)空變化

    這里給出干流水道口門站的水位過(guò)程。可見,在小潮期一個(gè)太陰日內(nèi),南渡江河口的水位經(jīng)歷了一次大的落-漲-落潮過(guò)程,但在大的漲落潮過(guò)程中,又發(fā)育有次一級(jí)的漲落潮過(guò)程,如在圖中26小時(shí)附近,有一次一級(jí)的落-漲過(guò)程。

    圖3 口門站小潮期的水位過(guò)程

    從站2的水位與流速過(guò)程的對(duì)比來(lái)看,流速變化幾乎滯后于水位變化90°,顯示出河口以駐波為主的潮波特征。由于南渡江的河口長(zhǎng)度較短,遠(yuǎn)小于主要分潮波(M2,S2,K1,O1)的1/4波長(zhǎng),駐波系統(tǒng)發(fā)育。從漲落潮歷時(shí)來(lái)看,漲潮歷時(shí)約7 h,而落潮歷時(shí)長(zhǎng)達(dá)18 h,漲潮歷時(shí)遠(yuǎn)短于落潮歷時(shí),與此相對(duì)應(yīng),漲潮最大流速與平均流速都大于落潮流速,表現(xiàn)出漲潮優(yōu)勢(shì)的流速不對(duì)稱類型。

    圖4 站2小潮期的流速、水體分層與鹽度過(guò)程

    與流速過(guò)程相對(duì)應(yīng),站2的鹽度表現(xiàn)出相應(yīng)的變化。整個(gè)觀測(cè)期內(nèi),水體的分層系數(shù)介于0至1之間,以部分混合的狀態(tài)為主(Prandle et al,1985)。落潮時(shí)鹽度不斷減小,最小鹽度出現(xiàn)在落憩時(shí)。與部分混合河口的一般規(guī)律不同(Dyer et al,1997),站2的水體分層系數(shù)在落急時(shí)較大,落憩時(shí)達(dá)到極小,符合鹽水楔河口的特征(Geyer etal,1989)。漲潮過(guò)程中,底部鹽度率先增大,水體分層增強(qiáng),水體中上層鹽度隨后增大,到漲急附近時(shí)水體混合充分,分層系數(shù)達(dá)到0。這種狀態(tài)一直持續(xù)到漲憩。在隨后的落潮過(guò)程中,表層鹽度先于底層鹽度減小,水體分層系數(shù)不斷增大??傮w而言,由于站2位于河口口門,潮流流速大,水體呈現(xiàn)出部分混合到充分混合的狀態(tài),落潮時(shí)為部分混合,漲潮時(shí)變?yōu)槌浞只旌?。這是由于潮汐應(yīng)變(Tidal straining)效應(yīng) (Simpson etal,1990)的影響,落潮時(shí)層化加強(qiáng),漲潮時(shí)混合加強(qiáng)。

    與站2相比,站3的流速明顯減小,其表層自向海流轉(zhuǎn)為向陸流的轉(zhuǎn)流時(shí)刻延遲,而底層的轉(zhuǎn)流時(shí)刻提前。與站2不同,在次一級(jí)的落-漲過(guò)程中(28-30 h),表層流仍為向海流,但中層與底層流則變?yōu)橄蜿懥?。可見,?在落潮期的垂向流速剪切要大于站2,而在漲潮期流速的垂向剪切較小。

    與流速較弱相對(duì)應(yīng),站3的水體分層系數(shù)在第一個(gè)落潮期內(nèi)幾乎都大于1,呈現(xiàn)出高度分層的狀態(tài)(Prandle etal,1985)。在落潮過(guò)程中,表層鹽度有所降低,而底層鹽度則未見減小,特別是在13-19 h之間,底層鹽度基本未變,這與該時(shí)段內(nèi)的底層流速近乎為0相對(duì)應(yīng)。站3的水深約為6 m,為人為挖沙形成的深坑(圖2),落潮時(shí)水體分層強(qiáng),抑制了表層與底層的動(dòng)量交換,流速的垂向剪切較大,底層幾乎為靜水。值得注意的是,在落潮后期,站3底層的鹽度大于站2,表明存在局部的底層鹽度逆梯度,這是由于挖沙坑對(duì)底部高鹽水的捕集作用所致。在漲潮期,站3的底層鹽度減小,表層鹽度增大,水體分層減弱,在漲急時(shí)水體的分層系數(shù)達(dá)到最小。我們推測(cè):漲急時(shí)刻在站3處可能發(fā)生了水躍(jump),產(chǎn)生類似于sill(檻)的動(dòng)力過(guò)程,垂向流速增大,水體中的動(dòng)能從上部傳到下部(Chao etal,1991),導(dǎo)致垂向混合作用顯著加強(qiáng)。在隨后的漲潮階段,鹽度增加,而水體分層系數(shù)則經(jīng)歷了從增加到減小的過(guò)程,到漲憩附近分層系數(shù)達(dá)到極小。漲潮過(guò)程中鹽躍層(垂向鹽度梯度最大處)有不斷向上抬升的現(xiàn)象,表明漲潮過(guò)程中底邊界層不斷增厚,這與大多數(shù)河口的底邊界層發(fā)育規(guī)律一致(Ralston etal,2010)。21-23 h之間表層鹽度減小明顯,是由于局部短時(shí)降雨所致。在隨后的落潮過(guò)程中,水體鹽度減小,特別是表層鹽度。與此同時(shí),水體分層則不斷增強(qiáng),這與部分混合河口的典型過(guò)程一致。可見,由于站3位于河口的中段,水體呈現(xiàn)出部分混合到高度分層的狀態(tài),以部分混合的狀態(tài)為主。

    圖5 站3小潮期的流速、水體分層與鹽度過(guò)程

    圖6 站4小潮期的流速、水體分層與鹽度過(guò)程

    再向上游,到站4時(shí),流速進(jìn)一步減小,表、中、底層的流速相位不同步。值得注意的是,水體中的底層流在大多數(shù)時(shí)刻都表現(xiàn)為向陸流,表明在站4的底層,落潮期由鹽度梯度產(chǎn)生的斜壓力要大于由水位梯度產(chǎn)生的正壓力。站4也位于一深坑位置(圖2),在落潮期水體中上部的鹽度減小,而底層鹽度幾乎不變,表現(xiàn)在水體分層上,即分層系數(shù)不斷增大,到落憩時(shí)達(dá)到最大。漲潮時(shí)由于潮汐應(yīng)變效應(yīng)以及由于地形所產(chǎn)生的動(dòng)力過(guò)程,高鹽水覆蓋在低鹽水之上,導(dǎo)致水體不穩(wěn)定,垂向混合加強(qiáng)。該時(shí)刻底層鹽度減小,表層鹽度增大,水體的分層系數(shù)減小,到漲急時(shí)分層最小。在隨后的落潮過(guò)程中,表層鹽度減小,而底層鹽度則增大,與之相對(duì)應(yīng),水體的分層系數(shù)不斷增大。落潮時(shí)底層鹽度增大的原因在于該時(shí)刻底層流為向陸流,高鹽水自海向陸輸運(yùn)??傮w來(lái)看,站4的水體分層系數(shù)大,處于高度分層狀態(tài)。

    2.2 大潮期的流速、鹽度、水體分層的時(shí)空變化

    在大潮期,水位過(guò)程表現(xiàn)為兩漲兩落,其中后一次的落潮過(guò)程中水位變化的幅度較小。

    與小潮期不同,在大潮期的落潮流速要大于漲潮流速,流速不對(duì)稱表現(xiàn)為落潮優(yōu)勢(shì)型。河口區(qū)的流速不對(duì)稱受多種因素控制,F(xiàn)riedrichs與Madsen(1992)對(duì)美國(guó)河口的流速不對(duì)稱進(jìn)行了綜合研究,提出了判據(jù)其中b0為斷面的平均寬度,Δb為斷面高水位的寬度與低水位時(shí)的斷面寬度差。該判據(jù)反映的是潮波傳播過(guò)程中的底磨擦效應(yīng)與潮灘地形效應(yīng)的對(duì)比情況,當(dāng)r>0時(shí),底磨擦效應(yīng)強(qiáng),漲潮優(yōu)勢(shì)發(fā)育,反之則落潮優(yōu)勢(shì)發(fā)育。但在日潮為主的河口區(qū),Ranasinghe與Pat-tiaratchi(2000)研究發(fā)現(xiàn)其流速不對(duì)稱主要受外海潮波的漲落潮歷時(shí)不等所控制。對(duì)于南渡江河口,在大潮時(shí)r顯然要大于小潮時(shí),根據(jù)Friedrichs與Madsen(1992)的判據(jù)應(yīng)為漲潮優(yōu)勢(shì)更發(fā)育,而結(jié)果則為大潮時(shí)落潮優(yōu)勢(shì)更明顯,表明南渡江河口的流速不對(duì)稱主要受外海潮波的漲落潮歷時(shí)不等所控制。大潮時(shí)落潮歷時(shí)要短于漲潮歷時(shí),而小潮時(shí)落潮歷時(shí)長(zhǎng)于漲潮歷時(shí)。

    圖7 口門站大潮期的水位過(guò)程

    與小潮時(shí)相比,大潮時(shí)站2的水體分層系數(shù)大大減小,水體分層系數(shù)在0-0.5之間變動(dòng),表明大潮期的混合作用更強(qiáng),這與世界上許多河口相類似(Haasetal,1977)。在落潮期表層鹽度減小,水體分層增強(qiáng),在落急時(shí)分層系數(shù)達(dá)到極大。表明落潮時(shí)潮流的混合作用小于由于潮汐應(yīng)變導(dǎo)致分層加強(qiáng)的效應(yīng)。在落潮后期,鹽度進(jìn)一步減小,分層系數(shù)也有所減小,但在落憩附近分層又增大。漲潮時(shí)鹽度增高,水體混合加強(qiáng),呈現(xiàn)為充分混合狀態(tài)。在隨后的落潮期鹽度又不斷減小,水體分層也隨之增強(qiáng)??傮w來(lái)看,大潮期站2為部分混合至充分混合狀態(tài),潮汐應(yīng)變效應(yīng)顯著,水體分層落潮期遠(yuǎn)大于漲潮期。

    圖8 站2大潮期的流速、水體分層與鹽度過(guò)程

    圖9 站3大潮期的流速、水體分層與鹽度過(guò)程

    在站3,大潮期流速顯著增大,底層流的漲潮歷時(shí)要長(zhǎng)于表層與中層。水體分層系數(shù)(0.2~1.5之間)比小潮期減小,混合作用加強(qiáng)。落潮時(shí)鹽度減小,水體分層增強(qiáng),鹽躍層在垂向上不斷下移,表明底邊界層受分層作用的抑制不斷變薄。這種現(xiàn)象在小潮期的落潮過(guò)程中不明顯。漲潮期水體垂向混合加強(qiáng),分層系數(shù)減小,在漲憩時(shí)鹽度達(dá)到最大,分層系數(shù)達(dá)到最小。在隨后的落潮過(guò)程中,鹽度減小,分層增強(qiáng),呈現(xiàn)出部分混合型河口的一般特征。

    圖10 站4大潮期的流速、水體分層與鹽度過(guò)程

    大潮期站4的表層流以向海流動(dòng)為主,而底層流以向陸流動(dòng)為主。在小的落潮過(guò)程中(26 h附近),底層出現(xiàn)向海流的時(shí)段很短。在漲潮期,中層流的流速大于表層(21-24 h),呈現(xiàn)出典型的漲潮期河口垂向流速分布中次表層流速最大的特征(Unclesetal,2002)。在落潮期,水體上部的鹽度減小劇烈,且鹽躍層不斷下移,水體的分層系數(shù)在落潮期不斷增大。漲潮時(shí)底層鹽度減小,表層鹽度增大,水體混合加強(qiáng),分層減弱。在隨后的次一級(jí)的落-漲過(guò)程中,鹽度呈現(xiàn)出先減后增的過(guò)程,水體分層則先增大后減小。

    2.3 連續(xù)38天觀測(cè)所揭示的鹽度變化規(guī)律

    連續(xù)38天觀測(cè)站觀測(cè)期只進(jìn)行了溫鹽觀測(cè),為分析觀測(cè)期鹽度隨大小潮的變化情況,采用EFDC模型模擬了觀測(cè)期間觀測(cè)點(diǎn)的水位與流速變化,將模型計(jì)算的流速(這里給出的是垂向平均流速)與觀測(cè)得到的鹽度進(jìn)行綜合分析。

    從圖中可以看出,觀測(cè)站在大潮期的鹽度變化大,高鹽度值基本出現(xiàn)于漲憩階段,而低鹽度值出現(xiàn)于漲急之前。從一個(gè)潮周期內(nèi)的鹽度變化幅度來(lái)看,觀測(cè)站大潮期的最大鹽度要大于小潮期,最小鹽度則要小于小潮期。從圖11b可以看出,大潮期落潮時(shí)鹽度減小緩慢,而在漲潮期的漲急后鹽度迅速增大,這與落潮期流速較小,而漲潮期流速較大有關(guān),另一方面,漲急時(shí)刻由于流速大,混合作用強(qiáng),也是表層鹽度迅速增大的原因。這與站3(圖9)在大潮期表層鹽度的變化相一致。在小潮期鹽度變化與流速的相關(guān)性不是很強(qiáng),但基本上呈現(xiàn)出漲潮時(shí)鹽度增大,落潮時(shí)鹽度減小,鹽度的極大值出現(xiàn)于漲急至漲憩期間的規(guī)律。

    圖11 連續(xù)38天觀測(cè)站的軸向流速與鹽度變化過(guò)程

    3 討論與結(jié)論

    3.1 枯季時(shí)鹽水入侵到底是大潮時(shí)強(qiáng)還是小潮時(shí)強(qiáng)

    趙軍鵬 (2011)發(fā)現(xiàn):根據(jù) Lerczak等(2006)對(duì)鹽分通量的分解方法,將FVCOM模型計(jì)算結(jié)果進(jìn)行斷面鹽分輸運(yùn)通量分析,分解成平流項(xiàng) (Advection)、潮汐振蕩項(xiàng) (tidal oscillatory)及平均流剪切項(xiàng)(steady shear)。平流項(xiàng)為由徑流作用及潮汐的Stokes效應(yīng)產(chǎn)生的鹽分輸運(yùn),以向海為主,平均流剪切項(xiàng)主要由河口環(huán)流產(chǎn)生,而潮汐振蕩項(xiàng)則由流速與鹽度的潮內(nèi)(intratidal)變化的相互作用產(chǎn)生,包括潮泵(tidal pumping)與潮汐捕集(tidal trapping)效應(yīng)。平均流剪切與潮汐振蕩是鹽分向陸輸運(yùn)的主要?jiǎng)恿C(jī)制。在南渡江干流,自海向陸,隨著潮汐動(dòng)力的不斷減弱,徑流動(dòng)力的不斷增強(qiáng),鹽分向陸輸運(yùn)的貢獻(xiàn)項(xiàng)中,潮汐振蕩項(xiàng)不斷減小,而由重力環(huán)流所造成的平均流剪切項(xiàng)則不斷增大。在口門附近,鹽分向陸輸運(yùn)的主要貢獻(xiàn)項(xiàng)為潮汐振蕩項(xiàng),在站3附近(也為38天連續(xù)觀測(cè)點(diǎn)所在位置)的斷面上,鹽分向陸輸運(yùn)中潮汐振蕩項(xiàng)仍占相當(dāng)重要的地位。一般而言,潮汐振蕩項(xiàng)與潮流速度的大小成正相關(guān)(Lerczak et al,2006),潮流速度越大,由潮汐振蕩項(xiàng)所形成的鹽分向陸輸運(yùn)通量越大,這是造成38天連續(xù)觀測(cè)站中最大鹽度在大潮期要大于小潮期的動(dòng)力機(jī)制所在。再向陸,到站4所在的斷面時(shí),從圖10和圖6的對(duì)比可以看出,站4大潮期的鹽度要小于小潮期。這是由于在站4附近所在斷面,鹽分向陸輸運(yùn)主要通過(guò)平均流剪切項(xiàng)來(lái)完成,而平均流剪切項(xiàng)又主要由重力環(huán)流所造成,重力環(huán)流在小潮期要強(qiáng)于大潮期(Haasetal,1977)。如果以某一量值(如1ppt)的底層鹽度向陸延伸的距離作為鹽水入侵長(zhǎng)度的定義,可以斷定,南渡江鹽水入侵長(zhǎng)度在小潮時(shí)要大于大潮時(shí)。這與部分混合的河口(如珠江口磨刀門 (Gong et al,2011);美國(guó) Hudson河口(Lerczak etal,2006))的規(guī)律相一致,但與充分混和河口(如長(zhǎng)江口)的變化規(guī)律不同(Wu etal,2010)。對(duì)于充分混合的河口,鹽水入侵的長(zhǎng)度與潮流流速振幅的一次方至二次方成正比(Mac-Cready etal,2007)。

    3.2 人工挖沙所形成的深坑對(duì)鹽水入侵產(chǎn)生怎樣的影響

    從對(duì)站3在小潮期鹽度變化的分析來(lái)看,挖沙形成的深坑對(duì)高鹽水具有捕集作用。漲潮時(shí)(特別是大潮期)通過(guò)Sill動(dòng)力過(guò)程進(jìn)入深坑的高鹽水,在落潮期由于水體分層,垂向混合作用小,再加上深坑處的水平流速近于0,不能被向海輸運(yùn)而富集于深坑處。這些富集的高鹽水在下一次的漲潮過(guò)程中被卷吸到水體上部隨漲潮流向陸輸運(yùn),從而在一次漲落潮過(guò)程中,深坑所在處落潮時(shí)鹽度低,漲潮時(shí)鹽度高,類似于潮泵(tidal pumping)過(guò)程,或Fisher等 (Fisher etal,1979;Geyer etal,1992;)所描述的河口邊灘小灣對(duì)高鹽水的捕集效應(yīng),成為鹽分向陸輸運(yùn)的一個(gè)重要機(jī)制,對(duì)鹽水入侵產(chǎn)生加劇作用。這一效應(yīng)可能是南渡江近年來(lái)鹽水入侵加劇的一個(gè)重要?jiǎng)恿σ蛩?,有待今后進(jìn)一步地定量研究。

    本研究的主要結(jié)論如下:

    (1)南渡江河口的潮波為駐波,流速不對(duì)稱主要受外海潮波的漲落潮歷時(shí)不等所控制,在2009年2月觀測(cè)的大潮期表現(xiàn)為落潮優(yōu)勢(shì),小潮期表現(xiàn)為漲潮優(yōu)勢(shì)。

    (2)在小潮期水體分層自海向陸不斷增強(qiáng),站2為充分混合至部分混合狀態(tài),站3為部分混合至高度分層狀態(tài),而站4則以高度分層狀態(tài)為主。各站的鹽度在落潮時(shí)減小,漲潮時(shí)增大,而其分層系數(shù)則在落潮時(shí)增大,在漲潮時(shí)減小。站3與站4的底層鹽度在落潮時(shí)變化很小。

    (3)大潮期河口的混合作用加強(qiáng),各站的分層系數(shù)都相應(yīng)減小。各站在落潮時(shí)鹽度減小,水體分層增強(qiáng),落憩時(shí)分層最強(qiáng)。漲潮時(shí)鹽度增大,水體充分混合。站3和站4在落潮期鹽躍層不斷下移。漲急時(shí)混合增強(qiáng),其后鹽躍層上移,鹽度增大。

    (4)結(jié)合數(shù)值模擬的流速數(shù)據(jù),對(duì)38天觀測(cè)站鹽度數(shù)據(jù)的分析表明,觀測(cè)站的表層鹽度在大潮期的變化幅度大,最大鹽度高于小潮期,最小鹽度低于小潮期。鹽度降低的過(guò)程緩慢,而增加的過(guò)程迅速,這與河口流速不對(duì)稱以及鹽度變化的動(dòng)力機(jī)制密切相關(guān)。

    (5)南渡江枯季時(shí)鹽水入侵長(zhǎng)度在小潮期要大于大潮期。這與其鹽分輸運(yùn)的動(dòng)力機(jī)制有關(guān)。重力環(huán)流所造成的平均流剪切項(xiàng)是鹽分向陸輸運(yùn)的重要貢獻(xiàn)項(xiàng),該項(xiàng)在河口的向陸部分小潮期要強(qiáng)于大潮期。

    (6)人工挖沙形成的深坑對(duì)漲潮帶入的高鹽水具有捕集作用,深坑處一個(gè)潮周期內(nèi)平均鹽分向陸輸運(yùn),可能是近年來(lái)南渡江河口鹽水入侵加劇的重要?jiǎng)恿C(jī)制之一。

    Chao S Y,Paluszkiewicz T,1991.The hydraulics ofdensity currents over estuarine sills.Journal of Geophysical Research,96 (C4):7065-7076.

    Dyer K R,1997.Estuaries:a physical introduction.2nd Edition,John Wiley&Sons Ltd.,195.

    Emery W J,Thomas R E,2001.Data analysis methods in physical oceanography.Elsevier Science,658.

    Fisher H B,ListE J,Koh R C Y,etal,1979.Mixing in inland and coastal waters.Academic Press,483.

    Friedrichs C T,Madsen O S,1992.Non-lineardiffusion ofthe tidalsignal in frictionally dominated embayments.Journal of Geophysical Research,97:5637-5650.

    Geyer W R,F(xiàn)armer D M,1989.Tide‐induced variation ofthe dynamics of a salt wedge estuary.Journal of Physical Oceanography,19(8):1060-1072.

    Geyer W R,Signell R P,1992.A reassessmentofthe role oftidaldispersion in estuaries and bays,Estuaries,15(2):97-108.doi:10.2307/1352684.

    Gong W,Shen J,2011.The response of salt intrusion to changes in river discharge and tidal mixing during the dry season in the Modaomen Estuary,China.Continental ShelfResearch,31:769-788.

    Haas L W,1977.The effect of spring-neap tidal cycle on the vertical salinity structure ofthe James,York and Rappahannock Rivers,Virginia,USA.Estuarine,CoastalMarine Science,5:485-496.

    Hamrick J M,1992.A Three-Dimensional Environmental Fluid Dynamics Computer Code:Theoretical and Computational Aspects.Special Reportin Applied Marine Science and Ocean Engineering.No.317 College ofWilliam and Mary,VIMS,63.

    Lerczak J A,Geyer W R,Chant R J,2006.Mechanisms Driving the Time-Dependent Salt Flux in a Partially Stratified Estuary.Journal ofPhysical Oceanography,36:2296-2311.

    Mac Creacly P,2007.Estuarine adjustment.Journal of Physical O-ceanography,37:2133-2145.

    Prandle D,1985.On salinity regimes and the verticalstructure ofresidual flows in narrow tidalestuaries.Estuarine,Coastaland Shelf Science,20:615-635.

    Ralston D K,Geyer W R,Lerczak J A,2010.Structure,variability,and saltflux in a strongly forced saltwedge estuary.JournalofGeophysicalResearch,115,C06005.doi:10.1029/2009JC005806.

    Ranasinghe R,Pattiaratchi C,2000.Tidalinletvelocity asymmetry in diurnalregimes.ContinentalShelfResearch,29(17):2347-2366.

    Simpson J H,Brown J,Matthews J,et al,1990.Tidal straining,density currents,and stirring in controlofestuarine stratification.Estuaries,13(2):125-132.

    Uncles R J,2002.Estuarine physical processes research:some recent studies and progress.Estuarine,Coastal and Shelf Science,55(6):829-856.

    Wu H,Zhu J R,2010.Advection scheme with 3rd high-order spatialinterpolation atthe middle temporalleveland its application to saltwater intrusion in the Changjiang Estuary.Ocean Modelling,33(1-2):33-51.

    趙軍鵬,2011.海南島南渡江鹽水入侵及其對(duì)人類活動(dòng)的影響.碩士學(xué)位論文.廣州:中山大學(xué).

    猜你喜歡
    南渡江小潮落潮
    潮 汐
    希 望
    海南有條南渡江
    新教育(2018年17期)2018-02-23 11:35:10
    希望
    中資出海并購(gòu)“落潮”
    新一季流行色已發(fā)布?快來(lái)入手同色系數(shù)碼3C小潮物!
    為什么會(huì)有漲潮和落潮
    海南行駛在河海上的列車
    南渡江三角洲地貌演變研究分析
    生態(tài)補(bǔ)償機(jī)制實(shí)驗(yàn)帶的創(chuàng)建路徑初探
    亚洲四区av| 综合色丁香网| 一区二区三区精品91| 亚洲经典国产精华液单| 国产精品麻豆人妻色哟哟久久| 国产在视频线精品| a级毛片免费高清观看在线播放| 两个人的视频大全免费| 国产伦在线观看视频一区| 亚洲美女搞黄在线观看| 春色校园在线视频观看| 18禁在线播放成人免费| 日日撸夜夜添| 人妻夜夜爽99麻豆av| 三级经典国产精品| 亚洲精品一区蜜桃| 日韩不卡一区二区三区视频在线| 99热这里只有是精品在线观看| 麻豆成人午夜福利视频| 午夜福利在线观看免费完整高清在| 日本av手机在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 我的老师免费观看完整版| 中文字幕久久专区| 少妇人妻精品综合一区二区| 91在线精品国自产拍蜜月| 十分钟在线观看高清视频www | 欧美bdsm另类| 美女脱内裤让男人舔精品视频| 国产一区有黄有色的免费视频| 中文字幕免费在线视频6| 卡戴珊不雅视频在线播放| 国内揄拍国产精品人妻在线| 高清午夜精品一区二区三区| 亚洲人成网站在线播| 日韩强制内射视频| 欧美日韩视频精品一区| 看十八女毛片水多多多| 久久久久久久大尺度免费视频| 美女福利国产在线 | 午夜激情久久久久久久| 亚洲精品亚洲一区二区| 九草在线视频观看| 亚洲三级黄色毛片| 女性被躁到高潮视频| 激情五月婷婷亚洲| 亚洲av.av天堂| 91久久精品电影网| 亚洲最大成人中文| 人人妻人人看人人澡| 国产一级毛片在线| 一区二区三区免费毛片| 亚洲国产精品一区三区| 亚洲第一区二区三区不卡| 久久亚洲国产成人精品v| 亚洲av不卡在线观看| 国产精品99久久99久久久不卡 | 欧美日韩视频精品一区| 亚洲精品亚洲一区二区| 少妇人妻精品综合一区二区| 又粗又硬又长又爽又黄的视频| 亚洲av欧美aⅴ国产| 最近中文字幕2019免费版| 高清视频免费观看一区二区| 久热这里只有精品99| 精品亚洲乱码少妇综合久久| 中文乱码字字幕精品一区二区三区| 美女xxoo啪啪120秒动态图| 日韩三级伦理在线观看| 亚洲精品日韩av片在线观看| 欧美日韩综合久久久久久| 人人妻人人爽人人添夜夜欢视频 | 视频区图区小说| 男人舔奶头视频| 制服丝袜香蕉在线| 亚洲经典国产精华液单| 国产精品.久久久| 欧美高清成人免费视频www| 精品一品国产午夜福利视频| 美女中出高潮动态图| 午夜免费观看性视频| 高清在线视频一区二区三区| 国产一级毛片在线| 色视频www国产| 欧美日韩国产mv在线观看视频 | 午夜免费观看性视频| tube8黄色片| 亚洲精品乱码久久久v下载方式| 天堂俺去俺来也www色官网| 久久人人爽人人片av| 国内精品宾馆在线| 丰满迷人的少妇在线观看| 亚洲精品色激情综合| av国产免费在线观看| 精品久久久噜噜| 亚洲人成网站在线播| 干丝袜人妻中文字幕| 如何舔出高潮| 国产精品伦人一区二区| 亚洲av成人精品一二三区| 寂寞人妻少妇视频99o| 1000部很黄的大片| 伦理电影大哥的女人| 午夜福利网站1000一区二区三区| 国产日韩欧美亚洲二区| 在线观看免费视频网站a站| 精品国产一区二区三区久久久樱花 | 免费大片黄手机在线观看| av免费观看日本| 国产精品女同一区二区软件| 亚洲成人一二三区av| 国产在线视频一区二区| 五月伊人婷婷丁香| av在线播放精品| 女的被弄到高潮叫床怎么办| 亚洲,欧美,日韩| 亚洲精品日韩av片在线观看| 亚洲欧美日韩卡通动漫| 亚洲美女搞黄在线观看| 欧美日韩综合久久久久久| 最近最新中文字幕免费大全7| 一区二区三区乱码不卡18| 中文字幕免费在线视频6| 99热这里只有是精品在线观看| 午夜激情福利司机影院| 老司机影院成人| 精品视频人人做人人爽| 久久av网站| 亚洲色图av天堂| 在线天堂最新版资源| 国产av精品麻豆| 成人美女网站在线观看视频| 精品视频人人做人人爽| 99国产精品免费福利视频| 五月开心婷婷网| 天堂中文最新版在线下载| 欧美日韩在线观看h| 亚洲国产精品专区欧美| 青春草亚洲视频在线观看| 一级毛片我不卡| 青春草国产在线视频| 老师上课跳d突然被开到最大视频| 欧美老熟妇乱子伦牲交| 看十八女毛片水多多多| 亚洲av欧美aⅴ国产| 激情 狠狠 欧美| 女的被弄到高潮叫床怎么办| 交换朋友夫妻互换小说| 国精品久久久久久国模美| 国产 一区 欧美 日韩| 在线看a的网站| 狂野欧美激情性xxxx在线观看| 欧美 日韩 精品 国产| 九九在线视频观看精品| 亚洲精品国产色婷婷电影| 欧美成人精品欧美一级黄| 人体艺术视频欧美日本| 男女国产视频网站| 中文在线观看免费www的网站| 午夜福利网站1000一区二区三区| 久久国产亚洲av麻豆专区| 久久久a久久爽久久v久久| 国产永久视频网站| 精品久久久久久久末码| 另类亚洲欧美激情| 一本—道久久a久久精品蜜桃钙片| 人妻少妇偷人精品九色| 女性生殖器流出的白浆| 日韩强制内射视频| 亚洲国产毛片av蜜桃av| 黄片无遮挡物在线观看| 我要看黄色一级片免费的| 午夜免费男女啪啪视频观看| 草草在线视频免费看| 在现免费观看毛片| 五月开心婷婷网| 又粗又硬又长又爽又黄的视频| 久久99热6这里只有精品| 成人亚洲欧美一区二区av| 精品午夜福利在线看| 免费黄色在线免费观看| 欧美区成人在线视频| 国产精品三级大全| 精品亚洲乱码少妇综合久久| 国产真实伦视频高清在线观看| 最后的刺客免费高清国语| 中文字幕久久专区| 亚洲欧美成人综合另类久久久| 中文乱码字字幕精品一区二区三区| 熟女av电影| 日韩制服骚丝袜av| 色视频www国产| 国产高清不卡午夜福利| 精品少妇久久久久久888优播| 国产日韩欧美亚洲二区| 日日啪夜夜撸| 99久国产av精品国产电影| 人人妻人人看人人澡| 久久久久精品久久久久真实原创| 国产淫片久久久久久久久| 色5月婷婷丁香| 国产成人免费观看mmmm| 日韩一区二区三区影片| 亚洲第一区二区三区不卡| 一级二级三级毛片免费看| 欧美成人午夜免费资源| 国产黄频视频在线观看| 夜夜骑夜夜射夜夜干| 18+在线观看网站| 亚洲欧洲国产日韩| 欧美精品亚洲一区二区| 黑丝袜美女国产一区| a级毛片免费高清观看在线播放| 亚洲人成网站在线播| 永久免费av网站大全| 午夜免费鲁丝| 国产高清三级在线| 日韩精品有码人妻一区| 久久99热6这里只有精品| 精品国产三级普通话版| 日韩免费高清中文字幕av| 国产成人一区二区在线| 美女国产视频在线观看| 久久精品夜色国产| 美女脱内裤让男人舔精品视频| 久久影院123| 亚洲av在线观看美女高潮| 亚洲精品乱码久久久久久按摩| 国产精品女同一区二区软件| 国产淫语在线视频| 我的女老师完整版在线观看| 91精品国产九色| 交换朋友夫妻互换小说| 一二三四中文在线观看免费高清| 亚洲欧美成人精品一区二区| 亚洲自偷自拍三级| videos熟女内射| 欧美日韩视频精品一区| 人妻系列 视频| 日本-黄色视频高清免费观看| 我要看日韩黄色一级片| 青春草视频在线免费观看| 少妇被粗大猛烈的视频| 欧美xxxx性猛交bbbb| 免费av不卡在线播放| 成人漫画全彩无遮挡| 国产成人精品一,二区| 精品人妻视频免费看| 国产老妇伦熟女老妇高清| 成人免费观看视频高清| 亚洲精品久久午夜乱码| 欧美高清成人免费视频www| 亚洲精品乱久久久久久| 高清在线视频一区二区三区| 99热国产这里只有精品6| 日韩 亚洲 欧美在线| 久久国产精品大桥未久av | 亚洲精品成人av观看孕妇| 插逼视频在线观看| 美女内射精品一级片tv| 国产一区二区三区综合在线观看 | 国产在线视频一区二区| 国产视频首页在线观看| 日韩一区二区视频免费看| 美女国产视频在线观看| 精品99又大又爽又粗少妇毛片| 久久精品人妻少妇| 99久国产av精品国产电影| 伦理电影大哥的女人| 在线观看国产h片| tube8黄色片| 日韩伦理黄色片| 三级国产精品片| 日本免费在线观看一区| 亚洲色图综合在线观看| 亚洲性久久影院| 国产伦精品一区二区三区视频9| 亚洲av.av天堂| 久久av网站| 男女免费视频国产| 日韩精品有码人妻一区| 日本av手机在线免费观看| 欧美极品一区二区三区四区| 国产精品.久久久| 日韩制服骚丝袜av| 一区二区三区免费毛片| 国产在线一区二区三区精| 欧美最新免费一区二区三区| 最新中文字幕久久久久| 天美传媒精品一区二区| 妹子高潮喷水视频| 秋霞伦理黄片| 男女啪啪激烈高潮av片| 一边亲一边摸免费视频| 免费观看的影片在线观看| 亚洲av成人精品一区久久| 交换朋友夫妻互换小说| 夫妻午夜视频| 国产 精品1| 伦精品一区二区三区| 亚洲国产av新网站| 国模一区二区三区四区视频| 色5月婷婷丁香| 欧美国产精品一级二级三级 | 国产av一区二区精品久久 | 国产深夜福利视频在线观看| 日韩伦理黄色片| 啦啦啦啦在线视频资源| 日本与韩国留学比较| 狂野欧美激情性bbbbbb| 日韩av在线免费看完整版不卡| 美女cb高潮喷水在线观看| 网址你懂的国产日韩在线| 国产黄片视频在线免费观看| 久久av网站| 丰满乱子伦码专区| 精品视频人人做人人爽| 有码 亚洲区| 国产色爽女视频免费观看| 菩萨蛮人人尽说江南好唐韦庄| 热re99久久精品国产66热6| 久久久亚洲精品成人影院| av不卡在线播放| 女人久久www免费人成看片| 欧美日韩视频高清一区二区三区二| 欧美97在线视频| 少妇人妻久久综合中文| 中文字幕人妻熟人妻熟丝袜美| 成人国产麻豆网| 晚上一个人看的免费电影| 欧美精品一区二区免费开放| 性色av一级| 天天躁夜夜躁狠狠久久av| 国内精品宾馆在线| 精品国产露脸久久av麻豆| 亚洲成人手机| 黑丝袜美女国产一区| 国产精品精品国产色婷婷| 夜夜爽夜夜爽视频| 亚洲国产av新网站| 久热久热在线精品观看| 五月伊人婷婷丁香| 亚洲人成网站在线播| 国产 一区精品| 国产精品一二三区在线看| 中文字幕av成人在线电影| 久久国内精品自在自线图片| 免费观看av网站的网址| 国模一区二区三区四区视频| 欧美日韩综合久久久久久| 王馨瑶露胸无遮挡在线观看| 国产有黄有色有爽视频| 亚洲四区av| 亚洲国产高清在线一区二区三| 日韩欧美精品免费久久| 日本免费在线观看一区| 老司机影院成人| 亚洲欧美一区二区三区黑人 | 大又大粗又爽又黄少妇毛片口| av又黄又爽大尺度在线免费看| 国产一区有黄有色的免费视频| 18禁裸乳无遮挡动漫免费视频| 中文在线观看免费www的网站| 亚洲av男天堂| 久久精品国产鲁丝片午夜精品| 久久99热这里只有精品18| 国产91av在线免费观看| 亚洲美女视频黄频| 观看免费一级毛片| 高清av免费在线| 日韩电影二区| 97在线视频观看| 国产伦精品一区二区三区四那| 男女边吃奶边做爰视频| 少妇的逼水好多| 2018国产大陆天天弄谢| 国产精品麻豆人妻色哟哟久久| 十分钟在线观看高清视频www | 亚洲av不卡在线观看| 精品久久久久久久久亚洲| 亚洲国产成人一精品久久久| 最近的中文字幕免费完整| 日韩人妻高清精品专区| 天堂俺去俺来也www色官网| 自拍偷自拍亚洲精品老妇| 国产69精品久久久久777片| 亚洲欧美日韩无卡精品| 午夜精品国产一区二区电影| 亚洲国产精品专区欧美| 夜夜爽夜夜爽视频| 国产视频首页在线观看| 不卡视频在线观看欧美| 在线观看免费高清a一片| 成人亚洲精品一区在线观看 | 80岁老熟妇乱子伦牲交| 国产精品精品国产色婷婷| 亚洲久久久国产精品| 超碰av人人做人人爽久久| 欧美+日韩+精品| 日本爱情动作片www.在线观看| 欧美激情极品国产一区二区三区 | 免费观看在线日韩| 亚洲aⅴ乱码一区二区在线播放| 在线精品无人区一区二区三 | 肉色欧美久久久久久久蜜桃| av女优亚洲男人天堂| 91午夜精品亚洲一区二区三区| 免费在线观看成人毛片| 观看av在线不卡| 晚上一个人看的免费电影| 看非洲黑人一级黄片| 特大巨黑吊av在线直播| 99久久精品热视频| 在线观看一区二区三区激情| 国产精品三级大全| 欧美+日韩+精品| 91久久精品国产一区二区成人| 国产白丝娇喘喷水9色精品| av国产久精品久网站免费入址| av.在线天堂| 老司机影院成人| 亚洲精品乱码久久久久久按摩| 性色av一级| 精品国产乱码久久久久久小说| 伊人久久国产一区二区| 日韩大片免费观看网站| 午夜视频国产福利| 视频中文字幕在线观看| 久久久久久久久久人人人人人人| 男人爽女人下面视频在线观看| 亚洲精品亚洲一区二区| 特大巨黑吊av在线直播| 最近最新中文字幕免费大全7| 国产高清三级在线| 精品久久久久久久久av| 91狼人影院| 日韩成人av中文字幕在线观看| 久久久久久久久久久免费av| 能在线免费看毛片的网站| 久久国产精品大桥未久av | 亚洲国产精品成人久久小说| 亚洲av欧美aⅴ国产| 国产精品一二三区在线看| 国产一区二区在线观看日韩| 老女人水多毛片| 精品久久久久久电影网| 亚洲av不卡在线观看| 一区在线观看完整版| 直男gayav资源| 最新中文字幕久久久久| 亚洲国产成人一精品久久久| 一区在线观看完整版| 久久精品久久久久久久性| 在线观看av片永久免费下载| av不卡在线播放| 十八禁网站网址无遮挡 | 欧美精品国产亚洲| 日韩,欧美,国产一区二区三区| 高清av免费在线| 久久久a久久爽久久v久久| 久久久久精品久久久久真实原创| 美女高潮的动态| 纵有疾风起免费观看全集完整版| 日韩欧美精品免费久久| h日本视频在线播放| 人人妻人人添人人爽欧美一区卜 | 亚洲av电影在线观看一区二区三区| 亚洲久久久国产精品| 99久久综合免费| 欧美xxⅹ黑人| 美女高潮的动态| 国产一区二区在线观看日韩| 亚洲精品456在线播放app| 免费人妻精品一区二区三区视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美清纯卡通| 人人妻人人澡人人爽人人夜夜| 超碰av人人做人人爽久久| 日日摸夜夜添夜夜添av毛片| 久久久久久久久大av| 色视频www国产| 人妻夜夜爽99麻豆av| 日韩欧美精品免费久久| 岛国毛片在线播放| av卡一久久| 我的老师免费观看完整版| 久久精品久久久久久久性| 视频区图区小说| 精品一区二区三卡| 一级毛片久久久久久久久女| 成人免费观看视频高清| 高清欧美精品videossex| 午夜福利在线在线| 伊人久久国产一区二区| 激情五月婷婷亚洲| h日本视频在线播放| 在线免费十八禁| 免费高清在线观看视频在线观看| 免费黄网站久久成人精品| 大片电影免费在线观看免费| 久久久久久久久久成人| 夜夜骑夜夜射夜夜干| 蜜臀久久99精品久久宅男| 欧美成人精品欧美一级黄| 我的女老师完整版在线观看| 欧美老熟妇乱子伦牲交| 国产高清有码在线观看视频| 欧美xxⅹ黑人| 性色av一级| 成人国产av品久久久| 国产亚洲最大av| 黄色一级大片看看| 国产人妻一区二区三区在| 精品久久久久久久久亚洲| 久久国产精品男人的天堂亚洲 | 精品国产一区二区三区久久久樱花 | 日本欧美视频一区| 一级毛片久久久久久久久女| 日本免费在线观看一区| 免费看av在线观看网站| 十分钟在线观看高清视频www | 久久久久久久久久人人人人人人| 97精品久久久久久久久久精品| 在线观看一区二区三区| 久久这里有精品视频免费| 高清不卡的av网站| 国产高潮美女av| 国产一级毛片在线| 又粗又硬又长又爽又黄的视频| 免费看光身美女| 毛片女人毛片| 免费人妻精品一区二区三区视频| 国产 一区 欧美 日韩| 一级二级三级毛片免费看| 国产熟女欧美一区二区| 国产精品一二三区在线看| 十分钟在线观看高清视频www | 美女脱内裤让男人舔精品视频| 99热全是精品| 十分钟在线观看高清视频www | 国产熟女欧美一区二区| 国产亚洲一区二区精品| 亚洲精品色激情综合| 亚洲在久久综合| 国产精品嫩草影院av在线观看| 亚洲天堂av无毛| 免费少妇av软件| 久久精品久久精品一区二区三区| 边亲边吃奶的免费视频| 内地一区二区视频在线| 国产精品爽爽va在线观看网站| 成人亚洲欧美一区二区av| 三级国产精品片| 一本色道久久久久久精品综合| 国产精品99久久久久久久久| 久久人人爽人人爽人人片va| 久久鲁丝午夜福利片| 国产精品久久久久久精品电影小说 | 国内精品宾馆在线| 五月伊人婷婷丁香| 人妻 亚洲 视频| 男男h啪啪无遮挡| 久久久午夜欧美精品| 观看免费一级毛片| 久久99热这里只频精品6学生| 大话2 男鬼变身卡| 18禁动态无遮挡网站| 日韩中文字幕视频在线看片 | 99久国产av精品国产电影| 亚洲成人中文字幕在线播放| 夜夜爽夜夜爽视频| 国产成人a∨麻豆精品| 男女边摸边吃奶| videos熟女内射| 久久久久久久精品精品| 欧美精品一区二区大全| 内地一区二区视频在线| 中国国产av一级| 制服丝袜香蕉在线| 成年av动漫网址| 在线观看国产h片| 国产免费福利视频在线观看| 99热这里只有精品一区| av一本久久久久| 欧美日韩一区二区视频在线观看视频在线| 免费看日本二区| 多毛熟女@视频| 成人亚洲欧美一区二区av| 亚洲成人一二三区av| 国产男女内射视频| 黄色日韩在线| 妹子高潮喷水视频| 一级a做视频免费观看| 99精国产麻豆久久婷婷| 美女视频免费永久观看网站| 国产亚洲午夜精品一区二区久久| 最近中文字幕2019免费版| 91精品国产九色| 精品一区在线观看国产| 高清午夜精品一区二区三区| 女的被弄到高潮叫床怎么办| 亚洲精品国产色婷婷电影| 永久网站在线| 国产久久久一区二区三区| 久久综合国产亚洲精品| 男的添女的下面高潮视频| 丝瓜视频免费看黄片| 青春草亚洲视频在线观看| 成人毛片60女人毛片免费| 精品亚洲乱码少妇综合久久| 校园人妻丝袜中文字幕| 日本-黄色视频高清免费观看| 韩国av在线不卡| 精品久久久久久电影网| 麻豆成人av视频|