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

    中國歷史探空資料部分等壓面位勢高度訂正

    2015-07-05 15:14:00熊安元胡開喜梁秀慧王志文楊燕茹鄒鳳玲
    應(yīng)用氣象學(xué)報 2015年3期
    關(guān)鍵詞:位勢探空氣層

    阮 新 熊安元 胡開喜 梁秀慧 王志文楊燕茹 鄒鳳玲

    1)(國家氣象信息中心,北京 100081)2)(河北省氣象局,石家莊 050021) 3)(北京市氣象局,北京 100089)

    ?

    中國歷史探空資料部分等壓面位勢高度訂正

    阮 新1)*熊安元1)胡開喜1)梁秀慧2)王志文3)楊燕茹1)鄒鳳玲1)

    1)(國家氣象信息中心,北京 100081)2)(河北省氣象局,石家莊 050021)3)(北京市氣象局,北京 100089)

    對中國高空探測歷史資料進(jìn)行靜力學(xué)檢查發(fā)現(xiàn),1963年前900,800,80 hPa和60 hPa位勢高度記錄錯誤率異常偏高,20世紀(jì)70年代中期個別站250 hPa和70 hPa位勢高度記錄錯誤率也異常偏高。通過對中國高空觀測規(guī)范歷史沿革的分析和資料驗證,指出受當(dāng)時計算條件所限,這些等壓面上的位勢高度是在其他規(guī)定等壓面的位勢高度查表計算完成后,利用時間高度曲線通過人工插值估算獲得,主觀因素影響較大,易產(chǎn)生較大誤差,因而錯誤率異常偏高。該文提出依據(jù)靜力學(xué)原理重新計算這些等壓面位勢高度記錄的訂正方案,并進(jìn)行效果檢驗。結(jié)果表明:該方案能夠比較精確地還原臺站觀測時對位勢高度的正確計算,錯誤資料訂正后,相關(guān)高度靜力學(xué)余差序列趨向均一,訂正方案合理、有效。

    探空資料; 位勢高度; 訂正

    引 言

    探空資料可以反映大氣的垂直結(jié)構(gòu),在大氣污染擴(kuò)散研究、數(shù)值天氣預(yù)報、校準(zhǔn)衛(wèi)星資料、氣候變化研究等方面具有重要作用[1]。探空資料也被廣泛應(yīng)用于各種極端天氣氣候事件分析[2-7]。相對于衛(wèi)星資料而言,探空資料具有時間序列長、垂直分辨率高等優(yōu)點,是監(jiān)測高空大氣變化的重要數(shù)據(jù)來源[8]。而對于20世紀(jì)70年代以前,探空資料則是高空大氣的唯一信息來源[9]。位勢高度場對于分析大氣環(huán)流演變、理解全球及區(qū)域氣候變化和極端事件的動力特征具有重要意義,并對氣候變化研究、氣候模擬的檢驗以及再分析資料的評估和改進(jìn)具有重要作用[8,10-11]。

    中國探空觀測開始于20世紀(jì)50年代初,此后迅速發(fā)展,50年代中期約有60個站,到1960年探空觀測站網(wǎng)已達(dá)到約100個站,并逐漸積累了完整、連續(xù)且覆蓋全國的歷史探空資料。中國探空觀測50年代初始階段采用美國和芬蘭等國的設(shè)備,50年代中期開始引進(jìn)、仿制蘇聯(lián)的觀測設(shè)備,并采用蘇聯(lián)的觀測規(guī)范。隨著時間的推移和我國氣象科學(xué)事業(yè)的發(fā)展,先后經(jīng)歷了從蘇式P3-049型探空儀到我國自主研制的59型探空儀,以及目前使用的L波段高空探測系統(tǒng)的多次換型。觀測規(guī)范也經(jīng)歷了多次修訂。這些變化給高空觀測資料的質(zhì)量帶來各種各樣的影響,引入不同種類的誤差。這些誤差可能會影響基于歷史探空資料的科學(xué)研究結(jié)果的準(zhǔn)確性。

    觀測儀器的變化是導(dǎo)致高空觀測資料出現(xiàn)系統(tǒng)偏差,進(jìn)而影響資料序列均一性的主要原因。大量研究結(jié)果表明,我國高空探測資料序列由于歷史上探空儀的多次更新而導(dǎo)致不連續(xù)。眾多研究表明,中國探空溫度和濕度資料序列由于儀器變更等原因?qū)е虏痪籟12-16]。馬穎等[17]利用中國70個探空站1個月的對比觀測資料,比較了59 型探空儀和L波段高空探測系統(tǒng)使用的電子探空儀溫度和位勢高度記錄的差異,分析了探空儀換型對資料一致性的影響。陶士偉等[18]以T213模式6 h預(yù)報場為背景場對L波段探空儀溫度資料的系統(tǒng)性偏差和隨機(jī)誤差進(jìn)行了評估。但這些研究主要針對探空資料中由于觀測原因造成的系統(tǒng)偏差和隨機(jī)誤差,且多限于較短時間的探空資料,僅翟盤茂[19]對中國歷史探空資料中的粗大誤差整體情況進(jìn)行了分析,并對20世紀(jì)60年代探空儀器變化給溫度和濕度時間序列帶來的不均一性進(jìn)行研究,指出不同時期的中國歷史探空資料因來源不同,錯誤率存在明顯差異。探空資料序列這些不連續(xù)將嚴(yán)重影響氣候變化事實的分析結(jié)果。

    高空觀測中的記錄整理方法(如輻射誤差訂正方法、位勢高度計算參數(shù)等)變化同樣是引起觀測資料序列不均一的重要原因。20世紀(jì)80年代以來,世界氣象組織(WMO)委托英國氣象局利用歐洲中期數(shù)值預(yù)報中心6 h預(yù)報場評估各國探空觀測的系統(tǒng)偏差和隨機(jī)誤差,分析全球高空記錄的一致性。20世紀(jì)90年代,該項工作曾報告中國探空觀測存在較大系統(tǒng)偏差和隨機(jī)誤差[20]。中國氣象局于2000—2001年對當(dāng)時59型探空儀的觀測誤差訂正方法進(jìn)行修正[21]。陶士偉等[22]依托中國國家氣象中心全球資料同化系統(tǒng)對訂正效果進(jìn)行了評估。

    本研究發(fā)現(xiàn),由于歷史上對原始觀測資料計算上的錯誤,導(dǎo)致1963年以前(含1963年)的中國歷史探空資料規(guī)定等壓面位勢高度記錄錯誤比例異常偏高,進(jìn)一步的研究發(fā)現(xiàn)錯誤記錄集中在900,800,80 hPa和60 hPa 4個規(guī)定等壓面,比同期其他規(guī)定等壓面以及1964年以后相同規(guī)定等壓面位勢高度的錯誤率高出約兩個量級,且這一問題是全國性的。此外,1976—1977年云南省的昆明(56778)、河口(56989)、麗江(56651)和騰沖(56739)4個站在250 hPa和70 hPa兩個規(guī)定等壓面依然存在位勢高度記錄錯誤比例異常偏高的問題。在研究中國探空觀測歷史沿革的基礎(chǔ)上,深入分析了這一階段特定規(guī)定等壓面位勢高度錯誤記錄異常偏多的原因,提出了訂正方案,并從多個角度對訂正效果進(jìn)行檢驗。結(jié)果表明:依據(jù)本文的訂正方案重新計算的位勢高度資料與正確的原始記錄非常接近;錯誤資料訂正后,相關(guān)氣層靜力學(xué)余差序列趨向均一,驗證了訂正方案的可行性。

    1 中國早期探空觀測位勢高度資料存在的問題

    靜力學(xué)檢查是針對探空資料中位勢高度和溫度在觀測完成后發(fā)生錯誤最有效的檢查方法。高空觀測中規(guī)定等壓面的高度由測站海拔位勢高度上累加規(guī)定等壓面之間厚度得到。臺站在記錄整理完成后,不僅報告觀測到的溫、壓、濕數(shù)據(jù),同時報告規(guī)定等壓面的位勢高度,這就構(gòu)成了探空資料質(zhì)量控制的一個重要冗余信息。靜力學(xué)檢查以在一定的誤差范圍內(nèi),能否用資料中的高度、溫度和濕度記錄還原臺站觀測時對規(guī)定等壓面之間厚度的計算作為判據(jù),檢查溫度和位勢高度記錄是否存在粗大誤差[23]。

    靜力學(xué)余差是指兩個等壓面之間的報告厚度與依據(jù)兩個等壓面的溫、壓、濕計算所得厚度的差[24],這一差值是因為等壓面資料無法還原當(dāng)時的精確觀測過程產(chǎn)生的。它能夠反映探空資料是否在觀測完成后被引入了粗大誤差,如果資料正確,靜力學(xué)余差應(yīng)在一定范圍內(nèi)變化。如果靜力學(xué)余差超過了通常的變動范圍,則提示資料可能在觀測完成之后發(fā)生了錯誤,這個通常的變動范圍即靜力學(xué)檢查的界限值。對于同一觀測站,一般在同一高度上,氣層越厚靜力學(xué)余差越大。靜力學(xué)檢查通過相鄰兩個氣層靜力學(xué)余差的大小和關(guān)系判斷資料是否存在錯誤。最典型的判斷是兩個相鄰厚度層靜力學(xué)余差都超過了界限值,如果它們符號相反,大小接近,則判斷這兩個氣層共用等壓面的位勢高度記錄為錯誤記錄;如果符號相同,大小接近則判斷是它們共用的等壓面的溫度為錯誤記錄。靜力學(xué)檢查具體方法可參閱文獻(xiàn)[23-26]。

    本文用靜力學(xué)檢查對1951—2012年中國探空觀測所有規(guī)定等壓面資料進(jìn)行了質(zhì)量控制,包括1000,925,900,850,800,700,600,500,400,300,250,200,150,100,80,70,60,50,40,30,20,15,10,7 hPa和5 hPa。其中,900,800,80 hPa和60 hPa資料1975年12月截止,250 hPa和70 hPa資料始于1976年1月。圖1是中國探空觀測規(guī)定等壓面位勢高度逐年錯誤率,可以看出,1963年以前中國探空位勢高度資料的錯誤率與其他時段存在量級上的差異:1963年以前,規(guī)定等壓面位勢高度資料錯誤率明顯偏高,部分年份高達(dá)6%;1964年以后除20世紀(jì)70年代中期錯誤率稍高外,其他年份均未超過1‰;20世紀(jì)70年代末開始使用計算機(jī)以后,位勢高度的錯誤率降低到約0.1‰,甚至更低。圖2給出了各規(guī)定等壓面的位勢高度錯誤率,位勢高度的錯誤主要集中在900,800,80 hPa和60 hPa 4個規(guī)定等壓面,錯誤率達(dá)到2%~11%,而其他等壓面錯誤率不到1‰。

    圖1 1951—2012年中國探空資料規(guī)定等壓面位勢高度逐年錯誤率Fig.1 Annual erroneous ratio of mandatory-level geopotential heights of Chinese radiosonde observations from 1951 to 2012

    圖2 1951—2012年中國探空資料規(guī)定等壓面位勢高度錯誤率Fig.2 Erroneous ratio of each mandatory-level geopotential heights of Chinese radiosonde observations from 1951 to 2012

    進(jìn)一步研究發(fā)現(xiàn),幾乎全國所有探空站與這些等壓面相關(guān)氣層的靜力學(xué)余差在1963年以前(含1963年)都遠(yuǎn)大于1963年以后。以江西省贛州站(57993)為例,其1963年12月以前和1964年1月—1975年12月兩個階段各規(guī)定等壓面之間氣層靜力學(xué)余差的標(biāo)準(zhǔn)差如圖3所示,該站1963年以前900~1000 hPa,850~900 hPa,800~850 hPa,700~850 hPa和80~100 hPa,60~80 hPa靜力學(xué)余差的標(biāo)準(zhǔn)差明顯比1963年以后大,而其他氣層的靜力學(xué)余差標(biāo)準(zhǔn)差的差異不明顯。

    圖3 1963年12月前后江西省贛州站各相鄰等壓面之間氣層靜力學(xué)余差標(biāo)準(zhǔn)差Fig.3 Standard deviation of hydrostatic residue of each mandatory level at Ganzhou Station of Jiangxi before and after Dec 1963

    此外,1976—1977年昆明、河口、麗江和騰沖4個站250 hPa和70 hPa規(guī)定等壓面的位勢高度記錄錯誤比例也異常偏高,接近50%。這兩個規(guī)定等壓面從1976年1月開始觀測,以昆明站為例,比較了其1976年1月—1977年5月和1977年6月—1999年12月兩個階段各氣層靜力學(xué)余差的變化范圍(圖4)。結(jié)果表明:1977年5月以前250~300 hPa, 200~300 hPa, 70~100 hPa和50~70 hPa 4個氣層靜力學(xué)余差的變化范圍遠(yuǎn)大于1977年5月以后,其他氣層的靜力學(xué)余差在前后兩個階段差異不明顯。河口、麗江和騰沖3個站存在同樣問題。

    圖4 1977年5月前后云南省昆明站各相鄰等壓面之間氣層靜力學(xué)余差標(biāo)準(zhǔn)差Fig.4 Standard deviation of hydrostatic residue of each mandatory level at Kunming Station of Yunnan before and after May 1977

    靜力學(xué)余差標(biāo)準(zhǔn)差偏大,表明這些氣層在相應(yīng)階段的靜力學(xué)余差的變化范圍偏大,圖5和圖6給出了更直觀的展示。

    受限于20世紀(jì)50年代和60年代中國高空觀測的技術(shù)條件,達(dá)到80 hPa或60 hPa的探空不多,且觀測誤差較大,因此,本文后續(xù)的實例分析主要以900 hPa,800 hPa以及1976年后250 hPa和70 hPa為例。

    2 問題原因

    靜力學(xué)檢查所判定的錯誤一般是觀測完成之后引入的[26]。這類錯誤通常是在觀測完成后的數(shù)據(jù)抄錄、編碼解碼和傳輸?shù)拳h(huán)節(jié)由于失誤等原因造成的,多為零散發(fā)生。但中國早期探空資料中900 hPa和800 hPa以及個別站250 hPa和70 hPa的位勢高度記錄的錯誤在特定階段集中發(fā)生。對1963年以前存在900,800,80 hPa和60 hPa等規(guī)定等壓面位勢高度錯誤的探空進(jìn)行分析發(fā)現(xiàn),這些探空的相關(guān)靜力學(xué)余差存在共同點,即這些等壓面上下相鄰氣層的靜力學(xué)余差較大,超過界限值,且大小接近,符號相反。以800 hPa的位勢高度錯誤為例,這些探空800~850 hPa和700~800 hPa兩個氣層的靜力學(xué)余差均超過界限值,且大小接近、符號相反,因此,靜力學(xué)檢查判斷是它們共用的800 hPa高度錯誤,而跨過800 hPa厚度更大的700~850 hPa的靜力學(xué)余差卻很小。昆明、河口、麗江和騰沖4個站1976—1977年存在250 hPa和70 hPa位勢高度錯誤探空的相關(guān)靜力學(xué)余差也存在同樣問題。

    圖5 1963年12月前后江西省贛州站900 hPa和800 hPa相關(guān)氣層靜力學(xué)余差變化Fig.5 Different range of hydrostatic residue for 900 hPa and 800 hPa at Ganzhou Station of Jiangxi before and after Dec 1963

    圖6 1977年5月前后云南省昆明站250 hPa和70 hPa相關(guān)氣層靜力學(xué)余差變化Fig.6 Different range of hydrostatic residue for 250 hPa and 70 hPa at Kunming Station of Yunnan before and after May 1977

    圖7給出了贛州站1963年12月前后兩個階段各跨層靜力學(xué)余差標(biāo)準(zhǔn)差的統(tǒng)計結(jié)果,從這兩個階段跨層靜力學(xué)余差的變化范圍看,該站850~1000 hPa和700~850 hPa的靜力學(xué)余差在1963年前后差異不大。圖8給出了昆明站1977年5月前后各跨層靜力學(xué)余差標(biāo)準(zhǔn)差的統(tǒng)計結(jié)果,該站200~300 hPa和50~100 hPa的靜力學(xué)余差的變化范圍在1977年5月前后也沒有明顯的差異。出現(xiàn)這一問題是在特定的時間段內(nèi)900,800,80 hPa和60 hPa位勢高度記錄的獲取方式與其他規(guī)定等壓面不同造成的。

    圖7 1963年12月前后江西省贛州站各跨層靜力學(xué)余差標(biāo)準(zhǔn)差Fig.7 Standard deviation of hydrostatic residue at Ganzhou Station of Jiangxibefore and after Dec 1963

    圖8 1977年5月前后云南省昆明站各跨層靜力學(xué)余差標(biāo)準(zhǔn)差Fig.8 Standard deviation of hydrostatic residueat at Kunming Station of Yunnan before and after May 1977

    高空觀測中記錄整理方式一般如下:首先將數(shù)據(jù)采樣點點繪并連接得到時間氣壓曲線、時間溫度曲線和時間濕度曲線,再通過氣壓與時間的對應(yīng)關(guān)系從時溫曲線和時濕曲線上讀取某個氣壓所對應(yīng)的溫度和濕度,得到該等壓面的溫度和濕度。然后按照定積分的中值定理,在時溫曲線和時濕曲線上用等面積法求得兩個規(guī)定等壓面之間氣層的平均溫度和平均濕度,再用靜力學(xué)公式計算就能夠得到相應(yīng)氣層的厚度,累加得到各等壓面高度。各規(guī)定等壓面高度點繪成時間高度曲線,可進(jìn)一步獲取規(guī)定高度風(fēng)等信息。

    高空觀測過程中需要進(jìn)行大量計算。中國高空觀測在20世紀(jì)70年代末開始使用計算機(jī)前,由于缺乏計算條件,規(guī)定等壓面之間厚度通過查表得到。1963年以前觀測規(guī)范[27-30]給出的規(guī)定等壓面間厚度表里只有850~1000 hPa,700~850 hPa等氣層的厚度,900~1000 hPa和800~850 hPa等氣層的厚度無法通過查表直接得到,而是在查表得到1000,850 hPa和700 hPa的高度后,在時間高度曲線上通過內(nèi)插的方法人工讀取,主觀因素的影響較大,容易發(fā)生錯誤。反映在歷史探空資料各層的靜力學(xué)余差上,便出現(xiàn)了圖3~6所揭示的現(xiàn)象,部分相鄰規(guī)定等壓面間靜力學(xué)余差明顯大,而跨過這些等壓面厚度更大的氣層的靜力學(xué)余差反而小。

    1963年7月中央氣象局頒布了新的《高空氣象觀測規(guī)范(使用P3-049 型探空儀)》[31],明確了規(guī)定等壓面包括1000,900,850,800,700,600,500,400,300,200,150,100,80,60,50,40,30,20,15 hPa和10 hPa,在規(guī)定等壓面厚度表中增加900~1000 hPa,850~900 hPa,800~850 hPa,700~800 hPa 等氣層的厚度。1964年8月印發(fā)的《高空氣象觀測規(guī)范說明及問題綜合解答(使用P3-049 型探空儀)》[32]中再次強(qiáng)調(diào)各規(guī)定等壓面的高度必須通過計算得到。絕大多數(shù)臺站從1964年1月起,900~1000 hPa,850~900 hPa,800~850 hPa 和700~800 hPa等氣層的靜力學(xué)余差變化范圍減小到與此后一致的水平,但北京站和貴陽站分別到1964年2月和1965年1月才開始按照新規(guī)范獲取這兩個等壓面的位勢高度。

    1976年觀測規(guī)范[33]新增加了250 hPa和70 hPa兩個規(guī)定等壓面,昆明、河口、麗江和騰沖4個站1976—1977年250 hPa和70 hPa兩個規(guī)定等壓面的位勢高度出現(xiàn)了類似早期900 hPa和800 hPa的問題,最長持續(xù)19個月。產(chǎn)生這一問題的具體原因已經(jīng)無從查考,可能與當(dāng)?shù)赜^測規(guī)范沒有及時更新貫徹,觀測人員習(xí)慣性地沿用1963年以前對900 hPa和800 hPa位勢高度記錄的處理方法有關(guān)。

    3 訂正方案

    由于臺站原始探空觀測信息已難以獲取,僅規(guī)定等壓面的溫濕度不足以反映大氣垂直結(jié)構(gòu)的所有細(xì)節(jié)。特性層捕捉了大氣層結(jié)中等溫層和逆溫層的起始位置、終止位置以及其他溫度和濕度梯度的顯著轉(zhuǎn)折點[31-35],將特性層的溫、濕度記錄引入厚度計算能夠較精細(xì)反映真實的大氣垂直結(jié)構(gòu),得出比較準(zhǔn)確的計算高度。中國早期探空觀測規(guī)范厚度表包含的1000,850,300 hPa和100 hPa等規(guī)定等壓面的位勢高度通過查表計算直接求得,得到的高度數(shù)值較為準(zhǔn)確。因此,可在這些規(guī)定等壓面高度的基礎(chǔ)上,充分利用特性層資料,采用靜力學(xué)公式對900,800,250,80,70 hPa和60 hPa的位勢高度進(jìn)行重新計算,完成對位勢高度的訂正。

    根據(jù)壓高公式[35],等壓面第(i-1)層和i之間的厚度為

    (1)

    假設(shè)兩個規(guī)定等壓面第(i-1)層和第i層之間有m個特性層,則這兩個規(guī)定等壓面之間的氣層被分隔為(m+1)個薄的氣層。為表述方便,將上下兩個規(guī)定等壓面和特性層統(tǒng)一重新編號為0到(m+1)。這些薄的氣層內(nèi)的溫濕度隨高度接近線性變化,該氣層內(nèi)的平均虛溫可以由上下層的虛溫求算術(shù)平均得到,并通過式(1)求得各個薄層的厚度,在規(guī)定等壓面第(i-1)層的高度上累加各個薄層的厚度,可以準(zhǔn)確地計算出規(guī)定第i層等壓面的位勢高度。

    結(jié)合特性層資料,規(guī)定等壓面第(i-1)層和第i層之間厚度的計算公式為

    (2)

    我國2000年以前厚度計算公式中,干空氣比氣體常數(shù)Rd和重力加速度g的數(shù)值與目前所使用的數(shù)值略有不同。為了與同期其他規(guī)定等壓面位勢高度的計算保持一致,本文的厚度計算中這兩個參數(shù)也采用當(dāng)時的數(shù)值,即Rd取287.04 J/(kg·K),g取9.80 m/s2。

    此外,中國探空觀測對于400 hPa 以上的規(guī)定等壓面位勢高度計算時不進(jìn)行虛溫訂正[31-35],因此,在重新計算250 hPa和70 hPa的位勢高度時,公式為

    (3)

    對這些規(guī)定等壓面位勢高度通過重新計算來訂正,要首先保證被訂正的等壓面溫度和濕度記錄、下面相鄰規(guī)定等壓面位勢高度、溫度和濕度記錄,以及兩個規(guī)定等壓面之間特性層記錄的正確。因此,首先采用氣候?qū)W界限值檢查、溫度的超絕熱遞減率檢查、靜力學(xué)檢查等方法對需要參與訂正計算的探空資料進(jìn)行質(zhì)量控制,只在這些資料正確的情況下進(jìn)行相應(yīng)探空的訂正計算。

    4 對訂正方案的效果檢驗

    4.1 獨立樣本檢驗

    為了驗證訂正方案的有效性,首先采用贛州站和昆明站資料進(jìn)行效果檢驗。檢驗方法的設(shè)計思路是如果采用上述訂正方案重新計算的某等壓面的位勢高度(簡稱為計算高度),與觀測時根據(jù)氣壓、溫度和濕度曲線正確計算得出的原始位勢高度記錄數(shù)值相差不大,那么訂正方案是可行的。由圖5和圖6可知,1964—1975年贛州站探空資料900 hPa和800 hPa的位勢高度不存在系統(tǒng)性計算錯誤問題,而1977年6月以后昆明站250 hPa和70 hPa的位勢高度資料也不存在類似問題。為此,分別采用1964年1月—1975年12月贛州站和1977年6月—1998年12月昆明站探空資料對訂正方案進(jìn)行檢驗。

    圖9和圖10分別是1964年1月—1975年12月贛州站900 hPa和800 hPa以及1977年6月—1998年12月昆明站250 hPa和70 hPa的計算高度與原始高度記錄差值的概率分布。由圖9可以看到,贛州站900 hPa和800 hPa的計算高度的誤差均在±5 gpm范圍內(nèi),其中分別有56.7%和44.4%的計算高度與原始高度記錄完全一致。對昆明站250 hPa位勢高度重新計算也取得了很好的效果,所有計算高度的誤差在±5 gpm范圍內(nèi),其中,79.3%的計算高度誤差在±2 gpm范圍內(nèi)。但對昆明站70 hPa位勢高度的重新計算則誤差稍大,有近20%的計算高度誤差超過±5 gpm,但均在±19 gpm范圍內(nèi)。70~100 hPa氣層的厚度約為2300 gpm,遠(yuǎn)超過900~1000 hPa,800~850 hPa和250~300 hPa的厚度,這是用100 hPa的位勢高度重新計算70 hPa位勢高度誤差稍大的原因之一;此外,高空觀測規(guī)范規(guī)定對流層頂以上的特性層選取標(biāo)準(zhǔn)與對流層相比要寬泛一些[26-34],只捕捉非常顯著的溫度拐點,這是采用100 hPa高度重新計算70 hPa高度誤差稍大的另一個原因。

    圖9 江西省贛州站計算高度與原始高度差的概率密度分布Fig.9 Probability density distribution of difference between calculated and original geopotential heights at Ganzhou Station of Jiangxi

    圖10 云南省昆明站計算高度與原始高度差的概率密度分布Fig.10 Probability density distribution of difference between calculated and original geopotential heights at Kunming Station of Yunnan

    4.2 訂正效果檢驗

    對于確實需要訂正的資料無法獲取其應(yīng)有的位勢高度真值,只能采用其他方法進(jìn)行訂正效果檢驗。以800 hPa位勢高度的訂正為例,采用在850 hPa位勢高度的基礎(chǔ)上累加根據(jù)溫、壓、濕計算出的800~850 hPa的厚度的方法,得到800 hPa訂正位勢高度。如果訂正合理,700~800 hPa之間氣層的靜力學(xué)余差應(yīng)減小。圖11是贛州站1963年12月以前900 hPa和800 hPa位勢高度訂正后850~900 hPa和700~800 hPa的靜力學(xué)余差??梢钥闯?,與圖5相比,訂正后1963年以前700~800 hPa的靜力學(xué)余差明顯減小,并與1963年以后相應(yīng)氣層的靜力學(xué)余差表現(xiàn)一致,證明訂正方案可行。圖12是昆明站1977年5月以前250 hPa和70 hPa位勢高度訂正后相關(guān)氣層的靜力學(xué)余差與后期的比較。

    圖11 江西省贛州站900 hPa和800 hPa相關(guān)氣層靜力學(xué)余差(1963年12月以前為訂正后位勢高度)Fig.11 Hydrostatic residue of 900 hPa and 800 hPa at Ganzhou Station of Jiangxi (with corrected geopotential heights before Dce 1963)

    圖12 云南省昆明站250 hPa和70 hPa位勢高度相關(guān)氣層靜力學(xué)余差(1977年5月以前為訂正后位勢高度)Fig.12 Hydrostatic residue of 250 hPa and 70 hPa at Kunming Station of Yunnan (with corrected geopotential heights before May 1977)

    5 結(jié) 論

    研究表明:

    1) 1963年以前的中國探空資料中900,800,80 hPa和60 hPa等規(guī)定等壓面位勢高度在靜力學(xué)檢查中發(fā)現(xiàn)大量錯誤,錯誤率達(dá)到2%~11%。原因是當(dāng)時觀測規(guī)范中尚不包含根據(jù)靜力學(xué)公式制作的這些規(guī)定等壓面在各種溫、濕條件下的厚度表。與其他規(guī)定等壓面不同,這些等壓面相關(guān)氣層的厚度不是查表得到,而是在時間高度曲線上通過內(nèi)插的方法人工讀取,受主觀因素影響,誤差較大。

    2) 1976—1977年云南省昆明、河口、麗江和騰沖4個站250 hPa和70 hPa兩個規(guī)定等壓面相關(guān)的位勢高度記錄錯誤率也異常偏高,接近50%,相關(guān)氣層靜力學(xué)余差遠(yuǎn)大于其他氣層。

    3) 本文依據(jù)靜力學(xué)原理,給出了結(jié)合特性層資料重新計算相應(yīng)規(guī)定等壓面位勢高度的誤差訂正方案。檢驗結(jié)果表明:依據(jù)訂正方案重新計算的位勢高度80%以上與正常的原始記錄誤差不超過±5 gpm,能夠比較精確地還原臺站觀測時對位勢高度的正確計算;錯誤資料訂正后,相關(guān)氣層靜力學(xué)余差序列趨向均一,表明訂正方案合理、有效。

    [1] Office of the Federal Coordinator for Meteorological Service and Supporting Research.Federal Meteorological Handbook No.3, Rawinsonde and Pibal Observation. Washington D C,1997.

    [2] 祁秀香,鄭永光.2007年夏季我國深對流活動時空分布特征.應(yīng)用氣象學(xué)報,2009,20(3):286-294.

    [3] 孫虎林,羅亞麗,張人禾.2009年6月3—4日黃淮地區(qū)強(qiáng)颮線天氣過程成熟階段特征分析.大氣科學(xué),2011,35(1):1-16.

    [4] 江吉喜,項續(xù)康,范梅珠.青藏高原夏季中尺度強(qiáng)對流系統(tǒng)的時空分布.應(yīng)用氣象學(xué)報,1996,7(4):473-478.

    [5] 張騰飛,魯亞斌,張杰,等.2000年以來云南4次強(qiáng)降雪過程的對比分析.應(yīng)用氣象學(xué)報,2007,18(1):64-72.

    [6] 王迎春,錢婷婷,鄭永光.北京連續(xù)降雪過程分析.應(yīng)用氣象學(xué)報,2004,15(1):58-65.

    [7] 葉晨,王建捷,張文龍.北京2009年“1101”暴雪的形成機(jī)制.應(yīng)用氣象學(xué)報,2011,22(4):398-409.

    [8] Marshall G J.Trends in Antarctic geopotential height and temperature:A comparison between radiosonde and NCEP-NCAR reanalysis data.JClimate,2002,15(6):659-673.

    [9] Durre I,Vose R S,Wuertz D B.Overview of the integrated global radiosonde archive.JClimate,2006,19(1):53-68.

    [10] Klein W H,Dai Y.Reconstruction of monthly mean 700-mb height from surface data by reverse specification.JClimate,1998,11(8):2136-2146.

    [11] Gong D Y,Drange H,Gao Y.Reconstruction of Northern Hemisphere 500 hPa geopotential heights back to the late 19th century.TheoreticalandAppliedClimatology,2007,90(1-2):83-102.

    [12] Zhai P M,Eskridge R E.Analyses of inhomogeneities in radiosonde temperature and humidity time Series.JClimate,1996,9(4):884-894.

    [13] Guo Y J,Thorne P W,McCarthy M P,et al.Radiosonde temperature trends and their uncertainties over eastern China.IntJClimatol,2008,28:1269-1281.

    [14] Guo Y J,Ding Y H.Long-term free-atmosphere temperature trends in China derived from homogenized in situ radiosonde temperature series.JClimate,2009,22(4):1037-1051.

    [15] Dai A,Wang J,Thorne P W,et al.A new approach to homogenize daily radiosonde humidity data.JClimate,2011,24(4):965-991.

    [16] Zhao T,Dai A,Wang J.Trends in tropospheric humidity from 1970 to 2008 over China from a homogenized radiosonde dataset.JClimate,2012,25(13):4549-4567.

    [17] 馬穎,姚雯,黃炳勛.59型和L波段探空儀溫度和位勢高度記錄對比.應(yīng)用氣象學(xué)報,2010,21(2):214-220.

    [18] 陶士偉,張躍堂.全球觀測資料質(zhì)量監(jiān)視評估.氣象,2006,32(6):53-58.

    [19] 翟盤茂.中國歷史探空資料中的一些過失誤差及偏差問題.氣象學(xué)報,1997,55(5):563-572.

    [20] Oakley T.Report by the Rapporteur on Radiosonde Compatibility Monitoring.WMO/TD-No.587,1993.

    [21] Guo Y,Huang B,Hu D,et al.Correction for Bias of Chinese Upper-air Measurements.WMO Teco,2002.

    [22] 陶士偉,陳曉紅,龔建東.L波段探空儀溫度資料誤差分析.氣象,2006,32(10):46-51.

    [23] Collins W G,Gandin L S.Comprehensive hydrostatic quality control at National Meteorological Center.MonWeaRev,1990,118(12):2752-2767.

    [24] Collins W G.The operational complex quality control of radiosonde height and temperatures at the national centers for environmental prediction.Part I,description of the method.JAppliedMeteor,2001,40(2):137-151.

    [25] Gandin L S.Complex quality control of meteorological observations.MonWeaRev,1990,116(5):1137-1155.

    [26] 翟盤茂,周琴芳.高空資料質(zhì)量控制的CHQC方法.氣象科技,1993,11(3):92-97.

    [27] 中央人民政府人民革命軍事委員會氣象局.蘇聯(lián)預(yù)報工作規(guī)范(第1部),天氣圖和高空圖的繪制.1953.

    [28] 中央氣象局.蘇聯(lián)部長會議水文氣象總部水文氣象站哨規(guī)范(第四分冊),測站高空氣象觀測.1956.

    [29] 中央氣象局.蘇式(P3-049)型探空觀測規(guī)范.1959.

    [30] 中央氣象局.蘇式(P3-049)型探空觀測規(guī)范(重印本).1959.

    [31] 中央氣象局.高空氣象觀測規(guī)范(使用P3-049型探空儀).1963.

    [32] 中央氣象局觀象臺.高空氣象觀測規(guī)范(使用P3-049型探空儀)說明及問題綜合解答.1964.

    [33] 中央氣象局.高空氣象觀測規(guī)范(使用59型探空儀).1976.

    [34] 中國氣象局監(jiān)測網(wǎng)絡(luò)司.高空氣象探測手冊(59-701微機(jī)處理系統(tǒng)部分).北京:氣象出版社,2001.

    [35] 中國氣象局監(jiān)測網(wǎng)絡(luò)司.L波段(I型)高空氣象探測系統(tǒng)業(yè)務(wù)操作手冊.北京:氣象出版社,2005.

    Correcting Geopotential Height Errors of Some Mandatory Levels of Chinese Historic Radiosonde Observations

    Ruan Xin1)Xiong Anyuan1)Hu Kaixi1)Liang Xiuhui2)Wang Zhiwen3)Yang Yanru1)Zou Fengling1)

    1)(NationalMeteorologicalInformationCenter,Beijing100081)2)(HebeiProvincialMeteorologicalBureau,Shijiazhuang050021)3)(BeijingMeteorologicalBureau,Beijing100089)

    Radiosonde observations are crucial in weather forecast and upper-air climate research. Due to their high vertical resolution, they are also important for calibration and validation of satellite temperature and water vapor retrievals. Quality control of Chinese radiosonde observations of 1951-2012 is conducted with hydrostatic check, and the result reveals that ratios of erroneous geopotential height (GPH) of 900 hPa, 800 hPa, 80 hPa and 60 hPa before 1963 are 2%-11%, which are much higher than other levels. The ratios of erroneous GPH of 250 hPa and 70 hPa in mid of the 1970s of a few stations in Yunnan Province are also much higher than other levels. Changes with observation code and practice in Chinese radiosonde observations are thoroughly investigated. Then the possible cause of the phenomenon is analyzed and validated with observations.

    The height of each pressure level in sounding is calculated by accumulating the thickness of layers between each pressure levels below on the base of station elevation. And the thickness is determined by hydrostatic equation automatically nowadays, but due to the limitation of calculation capacity before the 1980s, they are available in ready-made GPH-table. However, 900 hPa, 800 hPa, 80 hPa and 60 hPa are not included in the GPH-table possibly because they are not directly used in weather charts then. So GPHs of levels in question are obtained through interpolation on time-GPH chart after GPH of other levels are available on the chart through looking up the GPH-table. This practice makes GPH of levels in question vulnerable to human mistakes, thus many records contain errors. The explanation to the cause of the problem is validated with observations.

    Furthermore, the method to correct the erroneous GPH record is proposed, which is to recalculate it based on hydrostatic equation with correct observations of other levels from the same sounding. The correcting method is validated since the originally correct GPH records are reproduced accurately. For 900 hPa and 800 hPa of Ganzhou Station, differences between reproduced GPHs and original ones are less than 5 gpm, and about 50% of reproduced GPHs are exactly the same as original ones. Differences for 250 hPa of Kunming Station are all less than 5 gpm either, but differences are slightly larger for 70 hPa, possibly because the layer between 70-100 hPa is as large as 2300 gpm approximately, and significant levels are not enough to reveal all the details of temperature profile. Besides, the correcting method is validated from another aspect, since after the erroneous GPH are corrected, the relevant hydrostatic residue series become consistent.

    radiosonde observations; geopotential height; correction

    10.11898/1001-7313.20150301

    國家重點基礎(chǔ)研究發(fā)展計劃(2010CB951602),公益性行業(yè)(氣象)科研專項(GYHY201406001),國家科技支撐計劃(2012BAC22B05),中國氣象局氣候變化專項(CCSF201330)

    阮新,熊安元,胡開喜,等. 中國歷史探空資料部分等壓面位勢高度錯誤訂正. 應(yīng)用氣象學(xué)報,2015,26(3):257-267.

    2014-09-10收到, 2015-03-02收到再改稿。

    * email: ruanx@cma.gov.cn

    猜你喜歡
    位勢探空氣層
    探空氣球為高空氣象“把脈”
    含Hardy位勢的非線性Schr?dinger-Poisson方程正規(guī)化解的多重性
    一類帶強(qiáng)制位勢的p-Laplace特征值問題
    基于孔、裂隙理論評價致密氣層
    含變號位勢的ρ-Kirchhoff型方程組無窮多個高能量解的存在性
    含位勢的非線性雙調(diào)和方程解的存在性
    涪陵頁巖氣層構(gòu)造變形量化研究獲新進(jìn)展
    錄井工程(2017年1期)2017-07-31 17:44:42
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對比分析
    船底凹槽氣層流動特征及影響因素研究
    蘇里格氣田致密砂巖氣層識別難點及方法評述
    亚洲综合色网址| 久久精品国产a三级三级三级| 成人毛片a级毛片在线播放| 高清在线视频一区二区三区| 日韩不卡一区二区三区视频在线| 亚洲精品美女久久久久99蜜臀 | 成年美女黄网站色视频大全免费| av在线老鸭窝| 国产精品嫩草影院av在线观看| 亚洲内射少妇av| 这个男人来自地球电影免费观看 | 亚洲,欧美精品.| 精品亚洲乱码少妇综合久久| 男人爽女人下面视频在线观看| 欧美国产精品一级二级三级| 国产片内射在线| www.熟女人妻精品国产| av一本久久久久| 国产视频首页在线观看| 看非洲黑人一级黄片| 亚洲av欧美aⅴ国产| 90打野战视频偷拍视频| 热99久久久久精品小说推荐| 老汉色∧v一级毛片| 五月天丁香电影| 亚洲av电影在线观看一区二区三区| 亚洲欧美成人综合另类久久久| 永久免费av网站大全| 国产成人精品福利久久| 免费少妇av软件| 国产精品久久久久久久久免| 黄色视频在线播放观看不卡| 黄片播放在线免费| 亚洲av中文av极速乱| h视频一区二区三区| 纯流量卡能插随身wifi吗| 老司机影院成人| 久久国产亚洲av麻豆专区| 婷婷色综合大香蕉| 亚洲av欧美aⅴ国产| 叶爱在线成人免费视频播放| 久久久欧美国产精品| 十八禁高潮呻吟视频| 国产亚洲一区二区精品| 中文字幕精品免费在线观看视频| 久久影院123| 亚洲av福利一区| 久久精品国产自在天天线| 天堂中文最新版在线下载| 亚洲男人天堂网一区| 伊人久久国产一区二区| 女的被弄到高潮叫床怎么办| 亚洲av免费高清在线观看| 国产综合精华液| 亚洲精品国产一区二区精华液| av免费在线看不卡| 久久99一区二区三区| 亚洲经典国产精华液单| 久久久欧美国产精品| 欧美成人午夜免费资源| 人妻少妇偷人精品九色| 久久97久久精品| 亚洲欧美成人综合另类久久久| 成人国产av品久久久| 国产免费视频播放在线视频| 精品久久蜜臀av无| 精品人妻偷拍中文字幕| 久久鲁丝午夜福利片| 免费黄频网站在线观看国产| 青春草国产在线视频| 老司机影院成人| 91精品国产国语对白视频| 亚洲国产欧美在线一区| 视频区图区小说| 国产精品国产三级国产专区5o| 久久久久久久久免费视频了| 亚洲av男天堂| 91在线精品国自产拍蜜月| 电影成人av| 久久久久精品人妻al黑| 午夜影院在线不卡| 一级毛片我不卡| 国产精品人妻久久久影院| 亚洲av日韩在线播放| 欧美人与善性xxx| 女人高潮潮喷娇喘18禁视频| 日韩一区二区三区影片| 中文字幕精品免费在线观看视频| 黑人欧美特级aaaaaa片| 中文字幕人妻熟女乱码| 免费观看在线日韩| 欧美激情极品国产一区二区三区| 黄片小视频在线播放| 乱人伦中国视频| a级毛片黄视频| 精品第一国产精品| 亚洲熟女精品中文字幕| 久久精品人人爽人人爽视色| 男女无遮挡免费网站观看| kizo精华| 女人高潮潮喷娇喘18禁视频| 女性被躁到高潮视频| 99久久综合免费| 人妻人人澡人人爽人人| 国产男女内射视频| 亚洲欧洲国产日韩| 日韩人妻精品一区2区三区| 成年动漫av网址| 狠狠精品人妻久久久久久综合| 婷婷色综合www| 我的亚洲天堂| 欧美另类一区| 欧美精品国产亚洲| 久久久久久伊人网av| 日韩不卡一区二区三区视频在线| 91午夜精品亚洲一区二区三区| 日本-黄色视频高清免费观看| 国产伦理片在线播放av一区| 女的被弄到高潮叫床怎么办| 亚洲五月色婷婷综合| 亚洲欧美中文字幕日韩二区| 久久韩国三级中文字幕| 国产成人av激情在线播放| 综合色丁香网| 青青草视频在线视频观看| 丰满乱子伦码专区| 国产人伦9x9x在线观看 | 精品久久久精品久久久| 久久 成人 亚洲| 亚洲,一卡二卡三卡| 久久 成人 亚洲| 亚洲国产最新在线播放| 精品一区二区三区四区五区乱码 | 香蕉国产在线看| 亚洲av中文av极速乱| 最新的欧美精品一区二区| 久久精品亚洲av国产电影网| 日本av免费视频播放| 国产免费福利视频在线观看| 国产成人精品无人区| 亚洲精品aⅴ在线观看| 一级毛片黄色毛片免费观看视频| 美女脱内裤让男人舔精品视频| 免费女性裸体啪啪无遮挡网站| 黄频高清免费视频| 亚洲av中文av极速乱| 欧美日韩成人在线一区二区| 精品国产一区二区久久| 亚洲欧美成人综合另类久久久| 日韩伦理黄色片| 精品一区二区三区四区五区乱码 | 高清不卡的av网站| 丝袜美足系列| 我要看黄色一级片免费的| 人妻人人澡人人爽人人| 亚洲av免费高清在线观看| 黄色视频在线播放观看不卡| a 毛片基地| 亚洲精品成人av观看孕妇| 日韩中文字幕欧美一区二区 | av在线观看视频网站免费| 欧美日韩视频精品一区| 欧美日韩一级在线毛片| 一区二区av电影网| 日本欧美国产在线视频| 亚洲精品一区蜜桃| 精品酒店卫生间| 日本免费在线观看一区| 一二三四在线观看免费中文在| 亚洲av日韩在线播放| 永久网站在线| 男的添女的下面高潮视频| 国产精品秋霞免费鲁丝片| 午夜久久久在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲情色 制服丝袜| 黄色怎么调成土黄色| 国产成人免费无遮挡视频| av卡一久久| 国产精品久久久久成人av| 高清欧美精品videossex| 精品亚洲成a人片在线观看| 成年av动漫网址| 亚洲美女搞黄在线观看| 国产xxxxx性猛交| 午夜免费男女啪啪视频观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲经典国产精华液单| 美女脱内裤让男人舔精品视频| 成人毛片60女人毛片免费| 青草久久国产| 日韩中字成人| 一本色道久久久久久精品综合| 97在线视频观看| 啦啦啦在线免费观看视频4| videossex国产| 国产精品久久久久久精品电影小说| 少妇熟女欧美另类| 亚洲国产成人一精品久久久| 亚洲欧美色中文字幕在线| 亚洲国产精品999| 精品国产露脸久久av麻豆| 午夜福利网站1000一区二区三区| 91国产中文字幕| 黑丝袜美女国产一区| 三级国产精品片| 国产一级毛片在线| 午夜福利,免费看| 国产高清不卡午夜福利| 这个男人来自地球电影免费观看 | 欧美日韩一区二区视频在线观看视频在线| 中文天堂在线官网| 国产极品粉嫩免费观看在线| 精品久久蜜臀av无| 成人漫画全彩无遮挡| 久久韩国三级中文字幕| 黄色视频在线播放观看不卡| 晚上一个人看的免费电影| 男女午夜视频在线观看| 国产又色又爽无遮挡免| 天堂中文最新版在线下载| 99热网站在线观看| 日日摸夜夜添夜夜爱| 九九爱精品视频在线观看| 欧美精品高潮呻吟av久久| 国产伦理片在线播放av一区| 九草在线视频观看| 高清黄色对白视频在线免费看| 搡女人真爽免费视频火全软件| 在线观看三级黄色| 国产成人免费无遮挡视频| 国产老妇伦熟女老妇高清| 久久久久久免费高清国产稀缺| 国产亚洲av片在线观看秒播厂| 日韩欧美一区视频在线观看| 亚洲中文av在线| 国产免费福利视频在线观看| 如日韩欧美国产精品一区二区三区| 久久久久精品久久久久真实原创| 精品久久蜜臀av无| 欧美日韩视频高清一区二区三区二| 丰满少妇做爰视频| 97精品久久久久久久久久精品| 亚洲精品国产av蜜桃| 在线观看免费视频网站a站| 1024香蕉在线观看| 久久国产精品大桥未久av| 91aial.com中文字幕在线观看| 爱豆传媒免费全集在线观看| 亚洲情色 制服丝袜| 国产精品熟女久久久久浪| 蜜桃国产av成人99| 男人添女人高潮全过程视频| 伦理电影大哥的女人| 日韩人妻精品一区2区三区| 国产麻豆69| 国产人伦9x9x在线观看 | 成年女人在线观看亚洲视频| 老司机影院毛片| 日本av免费视频播放| 大片电影免费在线观看免费| 嫩草影院入口| 欧美变态另类bdsm刘玥| 久久毛片免费看一区二区三区| 国产精品.久久久| 在线观看免费日韩欧美大片| 久久久久久久久免费视频了| 不卡av一区二区三区| 男女午夜视频在线观看| 黄色配什么色好看| 亚洲精品久久久久久婷婷小说| 人妻一区二区av| 午夜福利在线观看免费完整高清在| 日本色播在线视频| 亚洲欧洲精品一区二区精品久久久 | 久久精品aⅴ一区二区三区四区 | 观看av在线不卡| 男女高潮啪啪啪动态图| 亚洲欧美一区二区三区久久| 成人二区视频| 日韩欧美一区视频在线观看| 国产男女超爽视频在线观看| 免费观看a级毛片全部| 黑人猛操日本美女一级片| 亚洲婷婷狠狠爱综合网| 精品少妇一区二区三区视频日本电影 | 波多野结衣一区麻豆| 日韩av免费高清视频| 热re99久久国产66热| 国产黄色视频一区二区在线观看| 国产亚洲av片在线观看秒播厂| www.av在线官网国产| 亚洲一区二区三区欧美精品| 99久久中文字幕三级久久日本| 黄色毛片三级朝国网站| 天堂俺去俺来也www色官网| 免费少妇av软件| 你懂的网址亚洲精品在线观看| 日本色播在线视频| 精品一区在线观看国产| 国产免费视频播放在线视频| 三级国产精品片| 国产又爽黄色视频| 最近中文字幕2019免费版| 国产欧美日韩综合在线一区二区| 丰满乱子伦码专区| 美国免费a级毛片| 九九爱精品视频在线观看| 国产成人一区二区在线| 亚洲av欧美aⅴ国产| 十分钟在线观看高清视频www| 看免费av毛片| 啦啦啦啦在线视频资源| 国产伦理片在线播放av一区| 一级片免费观看大全| 亚洲国产欧美在线一区| 欧美日韩国产mv在线观看视频| 国产一级毛片在线| 亚洲国产精品999| 一级片免费观看大全| 18禁裸乳无遮挡动漫免费视频| av卡一久久| 熟妇人妻不卡中文字幕| 在线精品无人区一区二区三| 亚洲av中文av极速乱| 国产精品三级大全| av女优亚洲男人天堂| 亚洲图色成人| 黄色配什么色好看| 极品人妻少妇av视频| 欧美日韩国产mv在线观看视频| 十八禁高潮呻吟视频| 国产精品二区激情视频| 这个男人来自地球电影免费观看 | 亚洲国产欧美在线一区| 黄色视频在线播放观看不卡| 国产色婷婷99| 国产精品不卡视频一区二区| 亚洲一码二码三码区别大吗| 国产男女内射视频| 老司机亚洲免费影院| 韩国高清视频一区二区三区| 久久久欧美国产精品| 九草在线视频观看| freevideosex欧美| 色视频在线一区二区三区| 亚洲精品中文字幕在线视频| 看免费成人av毛片| 精品一区二区三区四区五区乱码 | 高清在线视频一区二区三区| 亚洲欧美精品自产自拍| 高清av免费在线| 最近的中文字幕免费完整| 国产毛片在线视频| 亚洲色图 男人天堂 中文字幕| 中文字幕av电影在线播放| 老熟女久久久| 大话2 男鬼变身卡| 国产视频首页在线观看| 青春草国产在线视频| 高清视频免费观看一区二区| 日本色播在线视频| 久久 成人 亚洲| 97在线人人人人妻| 黄频高清免费视频| 2018国产大陆天天弄谢| 最近的中文字幕免费完整| 久久亚洲国产成人精品v| 两性夫妻黄色片| 观看美女的网站| 欧美日韩亚洲高清精品| 伦理电影大哥的女人| 人人妻人人澡人人看| 婷婷色综合www| 考比视频在线观看| 美女国产高潮福利片在线看| 电影成人av| 亚洲人成77777在线视频| 精品国产乱码久久久久久小说| 色视频在线一区二区三区| 亚洲成人一二三区av| 久久鲁丝午夜福利片| 亚洲天堂av无毛| 亚洲国产av影院在线观看| 精品少妇久久久久久888优播| 999久久久国产精品视频| 亚洲久久久国产精品| 成人国产麻豆网| 国产综合精华液| 美女视频免费永久观看网站| 尾随美女入室| 国产精品一国产av| 精品第一国产精品| 在线免费观看不下载黄p国产| 精品一区二区三卡| 天天操日日干夜夜撸| 中文字幕最新亚洲高清| 最近的中文字幕免费完整| 精品一品国产午夜福利视频| 交换朋友夫妻互换小说| 成人国产av品久久久| 成人国产av品久久久| 少妇的丰满在线观看| 曰老女人黄片| 丝袜喷水一区| 免费观看无遮挡的男女| 精品国产乱码久久久久久小说| 美女福利国产在线| 国产精品蜜桃在线观看| 两个人免费观看高清视频| 老鸭窝网址在线观看| 日韩中文字幕视频在线看片| 天天躁日日躁夜夜躁夜夜| 丝瓜视频免费看黄片| 99国产综合亚洲精品| 多毛熟女@视频| 亚洲精品第二区| 亚洲视频免费观看视频| 在线观看免费日韩欧美大片| 免费观看无遮挡的男女| 久久久久久人人人人人| 日韩,欧美,国产一区二区三区| 国产一区二区三区av在线| av卡一久久| 日韩精品免费视频一区二区三区| 在线看a的网站| 搡老乐熟女国产| 婷婷色av中文字幕| 欧美日韩综合久久久久久| 精品久久蜜臀av无| 亚洲精品第二区| 亚洲人成网站在线观看播放| 狠狠婷婷综合久久久久久88av| 精品国产一区二区久久| 一级毛片我不卡| 热99久久久久精品小说推荐| av在线观看视频网站免费| av电影中文网址| av在线app专区| 久久免费观看电影| 制服人妻中文乱码| 极品少妇高潮喷水抽搐| 亚洲欧美一区二区三区黑人 | 人妻系列 视频| 在线观看免费视频网站a站| 亚洲精品视频女| 免费av中文字幕在线| 久久午夜福利片| 巨乳人妻的诱惑在线观看| 中文欧美无线码| 亚洲色图综合在线观看| 久久精品熟女亚洲av麻豆精品| 18禁裸乳无遮挡动漫免费视频| 亚洲欧美精品综合一区二区三区 | 波野结衣二区三区在线| 国产精品秋霞免费鲁丝片| 亚洲四区av| 亚洲国产日韩一区二区| 精品国产一区二区久久| 久久韩国三级中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 高清黄色对白视频在线免费看| 国产精品 欧美亚洲| 久热久热在线精品观看| 女性生殖器流出的白浆| 最近最新中文字幕免费大全7| 亚洲欧美一区二区三区国产| 日韩熟女老妇一区二区性免费视频| 亚洲av成人精品一二三区| 丝袜美足系列| 日韩一区二区视频免费看| 在线观看免费日韩欧美大片| 91在线精品国自产拍蜜月| 在线观看一区二区三区激情| 成人毛片60女人毛片免费| 国产精品亚洲av一区麻豆 | 欧美日韩亚洲国产一区二区在线观看 | 2018国产大陆天天弄谢| a级毛片黄视频| 亚洲,欧美精品.| 色婷婷av一区二区三区视频| 9191精品国产免费久久| av在线老鸭窝| 亚洲精品久久久久久婷婷小说| av卡一久久| 国产精品麻豆人妻色哟哟久久| 2022亚洲国产成人精品| 我要看黄色一级片免费的| 另类亚洲欧美激情| 人人妻人人澡人人看| 亚洲精品第二区| 国产一区二区激情短视频 | 久久久a久久爽久久v久久| a级毛片在线看网站| 性色avwww在线观看| 国产又爽黄色视频| 国产熟女午夜一区二区三区| 欧美日韩av久久| 捣出白浆h1v1| 亚洲国产欧美在线一区| 欧美日韩亚洲高清精品| 狠狠婷婷综合久久久久久88av| 精品一区在线观看国产| 国产野战对白在线观看| 国产黄色免费在线视频| 亚洲av福利一区| 黑人巨大精品欧美一区二区蜜桃| 国产熟女欧美一区二区| 亚洲内射少妇av| 狠狠精品人妻久久久久久综合| 满18在线观看网站| 啦啦啦在线免费观看视频4| 国产午夜精品一二区理论片| 日韩精品有码人妻一区| 久久久久精品性色| 国产精品av久久久久免费| 人人澡人人妻人| 国产淫语在线视频| 99久久中文字幕三级久久日本| 国产精品国产三级国产专区5o| 高清黄色对白视频在线免费看| 人人妻人人添人人爽欧美一区卜| 美女午夜性视频免费| a级片在线免费高清观看视频| 男人操女人黄网站| 中文乱码字字幕精品一区二区三区| 九草在线视频观看| 亚洲国产av新网站| 18禁动态无遮挡网站| 国产免费福利视频在线观看| 热re99久久国产66热| 久久精品国产亚洲av涩爱| 免费不卡的大黄色大毛片视频在线观看| 最近的中文字幕免费完整| 国产精品av久久久久免费| 成人影院久久| 国产精品久久久久久精品古装| 国产精品 国内视频| 国产日韩一区二区三区精品不卡| 国产精品久久久久久久久免| 91午夜精品亚洲一区二区三区| 嫩草影院入口| 久久久久精品人妻al黑| 免费不卡的大黄色大毛片视频在线观看| 精品国产超薄肉色丝袜足j| 天天躁夜夜躁狠狠久久av| 亚洲国产欧美网| 欧美人与善性xxx| 精品午夜福利在线看| 在现免费观看毛片| 在线精品无人区一区二区三| 亚洲国产av新网站| 三级国产精品片| 国产极品天堂在线| 不卡av一区二区三区| 成人亚洲欧美一区二区av| 人人澡人人妻人| 91成人精品电影| 精品一区二区免费观看| 亚洲欧美精品综合一区二区三区 | 美女中出高潮动态图| 99久久中文字幕三级久久日本| 欧美成人精品欧美一级黄| 久久99蜜桃精品久久| 欧美人与性动交α欧美软件| 看免费成人av毛片| 久久精品久久精品一区二区三区| 国产精品国产av在线观看| 丰满迷人的少妇在线观看| tube8黄色片| 电影成人av| 丝袜脚勾引网站| 色94色欧美一区二区| 美女xxoo啪啪120秒动态图| 少妇被粗大的猛进出69影院| 最近最新中文字幕免费大全7| av在线app专区| 26uuu在线亚洲综合色| 香蕉精品网在线| 国产精品 国内视频| 日日摸夜夜添夜夜爱| 18+在线观看网站| 色吧在线观看| videos熟女内射| 人人妻人人澡人人看| 国产爽快片一区二区三区| 交换朋友夫妻互换小说| 国产免费一区二区三区四区乱码| 欧美日本中文国产一区发布| 人人妻人人添人人爽欧美一区卜| 人妻 亚洲 视频| 成年女人毛片免费观看观看9 | 久久免费观看电影| 人妻少妇偷人精品九色| 大陆偷拍与自拍| 成人影院久久| 深夜精品福利| 中国国产av一级| 人妻一区二区av| 亚洲,一卡二卡三卡| 少妇的逼水好多| 国产日韩一区二区三区精品不卡| 欧美黄色片欧美黄色片| 欧美精品av麻豆av| 黄片小视频在线播放| 精品一区在线观看国产| 国产又色又爽无遮挡免| 色哟哟·www| av视频免费观看在线观看| www.精华液| 赤兔流量卡办理| 亚洲内射少妇av| 午夜福利一区二区在线看| 亚洲国产欧美在线一区| h视频一区二区三区|