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

    天津地區(qū)霧天不同高度湍流輸送特征的實(shí)驗(yàn)研究

    2016-10-25 08:09:23葉鑫欣魏偉李航張宏升
    關(guān)鍵詞:平流霧逆溫層比濕

    葉鑫欣 魏偉 李航 張宏升

    ?

    天津地區(qū)霧天不同高度湍流輸送特征的實(shí)驗(yàn)研究

    葉鑫欣 魏偉 李航 張宏升?

    北京大學(xué)物理學(xué)院大氣與海洋科學(xué)系, 氣候與海?氣實(shí)驗(yàn)室, 北京100871; ? 通信作者, E-mail: hsdq@pku.edu.cn

    利用天津 255 米氣象塔層大氣邊界層觀測資料, 分析霧日各氣象要素的特征, 研究湍流輸送規(guī)律。結(jié)果表明: 霧前, 大氣濕度較大, 逆溫層高度約為 100m, 輻射霧發(fā)生前的逆溫強(qiáng)于平流霧; 霧中, 逆溫層持續(xù)變強(qiáng)、增厚; 霧后, 逆溫層出現(xiàn)抬升, 大氣呈近中性偏不穩(wěn)定的層結(jié)特征。輻射霧過程的逆溫現(xiàn)象比平流霧明顯; 輻射霧霧頂較低, 平流霧較高; 輻射霧的消散主要受溫度影響, 平流霧主要受風(fēng)速影響。霧天氣過程改變了大氣層結(jié)結(jié)構(gòu), 夜間可能呈現(xiàn)不穩(wěn)定層結(jié)狀態(tài); 霧前和霧中不同高度的湍流垂直輸送微弱, 霧過程后期的水平輸送突然增強(qiáng)。輻射霧霧前不同高度的平均動(dòng)能數(shù)值較小, 霧中呈增大趨勢, 霧消散階段逐漸增大, 湍流動(dòng)能的增大是影響輻射霧消散的重要因素。霧前不同高度平均動(dòng)能與湍流動(dòng)能比值的突然增大可能是霧發(fā)生的湍流信號(hào)之一, 比值劇增之后降至霧前水平則為霧消散的信號(hào)。

    大氣邊界層; 霧天氣; 湍流動(dòng)能; 天津

    霧是由于空氣中懸浮著大量的微小水滴或冰晶, 使水平能見度低于 1km 的天氣現(xiàn)象。霧是對人類交通活動(dòng)影響最大的天氣現(xiàn)象之一: 霧時(shí)的能見度降低會(huì)導(dǎo)致很多交通工具(如飛機(jī)等)無法正常使用, 甚至造成意外事故, 還會(huì)給人們的日常生活及農(nóng)作物生長帶來不利影響, 也會(huì)對人的健康產(chǎn)生很大影響[1-2]。揭示霧的發(fā)生、發(fā)展和消散規(guī)律, 正確認(rèn)識(shí)霧生消過程中的湍流輸送特征, 對提高霧的生消預(yù)報(bào)和人工影響霧實(shí)驗(yàn)方案的設(shè)計(jì)是非常必要的。

    由于獲取大氣邊界層觀測資料和快速響應(yīng)資料較困難以及霧天大氣邊界層結(jié)構(gòu)較復(fù)雜等原因, 霧的研究多集中在大尺度環(huán)流背景[3-5]、中小尺度地形影響、霧天大氣邊界層特點(diǎn)[6-10]、霧中微觀物理量演變特征[11]、垂直湍流輸送特征[12-13]、霧與環(huán)境的相互影響[14]以及不同地區(qū)霧過程的特點(diǎn)[15-17]等方面, 霧中邊界層特征研究和微觀物理量(如霧滴譜等)的研究則集中在輻射霧方面, 霧中湍流相似性理論驗(yàn)證及輸送規(guī)律的研究較少涉及。

    根據(jù)生成條件可以將霧分為輻射霧、平流霧、山地霧和鋒面霧等。霧發(fā)生、發(fā)展和消散過程與湍流運(yùn)動(dòng)密切相關(guān)[18-19]。湍流運(yùn)動(dòng)既可促進(jìn)霧的發(fā)展, 也可導(dǎo)致霧的消散。Zhou 等[20]認(rèn)為湍流和輻射及重力沉降之間存在微妙的平衡, 由于輻射冷卻而產(chǎn)生的液態(tài)水會(huì)隨著湍流輸送作用的增強(qiáng), 上升并蒸發(fā)消散。Nakanishi[21]研究了湍流輸送對霧的形成和消散機(jī)制。吳彬貴[22]指出, 霧過程中平均動(dòng)能較小, 但湍流擾動(dòng)活躍, 湍流動(dòng)量輸送以垂直方向?yàn)橹? 平均動(dòng)能和湍流動(dòng)能在霧前出現(xiàn)異常增強(qiáng), 可能是平流霧的啟動(dòng)信號(hào)。輻射霧過程中, 低空往往會(huì)伴隨著逆溫層的出現(xiàn), 近地面大氣層結(jié)常常呈現(xiàn)不穩(wěn)定特征[23]。Argentini 等[24]評估了小風(fēng)弱切變條件下城市平流霧過程。

    本文針對霧天氣過程與湍流相互作用研究的相對不足, 利用天津255米氣象塔層大氣邊界層探測資料, 對比分析城市地區(qū)不同霧天氣過程氣象要素廓線演變規(guī)律和大氣湍流特征, 探討湍流輸送對霧發(fā)生、發(fā)展和消散的影響機(jī)制, 試圖找到預(yù)測霧發(fā)生的湍流信號(hào)。

    1 數(shù)據(jù)獲取與處理

    天津市位于華北平原東部, 鄰近渤海。天津255米氣象鐵塔位于天津南郊(39.06°N, 117.1°E), 海拔3.3 m。該氣象鐵塔周邊500 m內(nèi)沒有高層建筑, 300 m內(nèi)建筑物高度控制在3~30 m之間。氣象鐵塔東側(cè)為商業(yè)區(qū)(建筑物高10~15 m), 南側(cè)和西側(cè)為居民區(qū)(建筑物高15~20 m), 北側(cè)約100 m處有20~30 m高的建筑, 基本上代表了城區(qū)復(fù)雜下墊面條件[25]。

    氣象鐵塔設(shè)置 15 層觀測平臺(tái), 高度分別為 5, 10, 20, 30, 40, 60, 80, 100, 120, 140, 160, 180, 200, 220 和 250 m。每個(gè)觀測平臺(tái)在南北兩個(gè)方向上各設(shè)有伸展臂, 長度約4 m, 進(jìn)行同步觀測。每層平臺(tái)均安裝有風(fēng)速、風(fēng)向和溫濕度傳感器。風(fēng)、溫、濕梯度觀測儀分別采用長春氣象儀器廠生產(chǎn)的風(fēng)速風(fēng)向儀和溫度濕度儀。風(fēng)速、風(fēng)向、溫度和濕度觀測為全天候、連續(xù)和自動(dòng)觀測, 采樣間隔為20 s。

    采用美國 CAMPBELL 公司生產(chǎn)的 CSAT3 型超聲風(fēng)溫儀, 安裝在南伸展臂頂端, 開口向東, 采樣頻率為10 Hz。2006年1月1日, 40 m高度的超聲風(fēng)溫儀正式開始自動(dòng)連續(xù)觀測。2007 年, 在 40 m 高度同方位加裝一臺(tái) CAMPBELL 公司生產(chǎn)的KH20 濕度快速響應(yīng)儀, 采樣頻率為 10 Hz。能見度觀測采用Belfort Model 600能見度監(jiān)測儀, 安裝在距地面2 m高度處, 采樣間隔為1 min。2010年, 氣象塔觀測設(shè)備全部更新后, 塔上 15 層平臺(tái)風(fēng)向風(fēng)速、溫度和濕度觀測儀器仍采用長春氣象儀器廠設(shè)備; 40 m平臺(tái)有CAMPBELL公司 CSAT3 型超聲風(fēng)溫儀和LI-7500二氧化碳水汽分析儀, 能見度觀測沿用Belfort Model 600能見度監(jiān)測儀。地面自動(dòng)站溫濕度傳感器采用芬蘭 V?is?l? 公司生產(chǎn)的HMP45D型傳感器。

    考慮到霧天氣邊界層氣象要素演變和城市下墊面變化的影響。本文針對 2010—2011 年秋冬季 4次霧天氣過程(表 1)分別分析霧發(fā)生前、發(fā)展中和消散過程的邊界層結(jié)構(gòu)特點(diǎn)。每次霧天氣過程包括霧形成前和霧消散后數(shù)小時(shí)的資料。所有資料都進(jìn)行嚴(yán)格的質(zhì)量控制, 包括湍流數(shù)據(jù)預(yù)處理、湍流宏觀統(tǒng)計(jì)量計(jì)算、微觀特征量計(jì)算等。其中, 大氣湍流數(shù)據(jù)預(yù)處理包括剔除野點(diǎn)、趨勢項(xiàng)去除、數(shù)字濾波、坐標(biāo)旋轉(zhuǎn)等。湍流數(shù)據(jù)經(jīng)過預(yù)處理后, 去除平均風(fēng)速小于0.5 m/s的資料組, 取整點(diǎn)后30 min數(shù)據(jù)滑動(dòng)進(jìn)行湍流分析。

    表1 霧過程個(gè)例

    2 背景場分析

    圖1給出2010—2011年10, 11和12月的平均溫濕風(fēng)廓線。從整體上看, 10月氣溫在0~100 m高度表現(xiàn)為下高上低, 20 m 以下高度溫度梯度較大, 20~100 m 高度之間溫度梯度較小, 部分層結(jié)呈現(xiàn)近中性。100~120 m高度和160~220 m高度常出現(xiàn)逆溫層。查詢氣溫日變化(圖略)可知: 最高氣溫大約出現(xiàn)在北京時(shí) 15:00 左右; 最低氣溫一般出現(xiàn)在05:00, 較高層的最低氣溫可延遲 0.5~1 小時(shí)出現(xiàn), 高層氣溫隨時(shí)間的變化較低層小一些。11 和 12 月與 10 月的溫度變化規(guī)律大體上一致, 逆溫層高度和最高/低氣溫出現(xiàn)的時(shí)間和量值略有差別。

    平均風(fēng)速廓線結(jié)果表明: 10月風(fēng)速呈低層小、高層大的趨勢, 并且白天(08:00—18:00)高低空風(fēng)速差較夜間小(白天差值約 3 m/s, 夜間約 7 m/s, 圖略)。低空, 風(fēng)速隨時(shí)間變化較小; 高空, 風(fēng)速隨時(shí)間的變化較大。80~100 m 高度是上述高低空風(fēng)速變化趨勢的大致拐點(diǎn)。11 和 12 月風(fēng)速與 10 月整體趨勢較一致, 主要差別在于高低空風(fēng)速差值較小, 高低空風(fēng)速變化趨勢的拐點(diǎn)大致在 100 m高度。

    平均相對濕度廓線表明: 10月平均相對濕度在50%~70%之間。低空相對濕度晝夜差值較大, 夜間22:00 至次日 08:00, 相對濕度在 65%以上, 午時(shí), 相對濕度一般在 40%~50%; 高空相對濕度晝夜差值較小, 一般在 5%左右(圖略)。11和12月相對濕度變化趨勢與10月基本上一致。

    3 霧天氣過程溫濕風(fēng)演變分析

    從圖2的能見度和相對濕度變化情況可以確定, 2010年 10月22日22:00—10月23日06:00和10 月 24 日 02:00—06:00 期間出現(xiàn)兩次霧過程(個(gè)例 001 和 002), 最低能見度分別為 0.54 和 0.36 km。查詢歷史天氣可知, 24 日下午 12:00—17:00時(shí)段內(nèi)出現(xiàn)小雨天氣。

    3.1 個(gè)例001分析

    綜合圖3中相對濕度、比濕、氣溫和風(fēng)速的時(shí)間?高度分布情況, 分析 22 日 22:00—23 日 06:00的霧過程的邊界層特征以及霧體的特點(diǎn), 可以得到以下結(jié)果。

    霧發(fā)生前數(shù)小時(shí), 整層大氣比濕保持在一個(gè)較大的數(shù)值, 約 9 g/kg; 相對濕度同時(shí)出現(xiàn)明顯增大, 達(dá)到 85%以上, 低層大氣甚至達(dá)到 90%。150 m 左右高度出現(xiàn)逆溫層, 并逐漸發(fā)展、變厚、下壓。低空風(fēng)速持續(xù)低于1 m/s, 有利于霧的形成, 并且霧發(fā)生前十幾個(gè)小時(shí)風(fēng)速都較小。風(fēng)向由霧前的北風(fēng)逐漸轉(zhuǎn)為東偏南風(fēng)(圖略)。由于天津地區(qū)東臨渤海, 東風(fēng)帶來的暖濕氣流是霧形成的水汽條件, 也是此次霧形成的可能原因之一。

    在此次霧發(fā)生發(fā)展過程中, 大氣的相對濕度幾乎處于飽和狀態(tài), 比濕則比霧前有減小的趨勢。在整個(gè)霧過程中, 比濕與相對濕度存在一定的反對應(yīng)關(guān)系, 即相對濕度明顯變大的時(shí)刻和高度, 比濕一般都呈現(xiàn)減小的趨勢, 但兩者都比平時(shí)的水平高, 這與前人關(guān)于輻射霧的結(jié)論[25]一致。霧發(fā)生前中期, 逆溫層持續(xù)變強(qiáng)、增厚, 并下壓逐漸接地; 而霧過程的后期, 逆溫層底出現(xiàn)抬升的情況, 大氣呈現(xiàn)弱不穩(wěn)定層結(jié)特征。整個(gè)霧發(fā)生過程中, 霧區(qū)低層風(fēng)速較小(約 1 m/s), 直到霧臨近消散, 才有緩慢增加的趨勢。根據(jù)相對濕度低于 80%可以判定為非霧區(qū), 不難得出結(jié)論: 此次霧過程的大部分時(shí)間段, 霧頂高度在 100~150 m 之間, 整個(gè)霧層較淺, 只在霧發(fā)展的后期, 霧頂才略有抬升。

    霧消散前, 隨著逆溫層底的抬升, 逆溫強(qiáng)度也在減弱, 直至霧消散前 1~2 小時(shí), 大氣呈現(xiàn)弱不穩(wěn)定、近中性的特征, 對流增強(qiáng), 加快了霧的消散速度。06:00 以后, 隨著太陽輻射逐漸出現(xiàn), 大氣緩慢升溫, 湍流交換變強(qiáng), 霧加速消散。并且, 從 23日 06:00 開始, 高層的極大風(fēng)速逐漸下壓, 低層的風(fēng)速也在緩慢增加, 空氣流動(dòng)加強(qiáng), 促使霧的消散??梢哉J(rèn)為, 溫度和風(fēng)速是影響本次霧消散的重要因素。

    3.2 個(gè)例002分析

    2010年10月24日02:00—06:00的霧過程(個(gè)例002)與22日22:00—23日06:00的霧過程(個(gè)例001)時(shí)間間隔比較短, 為了更好地對比不同霧日的特點(diǎn), 兩次霧過程均呈現(xiàn)在圖3中。對個(gè)例002的相對濕度、比濕、氣溫和風(fēng)速的變化進(jìn)行分析, 可以看出002號(hào)霧過程具有以下特點(diǎn)。

    霧前 3~4 個(gè)小時(shí), 大氣的相對濕度也明顯增大, 達(dá)到80%以上, 低層大氣的相對濕度霧前接近90%。霧前逆溫層開始出現(xiàn)并逐漸發(fā)展、增厚、下壓。低空風(fēng)速持續(xù)較低, 有些時(shí)段風(fēng)速低于 1 m/s, 甚至無風(fēng)。霧前風(fēng)向雖然持續(xù)變化, 但是都有偏東的風(fēng)向, 為此次霧的形成提供了水汽。

    此次霧發(fā)生發(fā)展過程中, 大氣的相對濕度幾乎處于飽和狀態(tài), 相比于霧前比濕則略微減小。在霧發(fā)生一小時(shí)后, 逆溫層底出現(xiàn)抬升的情況, 大氣呈現(xiàn)弱不穩(wěn)定層結(jié)特征, 逆溫層減弱。與個(gè)例001不同的是, 此次霧過程逆溫層較弱, 逆溫層底的下壓、抬升也不明顯, 不易通過逆溫層判斷霧頂?shù)母叨?。但? 霧前以及霧中逆溫層的出現(xiàn)和發(fā)展依然是霧過程的重要特點(diǎn)。從02:00開始, 霧區(qū)風(fēng)速持續(xù)增大, 達(dá)到6 m/s左右; 06:00以后的霧消散期間, 在200~250 m高空出現(xiàn)風(fēng)的極大值。本次霧過程期間, 風(fēng)速一直較大, 霧中一段時(shí)間有東偏南風(fēng)向, 暖濕氣流提供了水汽輸送??梢耘卸ù舜戊F為平流霧。此次霧過程的大部分時(shí)間段, 霧頂高度在 250 m以上。霧發(fā)展的初期高度較低, 之后逐漸升高, 直到大于250 m。

    霧消散前, 大氣呈現(xiàn)弱不穩(wěn)定, 近中性, 對流增強(qiáng), 利于霧的消散。06:00 以后, 隨著太陽輻射逐漸增強(qiáng), 大氣緩慢升溫, 湍流交換變強(qiáng), 加快了霧的消散速度。同時(shí)風(fēng)速在高層達(dá)到一個(gè)極大值, 空氣流動(dòng)加強(qiáng), 促使霧消散??梢哉J(rèn)為風(fēng)速是影響此次霧消散的重要因素。

    3.3 普適性分析

    為了驗(yàn)證上述結(jié)果的普適性, 我們另外選取兩次霧天氣個(gè)例進(jìn)行分析。分別是 2010 年 11 月 19日02:00—08:00的輻射霧過程(記為003)和2011 年11 月 28 日 04:00—10:00 的平流霧過程(記為004)。

    將003 與 001 兩次輻射霧過程進(jìn)行對比(圖略), 相同點(diǎn)是: 1) 霧前比濕、相對濕度較平均狀態(tài)大, 風(fēng)速較平均狀態(tài)小, 利于霧的形成; 霧前出現(xiàn)逆溫層, 并且持續(xù)發(fā)展、增厚; 2) 霧中相對濕度持續(xù)增大, 但比濕略有減小, 保持較小風(fēng)速, 利于霧的維持, 霧初期逆溫層變強(qiáng)、增厚、下壓, 霧后半段逆溫層底抬升, 逆溫層減弱; 3) 溫度變化與 001 號(hào)霧過程比較一致, 溫度是導(dǎo)致霧消散的主要原因。不同點(diǎn)是: 1) 003 號(hào)霧過程的霧前比濕小于 001 號(hào)霧過程的9 g/kg, 相對濕度也較小, 003號(hào)霧過程能見度只略小于1 km; 2) 霧中相對濕度略小于90%, 與001 號(hào)霧過程中 95%以上的高濕狀態(tài)有一定區(qū)別, 霧頂高度在 200 m 左右, 略低于 001 號(hào)霧過程; 3) 霧消散階段風(fēng)速增大現(xiàn)象不如001 號(hào)霧過程明顯。

    將004與002兩次平流霧過程進(jìn)行對比(圖略), 相同點(diǎn)是: 1) 霧前比濕、相對濕度較平均狀態(tài)大, 逆溫層和小風(fēng)速均利于霧的形成; 2) 霧中相對濕度持續(xù)增大, 比濕略有減小, 風(fēng)速有增大趨勢, 利于平流霧的發(fā)生, 霧初期逆溫層變強(qiáng)、增厚, 并下壓逐漸接地, 霧后半段逆溫層底抬升, 逆溫層減弱; 3) 004 號(hào)霧過程霧消散時(shí)的溫度和風(fēng)速變化與 002 號(hào)霧過程趨勢一致, 風(fēng)速是導(dǎo)致霧消散的主要原因。不同點(diǎn)是: 1) 004號(hào)霧過程霧前比濕和風(fēng)速略有差異, 2) 霧中風(fēng)速尤其是低層的風(fēng)速比 002 號(hào)霧過程中小, 平流不強(qiáng), 初期霧頂高度在 150 m 左右, 很快霧頂就上移到250 m以上。

    4 秋冬季霧天的湍流輸送特征研究

    4.1 個(gè)例001分析

    圖4給出三層高度(40, 120和220 m)水平熱量輸送和垂直熱量輸送??梢钥闯? 不論是霧前還是霧持續(xù)過程中, 垂直輸送都非常弱, 水平輸送則 一直存在, 尤其是霧過程的后半段, 40 和120 m高度還出現(xiàn)較強(qiáng)的水平熱量輸送, 有利于霧的消散。檢查三層高度穩(wěn)定度參數(shù)(圖略), 2010年10月22日白天, 大氣層結(jié)都以不穩(wěn)定層結(jié)為主。霧前一小時(shí)直至霧開始消散時(shí)間段, 三層穩(wěn)定度都呈現(xiàn)明顯的穩(wěn)定層結(jié)與不穩(wěn)定間層結(jié)振蕩的趨勢, 前半時(shí)間段以穩(wěn)定層結(jié)為主, 后半段時(shí)間以不穩(wěn)定間層結(jié)為主。從霧前后層結(jié)穩(wěn)定度參數(shù)的演變看, 本次霧天氣的存在改變了大氣原有的層結(jié)結(jié)構(gòu), 降低了夜間的層結(jié)穩(wěn)定度, 使夜間部分時(shí)段呈現(xiàn)不穩(wěn)定層結(jié)狀態(tài)。這一點(diǎn)與前人研究結(jié)果相似: 輻射霧過程中由穩(wěn)定層結(jié)演變?yōu)椴环€(wěn)定層結(jié)[26-27]。

    然而, 水汽輸送結(jié)果(圖5)顯示, 不論是霧前還是霧持續(xù)過程中, 垂直水汽輸送都非常弱, 以水平輸送為主, 尤其是霧過程的前段和后段, 120 m高度和220 m高度出現(xiàn)較強(qiáng)的水平水汽輸送, 說明本次霧的水汽來源以平流輸送為主。即使在霧結(jié)束階段, 水平輸送量也能達(dá)到垂直輸送的3倍, 甚至更多。值得注意的是, 23日02:00至霧開始消散的這段時(shí)間, 三層水汽水平輸送都有增大的趨勢, 但是40 m高度是正方向, 而120 m高度和220 m高度是負(fù)方向。

    圖 6(a)是三層高度平均動(dòng)能隨時(shí)間的變化。整體來看, 16:00 以前, 三層高度平均動(dòng)能都不大。16:00 以后, 平均動(dòng)能有一個(gè)快速的升高過程, 并很快達(dá)到峰值, 40 m高度在4 m2/s2左右, 120 m高度在7 m2/s2左右, 220 m高度達(dá)到8 m2/s2以上。之后又有降低的趨勢, 其中40 m高度和120 m高度平均動(dòng)能的降低一直持續(xù)到霧前的22:00, 基本上達(dá)到最小值。但是, 220 m 高度的平均動(dòng)能卻在19:00 開始再次上升, 并且上升趨勢持續(xù)到霧前, 并在霧后1小時(shí)左右達(dá)到極大值。霧的形成需要在較低的平均動(dòng)能背景下, 此次霧過程由于120 m及以上高度的平均動(dòng)能過大, 不利于霧的形成和保持。我們可以推測霧頂?shù)母叨葢?yīng)該在 100 m 左右,只在中后期短時(shí)內(nèi)霧頂高度可能有所提升。霧消散前, 平均動(dòng)能較霧中期的平均值小, 在開始消散后的 2 小時(shí)左右才達(dá)到較大的值。可以認(rèn)為, 平均動(dòng)能對此次的霧消散有影響, 但不是決定性因素。

    圖 6(b)是三層高度湍流動(dòng)能隨時(shí)間的變化??梢钥闯? 22 日白天三層湍流動(dòng)能的波動(dòng)都較大, 波動(dòng)趨勢基本上一致。霧前的幾個(gè)小時(shí), 三層湍流動(dòng)能均開始降低; 霧中期三層湍流動(dòng)能均在波動(dòng)中穩(wěn)步上升, 在霧臨近消散前達(dá)到一個(gè)較大的數(shù)值。湍流動(dòng)能減小到一定數(shù)值是霧發(fā)生的條件之一, 可以認(rèn)為, 湍流動(dòng)能的增大是影響霧消散的重要因素。

    圖7是40, 120和220 m三層高度平均動(dòng)能與湍流動(dòng)能的比值演變情況。不難發(fā)現(xiàn), 霧前數(shù)小時(shí), 比值突然增大到 30 倍左右, 并保持一定的時(shí)間, 直到霧發(fā)生。霧持續(xù)的中間時(shí)段, 比值有一個(gè)劇增, 可達(dá) 100 倍以上。之后, 比值降低到霧前的水平, 霧開始消散, 消散期間, 比值較小。在霧頂高度以下的40 m和120 m高度, 這一特征更為明顯。吳彬貴等[28]針對平流霧的研究指出: 霧前的湍流動(dòng)能和平均動(dòng)能異常, 且前者大于后者, 可能是平流霧的啟動(dòng)信號(hào)。根據(jù)對三層高度平均動(dòng)能與湍流動(dòng)能比值的分析, 可將二者比值的突然增大看做輻射霧發(fā)生的一個(gè)信號(hào), 而霧過程中比值在劇增之后降低到霧前的水平, 可以作為輻射霧即將消散的信號(hào)(通過下文的分析可以看出, 這一比值也可以作為平流霧發(fā)生和消散的信號(hào))。

    4.2 個(gè)例002分析

    參照圖5~7后半時(shí)段的結(jié)果, 我們對002個(gè)例中平流霧過程的湍流輸送特征進(jìn)行分析。

    分析圖 4 中三層熱通量垂直輸送的方向和強(qiáng)度, 可以看出: 04:00-08:00, 40 m高度向上輸送, 120 m和220 m高度向下輸送, 120 m高度熱量垂直輸送強(qiáng)度大于220 m和40 m高度; 08:00-12:00, 三層高度都向上輸送, 40 m 高度的輸送強(qiáng)度大于120 m和220 m高度的輸送強(qiáng)度。觀察溫度的梯度情況(圖 3), 可以看出 120 m 高度的溫度梯度較大, 而220 m和40 m高度較弱, 這可能是導(dǎo)致輸送強(qiáng)度不同的重要原因。08:00-12:00, 在太陽輻射的加熱作用下, 整層大氣呈現(xiàn)弱不穩(wěn)定層結(jié), 熱量向上輸送, 但強(qiáng)度較弱。

    分析圖4中水平方向三層高度熱量輸送的情況可知, 23 日 08:00-14:00, 40 m 高度以負(fù)方向輸送為主, 平均值在0.1 K · m/s左右, 120 m和220 m高度則以正方向輸送為主, 輸送強(qiáng)度比 40 m 高度大, 其中 220 m 高度的峰值達(dá)到 0.4 K · m/s, 平均值在0.2 K · m/s 左右。14:00 以后一直到霧發(fā)生的前兩個(gè)小時(shí), 從整體上看, 水平輸送很微弱, 大部分時(shí)段幾乎為零。04:00 至霧消散階段, 三層的水平熱量輸送相對之前都增強(qiáng), 其中 40 m 高度基本上保持負(fù)方向輸送, 數(shù)值持續(xù)增大; 120 m 和 220 m 高度則表現(xiàn)出正負(fù)方向交替, 峰值和平均值都較大。

    檢查三層高度穩(wěn)定度參數(shù)(圖略), 2010年10月23 日 20:00-24 日 02:00 (霧前), 三層高度都呈現(xiàn)明顯的穩(wěn)定層結(jié)與不穩(wěn)定間層結(jié)振蕩的趨勢, 前半時(shí)段以穩(wěn)定層結(jié)為主, 后半時(shí)段以不穩(wěn)定層結(jié)為主。霧持續(xù)期間, 40 m高度由不穩(wěn)定逐漸過渡到近中性并持續(xù)到霧消散, 而120 m和220 m高度則基本上保持穩(wěn)定層結(jié)。這與持續(xù)的較大風(fēng)速有關(guān), 也受到溫度廓線的影響。直到 08:00 以后, 在更大的風(fēng)速影響下, 三層高度都調(diào)整到近中性并保持到12:00以后開始降水。

    從三層高度水汽輸送對比(圖 5)可以看出, 霧前數(shù)小時(shí)直至霧前半時(shí)段, 垂直輸送一直都非常弱, 以水平輸送為主。霧后半時(shí)段, 各層都出現(xiàn)較大的垂直輸送和水平輸送, 但這可能與霧滴附著在儀器表面有關(guān)。整體上以水平輸送為主, 平均值約為垂直輸送的3倍, 甚至更多。

    通過比較 3 個(gè)高度層的平均動(dòng)能(圖 6(a))可以看出, 三層高度的平均動(dòng)能在霧前都有一個(gè)相當(dāng)長時(shí)間的低值區(qū), 霧的出現(xiàn)也伴隨著平均動(dòng)能的迅速增大。當(dāng)平均動(dòng)能增大到一定數(shù)值之后, 霧開始消散。由此也可以判定, 霧的形成是由比較強(qiáng)的平流引起的, 而霧的消散同樣是大風(fēng)引起的, 這非常符合平流霧的定義和特征。

    從圖6(b)后半段002個(gè)例三層高度湍流動(dòng)能來看, 霧前較長時(shí)間湍流動(dòng)能很小, 隨著霧的發(fā)生迅速增大。當(dāng)湍流動(dòng)能增大到一定數(shù)值后, 霧開始消散。這也是典型的平流霧的特征。值得注意的是, 在霧的消散階段, 可以發(fā)現(xiàn)湍流動(dòng)能呈現(xiàn)上低下高的狀態(tài), 即 40 m 高度湍流動(dòng)能 > 120 m 高度湍流動(dòng)能 > 220 m高度湍流動(dòng)能。

    從圖7后半段三層高度平均動(dòng)能與湍流動(dòng)能的比值演變情況可以看到, 23 日傍晚一直到霧前, 比值都有變大的趨勢, 三層高度的比值平均增大 3 倍左右。霧中期, 40 m高度和220 m高度的比值明顯有一個(gè)突增, 分別可達(dá)40和150, 之后又迅速降低到霧前的水平, 很快霧開始消散, 消散期間的比值也基本上保持在霧前較大的水平。120 m高度變化不明顯, 但也可以看到這個(gè)趨勢。比值的突然增大可以看做霧發(fā)生的一個(gè)信號(hào); 比值在劇增之后降低到霧前的水平, 可以作為霧即將消散的信號(hào)。

    4.3 普適性分析

    同樣將003號(hào)與001號(hào)霧過程、004號(hào)與002號(hào)霧過程進(jìn)行對比分析(圖略), 得到以下一致性的結(jié)論。1) 綜合三層熱量輸送情況來看, 不論是霧前還是霧持續(xù)過程中, 垂直輸送都非常弱, 水平輸送則一直存在, 可見霧中的熱量輸送以水平輸送為主。2) 不論是霧前還是霧持續(xù)過程中, 垂直水汽輸送都非常弱, 以水平水汽輸送為主。3) 霧過程中, 三層高度的平均動(dòng)能都有增大趨勢; 霧消散階段, 三層的平均動(dòng)能先略有減小然后逐漸增大??梢钥闯? 霧的形成需要在較低的平均動(dòng)能背景下。霧消散前, 平均動(dòng)能數(shù)值增大到較大的數(shù)值, 對此次的霧消散有很大的影響。4) 霧前幾小時(shí), 三層高度的湍流動(dòng)能一致降低; 霧持續(xù)過程中, 雖然各層湍流動(dòng)能先在短時(shí)間內(nèi)降低, 但之后穩(wěn)步上升, 在達(dá)到一定的數(shù)值后, 霧開始消散。湍流動(dòng)能減小到一定的數(shù)值是霧發(fā)生的條件之一, 湍流動(dòng)能的增大是影響霧消散的重要因素。5) 霧頂以下高度的平均動(dòng)能與湍流動(dòng)能比值在霧消散前均增大。

    5 結(jié)論

    本文采用天津市255米氣象鐵塔不同高度的常規(guī)氣象和湍流觀測資料, 對秋冬季節(jié)4次霧過程的大氣邊界層結(jié)構(gòu)和湍流輸送特征進(jìn)行分析。得到以下結(jié)論。

    1) 霧前數(shù)小時(shí), 大氣邊界層整層的比濕和相對濕度數(shù)值都較大; 逆溫層發(fā)生并逐漸發(fā)展, 其中輻射霧霧前的逆溫強(qiáng)度大于平流霧霧前的逆溫強(qiáng)度; 霧前風(fēng)速較小, 平流霧發(fā)生前1~2小時(shí)呈增大趨勢。霧中逆溫層厚度和強(qiáng)度持續(xù)發(fā)展, 大氣呈弱不穩(wěn)定層結(jié)特征, 輻射霧的逆溫強(qiáng)度高于平流霧。輻射霧霧區(qū)的低風(fēng)速現(xiàn)象持續(xù)到霧消散時(shí)才略有增加; 平流霧后期則往往出現(xiàn)風(fēng)速的增大。輻射霧霧頂高度(150~200 m)略低于平流霧(大于250 m)。輻射霧的消散主要受溫度的影響, 而平流霧消散的主要影響因素是風(fēng)速。

    2) 霧前和霧過程中, 熱量和水汽的垂直輸送較弱, 水平輸送相對較強(qiáng)。逆溫層的發(fā)展導(dǎo)致霧前期的大氣層結(jié)狀態(tài)以穩(wěn)定層結(jié)為主。霧天氣的存在降低了夜間的大氣層結(jié)穩(wěn)定度, 部分時(shí)段呈現(xiàn)弱不穩(wěn)定層結(jié)狀態(tài)。霧過程中, 不同高度的平均動(dòng)能有增大趨勢; 霧消散前, 平均動(dòng)能數(shù)值有明顯增加。湍流動(dòng)能的增加是霧消散的重要因素之一, 平均動(dòng)能與湍流動(dòng)能比值的突然增大可以視為霧發(fā)生的湍流信號(hào), 而霧過程中出現(xiàn)的比值數(shù)值的迅速減小可以作為霧即將消散的湍流信號(hào)。

    參考文獻(xiàn)

    [1]吳兌, 鄧雪嬌, 范紹佳, 等. 南嶺大瑤山霧區(qū)鋒面降水的雨水化學(xué)成分研究. 中山大學(xué)學(xué)報(bào): 自然科學(xué)版, 2005, 44(6): 105?109

    [2]Elsom D. Smog alert: managing urban air quality. London: Earthscan, 1996

    [3]王麗榮, 連志鶯. 河北省中南部一次大霧天氣過程分析. 氣象, 2005, 31(4): 65?68

    [4]劉秀梅, 葉鳳娟, 鄒濤, 等. 天津近海沿岸霧的分析. 天津航海, 2005(2): 50?51

    [5]羅學(xué)義, 張雪梅, 蘭博文. 哈爾濱市區(qū)霧的特征分析及預(yù)報(bào)指標(biāo)研究. 自然災(zāi)害學(xué)報(bào), 2005, 14(2): 47?49

    [6]Fuzzi S, Facchini M C, Orsi G, et al. The Po valley fog experiment 1989. Tellus B, 1992, 44(5): 448?468

    [7]Frank G, Martinsson B, Cederfelt S I, et al. Droplet formation and growth in polluted fogs. Contributions to Atmospheric Physics, 1998, 71(1): 65?85

    [8]Leipper D F. Fog on the US west coast: a review. Bulletin of the American Meteorological Society, 1994, 75(2): 229?240

    [9]Kloesel K A. A 70-year history of marine stratocumulus cloud field experiments off the coast of California. Bulletin of the American Meteorological Society, 1992, 73(10): 1581?1585

    [10]Gultepe I, Hansen B, Cober S G, et al. The fog remote sensing and modeling field project. Bulletin of the American Meteorological Society, 2009, 90(3): 341?359

    [11]黃玉生, 黃玉仁, 李子華, 等. 西雙版納冬季霧的微物理結(jié)構(gòu)及演變過程. 氣象學(xué)報(bào), 2000, 58(6): 715?725

    [12]唐浩華, 范紹佳, 吳兌, 等. 南嶺山地濃霧的微物理結(jié)構(gòu)及演變過程. 中山大學(xué)學(xué)報(bào): 自然科學(xué)版, 2002, 41(4): 92?96

    [13]吳兌, 趙博, 鄧雪嬌, 等. 南嶺山地高速公路霧區(qū)惡劣能見度研究. 高原氣象, 2007, 26(3): 649?654

    [14]陳祥輝, 馮民學(xué). 滬寧高速公路大霧監(jiān)測與預(yù)報(bào). 中國交通信息產(chǎn)業(yè), 2003(2): 100?101

    [15]李子華, 仲良喜, 俞香仁. 西南地區(qū)和長江下游霧的時(shí)空分布和物理結(jié)構(gòu). 地理學(xué)報(bào), 1992, 47(3): 242?251

    [16]李子華, 吳君. 重慶市區(qū)冬季霧滴譜特征. 南京氣象學(xué)院學(xué)報(bào), 1995, 18(1): 46?51

    [17]羅清泉, 鮮學(xué)福. 重慶市霧水的離子組分特征. 西南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2005, 27(3): 393?396

    [18]吳彬貴, 張宏升, 張長春, 等. 華北地區(qū)平流霧過程湍流輸送及演變特征. 大氣科學(xué), 2010, 34(2): 440?448

    [19]Oliver D A, Lewellen W S, Williamson G G. The interaction between turbulent and radiative transport in the development of fog and low-level stratus. Journal of the Atmospheric Sciences, 1978, 35(2): 301?316

    [20]Zhou B, Ferrier B S. Asymptotic analysis of equilibrium in radiation fog. Journal of Applied Meteorology and Climatology, 2008, 47(6): 1704? 1722

    [21]Nakanishi M. Large-eddy simulation of radiation fog. Boundary-Layer Meteorology, 2000, 94(3): 461?493

    [22]吳彬貴. 天津地區(qū)霧天湍流輸送特征的實(shí)驗(yàn)研究[D]. 北京: 北京大學(xué), 2011

    [23]Nishikawa T, Maruyama S, Sakai S. Radiative heat transfer and hydrostatic stability in nocturnal fog. Boundary-Layer Meteorology, 2004, 113(2): 273?286

    [24]Argentini S, Mastrantonio G, Lena F. Case studies of the wintertime convective boundary-layer structure in the urban area of Milan, Italy. Boundary-Layer Meteorology, 1999, 93(2): 253?267

    [25]吳彬貴, 張宏升, 張長春, 等. 天津城區(qū)大霧過程近地面層特征研究. 氣候與環(huán)境研究, 2010, 15(2): 179?190

    [26]黃建平, 朱詩武, 朱彬. 輻射霧的大氣邊界特征. 南京氣象學(xué)院學(xué)報(bào), 1998, 21(2): 258?265

    [27]李子華. 中國近40年來霧的研究. 氣象學(xué)報(bào), 2001, 59(5): 616?624

    [28]吳彬貴, 張宏升, 王兆宇, 等. 平流霧過程湍流微結(jié)構(gòu)與能量輸送的分析研究. 北京大學(xué)學(xué)報(bào): 自然科學(xué)版, 2011, 47(2): 295?301

    Study of Turbulence Transfer at Different Levels during Fog Periods in Tianjin

    YE Xinxin, WEI Wei, LI Hang, ZHANG Hongsheng?

    Laboratory for Climate and Ocean-Atmosphere Studies, Department of Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing 100871; ? Corresponding author, E-mail: hsdq@pku.edu.cn

    Based on the observational atmospheric turbulent data obtained from 255-m tower in Tianjin, the characteristics of different atmospheric variables and turbulent transfer during the fog periods were studied. The results show that before the fog, there exists high humidity, and the height of inversion reaches to 100 m.The inversion of radiation fog is stronger than that of the advection fog. During the fog, the inversion continues strengthening. With the lifting of inversion, the neutral and unstable stratification occurres which means the dissipation of fogs. Meanwhile, the development of inversion in the radiation fog is more obvious than that in advection fogs, the height of radiation fogs is lower than that of advection fogs, and the main reasons for the dissipation of radiation and advection fogs are temperature and wind speed, respectively. The stratification is changed by the fog period, resulting in an unstable nocturnal atmosphere. The vertical transfer is weak and the horizontal transfer strengthened in the later stage. The increasing in the ratio of mean kinetic energy to turbulent kinetic energy before the fog can be treated as the signal of the occurrence of fog and the decreasing marks the dissipation of fog.

    atmospheric boundary layer; fog; turbulent kinetic energy; Tianjin

    10.13209/j.0479-8023.2015.141

    P404

    2015-03-31;

    2015-04-29; 網(wǎng)絡(luò)出版日期: 2016-03-21

    環(huán)境保護(hù)部公益性行業(yè)科研專項(xiàng)(201409001, 201309009)、國家自然科學(xué)基金(41475007)和天津市自然科學(xué)基金(13JCYBJC20000)資助

    猜你喜歡
    平流霧逆溫層比濕
    “霧”入仙境
    巴彥淖爾市臨河區(qū)逆溫層特征及對空氣污染物濃度的影響
    雅安地區(qū)近50 年濕度變化特征分析
    基于探空資料的1961—2018年新疆高空大氣比濕氣候特征分析
    南方比濕特征及其與暴雨的關(guān)系
    黃河源地區(qū)近地面逆溫層特征及形成原因分析
    平流霧下的廣西金城江城區(qū)
    平流霧罩,海上蓬萊膠州灣
    四川盆地逆溫層特征對空氣污染的影響
    2017年朝陽市地面比濕特征分析
    国产色婷婷99| 成人一区二区视频在线观看| 国产精品.久久久| 欧美+亚洲+日韩+国产| 免费电影在线观看免费观看| 最后的刺客免费高清国语| 女人被狂操c到高潮| 成人欧美大片| 国产精品久久久久久久电影| 99热网站在线观看| 午夜免费激情av| 内地一区二区视频在线| 九草在线视频观看| 亚洲精品成人久久久久久| 日韩强制内射视频| 少妇熟女欧美另类| 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 美女国产视频在线观看| 亚洲精品久久国产高清桃花| 成年av动漫网址| 日本熟妇午夜| 日韩欧美精品v在线| 九九热线精品视视频播放| 精品欧美国产一区二区三| 亚洲av电影不卡..在线观看| 国产免费男女视频| 看免费成人av毛片| 国产在线精品亚洲第一网站| 免费大片18禁| 丰满乱子伦码专区| 亚洲欧美精品专区久久| 天堂中文最新版在线下载 | 一级黄色大片毛片| 国产精品乱码一区二三区的特点| 99热精品在线国产| 十八禁国产超污无遮挡网站| 精品久久久久久久久久久久久| 亚洲欧美精品综合久久99| 午夜福利成人在线免费观看| 一级毛片我不卡| 亚洲av熟女| or卡值多少钱| 中文字幕熟女人妻在线| 听说在线观看完整版免费高清| 深爱激情五月婷婷| 亚洲精品粉嫩美女一区| 男女边吃奶边做爰视频| 日韩欧美一区二区三区在线观看| 亚州av有码| 久久精品综合一区二区三区| 观看美女的网站| 99在线人妻在线中文字幕| 成人美女网站在线观看视频| 国产成人a∨麻豆精品| 亚洲精品久久国产高清桃花| 国产单亲对白刺激| 看片在线看免费视频| 熟女电影av网| 亚洲精品日韩在线中文字幕 | 少妇裸体淫交视频免费看高清| 性插视频无遮挡在线免费观看| 午夜激情欧美在线| 夜夜夜夜夜久久久久| 又爽又黄无遮挡网站| 精品99又大又爽又粗少妇毛片| 亚洲丝袜综合中文字幕| 天天一区二区日本电影三级| 日韩欧美 国产精品| 欧美日本亚洲视频在线播放| 婷婷色综合大香蕉| 九九久久精品国产亚洲av麻豆| 深夜精品福利| 日韩成人av中文字幕在线观看| 麻豆乱淫一区二区| 国产又黄又爽又无遮挡在线| 日韩国内少妇激情av| 日日啪夜夜撸| 波多野结衣高清作品| а√天堂www在线а√下载| 麻豆成人午夜福利视频| 国产成人aa在线观看| 欧美成人免费av一区二区三区| 成人一区二区视频在线观看| 亚洲欧美成人精品一区二区| 色5月婷婷丁香| 国产精品久久久久久久电影| 不卡一级毛片| 亚洲自拍偷在线| 久久精品国产亚洲av香蕉五月| 国国产精品蜜臀av免费| 91在线精品国自产拍蜜月| 女人十人毛片免费观看3o分钟| 日本欧美国产在线视频| 久久久色成人| 春色校园在线视频观看| 免费观看的影片在线观看| 极品教师在线视频| 国产精品电影一区二区三区| 亚洲欧美日韩东京热| 欧美变态另类bdsm刘玥| videossex国产| 简卡轻食公司| 级片在线观看| 看黄色毛片网站| 日韩成人伦理影院| 日韩,欧美,国产一区二区三区 | 少妇人妻精品综合一区二区 | 日韩视频在线欧美| 日韩在线高清观看一区二区三区| 中文资源天堂在线| 伦精品一区二区三区| 国产午夜精品久久久久久一区二区三区| 综合色丁香网| 亚洲美女搞黄在线观看| 亚洲一区二区三区色噜噜| 国产一区二区在线av高清观看| 一进一出抽搐动态| 少妇熟女aⅴ在线视频| 国产高清激情床上av| 最后的刺客免费高清国语| 日韩大尺度精品在线看网址| 亚洲精品日韩在线中文字幕 | 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看 | 中文资源天堂在线| 免费av毛片视频| 精品欧美国产一区二区三| 中文字幕av在线有码专区| 国产午夜精品久久久久久一区二区三区| 狠狠狠狠99中文字幕| 久久精品国产亚洲av天美| 亚洲中文字幕日韩| 麻豆国产97在线/欧美| 一本久久中文字幕| 欧美xxxx性猛交bbbb| 欧美一区二区国产精品久久精品| 成人国产麻豆网| 国产精品,欧美在线| 偷拍熟女少妇极品色| 日韩欧美精品免费久久| 免费在线观看成人毛片| 亚洲人成网站在线播| 亚洲三级黄色毛片| 特级一级黄色大片| 亚洲人成网站在线播放欧美日韩| 中文字幕av成人在线电影| 又粗又爽又猛毛片免费看| 亚洲精品亚洲一区二区| 国产精品爽爽va在线观看网站| 免费看av在线观看网站| 国产老妇女一区| 国产成人精品婷婷| 99热这里只有是精品在线观看| 少妇熟女aⅴ在线视频| 高清日韩中文字幕在线| 搡老妇女老女人老熟妇| 精品久久久久久久末码| 一边亲一边摸免费视频| 国产成人a∨麻豆精品| 少妇人妻精品综合一区二区 | 不卡视频在线观看欧美| 国产美女午夜福利| 美女高潮的动态| 国内精品久久久久精免费| 亚洲无线观看免费| 99久久久亚洲精品蜜臀av| 亚洲综合色惰| 一区二区三区四区激情视频 | 波多野结衣高清无吗| 97超视频在线观看视频| 亚洲精品久久国产高清桃花| 欧美高清性xxxxhd video| 国产伦在线观看视频一区| 亚洲精品乱码久久久久久按摩| 丰满人妻一区二区三区视频av| 性插视频无遮挡在线免费观看| 特大巨黑吊av在线直播| 国产伦一二天堂av在线观看| 国产av在哪里看| 最近视频中文字幕2019在线8| 最近2019中文字幕mv第一页| 波多野结衣巨乳人妻| 中文亚洲av片在线观看爽| 变态另类丝袜制服| 麻豆成人午夜福利视频| 国产av麻豆久久久久久久| 亚洲中文字幕一区二区三区有码在线看| 一级毛片aaaaaa免费看小| 久久久精品大字幕| 精品久久久久久成人av| 高清在线视频一区二区三区 | 亚洲天堂国产精品一区在线| 小说图片视频综合网站| 国内揄拍国产精品人妻在线| 欧美激情在线99| 亚洲成人久久爱视频| 国产在线精品亚洲第一网站| 日韩成人伦理影院| 国产av在哪里看| 国产黄色小视频在线观看| 美女 人体艺术 gogo| 联通29元200g的流量卡| 亚洲中文字幕一区二区三区有码在线看| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲内射少妇av| 中国美女看黄片| 国产精品久久久久久久久免| 国产在视频线在精品| 最近最新中文字幕大全电影3| 亚洲av二区三区四区| 国产精品无大码| 久久综合国产亚洲精品| 一级毛片我不卡| 亚洲人成网站在线播| 亚洲国产精品合色在线| 91aial.com中文字幕在线观看| 99久久中文字幕三级久久日本| 91在线精品国自产拍蜜月| 国产黄色视频一区二区在线观看 | 亚洲成a人片在线一区二区| 亚洲人成网站在线播| 一个人看视频在线观看www免费| 天堂中文最新版在线下载 | 可以在线观看毛片的网站| 看片在线看免费视频| 老熟妇乱子伦视频在线观看| 日本免费一区二区三区高清不卡| 美女大奶头视频| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 99久国产av精品国产电影| 性欧美人与动物交配| 小蜜桃在线观看免费完整版高清| 国产极品精品免费视频能看的| 26uuu在线亚洲综合色| 少妇丰满av| 久久久精品欧美日韩精品| 国产成人影院久久av| 午夜激情福利司机影院| 国产在视频线在精品| 国产伦精品一区二区三区四那| 成年女人看的毛片在线观看| 亚洲精品色激情综合| 国产精华一区二区三区| 成人美女网站在线观看视频| 亚洲人成网站在线观看播放| 精品久久久久久久久久免费视频| 国产精品人妻久久久久久| 免费av不卡在线播放| 亚洲欧美精品自产自拍| 我的女老师完整版在线观看| 日本爱情动作片www.在线观看| 久久99热6这里只有精品| 少妇的逼水好多| 神马国产精品三级电影在线观看| 欧美性猛交黑人性爽| 非洲黑人性xxxx精品又粗又长| 18+在线观看网站| 村上凉子中文字幕在线| 久久综合国产亚洲精品| 极品教师在线视频| 久久久久国产网址| 国产蜜桃级精品一区二区三区| 中国美女看黄片| 成人美女网站在线观看视频| 欧美+亚洲+日韩+国产| 欧美性感艳星| 国产精品伦人一区二区| 国产精品乱码一区二三区的特点| 中文欧美无线码| а√天堂www在线а√下载| 国产真实乱freesex| 99在线人妻在线中文字幕| 天天躁日日操中文字幕| 国产精品蜜桃在线观看 | 大型黄色视频在线免费观看| 亚洲av一区综合| 午夜久久久久精精品| 午夜福利在线在线| 成人午夜高清在线视频| 99久久人妻综合| 成人特级黄色片久久久久久久| 天堂√8在线中文| 日本撒尿小便嘘嘘汇集6| 在线观看免费视频日本深夜| 日韩欧美一区二区三区在线观看| 综合色丁香网| 女的被弄到高潮叫床怎么办| 边亲边吃奶的免费视频| 国产精品爽爽va在线观看网站| 特大巨黑吊av在线直播| 夫妻性生交免费视频一级片| 男女下面进入的视频免费午夜| 97热精品久久久久久| 国产精品久久久久久久久免| 特级一级黄色大片| 婷婷色综合大香蕉| 久久久色成人| 国产午夜福利久久久久久| 91久久精品国产一区二区成人| 尾随美女入室| АⅤ资源中文在线天堂| 久久久久久久亚洲中文字幕| 丰满的人妻完整版| 特大巨黑吊av在线直播| 日韩中字成人| 特大巨黑吊av在线直播| 亚洲精品成人久久久久久| 欧美色欧美亚洲另类二区| ponron亚洲| 最近手机中文字幕大全| 网址你懂的国产日韩在线| 我要搜黄色片| 亚洲,欧美,日韩| 国产精品日韩av在线免费观看| 精品久久久久久久久久免费视频| 亚洲自偷自拍三级| 十八禁国产超污无遮挡网站| 丰满人妻一区二区三区视频av| 高清毛片免费看| av免费观看日本| 久久99热6这里只有精品| 久久精品人妻少妇| 国产精品美女特级片免费视频播放器| 99久久久亚洲精品蜜臀av| 99久久精品国产国产毛片| 人人妻人人澡人人爽人人夜夜 | 村上凉子中文字幕在线| 干丝袜人妻中文字幕| 又粗又爽又猛毛片免费看| 麻豆国产av国片精品| 亚洲人成网站在线播放欧美日韩| 一本久久精品| 欧美区成人在线视频| 亚洲五月天丁香| 国内精品一区二区在线观看| 欧美激情在线99| 天天一区二区日本电影三级| 欧美zozozo另类| 伊人久久精品亚洲午夜| 日本爱情动作片www.在线观看| 国产午夜精品论理片| 成年女人永久免费观看视频| 国产激情偷乱视频一区二区| 国产精品爽爽va在线观看网站| 天堂av国产一区二区熟女人妻| 长腿黑丝高跟| 成人综合一区亚洲| 九色成人免费人妻av| 国产综合懂色| 亚洲国产欧美在线一区| 最近最新中文字幕大全电影3| 免费av毛片视频| 亚洲国产精品国产精品| 国产成人91sexporn| 国产久久久一区二区三区| 婷婷色综合大香蕉| 欧美高清性xxxxhd video| 18禁黄网站禁片免费观看直播| a级毛色黄片| 亚洲精品久久久久久婷婷小说 | 成熟少妇高潮喷水视频| 亚洲欧美成人综合另类久久久 | 欧美最新免费一区二区三区| 亚洲va在线va天堂va国产| 波野结衣二区三区在线| 亚洲五月天丁香| 国产精品1区2区在线观看.| 夜夜爽天天搞| av天堂中文字幕网| 毛片女人毛片| 99久久久亚洲精品蜜臀av| 国产淫片久久久久久久久| 全区人妻精品视频| 国产片特级美女逼逼视频| 伊人久久精品亚洲午夜| 哪里可以看免费的av片| 少妇人妻精品综合一区二区 | 一级黄片播放器| 一个人看的www免费观看视频| 国产精品一区二区三区四区久久| 中文字幕制服av| 日本与韩国留学比较| 国产片特级美女逼逼视频| 人人妻人人澡人人爽人人夜夜 | 久久久久久久久久成人| 99在线人妻在线中文字幕| 天堂影院成人在线观看| avwww免费| 国产麻豆成人av免费视频| 小说图片视频综合网站| 岛国在线免费视频观看| 中文字幕久久专区| 国产精品久久久久久精品电影| 美女高潮的动态| 在线a可以看的网站| 狂野欧美激情性xxxx在线观看| 精品久久久久久久久久免费视频| 黄片wwwwww| 亚洲成a人片在线一区二区| 人人妻人人澡欧美一区二区| 99国产精品一区二区蜜桃av| 国产精品一区二区在线观看99 | 国产一区二区在线观看日韩| 亚洲国产精品国产精品| 亚洲国产欧洲综合997久久,| 中文字幕av在线有码专区| 乱人视频在线观看| 欧美激情久久久久久爽电影| 免费人成视频x8x8入口观看| 国产老妇女一区| 一本久久中文字幕| 欧美一区二区精品小视频在线| 美女 人体艺术 gogo| 自拍偷自拍亚洲精品老妇| 麻豆av噜噜一区二区三区| av国产免费在线观看| 国产精品久久电影中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲91精品色在线| 波多野结衣巨乳人妻| 91久久精品国产一区二区三区| 久久久久免费精品人妻一区二区| 欧美在线一区亚洲| 亚洲电影在线观看av| 在现免费观看毛片| 此物有八面人人有两片| 成人二区视频| 免费av毛片视频| 欧美一级a爱片免费观看看| 成年女人永久免费观看视频| 国产免费男女视频| 欧美xxxx黑人xx丫x性爽| 夜夜看夜夜爽夜夜摸| 国产真实乱freesex| 日日啪夜夜撸| 99热这里只有精品一区| 国产精品人妻久久久影院| 久久这里只有精品中国| 久久精品国产自在天天线| 美女 人体艺术 gogo| 中文字幕av成人在线电影| 久久久久网色| 寂寞人妻少妇视频99o| 插阴视频在线观看视频| a级一级毛片免费在线观看| 国产成年人精品一区二区| 高清毛片免费看| 日韩制服骚丝袜av| 久久精品国产亚洲av天美| 亚洲欧美日韩无卡精品| 国产日本99.免费观看| 亚洲精品久久久久久婷婷小说 | 男人舔奶头视频| 99热6这里只有精品| 日韩人妻高清精品专区| 可以在线观看的亚洲视频| 天堂中文最新版在线下载 | 久久久国产成人免费| 三级男女做爰猛烈吃奶摸视频| 成人特级av手机在线观看| 在线播放国产精品三级| 国产成人91sexporn| 亚洲人成网站在线播| 免费看a级黄色片| 欧美xxxx性猛交bbbb| 人妻夜夜爽99麻豆av| av在线播放精品| av在线亚洲专区| 97在线视频观看| 99精品在免费线老司机午夜| 国产成人福利小说| 色尼玛亚洲综合影院| 一个人看的www免费观看视频| 日韩av在线大香蕉| 中国美白少妇内射xxxbb| 亚洲内射少妇av| 在线播放国产精品三级| 久久欧美精品欧美久久欧美| 99在线视频只有这里精品首页| 成人亚洲欧美一区二区av| 18禁在线播放成人免费| 国产精品不卡视频一区二区| 久久精品综合一区二区三区| 亚洲在线自拍视频| 国产精品人妻久久久影院| 免费观看精品视频网站| 99在线人妻在线中文字幕| 美女被艹到高潮喷水动态| 日韩欧美三级三区| 老女人水多毛片| 亚洲av男天堂| 国产精品永久免费网站| 精品久久久久久久久亚洲| 麻豆精品久久久久久蜜桃| 日本色播在线视频| 国产亚洲av嫩草精品影院| 欧美区成人在线视频| 久久午夜福利片| .国产精品久久| 午夜老司机福利剧场| 久久鲁丝午夜福利片| 少妇猛男粗大的猛烈进出视频 | 久久欧美精品欧美久久欧美| 午夜激情福利司机影院| 此物有八面人人有两片| 少妇人妻精品综合一区二区 | 噜噜噜噜噜久久久久久91| 国产成人精品一,二区 | 国产v大片淫在线免费观看| 少妇丰满av| 亚洲第一电影网av| 哪个播放器可以免费观看大片| 日本三级黄在线观看| 免费看a级黄色片| 日韩一区二区三区影片| 国产精品久久久久久亚洲av鲁大| 桃色一区二区三区在线观看| 深夜a级毛片| 欧美3d第一页| www日本黄色视频网| 欧美成人免费av一区二区三区| 亚洲乱码一区二区免费版| 国产日本99.免费观看| 亚洲欧美精品专区久久| 欧美另类亚洲清纯唯美| 99久久中文字幕三级久久日本| 真实男女啪啪啪动态图| 我的女老师完整版在线观看| 又黄又爽又刺激的免费视频.| 亚洲av男天堂| 中文字幕人妻熟人妻熟丝袜美| 亚洲自偷自拍三级| 成人午夜精彩视频在线观看| 国产高清视频在线观看网站| 麻豆一二三区av精品| 久久精品国产清高在天天线| 欧美一级a爱片免费观看看| 综合色av麻豆| 国产伦理片在线播放av一区 | 人妻系列 视频| 在线观看一区二区三区| 久久精品国产亚洲av天美| 精品一区二区免费观看| 国产大屁股一区二区在线视频| 国产精品国产三级国产av玫瑰| 一卡2卡三卡四卡精品乱码亚洲| 久久精品国产亚洲av香蕉五月| 人妻系列 视频| 久久这里只有精品中国| 人妻少妇偷人精品九色| 全区人妻精品视频| 国产亚洲精品久久久com| 日本免费一区二区三区高清不卡| 中文字幕精品亚洲无线码一区| 97人妻精品一区二区三区麻豆| 高清在线视频一区二区三区 | 成人三级黄色视频| 国产91av在线免费观看| 我的女老师完整版在线观看| 国产亚洲av片在线观看秒播厂 | 韩国av在线不卡| 一进一出抽搐gif免费好疼| 成人性生交大片免费视频hd| 午夜福利在线在线| 99热这里只有是精品在线观看| 日韩中字成人| 别揉我奶头 嗯啊视频| 免费观看a级毛片全部| 成人鲁丝片一二三区免费| 国产午夜精品久久久久久一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| 亚洲不卡免费看| 国产精品,欧美在线| 国产亚洲av片在线观看秒播厂 | 久久草成人影院| 国产成人91sexporn| 国产精品电影一区二区三区| 在线观看美女被高潮喷水网站| 人人妻人人澡人人爽人人夜夜 | 久久亚洲国产成人精品v| 国产黄色小视频在线观看| 一级黄色大片毛片| 99久国产av精品国产电影| 久久久久久九九精品二区国产| 99热全是精品| 成人三级黄色视频| 久久精品国产亚洲av天美| АⅤ资源中文在线天堂| 中文字幕av在线有码专区| 中国国产av一级| 成人高潮视频无遮挡免费网站| 久久韩国三级中文字幕| 三级男女做爰猛烈吃奶摸视频| 久久九九热精品免费| 国产精品久久久久久亚洲av鲁大| 久久精品国产99精品国产亚洲性色| 免费观看a级毛片全部| 国产在视频线在精品| 丰满乱子伦码专区| 亚洲国产日韩欧美精品在线观看| 日本免费a在线| 欧美最新免费一区二区三区| 久99久视频精品免费| 国产一区二区三区av在线 | 国产不卡一卡二| 久久久久久久亚洲中文字幕| 狂野欧美激情性xxxx在线观看| 国产毛片a区久久久久| 国产成人影院久久av| av天堂在线播放| 九九爱精品视频在线观看| 久久久久久伊人网av|