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

    晨昏軌道微波溫度計(jì)資料同化對(duì)降水定量預(yù)報(bào)的影響及其對(duì)三軌衛(wèi)星系統(tǒng)的意義

    2016-03-23 00:34:33鄒曉蕾秦正坤翁富忠
    大氣科學(xué) 2016年1期

    鄒曉蕾秦正坤翁富忠

    1馬里蘭大學(xué)地球系統(tǒng)多學(xué)科合作研究中心,美國(guó)College Park 207402南京信息工程大學(xué)大氣科學(xué)學(xué)院,南京2100443美國(guó)國(guó)家海洋大氣管理局衛(wèi)星研究和應(yīng)用中心,美國(guó)College Park 20740

    ?

    晨昏軌道微波溫度計(jì)資料同化對(duì)降水定量預(yù)報(bào)的影響及其對(duì)三軌衛(wèi)星系統(tǒng)的意義

    鄒曉蕾1秦正坤2翁富忠3

    1馬里蘭大學(xué)地球系統(tǒng)多學(xué)科合作研究中心,美國(guó)College Park 20740
    2南京信息工程大學(xué)大氣科學(xué)學(xué)院,南京210044
    3美國(guó)國(guó)家海洋大氣管理局衛(wèi)星研究和應(yīng)用中心,美國(guó)College Park 20740

    摘 要極軌衛(wèi)星的高級(jí)微波溫度計(jì)(Advanced Microwave Sounding Unit-A,簡(jiǎn)稱AMSU-A)輻射資料對(duì)提高降水定量預(yù)報(bào)的水平有重要作用。但是極軌衛(wèi)星的軌道特征導(dǎo)致乘載其上的微波溫度計(jì)資料在區(qū)域同化系統(tǒng)中存在嚴(yán)重缺測(cè)。本研究重點(diǎn)分析了晨昏軌道衛(wèi)星上微波溫度計(jì)資料同化對(duì)墨西哥灣沿岸定量降水預(yù)報(bào)的重要影響。研究選取了早晨星NOAA-15、上午星MetOp-A和下午星NOAA-18,利用美國(guó)NCEP(National Centers for Environmental Prediction)的業(yè)務(wù)同化系GSI(Gridpoint Statistical Interpolation)資料同化系統(tǒng),進(jìn)行了加和不加NOAA-15 AMSU-A資料的兩組資料同化和預(yù)報(bào)試驗(yàn),來(lái)闡明晨昏軌道衛(wèi)星上微波溫度計(jì)資料同化對(duì)墨西哥灣沿岸降水預(yù)報(bào)的重要影響。試驗(yàn)結(jié)果分析表明如果僅同化NOAA-18和MetOp-A資料,在協(xié)調(diào)世界時(shí)00:00和12:00的同化時(shí)間,在墨西哥灣和美國(guó)西部大陸就是衛(wèi)星觀測(cè)資料缺測(cè)區(qū),而早晨星NOAA-15資料正好可以填補(bǔ)這個(gè)資料空缺。模式預(yù)報(bào)也表明,同化NOAA-15的AMSU-A資料可以對(duì)墨西哥灣降水有持續(xù)的正影響。這一研究證明了保持有搭載著AMSU-A或者相似儀器的早晨星,對(duì)區(qū)域降水預(yù)報(bào)的重要性。由于目前NOAA-15是唯一的一顆正在運(yùn)行的、已遠(yuǎn)超過(guò)其正常運(yùn)行期的早晨星,通過(guò)技術(shù)手段維持NOAA-15的AMSU-A儀器更超長(zhǎng)期運(yùn)行也就特別重要。

    關(guān)鍵詞資料同化 晨昏軌道 極軌衛(wèi)星 降水預(yù)報(bào)

    資助項(xiàng)目 江蘇省自然科學(xué)基金項(xiàng)目BK2011039,國(guó)家自然科學(xué)基金項(xiàng)目41475103

    Found by the Natural Science Foundation of Jiangsu Province (Grant BK2011039), National Natural Science Foundation of China (Grant 41475103)

    Impact of Dawn–Dusk Satellite AMSU-A Data on Quantitative Precipitation Forecasts and the Implications for Three-Orbit Constellation

    ZOU Xiaolei1, QIN Zhengkun2, and WENG Fuzhong3

    1 Earth System Science Interdisciplinary Center, University of Maryland, College Park 20740, USA
    2 College of Atmospheric Science, Nanjing University of Information Science and Technology, Nanjing 210044
    3 National Oceanic and Atmospheric Administration (NOAA) Center for Satellite Applications and Research, College Park 20740, USA

    Abstract This study demonstrates the importance of an early morning orbit satellite for improved coastal quantitative precipitation forecasts (QPFs) near the Gulf of Mexico. The Advanced Microwave Sounding Unit-A (AMSU-A) radiance observations from the early morning satellite NOAA-15, the mid-morning satellite MetOp-A, and the afternoon satellite NOAA-18, are assimilated using the National Centers for Environmental Prediction (NCEP) Gridpoint StatisticalInterpolation (GSI). For the Gulf coast precipitation case selected in this study, two pairs of data assimilation and forecasting experiments are carried out to compare the differences in QPFs with and without assimilating AMSU-A data from NOAA-15. It is shown that the two orbits provided by NOAA-18 and MetOp-A render both the Gulf of Mexico and the western continent of the United States as two data-void areas at 0000 UTC and 1200 UTC. The NOAA-15 orbit fills these data gaps. Adding NOAA-15 AMSU-A data into GSI data assimilation results in a consistently positive impact on the QPFs near the Gulf coast. It is thus suggested to have a continual availability of an early morning orbiting satellite with an AMSU-A or AMSU-A-like instrument onboard. A partial solution would be a sustained effort to maintain the NOAA-15 AMSU-A for longer-lived operation.

    Keywords Data assimilation, Dawn–dusk orbit, Polar orbit satellite, quantitative precipitation forecast

    1 引言

    極軌衛(wèi)星微波溫度計(jì)資料同化一直是數(shù)值天氣預(yù)報(bào)改進(jìn)的重要方法之一。自從1978年美國(guó)的電視和紅外輻射觀測(cè)衛(wèi)星TIROS-N(Television and InfraRed Observation Satellite)發(fā)射以來(lái),國(guó)際上發(fā)射了一系列搭載微波溫度計(jì)的極軌氣象衛(wèi)星。在1978~2004年的氣象衛(wèi)星發(fā)展早期,搭載在8個(gè)NOAA(National Oceanic and Atmospheric Administration)系列極軌衛(wèi)星(NOAA-6、NOAA-7、NOAA-8、NOAA-9、NOAA-10、NOAA-11、NOAA-12和NOAA-14)上的微波溫度探測(cè)計(jì)(Microwave Sounding Unit,簡(jiǎn)稱MSU)只有4個(gè)大氣溫度廓線探測(cè)通道。從1998年開始,AMSU-A (Advanced Microwave Sounding Unit-A)開始替代MSU,搭載在相繼發(fā)射的5個(gè)NOAA系列氣象衛(wèi)星(NOAA-15、NOAA-16、NOAA-17、NOAA-18 和NOAA-19)和兩個(gè)EUMESAT(European Organisation for the Exploitation of Meteorological Satellites)衛(wèi)星(Meteorological Operational Satellite Program of Europe-A/B,簡(jiǎn)稱MetOp-A和MetOp-B)。AMSU-A比MSU多11個(gè)通道,其中AMSU-A的通道3、5、7和9與MSU的4個(gè)通道頻率相似。AMSU-A另外還包括了3個(gè)窗區(qū)通道,通道1、2和15,其中心頻率分別位于23.8、31.4和89 GHz。窗區(qū)通道主要用于探測(cè)云和地表參數(shù)。額外增加的8個(gè)中高層通道也使得AMSU-A能夠更好的探測(cè)對(duì)流層和平流層的大氣溫度廓線特征。

    每個(gè)極軌衛(wèi)星都是沿太陽(yáng)同步軌道運(yùn)行,在800 km左右的高度圍繞地球旋轉(zhuǎn),所以極軌衛(wèi)星能夠提供一天兩次衛(wèi)星局地越赤道時(shí)刻的全球輻射觀測(cè)資料。一次全球資料一般由14個(gè)軌道資料組成,每個(gè)軌道資料耗時(shí)約為100 min。每個(gè)軌道有固定的局地越赤道時(shí)間,降軌越赤道時(shí)間為09:30(當(dāng)?shù)貢r(shí)間)左右的衛(wèi)星稱為上午星,例如NOAA-17、MetOp-A和MetOp-B等,升軌越赤道時(shí)間為下午01:30(當(dāng)?shù)貢r(shí)間)左右的衛(wèi)星稱為下午星,比如TIROS-N、NOAA-7、NOAA-9、NOAA-11、NOAA-14、NOAA-16、NOAA-18和NOAA-19,晨昏軌道極軌衛(wèi)星稱為早晨星,主要有NOAA-6、NOAA-8、NOAA-10、NOAA-12、NOAA-15,其降軌越赤道時(shí)間一般為當(dāng)?shù)貢r(shí)間05:30或者07:30。利用早晨星、上午星和下午星的共同觀測(cè),極軌氣象衛(wèi)星能夠以4 h間隔提供一天6次全球大氣溫度廓線觀測(cè)資料。但是目前NOAA-15極軌衛(wèi)星已經(jīng)是唯一的早晨軌道極軌衛(wèi)星,而且超長(zhǎng)時(shí)間運(yùn)行了15年以上,遠(yuǎn)遠(yuǎn)超過(guò)了衛(wèi)星的設(shè)計(jì)壽命。搭載在NOAA-15上的AMSU-A儀器,除了通道11和14之外,大部分通道仍然保持了良好的業(yè)務(wù)運(yùn)行狀態(tài)。除了NOAA-15之外,還有兩個(gè)上午星MetOp-A和MetOp-B,和兩個(gè)下午星NOAA-18和NOAA-19能夠?yàn)槎鄠€(gè)國(guó)家業(yè)務(wù)資料同化系統(tǒng)提供AMSU-A觀測(cè)資料。

    AMSU-A輻射資料同化多是利用三維變分或者四維變分資料同化系統(tǒng)進(jìn)行。變分系統(tǒng)中輻射資料的有效同化主要是得益于快速輻射傳輸模式(the fast Radiative Transfer Model,簡(jiǎn)稱RTM)的開發(fā)。基于給定的大氣狀態(tài)模擬廓線,RTM可以準(zhǔn)確模擬出AMSU-A儀器各個(gè)觀測(cè)頻率的輻射量(McMillin and Fleming, 1976; Saunders et al., 1999, 2007; Weng, 2007)。同化通過(guò)調(diào)整大氣模式模擬的大氣廓線,并作為RTM模式輸入量進(jìn)行模擬,從而使得RTM模擬的輻射量與觀測(cè)更為接近。這里的“接近”指的是觀測(cè)和模擬之間有最大相似度或誤差協(xié)方差最小。自從同化了極軌氣象衛(wèi)星的微波溫度計(jì)輻射資料后,全球數(shù)值天氣預(yù)報(bào)水平得到了快速穩(wěn)定的提高(Eyre et al., 1993; Andersson et al., 1994)。早在20世紀(jì)90年代,美國(guó)的數(shù)值天氣預(yù)報(bào)中心NCEP (National Centers for Environmental Prediction)和歐洲短期天氣預(yù)報(bào)中心ECMWF(European Center of Medium-Range Weather Forecasts)就實(shí)現(xiàn)了業(yè)務(wù)天氣預(yù)報(bào)中同化TIROS(Television Infrared Observation Satellite Program)的業(yè)務(wù)垂直探測(cè)計(jì)TOVS(Operational Vertical Sounder)的輻射資料,并顯著提高了全球數(shù)值天氣預(yù)報(bào)的業(yè)務(wù)水平。

    衛(wèi)星資料同化對(duì)全球數(shù)值預(yù)報(bào)的提高效果得到了普遍認(rèn)可,其中主要原因之一就是因?yàn)樾l(wèi)星資料的高空間覆蓋率,尤其是能夠?qū)Q蠛透呱降貐^(qū)等常規(guī)觀測(cè)資料缺乏的地區(qū)提供遙感觀測(cè)資料。然而,單顆極軌衛(wèi)星一天只能夠覆蓋地球兩次,而業(yè)務(wù)全球資料同化系統(tǒng)常規(guī)需要每天至少進(jìn)行四次同化(00:00、06:00、12:00和18:00,協(xié)調(diào)世界時(shí)),所以為了提高全球觀測(cè)資料的時(shí)間頻率, 在有上午星和下午星的基礎(chǔ)上,增加早晨星,能夠?yàn)槿蛸Y料同化系統(tǒng)提供以4 h間隔的一天6次的全球大氣溫度和水汽廓線觀測(cè)資料。三軌衛(wèi)星的任何一顆星缺測(cè)就會(huì)造成全球資料同化系統(tǒng)中觀測(cè)資料的缺失。極軌衛(wèi)星資料覆蓋率的影響在區(qū)域模式的資料同化中更為明顯。比如在6 h的同化窗口內(nèi),由于上午星和下午星在美國(guó)南部和東部海岸存在資料空白區(qū),這就容易導(dǎo)致對(duì)流天氣系統(tǒng)上游條件不能被很好的觀測(cè),進(jìn)而減弱了資料同化對(duì)區(qū)域短期(24~36 h)定量降水預(yù)報(bào)的改進(jìn)效果。由于國(guó)際上業(yè)務(wù)同化普遍使用的三種軌道極軌氣象衛(wèi)星中的早晨星僅剩下一顆超期服役的NOAA-15衛(wèi)星,為了進(jìn)一步明確早晨星的重要性,本研究以區(qū)域天氣預(yù)報(bào)為例,評(píng)估了同時(shí)同化早晨星、上午星和下午星的AMSU-A觀測(cè)資料和僅同化上午星和下午星AMSU-A資料對(duì)美國(guó)墨西哥灣沿岸降水定量預(yù)報(bào)的影響。研究利用NCEP的業(yè)務(wù)資料同化系統(tǒng)GSI同化系統(tǒng)(Wu et al., 2002; Purser et al., 2003a, 2003b)進(jìn)行,定量降水預(yù)報(bào)QPFs(Quantitative Precipitation Forecasts)是基于WRF-ARW(the Advanced Research WRF)模式。

    文章的結(jié)構(gòu)如下:第二節(jié)描述了AMSU-A觀測(cè)資料的主要特征,第三節(jié)簡(jiǎn)單介紹了WRF-ARW模式和GSI資料同化系統(tǒng),并描述了試驗(yàn)設(shè)計(jì)方法,第四節(jié)則是介紹了同化系統(tǒng)中的質(zhì)量控制方法、偏差訂正方法,還細(xì)致比較了同化前(O-B)和同化后(O-A)的AMSU-A觀測(cè)亮溫和模擬亮溫之間差值變化特征,第五節(jié)主要分析了同化NOAA-15的AMSU-A輻射資料對(duì)墨西哥灣沿岸降水預(yù)報(bào)的影響特征,主要是證明了NOAA-15的AMSU-A資料同化對(duì)NOAA-18、MetOp-A衛(wèi)星資料和常規(guī)資料的附加效果,第六節(jié)是總結(jié)和討論。

    2 試驗(yàn)設(shè)計(jì)

    2.1 資料同化系統(tǒng)和預(yù)報(bào)模式

    NCEP的業(yè)務(wù)資料同化系統(tǒng)GSI同化系統(tǒng)是一個(gè)三維變分同化系統(tǒng)(Wu et al., 2002; Purser et al., 2003a, 2003b)。通過(guò)引入遞歸濾波方法,GSI系統(tǒng)建立了格點(diǎn)上的背景誤差協(xié)方差,具有適應(yīng)實(shí)際情況、各項(xiàng)異性且非均一的特征,對(duì)中尺度數(shù)據(jù)同化和預(yù)報(bào)尤其重要,因?yàn)閷?duì)于中尺度系統(tǒng),無(wú)論是大氣狀態(tài)特征,還是觀測(cè)數(shù)據(jù)的密度和質(zhì)量都存在顯著的空間非均一性。GSI系統(tǒng)的用戶指南提供了關(guān)于系統(tǒng)安裝、編譯和運(yùn)行的詳細(xì)信息。目前GSI資料同化系統(tǒng)已經(jīng)成功的移植到南京信息工程大學(xué)資料同化研究與應(yīng)用中心的Linux平臺(tái)大型計(jì)算機(jī)中,所有的同化和預(yù)報(bào)試驗(yàn)結(jié)果都是依賴于該計(jì)算系統(tǒng)。

    GSI同化系統(tǒng)整合了通用輻射傳輸模式CRTM (the Community Radiative Transfer Model)(Weng, 2007; Han et al., 2007)。該模式在晴空條件下能夠準(zhǔn)確模擬AMSU-A的所有通道輻射亮溫。但是在有云情況下模式的模擬能力還存在一定的不足,所以GSI系統(tǒng)主要同化晴空條件下的AMSU-A資料。晴空條件下輻射資料的模擬需要大氣溫度廓線和水汽廓線,以及一些地表參數(shù)(比如地表溫度和表面風(fēng)速等),另外還有衛(wèi)星觀測(cè)資料的幾何參數(shù)(比如方位角、天頂角、大氣透過(guò)率等)。CRTM的Jacobian模塊也被包含在GSI同化系統(tǒng)中,從而能夠有效的計(jì)算模擬亮溫的梯度變化信息,實(shí)現(xiàn)衛(wèi)星輻射資料的快速同化。

    研究所用的WRF-ARW模式為3.0版本。模式水平分辨率為10 km,垂直為50層,模式層頂設(shè)在10 hPa,模擬區(qū)域的格點(diǎn)數(shù)為250×200×27。模式的初始條件和邊界條件都是來(lái)自于美國(guó)NCEP的1°×1°業(yè)務(wù)全球分析資料FNL(Final)。WRF模式選用了單通量三級(jí)微物理過(guò)程(Hong and Lim, 2006),Kain-Fritsch積云參數(shù)化方案(Kain and Fritsch, 1990; 1993; 2004)以及Yonsei 行星邊界層算法(Hong and Dudhia, 2003)。

    2.2 試驗(yàn)方案設(shè)計(jì)

    圖1給出了1998~2014年NOAA-15、NOAA-18、NOAA-19和MetOp-A衛(wèi)星的升軌和降軌的LECT(Local Equatorial Crossing Time)時(shí)間。從LECT時(shí)間的比較中可以看到,NOAA-15、MetOp-A 和NOAA-18這三顆衛(wèi)星就可以組成很好的早晨星、上午星和下午星的極軌衛(wèi)星觀測(cè)系統(tǒng)。在2008年,它們的越赤道時(shí)間分別為當(dāng)?shù)貢r(shí)間的早晨05:00、上午09:30和下午02:00,三軌極軌衛(wèi)星觀測(cè)的高覆蓋率對(duì)于中尺度區(qū)域短期天氣預(yù)報(bào)尤其重要。研究中為了避免區(qū)域天氣預(yù)報(bào)受到模式側(cè)邊界條件的影響,所以同化的循環(huán)時(shí)間設(shè)為6 h間隔,而且同化次數(shù)不超過(guò)3次。

    數(shù)值試驗(yàn)的設(shè)計(jì)主要是為了更好的評(píng)估NOAA-15衛(wèi)星的AMSU-A輻射資料同化對(duì)美國(guó)墨西哥灣24 h沿岸降水預(yù)報(bào)的附加改進(jìn)效果。研究設(shè)計(jì)了4個(gè)數(shù)值試驗(yàn),試驗(yàn)方案見表1。在第一個(gè)試驗(yàn)即控制試驗(yàn)中(CTRL),常規(guī)資料和NOAA-18 和MetOp-A的AMSU-A資料分別在協(xié)調(diào)世界時(shí)2008年5月22日的12:00、18:00和23日的00:00進(jìn)行同化。第二個(gè)試驗(yàn)(EXP1)是在控制試驗(yàn)的基礎(chǔ)上增加了NOAA-15的AMSU-A資料。為了進(jìn)一步明確AMSU-A資料同化影響的主要來(lái)源,我們還設(shè)計(jì)了另外兩個(gè)試驗(yàn)分析AMSU-A窗區(qū)通道的影響。第三個(gè)試驗(yàn)(EXP2)與第二個(gè)試驗(yàn)同化同樣種類的觀測(cè)資料,但是AMSU-A的窗區(qū)通道1和2在陸地上資料被剔除,在EXP3試驗(yàn)中,AMSU-A窗區(qū)通道1和2的陸地和海洋資料同時(shí)被剔除。

    表1 試驗(yàn)設(shè)計(jì)方案Table 1 Experiment design

    3 AMSU-A資料簡(jiǎn)介

    圖1 1998~2014年極軌衛(wèi)星升軌(實(shí)線)和降軌(虛線)局地越赤道時(shí)間(紅色為NOAA-15,藍(lán)色為NOAA-18,綠色為NOAA-19,桔黃色為MetOp-A)Fig. 1 Local equator crossing time (LECT) of polar orbit satellites NOAA-15 (red), NOAA-18 (blue), NOAA-19 (green) and, MetOp-A (orange) at the ascending node (solid) and descending node (dashed) from 1998 to 2014

    圖2 AMSU-A儀器1至15通道在美國(guó)標(biāo)準(zhǔn)大氣條件下的權(quán)重函數(shù)垂直分布圖[灰色橫線為WRF-ARW(Advanced Research WRF)的模式層隨高度分布情況]Fig. 2 Weighting functions for AMSU-A channels 1–15 overlapped onto the WRF-ARW (Advanced Research WRF) model levels (gray horizontal line)

    AMSU-A儀器可以測(cè)量微波頻段在氧氣吸收帶的輻射量,所以該儀器主要用于提供除了強(qiáng)降水天氣條件之外的近乎全天氣條件的大氣溫度廓線信息。AMSU-A是一種交軌微波探測(cè)儀,它的探測(cè)角度是在以星下點(diǎn)為中心的±48.7°以內(nèi)。該儀器有15個(gè)探測(cè)通道,頻率范圍是在23.8~89.0 GHz之間。儀器的瞬時(shí)視場(chǎng)(Field-of-View,簡(jiǎn)稱FOV)為3.3°,在星下點(diǎn)兩側(cè)分別有15個(gè)相鄰的視場(chǎng),每條掃描線共30個(gè)FOV。大氣的溫度廓線信息主要來(lái)自位于50~60 GHz的通道3~14的觀測(cè)信息,隨著通道從3~14的逐漸增加,觀測(cè)高度逐漸增高,AMSU-A最高可以觀測(cè)到地球以上大約42 km高度(大約2 hPa)大氣溫度信息。通道1、2和15主要用于獲得地表信息和云的特征。AMSU-A資料星下點(diǎn)處的空間分辨率約為48 km,隨著掃描角的增加,其分辨率以二次方的速度逐漸加粗,在沿軌道方向,最大掃描角(第1個(gè)FOV和第30個(gè)FOV)的視場(chǎng)覆蓋范圍可以比星下點(diǎn)視場(chǎng)大2倍,在交軌方向則會(huì)大到3倍以上。

    圖2給出了AMSU-A的15個(gè)通道權(quán)重函數(shù)垂直分布圖,其上疊加的橫線為WRF-ARW模式的50個(gè)sigma層次分布特征,模式層頂為10 hPa。從圖2中可以看出,12個(gè)AMSU-A溫度探測(cè)通道均勻的分布在整個(gè)地球大氣層。WRF-ARW模式層次設(shè)置主要是為了保證溫度探測(cè)通道4-9的有效同化和用于質(zhì)量控制的通道1和2的合理應(yīng)用,從而保證那些對(duì)于天氣尺度和中尺度天氣系統(tǒng)最重要的溫度擾動(dòng)能夠得到有效的觀測(cè)。

    在GSI資料同化系統(tǒng)中,AMSU-A資料的質(zhì)量控制主要是依賴于AMSU-A通道1和2的觀測(cè)和模擬亮溫進(jìn)行的。其他通道的觀測(cè)和模擬亮溫偏差(O-B)以及通道1~3的O-B與地表發(fā)射率的比值也是質(zhì)量控制的因子。關(guān)于質(zhì)量控制更詳細(xì)的描述請(qǐng)見Zou et al.(2013)。

    4 數(shù)值試驗(yàn)

    利用三維變分同化系統(tǒng)進(jìn)行的同化試驗(yàn)一般是取6 h的同化時(shí)間窗,分別在協(xié)調(diào)世界時(shí)00:00、06:00、12:00和18:00時(shí)刻同化常規(guī)資料和衛(wèi)星資料。對(duì)于極軌衛(wèi)星而言,衛(wèi)星資料出現(xiàn)的當(dāng)?shù)貢r(shí)間是幾乎固定。所以對(duì)于中尺度的區(qū)域數(shù)值預(yù)報(bào)而言,是否有極軌衛(wèi)星資料主要取決于衛(wèi)星的LECT和模式區(qū)域的時(shí)區(qū)之間的對(duì)應(yīng)關(guān)系。圖3給出了AMSU-A資料的空間覆蓋率隨時(shí)間變化特征。圖中分別給出了2008年5月22日協(xié)調(diào)世界時(shí)00:00、06:00、12:00和18:00時(shí)刻±3h內(nèi)的AMSU-A資料空間覆蓋情況。其中NOAA-18為下午星(當(dāng)?shù)貢r(shí)間14:00),MetOp-A為上午星(當(dāng)?shù)貢r(shí)間09:30),NOAA-15則是早晨星(當(dāng)?shù)貢r(shí)間07:00)。從圖3可以看出,NOAA-18和MetOp-A在06:00(協(xié)調(diào)世界時(shí),下同)和18:00對(duì)美國(guó)有很好的覆蓋率,但是在00:00和12:00,這兩顆衛(wèi)星在美國(guó)的中部和墨西哥灣幾乎沒(méi)有觀測(cè)。早晨星NOAA-15的觀測(cè)資料則能夠很好的填補(bǔ)上午星和下午星的這段資料空白,這就可以保證那些影響墨西哥灣和東海岸降水的上游天氣系統(tǒng)能夠得到有效的觀測(cè)。

    圖4 試驗(yàn)EXP1中2008年5月22日12:00(協(xié)調(diào)世界時(shí))時(shí)刻N(yùn)OAA-15衛(wèi)星AMSU-A資料1到10以及15通道同化前(O-B,左列)和同化后(O-A,右列)觀測(cè)和模擬亮溫偏差的空間分布特征(黑色點(diǎn)線表示沒(méi)有通過(guò)質(zhì)量控制的觀測(cè)資料)Fig. 4 Differences between observations and background fields (O-B, left panels), and between observations and analysis fields (O-A, right panels) for NOAA-15 AMSU-A channels 1–10 and channel 15 data assimilated at 1200 UTC 22 May 2008 from expt EXP1. AMSU-A data points not assimilated are indicated by black dots

    圖4 (續(xù))Fig. 4 (Continued)

    圖4 (續(xù))Fig. 4 (Continued)

    為了評(píng)估AMSU-A資料同化影響,我們首先分析觀測(cè)和模擬亮溫偏差在同化前(O-B)和同化后(O-A)的改變情況。圖4給出了EXP1的NOAA-15衛(wèi)星的AMSU-A通道1~10在2008年5 月22日12:00時(shí)刻的O-B和O-A的空間分布特征。并不是所有的資料都能夠進(jìn)入同化系統(tǒng),圖4中的黑色點(diǎn)是被剔除的觀測(cè)資料。圖4左列為O-B,可以看出多個(gè)通道的O-B可以超出±0.5K的范圍,但是同化后,絕大部分的O-A(右列)都可以減小到±0.2 K的范圍。

    圖5 試驗(yàn)EXP1中2008年5月22日12:00時(shí)刻N(yùn)OAA-15衛(wèi)星AMSU-A資料1~10以及15通道亮溫O-B(虛線)和O-A(實(shí)線)平均值隨視場(chǎng)的變化曲線Fig. 5 Means varying with FOV (Field-of-View) for the O-B (dashed) and O-A (solid) brightness temperature differences of NOAA-15 AMSU-A channels 1–10 and 15 at 1200 UTC 22 May 2008 from expt EXP1

    O-B和O-A的均值和標(biāo)準(zhǔn)差分別表示在圖5和圖6,通道1~3和通道15的O-B均值最高可以達(dá)到±2 K的范圍,但是O-A則基本可以控制的±0.5 K以內(nèi),通道4~19的O-B均值則是表現(xiàn)出明顯的負(fù)偏差,通道9的負(fù)偏差接近-0.8 K,這些通道同化后的O-A均值均接近±0.2 K,而且沒(méi)有明顯的負(fù)偏差和正偏差特征,相應(yīng)的標(biāo)準(zhǔn)差經(jīng)過(guò)同化后也得到了顯著的降低。另外均值和標(biāo)準(zhǔn)差隨FOV的變化幅度也得到了明顯的減弱。

    圖7給出了試驗(yàn)CTRL和EXP1在200、500 和850 hPa的大氣溫度分析增量的空間分布情況。CRTL和EXP1在200 hPa的溫度分析增量表現(xiàn)出反向變化的現(xiàn)象,CRTL試驗(yàn)的溫度分析增量主要為正值,只是在模式的北部有小部分區(qū)域表現(xiàn)為負(fù)增量,但是EXP1的溫度分析增量則主要是負(fù)值,只是在中部陸地地區(qū)略有一小塊弱的正值,EXP1的負(fù)增量在海洋上空尤其明顯。對(duì)于500和850 hPa高度,平均而言,EXP1的分析增量在24°N以北區(qū)域要比CTRL冷,但是在其南部則是比CTRL要暖。圖8則是給出了兩個(gè)試驗(yàn)在500和850 hPa高度相對(duì)濕度增量的空間分布,可以看出在同化了NOAA-15的AMSU-A資料后,對(duì)流層低層的濕度增加,兩個(gè)試驗(yàn)的相對(duì)濕度增量最大可以達(dá)到20%。在對(duì)流層中層,EXP1的墨西哥灣中部大氣的濕度也是增加,這就有助于沿岸對(duì)流系統(tǒng)的水汽輸送。

    圖9給出了EXP1試驗(yàn)中對(duì)流系統(tǒng)發(fā)生、發(fā)展的大尺度環(huán)境場(chǎng),包括每3h的500 hPa高度場(chǎng),風(fēng)矢量場(chǎng)和風(fēng)速特征,時(shí)間是2008年5月23日00:00~21:00。在5月23日00:00時(shí)刻,在墨西哥灣92°W附近的500 hPa高度,有一個(gè)明顯的高空槽形成,在隨后的15 h內(nèi),槽逐漸東移,同時(shí)有加深并向南擴(kuò)展的現(xiàn)象,在09:00和12:00,槽發(fā)展到最鼎盛時(shí)期,同時(shí)也伴隨著最大的對(duì)流性降水現(xiàn)象。

    圖10給出了觀測(cè)和模擬的2008年5月23日06:00~09:00時(shí)段內(nèi)的墨西哥灣沿岸3 h累積降水空間分布結(jié)構(gòu)。圖10a為NCEP多普勒雷達(dá)資料反演得到的降水分布,圖10b和c則分別給出了CTRL 和EXP1的模擬降水空間特征。可以看出CTRL的降水落區(qū)出現(xiàn)了明顯的偏差,EXP1則能夠成功的模擬觀測(cè)降水的大值區(qū)。從降水的時(shí)間演變圖(圖略)還可以看出EXP1可以很好的模擬降水區(qū)逐漸東移的現(xiàn)象,但是EXP1的3 h降水量比觀測(cè)偏大。圖11展示了3 h累積降水傳統(tǒng)TS(Threat Score)評(píng)分隨時(shí)間變化情況。傳統(tǒng)TS評(píng)分的定義(Anthes et al., 1989 )為:TS=H/(F+O-H),這里的F為預(yù)報(bào)的觀測(cè)站降水超過(guò)臨界值的站數(shù),O是表示觀測(cè)降水超過(guò)臨界值的站數(shù),H是預(yù)報(bào)降水和觀測(cè)降水同時(shí)超過(guò)臨界值的站數(shù)。TS評(píng)分主要表現(xiàn)了模式對(duì)給定降水臨界值的預(yù)報(bào)準(zhǔn)確程度。為了進(jìn)一步明確模式無(wú)雨區(qū)對(duì)降水評(píng)分的影響,我們分別給出了整個(gè)模式區(qū)域和沿岸降水區(qū)[見圖10a中的框區(qū),(28°~32°N,83°~92°W)]的TS評(píng)分結(jié)果。從TS評(píng)分的直方圖中可以清楚看到EXP1的評(píng)分明顯高于CTRL。同化NOAA-15衛(wèi)星AMSU-A資料的最大正影響出現(xiàn)在06:00~18:00時(shí)間,這個(gè)時(shí)段的對(duì)流系統(tǒng)最為活躍??刂圃囼?yàn)CTRL對(duì)于臨界值大于5 mm的降水預(yù)報(bào)技巧明顯低于1 mm以下的降水,但是EXP1的小雨和大雨的預(yù)報(bào)技巧則是相當(dāng)?shù)摹?/p>

    圖6 同圖5,但為標(biāo)準(zhǔn)差隨視場(chǎng)的變化曲線Fig. 6 As in Fig. 5, but for the standard deviation

    圖7 試驗(yàn)CTRL(左列)和EXP1(右列)在2008年5月22日12:00 時(shí)刻(a、b)200 hPa、(c、d)500 hPa和(e、f)850 hPa高度的溫度分析增量(單位:°C)空間分布Fig. 7 Analysis increments of temperature (color scale; units: °C) at (a, b) 200 hPa, (c, d) 500 hPa and (e, f) 850 hPa at 1200 UTC 22 May 2008 from expts CTRL1 (left column) and EXP1 (right column)

    圖8 2008年5月22日12:00 時(shí)刻試驗(yàn)EXP1和CTRL的相對(duì)濕度分析增量差值(EXP1減去CTRL)的空間分布圖:(a)500 hPa;(b)850 hPaFig. 8 Analysis difference of the atmospheric relative humidity between EXP1 and CTRL1 (EXP1 minus CTRL1) at (a) 500 hPa and (b) 850 hPa at 1200 UTC 22 May 2008

    圖10 2008年5月23日06:00~09:00的3 h累積降水(單位:mm)空間分布:(a)NCEP觀測(cè);(b)試驗(yàn)CTRL模擬;(c)試驗(yàn)EXP1模擬Fig. 10 Accumulative rainfall (units: mm) during 0600–0900 UTC 23 May 2008 from (a) NCEP observations, model forecasts (b) without and (c) with NOAA-15 data

    圖11 試驗(yàn)CTRL和EXP1在2008年5月23日模擬的3h累積降水的TS(Threat Score)評(píng)分隨時(shí)間變化直方圖。TS評(píng)分計(jì)算中設(shè)定的臨界值分別為1、5、10和15 mm。圖中CRTLs和EXP1s代表TS評(píng)分的檢驗(yàn)區(qū)域?yàn)檠匕秴^(qū)域,即圖10a中的長(zhǎng)方形框區(qū)Fig. 11 Threat scores of 3-h accumulative rainfall during 0000 UTC–2400 UTC 23 May 2008 from CTRL1 and EXP1 in the model domain (gray and cyan bars) and the small coastal area indicated in Fig. 10a (red and blue bars) at thresholds 1, 5, 10 and 15 mm

    圖12 試驗(yàn)EXP3(左列)和EXP1(右列)在2008年5月22日12:00時(shí)刻925 hPa高度的溫度分析場(chǎng)和對(duì)應(yīng)的分析增量:(a、b)溫度分析場(chǎng);(c、d)溫度分析增量。單位:KFig. 12 (a, b) Temperature analysis fields (units: K) and (c, d) analysis increments (units: K) at 925 hPa at 1200 UTC 22 May 2008 in EXP3 (left panels) and EXP1 (right panels)

    敏感性試驗(yàn)研究表明AMSU-A的兩個(gè)窗區(qū)通道資料的同化對(duì)沿岸定量降水預(yù)報(bào)改進(jìn)有重要的作用。圖12展示了EXP1和EXP3在925 hPa的大氣溫度分析場(chǎng)和分析增量,比較結(jié)果表明當(dāng)兩個(gè)AMSU-A窗區(qū)通道在海洋和陸地上的資料被剔除后,墨西哥灣東部的溫度增量明顯減小,同化窗區(qū)通道資料可以使得對(duì)流發(fā)展地區(qū)的溫度和水汽增加,從而更有助于對(duì)流活動(dòng)的發(fā)生和發(fā)展。EXP3 的3 h累積降水的TS評(píng)分也明顯低于EXP1(圖13)。相比而言,EXP2的TS評(píng)分結(jié)果則是和EXP1的結(jié)果比較相似,只是在15 mm以上的大雨區(qū),EXP2的TS評(píng)分結(jié)果比EXP1要低從TS評(píng)分結(jié)果上來(lái)看,大部分窗區(qū)通道的正影響是來(lái)自于海洋地區(qū)的資料,陸地地區(qū)的窗區(qū)資料主要對(duì)15 mm以上的降水預(yù)報(bào)有一定的正影響。

    5 結(jié)論與討論

    本文簡(jiǎn)單地評(píng)估了早晨星軌道的極軌衛(wèi)星AMSU-A資料同化對(duì)墨西哥灣沿岸降水預(yù)報(bào)的影響。研究表明在GSI資料同化系統(tǒng)中同化了早晨星NOAA-15的AMSU-A輻射資料,可以明顯改善WRF-ARW模式對(duì)墨西哥灣24 h沿岸降水預(yù)報(bào)技巧,也進(jìn)一步證明了極軌衛(wèi)星的三軌衛(wèi)星觀測(cè)系統(tǒng)對(duì)于區(qū)域降水預(yù)報(bào)改進(jìn)的重要作用。現(xiàn)有的極軌衛(wèi)星AMSU-A:資料包括NOAA-15、NOAA-18、 NOAA-19、MetOp-A和MetOp-B,都已經(jīng)成為多個(gè)國(guó)家業(yè)務(wù)預(yù)報(bào)中心常規(guī)同化的衛(wèi)星資料,這些衛(wèi)星的資料也隨時(shí)可以從美國(guó)的JPSS(Joint Polar Satellite System)資料中心下載,所以極軌衛(wèi)星AMSU-A資料同化對(duì)今后數(shù)值預(yù)報(bào)業(yè)務(wù)和研究的影響也會(huì)與日俱增,這也使得如何配置極軌衛(wèi)星系統(tǒng)顯得尤為重要。本研究初步證明了極軌衛(wèi)星三軌觀測(cè)系統(tǒng)維持的重要性,而早晨星的超期服役問(wèn)題正是維持三軌觀測(cè)系統(tǒng)所面臨的主要問(wèn)題,不間斷的發(fā)射這一類搭載了AMSU-A的早晨星對(duì)于維持全球極軌衛(wèi)星三軌觀測(cè)系統(tǒng)是值得考慮的方法。由于目前NOAA-15是唯一的一顆正在運(yùn)行的、已遠(yuǎn)超過(guò)其正常運(yùn)行期的早晨星,通過(guò)技術(shù)手段維持NOAA-15的AMSU-A超長(zhǎng)期運(yùn)行也是解決方法之一。

    圖13 同圖11,但 為試驗(yàn)EXP1、EXP2和EXP3的TS評(píng)分直方圖Fig. 13 Same as Fig. 11, but for EXP1, EXP2 and EXP3

    本文的研究結(jié)果是對(duì)早期基于ECMWF的數(shù)值天氣預(yù)報(bào)系統(tǒng)研究結(jié)果(Kelly, 2002;13th ITSC)的很好補(bǔ)充,早期研究結(jié)果表明NOAA-15、NOAA-16和NOAA-17這三個(gè)衛(wèi)星可以提供全球的三軌衛(wèi)星AMSU-A觀測(cè)資料,從500 hPa高度場(chǎng)異常相關(guān)系數(shù)檢驗(yàn)結(jié)果表明同化三軌AMSU-A資料可以獲得更好的中期數(shù)值預(yù)報(bào)水平。但是對(duì)于區(qū)域數(shù)值天氣預(yù)報(bào)而言,極軌衛(wèi)星AMSU-A資料同化對(duì)區(qū)域降水預(yù)報(bào)的影響有明顯的區(qū)域依賴性(主要是經(jīng)度依賴性)。在今后的研究中,我們還將進(jìn)一步選取更多個(gè)例和進(jìn)行更長(zhǎng)時(shí)間尺度的同化試驗(yàn)分析,另外一個(gè)較好的試驗(yàn)時(shí)段是從2002年6月24日NOAA-17發(fā)射到2003年10月30日NOAA-16的停止服役,在這一年左右的時(shí)間,NOAA-15、NOAA-16和NOAA-17同時(shí)運(yùn)行,這三個(gè)衛(wèi)星也可以提供更高質(zhì)量的三軌AMSU-A觀測(cè)資料,今后研究中我們將利用新的觀測(cè)資料進(jìn)行更多的評(píng)估和試驗(yàn)。

    參考文獻(xiàn)(References)

    Anthes R A, Kuo Y, Hsie E, et al. 1989. Estimation of skill and uncertainty in regional numerical models [J]. Quart. J. Roy. Meteor. Soc., 115, 763–806, doi: 10.1002/qj.49711548803.

    Eyre J R, Kelly G, McNally A P, et al. 1993. Assimilation of TOVS radiance information through one-dimensional variational analysis [J]. Quart. J. Roy. Meteor. Soc, 119: 1427–1463, doi:10.1002/qj.49711951411.

    Han Y, Weng F Z, Liu Q H, et al. 2007. A fast radiative transfer model for SSMIS upper atmosphere sounding channels [J]. J. Geophys. Res., 112: D11121, doi:10.1029/2006JD008208.

    Hong S Y, Dudhia J. 2003. Testing of a new non-local boundary layer vertical diffusion scheme in numerical weather prediction applications [C]. 20th Conference on Weather Analysis and Forecasting/16th Conference on Numerical Weather Prediction, Seattle, WA.

    Hong S Y, Lim J -O J. 2006. The WRF single-moment 6-class microphysics scheme (WSM6) [J]. J. Korean Meteor. Soc., 42, 129–151.

    Kain J S. 2004. The Kain-Fritsch convective parameterization: An update [J]. J. Appl. Meteor., 43, 170–181, doi: 10.1175/1520-0450(2004)043<0170: TKCPAU>2.0.CO;2

    Kain J S, Fritsch J M. 1990. A one-dimensional entraining/detraining plume model and its application in convective parameterization [J]. J. Atmos. Sci., 47, 2784–2802, doi: 10.1175/1520-0469(1990)047<2784: AODEPM>2.0.CO;2

    Kain, J. S., and J. M. Fritsch, 1993: Convective parameterization for mesoscale models: The Kain-Fritsch scheme. The Representation of Cumulus Convection in Numerical Models, Meteor. Monogr., 24, Amer. Meteor. Soc., 165–170.

    Kelly G, Anderson E, Bauer P, et al. 2003. Use of satellite radiances in the ECMWF system [C]// Proceedings of the 13thInternational TOVS Conference. Quebec, Canada.

    McMillin L, Fleming H. 1976. Atmospheric transmittance of an absorbing gas: A computationally fast and accurate transmittance model for absorbing gases with constant mixing ratios in inhomogeneous atmospheres [J]. Appl. Optics, 15: 358–363, doi:10.1364/AO.15.000358.

    Purser R J, Wu Wanshu, Parrish D F, et al. 2003a. Numerical aspects of the application of recursive filters to variational statistical analysis. Part I: Spatially homogeneous and isotropic Gaussian covariances [J]. Mon. Wea. Rev., 131: 1524–1535, doi:10.1175//1520-0493(2003)131<1524: NAOTAO>2.0.CO;2.

    Purser R J, Wu Wanshu, Parrish D F, et al. 2003b. Numerical aspects of the application of recursive filters to variational statistical analysis. Part II: Spatially inhomogeneous and anisotropic general covariances [J]. Mon. Wea. Rev., 131: 1536–1548, doi:10.1175//2543.1.

    Saunders R, Matricardi M, Brunel P. 1999. An improved fast radiative transfer model for assimilation of satellite radiance observations [J]. Quart. J. Roy. Meteor. Soc., 125: 1407–1425, doi:10.1002/qj.1999. 49712555615.

    Saunders R, Rayer P, Brunel P, et al. 2007. A comparison of radiative transfer models for simulating Atmospheric Infrared Sounder (AIRS) radiance [J]. J. Geophys. Res., 112: D01S90, doi:10.1029/ 2006JD007088.

    Weng F Z. 2007. Advances in radiative transfer modeling in support of satellite data assimilation [J]. J. Atmos. Sci., 64: 3799–3807, doi:10.1175/2007JAS2112.1.

    Wu Wanshu, Purser R J, Parrish D F. 2002. Three-dimensional variational analysis with spatially inhomogeneous covariances [J]. Mon. Wea. Rev., 130: 2905–2916, doi:10.1175/1520-0493(2002)130<2905:TDVAWS> 2.0.CO;2.

    Zou X, Weng F, Zhang B, et al. 2013. Impacts of assimilation of ATMS data in HWRF on track and intensity forecasts of 2012 four landfall hurricanes [J]. J. Geophys. Res. Atmos., 118: 11558–11576, doi:10.1002/2013JD020405.

    鄒曉蕾,秦正坤,翁富忠. 2016. 晨昏軌道微波溫度計(jì)資料同化對(duì)降水定量預(yù)報(bào)的影響及其對(duì)三軌衛(wèi)星系統(tǒng)的意義 [J]. 大氣科學(xué), 40 (1): 46?62. Zou Xiaolei, Qin Zhengkun, Weng Fuzhong. 2016. Impact of dawn–dusk satellite AMSU-A data on quantitative precipitation forecasts and the implications for three-orbit constellation [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 40 (1): 46?62, doi: 10.3878/j.issn.1006-9895.1508.15137.

    作者簡(jiǎn)介鄒曉蕾,女,1960年出生,博士,教授,主要從事資料同化研究。E-mail: xiaoleizou@nuist.edu.cn

    收稿日期2015-03-01;網(wǎng)絡(luò)預(yù)出版日期 2015-08-24

    doi:10.3878/j.issn.1006-9895.1508.15137

    文章編號(hào)1006-9895(2016)01-0046-17

    中圖分類號(hào)P456.7

    文獻(xiàn)標(biāo)識(shí)碼A

    国产又色又爽无遮挡免| 夫妻性生交免费视频一级片| 免费少妇av软件| 亚洲av成人不卡在线观看播放网 | 真人做人爱边吃奶动态| 一级,二级,三级黄色视频| 国产成人精品在线电影| 中文字幕亚洲精品专区| 亚洲自偷自拍图片 自拍| 50天的宝宝边吃奶边哭怎么回事| 午夜激情久久久久久久| 亚洲色图 男人天堂 中文字幕| 免费高清在线观看日韩| 19禁男女啪啪无遮挡网站| 日韩视频在线欧美| 亚洲五月色婷婷综合| 国产1区2区3区精品| 国产一区二区 视频在线| 久久久久久亚洲精品国产蜜桃av| 欧美亚洲 丝袜 人妻 在线| 夫妻性生交免费视频一级片| 一二三四社区在线视频社区8| 国产在视频线精品| 中文精品一卡2卡3卡4更新| 欧美日韩福利视频一区二区| 麻豆国产av国片精品| 天天添夜夜摸| 久久毛片免费看一区二区三区| 老熟女久久久| 狠狠精品人妻久久久久久综合| 18禁观看日本| 人妻 亚洲 视频| 美女午夜性视频免费| 精品人妻一区二区三区麻豆| 中文乱码字字幕精品一区二区三区| 国产免费视频播放在线视频| 国产精品一区二区免费欧美 | 日本猛色少妇xxxxx猛交久久| 成人国产av品久久久| 欧美激情高清一区二区三区| 国产精品香港三级国产av潘金莲 | 欧美日韩成人在线一区二区| 久久性视频一级片| 男女边摸边吃奶| 十八禁人妻一区二区| 亚洲一区二区三区欧美精品| 亚洲七黄色美女视频| 国产精品.久久久| 91麻豆精品激情在线观看国产 | 欧美国产精品一级二级三级| 日韩中文字幕视频在线看片| 婷婷丁香在线五月| 精品久久蜜臀av无| 日本vs欧美在线观看视频| 免费不卡黄色视频| 狂野欧美激情性xxxx| 视频区图区小说| 日韩一卡2卡3卡4卡2021年| 久久这里只有精品19| 国产成人一区二区三区免费视频网站 | 久久毛片免费看一区二区三区| 国产欧美日韩一区二区三 | 丁香六月欧美| 一级毛片我不卡| a 毛片基地| 一级a爱视频在线免费观看| 美女视频免费永久观看网站| 成年av动漫网址| 久久久久久久国产电影| 亚洲国产精品国产精品| 97在线人人人人妻| 在线观看免费视频网站a站| 99国产精品免费福利视频| 一本一本久久a久久精品综合妖精| 精品少妇黑人巨大在线播放| 男女之事视频高清在线观看 | av又黄又爽大尺度在线免费看| 女人精品久久久久毛片| 91老司机精品| 日韩免费高清中文字幕av| 久热爱精品视频在线9| 国产成人欧美| av天堂在线播放| 赤兔流量卡办理| 国产不卡av网站在线观看| 欧美日本中文国产一区发布| 女人久久www免费人成看片| 黄频高清免费视频| 只有这里有精品99| 亚洲av男天堂| 亚洲欧洲日产国产| 黄片播放在线免费| 天堂俺去俺来也www色官网| 亚洲国产毛片av蜜桃av| 久热爱精品视频在线9| 女人被躁到高潮嗷嗷叫费观| 国产精品秋霞免费鲁丝片| 中文欧美无线码| 观看av在线不卡| 色综合欧美亚洲国产小说| 天堂8中文在线网| 精品少妇内射三级| 欧美性长视频在线观看| 国产精品一二三区在线看| 97在线人人人人妻| 国产成人一区二区在线| 精品国产乱码久久久久久男人| 最新的欧美精品一区二区| 久久精品久久久久久噜噜老黄| 日本色播在线视频| 国产男人的电影天堂91| 国产成人av激情在线播放| 一区二区三区乱码不卡18| 少妇 在线观看| 日本五十路高清| 高清视频免费观看一区二区| 久久精品人人爽人人爽视色| 永久免费av网站大全| 免费观看a级毛片全部| 成年女人毛片免费观看观看9 | 十分钟在线观看高清视频www| 国产精品欧美亚洲77777| 久久综合国产亚洲精品| 午夜91福利影院| 欧美成人午夜精品| 久久精品国产a三级三级三级| 男女无遮挡免费网站观看| 女性生殖器流出的白浆| 国产女主播在线喷水免费视频网站| av福利片在线| 亚洲综合色网址| 9色porny在线观看| 一本大道久久a久久精品| 另类精品久久| 国产一区有黄有色的免费视频| 赤兔流量卡办理| 精品国产乱码久久久久久男人| 免费高清在线观看视频在线观看| 18禁裸乳无遮挡动漫免费视频| 成人影院久久| 91九色精品人成在线观看| 韩国高清视频一区二区三区| 乱人伦中国视频| 久久久欧美国产精品| 欧美xxⅹ黑人| 美女主播在线视频| 一级片'在线观看视频| 看免费av毛片| 97人妻天天添夜夜摸| 伊人久久大香线蕉亚洲五| 日韩 欧美 亚洲 中文字幕| 90打野战视频偷拍视频| 国产成人影院久久av| 国产97色在线日韩免费| 丰满迷人的少妇在线观看| 亚洲国产精品成人久久小说| 日韩制服丝袜自拍偷拍| 成年动漫av网址| 夫妻性生交免费视频一级片| 超色免费av| av网站在线播放免费| av国产久精品久网站免费入址| 少妇 在线观看| 巨乳人妻的诱惑在线观看| 欧美日韩综合久久久久久| 欧美日韩视频精品一区| 熟女少妇亚洲综合色aaa.| 国产免费视频播放在线视频| 高清av免费在线| 男女边吃奶边做爰视频| 亚洲欧美一区二区三区久久| 秋霞在线观看毛片| 久久久国产欧美日韩av| 一二三四社区在线视频社区8| 电影成人av| 成年女人毛片免费观看观看9 | 久久精品国产综合久久久| 亚洲精品成人av观看孕妇| 狠狠精品人妻久久久久久综合| 深夜精品福利| 国产成人精品久久二区二区91| 18禁黄网站禁片午夜丰满| 中国国产av一级| www日本在线高清视频| 国产欧美日韩一区二区三区在线| 成人午夜精彩视频在线观看| 老司机影院毛片| 69精品国产乱码久久久| 美女扒开内裤让男人捅视频| netflix在线观看网站| 婷婷色av中文字幕| 久久精品aⅴ一区二区三区四区| 国产精品av久久久久免费| 亚洲美女黄色视频免费看| 国产精品 欧美亚洲| av国产久精品久网站免费入址| 久久精品久久久久久久性| 国产精品久久久av美女十八| 国产欧美亚洲国产| 亚洲av日韩在线播放| 麻豆国产av国片精品| e午夜精品久久久久久久| 国产成人av激情在线播放| 久久国产精品影院| 欧美在线黄色| 最近中文字幕2019免费版| 欧美成狂野欧美在线观看| av一本久久久久| 国语对白做爰xxxⅹ性视频网站| 丝袜在线中文字幕| 国产淫语在线视频| 亚洲av电影在线进入| 制服诱惑二区| 久久国产精品人妻蜜桃| 亚洲精品一卡2卡三卡4卡5卡 | 欧美人与性动交α欧美软件| 777久久人妻少妇嫩草av网站| 亚洲欧美色中文字幕在线| 精品免费久久久久久久清纯 | 自线自在国产av| 国产精品久久久人人做人人爽| 男女下面插进去视频免费观看| 中文字幕色久视频| 9191精品国产免费久久| 黑人巨大精品欧美一区二区蜜桃| 精品少妇黑人巨大在线播放| 久久性视频一级片| 国产精品一区二区在线不卡| 亚洲精品一卡2卡三卡4卡5卡 | 色婷婷久久久亚洲欧美| 无遮挡黄片免费观看| 波野结衣二区三区在线| 观看av在线不卡| 丰满迷人的少妇在线观看| 九色亚洲精品在线播放| 国产亚洲精品第一综合不卡| 国产一区二区激情短视频 | 18在线观看网站| 国产成人欧美在线观看 | 又紧又爽又黄一区二区| 精品一区二区三卡| 久久人妻熟女aⅴ| 极品人妻少妇av视频| 亚洲 欧美一区二区三区| 亚洲欧美日韩另类电影网站| 考比视频在线观看| 纵有疾风起免费观看全集完整版| 日韩制服骚丝袜av| 中文乱码字字幕精品一区二区三区| 欧美在线一区亚洲| 国产精品 国内视频| 精品久久久久久电影网| 久久影院123| 欧美久久黑人一区二区| 国产老妇伦熟女老妇高清| 夜夜骑夜夜射夜夜干| 日韩制服丝袜自拍偷拍| 国产在线免费精品| 日本av免费视频播放| 中国国产av一级| 亚洲三区欧美一区| 久久亚洲国产成人精品v| 伊人亚洲综合成人网| 国产精品人妻久久久影院| 免费在线观看视频国产中文字幕亚洲 | 亚洲国产av新网站| 人妻人人澡人人爽人人| 国产日韩欧美亚洲二区| 亚洲人成77777在线视频| 在线观看免费午夜福利视频| 校园人妻丝袜中文字幕| 每晚都被弄得嗷嗷叫到高潮| 99国产精品免费福利视频| 青草久久国产| 国产日韩一区二区三区精品不卡| 一二三四在线观看免费中文在| 日本欧美视频一区| 亚洲一码二码三码区别大吗| 久久久久国产一级毛片高清牌| 欧美 亚洲 国产 日韩一| 欧美精品高潮呻吟av久久| 最近中文字幕2019免费版| av一本久久久久| 婷婷色av中文字幕| 精品国产乱码久久久久久男人| svipshipincom国产片| 老汉色∧v一级毛片| 国产精品九九99| 香蕉国产在线看| 精品第一国产精品| 亚洲精品成人av观看孕妇| 午夜福利一区二区在线看| 亚洲精品自拍成人| 爱豆传媒免费全集在线观看| 老司机深夜福利视频在线观看 | 好男人电影高清在线观看| 捣出白浆h1v1| 中文欧美无线码| 校园人妻丝袜中文字幕| 晚上一个人看的免费电影| 人人妻人人澡人人爽人人夜夜| 日本vs欧美在线观看视频| 99国产精品一区二区三区| 波多野结衣av一区二区av| 性色av乱码一区二区三区2| 精品久久蜜臀av无| 国产精品久久久久久精品电影小说| 日韩欧美一区视频在线观看| 麻豆av在线久日| 日本av手机在线免费观看| 国产无遮挡羞羞视频在线观看| 黄片小视频在线播放| 91国产中文字幕| 高潮久久久久久久久久久不卡| 欧美+亚洲+日韩+国产| 国产精品久久久久久人妻精品电影 | 一级黄色大片毛片| 国产精品国产av在线观看| 国产精品久久久久久人妻精品电影 | 亚洲欧美精品综合一区二区三区| 国产精品国产av在线观看| 99香蕉大伊视频| 丰满少妇做爰视频| 欧美av亚洲av综合av国产av| 日韩av不卡免费在线播放| 伦理电影免费视频| 婷婷成人精品国产| 久久中文字幕一级| 自拍欧美九色日韩亚洲蝌蚪91| 国产一级毛片在线| 国产av精品麻豆| 国产激情久久老熟女| 免费看av在线观看网站| 亚洲中文av在线| e午夜精品久久久久久久| 久久久久网色| 深夜精品福利| 一级毛片 在线播放| 国产片特级美女逼逼视频| 1024香蕉在线观看| 最黄视频免费看| 真人做人爱边吃奶动态| 在线观看一区二区三区激情| 麻豆av在线久日| 国产免费视频播放在线视频| 精品人妻1区二区| 人成视频在线观看免费观看| 精品第一国产精品| 好男人电影高清在线观看| 免费观看人在逋| 亚洲欧洲精品一区二区精品久久久| 婷婷成人精品国产| 久久国产精品大桥未久av| 夫妻午夜视频| 中文字幕亚洲精品专区| 精品国产乱码久久久久久男人| 欧美另类一区| 老汉色∧v一级毛片| 国产亚洲精品久久久久5区| 夫妻午夜视频| 美女视频免费永久观看网站| 亚洲av日韩精品久久久久久密 | 国产成人一区二区三区免费视频网站 | 国产1区2区3区精品| 国产xxxxx性猛交| 婷婷色麻豆天堂久久| 少妇猛男粗大的猛烈进出视频| 国产精品麻豆人妻色哟哟久久| 亚洲精品自拍成人| 99久久综合免费| 国产真人三级小视频在线观看| 两性夫妻黄色片| 精品久久久精品久久久| 日韩精品免费视频一区二区三区| 久久天躁狠狠躁夜夜2o2o | 欧美老熟妇乱子伦牲交| 免费在线观看黄色视频的| a级毛片黄视频| 精品少妇黑人巨大在线播放| 亚洲免费av在线视频| 麻豆乱淫一区二区| 成年人免费黄色播放视频| 国产男人的电影天堂91| 一级黄色大片毛片| 女性生殖器流出的白浆| 啦啦啦在线观看免费高清www| 一边摸一边做爽爽视频免费| 90打野战视频偷拍视频| 精品久久蜜臀av无| 精品国产超薄肉色丝袜足j| 一级毛片黄色毛片免费观看视频| 欧美中文综合在线视频| 一级黄片播放器| 成人亚洲欧美一区二区av| 2018国产大陆天天弄谢| 欧美日韩黄片免| 90打野战视频偷拍视频| 欧美国产精品一级二级三级| 晚上一个人看的免费电影| 天堂中文最新版在线下载| 国产高清videossex| 韩国精品一区二区三区| 日韩一本色道免费dvd| 久久 成人 亚洲| 日韩电影二区| 美女主播在线视频| 美女国产高潮福利片在线看| 熟女av电影| 日韩制服骚丝袜av| 国产精品一区二区精品视频观看| 国产精品.久久久| 男女边吃奶边做爰视频| 亚洲伊人色综图| 国产色视频综合| 亚洲国产精品成人久久小说| 美女高潮到喷水免费观看| 国产无遮挡羞羞视频在线观看| 亚洲人成77777在线视频| 久久久久久久久免费视频了| 无遮挡黄片免费观看| 免费观看a级毛片全部| 一级a爱视频在线免费观看| 最近中文字幕2019免费版| 黄片播放在线免费| 在线观看一区二区三区激情| 国产精品偷伦视频观看了| 最近手机中文字幕大全| 亚洲国产精品999| av又黄又爽大尺度在线免费看| 免费在线观看影片大全网站 | 国产成人啪精品午夜网站| 乱人伦中国视频| 午夜av观看不卡| 手机成人av网站| 午夜福利在线免费观看网站| 国产日韩一区二区三区精品不卡| 欧美97在线视频| 亚洲av男天堂| 热re99久久精品国产66热6| 最近最新中文字幕大全免费视频 | 久久鲁丝午夜福利片| 国产精品欧美亚洲77777| av线在线观看网站| 亚洲人成电影观看| 老司机亚洲免费影院| 国产在线免费精品| 久久久精品94久久精品| 日韩制服丝袜自拍偷拍| 国产黄色免费在线视频| 国产精品熟女久久久久浪| 国产一区二区激情短视频 | 黑人欧美特级aaaaaa片| 在线观看免费视频网站a站| 亚洲国产av影院在线观看| 国产视频一区二区在线看| 亚洲欧美中文字幕日韩二区| 伊人久久大香线蕉亚洲五| 久久精品亚洲av国产电影网| 国产精品久久久久久精品电影小说| 精品国产乱码久久久久久男人| 91精品国产国语对白视频| 免费在线观看完整版高清| 日韩av免费高清视频| av在线老鸭窝| 99精品久久久久人妻精品| 国产视频一区二区在线看| 亚洲精品自拍成人| 国产老妇伦熟女老妇高清| 欧美日韩综合久久久久久| 黄色毛片三级朝国网站| www.999成人在线观看| 熟女少妇亚洲综合色aaa.| 新久久久久国产一级毛片| 99香蕉大伊视频| 亚洲精品在线美女| 麻豆av在线久日| 亚洲欧美清纯卡通| 国产精品一区二区在线观看99| 久久人人爽av亚洲精品天堂| 精品视频人人做人人爽| 亚洲国产av影院在线观看| 欧美日韩亚洲高清精品| 18在线观看网站| 日韩 亚洲 欧美在线| 男女免费视频国产| 99九九在线精品视频| 一二三四社区在线视频社区8| 激情五月婷婷亚洲| 午夜福利乱码中文字幕| 日韩中文字幕视频在线看片| 精品人妻一区二区三区麻豆| 亚洲欧洲国产日韩| 欧美 日韩 精品 国产| 大香蕉久久成人网| 国产黄色视频一区二区在线观看| 男女之事视频高清在线观看 | av天堂在线播放| 国产欧美日韩综合在线一区二区| 亚洲精品在线美女| 丝袜人妻中文字幕| av不卡在线播放| 视频区欧美日本亚洲| 一级毛片女人18水好多 | 亚洲精品国产色婷婷电影| 国产xxxxx性猛交| 七月丁香在线播放| 999久久久国产精品视频| 欧美人与性动交α欧美精品济南到| 一级黄片播放器| 伊人亚洲综合成人网| 在现免费观看毛片| 国产男女超爽视频在线观看| 99国产精品99久久久久| 一二三四社区在线视频社区8| 9色porny在线观看| 免费看不卡的av| 欧美中文综合在线视频| 国产成人系列免费观看| 99久久精品国产亚洲精品| 999久久久国产精品视频| 国产精品成人在线| 97人妻天天添夜夜摸| 永久免费av网站大全| 99国产精品一区二区三区| 一级毛片黄色毛片免费观看视频| 咕卡用的链子| 久久久国产精品麻豆| 中文字幕av电影在线播放| 久久鲁丝午夜福利片| 黄色怎么调成土黄色| 尾随美女入室| 国产一区二区 视频在线| a级片在线免费高清观看视频| 一级毛片电影观看| 午夜福利乱码中文字幕| xxx大片免费视频| 国产成人系列免费观看| 国产精品人妻久久久影院| 无遮挡黄片免费观看| 一本色道久久久久久精品综合| 免费黄频网站在线观看国产| 两人在一起打扑克的视频| 人人妻,人人澡人人爽秒播 | 欧美久久黑人一区二区| 午夜福利,免费看| 天堂8中文在线网| 午夜福利,免费看| 黑人欧美特级aaaaaa片| 亚洲第一青青草原| 久久影院123| 黄片小视频在线播放| 两个人看的免费小视频| 亚洲,欧美精品.| 国产日韩一区二区三区精品不卡| 黄片小视频在线播放| 欧美精品av麻豆av| 国产欧美日韩一区二区三 | 中文欧美无线码| 免费看av在线观看网站| 欧美日韩视频精品一区| 亚洲av综合色区一区| 麻豆国产av国片精品| 男人舔女人的私密视频| 18禁黄网站禁片午夜丰满| 嫩草影视91久久| av片东京热男人的天堂| 一个人免费看片子| 在线精品无人区一区二区三| av国产精品久久久久影院| 久久女婷五月综合色啪小说| 国产片特级美女逼逼视频| 嫁个100分男人电影在线观看 | 女警被强在线播放| 国产日韩欧美亚洲二区| 国产亚洲欧美精品永久| 国产xxxxx性猛交| 夫妻午夜视频| 午夜福利,免费看| 一级毛片女人18水好多 | 中文字幕av电影在线播放| 成人亚洲欧美一区二区av| 一本综合久久免费| 最近最新中文字幕大全免费视频 | 男女免费视频国产| 亚洲欧美精品综合一区二区三区| 欧美精品一区二区大全| 一区二区三区激情视频| 婷婷色综合大香蕉| 欧美+亚洲+日韩+国产| 大码成人一级视频| 欧美老熟妇乱子伦牲交| 亚洲精品av麻豆狂野| 国产伦人伦偷精品视频| 亚洲五月色婷婷综合| 欧美人与性动交α欧美精品济南到| 日韩制服丝袜自拍偷拍| 女人精品久久久久毛片| 欧美精品高潮呻吟av久久| 国产精品三级大全| 免费在线观看完整版高清| 亚洲,一卡二卡三卡| 国产男女超爽视频在线观看| 99国产精品一区二区三区| 电影成人av| 国产亚洲一区二区精品| 18禁黄网站禁片午夜丰满| videos熟女内射| 午夜久久久在线观看| 久久久欧美国产精品| 亚洲色图综合在线观看| 99国产综合亚洲精品| 国产老妇伦熟女老妇高清| 中国美女看黄片| 每晚都被弄得嗷嗷叫到高潮| 女人精品久久久久毛片|