• <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ù)對比分析
    船底凹槽氣層流動特征及影響因素研究
    蘇里格氣田致密砂巖氣層識別難點及方法評述
    欧美日韩瑟瑟在线播放| 国产欧美日韩综合在线一区二区| 国产伦人伦偷精品视频| 中文字幕人妻丝袜一区二区| 亚洲av成人不卡在线观看播放网| 亚洲专区国产一区二区| 日韩免费高清中文字幕av| av天堂久久9| 精品久久蜜臀av无| 久99久视频精品免费| 又大又爽又粗| 免费高清视频大片| 51午夜福利影视在线观看| 最好的美女福利视频网| 精品国产亚洲在线| 欧美精品啪啪一区二区三区| 99久久99久久久精品蜜桃| 黄片播放在线免费| 级片在线观看| 伦理电影免费视频| 高清av免费在线| 成人国语在线视频| 18禁黄网站禁片午夜丰满| 99香蕉大伊视频| 欧美乱妇无乱码| 交换朋友夫妻互换小说| 亚洲av日韩精品久久久久久密| 亚洲成人国产一区在线观看| 99re在线观看精品视频| 精品久久久久久,| 日韩精品青青久久久久久| 麻豆成人av在线观看| 久久精品亚洲熟妇少妇任你| 欧美日韩av久久| www.精华液| 久久性视频一级片| 香蕉国产在线看| 桃红色精品国产亚洲av| 国产熟女午夜一区二区三区| 青草久久国产| 久久亚洲真实| 国内毛片毛片毛片毛片毛片| 国产一区二区三区在线臀色熟女 | 亚洲人成77777在线视频| 一区福利在线观看| 在线观看一区二区三区激情| 99国产极品粉嫩在线观看| 亚洲国产欧美日韩在线播放| 亚洲精品国产区一区二| 国产精品一区二区精品视频观看| 欧美日韩国产mv在线观看视频| 最新美女视频免费是黄的| 欧美丝袜亚洲另类 | 男人操女人黄网站| a级毛片在线看网站| 1024视频免费在线观看| 亚洲欧美一区二区三区黑人| 成人影院久久| 亚洲视频免费观看视频| 在线观看免费高清a一片| 午夜福利在线观看吧| 欧美日韩黄片免| 亚洲中文字幕日韩| 又紧又爽又黄一区二区| 亚洲精品国产区一区二| 91字幕亚洲| 亚洲少妇的诱惑av| 亚洲专区中文字幕在线| 亚洲av成人不卡在线观看播放网| 999精品在线视频| 亚洲中文字幕日韩| 日本黄色视频三级网站网址| 女人被躁到高潮嗷嗷叫费观| 国产精品影院久久| av天堂在线播放| 99re在线观看精品视频| 在线永久观看黄色视频| 99在线人妻在线中文字幕| 老汉色av国产亚洲站长工具| 黄色视频不卡| 老汉色av国产亚洲站长工具| 国产免费现黄频在线看| 美女高潮喷水抽搐中文字幕| 亚洲av美国av| 亚洲av成人av| 天堂动漫精品| 成人亚洲精品一区在线观看| 欧美激情高清一区二区三区| 亚洲精品中文字幕一二三四区| 夜夜爽天天搞| 中亚洲国语对白在线视频| 高清毛片免费观看视频网站 | www国产在线视频色| 成人亚洲精品av一区二区 | 久久久久久人人人人人| 手机成人av网站| 精品一区二区三区视频在线观看免费 | 亚洲第一欧美日韩一区二区三区| 国产成人一区二区三区免费视频网站| 亚洲人成电影免费在线| 嫁个100分男人电影在线观看| 视频在线观看一区二区三区| 一个人观看的视频www高清免费观看 | 在线看a的网站| 午夜福利在线观看吧| av网站免费在线观看视频| 色精品久久人妻99蜜桃| 中亚洲国语对白在线视频| 亚洲人成电影观看| 在线观看免费日韩欧美大片| 国产av又大| 波多野结衣一区麻豆| 欧美一级毛片孕妇| av网站在线播放免费| 久久久久久久久免费视频了| 欧美最黄视频在线播放免费 | 欧美精品亚洲一区二区| 欧美精品亚洲一区二区| 日韩精品中文字幕看吧| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品av麻豆狂野| 中文字幕av电影在线播放| 成人精品一区二区免费| 又黄又爽又免费观看的视频| 精品国产超薄肉色丝袜足j| 亚洲少妇的诱惑av| 日韩精品免费视频一区二区三区| 午夜a级毛片| 99在线人妻在线中文字幕| 亚洲一码二码三码区别大吗| 亚洲人成电影免费在线| 一级片'在线观看视频| 国产激情欧美一区二区| 中文字幕精品免费在线观看视频| 一级a爱视频在线免费观看| 午夜精品国产一区二区电影| 免费少妇av软件| 国产一卡二卡三卡精品| 亚洲精品国产精品久久久不卡| 99国产综合亚洲精品| 国产精品久久久久久人妻精品电影| 狠狠狠狠99中文字幕| 国产精品成人在线| 热99re8久久精品国产| 十分钟在线观看高清视频www| 国产精品久久电影中文字幕| 亚洲欧美日韩高清在线视频| netflix在线观看网站| 两性午夜刺激爽爽歪歪视频在线观看 | 成人手机av| 在线观看免费午夜福利视频| 欧美中文日本在线观看视频| 亚洲国产精品999在线| 亚洲人成77777在线视频| 久久狼人影院| 精品一品国产午夜福利视频| 成人免费观看视频高清| 免费在线观看完整版高清| 久久精品人人爽人人爽视色| 长腿黑丝高跟| 91精品三级在线观看| 俄罗斯特黄特色一大片| a在线观看视频网站| 亚洲人成伊人成综合网2020| 夫妻午夜视频| 久久影院123| 不卡av一区二区三区| 亚洲人成网站在线播放欧美日韩| 女人被躁到高潮嗷嗷叫费观| 天天影视国产精品| 国产一卡二卡三卡精品| 黄色视频不卡| 亚洲五月天丁香| 久久人妻熟女aⅴ| 日本免费a在线| 99精国产麻豆久久婷婷| 精品久久久久久久毛片微露脸| a级毛片在线看网站| 久99久视频精品免费| 91麻豆av在线| 精品国产乱子伦一区二区三区| 亚洲 欧美一区二区三区| 91在线观看av| 五月开心婷婷网| 在线天堂中文资源库| 国内毛片毛片毛片毛片毛片| 亚洲人成电影观看| 精品无人区乱码1区二区| 亚洲欧洲精品一区二区精品久久久| 丰满迷人的少妇在线观看| 久久精品亚洲精品国产色婷小说| 国产一区二区三区视频了| 亚洲美女黄片视频| 国产黄a三级三级三级人| 亚洲欧美一区二区三区黑人| 久久香蕉激情| 亚洲男人的天堂狠狠| 一区二区三区国产精品乱码| 国产精品98久久久久久宅男小说| 亚洲精品在线观看二区| 国产一区二区三区综合在线观看| 日韩av在线大香蕉| 夜夜夜夜夜久久久久| 久久国产亚洲av麻豆专区| 亚洲熟女毛片儿| 丁香欧美五月| 国产一区二区三区在线臀色熟女 | 国产亚洲精品综合一区在线观看 | 好男人电影高清在线观看| 日本 av在线| 性色av乱码一区二区三区2| 中文字幕人妻熟女乱码| 成人18禁在线播放| 欧美 亚洲 国产 日韩一| 国产成人欧美| 午夜免费观看网址| 18禁美女被吸乳视频| 很黄的视频免费| 一级毛片精品| 黄网站色视频无遮挡免费观看| 亚洲国产毛片av蜜桃av| 88av欧美| 一个人观看的视频www高清免费观看 | 午夜精品久久久久久毛片777| 大码成人一级视频| 在线观看舔阴道视频| 天天躁夜夜躁狠狠躁躁| 在线免费观看的www视频| 久久久久久久久久久久大奶| 女人被狂操c到高潮| 精品高清国产在线一区| 纯流量卡能插随身wifi吗| 国产精品二区激情视频| 一边摸一边抽搐一进一小说| 可以免费在线观看a视频的电影网站| 亚洲欧美一区二区三区久久| 丝袜在线中文字幕| 亚洲第一av免费看| 国产免费现黄频在线看| a在线观看视频网站| 国产精品 欧美亚洲| 99国产精品一区二区三区| 久久久久亚洲av毛片大全| 久久香蕉国产精品| 两性夫妻黄色片| 国产精品一区二区在线不卡| 日韩国内少妇激情av| 神马国产精品三级电影在线观看 | 日本黄色日本黄色录像| 日韩大码丰满熟妇| 精品久久久久久成人av| 99国产精品一区二区三区| 午夜影院日韩av| 国内毛片毛片毛片毛片毛片| 免费在线观看亚洲国产| 啦啦啦在线免费观看视频4| 欧美日韩亚洲高清精品| 国产黄a三级三级三级人| 大香蕉久久成人网| 国产av一区在线观看免费| 国产精品99久久99久久久不卡| 久久狼人影院| 久久久久久久久久久久大奶| 成人三级做爰电影| 两性午夜刺激爽爽歪歪视频在线观看 | 丰满迷人的少妇在线观看| 国产精品1区2区在线观看.| 一a级毛片在线观看| 在线看a的网站| 最新在线观看一区二区三区| 亚洲午夜理论影院| 国产成人欧美| 天堂√8在线中文| 后天国语完整版免费观看| 国产成人av激情在线播放| 国产精品98久久久久久宅男小说| 一a级毛片在线观看| 久久亚洲真实| 正在播放国产对白刺激| 变态另类成人亚洲欧美熟女 | 亚洲av美国av| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色a级毛片大全视频| 亚洲精品国产区一区二| 亚洲久久久国产精品| 国产极品粉嫩免费观看在线| 亚洲av成人一区二区三| 亚洲自拍偷在线| 国产有黄有色有爽视频| 女人精品久久久久毛片| 韩国av一区二区三区四区| 国产男靠女视频免费网站| 免费女性裸体啪啪无遮挡网站| 首页视频小说图片口味搜索| 欧美乱色亚洲激情| 午夜91福利影院| 国产男靠女视频免费网站| 免费观看人在逋| 91麻豆av在线| 老司机在亚洲福利影院| 日本 av在线| 久久精品成人免费网站| 9色porny在线观看| 12—13女人毛片做爰片一| 美女高潮到喷水免费观看| 女警被强在线播放| 国产亚洲精品综合一区在线观看 | 午夜免费成人在线视频| 视频区欧美日本亚洲| 欧美日韩一级在线毛片| 亚洲 欧美一区二区三区| 精品电影一区二区在线| 露出奶头的视频| 亚洲成a人片在线一区二区| 久久精品亚洲精品国产色婷小说| 亚洲国产毛片av蜜桃av| 国产成人精品无人区| 日韩一卡2卡3卡4卡2021年| 亚洲午夜理论影院| 亚洲国产精品999在线| 99国产精品一区二区蜜桃av| 一级片'在线观看视频| 欧美国产精品va在线观看不卡| 国产av一区二区精品久久| 日本精品一区二区三区蜜桃| 18美女黄网站色大片免费观看| 亚洲精品久久午夜乱码| 999久久久国产精品视频| 国产免费现黄频在线看| 在线av久久热| 国产乱人伦免费视频| 亚洲成人国产一区在线观看| 国产欧美日韩一区二区精品| 亚洲色图综合在线观看| 好男人电影高清在线观看| 黑丝袜美女国产一区| e午夜精品久久久久久久| 久久精品亚洲熟妇少妇任你| 美女午夜性视频免费| 国产精品久久视频播放| 久久人妻av系列| 老汉色∧v一级毛片| 五月开心婷婷网| 精品久久久久久电影网| 又黄又粗又硬又大视频| 80岁老熟妇乱子伦牲交| 嫁个100分男人电影在线观看| av国产精品久久久久影院| av天堂在线播放| 大型av网站在线播放| 久久热在线av| 三上悠亚av全集在线观看| 国产有黄有色有爽视频| 午夜久久久在线观看| 激情视频va一区二区三区| 美女高潮到喷水免费观看| 久久久久久亚洲精品国产蜜桃av| 成年版毛片免费区| 久久精品国产清高在天天线| 波多野结衣一区麻豆| 国产精品乱码一区二三区的特点 | 天天影视国产精品| 香蕉国产在线看| 99久久99久久久精品蜜桃| 精品国产一区二区三区四区第35| 91国产中文字幕| 99国产精品99久久久久| 亚洲伊人色综图| 18禁美女被吸乳视频| 久久久久久久久免费视频了| 亚洲一区二区三区色噜噜 | 热99re8久久精品国产| 色老头精品视频在线观看| 国产一卡二卡三卡精品| 视频区图区小说| 无限看片的www在线观看| 精品一区二区三卡| 久久精品亚洲av国产电影网| 免费观看精品视频网站| 神马国产精品三级电影在线观看 | 亚洲五月婷婷丁香| 国产成人免费无遮挡视频| 国产xxxxx性猛交| 制服诱惑二区| 搡老熟女国产l中国老女人| 久热这里只有精品99| 69精品国产乱码久久久| 精品久久久精品久久久| 亚洲国产看品久久| 欧美精品一区二区免费开放| 国产亚洲精品综合一区在线观看 | 亚洲激情在线av| 无人区码免费观看不卡| 欧美精品啪啪一区二区三区| 欧美老熟妇乱子伦牲交| e午夜精品久久久久久久| 一级a爱视频在线免费观看| 五月开心婷婷网| 国产精品永久免费网站| 老熟妇仑乱视频hdxx| 久久精品国产亚洲av高清一级| 久久国产精品人妻蜜桃| 久久欧美精品欧美久久欧美| 中亚洲国语对白在线视频| 亚洲 欧美一区二区三区| 午夜免费成人在线视频| 麻豆av在线久日| 欧美乱码精品一区二区三区| 国产高清激情床上av| 国产国语露脸激情在线看| 国产亚洲av高清不卡| 亚洲九九香蕉| 国产在线精品亚洲第一网站| 久久人人精品亚洲av| 欧美乱妇无乱码| 欧美国产精品va在线观看不卡| 久久久水蜜桃国产精品网| 欧美日韩亚洲国产一区二区在线观看| 999精品在线视频| 大型av网站在线播放| 国产精品野战在线观看 | 女人被狂操c到高潮| 午夜两性在线视频| 男女之事视频高清在线观看| 国产日韩一区二区三区精品不卡| 亚洲av成人不卡在线观看播放网| 交换朋友夫妻互换小说| 伦理电影免费视频| 嫩草影视91久久| 麻豆av在线久日| 久久国产亚洲av麻豆专区| 国产高清激情床上av| 美女 人体艺术 gogo| 9热在线视频观看99| 国产一区二区激情短视频| 国产亚洲欧美精品永久| 成人av一区二区三区在线看| 亚洲aⅴ乱码一区二区在线播放 | 国产成人精品久久二区二区免费| 天天影视国产精品| 午夜精品久久久久久毛片777| 两性午夜刺激爽爽歪歪视频在线观看 | 不卡一级毛片| 黑丝袜美女国产一区| 天堂中文最新版在线下载| 国产有黄有色有爽视频| 91大片在线观看| 男女午夜视频在线观看| 午夜精品在线福利| 欧美日韩精品网址| 欧美激情高清一区二区三区| 国产成人av激情在线播放| 日韩精品免费视频一区二区三区| bbb黄色大片| 天堂俺去俺来也www色官网| 国产伦人伦偷精品视频| 看免费av毛片| 亚洲精品久久成人aⅴ小说| 在线观看日韩欧美| 热re99久久精品国产66热6| 男女之事视频高清在线观看| 亚洲精品av麻豆狂野| 在线观看免费视频网站a站| 精品久久久久久电影网| 亚洲人成网站在线播放欧美日韩| 国产有黄有色有爽视频| 欧美最黄视频在线播放免费 | 欧美激情高清一区二区三区| 最近最新中文字幕大全电影3 | 国产精品久久久久久人妻精品电影| 丁香六月欧美| 亚洲激情在线av| 国产成人欧美在线观看| 久久香蕉精品热| 亚洲精品美女久久av网站| 亚洲人成伊人成综合网2020| 国产成人免费无遮挡视频| 欧美日韩瑟瑟在线播放| 国产视频一区二区在线看| 国产精品久久久av美女十八| 乱人伦中国视频| 免费在线观看影片大全网站| 国产成年人精品一区二区 | 午夜免费成人在线视频| 国产精品偷伦视频观看了| 亚洲色图av天堂| 叶爱在线成人免费视频播放| 国产1区2区3区精品| 亚洲男人天堂网一区| 午夜成年电影在线免费观看| 99久久人妻综合| 国产精品偷伦视频观看了| 水蜜桃什么品种好| 18禁美女被吸乳视频| 亚洲国产欧美网| 国产三级黄色录像| 一级毛片高清免费大全| 91精品三级在线观看| 国产精品一区二区在线不卡| 久久狼人影院| 欧美黄色淫秽网站| 三上悠亚av全集在线观看| 亚洲色图 男人天堂 中文字幕| 变态另类成人亚洲欧美熟女 | 99国产精品一区二区蜜桃av| 亚洲一区二区三区不卡视频| 成人亚洲精品一区在线观看| 亚洲精品av麻豆狂野| 99国产精品免费福利视频| 高清在线国产一区| 亚洲精品久久成人aⅴ小说| 成人三级做爰电影| 亚洲自拍偷在线| 性少妇av在线| 女人被躁到高潮嗷嗷叫费观| 成年人黄色毛片网站| 在线观看免费视频日本深夜| 1024视频免费在线观看| 国产精品乱码一区二三区的特点 | 欧美人与性动交α欧美软件| 淫秽高清视频在线观看| 国产精品电影一区二区三区| 亚洲熟妇熟女久久| 亚洲情色 制服丝袜| 18禁国产床啪视频网站| 免费av中文字幕在线| 国产精品香港三级国产av潘金莲| 亚洲成国产人片在线观看| 色在线成人网| 最新在线观看一区二区三区| 新久久久久国产一级毛片| 极品人妻少妇av视频| 精品电影一区二区在线| 免费在线观看黄色视频的| 成人特级黄色片久久久久久久| 久久精品成人免费网站| 精品福利观看| 免费搜索国产男女视频| 亚洲一区中文字幕在线| 99国产精品一区二区蜜桃av| 欧美国产精品va在线观看不卡| 性色av乱码一区二区三区2| 欧美乱码精品一区二区三区| 久久中文字幕人妻熟女| 咕卡用的链子| 亚洲九九香蕉| 国产免费男女视频| 少妇裸体淫交视频免费看高清 | 久久伊人香网站| 一本大道久久a久久精品| videosex国产| 91在线观看av| 国产精品自产拍在线观看55亚洲| 久久国产亚洲av麻豆专区| 久久久国产精品麻豆| 午夜91福利影院| 日韩欧美免费精品| 欧美人与性动交α欧美软件| 亚洲人成77777在线视频| 中文字幕高清在线视频| 国产一区二区三区在线臀色熟女 | 美女高潮到喷水免费观看| 可以在线观看毛片的网站| 欧美日本中文国产一区发布| 国产成人免费无遮挡视频| 视频区欧美日本亚洲| 亚洲av成人一区二区三| 亚洲性夜色夜夜综合| 欧美日韩国产mv在线观看视频| 午夜福利免费观看在线| 欧美乱码精品一区二区三区| 俄罗斯特黄特色一大片| 中文字幕精品免费在线观看视频| 精品免费久久久久久久清纯| 美女国产高潮福利片在线看| 夜夜爽天天搞| 国产亚洲av高清不卡| 久久精品亚洲av国产电影网| 亚洲一区二区三区色噜噜 | 国产精品电影一区二区三区| 亚洲中文av在线| 好看av亚洲va欧美ⅴa在| 少妇裸体淫交视频免费看高清 | 国产精品一区二区免费欧美| 亚洲精品一二三| 丰满迷人的少妇在线观看| a级片在线免费高清观看视频| 国产在线精品亚洲第一网站| 亚洲精品国产精品久久久不卡| 欧美乱码精品一区二区三区| 国产一区二区三区在线臀色熟女 | 亚洲中文av在线| 中文字幕高清在线视频| 亚洲九九香蕉| 又黄又爽又免费观看的视频| 最近最新免费中文字幕在线| 黄色成人免费大全| 亚洲五月婷婷丁香| 欧美日韩亚洲高清精品| 色婷婷av一区二区三区视频| 国产又爽黄色视频| 精品久久久久久成人av| 午夜免费观看网址| 怎么达到女性高潮| 久久精品人人爽人人爽视色| 香蕉久久夜色| 一个人观看的视频www高清免费观看 | 久久欧美精品欧美久久欧美| 欧美黄色淫秽网站| 久久狼人影院| 夫妻午夜视频| 午夜福利,免费看| 国产激情欧美一区二区| 久久中文字幕人妻熟女| 欧美日韩av久久|