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

    大氣邊界層物理研究進展

    2013-09-22 05:32:38劉輝志馮健武王雷洪鐘祥
    大氣科學 2013年2期
    關鍵詞:下墊面陣風邊界層

    劉輝志 馮健武 王雷 洪鐘祥

    中國科學院大氣物理研究所大氣邊界層物理和大氣化學國家重點實驗室,北京100029

    1 引言

    大氣邊界層是指離地球表面約1~2 km高度的低層大氣。由于大氣邊界層受地球表面的影響最大,該層大氣有著區(qū)別于上層自由大氣顯著不同的特征,例如各種氣象要素(氣溫、濕度和風速等)日變化較大、垂直梯度較大等。大氣邊界層是人類的生活和生產(chǎn)活動主要場所。由人類活動帶來污染物的排放、傳輸和轉(zhuǎn)化大部分發(fā)生在該層,因此大氣邊界層的環(huán)境問題直接影響到人類的健康和生存。大氣邊界層同時也是地球各個圈層相互作用的關鍵區(qū)域。大氣邊界層的變化直接影響到地圈、水圈、冰雪圈和生物圈與大氣圈的能量和物質(zhì)交換過程,同時對天氣和氣候產(chǎn)生重要的影響。天氣及氣候模式中大氣邊界層物理過程參數(shù)化方案的改進是提高其模擬性能的關鍵科學問題之一,也是當前大氣科學研究的基礎前沿問題。由于全球變化研究包括氣候異常、生態(tài)環(huán)境惡化、水資源短缺等問題以及可持續(xù)發(fā)展研究等的需要,大氣邊界層物理研究已成為大氣科學研究的前沿學科之一。

    中國科學院大氣物理研究所大氣邊界層物理和大氣化學國家重點實驗室(簡稱LAPC)定位于低層大氣中物理和化學過程的基礎研究。自實驗室建立來,我們一直將大氣邊界層物理的基本理論作為重點研究方向。近年來,LAPC在大氣邊界層物理和大氣化學聯(lián)網(wǎng)觀測研究,大氣邊界層實驗、理論與參數(shù)化研究,大氣化學過程與氣候變化的相互作用模擬研究,陸地生態(tài)系統(tǒng)與大氣間碳氮交換研究,空氣質(zhì)量多模式集合預報系統(tǒng)研制與應用等方面取得了重要進展。對于大氣邊界層物理方面,近年來主要研究方向包括:(1)城市復雜下墊面湍流相干結(jié)構和邊界層陣風機理,非均勻下墊面大氣邊界層結(jié)構和交換過程;(2)不同生態(tài)系統(tǒng)地—氣湍流物質(zhì)、能量交換規(guī)律及特征;(3)海洋大氣邊界層物理過程,數(shù)值模式中的大氣邊界層參數(shù)化。對于 2008年之前的工作,胡非等(2003)和張美根等(2008)已作了較全面的總結(jié),因此本文將扼要介紹近4年來(2009~2012)LAPC在大氣邊界層物理方面取得的最新研究進展。

    2 大氣邊界層結(jié)構和理論研究

    大氣邊界層氣象學是以湍流理論為基礎的,在湍流理論有了一定的發(fā)展之后才得以有邊界層氣象學的產(chǎn)生。大氣邊界層的基礎理論早在 20世紀50年代末已經(jīng)基本形成。到20世紀70年代末,對均勻下墊面大氣邊界層物理結(jié)構有了較全面的認識,大氣邊界層物理學開始作為氣象學的一門相對獨立的分支學科出現(xiàn)。目前,大氣邊界層領域主要圍繞解決大氣數(shù)值模式中邊界層和地表通量參數(shù)化的問題上展開研究。LAPC近年來在深入認識湍流機制和本質(zhì)的基礎上,在湍流相干結(jié)構和邊界層陣風機理方面取得了重要的進展。

    在我國華北地區(qū),春季常發(fā)生揚沙和沙塵暴天氣,通常伴隨有冷鋒過境及強風,嚴重影響著人們的生活和生產(chǎn)活動。對起沙揚沙過程基本物理規(guī)律的深入認識是預報和減少沙塵災害、改善空氣質(zhì)量的基礎。通過分析北京325米氣象塔觀測的風速梯度資料以及渦動相關資料(Zeng et al.,2010;Cheng et al.,2011,2012a,2012b,2012c),曾慶存等(Zeng et al.,2010)提出了自然界陣風產(chǎn)生的一種新機制,即自然界陣風產(chǎn)生的機制不僅有劇烈天氣的外部動力作用原因,同時也是邊界層湍流內(nèi)部自組織結(jié)構發(fā)展的結(jié)果。在強風天氣中,當風速足夠大,使得地面摩擦速度大于某一臨界值時,在當?shù)仄鹕硥m。但是大中尺度系統(tǒng)下沉氣流使沙塵只能積聚在大氣邊界層的低層(約200 m厚)。由于陣風存在明顯的相干結(jié)構(圖1);陣風風速波峰期(即陣風時)伴隨著下沉氣流,波谷期(即陣風歇時)伴隨著上升氣流,當陣風停歇時上升氣流才能將沙塵傳入邊界層上層,然后借助大中尺度的系統(tǒng)性上升氣流再往上傳輸進入對流層?;趯﹃囷L機理的深入認識,發(fā)展了一個利用頻譜和混沌時間序列突變點檢測理論來提取陣風信號的方法(Cheng et al.,2012a),提出了一個新的一個陣風擾動參數(shù)化方案(Cheng et al.,2012b),導出地表摩擦速度新的理論公式,并以此改進了空氣質(zhì)量模式中的起沙和揚沙過程的模擬(Cheng et al.,2012c)。

    在大氣邊界層物理結(jié)構研究方面,發(fā)現(xiàn)了臺風、寒潮大風、沙塵暴過程和一般平穩(wěn)天氣大邊界層湍流結(jié)構完全具有相同的概率密度分布型(Liu et al.,2010),即具有穩(wěn)定分布和標度律,這一結(jié)果為強天氣條件下大氣邊界層湍流的理論建模和參數(shù)化提供了科學參考。另外,我們還給出了兩層大氣模式中滾軸渦漩(大氣中一種重要的有組織結(jié)構)的解析解,揭示出下層不穩(wěn)定層結(jié)對流活動與上層穩(wěn)定層結(jié)中重力波的相互作用是滾軸渦漩形成的機制之一(Liu and Sang,2009)。此外,LAPC科研人員采用飛機觀測大氣湍流數(shù)據(jù),通過分析大氣邊界層頂?shù)男〕叨韧牧鹘Y(jié)構來研究層積云覆蓋邊界層頂部的夾卷層特征(代成穎, 2012),發(fā)現(xiàn)了:夾卷層內(nèi)一直存在湍流,只是湍流強度在逆溫層底部開始向上逐漸減弱;云頂風切變有助于產(chǎn)生湍流,能夠增強云頂過渡層內(nèi)大氣的夾卷混合,進而減弱逆溫強度,縮小逆溫高度與湍流層頂之間距離;風切變產(chǎn)生的湍流能夠促使云頂多界面分離,使夾卷層、稀薄云層和逆溫層厚度明顯增加;當云頂處存在強風切變時,局地強風切變還能夠在云頂上方產(chǎn)生間歇性湍流。研究結(jié)果可為數(shù)值模式中邊界層過程的參數(shù)化提供科學的參考依據(jù)。

    3 大氣邊界層物理觀測實驗

    均勻下墊面的大氣邊界層近地層特征在Monin-Obukhov相似性理論的指導下研究比較成熟,但地球表面最普通存在的非均勻和復雜下墊面大氣邊界層尚未得到很好的解決。近年來,大氣邊界層研究重點從均勻下墊面逐漸轉(zhuǎn)向非均勻非定常下墊面。另一方面,觀測技術的迅猛發(fā)展為非均勻下墊面的直接觀測研究提供可能性。為了揭示真實下墊面條件下大氣邊界層特征和規(guī)律,LAPC在發(fā)展大氣邊界層理論并為模式發(fā)展提供基礎性支持的同時,利用鐵塔、汽艇、雷達、飛機和衛(wèi)星遙感等手段,在我國各種典型地表下墊面開展了有針對性的綜合觀測試驗,例如城市(北京)、半干旱草原(吉林通榆、錫林浩特)、山地(青藏高原)、海洋(北冰洋)等下墊面。上述實驗獲得了區(qū)別于傳統(tǒng)水平均勻下墊面邊界層的新發(fā)現(xiàn),積累了大量寶貴的觀測資料。

    圖1 2002年3月20日不同高度、不同時間間隔下1 min平均水平風速和垂直風速變化曲線。豎線表示陣風的相干結(jié)構:虛線代表陣風峰期,實線代表陣風谷期;陰影部分代表受對流和湍流影響的陣風相干性;“×”代表與相干結(jié)構相反的情況。(引自Zeng et al.,2010)Fig.1 1-min averaged horizontal downwind velocity and vertical velocities (u; w) at various heights and time intervals on 20 March 2002.Vertical lines show the coherent structure of gust-wind, dashed lines are the moments of the peaks in strong gust-wind, and solid lines are wind speed valleys.Shaded areas represent the coherent structures of gust-wind, but disturbed by convection or turbulence.Those marked “×” are out of phase from the coherent structure.(From Zeng et al., 2010)

    3.1 青藏高原觀測實驗

    在青藏高原的觀測研究方面,LAPC近年來聯(lián)合青藏高原研究所在喜馬拉雅山北側(cè)珠峰北坡絨布河谷地區(qū)聯(lián)合進行了兩次大氣強化觀測實驗(2006、2007年 5~6月)。該實驗在一個由海拔5800 m水平截面和沿河谷中部地區(qū)的垂直截面組成的封閉型河谷中進行。河谷中部架設了自動氣象站和風溫廓線儀,觀測地面風、大氣輻射、熱力狀況和垂直風廓線等參數(shù)。根據(jù)大氣物質(zhì)守恒定律,通過分析研究穿越河谷中部垂直截面的大氣物質(zhì)通量可初步定量估計河谷內(nèi)外大氣物質(zhì)的交換量。此外,結(jié)合同期的大尺度環(huán)流資料,該實驗可為研究珠峰北坡絨布河谷地區(qū)局地環(huán)流導致的地氣間物質(zhì)交換過程以及觀測期間南亞夏季風活動對地氣交換過程的影響提供第一手觀測資料。

    上述的資料分析得出了許多有意義的結(jié)果和發(fā)現(xiàn)(Zou et al.,2009;Zhou et al.,2009;Ma et al.,2009;周立波等,2010;Zhou et al.,2011),揭示了以喜馬拉雅山區(qū)為代表的青藏高原大型山地“山谷風”和“冰川風”環(huán)流的耦合機理。喜馬拉雅山區(qū)特殊的局地環(huán)流系統(tǒng)是強烈的太陽輻射在復雜下墊面上形成的熱動力復合環(huán)流,主要是“山谷風”和“冰川風”相互耦合的產(chǎn)物(Zou et al.,2009)。該環(huán)流明顯有別于傳統(tǒng)山地氣象學理論以及世界其它山地中的觀測結(jié)果,耦合環(huán)流系統(tǒng)的物質(zhì)輸送能力遠遠大于過去學術界的傳統(tǒng)認識,表明了以喜馬拉雅山區(qū)為代表的青藏高原大型山地中以局地環(huán)流為媒介的低層大氣與自由大氣間的物質(zhì)交換過程的重要性。分析結(jié)果表明:喜馬拉雅山區(qū)的近地面湍流熱交換接近高原其他地區(qū),其感熱輸送為69.9 W m?2,潛熱輸送為 32.2 W m?2,總熱通量為102.1 W m?2,波文比為2.17。由于局地環(huán)流的存在,極大提高了該地區(qū)的地氣交換能力,其物質(zhì)交換能力可以達到6.4×1011m3d?1,約為整個封閉河谷體積的 36倍。該地區(qū)地氣間的物質(zhì)輸送與南亞夏季風活動密切相關(周立波等,2010)。南亞夏季風活動弱(季風中斷期)則物質(zhì)交換量大(9.7×106m3s?1),季風活動強(季風活躍期)則物質(zhì)交換量?。?.6×106m3s?1),表明南亞夏季風可能通過改變局地大氣輻射和熱力狀況影響到該地區(qū)地氣間的物質(zhì)交換過程(圖2)。

    3.2 半干旱區(qū)地氣相互作用觀測實驗

    為了研究半干旱區(qū)不同下墊面地氣相互作用過程,LAPC和東亞中心的科研人員在吉林省通榆縣建立了地氣交換過程長期觀測實驗站(簡稱通榆實驗站)。該實驗站是國際能量和水分循環(huán)加強觀測計劃(CEOP)的地面觀測基準站之一,目前已經(jīng)積累了多年連續(xù)的半干旱區(qū)地表過程觀測資料。該觀測實驗與過去在我國進行陸面過程大氣邊界層觀測實驗HEIFE及IMGRASS最大的區(qū)別在于:改變了過去僅在加強期(一般在夏季)進行湍流通量的加強觀測,而采取一年四季連續(xù)的觀測地氣間通量交換,研究地氣通量交換過程的日、月、季及年變化,為氣候模式中半干旱地區(qū)陸面過程參數(shù)化的改進提供了可靠的驗證數(shù)據(jù)。通過分析通榆實驗站的觀測資料,我們在半干旱區(qū)近地層地氣之間能量和物質(zhì)交換以及水循環(huán)特征研究方面取得了一些有意義的結(jié)果(涂剛等,2009;Liu et al.,2012a),并給出了該半干旱區(qū)地表基本物理參數(shù)(空氣動力學粗糙度、熱力學粗糙度和整體輸送系數(shù))的變化特征(馮健武等,2012)。地氣之間的水循環(huán)過程是陸面過程中關鍵過程之一,直接影響到地表水分收支、植被的生產(chǎn)力和碳循環(huán)等重要的生態(tài)過程,從而間接地影響天氣和氣候。通過分析通榆站觀測資料,發(fā)現(xiàn)了地表實際蒸散在生長季的日總量最大值在退化草地下墊面為2.0~4.5 mm d–1,在農(nóng)田下墊面為1.5~5.5 mm d–1。降雨量的季節(jié)分布是地表實際蒸散季節(jié)變化的一個重要的影響因子。在半小時時間尺度上,大氣水分需求是地表實際蒸散的主要控制因子。生物控制因子在季節(jié)尺度上對半干旱區(qū)地表實際蒸散的影響較顯著,主要是通過表面導度季節(jié)和年際變化來實現(xiàn)。地表實際蒸散在年尺度上取決于年總降雨量(Liu et al.,2012a)。

    在中德合作項目“不同放牧強度對內(nèi)蒙古草原生態(tài)系統(tǒng)物質(zhì)流的影響”支持下,我們于 2006~2009年在中國科學院內(nèi)蒙古草原生態(tài)系統(tǒng)定位站開展內(nèi)蒙古半干旱草原下墊面地氣能量和物質(zhì)交換過程的觀測實驗。為了揭示放牧對地氣交換過程的影響,我們在該地區(qū)三種典型的羊草樣地(冬季放牧、持續(xù)放牧和過度放牧)進行了對比實驗。通過對觀測資料進行分析(王雷等,2009;王雷等,2010;Wang et al.,2012),發(fā)現(xiàn)了不同的放牧強度在季節(jié)尺度上都顯著減少了地氣間的實際蒸散(圖3),主要原因是:(1)放牧使地表反照率增大,從而減小了有效能量;(2)放牧導致土壤濕度的降低,減少了有效能量分配為潛熱通量的比例(Wang et al.,2012)。此外,Gao et al.(2009a)對內(nèi)蒙古典型草原下墊面地氣之間能量和物質(zhì)交換特征進行了研究,發(fā)現(xiàn)了該地區(qū)近地層能量平衡四分量具有顯著的日和季節(jié)變化特征,能量的不平衡問題主要是由于能量平衡四分量的不同相位造成的。

    圖2 2006 年5月1日至7月31日南亞夏季風指數(shù)變化曲線(a),季風中斷(b)和活躍期間(c)南亞地區(qū)500 hPa天氣形勢變化。星號代表絨布河谷的位置。(引自Zou et al.,2009)Fig.2 Variations of the South Asia summer monsoon index from May 1 to July 31 (a), 500 hPa streamlines in the SASM break period (b), and 500 hPa streamlines in the SASM active period (c).The star denotes the Rongbuk Valley.(From Zou et al., 2009)

    3.3 城市邊界層觀測實驗

    近十年來,隨著城市化發(fā)展加快,北京城市規(guī)模不斷擴大,城市建筑物逐漸向高空發(fā)展。城市下墊面變得更加復雜,動力學粗糙度日益增大。下墊面熱力學和動力學性質(zhì)的改變使該地區(qū)下墊面與大氣進行能量和物質(zhì)交換的特征發(fā)生變化,從而使大氣邊界層對天氣和氣候產(chǎn)生不可忽略的影響。針對城市邊界層理論研究,大氣物理研究所早在上世紀70年代起建立了北京325 m氣象塔(簡稱氣象塔)。目前,該氣象塔已經(jīng)積累了二十多年城市邊界層氣象數(shù)據(jù),為城市邊界層的研究提供寶貴數(shù)據(jù)。為了獲取城市冠層不同高度處大氣湍流脈動數(shù)據(jù),該氣象塔從2000年開始在47 m、120 m和280 m高度處增加了三層超聲風速溫度儀,為城市邊界層結(jié)構和湍流交換特征研究提供連續(xù)的湍流觀測資料。

    在城市邊界層湍流通量研究方面,通過分析氣象塔的渦動相關資料,發(fā)現(xiàn)了 2008年奧運期間CO2通量在量值上明顯小于過去同期的 CO2通量(Song and Wang,2012),分析結(jié)果證實了奧運期間所采取的車輛限行等政策取得了較顯著的效果。另外,通過分析多年(2006~2009)CO2通量資料(Liu et al.,2012b),發(fā)現(xiàn)了CO2通量在量值上與車流量關系密切,CO2通量日變化曲線的兩個峰值出現(xiàn)在上下班期間車流量最多的時次。結(jié)合氣象塔風場資料,發(fā)現(xiàn)了CO2通量最大值出現(xiàn)的方位在兩條車流量最多的公路交叉處。同時通過對比2006~2009年CO2排放總量,發(fā)現(xiàn)了2008年CO2排放總量明顯小于其他年份(表1),分析結(jié)果可為減少溫室氣體的排放提供參考依據(jù)。

    圖3 2005和2006年生長季(a)凈輻射日平均值,(b)LE/LEeq,(c)實際蒸散(ET)和潛在蒸散(PET)的季節(jié)變化。(引自Wang et al.,2012)Fig.3 Seasonal variation of (a) daily net radiation (Rn), (b) LE/LEeq, (c) ET and PET of the four sites during the growing seasons of 2005 and 2006.(From Wang et al., 2012)

    在城市邊界層湍流統(tǒng)計特征研究方面,利用氣象塔的渦動相關資料,分析了城市邊界層的陣風和湍流積分尺度以及湍流的垂直結(jié)構(王丙蘭,2012),發(fā)現(xiàn)了陣風因子的大小與陣風持續(xù)時間、風速平均時間和平均風速的大小有關,陣風因子隨陣風持續(xù)時間的增大而減小,隨風速平均時間的增大而增大,隨平均風速的增大而減小,當平均風速足夠大時,陣風因子為常數(shù);湍流積分尺度隨樣本長度的增大而增大;各種穩(wěn)定度條件下的上沖下掃運動對動量通量的貢獻不同,在城市冠層內(nèi)(47 m),下掃運動對通量的貢獻要大于上沖運動,但是在140 m和280 m高度處,均是上沖運動的貢獻略大于下掃運動。

    表1 2006~2009年北京325 m氣象塔觀測的CO2排放總量估算值(引自Liu et al.,2012b)Table 1 Total annual CO2 flux estimated using EC data from Beijing 325 m meteorological tower (2006–2009)(From Liu et al., 2012b)

    3.4 海洋邊界層觀測實驗

    在海洋下墊面大氣邊界層研究方面,LAPC科研人員利用 1998年在北冰洋地區(qū)大氣邊界層觀測實驗(SHEBA)獲取的觀測資料,研究了該地區(qū)大氣邊界層高度的演變特征并取得了一些有意義的發(fā)現(xiàn)(Dai et al.,2011)。由于下墊面熱力學和動力學特征的差異,北冰洋地區(qū)的大氣邊界層結(jié)構特征有別于中低緯度的大氣邊界層。由于觀測資料十分缺乏,該地區(qū)的大氣邊界層演變特征尚不清楚。通過分析該實驗期間獲取的飛機和GPS探空數(shù)據(jù)(圖4),發(fā)現(xiàn)了:溫度梯度方法較適用于估算該地區(qū)的大氣邊界層高度,并且給出了穩(wěn)定邊界層和對流邊界層的閾值(分別為6.5 K/100 m 和 1.0 K/100 m);該地區(qū)大氣邊界層高度春季大于夏季,最大高度出現(xiàn)在5月份(1150 m)。這些研究結(jié)果在一定程度上彌補了高緯度海洋下墊面大氣邊界層研究的空白。

    4 大氣邊界層數(shù)值模擬

    近地層湍流通量參數(shù)化方案是氣候和中尺度數(shù)值模擬研究中一個十分重要的課題。數(shù)值模式中湍流通量參數(shù)化一般有兩種計算方案,即迭代方案和非迭代方案。迭代方案完整地保留了 Monin-Obukhov相似理論,但迭代過程需要耗費大量的CPU時間。為了讓數(shù)值模式的通量計算方案具備更高的計算精度,同時花費更少的 CPU時間,可采用非迭代方案計算湍流通量。非迭代方案是目前湍流通量參數(shù)化計算方案研究的熱點?;诙嘣貧w分析方法,發(fā)展了一套湍流通量參數(shù)化方案的非迭代方法(李煜斌等,2009;Li et al.,2010)。新方法具有計算量少、高精度和允許 z0/z0h為大值等優(yōu)點。

    氣候模式幾乎都是基于能量平衡的原理,近地面可利用能量,即凈輻射與土壤熱通量之差(有時還考慮植被冠層熱儲存等),應與所分配的感熱通量和潛熱通量之和相等。然而,觀測實驗中普遍存在能量不平衡的問題,此問題在上世紀 90年代的陸面過程觀測實驗中被證實并得到國內(nèi)外科學家們的重視。Gao et al.(2010)分析了青藏高原安多站的觀測資料,發(fā)現(xiàn)土壤表面溫度與土壤熱通量日變化曲線的波相位差存在0至π/4的差異,此相位差對地表能量平衡的影響十分顯著。

    在海面湍流參數(shù)化研究方面,發(fā)展了一個可以用于氣候模式的海面湍流通量參數(shù)化方案(Ban et al.,2010)。該方案能夠給出更接近觀測事實的通量模擬結(jié)果。將新的湍流通量參數(shù)化方案應用到氣候模式CAM3中,不僅使得CAM3給出了更為精確的模擬結(jié)果,而且提高了CAM3的模擬時效(圖5)。此外,Gao et al.(2009b)在海/湖面摩擦系數(shù)、粗糙度和拖曳系數(shù)參數(shù)化方案的改進方面做了一些研究,并利用實際觀測資料評估了新的參數(shù)化方案。新的參數(shù)化方案使模擬結(jié)果更貼近真實值。

    在風場模擬和風能資源評估方面,利用擬合對數(shù)風廓線法和Davenport土地類型系統(tǒng)劃分法,結(jié)合風能資源利用的特點,估算了內(nèi)蒙東部草原地表粗糙度(李軍等,2011)。分析結(jié)果表明:該地區(qū)地表粗糙度長度有明顯月份和季節(jié)變化,夏秋植被茂盛期和冬春植被枯萎期粗糙度長度分別為0.138 m和 0.088 m;細微的地形起伏也會對粗糙度的評估造成影響;兩類方法估算的粗糙度長度大致相當,且估算的風能密度只相差2%左右。

    5 第二代超聲風速溫度儀觀測系統(tǒng)研制

    LAPC第二代超聲風速溫度儀(UAT-2)是在第一代的基礎上研發(fā)的新一代大氣湍流觀測系統(tǒng)。近年來,LAPC在利用多臺UAT-2構建湍流觀測系統(tǒng)方面取得了重要的進展。該系統(tǒng)可實時、同步完成對多點三維風速和溫度的脈動測量。系統(tǒng)擁有以下特點:(1)采樣頻率可達100 Hz。采樣頻率是超聲風速溫度儀的主要技術指標之一。UAT-2觀測系統(tǒng)最高采樣頻率提高至100 Hz。采樣頻率的提高有兩點重要意義:儀器可以在單位時間內(nèi)獲取更多測風數(shù)據(jù),能更好地、真實反映大氣中三維風速和溫度脈動值;采樣頻率的提高,反映儀器整體設計水平和工藝,以及相關最新技術的應用。(2)GPS信號的引用和實時同步組網(wǎng)觀測。隨著超聲風速溫度儀應用的深入與發(fā)展,需要多臺超聲實時、同步組網(wǎng)觀測。只有多臺超聲采用高精度同源時鐘,才可以保證此關鍵技術的實現(xiàn)。為此,系統(tǒng)引進了GPS秒信號作為網(wǎng)絡的共同時基。GPS秒信號的引進有兩個重要意義:可以把超聲有線局域網(wǎng)實時、同步觀測的理念拓展至無線廣域網(wǎng)的實時、同步觀測;在大大降低組網(wǎng)成本的同時,可以提高系統(tǒng)的采樣頻率。(3)聲陣實時垂直度的獲取。超聲測量對其傳感器——聲陣的垂直度有嚴格的要求,特別是在各種通量測量中,更突顯其重要性。為此,本系統(tǒng)加裝了專門測量聲陣動態(tài)垂直度的傾角儀子系統(tǒng)。子系統(tǒng)按100 Hz的采樣頻率采集動態(tài)垂直度數(shù)據(jù),保證了與超聲數(shù)據(jù)的同步,便于數(shù)據(jù)后處理時對相關數(shù)據(jù)進行必要的訂正,提高數(shù)據(jù)的準確性。圖 6為六臺UAT-2平行對比試驗結(jié)果。多臺UAT-2在采樣頻率為100 Hz時,觀測的脈動值具有很好的一致性,表明該系統(tǒng)達到國際同行的先進水平。

    圖5 CAM3中的參數(shù)化方案和新的參數(shù)化方案(GLL)計算的海面湍流通量與實測資料的對比:(a)風切應力;(b)感熱通量和(c)潛熱通量。紅線和藍線代表回歸直線;k1和k2分別為CAM3方案和GLL方案與觀測值的比率。(引自Ban et al.,2010)Fig.5 Fluxes calculated with the flux subroutine of the original CAM3 and by the GLL scheme versus measurements: ID wind stress (a), covariance sensible heat flux (b), and covariance latent heat flux (c).The red and blue lines represent regression lines; k1 and k2 are the slopes of fluxes calculated by CAM3 and GLL versus measurement.(From Ban et al., 2010)

    圖6 2011年8月6臺UAT-2超聲風速溫度儀平行對比實驗中觀測的垂直風速脈動(w ')(羅衛(wèi)東提供)Fig 6 Vertical wind speed fluctuation (w') were compared with 6 ultrasonic anemometers on 6 August, 2011 (From Luo Weidong, personnel communication)

    6 結(jié)論和展望

    本文從理論研究、觀測實驗、數(shù)值模擬和儀器研制四個方面,扼要總結(jié)了近4年來中國科學院大氣物理研究所大氣邊界層和大氣化學國家重點實驗室(LAPC)在大氣邊界層物理領域的主要研究進展。由于水平的限制,可能有不少不妥之處或不全面的地方。作者認為未來一段時間內(nèi)大氣邊界層物理研究重點仍然集中在非均勻下墊面(復雜地形,例如高大山地、城市等)大氣邊界層結(jié)構和特征、天氣和氣候模式中大氣邊界層參數(shù)化方案的改進、強風條件下海洋大氣邊界層參數(shù)化的方案及大型風電場風電量的短期預報系統(tǒng)等。非均勻下墊面大氣邊界層的定量描述是目前研究的重點和難點。此外,因為地面單點觀測通量的代表性通常在公里尺度,而中尺度模式和遙感資料反演得到的面上湍流通量通常尺度為幾公里到幾十公里。如何將單點觀測到的湍流通量升尺度到面上的通量,或者將遙感反演得到的面上湍流通量降尺度到單點,與地面單點的觀測結(jié)果比較,從而校驗和改進遙感反演算法。因此,地面單點觀測、模式輸出以及遙感反演結(jié)果的匹配也是大氣邊界層物理未來幾年研究的重點。

    References)

    Ban J M, Gao Z Q, Lenschow D H.2010.Climate simulations with a new air–sea turbulent flux parameterization in the National Center for Atmospheric Research Community Atmosphere Model (CAM3) [J].J.Geophys.Res., 115: D01106, doi: 10.1029/2009JD012802.

    Cheng Xueling, Zeng Qingcun, Hu Fei.2011.Characteristics of gusty wind disturbances and turbulent fluctuations in windy atmospheric boundary layer behind cold fronts [J].J.Geophys.Res., 116: D06101, doi:10.1029/2010JD015081.

    Cheng Xueling, Zeng Qingcun, Hu Fei.2012a.Stochastic modeling the effect of wind gust on dust entrainment during sand storm [J].Chinese Science Bulletin, 57 (27): 3595–3602, doi: 10.1007/s11434-012-5230-z.Cheng Xueling, Zeng Qingcun, Hu Fei.2012b.Parameterizations of some important characteristics of turbulent fluctuations and gusty wind disturbances in the atmospheric boundary layer [J].J.Geophys.Res., 117:D08113, doi: 10.1029/2011JD017191.

    Cheng Xueling, Hu Fei, Zeng Qingcun.2012c.Simulation of wind gust structure in the atmospheric boundary layer with Lattice Boltzmann Method [J].Chinese Science Bulletin, 57: 1196–1203.

    代成穎.2012.大氣邊界層頂部夾卷層特征及邊界層高度研究 [D].中國科學院大氣物理研究所博士學位論文, 129pp.Dai Chengying.2012.Properties of the entrainment zone atop the atmospheric boundary layer and determination of the boundary layer height [D].Ph.D.dissertation (in Chinese), Institute of Atmospheric Physics, Chinese Academy of Sciences, 129pp.

    Dai Chengying, Gao Zhiqiu, Wang Qing, et al.2011.Analysis ofatmospheric boundary layer height characteristic over the Arctic Ocean using the aircraft and GPS sounding [J].Atmospheric and Oceanic Science Letters, 4 (2): 124–130.

    馮健武, 劉輝志, 王雷, 等.2012.半干旱區(qū)不同下墊面地表粗糙度和湍流通量整體輸送系數(shù)變化特征 [J].中國科學 (地球科學), 42 (1):24–33.Feng Jianwu, Liu Huizhi, Wang Lei, et al.2012.Seasonal and inter-annual variation of surface roughness length and bulk transfer coefficients in a semiarid area [J].Science China (Earth Science), 55 (2):254–261.

    Gao Z Q, Lenschow D H, He Z, et al.2009a.Seasonal and diurnal variations in moisture, heat and CO2fluxes over a typical steppe prairie in Inner Mongolia, China [J].Hydrology and Earth System Sciences, 13:987–998.

    Gao Zhiqiu, Wang Qing, Zhou Mingyu.2009b.Wave-dependence of friction velocity, roughness length, and drag coefficient over coastal and open water surfaces by using three databases [J].Advance in Atmospheric Sciences, 26 (5): 887–894.

    Gao Z Q, Horton R, Liu H P.2010.Impact of wave phase difference between soil surface heat flux and soil surface temperature on soil surface energy balance closure [J].J.Geophys.Res., 115: D16112: doi:10.1029/2009JD013278.

    胡非, 洪鐘祥, 雷孝恩.2003.大氣邊界層和大氣環(huán)境研究進展 [J].大氣科學, 27 (4): 712–728.Hu Fei, Hong Zhongxiang, Lei Xiaoen.2003.Recent progress of atmospheric boundary layer physics and atmospheric environment research in IAP [J].Chinese Journal of Atmospheric Sciences (in Chinese), 27 (4): 712–728.

    Li Y B, Gao Z Q, Lenschow D H, et al.2010.An improved approach for parameterizing surface-layer turbulent transfer coefficients in numerical models [J].Bound.-Layer Meteor., 137: 153–165.

    Liu Huizhi, Sang Jianguo.2009.Analytical model of roll vortices in the convective boundary layer [J].Bound.-Layer Meteor., 130: 43–55.

    Liu Huizhi, Feng Jianwu.2012a.Seasonal and interannual variations of evapotranspiration and energy exchange over different land surfaces in a semi-arid area of China [J].Journal of Applied Meteorology and Climatology, 51 (10): 1875–1888.

    Liu Huizhi, Feng Jianwu, J?rvi L, et al.2012b.Four-year (2006–2009) eddy covariance measurements of CO2flux over an urban area in Beijing [J].Atmospheric Chemistry and Physics, 12: 1–87, doi: 10.5194/acp-12–1-2012.

    Liu Lei, Hu Fei, Cheng Xueling.2010.Probability density functions of velocity increments in the atmospheric boundary layer [J].Bound.-Layer Meteor., 134: 243–255.

    李軍, 胡非, 劉磊, 等.2011.風能資源評估中地表粗糙度的研究 [J].資源科學, 2011, 33 (12): 2341–2348.Li Jun, Hu Fei, Liu Lei, et al.2011.On estimation of surface roughness for wind energy resources assessment[J].Resources Science (in Chinese), 33 (12): 2341–2348.

    李煜斌, 高志球, 袁仁民, 等.2009.湍流通量參數(shù)化方案的非迭代方法研究 [J].大氣科學, 33 (4): 760–770.Li Yubin, Gao Zhiqiu, Yuan Renmin, et al.2009.Non-iteration methods of turbulent flux parameterization [J].Chinese Journal of Atmospheric Sciences (in Chinese), 33 (4): 760–770.

    Ma Shupo, Zhou Libo, Zou Han, et al.2009.Preliminary estimation of moisture exchange in rongbuk valley on the northern slope of Mt.Qomolangma [J].Atmospheric and Oceanic Science Letters, 2 (1): 40–44.

    Song Tao, Wang Yuesi.2012.Carbon dioxide fluxes from an urban area in Beijing [J].Atmospheric Research, 106: 139–149.

    涂剛, 劉輝志, 董文杰.2009.半干旱區(qū)不同下墊面近地層湍流通量特征分析 [J].大氣科學, 33 (4): 719–725.Tu Gang, Liu Huizhi, Dong Wenjie.2009.Characteristics of the surface turbulent fluxes over degraded grassland and cropland in the semi-arid area [J].Chinese Journal of Atmospheric Sciences (in Chinese), 33 (4): 719–725.

    王雷, 劉輝志, Ketzer B, 等.2009.放牧強度對內(nèi)蒙古半干旱草原地氣間能量和物質(zhì)交換的影響 [J].大氣科學, 33 (6): 120–1211.Wang Lei, Liu Huizhi, Ketzer B, et al.2009.Influence of grazing intensity on energy and mass exchange between the surface and atmosphere over semi-arid grassland in Inner Mongolia [J].Chinese Journal of Atmospheric Sciences (in Chinese), 33 (6): 1201–1211.

    王雷, 劉輝志, Schaffrath D, 等.2010.內(nèi)蒙古羊草草原和大針茅草原下墊面水汽、二氧化碳通量輸送特征 [J].高原氣象, 2010, 29 (3):605–613.Wang Lei, Liu Huizhi, Schaffrath D, et al.2010.Transfer characteristics of water vapor and CO2flux over leymus Chinensis and stip a grandis steppe surfaces in Inner Mongolia [J].Plateau Meteorology,29 (3): 605–613.

    Wang Lei, Liu Huizhi, Ketzer B, et al.2012.Effect of grazing intensity on evapotranspiration in the semiarid grasslands of Inner Mongolia, China[J].Journal of Arid Environments, 83: 15–24.

    王丙蘭.2012.城市邊界層湍流統(tǒng)計特征和環(huán)流結(jié)構研究 [D].中國科學院大氣物理研究所博士學位論文, 146pp.Wang Binglan.2012.Statistic characteristics of urban boundary layer turbulence and circulation structure [D].Ph.D.dissertation (in Chinese), Institute of Atmospheric Physics, Chinese Academy of Sciences, 146pp.

    Zeng Qingcun, Cheng Xueling, Hu Fei.2010.Gustiness and coherent structure of strong winds and their role in dust emission and entrainment[J].Advance in Atmospheric Sciences, 27 (1): 1–13.

    張美根, 胡非, 鄒捍, 等.2008.大氣邊界層物理與大氣環(huán)境過程研究進展 [J].大氣科學, 32 (4): 923–934.Zhang Meigen, Hu Fei, Zou Han,et al.2008.An overview of recent studies on atmospheric boundary layer physics and atmospheric environment in LAPC [J].Chinese Journal of Atmospheric Sciences (in Chinese), 32 (4): 923–934.

    周立波, 鄒捍, 馬舒坡, 等.2010.喜馬拉雅山地區(qū)地氣間物質(zhì)交換及其與南亞夏季風的聯(lián)系 [J].氣候與環(huán)境研究, 15 (3): 289–294.Zhou Libo, Zou Han, Ma Shupo, et al.2010.The air mass exchange between the surface and free atmosphere and its relationship with the South Asian summer monsoon [J].Climatic and Environmental Research (in Chinese),15 (3): 289–294.

    Zhou Libo, Zou Han, Ma Shupo, et al.2009.Quantitative estimation of air mass exchange by along-valley wind in the rongbuk valley[J].Atmospheric and Oceanic Science Letters, 2 (1): 29–34.

    Zhou Libo, Zou Han, Ma Shupo, et al.2011.Vertical air mass exchange driven by the local circulation on the northern slope of Mount Everest [J].Advance in Atmospheric Science, 28 (1): 217–222.

    Zou Han, Ma Shupo, Zhou Libo, et al.2009.Measured turbulent heat transfer on the northern slope of Mt.Everest and its relation to the south Asian summer monsoon [J].J.Geophys.Res., 36: L09810, doi:10.1029/2008GL036984.

    猜你喜歡
    下墊面陣風邊界層
    不同下墊面對氣溫的影響
    陣風戰(zhàn)斗機
    軍事文摘(2021年19期)2021-10-10 13:29:36
    法國陣風戰(zhàn)斗機
    軍事文摘(2021年17期)2021-09-24 01:11:46
    基于HIFiRE-2超燃發(fā)動機內(nèi)流道的激波邊界層干擾分析
    北京與成都城市下墊面閃電時空分布特征對比研究
    陣風勁吹
    航空世界(2018年12期)2018-07-16 08:34:50
    流域下墊面變化對潮白河密云水庫上游徑流影響分析
    下墊面變化對徑流及洪水影響分析
    一類具有邊界層性質(zhì)的二次奇攝動邊值問題
    非特征邊界的MHD方程的邊界層
    看十八女毛片水多多多| 国产精品99久久99久久久不卡 | 最近中文字幕高清免费大全6| videossex国产| 免费在线观看成人毛片| 日韩视频在线欧美| 国产熟女欧美一区二区| 亚洲精品乱久久久久久| 黄色毛片三级朝国网站 | 国产一区有黄有色的免费视频| 国产欧美日韩一区二区三区在线 | 久久精品国产鲁丝片午夜精品| 欧美成人午夜免费资源| 黑人巨大精品欧美一区二区蜜桃 | 国产精品人妻久久久影院| 免费不卡的大黄色大毛片视频在线观看| 日韩中文字幕视频在线看片| 夫妻午夜视频| 久久久久久久精品精品| 亚洲人与动物交配视频| 国产欧美日韩精品一区二区| 日韩人妻高清精品专区| 国产男女内射视频| 日本爱情动作片www.在线观看| 少妇人妻久久综合中文| 久久久久久久国产电影| 人妻系列 视频| 久久精品久久久久久久性| 亚洲精品国产av蜜桃| 新久久久久国产一级毛片| 王馨瑶露胸无遮挡在线观看| 午夜老司机福利剧场| 在现免费观看毛片| 少妇人妻精品综合一区二区| 国产69精品久久久久777片| 久久精品国产亚洲av天美| 免费看光身美女| 亚洲婷婷狠狠爱综合网| 性色av一级| 日韩成人伦理影院| 免费人妻精品一区二区三区视频| 亚洲色图综合在线观看| 乱人伦中国视频| 国产欧美日韩精品一区二区| 国产成人一区二区在线| 亚洲欧洲精品一区二区精品久久久 | 看十八女毛片水多多多| 亚洲第一区二区三区不卡| 精品熟女少妇av免费看| 搡女人真爽免费视频火全软件| 精品久久久久久久久av| 六月丁香七月| 亚洲精品第二区| 精品亚洲成a人片在线观看| 人人妻人人爽人人添夜夜欢视频 | 一边亲一边摸免费视频| 最后的刺客免费高清国语| 麻豆精品久久久久久蜜桃| 涩涩av久久男人的天堂| 国产精品人妻久久久久久| 免费av不卡在线播放| 精品人妻偷拍中文字幕| 久久国产乱子免费精品| 男女边吃奶边做爰视频| 丰满少妇做爰视频| 国产精品一二三区在线看| 国产91av在线免费观看| 国产色爽女视频免费观看| 在线播放无遮挡| 国产伦精品一区二区三区四那| h视频一区二区三区| 日韩电影二区| 亚洲成人一二三区av| 久久久久久久大尺度免费视频| 国产精品久久久久久久久免| 男人狂女人下面高潮的视频| 中文字幕免费在线视频6| 国产乱人偷精品视频| 免费黄频网站在线观看国产| 在线观看一区二区三区激情| 精品一区二区三卡| 亚洲无线观看免费| av有码第一页| 在线看a的网站| 久久99精品国语久久久| 黑人猛操日本美女一级片| 午夜老司机福利剧场| 伊人久久国产一区二区| 国产亚洲av片在线观看秒播厂| 国产成人精品一,二区| 国产伦精品一区二区三区四那| 久久精品国产亚洲av涩爱| 亚洲久久久国产精品| 亚洲第一av免费看| 少妇的逼水好多| 免费av不卡在线播放| 天天躁夜夜躁狠狠久久av| 久久精品国产自在天天线| 国产男女内射视频| 亚洲高清免费不卡视频| 99热6这里只有精品| 日日啪夜夜爽| 高清视频免费观看一区二区| 高清黄色对白视频在线免费看 | 三级国产精品片| 国产在线男女| 日韩伦理黄色片| 亚洲高清免费不卡视频| 午夜影院在线不卡| 国产精品国产三级专区第一集| 国产精品女同一区二区软件| 国产91av在线免费观看| 黄色一级大片看看| 人体艺术视频欧美日本| 毛片一级片免费看久久久久| 一本大道久久a久久精品| 韩国高清视频一区二区三区| av在线app专区| 国产有黄有色有爽视频| 丰满人妻一区二区三区视频av| 我要看黄色一级片免费的| 国产乱人偷精品视频| 午夜免费观看性视频| 日韩成人伦理影院| 成人毛片60女人毛片免费| 最黄视频免费看| 激情五月婷婷亚洲| 人妻人人澡人人爽人人| 欧美成人精品欧美一级黄| 国产精品国产av在线观看| 夫妻性生交免费视频一级片| 欧美精品一区二区免费开放| 国国产精品蜜臀av免费| 国产精品嫩草影院av在线观看| 99热这里只有精品一区| 久久久欧美国产精品| 大话2 男鬼变身卡| 亚洲欧美精品自产自拍| 色婷婷久久久亚洲欧美| 韩国av在线不卡| tube8黄色片| 欧美激情极品国产一区二区三区 | 人妻少妇偷人精品九色| 日本91视频免费播放| 精品人妻熟女毛片av久久网站| 成人漫画全彩无遮挡| 色婷婷久久久亚洲欧美| 久久久久久久亚洲中文字幕| 少妇高潮的动态图| 各种免费的搞黄视频| 国产精品一区www在线观看| 91精品一卡2卡3卡4卡| 国产欧美日韩综合在线一区二区 | 久久久国产精品麻豆| 99久久中文字幕三级久久日本| 人妻少妇偷人精品九色| av福利片在线观看| 日韩强制内射视频| 新久久久久国产一级毛片| 国产永久视频网站| 美女xxoo啪啪120秒动态图| 女性被躁到高潮视频| 色视频www国产| 99久久精品一区二区三区| 22中文网久久字幕| 美女福利国产在线| 久久毛片免费看一区二区三区| 最黄视频免费看| 欧美 亚洲 国产 日韩一| 赤兔流量卡办理| 99热国产这里只有精品6| 我要看黄色一级片免费的| 蜜桃在线观看..| 日韩精品有码人妻一区| 亚洲精品成人av观看孕妇| 精品酒店卫生间| 亚洲av福利一区| 欧美xxxx性猛交bbbb| 少妇高潮的动态图| 啦啦啦在线观看免费高清www| 国产69精品久久久久777片| 黑丝袜美女国产一区| 高清av免费在线| 亚洲国产欧美日韩在线播放 | 国产中年淑女户外野战色| 日韩,欧美,国产一区二区三区| 日韩三级伦理在线观看| 性色av一级| 国产有黄有色有爽视频| 国产av国产精品国产| 亚洲第一av免费看| 日本wwww免费看| 国产精品久久久久久精品电影小说| 熟女av电影| 久久久久久久久大av| 精品人妻熟女av久视频| 欧美日韩综合久久久久久| 免费久久久久久久精品成人欧美视频 | 欧美精品一区二区免费开放| 高清欧美精品videossex| 高清视频免费观看一区二区| √禁漫天堂资源中文www| 99热这里只有精品一区| 日韩一本色道免费dvd| 国精品久久久久久国模美| 在线观看一区二区三区激情| 精品国产乱码久久久久久小说| 中文字幕精品免费在线观看视频 | 两个人免费观看高清视频 | 99热这里只有是精品50| 亚洲av成人精品一二三区| 狂野欧美白嫩少妇大欣赏| 女人久久www免费人成看片| 国产精品三级大全| 三上悠亚av全集在线观看 | 国产精品久久久久久精品古装| av.在线天堂| 免费大片18禁| 国产午夜精品一二区理论片| 爱豆传媒免费全集在线观看| 亚洲图色成人| 观看av在线不卡| 一区二区三区乱码不卡18| 国产成人午夜福利电影在线观看| 在线观看三级黄色| 精品人妻一区二区三区麻豆| 国产欧美日韩综合在线一区二区 | av.在线天堂| 嫩草影院新地址| 成人综合一区亚洲| 美女国产视频在线观看| 亚洲国产欧美日韩在线播放 | 久热久热在线精品观看| 国产极品天堂在线| 成年人午夜在线观看视频| 黄色一级大片看看| 欧美日韩视频高清一区二区三区二| 亚洲精品乱码久久久v下载方式| 欧美激情极品国产一区二区三区 | 蜜臀久久99精品久久宅男| 一本色道久久久久久精品综合| 久久人人爽人人片av| 精品国产一区二区三区久久久樱花| 9色porny在线观看| 日韩一区二区视频免费看| 精品久久久久久电影网| 纯流量卡能插随身wifi吗| 天天躁夜夜躁狠狠久久av| 欧美日韩国产mv在线观看视频| 亚洲性久久影院| 国产男人的电影天堂91| 97在线人人人人妻| 久久精品久久久久久噜噜老黄| 99热6这里只有精品| 欧美三级亚洲精品| 大片免费播放器 马上看| 免费看不卡的av| 亚洲av.av天堂| 亚洲国产欧美日韩在线播放 | 极品教师在线视频| 交换朋友夫妻互换小说| 黄色毛片三级朝国网站 | 日韩三级伦理在线观看| 青青草视频在线视频观看| 菩萨蛮人人尽说江南好唐韦庄| 国产高清有码在线观看视频| 内地一区二区视频在线| 国产精品一区二区三区四区免费观看| av福利片在线| 久久午夜福利片| 国产色爽女视频免费观看| 国产女主播在线喷水免费视频网站| 激情五月婷婷亚洲| 国产69精品久久久久777片| av有码第一页| 亚州av有码| 少妇人妻精品综合一区二区| 尾随美女入室| 国内揄拍国产精品人妻在线| 插阴视频在线观看视频| a级一级毛片免费在线观看| 久久99蜜桃精品久久| 国产成人免费无遮挡视频| 黑人巨大精品欧美一区二区蜜桃 | 精华霜和精华液先用哪个| 日韩 亚洲 欧美在线| 你懂的网址亚洲精品在线观看| 中文字幕av电影在线播放| 国产精品国产av在线观看| 国产男人的电影天堂91| 亚洲不卡免费看| 伊人亚洲综合成人网| 日本av手机在线免费观看| 国产在线一区二区三区精| 一级毛片黄色毛片免费观看视频| 久久久精品94久久精品| 国产精品久久久久久久久免| 夫妻性生交免费视频一级片| 一级a做视频免费观看| 妹子高潮喷水视频| 黑丝袜美女国产一区| 久久久国产一区二区| 欧美xxxx性猛交bbbb| 免费少妇av软件| 日韩一本色道免费dvd| 人人妻人人添人人爽欧美一区卜| 久久久久久久国产电影| 日本-黄色视频高清免费观看| 91久久精品国产一区二区三区| 91精品国产九色| 人体艺术视频欧美日本| 丰满乱子伦码专区| 特大巨黑吊av在线直播| 免费观看av网站的网址| 搡老乐熟女国产| 99热这里只有精品一区| 国产永久视频网站| 能在线免费看毛片的网站| 国产欧美亚洲国产| 免费观看无遮挡的男女| 性高湖久久久久久久久免费观看| 国产成人a∨麻豆精品| 国产欧美亚洲国产| 欧美少妇被猛烈插入视频| 午夜精品国产一区二区电影| 伦精品一区二区三区| 国产精品人妻久久久久久| 男女边摸边吃奶| 毛片一级片免费看久久久久| 在线观看一区二区三区激情| 日韩av不卡免费在线播放| 多毛熟女@视频| 亚洲不卡免费看| 亚洲经典国产精华液单| 国产一区亚洲一区在线观看| 狂野欧美激情性bbbbbb| 国产白丝娇喘喷水9色精品| a级毛片免费高清观看在线播放| 国产91av在线免费观看| 免费av不卡在线播放| 国产精品人妻久久久影院| 国产精品一区www在线观看| 国产伦在线观看视频一区| 中国三级夫妇交换| 成人二区视频| 99视频精品全部免费 在线| 一级a做视频免费观看| 在线免费观看不下载黄p国产| 日韩一区二区视频免费看| 91aial.com中文字幕在线观看| 美女福利国产在线| 久久精品熟女亚洲av麻豆精品| 日韩中字成人| 日韩亚洲欧美综合| 国产伦在线观看视频一区| 18禁裸乳无遮挡动漫免费视频| 97精品久久久久久久久久精品| 欧美激情极品国产一区二区三区 | 久久精品国产亚洲av天美| 日本-黄色视频高清免费观看| 天堂8中文在线网| 日本欧美视频一区| 一区在线观看完整版| 免费看av在线观看网站| 日韩中字成人| 国产在线一区二区三区精| 亚洲国产av新网站| 丝袜喷水一区| 自拍欧美九色日韩亚洲蝌蚪91 | 丰满乱子伦码专区| 国产伦理片在线播放av一区| 9色porny在线观看| 日本黄色日本黄色录像| av卡一久久| 国产高清不卡午夜福利| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 乱系列少妇在线播放| 久久久亚洲精品成人影院| 久久精品夜色国产| 国产午夜精品久久久久久一区二区三区| 在线观看www视频免费| 久久99精品国语久久久| 秋霞在线观看毛片| 91午夜精品亚洲一区二区三区| 在线免费观看不下载黄p国产| 亚洲精品久久久久久婷婷小说| 乱系列少妇在线播放| 欧美三级亚洲精品| 最近中文字幕高清免费大全6| 大又大粗又爽又黄少妇毛片口| 国产一区二区在线观看av| 性色av一级| 亚洲欧美精品自产自拍| 在线观看免费高清a一片| 久久av网站| 国国产精品蜜臀av免费| 成人美女网站在线观看视频| 国产精品秋霞免费鲁丝片| a级毛色黄片| 91久久精品电影网| 精品卡一卡二卡四卡免费| 秋霞伦理黄片| 99久久精品热视频| 亚洲成人一二三区av| 精品人妻偷拍中文字幕| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 亚洲av男天堂| 嘟嘟电影网在线观看| 纵有疾风起免费观看全集完整版| 18禁在线无遮挡免费观看视频| 国产毛片在线视频| 精品久久国产蜜桃| 国产亚洲av片在线观看秒播厂| 婷婷色综合大香蕉| 秋霞在线观看毛片| 成人毛片60女人毛片免费| 成人无遮挡网站| 内射极品少妇av片p| 我要看日韩黄色一级片| 一级毛片aaaaaa免费看小| 色婷婷av一区二区三区视频| 女的被弄到高潮叫床怎么办| 国产熟女欧美一区二区| 日韩欧美 国产精品| 少妇被粗大猛烈的视频| 男人狂女人下面高潮的视频| 午夜av观看不卡| 国产 精品1| 欧美激情国产日韩精品一区| 国产国拍精品亚洲av在线观看| 伊人久久国产一区二区| www.色视频.com| av线在线观看网站| 日韩伦理黄色片| 国产欧美日韩一区二区三区在线 | 亚洲精品aⅴ在线观看| av不卡在线播放| 你懂的网址亚洲精品在线观看| av卡一久久| 欧美日韩亚洲高清精品| 久久 成人 亚洲| 男女边吃奶边做爰视频| 国产日韩欧美在线精品| 一本—道久久a久久精品蜜桃钙片| 香蕉精品网在线| 国产精品.久久久| 欧美日韩视频高清一区二区三区二| 日日摸夜夜添夜夜添av毛片| 美女xxoo啪啪120秒动态图| 欧美变态另类bdsm刘玥| 有码 亚洲区| 秋霞伦理黄片| 日本猛色少妇xxxxx猛交久久| 亚洲三级黄色毛片| 七月丁香在线播放| 熟妇人妻不卡中文字幕| 视频区图区小说| 久久精品国产a三级三级三级| 草草在线视频免费看| 香蕉精品网在线| 国产熟女午夜一区二区三区 | 性色av一级| 亚洲精品日韩在线中文字幕| 亚洲精品日本国产第一区| 三级国产精品欧美在线观看| 插阴视频在线观看视频| 九草在线视频观看| 少妇裸体淫交视频免费看高清| 麻豆成人av视频| 男女边吃奶边做爰视频| av在线app专区| 久久精品国产自在天天线| 嘟嘟电影网在线观看| 黄片无遮挡物在线观看| 亚洲成人手机| 九九在线视频观看精品| 日日摸夜夜添夜夜爱| 尾随美女入室| 欧美日韩视频高清一区二区三区二| kizo精华| 国产日韩一区二区三区精品不卡 | 国产精品一区二区在线不卡| 成年女人在线观看亚洲视频| 国产精品免费大片| 欧美bdsm另类| 777米奇影视久久| av黄色大香蕉| 日韩大片免费观看网站| 亚洲精品成人av观看孕妇| 久久免费观看电影| 在线观看免费高清a一片| av有码第一页| 最近手机中文字幕大全| 欧美老熟妇乱子伦牲交| 国产女主播在线喷水免费视频网站| 22中文网久久字幕| 久久久精品94久久精品| 性色av一级| 国产一区二区三区综合在线观看 | 18禁在线播放成人免费| 国产老妇伦熟女老妇高清| 成人综合一区亚洲| 精品久久久久久久久亚洲| 日韩欧美精品免费久久| 在线观看一区二区三区激情| 日韩一本色道免费dvd| 一区二区三区四区激情视频| 国产日韩欧美亚洲二区| freevideosex欧美| 丰满少妇做爰视频| 欧美老熟妇乱子伦牲交| 国产男人的电影天堂91| 丝袜脚勾引网站| 久久久久国产网址| 久久久国产一区二区| 简卡轻食公司| 人人妻人人澡人人爽人人夜夜| 夫妻午夜视频| 成人美女网站在线观看视频| 少妇丰满av| 熟女人妻精品中文字幕| 波野结衣二区三区在线| 青青草视频在线视频观看| 亚洲丝袜综合中文字幕| 人人妻人人澡人人看| 国产日韩欧美视频二区| 国产精品一区二区性色av| 成人漫画全彩无遮挡| 中文字幕久久专区| 日本色播在线视频| 国产精品一区www在线观看| 五月开心婷婷网| 亚洲av电影在线观看一区二区三区| videos熟女内射| 国产视频内射| 亚洲精品一二三| www.色视频.com| 欧美3d第一页| 国产免费一级a男人的天堂| videossex国产| 精品人妻熟女毛片av久久网站| 国产精品免费大片| 91成人精品电影| 日韩伦理黄色片| 亚洲av综合色区一区| 日韩人妻高清精品专区| 我要看黄色一级片免费的| 三级经典国产精品| 国产精品女同一区二区软件| 国产色婷婷99| 日日撸夜夜添| 国产91av在线免费观看| 国产精品成人在线| 熟女人妻精品中文字幕| 大话2 男鬼变身卡| 国模一区二区三区四区视频| 亚洲欧美日韩东京热| 国产精品一区二区在线不卡| 中文字幕久久专区| 久久久久视频综合| 蜜桃久久精品国产亚洲av| 欧美精品一区二区免费开放| 国产黄色视频一区二区在线观看| 亚洲电影在线观看av| 久久午夜福利片| 精品人妻熟女毛片av久久网站| 在线观看国产h片| 人妻 亚洲 视频| 国产男人的电影天堂91| 中文字幕久久专区| 少妇的逼好多水| 中文字幕久久专区| 各种免费的搞黄视频| 十八禁高潮呻吟视频 | 午夜福利视频精品| 十八禁网站网址无遮挡 | 午夜影院在线不卡| av免费在线看不卡| 极品少妇高潮喷水抽搐| 亚洲国产av新网站| 免费黄网站久久成人精品| 看非洲黑人一级黄片| 18禁动态无遮挡网站| 免费观看a级毛片全部| 国产精品久久久久久精品古装| 如日韩欧美国产精品一区二区三区 | 国产美女午夜福利| 伦理电影大哥的女人| 国产精品免费大片| 精品一区二区三卡| 日本wwww免费看| 久久久久国产精品人妻一区二区| 免费观看av网站的网址| 春色校园在线视频观看| 日本免费在线观看一区| 国产精品嫩草影院av在线观看| 免费观看av网站的网址| 我要看黄色一级片免费的| 亚洲欧美一区二区三区黑人 | 一级毛片我不卡| 边亲边吃奶的免费视频| 亚洲精品日本国产第一区| 亚洲av日韩在线播放| 亚洲精品日本国产第一区| 欧美精品一区二区大全| 国产淫片久久久久久久久| 纵有疾风起免费观看全集完整版| 日韩亚洲欧美综合| 亚洲av日韩在线播放| 水蜜桃什么品种好| 日韩av不卡免费在线播放| 亚洲欧美一区二区三区国产| 日韩av不卡免费在线播放|