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

    顧及氣候類型的區(qū)域加權(quán)平均溫度模型

    2023-02-23 07:41:28劉金濤呂偉才肖星星盧???/span>
    導(dǎo)航定位學(xué)報(bào) 2023年6期
    關(guān)鍵詞:探空平均溫度緯度

    劉金濤,呂偉才,肖星星,盧???/p>

    顧及氣候類型的區(qū)域加權(quán)平均溫度模型

    劉金濤,呂偉才,肖星星,盧福康

    (安徽理工大學(xué) 空間信息與測(cè)繪工程學(xué)院/礦區(qū)環(huán)境與災(zāi)害協(xié)同監(jiān)測(cè)煤炭行業(yè)工程研究中心/礦山采動(dòng)災(zāi)害空天地協(xié)同監(jiān)測(cè)與預(yù)警安徽普通高校重點(diǎn)實(shí)驗(yàn)室,安徽 淮南 232001)

    針對(duì)由于中國(guó)土地面積廣闊,氣候類型復(fù)雜,大氣加權(quán)平均溫度(m)易受氣候類型、緯度和海拔影響的問題,提出一種顧及氣候類型的區(qū)域加權(quán)平均溫度模型:根據(jù)中國(guó)區(qū)域內(nèi)分布廣泛的5種氣候類型,利用中國(guó)境內(nèi)73個(gè)探空站2007—2015年的探空數(shù)據(jù),分區(qū)域建立顧及氣候類型及地面氣溫(s)、飽和水汽壓(s)的大氣加權(quán)平均溫度模型ET-TM模型;并用2016年的探空數(shù)據(jù)結(jié)合貝維斯(Bevis)模型和GPT-3模型共同驗(yàn)證ET-TM模型的精度。結(jié)果表明:ET-TM模型整體精度高于Bevis模型和GPT-3模型;分區(qū)域建立大氣加權(quán)平均溫度模型有利于減少緯度和高程對(duì)大氣加權(quán)平均溫度的影響;對(duì)比Bevis模型和GPT-3模型,ET-TM模型更適用于中國(guó)區(qū)域內(nèi)的大氣可降水量反演。

    大氣可降水量(PWV);大氣加權(quán)平均溫度(m);氣候類型;水汽反演;殘差擬合

    0 引言

    全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system,GNSS)是20世紀(jì)對(duì)人類影響最大的空間信息技術(shù)。近年來,利用GNSS導(dǎo)航衛(wèi)星進(jìn)行水汽反演有了巨大的進(jìn)步[1]。大氣加權(quán)平均溫度(atmosphericweighted mean temperature,m)是濕延遲向大氣可降水量(precipitable water vapor,PWV)轉(zhuǎn)換的重要參數(shù),其m精度的高低將直接影響到GNSS水汽反演的精度[2-3]。目前常用的m計(jì)算方法主要有2種[4]:一種是利用探空站的水汽壓和氣溫?cái)?shù)據(jù)積分得到,該方法精度較高[5];另一種是利用經(jīng)驗(yàn)?zāi)P陀?jì)算,如全球壓力和溫度模型(global pressure and temperature model,GPT)GPT-2模型[6]、GPT-3模型[7]等。

    文獻(xiàn)[8]首次提出了m的概念,并認(rèn)為北美地區(qū)的m和地面氣溫(s)具有很強(qiáng)的相關(guān)性。該模型是基于美國(guó)探空站建立的適合中緯度地區(qū)的模型;但由于對(duì)流層的區(qū)域不確定性,該模型應(yīng)用于其他地區(qū)勢(shì)必會(huì)產(chǎn)生較大的偏差。此后,諸多學(xué)者開始分析研究局地m變化規(guī)律。文獻(xiàn)[9]利用2015—2018年的探空數(shù)據(jù)建立了顧及s、飽和水汽壓(s)、海拔、緯度及殘差的BET模型,其精度優(yōu)于貝維斯(Bevis)模型、GPT3w-1模型和GPT3w-5模型。該模型的偏差(ias)和均方根誤差(ms)抗高程和緯度的能力優(yōu)于上述3種模型。文獻(xiàn)[10]利用2007—2014年的全球大地測(cè)量觀測(cè)系統(tǒng)(global geodetic observing system,GGOS)數(shù)據(jù)建立了一種顧及季節(jié)變化以及垂直遞減率的CTm格網(wǎng)模型,且該模型無需氣象因子的輸入。相較于GPT2w模型,該模型在我國(guó)地區(qū)表現(xiàn)出良好的精度和穩(wěn)定性。文獻(xiàn)[11]利用2014—2016年的探空數(shù)據(jù),在Bevis模型的基礎(chǔ)上加入了高程和季節(jié)變化,建立了中國(guó)西部地區(qū)的新m模型,其在中國(guó)西部地區(qū)表現(xiàn)出了良好的適用性,受到高程和四季變化的影響較小。文獻(xiàn)[12]利用2015—2017年的中國(guó)南部地區(qū)的探空數(shù)據(jù)分別建立了顧及s的m-SC1模型以及顧及s、高程、季節(jié)變化和緯度的m-SC2模型,在中國(guó)南部地區(qū)未表現(xiàn)出明顯的季節(jié)性變化,在空間域上有更好的適用性,有效提高了中國(guó)南部地區(qū)的GNSS水汽反演精度。因此利用長(zhǎng)時(shí)間序列的氣象資料建立了局地m模型,在提升局地GNSS-PWV反演方面有著重要意義。

    中國(guó)土地面積廣闊,氣候類型復(fù)雜,地勢(shì)整體表現(xiàn)為西高東低,依據(jù)氣候類型進(jìn)行m建模不僅僅可以降低氣候?qū)的影響,還可相應(yīng)地降低緯度和海拔對(duì)m的影響。因此本文顧及氣候類型對(duì)m的影響,構(gòu)建多因子m模型,并對(duì)其進(jìn)行精度評(píng)定。

    1 數(shù)據(jù)來源及原理

    選擇中國(guó)境內(nèi)的73個(gè)探空站(如圖1所示),從2007—2016共10 a的探空數(shù)據(jù),采樣時(shí)間分別為協(xié)調(diào)世界時(shí)(universal time coordinated,UTC)的0點(diǎn)和12點(diǎn),探空數(shù)據(jù)來自美國(guó)懷俄明大學(xué)。為證明探空數(shù)據(jù)的完整性,從73個(gè)探空站點(diǎn)中隨機(jī)選取8個(gè)站點(diǎn)計(jì)算其數(shù)據(jù)完整性,結(jié)果如表1所示。

    圖1 中國(guó)部分探空站點(diǎn)分布

    表1 部分站點(diǎn)每年的數(shù)據(jù)完整率 %

    2 大氣加權(quán)溫度(Tm)模型的建立

    2.1 相關(guān)性分析

    中國(guó)陸地面積廣袤,地形氣候復(fù)雜,影響m的因素較多。研究表明,大氣加權(quán)平均溫度與地面氣溫[14]和水汽壓[15-16]均具有一定的相關(guān)性。選取全國(guó)73個(gè)探空站利用積分法計(jì)算出相應(yīng)的m。以北京探空站為例,探究m與s、s之間的相關(guān)性(如圖2所示)。

    圖2 Tm與Ts、Es相關(guān)性分析

    從圖2中可以看出:m與s之間具有線性關(guān)系,且隨著s的增加,m也隨之增加,二者呈正相關(guān)關(guān)系,相關(guān)系數(shù)為0.838,相關(guān)性較高;m與s存在冪關(guān)系,相關(guān)系數(shù)為0.855。因此在m建模的過程中,在Bevis模型的基礎(chǔ)上考慮了s因素。

    2.2 Tm模型的建立

    熱帶季風(fēng)氣候、亞熱帶季風(fēng)氣候、溫帶季風(fēng)氣候、溫帶大陸性氣候和高山高原氣候是我國(guó)主要的5種氣候類型。按照各探空站所屬的氣候類型,將73個(gè)探空站進(jìn)行分類[17],分類結(jié)果如表2所示。利用2007—2015年共9 a的探空數(shù)據(jù)進(jìn)行建模,并用2016年的探空數(shù)據(jù)驗(yàn)證所建立模型的精度。

    表2 各氣候區(qū)域內(nèi)探空站數(shù)量

    利用73個(gè)探空站9 a的數(shù)據(jù),采用多元非線性擬合的方法,建立s、s-m多因子的大氣加權(quán)溫度模型。

    設(shè)m函數(shù)方程為

    式中、、、為待擬合參數(shù)。

    擬合出的m模型如表3所示。

    表3 各氣候區(qū)域內(nèi)的多因子Tm模型

    (續(xù))

    由表3可知,5種氣候區(qū)域的m模型中,m-s、s模型計(jì)算的m與探空站積分計(jì)算的m相關(guān)系數(shù)均在0.8以上,擬合效果較好,其中溫帶季風(fēng)氣候區(qū)域內(nèi)m模型相關(guān)系數(shù)達(dá)到了0.92。

    2.3 殘差分析

    研究表明,實(shí)測(cè)獲取的m與m模型之間的殘差具有一定周期性[18],因此提升m模型的精度不能只考慮氣象因素,還應(yīng)當(dāng)考慮到s、s-m模型隨時(shí)間的變化關(guān)系[19]。如圖3所示,利用快速傅里葉變換(fast Fourier transform,F(xiàn)FT)分析m的周期特性,其能量級(jí)數(shù)的峰值出現(xiàn)在181和362 d附近。

    本文使用具有年周期、半年周期的三角函數(shù)來擬合s、s-m模型產(chǎn)生的殘差,將三角函數(shù)的值視為季節(jié)性修正,避免了計(jì)算不同季節(jié)的模型參數(shù),并將最終的模型命名為ET-TM模型。

    圖3 探空站點(diǎn)時(shí)間序列殘差與FFT下的周期特性

    設(shè)m殘差模型為

    式中:oy為年積日;(1)為殘差年均值;(2)、(3)為年周期系數(shù);(4)、(5)為半年周期系數(shù)。

    具體結(jié)果如表4所示。

    表4 各氣候類型下的Tm殘差模型系數(shù)

    最終m模型為

    式中m0為s、s-m模型。在利用該模型計(jì)算m時(shí),只須輸入s、s和oy,即可得到相應(yīng)的大氣加權(quán)平均溫度。

    3 實(shí)驗(yàn)與結(jié)果分析

    為驗(yàn)證該模型計(jì)算m的精度,將2016年的探空數(shù)據(jù)積分得到的m作為真實(shí)值,驗(yàn)證該模型的精度。利用平均偏差ias和均方根誤差ms作為精度評(píng)價(jià)指標(biāo),有:

    3.1 基于探空數(shù)據(jù)的模型精度分析

    利用73個(gè)探空站2016年的探空數(shù)據(jù)來驗(yàn)證ET-TM模型的精度,并結(jié)合Bevis模型和GPT-3模型進(jìn)行對(duì)比,綜合分析ET-TM模型計(jì)算m的性能,并計(jì)算各個(gè)模型的ias和ms。具體結(jié)果如表5所示。

    由表5可知:除熱帶季風(fēng)氣候區(qū)和亞熱帶季風(fēng)氣候區(qū)外,Bevis模型年均ias均表現(xiàn)為負(fù)偏差,其年均ias分別為1.28、-1.26、-2.32、-5.64、1.22 K;GPT-3模型在5種氣候區(qū)域內(nèi)均表現(xiàn)為正偏差,年平均偏差分別為0.90、0.38、1.54、3.65、0.20 K。因此,Bevis模型和GPT-3模型計(jì)算的m在中國(guó)區(qū)域內(nèi)存在一定的系統(tǒng)誤差,同時(shí)Bevis模型和GPT-3模型在高海拔地區(qū)的高山高原氣候區(qū)偏差最大。ET-TM模型在中國(guó)5種氣候區(qū)域的偏差最小為-2.16 K,最大為2.19 K,年平均偏差分別為-0.04、0.03、0.11、0.13、0.35 K,波動(dòng)較小。Bevis模型的ms在5種氣候區(qū)域內(nèi)分別為3.12、4.32、4.81、6.80、2.98 K;GPT-3模型在5種氣候區(qū)域內(nèi)的ms分別為3.82、4.77、4.99、5.07、2.70 K;ET-TM模型的ms分別為2.67、3.33、3.44、2.68、2.26 K,相對(duì)于Bevis模型分別提高了0.45 K(14.4%)、0.99 K(22.9%)、1.37 K(28.5%)、4.12 K(60.6%)、0.72 K(24.16%),相對(duì)于GPT-3分別提高了1.15 K (30.1%)、1.44 K(30.2%)、1.55 K(30.1%)、2.39 K (47.1%)、0.44 K(16.30%)。由此可以說明,ET-TM模型在中國(guó)5種氣候區(qū)域內(nèi)的精度均優(yōu)于Bevis模型和GPT-3模型,且涵蓋因素更多、偏差更小。如圖4所示為Bevis模型、GPT-3模型、ET-TM模型的年均ias、ms分布情況,圖中黑點(diǎn)的面積代表其數(shù)值大小。

    表5 4種氣候區(qū)域內(nèi)3種模型精度對(duì)比 K

    由圖4可知:在熱帶季風(fēng)區(qū)域,雖然探空站點(diǎn)較少,但仍能看出ET-TM模型要優(yōu)于Bevis模型和GPT-3模型;Bevis模型在亞熱帶季風(fēng)區(qū)域和溫帶季風(fēng)區(qū)域的南部適用性較好,年均ias在5種氣候區(qū)域內(nèi)較小,在海拔較高的高山高原氣候區(qū)域內(nèi)系統(tǒng)誤差最大,同時(shí)在溫帶季風(fēng)氣候區(qū)的北部年均ias較大,由此可見Bevis模型易受到海拔和緯度的影響;GPT-3模型在亞熱帶季風(fēng)、溫帶季風(fēng)氣候區(qū)域內(nèi)年均ias較小,在高山高原氣候區(qū)域內(nèi)偏差較大,同樣受到緯度和海拔的影響,且海拔造成的影響明顯大于緯度因素;ET-TM模型考慮到氣候類型對(duì)m的影響,按照氣候類型對(duì)中國(guó)區(qū)域進(jìn)行m建模,這既減小了氣候的影響,又減小了緯度和海拔2種因素的影響,其年均ias在3種模型中最小。從年均ms來看,Bevis模型為單因子模型,僅使用s進(jìn)行了建模,在中國(guó)亞熱帶季風(fēng)區(qū)域內(nèi)年均ms較小,但在高海拔地區(qū)青藏高原區(qū)域,年均ms較大,最大為7.62 K,在溫帶季風(fēng)區(qū)域北端,其年均偏差ms部分超過5 K;GPT-3模型在建模時(shí)考慮了殘差的季節(jié)性變化,其年均ms整體小于Bevis模型,年均ms較大值也出現(xiàn)在高海拔和高緯度地區(qū);ET-TM模型年均ms在3種模型中最小,與Bevis模型相比,最大優(yōu)勢(shì)在于青藏高原地區(qū),ET-TM模型年均ms在青藏高原地區(qū)最大值為2.68 K,與GPT-3模型相比,ET-TM模型年均ms相對(duì)穩(wěn)定,波動(dòng)較小??傮w來說,ET-TM模型要優(yōu)于Bevis模型、GPT-3模型。

    為了更好地檢驗(yàn)ET-TM模型的精度及適用性,繪制熱帶季風(fēng)氣候、亞熱帶季風(fēng)氣候、溫帶季風(fēng)氣候、溫帶大陸氣候和高山高原氣候區(qū)域內(nèi)的不同模型的日均ias和日均ms散點(diǎn)圖并進(jìn)行統(tǒng)計(jì)分析,結(jié)果如圖5、圖6所示。

    由圖5可以看出:ET-TM模型ias在5種氣候類型區(qū)域內(nèi)變化較為平緩,相對(duì)穩(wěn)定,ias總體較小;Bevis模型和GPT-3模型在我國(guó)亞熱帶季風(fēng)區(qū)域內(nèi)適應(yīng)性較高,但其在建模過程當(dāng)中未顧及高程因素,造成在青藏高原等高海拔地區(qū)產(chǎn)生了較大的偏差。Bevis模型在溫帶季風(fēng)、溫帶大陸和高山高原氣候區(qū)內(nèi)日均ias表現(xiàn)為春夏季節(jié)偏差小,秋冬季節(jié)偏差大的特點(diǎn),存在一定的季節(jié)性變化,在3個(gè)氣候區(qū)域內(nèi)整體表現(xiàn)為負(fù)偏差,而在熱帶季風(fēng)氣候區(qū)域內(nèi),其偏差在春冬季節(jié)偏差較大,其他時(shí)間偏差相對(duì)較??;GPT-3和ET-TM模型均考慮了殘差的季節(jié)性變化,因此GPT-3模型和ET-TM模型受季節(jié)的影響較小,二者日均ias變化趨勢(shì)基本相同,但GPT-3模型計(jì)算得出的m整體上高于積分法計(jì)算的m,在5種氣候區(qū)域內(nèi)均表現(xiàn)為正偏差。

    由圖6可知:從ms角度來說,在亞熱帶季風(fēng)區(qū)域內(nèi)3種模型均未表現(xiàn)出明顯的季節(jié)性變化,但在溫帶季風(fēng)區(qū)域內(nèi),3種模型的日均ms在夏季偏小,其余時(shí)間段較大;在熱帶季風(fēng)氣候、溫帶大陸性氣候和高山高原氣候區(qū)域內(nèi),Bevis模型表現(xiàn)出一定的季節(jié)性差異,在熱帶季風(fēng)氣候區(qū)域和高山高原氣候區(qū)域內(nèi)較為明顯,全年表現(xiàn)為夏季ms較小,春秋季節(jié)ms較大;GPT-3模型的日均ms在熱帶季風(fēng)區(qū)域和溫帶季風(fēng)區(qū)域顯示出一定的季節(jié)性變化,在其他3個(gè)區(qū)域均未出現(xiàn)明顯的季節(jié)性變化,但GPT-3模型的日均ms波動(dòng)較大且不穩(wěn)定;ET-TM模型日均ms均小于Bevis模型和GPT-3模型,波動(dòng)范圍小,相對(duì)穩(wěn)定。綜合圖5和圖6可知:從亞熱帶季風(fēng)氣候到溫帶季風(fēng)區(qū)域,隨著緯度的增加,Bevis模型和GPT-3模型的ias和ms逐漸增大;同時(shí),由亞熱帶季風(fēng)到高山高原氣候區(qū),隨著海拔的增加,Bevis模型和GPT-3模型的ias和ms也逐漸變大。由此可見,緯度和高程對(duì)m模型有著極大的影響,但海拔因素的影響要遠(yuǎn)大于緯度因素的影響。ET-TM模型雖未考慮高程和緯度因素的影響,但按照氣候類型劃分區(qū)域建模的方法總體上減小了高程和緯度的變化。因此ET-TM模型在全年當(dāng)中受到季節(jié)性影響較小,全年相對(duì)平穩(wěn),同時(shí)受緯度和高程的影響也相對(duì)較小。

    3.2 Tm對(duì)GNSS-PWV反演的精度分析

    m模型構(gòu)建的目的在于應(yīng)用于GNSS-PWV反演,但在多數(shù)情況下,GNSS觀測(cè)站與探空站并不共址,因此很難全面可靠地研究m模型對(duì)GNSS-PWV反演的影響。目前已有文獻(xiàn)從理論上研究了m對(duì)GNSS-PWV反演的影響[20],本文利用該方法評(píng)估ET-TM模型對(duì)GNSS-PWV反演的影響,即

    式中:p為大氣可降水量PWV的ms;為該區(qū)域的大氣可降水量PWV;π為水汽轉(zhuǎn)換系數(shù)的ms;Tm為大氣加權(quán)平均溫度的年均ms。利用p/作為表征m對(duì)GNSS-PWV反演的影響。

    如表6所示為3種模型計(jì)算PWV的理論ms和相對(duì)誤差。

    表6 ET-TM模型、Bevis模型和GPT-3模型計(jì)算PWV的理論Rms和相對(duì)誤差 mm

    (續(xù))

    依據(jù)表6可知,ET-TM模型的PWV在5種氣候區(qū)域內(nèi)的最大值、最小值分別為0.423、0.086 mm,小于Bevis模型和GPT-3模型,整體上優(yōu)于Bevis模型和GPT-3模型;從PWV/的結(jié)果來看,ET-TM模型在5個(gè)區(qū)域內(nèi)的平均值分別為0.009、0.012、0.013、0.01、0.008 mm,相比于Bevis模型和GPT3模型波動(dòng)更小,展現(xiàn)了ET-TM模型的穩(wěn)定性。因此該模型可以為中國(guó)區(qū)域的GNSS-PWV反演提供高精度的m值。

    4 結(jié)束語

    1)本文利用2007—2015年共9 a的探空數(shù)據(jù),采用回歸分析方法,建立了適用于中國(guó)熱帶季風(fēng)氣候、亞熱帶季風(fēng)氣候、溫帶季風(fēng)氣候、溫帶大陸性氣候和高山高原氣候5種氣候區(qū)域下的ET-TM模型;并使用2016年的數(shù)據(jù)對(duì)建立的模型進(jìn)行了精度分析,得到5種氣候類型下的年均ias和年均ms。

    2)ET-TM模型考慮到氣候因素,同時(shí)利用與m相關(guān)性較高的s和s2個(gè)因子對(duì)m進(jìn)行建模,并對(duì)殘差序列進(jìn)行分析。利用該方法進(jìn)行建模,不僅減小了氣候條件和季節(jié)的影響,還在一定程度上減少了緯度和海拔2種因素的影響。

    3)通過對(duì)PWV、PWV/的分析,ET-TM模型的PWV在5種氣候區(qū)域內(nèi)最大為0.325 mm,最小為0.086 mm,PWV/在5種氣候區(qū)域內(nèi)最大為0.013,最小為0.008,均優(yōu)于Bevis模型和GPT-3模型。

    總體來看,ET-TM模型計(jì)算的m精度更高且更可靠,可以更好地滿足中國(guó)區(qū)域內(nèi)的GNSS-PWV反演需要,更好地服務(wù)于中國(guó)GNSS氣象學(xué)的研究工作。

    [1] 姚歡歡, 黨亞民, 許長(zhǎng)輝, 等. 不同氣候類型北斗反演大氣可降水量性能分析[J]. 測(cè)繪科學(xué), 2022, 47(10): 21-28.

    [2] 黃良珂, 吳丕團(tuán), 王浩宇, 等. 中國(guó)西南地區(qū)GPS大氣水汽轉(zhuǎn)換系數(shù)模型精化研究[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2019, 39(3): 256-261.

    [3] 侯曉玲, 熊永良. 川渝地區(qū)大氣加權(quán)平均溫度模型優(yōu)化[J]. 導(dǎo)航定位學(xué)報(bào), 2023, 11(4): 63-69.

    [4] 謝劭峰, 張繼洪, 王義杰, 等. 大氣加權(quán)平均溫度垂直遞減率時(shí)空特性分析[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2023, 43(6): 593-599.

    [5] 鄭磊. 顧及季節(jié)變化的西南地區(qū)大氣加權(quán)平均溫度模型[J]. 導(dǎo)航定位學(xué)報(bào), 2021, 9(4): 98-103.

    [6] LAGLER K, SCHINDELEGGER M, B?HM J, et al. Gpt2: Empirical slant delay model for radio space geodetic techniques[J]. Geophysical Research Letters, 2013, 40(6): 1069-1073.

    [7] Landskron Daniel, Bohm Johannes. Vmf3/gpt3: Refined discrete and empirical troposphere mapping functions[J]. Journal of Geodesy, 2018, 92(4): 349-360.

    [8] BEVIS M, BUSINGER S, HERRING T A, et al. GPS meteorology: Mapping zenith wet delays onto precipitable water[J]. Journal of Applied Meteorology, 1994, 33(3): 379-386.

    [9] 李浩杰, 劉立龍, 黃良珂, 等. 顧及多因子影響的中國(guó)區(qū)域T_m模型精化研究[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2022, 42(4): 393-397.

    [10] 黃良珂, 彭華, 劉立龍, 等. 顧及垂直遞減率函數(shù)的中國(guó)區(qū)域大氣加權(quán)平均溫度模型[J]. 測(cè)繪學(xué)報(bào), 2020, 49(4): 432-442.

    [11] 莫智翔, 黎杏, 黃良珂, 等. 顧及多因子影響的中國(guó)西部地區(qū)大氣加權(quán)平均溫度模型精化研究[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2021, 41(2): 145-151.

    [12] 廖發(fā)圣, 黃良珂, 劉立龍, 等. 中國(guó)南部地區(qū)大氣加權(quán)平均溫度模型精化研究[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2022, 42(1): 41-47.

    [13] 羅宇, 羅林艷, 呂冠儒. 加權(quán)平均溫度模型對(duì)GPS水汽反演的影響[J]. 測(cè)繪科學(xué), 2018, 43(9): 6-9, 15.

    [14] Ding M. A neural network model for predicting weighted mean temperature[J]. Journal of Geodesy, 2018, 92(10): 1187-1198.

    [15] 王皓, 趙興旺. 香港地區(qū)大氣加權(quán)平均溫度建模與研究[J]. 全球定位系統(tǒng), 2019, 44(1): 119-124.

    [16] 王瀚弘, 魏冠軍, 張幸, 等. 蘭州地區(qū)區(qū)域加權(quán)平均溫度模型構(gòu)建方法研究[J]. 導(dǎo)航定位學(xué)報(bào), 2022, 10(1): 90-96, 129.

    [17] 中華人民共和國(guó)氣候圖集編委會(huì). 中華人民共和國(guó)氣候圖集[M]. 第1版. 北京: 氣象出版社, 2002.

    [18] Yao Y, Zhang B, Xu C, et al. Improved one/multi-parameter models that consider seasonal and geographic variations for estimating weighted mean temperature in ground-based gps meteorology[J]. Journal of Geodesy, 2014, 88(3): 273-282.

    [19] 張迪, 袁林果, 黃良珂, 等. 澳大利亞區(qū)域大氣加權(quán)平均溫度建模[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2022, 47(7): 1146-1153.

    [20] Huang Liangke, Jiang Weiping, Liu Lilong, et al. A new global grid model for the determination of atmospheric weighted mean temperature in gps precipitable water vapor[J]. Journal of Geodesy, 2019, 93(2): 159-176.

    Regional weighted average temperature model considering climate types

    LIU Jintao, LYU Weicai, XIAO Xingxing, LU Fukang

    (School of Geomatics/Coal Industry Engineering Research Center of Mining Area Environmental and Disaster Cooperative Monitoring/Key Laboratory of Aviation-Aerospace-Ground Cooperative Monitoring and Early Warning of Coal Mining-Induced Disasters of Anhui Higher Education Institutes, Anhui University of Science and Technology, Huainan, AnHui 232001, China)

    Aiming at the problem that the atmospheric weighted mean temperature (m) is susceptible to the influence of climate type, latitude and altitude due to the vast land area and complex climate types in China, the paper proposed a regional weighted average temperature model considering climate types: according to the five climate types widely distributed in China, based on the sounding data of 73 sounding stations in China from the year of 2007 to 2015, an atmospheric weighted average temperature model ET-TM model considering climate types, surface air temperature (s) and saturation water vapor pressure (s) was established in different regions; and the accuracy of the ET-TM model was verified by the sounding data of 2016 combined with Bevis model and GPT-3 model. Results showed that the overall accuracy of ET-TM model would be higher than that of Bevis model and GPT-3 model; and the establishment of atmospheric weighted average temperature model by region would be helpful to reduce the influence of latitude and elevation on atmospheric weighted average temperature; furthermore, compared with Bevis model and GPT-3 model, the ET-TM model could be more suitable for the retrieval of atmospheric precipitable water in China.

    precipitable water vapor (PWV); atmospheric weighted mean temperature (m); climate type; water vapor retrieval; residual fitting

    劉金濤,呂偉才,肖星星,等. 顧及氣候類型的區(qū)域加權(quán)平均溫度模型[J]. 導(dǎo)航定位學(xué)報(bào), 2023, 11(6): 119-128.(LIU Jintao, LYU Weicai, XIAO Xingxing, et al. Regional weighted average temperature model considering climate types[J]. Journal of Navigation and Positioning, 2023, 11(6): 119-128.)DOI:10.16547/j.cnki.10-1096.20230615.

    P228

    A

    2095-4999(2023)06-0119-10

    2023-03-02

    2021年度安徽省重大科技專項(xiàng)(202103a05020026);2021年度安徽省重點(diǎn)研究與開發(fā)計(jì)劃項(xiàng)目(202104a 07020014);安徽省自然科學(xué)基金項(xiàng)目(2008085MD114)。

    劉金濤(1998—),男,山東青島人,碩士研究生,研究方向?yàn)镚NSS數(shù)據(jù)處理。

    呂偉才(1965—),女,山東青島人,博士,教授,研究方向?yàn)楝F(xiàn)代測(cè)量數(shù)據(jù)處理、GNSS測(cè)量與數(shù)據(jù)處理。

    猜你喜歡
    探空平均溫度緯度
    探空氣球?yàn)楦呖諝庀蟆鞍衙}”
    3月熱浪來襲悉尼或迎165年以來新紀(jì)錄
    南方地區(qū)圓拱形和鋸齒形大棚內(nèi)溫度四季差別探究*
    云南保山氣溫變化特征及其均生函數(shù)預(yù)測(cè)實(shí)驗(yàn)
    緯度
    齊魯周刊(2017年29期)2017-08-08 06:28:15
    TK-2GPS人影火箭探空數(shù)據(jù)與L波段探空數(shù)據(jù)對(duì)比分析
    徐州地區(qū)加權(quán)平均溫度模型研究
    基于時(shí)空緯度的國(guó)內(nèi)農(nóng)民工創(chuàng)業(yè)研究
    淺談凈舉力對(duì)探空氣球升速及施放高度的影響
    河南科技(2014年14期)2014-02-27 14:12:08
    鄭州探空數(shù)據(jù)庫(kù)的建設(shè)簡(jiǎn)介
    河南科技(2014年4期)2014-02-27 14:06:59
    欧美激情在线99| 成人一区二区视频在线观看| 日本与韩国留学比较| 亚洲欧洲国产日韩| 欧美日韩一区二区视频在线观看视频在线 | 成年人午夜在线观看视频 | 国产亚洲av片在线观看秒播厂 | 国产永久视频网站| 久久久久网色| 日韩强制内射视频| 男女国产视频网站| 欧美xxⅹ黑人| 日本一本二区三区精品| 一本久久精品| 美女黄网站色视频| 一区二区三区高清视频在线| 免费黄频网站在线观看国产| xxx大片免费视频| 人体艺术视频欧美日本| 日韩一区二区视频免费看| 水蜜桃什么品种好| 国国产精品蜜臀av免费| 精品少妇黑人巨大在线播放| 91精品国产九色| 免费观看性生交大片5| 国产午夜精品久久久久久一区二区三区| 一个人看视频在线观看www免费| 欧美成人午夜免费资源| 搡老妇女老女人老熟妇| 亚洲精品,欧美精品| 欧美97在线视频| 男人爽女人下面视频在线观看| 床上黄色一级片| 国内精品一区二区在线观看| 国产一级毛片在线| 禁无遮挡网站| 亚洲国产欧美在线一区| 日韩av在线免费看完整版不卡| 精品亚洲乱码少妇综合久久| 男人爽女人下面视频在线观看| av在线观看视频网站免费| 国产免费又黄又爽又色| 少妇裸体淫交视频免费看高清| 午夜激情福利司机影院| 五月伊人婷婷丁香| 免费大片黄手机在线观看| 熟妇人妻不卡中文字幕| 一本久久精品| 久久综合国产亚洲精品| 国产不卡一卡二| av天堂中文字幕网| 丝袜喷水一区| 成人鲁丝片一二三区免费| 听说在线观看完整版免费高清| 美女国产视频在线观看| 亚洲怡红院男人天堂| 中文字幕久久专区| 久久国内精品自在自线图片| 日韩成人av中文字幕在线观看| 国产精品久久久久久精品电影| .国产精品久久| av专区在线播放| 日韩大片免费观看网站| 又黄又爽又刺激的免费视频.| 国模一区二区三区四区视频| 又粗又硬又长又爽又黄的视频| 白带黄色成豆腐渣| 国产综合懂色| 国产在视频线精品| 国产亚洲91精品色在线| 欧美三级亚洲精品| 久久久久久久久中文| 国产色爽女视频免费观看| 免费观看的影片在线观看| 亚洲熟妇中文字幕五十中出| 成人毛片a级毛片在线播放| 亚洲欧美成人精品一区二区| 欧美日韩亚洲高清精品| 乱人视频在线观看| 日韩不卡一区二区三区视频在线| 日本-黄色视频高清免费观看| 欧美人与善性xxx| 亚洲精品一区蜜桃| 亚洲欧美中文字幕日韩二区| 一区二区三区免费毛片| 国产精品国产三级国产av玫瑰| 日本-黄色视频高清免费观看| 人妻一区二区av| 中文字幕久久专区| 青春草亚洲视频在线观看| 国产精品人妻久久久久久| 老师上课跳d突然被开到最大视频| 国产在视频线精品| 午夜福利视频精品| 一级a做视频免费观看| 久久久成人免费电影| 成人午夜精彩视频在线观看| 亚洲精品一区蜜桃| 亚洲人与动物交配视频| 久久午夜福利片| 国产精品女同一区二区软件| 欧美日韩精品成人综合77777| 欧美一级a爱片免费观看看| 亚洲欧美中文字幕日韩二区| 黄片wwwwww| 亚洲av国产av综合av卡| 免费无遮挡裸体视频| a级一级毛片免费在线观看| 精品99又大又爽又粗少妇毛片| 最近2019中文字幕mv第一页| 免费大片18禁| av免费在线看不卡| 欧美人与善性xxx| 国产一区有黄有色的免费视频 | 国产黄频视频在线观看| 亚洲成人av在线免费| 久久久欧美国产精品| 乱码一卡2卡4卡精品| 国产精品伦人一区二区| 日日啪夜夜爽| 一二三四中文在线观看免费高清| 国产午夜精品久久久久久一区二区三区| 午夜福利网站1000一区二区三区| 三级经典国产精品| av专区在线播放| 日本三级黄在线观看| 尤物成人国产欧美一区二区三区| 天堂网av新在线| 国产精品不卡视频一区二区| 午夜视频国产福利| 成年人午夜在线观看视频 | 精品人妻偷拍中文字幕| 欧美激情久久久久久爽电影| 秋霞伦理黄片| 午夜激情久久久久久久| 大香蕉久久网| 看免费成人av毛片| kizo精华| 免费少妇av软件| 国产欧美日韩精品一区二区| 伊人久久精品亚洲午夜| 波多野结衣巨乳人妻| 国产成人午夜福利电影在线观看| 老女人水多毛片| 久久久久久久久久成人| 欧美高清性xxxxhd video| 91久久精品国产一区二区成人| 国产伦精品一区二区三区四那| 国产综合懂色| 免费人成在线观看视频色| 亚洲在线自拍视频| 亚洲国产精品国产精品| 国产精品一区www在线观看| 日日摸夜夜添夜夜爱| 国产精品爽爽va在线观看网站| 国产成人精品婷婷| 少妇被粗大猛烈的视频| 日韩人妻高清精品专区| 神马国产精品三级电影在线观看| 国产精品人妻久久久影院| 国产高清有码在线观看视频| 少妇丰满av| 国产精品一区二区在线观看99 | 精品一区二区三区人妻视频| 亚洲精品一二三| 欧美日本视频| 九草在线视频观看| 国产探花极品一区二区| 久久精品国产亚洲av涩爱| 亚洲色图av天堂| 可以在线观看毛片的网站| 久久久久久久久久久丰满| 内地一区二区视频在线| 久久精品久久久久久久性| 女人被狂操c到高潮| 嘟嘟电影网在线观看| 中文在线观看免费www的网站| 国产精品久久视频播放| 美女黄网站色视频| 亚洲成人中文字幕在线播放| 亚洲最大成人av| 成人午夜高清在线视频| 大香蕉97超碰在线| 亚洲va在线va天堂va国产| 久久久色成人| 国产人妻一区二区三区在| 少妇人妻精品综合一区二区| 日韩一区二区三区影片| 亚洲精品国产av蜜桃| 麻豆精品久久久久久蜜桃| 亚洲精品成人av观看孕妇| 久久精品夜色国产| 特大巨黑吊av在线直播| 日韩欧美国产在线观看| 国产淫语在线视频| 18+在线观看网站| 欧美潮喷喷水| 亚洲av成人精品一区久久| 99久久人妻综合| 亚洲av不卡在线观看| 国产精品人妻久久久久久| 国产伦在线观看视频一区| 80岁老熟妇乱子伦牲交| 久久久久久久亚洲中文字幕| 免费黄频网站在线观看国产| 天堂网av新在线| 午夜免费男女啪啪视频观看| 亚洲不卡免费看| av卡一久久| 精品久久久久久久久亚洲| 亚洲色图av天堂| 99热这里只有精品一区| 最后的刺客免费高清国语| 亚洲av免费在线观看| 身体一侧抽搐| 干丝袜人妻中文字幕| 91在线精品国自产拍蜜月| 日韩精品有码人妻一区| 亚洲精品一二三| 特大巨黑吊av在线直播| av免费在线看不卡| 亚洲精品日韩在线中文字幕| 国产午夜精品一二区理论片| 嘟嘟电影网在线观看| 最近手机中文字幕大全| 国产国拍精品亚洲av在线观看| 男的添女的下面高潮视频| 久久久欧美国产精品| 日本三级黄在线观看| 观看免费一级毛片| 欧美一区二区亚洲| 日本-黄色视频高清免费观看| 久久亚洲国产成人精品v| 人人妻人人看人人澡| 国产伦一二天堂av在线观看| 欧美高清成人免费视频www| a级一级毛片免费在线观看| 亚洲精华国产精华液的使用体验| 又黄又爽又刺激的免费视频.| 中文字幕制服av| 99久国产av精品| 国产欧美日韩精品一区二区| 久久久久久久久中文| 国产欧美日韩精品一区二区| 国产伦理片在线播放av一区| 天堂av国产一区二区熟女人妻| 日韩欧美一区视频在线观看 | 久久久久久国产a免费观看| 久久久精品94久久精品| 精品久久国产蜜桃| 久久人人爽人人爽人人片va| 美女大奶头视频| 久久人人爽人人爽人人片va| 中国国产av一级| 日日摸夜夜添夜夜爱| 日本一二三区视频观看| 激情五月婷婷亚洲| 亚洲国产欧美人成| 亚洲综合色惰| 搡老妇女老女人老熟妇| 亚洲高清免费不卡视频| 毛片女人毛片| 亚洲最大成人手机在线| 欧美日韩视频高清一区二区三区二| 亚洲精品影视一区二区三区av| 国产成人a∨麻豆精品| 欧美成人a在线观看| 一级黄片播放器| 免费看av在线观看网站| 麻豆精品久久久久久蜜桃| 亚洲av不卡在线观看| 精品久久久噜噜| 高清毛片免费看| 成人无遮挡网站| 国产精品国产三级国产专区5o| 91狼人影院| 亚洲欧美一区二区三区黑人 | 国产精品嫩草影院av在线观看| 久久久久网色| 人妻少妇偷人精品九色| 熟妇人妻不卡中文字幕| 午夜福利视频1000在线观看| 禁无遮挡网站| 2021少妇久久久久久久久久久| 国产人妻一区二区三区在| 精品欧美国产一区二区三| 国产亚洲一区二区精品| 亚洲四区av| 最近视频中文字幕2019在线8| 国产成人精品婷婷| 纵有疾风起免费观看全集完整版 | 又粗又硬又长又爽又黄的视频| 噜噜噜噜噜久久久久久91| 国产高清国产精品国产三级 | 国产亚洲av嫩草精品影院| 国产精品久久久久久精品电影小说 | 亚洲美女视频黄频| 搡女人真爽免费视频火全软件| 一边亲一边摸免费视频| 亚洲婷婷狠狠爱综合网| 日韩三级伦理在线观看| 一边亲一边摸免费视频| 久久6这里有精品| 乱系列少妇在线播放| 免费看日本二区| 午夜精品一区二区三区免费看| 高清视频免费观看一区二区 | 久久精品综合一区二区三区| 久久午夜福利片| 99久久人妻综合| 建设人人有责人人尽责人人享有的 | 内地一区二区视频在线| 高清欧美精品videossex| 国产成年人精品一区二区| 久久精品人妻少妇| av福利片在线观看| 国产成人精品久久久久久| 国产永久视频网站| 久久精品人妻少妇| 我的老师免费观看完整版| 国产极品天堂在线| 日韩欧美三级三区| 天堂√8在线中文| 亚洲欧洲日产国产| av免费在线看不卡| 亚洲精品乱码久久久久久按摩| 婷婷色麻豆天堂久久| 国产有黄有色有爽视频| 精品少妇黑人巨大在线播放| 中国美白少妇内射xxxbb| 亚洲图色成人| 99久国产av精品| 免费看日本二区| 免费在线观看成人毛片| 最近手机中文字幕大全| 99久久人妻综合| 麻豆乱淫一区二区| 免费观看性生交大片5| 色哟哟·www| 国产精品一区二区三区四区久久| 精品久久久久久久久亚洲| 不卡视频在线观看欧美| 国产一区二区三区综合在线观看 | 熟妇人妻不卡中文字幕| 中文字幕亚洲精品专区| av专区在线播放| 日本-黄色视频高清免费观看| 亚洲国产av新网站| 国产视频首页在线观看| 国模一区二区三区四区视频| 国产人妻一区二区三区在| 亚洲国产精品国产精品| 中文字幕亚洲精品专区| 久久久精品免费免费高清| 成人亚洲精品一区在线观看 | 久久精品久久久久久噜噜老黄| 高清日韩中文字幕在线| 黄片wwwwww| 亚洲av一区综合| 国产亚洲最大av| 毛片一级片免费看久久久久| 日本黄色片子视频| 久久人人爽人人片av| 国产色爽女视频免费观看| 久久久久免费精品人妻一区二区| 久久久久久久久久久免费av| 日本与韩国留学比较| 国产 亚洲一区二区三区 | 国产成年人精品一区二区| 嫩草影院精品99| 国产人妻一区二区三区在| 午夜福利在线观看吧| 少妇裸体淫交视频免费看高清| 免费看美女性在线毛片视频| 2018国产大陆天天弄谢| 日韩欧美三级三区| 亚洲欧洲国产日韩| 97超视频在线观看视频| 少妇的逼好多水| 日韩欧美 国产精品| 欧美+日韩+精品| 久久久久久久久久人人人人人人| 最近中文字幕2019免费版| 亚洲精华国产精华液的使用体验| 小蜜桃在线观看免费完整版高清| 国产精品人妻久久久久久| 搞女人的毛片| 久久久久久久久久人人人人人人| 国产av码专区亚洲av| 成人毛片a级毛片在线播放| 亚洲怡红院男人天堂| 69人妻影院| 国产精品一二三区在线看| 亚洲美女搞黄在线观看| 午夜免费观看性视频| 欧美精品一区二区大全| 在线免费观看的www视频| 免费观看无遮挡的男女| 久久久a久久爽久久v久久| 亚洲欧美精品专区久久| 午夜免费激情av| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 老司机影院成人| 亚洲三级黄色毛片| 国产精品福利在线免费观看| 青青草视频在线视频观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲av中文av极速乱| 亚洲成人中文字幕在线播放| 精品一区在线观看国产| 禁无遮挡网站| 最近的中文字幕免费完整| 精品久久久精品久久久| 免费看日本二区| 国产探花极品一区二区| 国产成人免费观看mmmm| 女人十人毛片免费观看3o分钟| 干丝袜人妻中文字幕| 成人午夜精彩视频在线观看| 日本黄大片高清| 偷拍熟女少妇极品色| 老师上课跳d突然被开到最大视频| av在线观看视频网站免费| 欧美xxxx黑人xx丫x性爽| 97超碰精品成人国产| 26uuu在线亚洲综合色| 国产不卡一卡二| 简卡轻食公司| 免费看光身美女| 国产精品国产三级专区第一集| 久久久久久伊人网av| 嫩草影院精品99| 麻豆乱淫一区二区| 97超碰精品成人国产| 亚洲精品久久久久久婷婷小说| 午夜激情久久久久久久| 午夜福利在线观看吧| 99久久人妻综合| 国产有黄有色有爽视频| 日韩精品有码人妻一区| 六月丁香七月| 国产精品久久久久久久久免| 国产又色又爽无遮挡免| 全区人妻精品视频| 日韩制服骚丝袜av| 一夜夜www| 好男人在线观看高清免费视频| av专区在线播放| 午夜日本视频在线| 国产精品国产三级国产av玫瑰| 国产亚洲精品久久久com| 青春草视频在线免费观看| 黄色一级大片看看| 街头女战士在线观看网站| 国产精品一及| 久久6这里有精品| 欧美精品一区二区大全| 美女黄网站色视频| 男人爽女人下面视频在线观看| 80岁老熟妇乱子伦牲交| 男女边摸边吃奶| 亚洲精品国产成人久久av| 六月丁香七月| 尤物成人国产欧美一区二区三区| 蜜桃亚洲精品一区二区三区| 两个人视频免费观看高清| 国产单亲对白刺激| a级一级毛片免费在线观看| 51国产日韩欧美| 又黄又爽又刺激的免费视频.| 日韩欧美 国产精品| 色吧在线观看| 五月天丁香电影| 成人一区二区视频在线观看| 日韩欧美 国产精品| a级毛片免费高清观看在线播放| 亚洲va在线va天堂va国产| av在线老鸭窝| 精品不卡国产一区二区三区| 天堂√8在线中文| 午夜福利视频1000在线观看| 男女那种视频在线观看| 国产探花极品一区二区| 日韩制服骚丝袜av| 在现免费观看毛片| 非洲黑人性xxxx精品又粗又长| 免费人成在线观看视频色| 婷婷色麻豆天堂久久| 国产高清有码在线观看视频| 尤物成人国产欧美一区二区三区| 网址你懂的国产日韩在线| 国产男人的电影天堂91| 能在线免费观看的黄片| 欧美变态另类bdsm刘玥| 亚洲精品视频女| 校园人妻丝袜中文字幕| av黄色大香蕉| 国产高潮美女av| 亚洲欧洲日产国产| 欧美成人一区二区免费高清观看| 亚洲欧美成人精品一区二区| 久久鲁丝午夜福利片| 亚洲国产av新网站| 久久国内精品自在自线图片| av在线老鸭窝| 日韩av在线免费看完整版不卡| 超碰av人人做人人爽久久| 最近的中文字幕免费完整| 久久鲁丝午夜福利片| 日韩欧美三级三区| 一级黄片播放器| 可以在线观看毛片的网站| 你懂的网址亚洲精品在线观看| 黑人高潮一二区| 亚洲人与动物交配视频| 99视频精品全部免费 在线| 在线观看免费高清a一片| 日韩人妻高清精品专区| 中文字幕免费在线视频6| 一级毛片黄色毛片免费观看视频| 国产男人的电影天堂91| 又爽又黄无遮挡网站| 99视频精品全部免费 在线| 久久综合国产亚洲精品| 成年av动漫网址| 大话2 男鬼变身卡| 一级毛片久久久久久久久女| 老师上课跳d突然被开到最大视频| 男人舔奶头视频| 久久午夜福利片| 亚洲欧美成人精品一区二区| 亚洲综合色惰| 国产精品一二三区在线看| 人体艺术视频欧美日本| 91在线精品国自产拍蜜月| 黄色配什么色好看| 久久97久久精品| 成年版毛片免费区| 色播亚洲综合网| 黄色一级大片看看| 少妇猛男粗大的猛烈进出视频 | 国产69精品久久久久777片| 国产精品久久久久久久电影| 深爱激情五月婷婷| 少妇熟女aⅴ在线视频| 校园人妻丝袜中文字幕| 少妇熟女aⅴ在线视频| 日日摸夜夜添夜夜爱| 成人高潮视频无遮挡免费网站| 欧美一区二区亚洲| 麻豆成人av视频| 丰满乱子伦码专区| 黄色欧美视频在线观看| 国产成人免费观看mmmm| 日本一本二区三区精品| av免费观看日本| 午夜精品国产一区二区电影 | 亚洲最大成人av| 老司机影院成人| 大话2 男鬼变身卡| 久久久久网色| 伦精品一区二区三区| xxx大片免费视频| 亚洲精品国产av成人精品| 国产一级毛片七仙女欲春2| 午夜精品在线福利| kizo精华| 日日干狠狠操夜夜爽| 亚洲精品乱码久久久久久按摩| 国产中年淑女户外野战色| 亚洲,欧美,日韩| 听说在线观看完整版免费高清| 国产91av在线免费观看| 黄片无遮挡物在线观看| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 成人特级av手机在线观看| a级一级毛片免费在线观看| 日日啪夜夜撸| 如何舔出高潮| 热99在线观看视频| 欧美激情在线99| 岛国毛片在线播放| 少妇人妻一区二区三区视频| 亚洲精品影视一区二区三区av| 欧美一级a爱片免费观看看| 亚洲精品乱久久久久久| 久久草成人影院| 日产精品乱码卡一卡2卡三| 综合色av麻豆| 建设人人有责人人尽责人人享有的 | 99久久九九国产精品国产免费| 国产男人的电影天堂91| 日韩欧美一区视频在线观看 | or卡值多少钱| 特级一级黄色大片| 床上黄色一级片| 在线免费十八禁| 国产午夜福利久久久久久| 久久午夜福利片| 一级爰片在线观看| 国产国拍精品亚洲av在线观看| 国产久久久一区二区三区| xxx大片免费视频| 免费看不卡的av| 97在线视频观看| 国内精品一区二区在线观看| a级一级毛片免费在线观看| 精品久久久久久久久av| 亚洲丝袜综合中文字幕| 亚洲av一区综合| 色尼玛亚洲综合影院| 日韩av免费高清视频| 中文字幕免费在线视频6|