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

    河口微結(jié)構(gòu)觀測(cè)與小尺度動(dòng)力過程研究進(jìn)展

    2016-11-09 03:52:23劉守江
    關(guān)鍵詞:微結(jié)構(gòu)河口湍流

    任 杰,劉守江

    (1.中山大學(xué)a.海洋學(xué)院;b.近岸海洋科學(xué)與技術(shù)研究中心,廣州 510275;2.西華師范大學(xué) 國(guó)土資源學(xué)院,四川 南充 637009)

    河口微結(jié)構(gòu)觀測(cè)與小尺度動(dòng)力過程研究進(jìn)展

    任 杰1,劉守江2

    (1.中山大學(xué)a.海洋學(xué)院;b.近岸海洋科學(xué)與技術(shù)研究中心,廣州 510275;2.西華師范大學(xué)國(guó)土資源學(xué)院,四川南充 637009)

    河口微結(jié)構(gòu)直接與耗散和混合現(xiàn)象聯(lián)系在一起,在能量串級(jí)中扮演重要角色。河口微結(jié)構(gòu)的研究起步于20世紀(jì)90年代,主要集中在下面幾個(gè)方面:(1)微結(jié)構(gòu)觀測(cè)技術(shù)與設(shè)備取得了突破性的進(jìn)展;(2)河口能耗、混合過程與湍流相互作用;(3)河口湍動(dòng)能平衡及浮力能量估算;(4)河口微結(jié)構(gòu)的典型尺度分析。河口微結(jié)構(gòu)觀測(cè)與小尺度過程盡管已取得系列重要進(jìn)展,但在河口界面的微結(jié)構(gòu)觀測(cè)與小尺度過程對(duì)極端性氣候-水文事件的響應(yīng)研究尚有發(fā)展的空間。

    微結(jié)構(gòu);小尺度過程;湍流;河口

    0 引 言

    在早期(20世紀(jì)50-80年代),以Pritchard[1-3],Ippen[4],F(xiàn)ischer[5],Dyer[6,7]和Officer[8]等為代表,對(duì)與鹽水運(yùn)動(dòng)相關(guān)的基本河口現(xiàn)象與過程進(jìn)行了開創(chuàng)性和奠基性的研究工作,這些研究工作主要集中在河口密度環(huán)流、鹽水入侵長(zhǎng)度、河口鹽淡水混合類型、河口混合過程、河口環(huán)流與最大渾濁帶的關(guān)系等基本和重要的水動(dòng)力學(xué)問題。鹽度的層化與混合在過去的經(jīng)典研究中雖有涉及,但其一般多局限在平均流或更宏觀的尺度上。當(dāng)近20-30年來小尺度上的細(xì)微結(jié)構(gòu)的觀測(cè)成為可能時(shí),許多河口過程的小尺度結(jié)構(gòu)才逐漸成為河口海岸學(xué)中一個(gè)非?;钴S的研究領(lǐng)域。

    海洋物理要素場(chǎng)的小尺度結(jié)構(gòu),通常是指垂向尺度小于常規(guī)海洋學(xué)觀察層次間距的海水狀態(tài)參數(shù)的層次結(jié)構(gòu)。一般而言,海水狀態(tài)參數(shù)的垂向尺度小于1m的稱為微結(jié)構(gòu)(m icrostructure)[9],河口微結(jié)構(gòu)的垂向尺度尚無明確的定義,應(yīng)該比海洋微結(jié)構(gòu)尺度要小。在微結(jié)構(gòu)現(xiàn)象中,湍流與分子過程起著主要作用。近岸河口區(qū)域的微結(jié)構(gòu)研究明顯晚于海洋微結(jié)構(gòu)研究。近10來年,國(guó)際同行在近岸河口環(huán)境的微結(jié)構(gòu)觀測(cè)和小尺度過程研究上取得了一系列重要的成果[10-17]。

    河口微結(jié)構(gòu)在能量串級(jí)中扮演重要角色,直接與耗散和混合現(xiàn)象聯(lián)系在一起,因此研究河口微結(jié)構(gòu)對(duì)于區(qū)域能量和動(dòng)量傳輸有著相當(dāng)重要的意義。

    1 河口微結(jié)構(gòu)過程

    常規(guī)調(diào)查結(jié)果顯示,河口環(huán)境水體的溫度和鹽度等狀態(tài)參數(shù)的分布和變化總是那么光滑而連續(xù)。但實(shí)際現(xiàn)象并非如此,在2008年1月23日的珠江虎門口的微結(jié)構(gòu)觀測(cè)中發(fā)現(xiàn),河口動(dòng)力環(huán)境中最重要的兩個(gè)狀態(tài)參數(shù)溫度與鹽度的垂向分布存在著許多時(shí)空尺度較小的復(fù)雜結(jié)構(gòu),它們是一系列的、由許多近乎均勻的水層和較薄的強(qiáng)梯度水層相間疊置的階梯狀結(jié)構(gòu)(圖1)。這種現(xiàn)象并不是珠江河口的特例,在世界其它河口的精密調(diào)查中也發(fā)現(xiàn)了類似的結(jié)構(gòu)[17,18]。

    圖1 珠江虎門口溫、鹽微結(jié)構(gòu)剖面(a)及垂向梯度分布(b)(2008年1月23日21∶00)

    在這種小尺度的垂向結(jié)構(gòu)中,常伴隨著顯著的流速鉛直切變[19]。我們?cè)谕握{(diào)查中發(fā)現(xiàn),在鹽度強(qiáng)梯度薄層中,流速鉛直切變模量也相應(yīng)出現(xiàn)峰值,而在各階梯層內(nèi),流速的平均切變模量明顯要小得多,兩者可以相差好幾個(gè)數(shù)量級(jí)(圖2)。計(jì)算表明,在強(qiáng)梯度薄層中,理查森數(shù)的量級(jí)為 1,呈穩(wěn)定層結(jié)狀態(tài);但在各階梯層內(nèi)理查森數(shù)較小,水體大體上處于臨界(中性)穩(wěn)定狀態(tài)。

    圖2 珠江崖門口近底層流速鉛直切變模量垂向分布((a)為瞬時(shí)流;(b)為平均流)

    在自然條件下,水體狀態(tài)參數(shù)鉛直分布的微結(jié)構(gòu),作為一些物理行為的表象,本質(zhì)上往往與多種重要的河口動(dòng)力過程或效應(yīng)聯(lián)合作用緊密聯(lián)系在一起,如剪切不穩(wěn)定(shear instabilities)、內(nèi)波剪切(internal shears)、潮汐應(yīng)變(tidal straining)、對(duì)流、密度倒置(density overturn)等,包括一些生物學(xué)過程也可能對(duì)微結(jié)構(gòu)有所貢獻(xiàn)[20]。在這個(gè)問題上,還需要進(jìn)一步的觀測(cè)、實(shí)驗(yàn)和理論研究。然而,微結(jié)構(gòu)的存在,可使聲散射增強(qiáng),引起局地聲強(qiáng)的急劇改變;對(duì)于估計(jì)通過躍層的鉛直熱交換速率,了解湍流和混合擴(kuò)散過程,都有重要的意義。

    與微結(jié)構(gòu)相關(guān)的能耗問題在珠江河口富有特色,表現(xiàn)為:(1)珠江河口能量消耗的空間分布很不均勻,存在若干高能耗區(qū)。不同部位潮平均能耗率相差可達(dá)4—5個(gè)數(shù)量級(jí)。其中,地貌尺度的動(dòng)力結(jié)構(gòu)對(duì)于整個(gè)河口消能至關(guān)重要[21-22]??臻g差異顯著的能耗現(xiàn)象與微結(jié)構(gòu)顯然有著不可分割的聯(lián)系。(2)珠江崖門河口的垂向積分尺度在觀測(cè)期間平均約5.2cm(圖3),為河口典型的流速微結(jié)構(gòu)形式,拖曳系數(shù)、擴(kuò)散系數(shù)、粗糙長(zhǎng)度等湍流特征參數(shù)都表現(xiàn)出空間-時(shí)間高度相關(guān)性,這與流速微結(jié)構(gòu)的時(shí)空分布密切聯(lián)系。然而微結(jié)構(gòu)與上述過程的深入關(guān)系則需要對(duì)這些小尺度行為作進(jìn)一步的研究。

    圖3 珠江河口崖門的積分尺度隨潮變化過程

    大量的觀測(cè)表明,這種狀態(tài)參數(shù)鉛直分布的細(xì)微結(jié)構(gòu),是河口中普遍存在的重要特征。在河口地區(qū),如圖1所示的這種流動(dòng)梯度或能源來源是非常豐富的,如潮汐變化,固體邊界以及風(fēng)的作用都可以產(chǎn)生梯度,從而提供湍流的垂向混合與水體層結(jié),形成系列復(fù)雜的微結(jié)構(gòu)。

    2 研究進(jìn)展

    2.1 微結(jié)構(gòu)觀測(cè)

    微結(jié)構(gòu)的觀測(cè)經(jīng)歷了一個(gè)較長(zhǎng)的過程,Stewart和Grant[23]率先在一個(gè)潮汐河道里進(jìn)行了湍流觀測(cè)。隨后他們進(jìn)一步應(yīng)用湍流譜分析方法討論了Kolmogorov慣性子區(qū)與慣性-平流子區(qū)(inertial-convective subrange)的動(dòng)力特征[24]。垂向剪切、剖面電導(dǎo)與溫度等這種小尺度脈動(dòng)的直接測(cè)量?jī)x器首先在美國(guó)[25-27]、加拿大[28]、俄羅斯[29,30]、德國(guó)[31](Prandke et al.,1985)和澳大利亞[32]問世,自此以后,不同的研究團(tuán)隊(duì)開始提高和發(fā)展新型微結(jié)構(gòu)儀器,如AMP和CHAMELEON[33],BAKLAN和GRIF[34],F(xiàn)LY[35]和EPSONDE[36]。近來隨著一些商業(yè)剖面儀的面市,如MSS[37],TURBOMAP[38]和PME[39],越來越多的研究者開始關(guān)注湍流微結(jié)構(gòu)動(dòng)力及其對(duì)水生生態(tài)系統(tǒng)的影響,尤其是近20—30年來,海洋微結(jié)構(gòu)研究成為海洋物理學(xué)研究中的一個(gè)非常活躍的領(lǐng)域。

    微結(jié)構(gòu)在近岸河口區(qū)域的存在,同海洋區(qū)域一樣普遍和重要,但其研究起步已是明顯晚于海洋微結(jié)構(gòu)研究數(shù)十年。過去河口海岸學(xué)者較多地關(guān)注河口過程、海岸行為的大尺度現(xiàn)象[40,41],直至近10來年,國(guó)際同行開始把目光更多地投注在近岸河口環(huán)境的微結(jié)構(gòu)觀測(cè)和小尺度過程上[10,12-17],并取得了一系列重要的進(jìn)展。

    盡管微結(jié)構(gòu)測(cè)量的歷史可以追溯到30多年前,微結(jié)構(gòu)變量的估算在伴隨著微結(jié)構(gòu)觀測(cè)的同時(shí)也就開始了,但像一些重要的微結(jié)構(gòu)特征量(如湍流動(dòng)能耗散,溫、鹽脈動(dòng),渦動(dòng)擴(kuò)散)在不同自然環(huán)境水體里的精確估算,特殊界面(躍層、鋒面)的微結(jié)構(gòu)特征,及微結(jié)構(gòu)資料處理技術(shù)(比如受波浪污染的微結(jié)構(gòu)數(shù)據(jù)的湍流分離)等問題還有待進(jìn)一步提高[42]。

    2.2 河口小尺度動(dòng)力過程研究

    河口混合過程的研究始于1960s[43,44],由于觀測(cè)技術(shù)的限制,過去的大部分研究都聚焦在平均流尺度的觀測(cè)和結(jié)構(gòu)上,及其對(duì)相關(guān)大尺度過程的解釋。湍流與動(dòng)量輸運(yùn)、鹽度混合之間的相互關(guān)系,及混合過程的動(dòng)力學(xué)機(jī)制是深入理解河口過程的關(guān)鍵問題[45]。隨著近年微結(jié)構(gòu)觀測(cè)技術(shù)的發(fā)展(如剪切剖面儀的應(yīng)用),湍流耗散可以直接測(cè)量,河口層化混合過程與湍流相互作用的研究開始了嘗試性的研究工作[46]。

    在20世紀(jì)的50-70年代,新的海洋觀測(cè)技術(shù)的發(fā)展使海洋要素的厘米尺度觀測(cè)成為可能,由此引發(fā)了對(duì)粘性和擴(kuò)散耗散尺度的直接觀測(cè)及湍流混合過程的深入分析。正如Gargett[47]指出,小尺度混合過程的研究有三種途徑:(1)根據(jù)大尺度行為的潛在原因進(jìn)行推論;(2)小尺度混合過程的直接觀測(cè);(3)以水槽實(shí)驗(yàn)和理論模型為基礎(chǔ)的過程研究。在20世紀(jì)60年代前,因湍流不能直接觀測(cè),推論法扮演了主要角色,以至于在較長(zhǎng)的一段時(shí)間內(nèi)我們對(duì)實(shí)際的微結(jié)構(gòu)認(rèn)識(shí)很不夠[48]。早期的河口微結(jié)構(gòu)測(cè)量聚焦在能量尺度(energetic scales),像湍流脈動(dòng)引起的粘性和擴(kuò)散等小尺度過程最先在深海展開,后來演變成了現(xiàn)在眾所周知的微結(jié)構(gòu)觀測(cè)。

    湍流垂向混合是河口動(dòng)量與鹽分平衡最重要的動(dòng)力過程之一,20世紀(jì)50-60年代,河口環(huán)流的經(jīng)典理論[2,3,49-53]認(rèn)為,鹽躍層(halocline)的湍流混合是河口動(dòng)力平衡項(xiàng)中的積分部分(integral part),量級(jí)最小。然而Peters[54]在Hudson河口進(jìn)行的一系列微結(jié)構(gòu)調(diào)查發(fā)現(xiàn),潮汐河口的湍流混合呈現(xiàn)出一種變化模式:在弱層化底層潮混合強(qiáng)烈,小潮期鹽躍層界面湍流混合弱,大潮落潮期則整個(gè)水體混合強(qiáng)烈,對(duì)于后者,垂向湍動(dòng)鹽度通量的典型值為0.001kg·m-2·s-1,這個(gè)數(shù)值足以消除長(zhǎng)達(dá)4h的高鹽水(haline)分層,由此說明了河口鹽躍層界面位置(cross halocline)湍流混合的重要性[55],因此河口動(dòng)力平衡分析中,小尺度的垂向混合過程在鹽度躍層的影響是不容忽視的。

    在物理海洋學(xué)中,動(dòng)量、物質(zhì)及熱的垂向湍流通量估算方法一般根據(jù)簡(jiǎn)化的穩(wěn)態(tài)湍動(dòng)能(TKE)方程(假定空間各向同性)。湍量動(dòng)能 q2的非穩(wěn)態(tài)平衡方程和標(biāo)量(如溫度、鹽度)脈動(dòng)(scalar fluctuation)方差Θ的平衡方程如下[56]:

    上式中因擴(kuò)散與平流項(xiàng)(diffusion and advection terms)通常比主要項(xiàng)小得多而被忽略。,穩(wěn)態(tài)條件下湍流動(dòng)能的剪切生成(shear production)G主要與粘性耗散(vicous dissipation)ε和浮力通量(buoyancy flux)(穩(wěn)態(tài)分層流)B平衡。ε以熱的形式損失能量,為匯項(xiàng)。在穩(wěn)定分層流下重力做負(fù)功,能量損失,B是匯項(xiàng),而在不穩(wěn)定流條件下重力做正功,B為源項(xiàng)[17]。而標(biāo)量脈動(dòng)(溫、電導(dǎo)、鹽等)的能量生成Gθ主要與標(biāo)量耗散χθ平衡。

    湍流混合涉及一個(gè)重要參量,即浮力通量的估算,垂向浮力通量是難于直接測(cè)量的參數(shù),因此常常通過一些容易直接觀測(cè)的變量去估算浮力通量。在過去的十多年,大量的學(xué)者對(duì)河口密度與流速的微結(jié)構(gòu)進(jìn)行了直接觀測(cè)[12,57-60],并間接估算浮力通量B,主要方法有:(1)Osborn[61]方法(1980),假定混合效率γ的上限為0.2,則浮力通量界限為 B≤0.2ε;(2)Ivey and Imberger[62]方法(1991),定義通量里查森數(shù) Rf,Rf=1/(1 +ε/B),Rf與湍流Froude數(shù)Frt=(LO/LE)2/3和湍流Reynolds數(shù)Fet=(LE/LK)4/3有關(guān),其中,LO=(ε/N3)1/2是 Ozmidov尺度,是穩(wěn)定分層流的最大湍渦(turbulent eddies),LK=(ν3/ε)1/4是 Kolmogorov尺度,是可能的最小尺度,LE=(ρ2)1/2?ρ/?z是Ellison尺度,是湍渦的垂向位移(vertical displacement),湍流動(dòng)能在此轉(zhuǎn)換為勢(shì)能。(3)Simpson et al.[63]方法(2005),Simpson認(rèn)為水平密度梯度的應(yīng)變(straining)導(dǎo)致漲潮期間動(dòng)力不穩(wěn)定,而落潮期間動(dòng)力穩(wěn)定,從而有,式中ξ為水面高度,?為深度平均流速。(4)Osborn and Cox方法(1972),對(duì)于穩(wěn)態(tài)和橫向均勻流,根據(jù)標(biāo)量脈動(dòng)守恒(conservation of scalar fluctuations)計(jì)算浮力通量[64]:,上述4種方法分別從平均流尺度與湍流尺度對(duì)浮力通量進(jìn)行估算。

    河口混合是通過不穩(wěn)定的小尺度運(yùn)動(dòng)來實(shí)現(xiàn)的,這種重力不穩(wěn)定機(jī)制稱之為倒置(overturn),倒置失穩(wěn)(collapse of overturn)會(huì)耗散機(jī)械能(mechanical energy),并使一些物理要素場(chǎng)(溫、鹽)發(fā)生擴(kuò)散輸運(yùn)。Thorpe[65]提出了一種發(fā)生密度倒置的觀測(cè)資料分析技術(shù),他定義了一個(gè)Thorpe尺度,即將一個(gè)包含密度倒置的瞬時(shí)密度剖面進(jìn)行重新排序,重構(gòu)為一個(gè)單調(diào)的穩(wěn)定剖面,然后計(jì)算垂向位移d’,則垂向位移的均方根即為 Thorpe尺度 LTh=rms(d’),而反映密度倒置的垂向尺寸用最大LThmax表示。Thorpe尺度也可直接估算湍流動(dòng)能的耗散,也可以用來量化垂向水層的混合率[66]。從密度倒置資料估算能耗率時(shí),假定Thorpe尺度與Ozmidov尺度同一量級(jí),LO=(e/N3)1/2,則浮力與慣性力相等,層化水流中存在最大垂向尺度運(yùn)動(dòng)。Dillon[67]認(rèn)為Thorpe尺度與Ozmidov尺度強(qiáng)烈相關(guān),并成比例,有LO=0.8LT,則能耗率可由Thorpe尺度估算:ε=(0.8LT)2N3,而渦動(dòng)擴(kuò)散系數(shù)可由Ozm idov尺度表示:Kν=0.2LO2N[28],也可表示為Kν=0.1LT

    2N[18]。在河口鹽度層化界面會(huì)表現(xiàn)出較大的Thorpe尺度分布,表明索普尺度與垂向混合存在密切關(guān)系。

    3 討 論

    河口水體層化控制垂向混合(影響熱、鹽、動(dòng)量、營(yíng)養(yǎng)物質(zhì)垂向通量)的強(qiáng)度,是影響河口動(dòng)量、能量交換與物質(zhì)輸運(yùn)的重要因子,因此河口水體層化、混合過程是淺海和河口物理學(xué)研究的關(guān)鍵問題之一[68]。盡管其間涉及的湍流是一個(gè)尚未完全解決的科學(xué)問題,甚至有關(guān)它的一些最基本的假設(shè)也受到理論和實(shí)驗(yàn)的質(zhì)疑[69]。但不進(jìn)入湍流微結(jié)構(gòu)的觀測(cè)與小尺度過程分析,是難以真正深入認(rèn)識(shí)河口動(dòng)力過程的。

    盡管在河口層化、混合過程的理論研究與微結(jié)構(gòu)觀測(cè)方面,國(guó)際上已初步取得了可喜的成果,但同時(shí)發(fā)現(xiàn),此領(lǐng)域尚有充分發(fā)展的空間。(1)河口存在豐富的界面過程,如鹽躍層、風(fēng)驅(qū)水表混合層、底剪切Ekman層、鋒面等,這些重要的河口物理過程的研究工作過去一直在大尺度過程上開展的,其湍流微結(jié)構(gòu)特征與小尺度動(dòng)力機(jī)制層面的理論成果也是令人期待的;(2)其次,河口小尺度過程對(duì)極端性氣候—水文事件的響應(yīng)研究較少(這與災(zāi)變事件的現(xiàn)場(chǎng)觀測(cè)存在極大難度有一定關(guān)系)。以往的研究多集中在利用湍流動(dòng)能平衡方程分析正常氣象條件下的微結(jié)構(gòu)特征,而缺少風(fēng)暴潮、大洪水等事件條件下的微結(jié)構(gòu)形態(tài)探討。實(shí)際上極端事件也是河口重要的動(dòng)力過程之一,它們所起的作用可能產(chǎn)生超過自身尺度數(shù)倍的響應(yīng)效應(yīng),如數(shù)天尺度的洪水作用,其漫灘淤積幅度可能比正常淤積強(qiáng)度大數(shù)個(gè)量級(jí)。因此這些過程下發(fā)生的微結(jié)構(gòu)更值得我們關(guān)注;(3)當(dāng)前涉及的河口微結(jié)構(gòu)觀測(cè)與理論分析很少與尺度問題聯(lián)系起來,它們與傳統(tǒng)的研究手段近似,僅僅局限在微結(jié)構(gòu)資料對(duì)應(yīng)的小尺度上,而與大尺度的聯(lián)系及兩者之間的解藕或藕合問題,恰恰是復(fù)雜的非線性的河口系統(tǒng)面臨的主要問題之一。

    [1]PRITCHARD DW.Estuarine hydrography[J].Advances in Geophysics,1952,1:243-280.

    [2]PRITCHARD D W.A study of the salt balance in a coastal p lain estuary[J].J.Marine Research,1954,13:133-144.

    [3]PRITCHARD DW.What is an estuary:physical viewpoint[C]//LAUFG H.Estuaries.American Association for the Advancement of Science,Washington D C,publ.no.83.1967.

    [4]IPPEN A T.Estuary and coastline hydrodynamics[M].New York:McGraw-Hill Book Co.1966.

    [5]H B FISCHER,E J LIST,R C Y KOH,et al.Mixing in inland and coastal waters[M].New York:Academic Press,1979.

    [6]DYER K R.Estuarine hydrography and sedimentation[M].Cambridge Uni.Press.1979.

    [7]KESTNER F JT,DYER K R.Estuaries:a physical introduction[J].Geographical Journal,1974,140(2):316-318.

    [8]OFFICER C B.Physical oceanography of estuaries[M].New York:John W iley and Sons.1976.

    [9]MUNK W H,ANDERSON E R.Notes on the theory of the thermocline[J].J.Marine Research,1948,7:276-295.

    [10]STACEY M T.Observations of turbulence in a partially stratified estuary[J].J.Phys Oceanogr,1999,29:1950-1970.

    [11]STACEY M T,F(xiàn)RAM JP,CHOW F K.Role of tidally periodic density stratification in the creation of estuarine subtidal circulation[J].J.Geophys.Res-Atmo.2008,113(C8):179-185.DOI:10.1029/2007JC004581.

    [12]TROWBRIDGE JH.A simple description of the deepening and structure of a stably stratified flow driven by a surface stress[J].Journal of Geophysical Research Oceans,1992,97(C10):15529-15543.

    [13]LU Y,LUECK R,HUANG D.Turbulence characteristics in a tidal channel[J].J.Phys Oceanogr,2000,30:855-867.

    [14]PETERSH,BOKHORST R.Microstructure observations of turbulentmixing in a partially mixed estuary.Part I:Dissipation rate[J].J.Phys Oceanogr,2000,30:1232-1244.

    [15]PETERSH,BOKHORST R.Microstructure observations of turbulentmixing in a partially mixed estuary.Part II:Salt flux and stress[J].J.Phys Oceanogr,2001,31:1105-1119.

    [16]DANIELA D I,KIR K.Variations of turbulent flow with river discharge in the A ltamaha River Estuary,Georgia[J].J.Geophysical Research,2007,112:1-18.

    [17]REN J,WU J X.Sediment trapping by haloclines of a river p lume in the Pearl River Estuary[J].Continental Shelf Research,2014,82:1-8.

    [18]YOSHIE K,PETER M,HAU G.The estimation of the dissipation rate from CTD profiles in the Greenl and Sea,Nelkin,Norwegian Ocean climate project technical report No.3[R].Edited by Solfrid Stre H jllo,2002:45-48

    [19]PETERSH,GREGG M C,SANFORD T B.Detail and scaling of turbulent overturns in the Pacific Equatorial Undercurrent[J]. J.Geophysical Research,1995,18:349-368.

    [20]DAVID D F,GREGORY B C,THOMASR O.Temperature and velocity m icrostructure caused by swimming fish[J].Limnology and Oceanography,1987,32(4):978-983.

    [21]劉歡.珠江河口地形致高能耗區(qū)的消能特點(diǎn)和機(jī)理研究[D].廣州:中山大學(xué).2009.

    [22]劉歡,吳超羽,徐煒銘,等.珠江河口底邊界層湍流特征量研究[J].海洋工程,2009,27:62-69.

    [23]STEWART R W,GRANT H L.Determination of the rate of dissipation of turbulent energy near the sea surface in the presence of waves[J].J.Geophysical Research,1962,67:3177-3180.

    [24]GRANT H L,HUGHES B A,VOGEL W M,et al.The spectrum of temperature fluctuations in turbulent flow[J].J.Fluid Mechanics,1968,34:423-442.

    [25]OSBORN T R.The design and performance of free-fallmicrostructure instruments at the Institute of Oceanography[R].University of British Columbia.IOUBC Manuscript Report No.30,1978.

    [26]GREGG M C,NODLAND W E,AAGAARD E E,et al.Use of fiber-optic cable with a free-fall microstructure profiler[C]// Mar.Technol.Soc.OCEANS’82 Conference Record,Washington D C:Mar.Technol,Soc.1982:260-265.

    [27]DEWEY R K,GARGET A E,OAKEY N S.A m icrostructure instrument for profiling oceanic turbulence in coastal bottom boundary layers[J].J.Atmos.Oceanic Technol,1987,4:288-297.

    [28]OAKEY N S.Determination of the rate of dissipation of turbulent energy from simultaneous temperature and velocity shearmicrostructure measurements[J].J.Phys.Oceanogr,1982,12:256-271.

    [29]MONIN A S,OZMIDOV R V.Turbulence in the Ocean[M]//Environmental Fluid Mechanics.D.Reidel Publishing Company. 1985,3.

    [30]ARVAN B M,KUSHNIKOV V V,NABATOV V N,et al.Free-falling microstructure sonde“Baklan”[M].//P P,ed.Method and Technology of Hydrophysical and Geophysical Studied of the World Ocean.Shirshov Institute of Oceanology,Moscow.1985:8-12

    [31]PRANDKE H,KRUEGER S,ROEDER W.Aufbau und Funktion einer frei fallenden Sonde zur Untersuchung der Mikrostruktur derthermohalinen Schichtung im Meer[J].Acta Hydrophysica,1985,29:65-210.

    [32]CARTER G D,IMBERGER J.Vertically risingmicrostructure profiler[J].J.Atmospheric and Oceanic Technology,1986,3:46.

    [33]MOUM JN,GREGG M,LICEN C,et al.Comparison of turbulence kinetic energy dissipation rate estimates from two ocean microstructure profiler[J].J.Atmospheric and Oceanic Technology,1995,12:346-365.

    [34]PAKA V T,NAVATOV V N,LOZOVATSKY I,et al.Ocean microstructure measurements by BAKLAN and GGIF[J].J.Atmospheric and Oceanic Technology,1999,16:1519-1532.

    [35]SIMPSON JH,CRAWFORD W R,RIPPETH T P,et al.Vertical structure of turbulent dissipation in shelf seas[J].J.Physical Oceanography,1996,26:1580-1590.

    [36]OAKEY N S.EPSONDE:an instrument tomeasure turbulence in the deep ocean[J].IEEE J.Oceanic Eng,1988,13:124-128.

    [37]PRANDKE H,STIPS A.Investigation of microstructure and turbulence in marine and limnic waters using the MST profiler[R].Ispra Joint Research Center,European Commission.Technical Note No.I.1996,96:87.

    [38]WOLK F,YAMAZAKIH,SEURONT L,et al.A new free-fall profiler formeasuring bio-physicalm icrostructure[J].J.Atmospheric and Oceanic Technology,2002,19:780-793.

    [39]STEVENSC,SM ITH M,ROSS A.SCAMP:measuring turbulence in estuaries,lakes,and coastal waters[J].NIWA-Waterand Atmosphere,1999,7(2):20-21.

    [40]DEVRIEND H J.Mathmatical modeling and large-scale coastal behaviour:Part1.Physical processes[J].J.Hydraulic Res,1991,29(6):727-740.

    [41]DEVRIEND H J,CAPOBIANCO M,CHESHER T,et al.Approaches to long-term modelling of coastalmorphology:a review[J].Coastal Eng.,1993,21(1-3):225-269.

    [42]ROGET E,LOZOVATSKY L,SANCHEZ X,et al.M icrostructure measurements in natural waters:Methodology and applications[J].Progress in Oceanography,2006,70:126-148.

    [43]BOWDON K F.Circulation and mixing in the Mersey Estuary[M].IASH Comm.Surface Waters,Publ.,1960,51:352-360.

    [44]BOWDEN K F,HOWE M R.Observation of turbulence in a tidal current[J].J.Fluid Mechanics,1963,17:271-284.

    [45]SIMPSON JH,WILLIAMSE,BRASSEUR L H,et al.The impact of tidal straining on the cycle of turbulence in a partially stratified estuary[J].Continental Shelf Research,2005,25:51-64.

    [46]MOUM JN,RIPPETH T P.Do observations adequately resolve the natural variability of oceanic turbulence[J].J.Marine Systems.2009,77(4):409-417.

    [47]GARGETT A E,OSBORN T R,NASMYTH PW.Local isotropy and the decay of turbulence in a stratified fluid[J].J.Fluid Mechanics,1984,144:231-280.

    [48]GREGG M C.The study ofmixing in the ocean:A brief history[J].Oceanography,1991,4:39-45.

    [49]PRITCHARD D W.Estuarine circulation patterns[J].Proc.Am.Soc Civ.Eng.,1955,81(717):1-11.

    [50]PRITCHARD D W.The dynamic structure of a coastal plain estuary[J].J.Marine Research,1956,15:33-42.

    [51]RATTRAY M,HANSEN D V.A similarity solution for the circulation in an estuary[J].J.Mar.Res.,1962,20:121-133

    [52]CAMERON W M,PRITCHARD D W.Estuaries[M]//HILL M N The Sea:Ideas and Observations on Progress in the Study of the Seas(ed.).New York:W iley-Interscience,1963:306-324.

    [53]HANSEN D V,RATTRAY M.Gravitational circulation in straits and estuaries[J].J.Mar.Res.,1965,23:104-122.

    [54]PETERSH,BOKHORSTR R.Observations of stratified turbulentmixing in an estuary:Neap-to-spring variations during high river flow[J].Estuarine,Coastal Shelf Sci.,1997,45:69-88.

    [55]Peters H.Spatial and temporal variability of turbulentmixing in an estuary[J].J.Mar.Res.,1999,57:805-845.

    [56]KUNDU P K.Fluid Mechanics[M].Academic Press.,1990:638.

    [57]MOUM JN.Profiler measurements of vertical velocity microstructure in the ocean[J].J.Atmos.Oceanic Technol.,1990,7:323-333.

    [58]FLEURY M,LUECK R G.Direct heat flux estimates using a towed vehicle[J].J.Physical Oceanography,1994,24:801-818.

    [59]WOLK F,LUECK R.Heat flux and mixing efficiency in the surface mixed layer[J].J.Geophys.Res.,2001,106:19547-19562.

    [60]SHAW W,TROWBRIDGE J,W ILLIAMSA.Budgets of turbulent kinetic energy and scalar variance in the continental shelf bottom boundary layer[J].J.Geophys.Res.,2001,106:9551-9564.

    [61]OSBORN T R.Estimates of the local rate of vertical diffusion from dissipation measurements[J].J.Phys.Oceanogr,1980,10:83-89.

    [62]IVEY G N,IMBERGER J.On the nature of turbulence in a stratified fluid.Part I:The energetics ofmixing[J].J.Physical O-ceanography,1991,21:650-658.

    [63]SIMPSON JH,WILLIAMSE,BRASSEUR L H,et al.The impact of tidal straining on the cycle of turbulence in a partially stratified estuary[J].Continental Shelf Research,2005,25:51-64.

    [64]TENNEKESH,Lum ley JL.A First Course in Turbulence[M].Cambridge,MA:The MIT Press,1972:300.

    [65]THORPE S A.Turbulence and mixing in a Scottish loch[J].Philos.Trans.Roy.Soc.London,1977,286 A:125-181.

    [66]W EST A,LORKE A.Small-scale hydrodynamics in lakes[J].Annual Review of Fluid Mechanics,2003,35:373-412.

    [67]DILLON T M,CALDWELL D R.The Batchelor spectrum and dissipation in the upper ocean[J].J.Geophys.Res.,1980,85(C4):1910-1916.

    [68]SIMPSON JH,BROWN J,MATTHEWS JP,et al.Tidal straining,density currents and stirring in the control of estuarine stratification[J].Estuaries,1990,13(2),125-132.

    [69]NELKIN M.In What sense is turbulence an unsolved Problem[J].Science New Series,1992,255(5044):566-570.

    A Literature Review on Estuarine M icrostructure and Small Scale Dynam ic Processes

    Ren Jie1,Liu Shoujiang2
    (1.Institute of Estuarine and Coastal Research,Marine school of SunYat-Sen University,Guangzhou 510275,China;2.Land&Resources School,China West Normal University,Nanchong Sichuan 637009,China)

    Estuarine microstructure,contacted directly with dissipation and mixing processes,plays an important role on energy cascade.Early research on estuarinem icrostructure beganin the 90’s of last century,and great progress has been made on the following aspects:(1)obvious progress of observation technology and instrument;(2)interaction of estuarine energy dissipation rate,m ixing processes and turbulence;(3)estuarine turbulent kinetic energy balance and estimate buoyant flux;(4)typical scales analysis on estuarine microstructure.Although considerable achievements have been madein estuarine m icrostructure and small scale processes,but m icrostructure measurement of estuarine interface and extreme hydrometeorological incident need more research.

    m icrostructure;small scale process;turbulent;estuary

    P737.121

    A

    10.16246/j.issn.1673-5072.2016.01.016

    1673-5072(2016)01-0107-07

    2015-12-24

    國(guó)家自然科學(xué)基金項(xiàng)目(41476072)

    任 杰(1975—),男,四川南部人,副教授,博士,主要從事近岸河口動(dòng)力沉積過程研究。

    任 杰,E-mail:renjie@mail.sysu.edu.cn

    猜你喜歡
    微結(jié)構(gòu)河口湍流
    重氣瞬時(shí)泄漏擴(kuò)散的湍流模型驗(yàn)證
    金屬微結(jié)構(gòu)電鑄裝置設(shè)計(jì)
    用于視角偏轉(zhuǎn)的光學(xué)膜表面微結(jié)構(gòu)設(shè)計(jì)
    粘結(jié)型La0.8Sr0.2MnO3/石墨復(fù)合材料的微結(jié)構(gòu)與電輸運(yùn)性質(zhì)
    他們?yōu)槭裁催x擇河口
    河口,我們的家
    特殊的河口水
    河口
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    两人在一起打扑克的视频| 黄色女人牲交| 国产精品电影一区二区三区| 别揉我奶头~嗯~啊~动态视频| 免费搜索国产男女视频| 国产精品av久久久久免费| 久久精品aⅴ一区二区三区四区| 纯流量卡能插随身wifi吗| 女同久久另类99精品国产91| 正在播放国产对白刺激| 亚洲人成电影观看| 18禁裸乳无遮挡免费网站照片 | 视频区欧美日本亚洲| 国产免费av片在线观看野外av| 亚洲精品在线观看二区| 亚洲精品成人av观看孕妇| 亚洲av片天天在线观看| 日韩有码中文字幕| 国产伦人伦偷精品视频| 精品乱码久久久久久99久播| 亚洲欧美日韩无卡精品| ponron亚洲| 亚洲精品久久午夜乱码| 欧美黄色淫秽网站| 一级,二级,三级黄色视频| 久久久久久久午夜电影 | 久久亚洲真实| 久久中文字幕一级| 亚洲成人国产一区在线观看| 亚洲欧美日韩无卡精品| 999久久久国产精品视频| 久久国产精品影院| 无限看片的www在线观看| 一级毛片女人18水好多| 操美女的视频在线观看| 高清黄色对白视频在线免费看| 精品人妻在线不人妻| 日本 av在线| 又紧又爽又黄一区二区| 50天的宝宝边吃奶边哭怎么回事| 性少妇av在线| 青草久久国产| 日韩av在线大香蕉| 欧美日韩中文字幕国产精品一区二区三区 | 精品一区二区三区视频在线观看免费 | 午夜老司机福利片| 国产午夜精品久久久久久| 久久国产精品男人的天堂亚洲| 国产99白浆流出| av电影中文网址| 脱女人内裤的视频| 欧美日韩亚洲国产一区二区在线观看| 国产色视频综合| 另类亚洲欧美激情| 午夜激情av网站| 亚洲在线自拍视频| 日韩欧美在线二视频| 在线av久久热| 国产精品秋霞免费鲁丝片| 亚洲aⅴ乱码一区二区在线播放 | 91麻豆精品激情在线观看国产 | 日韩欧美三级三区| 国产亚洲精品第一综合不卡| 黑人巨大精品欧美一区二区mp4| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产美女av久久久久小说| 国产一区在线观看成人免费| 亚洲欧美日韩无卡精品| 久久久国产欧美日韩av| 一级a爱视频在线免费观看| 老熟妇乱子伦视频在线观看| 欧美精品啪啪一区二区三区| 热re99久久国产66热| 午夜免费鲁丝| 窝窝影院91人妻| 丁香六月欧美| 在线播放国产精品三级| 成人18禁高潮啪啪吃奶动态图| 国产在线精品亚洲第一网站| 脱女人内裤的视频| 日日爽夜夜爽网站| 国产av一区在线观看免费| 亚洲精品国产色婷婷电影| 成人影院久久| 中文字幕av电影在线播放| 久久精品亚洲精品国产色婷小说| 黑人欧美特级aaaaaa片| 国内久久婷婷六月综合欲色啪| 国产亚洲欧美98| 18禁黄网站禁片午夜丰满| 欧美色视频一区免费| 在线观看免费日韩欧美大片| 免费搜索国产男女视频| 在线看a的网站| 亚洲国产精品一区二区三区在线| 51午夜福利影视在线观看| 69av精品久久久久久| 日本免费a在线| videosex国产| 99精品欧美一区二区三区四区| 亚洲精品久久成人aⅴ小说| 午夜影院日韩av| 午夜精品国产一区二区电影| 黄网站色视频无遮挡免费观看| 亚洲一区中文字幕在线| 好看av亚洲va欧美ⅴa在| 久久青草综合色| 搡老乐熟女国产| 法律面前人人平等表现在哪些方面| 两人在一起打扑克的视频| 一进一出抽搐动态| 欧美乱妇无乱码| 一级片免费观看大全| 日本a在线网址| 午夜福利免费观看在线| 亚洲av美国av| 美女高潮到喷水免费观看| 19禁男女啪啪无遮挡网站| 两个人免费观看高清视频| 美女大奶头视频| 久久国产乱子伦精品免费另类| 国产亚洲av高清不卡| 亚洲欧美激情在线| 999久久久精品免费观看国产| 不卡一级毛片| 午夜精品国产一区二区电影| 天天躁夜夜躁狠狠躁躁| 嫩草影视91久久| 国产在线观看jvid| 国产不卡一卡二| 久久天堂一区二区三区四区| 久久久国产欧美日韩av| 亚洲国产毛片av蜜桃av| 国产又色又爽无遮挡免费看| 午夜免费观看网址| 国产黄a三级三级三级人| 日韩三级视频一区二区三区| 精品人妻在线不人妻| 久久亚洲精品不卡| 日韩三级视频一区二区三区| 99国产精品免费福利视频| 亚洲中文av在线| 久久中文字幕一级| 久久久久久久久久久久大奶| 亚洲国产欧美一区二区综合| 日韩人妻精品一区2区三区| 中亚洲国语对白在线视频| 国产精品成人在线| 久久久国产成人免费| 极品人妻少妇av视频| e午夜精品久久久久久久| 国产欧美日韩综合在线一区二区| 香蕉国产在线看| 欧美日韩一级在线毛片| 美女 人体艺术 gogo| 国产免费av片在线观看野外av| 女人被躁到高潮嗷嗷叫费观| 一进一出抽搐gif免费好疼 | 露出奶头的视频| 国产有黄有色有爽视频| 男女下面插进去视频免费观看| 一边摸一边抽搐一进一小说| 精品熟女少妇八av免费久了| 久久性视频一级片| 国产精品美女特级片免费视频播放器 | 三级毛片av免费| 国产精品秋霞免费鲁丝片| 色婷婷av一区二区三区视频| aaaaa片日本免费| 欧美丝袜亚洲另类 | 国产精品国产高清国产av| 欧美日本亚洲视频在线播放| 国产av精品麻豆| 人人妻人人澡人人看| 久久久国产成人免费| 国产精华一区二区三区| 成年人黄色毛片网站| 99精品久久久久人妻精品| 国产精品美女特级片免费视频播放器 | 啦啦啦在线免费观看视频4| 日韩大尺度精品在线看网址 | 亚洲欧洲精品一区二区精品久久久| 脱女人内裤的视频| 性欧美人与动物交配| 91麻豆av在线| 级片在线观看| 午夜日韩欧美国产| 精品少妇一区二区三区视频日本电影| 国产99白浆流出| 黑人巨大精品欧美一区二区蜜桃| 久久人人精品亚洲av| 亚洲第一青青草原| 大陆偷拍与自拍| www.www免费av| 午夜福利在线观看吧| 免费看十八禁软件| 一进一出抽搐动态| 欧美成人免费av一区二区三区| 黑人猛操日本美女一级片| 麻豆一二三区av精品| 欧美日韩av久久| 亚洲精品国产区一区二| 黑人巨大精品欧美一区二区蜜桃| 国产一区在线观看成人免费| 少妇被粗大的猛进出69影院| 在线天堂中文资源库| 高清黄色对白视频在线免费看| 欧美成狂野欧美在线观看| 97超级碰碰碰精品色视频在线观看| 丰满人妻熟妇乱又伦精品不卡| 脱女人内裤的视频| 老司机在亚洲福利影院| 热re99久久精品国产66热6| 亚洲国产精品sss在线观看 | 久久久久久大精品| 日日干狠狠操夜夜爽| 久久久久久久精品吃奶| 在线观看午夜福利视频| 久久久国产成人精品二区 | 亚洲三区欧美一区| 午夜精品国产一区二区电影| 男女午夜视频在线观看| 男女做爰动态图高潮gif福利片 | 老熟妇乱子伦视频在线观看| 超碰97精品在线观看| 女人高潮潮喷娇喘18禁视频| 在线观看一区二区三区| 亚洲av日韩精品久久久久久密| 男女之事视频高清在线观看| 国产黄色免费在线视频| 国产极品粉嫩免费观看在线| 国产单亲对白刺激| 国产高清videossex| 一级毛片精品| 色婷婷久久久亚洲欧美| 一区二区三区激情视频| 新久久久久国产一级毛片| 亚洲国产看品久久| 天天添夜夜摸| 国产成人影院久久av| 免费不卡黄色视频| 中文字幕av电影在线播放| 丝袜美足系列| 一级黄色大片毛片| 精品人妻1区二区| 欧美黄色片欧美黄色片| 大型av网站在线播放| 在线观看日韩欧美| 日本三级黄在线观看| 操出白浆在线播放| 精品久久久久久成人av| 9191精品国产免费久久| 黄网站色视频无遮挡免费观看| 成人国产一区最新在线观看| 婷婷六月久久综合丁香| 久久中文字幕一级| tocl精华| 伦理电影免费视频| 亚洲va日本ⅴa欧美va伊人久久| 久久香蕉精品热| 国产无遮挡羞羞视频在线观看| 一级黄色大片毛片| 亚洲精品中文字幕在线视频| 久久久久久免费高清国产稀缺| 日本欧美视频一区| 免费高清在线观看日韩| 欧美色视频一区免费| 国产片内射在线| 欧美日韩精品网址| av视频免费观看在线观看| 日日干狠狠操夜夜爽| 人人妻人人澡人人看| 50天的宝宝边吃奶边哭怎么回事| 国产精品国产高清国产av| 精品电影一区二区在线| 中文字幕人妻丝袜一区二区| 麻豆av在线久日| 亚洲avbb在线观看| 精品国产国语对白av| 国产99白浆流出| 午夜福利一区二区在线看| 男人的好看免费观看在线视频 | 亚洲aⅴ乱码一区二区在线播放 | 国产黄色免费在线视频| 国产亚洲欧美98| 看免费av毛片| 精品少妇一区二区三区视频日本电影| 亚洲成人免费av在线播放| 91精品三级在线观看| 日韩大码丰满熟妇| 国产亚洲欧美精品永久| 国产熟女xx| 久久婷婷成人综合色麻豆| 精品人妻在线不人妻| 极品教师在线免费播放| 欧美日韩一级在线毛片| 黄色丝袜av网址大全| 好看av亚洲va欧美ⅴa在| 欧美日韩黄片免| 久久精品aⅴ一区二区三区四区| 成在线人永久免费视频| 欧美日韩福利视频一区二区| 伦理电影免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久电影中文字幕| 黄色怎么调成土黄色| 国产免费男女视频| 亚洲美女黄片视频| 51午夜福利影视在线观看| 免费在线观看亚洲国产| 人人澡人人妻人| 亚洲欧美精品综合一区二区三区| 日本欧美视频一区| 一夜夜www| 首页视频小说图片口味搜索| 精品高清国产在线一区| 超碰97精品在线观看| 亚洲精品久久成人aⅴ小说| 男女之事视频高清在线观看| 桃色一区二区三区在线观看| 9191精品国产免费久久| 国产av在哪里看| 日韩视频一区二区在线观看| 日韩有码中文字幕| 伦理电影免费视频| av免费在线观看网站| 国产成人系列免费观看| 亚洲精品中文字幕一二三四区| 日本vs欧美在线观看视频| 国产1区2区3区精品| 亚洲成人精品中文字幕电影 | 亚洲一卡2卡3卡4卡5卡精品中文| 我的亚洲天堂| 少妇粗大呻吟视频| 国产三级黄色录像| 午夜精品久久久久久毛片777| 最近最新中文字幕大全电影3 | 在线永久观看黄色视频| 一级片免费观看大全| 国产精品 欧美亚洲| 国产伦一二天堂av在线观看| 丁香欧美五月| 美女高潮到喷水免费观看| 久99久视频精品免费| av电影中文网址| 伊人久久大香线蕉亚洲五| 亚洲精品美女久久久久99蜜臀| av超薄肉色丝袜交足视频| 狂野欧美激情性xxxx| 午夜福利影视在线免费观看| 亚洲精品在线观看二区| 国产一区二区三区在线臀色熟女 | 欧美成人午夜精品| 亚洲视频免费观看视频| 免费人成视频x8x8入口观看| 日本wwww免费看| 成人精品一区二区免费| 波多野结衣高清无吗| 国产精品爽爽va在线观看网站 | 在线十欧美十亚洲十日本专区| 国产三级在线视频| 精品久久久久久久毛片微露脸| 国产精品美女特级片免费视频播放器 | 国产精品99久久99久久久不卡| 国产一区在线观看成人免费| 18禁美女被吸乳视频| 亚洲一区高清亚洲精品| 日韩免费高清中文字幕av| 婷婷精品国产亚洲av在线| 亚洲av电影在线进入| 精品一区二区三区四区五区乱码| 看黄色毛片网站| 日韩欧美国产一区二区入口| 90打野战视频偷拍视频| 国产麻豆69| 欧美乱码精品一区二区三区| 久久久久亚洲av毛片大全| 一个人免费在线观看的高清视频| 黄色丝袜av网址大全| 欧美日本亚洲视频在线播放| 日韩大码丰满熟妇| xxx96com| 久久婷婷成人综合色麻豆| 俄罗斯特黄特色一大片| 桃色一区二区三区在线观看| 欧美日韩瑟瑟在线播放| 日韩欧美一区视频在线观看| 又大又爽又粗| 怎么达到女性高潮| 久久久国产欧美日韩av| а√天堂www在线а√下载| 无人区码免费观看不卡| 好看av亚洲va欧美ⅴa在| 亚洲欧美日韩无卡精品| 亚洲国产欧美一区二区综合| 久久人人爽av亚洲精品天堂| 97碰自拍视频| 亚洲欧美一区二区三区黑人| www.精华液| 亚洲中文字幕日韩| 国产精品爽爽va在线观看网站 | 91在线观看av| avwww免费| 90打野战视频偷拍视频| 日本vs欧美在线观看视频| 久久精品亚洲熟妇少妇任你| 日本欧美视频一区| 99riav亚洲国产免费| 日日摸夜夜添夜夜添小说| 亚洲欧美日韩无卡精品| 一边摸一边抽搐一进一小说| 亚洲男人的天堂狠狠| 成人特级黄色片久久久久久久| 亚洲av成人av| 免费在线观看亚洲国产| 亚洲中文字幕日韩| 国产精品自产拍在线观看55亚洲| 精品日产1卡2卡| 亚洲国产毛片av蜜桃av| 国产亚洲欧美98| 丰满饥渴人妻一区二区三| 欧美激情极品国产一区二区三区| 日本vs欧美在线观看视频| 国产乱人伦免费视频| 久久精品影院6| 亚洲成国产人片在线观看| 水蜜桃什么品种好| 国产欧美日韩精品亚洲av| 无人区码免费观看不卡| 99国产精品一区二区三区| 婷婷六月久久综合丁香| 午夜精品久久久久久毛片777| 中文字幕精品免费在线观看视频| a级片在线免费高清观看视频| 69精品国产乱码久久久| 国产激情欧美一区二区| 女人高潮潮喷娇喘18禁视频| 午夜福利在线观看吧| 久久青草综合色| 亚洲自拍偷在线| 欧美精品一区二区免费开放| 99久久99久久久精品蜜桃| 一a级毛片在线观看| 18禁裸乳无遮挡免费网站照片 | 久久久久久久午夜电影 | 桃色一区二区三区在线观看| 午夜激情av网站| 国产一区二区在线av高清观看| 久久久久国产一级毛片高清牌| 精品国产超薄肉色丝袜足j| av在线天堂中文字幕 | 亚洲色图av天堂| 国产午夜精品久久久久久| 色综合欧美亚洲国产小说| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 国产在线精品亚洲第一网站| 美女大奶头视频| 成年人黄色毛片网站| 麻豆一二三区av精品| 在线观看一区二区三区| 欧美色视频一区免费| 老司机深夜福利视频在线观看| 亚洲成av片中文字幕在线观看| 久久久久国产一级毛片高清牌| 国产成人欧美| 交换朋友夫妻互换小说| 一区二区三区国产精品乱码| 午夜免费观看网址| 在线av久久热| bbb黄色大片| www日本在线高清视频| videosex国产| 91大片在线观看| 亚洲国产毛片av蜜桃av| 男人操女人黄网站| 欧美日韩精品网址| 淫秽高清视频在线观看| 国产aⅴ精品一区二区三区波| 日本五十路高清| 国产色视频综合| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇中文字幕五十中出 | 午夜视频精品福利| 电影成人av| 亚洲av成人一区二区三| 久久欧美精品欧美久久欧美| 91av网站免费观看| 黄色怎么调成土黄色| 在线观看www视频免费| 一个人免费在线观看的高清视频| xxx96com| 成人三级黄色视频| 极品人妻少妇av视频| 久久久久久免费高清国产稀缺| 精品免费久久久久久久清纯| 国产精品电影一区二区三区| 黄片小视频在线播放| 成熟少妇高潮喷水视频| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 神马国产精品三级电影在线观看 | 一夜夜www| 黄色毛片三级朝国网站| 18禁国产床啪视频网站| 精品欧美一区二区三区在线| 亚洲 国产 在线| 国产色视频综合| 老司机在亚洲福利影院| 国产精品av久久久久免费| 免费在线观看日本一区| 午夜福利一区二区在线看| 亚洲成国产人片在线观看| 嫩草影视91久久| 一区二区三区精品91| 一边摸一边抽搐一进一出视频| 久久久精品欧美日韩精品| 黄片播放在线免费| 男女做爰动态图高潮gif福利片 | 欧美日本中文国产一区发布| 夫妻午夜视频| 亚洲美女黄片视频| 免费在线观看黄色视频的| 一级黄色大片毛片| 天堂中文最新版在线下载| 看片在线看免费视频| 男女床上黄色一级片免费看| svipshipincom国产片| 在线视频色国产色| 色婷婷av一区二区三区视频| 国产精品综合久久久久久久免费 | 日韩欧美一区二区三区在线观看| 狂野欧美激情性xxxx| 亚洲全国av大片| 成年人黄色毛片网站| 久久国产精品影院| 亚洲av熟女| 一本大道久久a久久精品| 90打野战视频偷拍视频| 怎么达到女性高潮| 日韩欧美三级三区| 国产av一区在线观看免费| 久久亚洲精品不卡| 亚洲欧美精品综合久久99| 亚洲伊人色综图| bbb黄色大片| 国产高清视频在线播放一区| 中文欧美无线码| 亚洲中文字幕日韩| 一级,二级,三级黄色视频| 日韩精品中文字幕看吧| 中文字幕另类日韩欧美亚洲嫩草| 丁香六月欧美| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 免费av毛片视频| 黑人猛操日本美女一级片| 80岁老熟妇乱子伦牲交| 亚洲国产中文字幕在线视频| 免费少妇av软件| 中出人妻视频一区二区| 最好的美女福利视频网| 老司机午夜十八禁免费视频| 女人被躁到高潮嗷嗷叫费观| 成年人黄色毛片网站| 妹子高潮喷水视频| 亚洲人成电影观看| 欧美人与性动交α欧美软件| 香蕉丝袜av| 久久精品亚洲精品国产色婷小说| 久久久久久大精品| 丝袜美足系列| 国产精品爽爽va在线观看网站 | 少妇粗大呻吟视频| 欧美成狂野欧美在线观看| 亚洲三区欧美一区| 日韩大码丰满熟妇| 国产视频一区二区在线看| 中国美女看黄片| 女同久久另类99精品国产91| 免费看十八禁软件| 国产视频一区二区在线看| 这个男人来自地球电影免费观看| 亚洲精品成人av观看孕妇| 中国美女看黄片| 又紧又爽又黄一区二区| 丝袜人妻中文字幕| 一边摸一边做爽爽视频免费| 色精品久久人妻99蜜桃| 一级a爱视频在线免费观看| 香蕉国产在线看| 女人被躁到高潮嗷嗷叫费观| 99精品久久久久人妻精品| 纯流量卡能插随身wifi吗| 国产极品粉嫩免费观看在线| 乱人伦中国视频| 极品人妻少妇av视频| www.自偷自拍.com| 极品人妻少妇av视频| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| 国产xxxxx性猛交| 亚洲男人的天堂狠狠| 国产欧美日韩一区二区三| 精品第一国产精品| 国产欧美日韩精品亚洲av| 国产成人免费无遮挡视频| 亚洲狠狠婷婷综合久久图片| 无限看片的www在线观看| 亚洲av熟女| 亚洲欧美日韩另类电影网站| 大香蕉久久成人网| 麻豆久久精品国产亚洲av | 一二三四在线观看免费中文在| 国产精品日韩av在线免费观看 |