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

    毫秒脈沖星觀測反演臺站原子鐘行為研究

    2023-12-02 07:05:50王嘯威閆振沈志強王玲玲吳亞軍趙融冰黃志鵬王睿
    時間頻率學報 2023年3期
    關(guān)鍵詞:脈沖星氫原子天馬

    王嘯威,閆振,3,4,*,沈志強,3,王玲玲,吳亞軍,3,趙融冰,3,黃志鵬,5,王睿

    毫秒脈沖星觀測反演臺站原子鐘行為研究

    王嘯威1,2,閆振1,2,3,4,*,沈志強1,2,3,王玲玲1,3,吳亞軍1,2,3,趙融冰1,2,3,黃志鵬1,2,5,王睿1,2

    (1. 中國科學院 上海天文臺,上海 200030;2. 中國科學院大學 天文與空間科學學院,北京 100049;3. 上海市引力波探測前沿科學研究基地,上海 200240;4. 中國科學院 射電天文重點實驗室,南京 210008;5. 上??萍即髮W 物質(zhì)科學與技術(shù)學院,上海 201210)

    氫原子鐘作為高精度時間和頻率系統(tǒng)的核心部件,在射電天文觀測、高精度時間計量、空間導(dǎo)航等領(lǐng)域獲得廣泛應(yīng)用。盡管氫原子鐘具有極高的短期穩(wěn)定性,但是其長期穩(wěn)定度會受到各種環(huán)境因素的影響。毫秒脈沖星因其超乎尋常的長期穩(wěn)定自轉(zhuǎn)頻率特性被譽為“掛在天上的鐘”,在時間保持、深空導(dǎo)航等方面具有重要的潛在應(yīng)用價值。本文基于上海天馬望遠鏡毫秒脈沖星計時觀測資料對有關(guān)問題進行了研究,成功反演出臺站氫原子鐘的歷史行為,包括時鐘跳躍、換鐘等。通過分段擬合其計時殘差中偏離白噪聲的長期趨勢成功獲得氫原子鐘速。綜合考慮上述歷史行為后獲得了基于該毫秒脈沖星的時鐘改正文件,將其用于天馬望遠鏡其他多顆毫秒脈沖星計時分析并取得了非常理想的效果。

    毫秒脈沖星;PSR J1713+0747;脈沖星計時;氫原子鐘;鐘差改正

    0 引言

    氫原子鐘利用基態(tài)氫原子的超精細能級躍遷,輸出短期穩(wěn)定度極高的頻率信號(天穩(wěn)定度可達10-15),因而在導(dǎo)航、空間探測、守時授時和天文觀測等領(lǐng)域有著廣泛的應(yīng)用。不過,在實際應(yīng)用中氫原子鐘的性能會受到各種環(huán)境因素的影響,如溫度、濕度、振動和磁場等。其中溫度變化的影響尤為顯著,因為它會顯著影響氫原子鐘諧振腔輸出頻率的準確度和穩(wěn)定度[1]。因此,氫原子鐘的長期穩(wěn)定性相對較低,從而也會影響整個時頻系統(tǒng)的長期可靠性。如果能夠找到一種長期穩(wěn)定的信號源,并通過一定技術(shù)手段實現(xiàn)對氫原子鐘的校準,那將對提高其長期穩(wěn)定性大有裨益。

    毫秒脈沖星是宇宙空間中一類經(jīng)過“吸積再加速”的特殊年老脈沖星[2],其自轉(zhuǎn)周期比普通脈沖星更短(約為1.4~30 ms),周期變化率也比普通脈沖星小約5個數(shù)量級(僅有10-19~10-21s s-1)。和普通脈沖星相比,毫秒脈沖星的計時特性穩(wěn)定得多,因為它們極少像普通脈沖星一樣時常表現(xiàn)出一些自轉(zhuǎn)不穩(wěn)定性行為,主要包括周期躍變(glitch)和計時噪聲(timing noise)。所謂周期躍變是一種脈沖星自轉(zhuǎn)突然加快的現(xiàn)象[3],而計時噪聲表現(xiàn)為脈沖星自轉(zhuǎn)參數(shù)發(fā)生連續(xù)的、長時標的“隨機游走”。在目前已知所有的毫秒脈沖星中,只有PSR J1824-2452A[4]和PSR J0613-0200[5]觀測到發(fā)生過周期躍變。正是由于毫秒脈沖星的自轉(zhuǎn)周期極其穩(wěn)定,對其長期周期測量的準確度可優(yōu)于10-15水平[6],因而它可以被視為宇宙空間中的一種天然的具有長期穩(wěn)定性的信號源,被譽為“掛在天上的鐘”。

    脈沖星計時觀測是完成上述工作的基礎(chǔ)。所謂脈沖星計時觀測又稱脈沖到達時間(time of arrival,簡稱TOA)測量,是指精確測量脈沖星信號到達觀測天線的時間。臺站時間由氫原子鐘、全球定位系統(tǒng)(Global Positioning System,簡稱GPS)接收機、比對計數(shù)器等組成的時頻系統(tǒng)維持。利用觀測得到的脈沖星積分輪廓與其標準輪廓進行互相關(guān)精確獲得到達天線時間,然后再考慮行星歷表改正、相對論效應(yīng)等將其轉(zhuǎn)化為太陽系質(zhì)心處的實際的TOA,進而與脈沖星計時模型給出預(yù)報的TOA進行對比并得到計時殘差,通過最小二乘法對計時殘差進行多次擬合迭代可以得到更精確的脈沖星計時模型參數(shù)和擬合后的計時殘差,進而可以借此研究脈沖星基本性質(zhì)和內(nèi)部結(jié)構(gòu)等。如果對一批計時精度較高的毫秒脈沖星進行長期計時觀測,就可以組成脈沖星計時陣(Pulsar Timing Array,簡稱PTA)用以探測低頻引力波。目前已有國內(nèi)外研究人員嘗試利用脈沖星計時數(shù)據(jù)建立“綜合脈沖星時”、檢測國際原子時中可能存在的系統(tǒng)誤差以及臺站時鐘的不穩(wěn)定性等[10-12]。

    本文將介紹利用上海天馬射電望遠鏡(下面簡稱天馬望遠鏡)毫秒脈沖星長期計時觀測數(shù)據(jù)對臺站原子鐘歷史行為反演以及修正方面的研究進展?;诤撩朊}沖星觀測數(shù)據(jù)成功獲得臺站氫原子鐘的鐘差數(shù)據(jù)并反演出其歷史行為,然后通過擬合其計時殘差的長期變化趨勢在一定程度上改正了鐘差,從而提高了天馬望遠鏡毫秒脈沖星計時觀測的計時精度。

    1 毫秒脈沖星計時觀測

    利用上海天馬射電望遠鏡開展了毫秒脈沖星的長期計時觀測,并且按照脈沖星計時觀測的標準流程對脈沖信號進行在線消色散和周期折疊并與標準脈沖輪廓互相關(guān)得到脈沖到達時間,最后利用TEMPO2軟件包[13]對脈沖到達時間進行分析得到計時殘差。

    1.1 觀測簡介

    上海天馬射電望遠鏡是一臺口徑為65 m的全可動大型射電望遠鏡。因位于上海市松江區(qū)的天馬山附近,考慮到臺站射電干擾情況和毫秒脈沖星輻射特點,天馬望遠鏡在S和C波段對毫秒脈沖星開展長期計時觀測。S和C波段的中心頻率分別為2 250和4 820 MHz,有效帶寬分別為100和800 MHz。利用DIBAS數(shù)字終端[14]將其總帶寬分別劃分為1 MHz和2 MHz的子通道,采用非相干消色散在線折疊模式對觀測數(shù)據(jù)進行在線的消色散和周期折疊處理,然后輸出子積分時間為30 s的8-bit量化的PSRFITS格式文件并記錄到磁盤陣列。由于天馬望遠鏡是一臺既服務(wù)于月球探測、火星探測等國家重大工程項目又用于天文觀測研究的多功能綜合望遠鏡,因而受到望遠鏡觀測時間限制,在2022年6月以前對每顆毫秒脈沖星的單次觀測時長約為5~30 min(視脈沖星的流量密度大小而定)。在2022年6月完成一個階段的有關(guān)國家重大探測工程觀測任務(wù)以后,為了提高觀測所得的TOA的精度,幾顆有望取得更高計時精度的毫秒脈沖星的單次觀測時長增加至2~3 h,觀測間隔也進一步縮短(約為1~2周每次)。表1的1~7列分別展示了每顆毫秒脈沖星的名稱、自轉(zhuǎn)周期、1.4 GHz處的平均流量密度、觀測波段、觀測起止時間、時間跨度以及TOA的個數(shù)。

    表1 7顆毫秒脈沖星及計時觀測的基本參數(shù)

    注:1)此處為1.4 GHz的平均流量密度,數(shù)據(jù)來自Australia Telescope National Facility(ATNF)脈沖星數(shù)據(jù)庫(1.67版),https://www.atnf.csiro.au/research/pulsar/psrcat/;

    2)MJD為約化儒略日。

    1.2 脈沖到達時間測量

    觀測數(shù)據(jù)的處理主要采用被廣泛使用的PSRCHIVE等軟件包[15]。數(shù)據(jù)處理的過程分為以下3步:①使用CLFD軟件包[16]對觀測數(shù)據(jù)進行自動的射頻干擾(radio frequency interference,簡稱RFI)消除,并通過PSRCHIVE軟件包中的pazi命令進行交互式檢查并手動消除殘余的RFI;② 把每次觀測數(shù)據(jù)沿頻率軸和時間軸的方向進行折疊,從而得到一系列的高信噪比積分輪廓;③ 將這些積分輪廓與一個事先準備好的標準模版進行互相關(guān)從而得到每次觀測的TOA。而標準模版的生成與TOA的獲得過程如下:① 選取多次積分輪廓信噪比較高的觀測,利用psradd命令將積分輪廓進行疊加;② 用paas命令擬合疊加得到的積分輪廓生成標準模版;③利用pat命令使用Fourier Domain Monte-Carlo(FDM)算法將每次觀測得到的積分輪廓與標準輪廓進行互相關(guān)得到TOA及其誤差。表1的第8列展示了每顆毫秒脈沖星計時觀測得到的TOA誤差的中位數(shù)。

    1.3 計時分析

    本文的脈沖星計時分析使用了TEMPO2軟件包。除了之前已獲得的TOA數(shù)據(jù)外,脈沖星計時分析還需要有一個初始星表,其中包括脈沖星的自轉(zhuǎn)、位置、自行、視差、色散量和雙星軌道(若脈沖星處于雙星系統(tǒng)中)等脈沖星計時模型參數(shù)。PSR J1012+5307的初始星表選用北美納赫茲引力波天文臺脈沖星計時陣發(fā)布的12.5年的計時觀測得到的結(jié)果[17],其余6顆脈沖星(PSR J1022+1001、J1643-1224、J1713+0747、J1744-1134、J1909-3744和J1939+2134)的初始星表則選擇澳大利亞Parkes脈沖星計時陣(Parkes Pulsar Timing Array,簡稱PPTA)發(fā)布的第二批數(shù)據(jù)中的結(jié)果[18]。此外,計時分析中的太陽系行星歷表選用DE436。理論上脈沖到達太陽系質(zhì)心處的時間可以用如下公式表示[19]:

    在不進行任何參數(shù)擬合的情況下,通過TEMPO2計時分析得到的每顆毫秒脈沖星的計時殘差如圖1所示。圖中S和C波段的觀測數(shù)據(jù)分別用黑色十字和灰色空心三角表示,每幅圖的左下角標注了脈沖星名和以均方根rms表示的計時精度。從圖1中可以明顯看出計時殘差的變化趨勢在不同源、不同觀測波段間都有很強的相關(guān)性。由于毫秒脈沖星普遍擁有非常低的計時噪聲,因此在給定足夠準確的脈沖星計時模型的情況下,得到的計時殘差應(yīng)該是一種白噪聲。脈沖星計時模型的主要誤差包括脈沖星自轉(zhuǎn)不穩(wěn)定性、脈沖信號傳播路徑中的色散延遲誤差、太陽系行星歷表誤差、引力波和觀測站原子鐘的誤差[20],而這些誤差中只有測站原子鐘的誤差會在觀測不同方向的脈沖星時引入單極性的誤差。同時,天馬望遠鏡的毫秒脈沖星計時觀測中并未對提供時間參考的本地氫原子鐘的鐘差進行改正,因此圖1中計時殘差的變化趨勢可以認為是本地氫原子鐘的鐘差。

    圖1 7顆毫秒脈沖星的計時殘差

    2 利用毫秒脈沖星計時觀測進行鐘差改正

    在脈沖星計時模型足夠準確的情況下,對一顆時間噪聲很小的毫秒脈沖星進行計時觀測和分析得到的計時殘差應(yīng)該是理想的白噪聲。因此可以利用毫秒脈沖星的到達時間觀測對計時殘差中存在的原子鐘鐘差項的長期變化趨勢進行擬合,從而達到改正鐘差的效果。

    2.1 參考毫秒脈沖星選取

    利用單顆毫秒脈沖星的計時觀測數(shù)據(jù)進行觀測站原子鐘誤差的改正時,選取的毫秒脈沖星需要滿足以下條件:① 有足夠的觀測精度(TOA的誤差?。虎?計時噪聲小,即自轉(zhuǎn)穩(wěn)定性高;③ 覆蓋足夠長的觀測時間且觀測相對頻繁。因此,在表1所列出的7顆毫秒脈沖星中選取PSR J1713+0747作為參考毫秒脈沖星。

    PSR J1713+0747是一顆自轉(zhuǎn)周期約為4.57 ms的毫秒脈沖星,由Arecibo望遠鏡巡天發(fā)現(xiàn)[21]。PSR J1713+0747處在一個軌道周期約68 d的雙星系統(tǒng)中,其伴星是一顆質(zhì)量約為0.28 M⊙的白矮星[22]。PSR J1713+0747作為目前計時最為精確的毫秒脈沖星之一,國際上的脈沖星計時陣都對其有長期的計時觀測[23-25]。天馬望遠鏡在S波段對PSR J1713+0747進行了長期的計時觀測,總的時間跨度為5.5年。PSR J1713+0747的流量較強(在1.4 GHz的平均流量密度約為9.1 mJy),對其觀測可以獲得比較高的信噪比,觀測得到的TOA的誤差中位數(shù)為1.312 μs。在進行計時分析過程中發(fā)現(xiàn)PSR J1713+0747在約化儒略日(Modified Julian Date,簡稱MJD)59 320—59 321(2020年4月16日至17日)間發(fā)生了一次輪廓變化事件,隨后其平均脈沖輪廓在幾個月內(nèi)逐漸恢復(fù)到變化前的狀態(tài)[26-27]。PSR J1713+0747平均脈沖輪廓的變化也使得觀測到的TOA存在一個跳變和逐漸恢復(fù)的過程。參考J. Singha等人擬合此次輪廓變化對TOA的影響時采用的跳變加指數(shù)恢復(fù)的模型[28],利用Enterprise軟件包對輪廓變化引起的TOA變化進行擬合,得到在S波段上TOA的跳變?yōu)?7.5 μs,恢復(fù)時標約219 d。本文中展示的PSR J1713+0747的計時殘差均為修正了輪廓變化事件后的計時殘差。

    2.2 鐘差模型擬合與分析

    圖2中帶誤差棒的數(shù)據(jù)點為PSR J1713+0747的計時殘差(已改正輪廓變化帶來的影響)。利用上述方法成功反演出天馬望遠鏡時頻系統(tǒng)氫原子鐘的歷史行為,這些行為在如圖2所示的計時殘差圖上表現(xiàn)為跳變、轉(zhuǎn)折等。天馬望遠鏡毫秒脈沖星計時觀測中的時間參考由一臺氫原子鐘作為主鐘。圖2發(fā)生在2018年4月份的J1713+0747計時殘差的轉(zhuǎn)折對應(yīng)一次換鐘事件,2018年4月26日(MJD 58 234)天馬時頻系統(tǒng)主鐘由國產(chǎn)SOHM-4型88號鐘更換為國產(chǎn)SOHM-4型118號鐘(本編號為上海天文臺氫原子鐘的生產(chǎn)編號);2018年9月份的該源計時殘差轉(zhuǎn)折對應(yīng)另外一次換鐘事件,2018年9月15日(MJD 58 379)后時頻系統(tǒng)主鐘由國產(chǎn)SOHM-4型118號鐘更換為俄羅斯產(chǎn)VCH-1003M型氫原子鐘。

    注:圖2(a)黑色虛線和灰色區(qū)域分別代表換鐘的時間節(jié)點和跳變發(fā)生的時間區(qū)間

    采用加權(quán)的最小二乘法對PSR J1713+0747的計時殘差進行分段線性擬合,權(quán)重為TOA誤差的倒數(shù),分段點取兩次換鐘的時間節(jié)點。此外在第2次換鐘后還加入兩個跳變,跳變的時間點取相鄰兩次觀測時間的中點。擬合所得的鐘差模型表示為圖2(a)的黑色虛線。擬合得到的每一段的斜率可以表示不同時間段內(nèi)使用的氫原子鐘的鐘速,分別為-2.25×10-13s s-1(MJD 57 950—58 234)、5.05×10-13s s-1(MJD 58 235—58 379)和7.99×10-13s s-1(MJD 58 380—59 953)。2020年3月至12月底這一段時間內(nèi)觀測數(shù)據(jù)較為稀疏,這是因為在這一段時間內(nèi)集中開展了天線維保和“天問一號”火星探測觀測任務(wù)。但是,在這一段時間發(fā)現(xiàn)了兩次脈沖星計時殘差跳變事件,擬合得到的跳變幅度分別為40.3 μs和15.5 μs。通過查閱維保記錄發(fā)現(xiàn)第1次跳變是由于2020年4月19日至5月15日(MJD 58958—58984)期間維修鐘房空調(diào)斷電時氫原子鐘發(fā)生斷電重啟。在2020年10月中旬至12月下旬這段時間內(nèi),由于空調(diào)故障鐘房溫度出現(xiàn)較大波動,引起氫原子鐘跳變,其正對應(yīng)圖2(a)中擬合得到的第2次跳變。此外,通過查閱相關(guān)日志還得到了臺站氫原子鐘原始的鐘差數(shù)據(jù)(臺站氫原子鐘和GPS鐘的差值)并以黑色圓點表示在圖2(b)中。通過對比,可以看出擬合得到的鐘差模型很好地反演了氫原子鐘鐘差的實際走勢。同時,實測的鐘差數(shù)據(jù)存在一些由于測量誤差(或錯誤)導(dǎo)致的異常數(shù)據(jù)點,而通過脈沖星計時得到的計時殘差數(shù)據(jù)點異常點極少。

    2.3 鐘差改正后計時分析

    利用擬合所得的鐘差模型按時間順序生成一系列鐘差數(shù)據(jù)并保存為一個鐘差改正文件后,再使用TEMPO2軟件包進行計時分析,就可以利用這個鐘差改正文件來獲得經(jīng)過鐘差改正的計時殘差。圖3展示了7顆毫秒脈沖星經(jīng)過鐘差改正后的計時殘差情況。

    從圖3可以看出圖1中所表現(xiàn)出的不同源之間的具有相關(guān)性的計時殘差變化趨勢已經(jīng)消失。除PSR J1939+2134外,其余6顆毫秒脈沖星的計時殘差相比鐘差改正前均有明顯改善,從而證明了鐘差改正的有效性。其中,計時精度最高的PSR J1909-3744通過C波段0.6年的觀測可以獲得170 ns的計時精度。PSR J1939+2134的計時殘差在經(jīng)過鐘差改正后依然明顯偏離了白噪聲,但是先前研究已經(jīng)證明該脈沖星的計時殘差中存在明顯的紅噪聲。進一步結(jié)合天馬望遠鏡S波段以及PPTA在波長10 cm頻段獲得的PSR J1939+2134的TOA數(shù)據(jù)(PPTA的數(shù)據(jù)只取2015年后)進行計時分析,得到的計時殘差如圖4所示??梢钥闯鲭m然天馬的觀測得到的TOA誤差相對較大,但在與PPTA觀測重合的時間段內(nèi)兩者的計時殘差變化趨勢符合較好。因此,圖3中PSR J1939+2134計時殘差顯示出的紅噪聲可以認為是其本身的某種自轉(zhuǎn)不穩(wěn)定性導(dǎo)致的。

    圖4 PPTA及天馬望遠鏡鐘差改正后數(shù)據(jù)的PSR J1939+2134的計時殘差

    3 結(jié)論及展望

    總之,氫原子鐘的長期不穩(wěn)定性會導(dǎo)致用其作為時間參考時在脈沖星計時模型中引入額外的鐘差誤差。而毫秒脈沖星極高的自轉(zhuǎn)穩(wěn)定性使其可以達到很高的計時精度,因此利用對毫秒脈沖星的長期計時觀測就可以反推出計時觀測中使用的氫原子鐘的鐘差。本文的研究表明:① 對于單顆計時精度高的毫秒脈沖星,在給定足夠精確的脈沖星計時模型初始參數(shù)時,可以通過計時殘差來顯示觀測中作為時間參考的氫原子鐘的鐘差的短期跳變及長期變化趨勢,從而通過擬合達到改正鐘差的效果;② 利用對PSR J1713+0747的計時觀測進行鐘差改正后,天馬望遠鏡的毫秒脈沖星計時精度得到有效提高,其中PSR J1909-3744在C波段0.6 年的計時精度可以達到170 ns。

    我國有一定規(guī)模的脈沖星觀測始于20世紀90年代的烏魯木齊天文觀測站(現(xiàn)為新疆天文臺)25 m射電望遠鏡脈沖星觀測系統(tǒng)的建立[29]。和國際發(fā)達國家相比,可謂起步較晚。隨著國家經(jīng)濟實力的提高,國內(nèi)先后建起了密云50 m[30]、昆明40 m[31]、上海65 m、昊平40 m[32]、貴州500 m口徑球面射電望遠鏡[33]等大型射電望遠鏡。另外,奇臺110 m、景東120 m、修水 120 m等百米級別望遠鏡也已經(jīng)批準立項或正在順利推進。這些望遠鏡建設(shè)將極大提升我國脈沖星的觀測研究能力。各望遠鏡揚長避短,有機組成中國脈沖星計時陣列,這將不僅會使得國內(nèi)在脈沖星物理、納赫茲引力波探測等基礎(chǔ)領(lǐng)域占領(lǐng)國際制高點,而且還將在脈沖星時間標準建立與維持、深空自主導(dǎo)航等應(yīng)用領(lǐng)域開辟廣闊前景。

    [1] PARKER T E. Environmental factors and hydrogen maser frequency stability[J]. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 1999, 46(3): 745-751.

    [2] ALPAR M A, CHENG A F, RUDERMAN M A, et al. A new class of radio pulsars[J]. Nature, 1982(300): 728-730.

    [3] ESPINOZA C M, LYNE A G, STAPPERS B W, et al. A study of 315 glitches in the rotation of 102 pulsars[J]. Monthly Notices of the Royal Astronomical Society, 2011, 414(2): 1679-1704.

    [4] COGNARD I, BACKER D C. A microglitch in the millisecond pulsar PSR B1821-24 in M28[J]. The Astrophysical Journal, 2004, 612(2): 125-127.

    [5] MCKEE J W, JANSSEN G H, STAPPERS B W, et al. A glitch in the millisecond pulsar J0613-0200[J]. Monthly Notices of the Royal Astronomical Society, 2016, 461(3): 2809-2817.

    [6] VERBIEST J P W, BAILES M, VAN STRATEN W, et al. Precision timing of PSR J0437-4715: an accurate pulsar distance, a high pulsar mass, and a limit on the variation of Newton’s gravitational constant[J]. The Astrophysical Journal, 2008, 627(1): 675-680.

    [7] MATSAKIS D N, TAYLOR J H, EUBANKS T M. A statistic for describing pulsar and clock stabilities[J]. Astronomy and Astrophysics, 1997(326): 924-928.

    [8] VERBIEST J P W, BAILES M, COLES W A, et al. Timing stability of millisecond pulsars and prospects for gravitational-wave detection[J]. Monthly Notices of the Royal Astronomical Society, 2009, 400(2): 951-968.

    [9] 楊廷高, 仲崇霞. 脈沖星時穩(wěn)定度及可能應(yīng)用[J]. 時間頻率學報, 2004, 27(2): 129-137.

    [10] HOBBS G, COLES W, MANCHESTER R N, et al. Development of a pulsar-based time scale[J]. Monthly Notices of the Royal Astronomical Society, 2012, 427(4): 2780-2787.

    [11] HOBBS G, GOU L, CABALLERO R N, et al. A pulsar-based time-scale from the international pulsar timing array[J]. Monthly Notices of the Royal Astronomical Society, 2020, 491(4): 5951-5965.

    [12] 李志玄. 利用PSR J0437-47l5研究引力常數(shù)變化率和臺站時鐘跳變[D]. 云南: 中國科學院大學(中國科學院云南天文臺), 2019.

    [13] HOBBS G, EDWARDS R T, MANCHESTER R N. TEMPO2, a new pulsar-timing package-I. an overview[J]. Monthly Notices of the Royal Astronomical Society, 2006, 369(2): 655-672.

    [14] YAN Z, SHEN Z Q, WU Y J, et al. Pulsar research with the newly built Shanghai Tian Ma Radio Telescope[J]. URSI Radio Science Bulletin, 2018, 2018(366): 10-18.

    [15] HOTAN A, VAN STRATEN W, MANCHESTER R. Psrchive and psrfits: an open approach to radio pulsar data storage and analysis[J]. Publications of the Astronomical Society of Australia, 2004, 21(3): 302-309.

    [16] MORELLO V, BARR E D, COOPER S, et al. The high time resolution universe survey-XIV. Discovery of 23 pulsars through GPU-accelerated reprocessing[J]. Monthly Notices of the Royal Astronomical Society, 2019, 483(3): 3673-3685.

    [17] ALAM M F, ARZOUMANIAN Z, BAKER P T, et al. The NANOGrav 12.5 yr data set: wideband timing of 47 millisecond pulsars[J]. The Astrophysical Journal Supplement Series, 2021, 252(4): 54.

    [18] REARDON D J, SHANNON R M, CAMERON A D, et al. The Parkes Pulsar Timing Array second data release: timing analysis[J]. Monthly Notices of the Royal Astronomical Society, 2021, 507(2): 2137-2153.

    [19] NORBERT W. A timing formula for main-sequence star binary pulsars[J]. Monthly Notices of the Royal Astronomical Society, 1998, 298(1): 66-77.

    [20] EDWARDS R T, HOBBS G, MANCHESTER R N. TEMPO2, a new pulsar timing package-II. The timing model and precision estimates[J]. Monthly Notices of the Royal Astronomical Society, 2006, 372(4): 1549-1574.

    [21] FOSTER R S, WOLSZCZAN A, CAMILO F. A new binary millisecond pulsar[J]. Astrophysical Journal Letters, 1993(410): 91.

    [22] SPLAVER E M, NICE D J, STAIRS I H, et al. Masses, parallax, and relativistic timing of the PSR J1713+0747 binary system[J]. The Astrophysical Journal, 2005, 620(1): 405-415.

    [23] DESVIGNES G, CABALLERO R N, LENTATI L, et al. High-precision timing of 42 millisecond pulsars with the European Pulsar Timing Array[J]. Monthly Notices of the Royal Astronomical Society, 2016, 458(3): 3341-3380.

    [24] MCLAUGHLIN M A. The North American Nanohertz Observatory for gravitational waves[J]. Classical and Quantum Gravity, 2013, 30(22): 224008.

    [25] HOBBS G. The Parkes pulsar timing array[J]. Classical and Quantum Gravity, 2013, 30(22): 224007.

    [26] XU H, HUANG Y X, BURGAY M, et al. A sustained pulse shape change in PSR J1713+0747 possibly associated with timing and DM events[J]. The Astronomer’s Telegram, 2021, 14642: 1.

    [27] MEYERS B, CHIME/PULSAR COLLABORATION. Confirmation of a change in the emission properties of PSR J1713+0747[J]. The Astronomer’s Telegram, 2021, 14652: 1.

    [28] SINGHA J, SURNIS M P, JOSHI B C, et al. Evidence for profile changes in PSR J1713+0747 using the uGMRT[J]. Monthly Notices of the Royal Astronomical Society, 2021, 507(1): 57-61.

    [29] 王娜, 張晉, 艾力·玉,等. 烏魯木齊天文站的脈沖星脈沖到達時間觀測[C] // 1999年全國時間頻率學術(shù)報告會, 青島: 1999年全國時間頻率學術(shù)報告委員會.

    [30] JIN C, CAO Y, CHEN H, et al. The Miyun 50 m Pulsar Radio Telescope[J]. Chinese Journal of Astronomy and Astrophysics Supplement, 2006, 6(S2): 319-323.

    [31] HAO L F, WANG M, YANG J. VLBI observations with the Kunming 40-meter radio telescope[J]. Research in Astronomy and Astrophysics, 2010, 10(8): 805-814.

    [32] LUO J T, GAO Y P, YANG T G, et al. Pulsar timing observations with Haoping radio telescope[J]. Research in Astronomy and Astrophysics, 2020, 20(7): 111.

    [33] NAN R D. Five hundred meter aperture spherical radio telescope (FAST)[J]. Science in China: Physics, Mechanics and Astronomy, 2006, 49(2): 129-148.

    Retrieval of the hydrogen maser’s behavior at local station using observations of millisecond pulsars

    WANG Xiao-wei1, 2, YAN Zhen1, 2, 3, 4,*, SHEN Zhi-qiang1, 2, 3, WANG Ling-ling1, 3, WU Ya-jun1, 2, 3, ZHAO Rong-bing1, 2, 3, HUANG Zhi-peng1, 2, 5, WANG Rui1, 2

    (1. Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai 200030, China;2. University of Chinese Academy of Sciences, Beijing 100049, China;3. Shanghai Frontier Science Center of Gravitational Wave Detection, Shanghai 200240, China; 4. Key Laboratory of Radio Astronomy, Chinese Academy of Sciences, Nanjing 210008, China;5. School of Physical Science and Technology, Shanghai Tech University, Shanghai 201210, China)

    As the core component of a high-precision time and frequency system, the hydrogen maser clock has been widely used in radio astronomy observation, high-precision time measurement, space navigation and other fields. Although the hydrogen maser clock is extremely stable in the short term, its long-term stability could be affected by various environmental factors. The millisecond pulsar (MSP) is known as the “clock on the sky” because of its extraordinary long-term stable rotation frequency characteristics and its important potential application value in time maintenance, deep space navigation and other fields. Deep studies in these fields were done based on the timing observation data on MSPs with the Shanghai Tianma Radio Telescope (TMRT), and the historical behavior of the hydrogen atomic clock at the station, including clock jumping and changing, were successfully retrieved. The hydrogen maser clock speed is successfully obtained by fitting the long-term trend of deviation from white noise in timing residuals. Considering the above historical behavior, we obtained the clock-correction file based on pulsar timing. We applied this file to the timing analysis of other MSPs in the TMRT and achieved very good results.

    millisecond pulsar; PSR J1713+0747; pulsar timing; hydrogen maser clock; clock correction

    2023-02-25;

    2023-05-29

    國家重點研發(fā)計劃資助項目(2022YFA1603104);上海市自然科學基金資助項目(20ZR1467600)

    王嘯威, 閆振, 沈志強, 等. 毫秒脈沖星觀測反演臺站原子鐘行為研究[J]. 時間頻率學報, 2023, 46(3): 217-226.

    10.13875/j.issn.1674-0637.2023-03-0217-10

    猜你喜歡
    脈沖星氫原子天馬
    “中國天眼”已發(fā)現(xiàn)740余顆新脈沖星
    軍事文摘(2023年12期)2023-06-12 07:51:00
    發(fā)現(xiàn)脈沖星的女天文學家——貝爾
    科學(2022年4期)2022-10-25 02:43:42
    一維模型氫原子在啁啾激光場中的閾上電離
    “天馬”鬛羚有絕技
    關(guān)于氫原子躍遷問題的難點突破
    基于虛擬觀測值的X射線單脈沖星星光組合導(dǎo)航
    對氫原子能量吸收的再認識お
    衢州市常山縣天馬一小
    氫原子n=2能級的二級斯塔克效應(yīng)
    如浮云伴天馬,共演西游神話
    久久精品亚洲熟妇少妇任你| 9色porny在线观看| 一级毛片女人18水好多| 日韩中文字幕欧美一区二区| 日韩精品免费视频一区二区三区| 高清黄色对白视频在线免费看| 免费在线观看黄色视频的| 国产人伦9x9x在线观看| 亚洲熟女毛片儿| 国产高清激情床上av| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 国产精品野战在线观看 | 国产三级黄色录像| 欧美人与性动交α欧美精品济南到| 亚洲avbb在线观看| 国产乱人伦免费视频| 成人三级黄色视频| 免费人成视频x8x8入口观看| 亚洲欧美精品综合久久99| 咕卡用的链子| 精品久久久久久,| av在线天堂中文字幕 | 在线观看免费视频日本深夜| 黄片小视频在线播放| 亚洲国产欧美一区二区综合| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 亚洲美女黄片视频| 亚洲国产精品一区二区三区在线| 亚洲国产中文字幕在线视频| 亚洲五月色婷婷综合| 男人舔女人下体高潮全视频| 亚洲一码二码三码区别大吗| 高清在线国产一区| 精品国产亚洲在线| 99精品欧美一区二区三区四区| 十八禁人妻一区二区| 国产av一区在线观看免费| 欧美日韩av久久| 亚洲成人免费av在线播放| 99热只有精品国产| 中文字幕另类日韩欧美亚洲嫩草| 性欧美人与动物交配| 成熟少妇高潮喷水视频| 老熟妇乱子伦视频在线观看| 精品国产一区二区三区四区第35| 中亚洲国语对白在线视频| 免费久久久久久久精品成人欧美视频| 亚洲欧美激情在线| 啦啦啦 在线观看视频| 午夜精品在线福利| 看黄色毛片网站| 91字幕亚洲| 国产精品日韩av在线免费观看 | 国产精品电影一区二区三区| 亚洲精品在线观看二区| 精品免费久久久久久久清纯| 亚洲男人的天堂狠狠| 老熟妇乱子伦视频在线观看| 80岁老熟妇乱子伦牲交| tocl精华| 精品卡一卡二卡四卡免费| 精品电影一区二区在线| 久久精品人人爽人人爽视色| 可以免费在线观看a视频的电影网站| 好看av亚洲va欧美ⅴa在| 亚洲伊人色综图| 国产精品爽爽va在线观看网站 | 女生性感内裤真人,穿戴方法视频| 欧美色视频一区免费| 免费在线观看日本一区| 91字幕亚洲| 国产精品爽爽va在线观看网站 | 叶爱在线成人免费视频播放| 亚洲欧美精品综合久久99| 欧美激情久久久久久爽电影 | 国产成人精品久久二区二区91| 亚洲一区二区三区欧美精品| 亚洲男人的天堂狠狠| av免费在线观看网站| 国产av一区二区精品久久| 母亲3免费完整高清在线观看| www.精华液| 欧美精品亚洲一区二区| 欧美黑人欧美精品刺激| 精品无人区乱码1区二区| 黑人操中国人逼视频| 好男人电影高清在线观看| 久久精品人人爽人人爽视色| 中出人妻视频一区二区| 亚洲三区欧美一区| 欧美日韩中文字幕国产精品一区二区三区 | 欧美黑人精品巨大| 国产av又大| 精品国内亚洲2022精品成人| 亚洲精品在线观看二区| 久久久久久久久久久久大奶| 级片在线观看| 男人舔女人的私密视频| 18禁美女被吸乳视频| 校园春色视频在线观看| 亚洲 欧美 日韩 在线 免费| 女性生殖器流出的白浆| 狂野欧美激情性xxxx| 又黄又爽又免费观看的视频| 免费在线观看完整版高清| 香蕉国产在线看| 久久 成人 亚洲| 亚洲伊人色综图| 黄频高清免费视频| 欧美大码av| 亚洲精品中文字幕在线视频| 欧美日韩亚洲综合一区二区三区_| 久久久国产成人精品二区 | 大香蕉久久成人网| 亚洲欧美一区二区三区久久| 国产主播在线观看一区二区| 久9热在线精品视频| 一级毛片精品| 亚洲欧美日韩无卡精品| 日韩高清综合在线| 91在线观看av| 午夜免费观看网址| 亚洲人成伊人成综合网2020| 亚洲欧美精品综合一区二区三区| 黄色视频不卡| 操美女的视频在线观看| 又紧又爽又黄一区二区| 久久香蕉国产精品| 精品欧美一区二区三区在线| 999精品在线视频| 国产免费现黄频在线看| 美女扒开内裤让男人捅视频| 99久久久亚洲精品蜜臀av| 午夜福利在线免费观看网站| 国产精品永久免费网站| 男人舔女人的私密视频| 久久久久久久久久久久大奶| 国产1区2区3区精品| 神马国产精品三级电影在线观看 | av免费在线观看网站| 岛国视频午夜一区免费看| 国产极品粉嫩免费观看在线| 欧美不卡视频在线免费观看 | 窝窝影院91人妻| 精品人妻1区二区| 男女下面插进去视频免费观看| 国产高清国产精品国产三级| 亚洲伊人色综图| 欧美在线一区亚洲| 国产精品野战在线观看 | 最近最新中文字幕大全电影3 | 亚洲 欧美一区二区三区| 国产精品影院久久| 高清毛片免费观看视频网站 | 1024香蕉在线观看| 日韩一卡2卡3卡4卡2021年| 在线十欧美十亚洲十日本专区| 国产精品免费一区二区三区在线| 日本撒尿小便嘘嘘汇集6| 亚洲一区中文字幕在线| 啦啦啦 在线观看视频| 狂野欧美激情性xxxx| 精品少妇一区二区三区视频日本电影| 亚洲av电影在线进入| 亚洲精品国产一区二区精华液| 少妇 在线观看| 国产午夜精品久久久久久| 国产成人欧美在线观看| 在线观看66精品国产| 日本黄色日本黄色录像| 亚洲avbb在线观看| 国产深夜福利视频在线观看| 视频在线观看一区二区三区| 国产亚洲精品综合一区在线观看 | 9色porny在线观看| 巨乳人妻的诱惑在线观看| 99riav亚洲国产免费| 热re99久久国产66热| 亚洲第一av免费看| 在线观看免费午夜福利视频| 一二三四在线观看免费中文在| 国产欧美日韩精品亚洲av| 国产区一区二久久| 999精品在线视频| 日韩欧美免费精品| 精品久久久久久久久久免费视频 | 精品熟女少妇八av免费久了| 日韩国内少妇激情av| 国产91精品成人一区二区三区| 日韩有码中文字幕| 一级黄色大片毛片| 俄罗斯特黄特色一大片| 国产人伦9x9x在线观看| 亚洲欧美日韩无卡精品| bbb黄色大片| 日韩精品中文字幕看吧| 亚洲成av片中文字幕在线观看| 久久久久久亚洲精品国产蜜桃av| 亚洲精品中文字幕在线视频| 久久精品91无色码中文字幕| 美国免费a级毛片| 91麻豆精品激情在线观看国产 | 久久狼人影院| 日本一区二区免费在线视频| 夜夜看夜夜爽夜夜摸 | 精品国内亚洲2022精品成人| 国产精品国产av在线观看| 人人妻人人澡人人看| 欧美日韩av久久| 国产精品九九99| 亚洲国产中文字幕在线视频| 国产野战对白在线观看| 久久性视频一级片| 老司机在亚洲福利影院| 天天添夜夜摸| 夜夜躁狠狠躁天天躁| 亚洲精品国产精品久久久不卡| 欧美丝袜亚洲另类 | 97碰自拍视频| a在线观看视频网站| 香蕉久久夜色| 啦啦啦 在线观看视频| 亚洲男人的天堂狠狠| 欧美日韩黄片免| 日韩欧美免费精品| 亚洲中文日韩欧美视频| 久热这里只有精品99| 99热国产这里只有精品6| 9热在线视频观看99| 亚洲中文av在线| 亚洲国产中文字幕在线视频| 亚洲国产精品sss在线观看 | 精品一区二区三区视频在线观看免费 | 国产精品秋霞免费鲁丝片| 九色亚洲精品在线播放| 黄色视频,在线免费观看| 国产精品一区二区免费欧美| 久久精品国产亚洲av香蕉五月| 国产国语露脸激情在线看| 亚洲欧美精品综合一区二区三区| 狠狠狠狠99中文字幕| 欧美最黄视频在线播放免费 | 亚洲一区高清亚洲精品| 三上悠亚av全集在线观看| 最近最新免费中文字幕在线| 757午夜福利合集在线观看| 在线天堂中文资源库| 精品久久久久久成人av| 天天躁狠狠躁夜夜躁狠狠躁| 国产国语露脸激情在线看| 在线观看免费视频日本深夜| 中文字幕另类日韩欧美亚洲嫩草| 美女福利国产在线| 一二三四社区在线视频社区8| 女人被狂操c到高潮| 超碰97精品在线观看| 69av精品久久久久久| 一区二区三区国产精品乱码| 免费高清视频大片| 成人影院久久| 亚洲片人在线观看| 精品午夜福利视频在线观看一区| 一进一出抽搐gif免费好疼 | 免费看a级黄色片| 午夜精品在线福利| 亚洲三区欧美一区| 久久精品91无色码中文字幕| 手机成人av网站| 90打野战视频偷拍视频| 淫妇啪啪啪对白视频| 日韩 欧美 亚洲 中文字幕| x7x7x7水蜜桃| 亚洲一区二区三区色噜噜 | 亚洲专区字幕在线| 成人免费观看视频高清| 欧美最黄视频在线播放免费 | 亚洲欧美精品综合久久99| 97碰自拍视频| 女警被强在线播放| 亚洲国产欧美日韩在线播放| 国产黄色免费在线视频| 手机成人av网站| 国产有黄有色有爽视频| 亚洲午夜理论影院| 国产一区二区三区综合在线观看| 9色porny在线观看| 欧美老熟妇乱子伦牲交| 欧美另类亚洲清纯唯美| 国产又色又爽无遮挡免费看| 婷婷丁香在线五月| 欧美av亚洲av综合av国产av| 亚洲一区高清亚洲精品| 成人国产一区最新在线观看| 午夜精品久久久久久毛片777| 一级a爱视频在线免费观看| 欧美激情高清一区二区三区| 免费一级毛片在线播放高清视频 | 男人操女人黄网站| 又紧又爽又黄一区二区| 1024视频免费在线观看| 激情在线观看视频在线高清| 亚洲性夜色夜夜综合| 怎么达到女性高潮| 丰满的人妻完整版| 伦理电影免费视频| av欧美777| 亚洲一区高清亚洲精品| 无遮挡黄片免费观看| 黄色片一级片一级黄色片| 999精品在线视频| 18禁裸乳无遮挡免费网站照片 | 国产成人啪精品午夜网站| 精品久久久久久,| 国产麻豆69| 99久久久亚洲精品蜜臀av| 国产精品成人在线| av免费在线观看网站| 最近最新中文字幕大全免费视频| www日本在线高清视频| 嫩草影院精品99| 91精品国产国语对白视频| 美女高潮喷水抽搐中文字幕| 黄片播放在线免费| 国产成年人精品一区二区 | 久久精品人人爽人人爽视色| 一级毛片高清免费大全| 亚洲熟妇熟女久久| 欧美最黄视频在线播放免费 | 欧美中文综合在线视频| 男女下面进入的视频免费午夜 | 一夜夜www| 美女福利国产在线| 国产一区二区三区在线臀色熟女 | 日本a在线网址| 丝袜美足系列| 亚洲人成77777在线视频| 久久国产亚洲av麻豆专区| 中文字幕人妻丝袜一区二区| 人人妻人人添人人爽欧美一区卜| 男女床上黄色一级片免费看| xxxhd国产人妻xxx| 免费高清视频大片| 国产成人免费无遮挡视频| 亚洲熟妇熟女久久| 免费一级毛片在线播放高清视频 | 精品熟女少妇八av免费久了| 大型黄色视频在线免费观看| 精品无人区乱码1区二区| 大型黄色视频在线免费观看| 亚洲人成伊人成综合网2020| 91成人精品电影| 97超级碰碰碰精品色视频在线观看| 日韩精品青青久久久久久| 69精品国产乱码久久久| 亚洲国产欧美日韩在线播放| 18禁黄网站禁片午夜丰满| 老熟妇仑乱视频hdxx| 欧美色视频一区免费| 身体一侧抽搐| 日本一区二区免费在线视频| 人人妻,人人澡人人爽秒播| 视频区图区小说| 三上悠亚av全集在线观看| 欧美日韩一级在线毛片| 久久欧美精品欧美久久欧美| 十分钟在线观看高清视频www| 神马国产精品三级电影在线观看 | 国产不卡一卡二| 欧美一级毛片孕妇| 亚洲第一av免费看| 老司机福利观看| 巨乳人妻的诱惑在线观看| 精品午夜福利视频在线观看一区| 欧美黑人精品巨大| 最新美女视频免费是黄的| 18禁国产床啪视频网站| 成人18禁高潮啪啪吃奶动态图| 国产av一区二区精品久久| 国产av一区在线观看免费| 国产成年人精品一区二区 | 亚洲狠狠婷婷综合久久图片| 一本综合久久免费| 男女下面插进去视频免费观看| 亚洲一区中文字幕在线| 亚洲一区二区三区欧美精品| 在线观看www视频免费| 男人的好看免费观看在线视频 | 在线观看免费午夜福利视频| 婷婷精品国产亚洲av在线| 国产精华一区二区三区| 成人永久免费在线观看视频| 少妇裸体淫交视频免费看高清 | 免费在线观看影片大全网站| 三级毛片av免费| 亚洲伊人色综图| 最新在线观看一区二区三区| 天天影视国产精品| 国产精品永久免费网站| 99精品欧美一区二区三区四区| 黄色怎么调成土黄色| 亚洲av成人av| 亚洲av成人一区二区三| 欧美日韩乱码在线| 国产乱人伦免费视频| 18美女黄网站色大片免费观看| 在线观看66精品国产| 国产精品偷伦视频观看了| 久久精品91蜜桃| 99在线视频只有这里精品首页| 亚洲欧美激情在线| 中文字幕人妻熟女乱码| 国产精品影院久久| 色综合欧美亚洲国产小说| 精品国产一区二区三区四区第35| 国产高清videossex| 国产成年人精品一区二区 | а√天堂www在线а√下载| 激情在线观看视频在线高清| 久久久久九九精品影院| 亚洲国产欧美日韩在线播放| 日本撒尿小便嘘嘘汇集6| 国产成人系列免费观看| 国内久久婷婷六月综合欲色啪| 国产精品国产av在线观看| 国产在线精品亚洲第一网站| 久久久国产精品麻豆| 9191精品国产免费久久| 中文字幕另类日韩欧美亚洲嫩草| 97人妻天天添夜夜摸| 视频在线观看一区二区三区| 亚洲成人久久性| 日本wwww免费看| 女人高潮潮喷娇喘18禁视频| 老司机福利观看| 一进一出好大好爽视频| 日本a在线网址| 久久精品亚洲熟妇少妇任你| 视频在线观看一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 日本撒尿小便嘘嘘汇集6| 一级a爱片免费观看的视频| 成人18禁在线播放| 中文字幕人妻丝袜一区二区| 亚洲av成人一区二区三| 69av精品久久久久久| 中文亚洲av片在线观看爽| 一级作爱视频免费观看| 国产三级在线视频| 免费av中文字幕在线| 欧美黑人欧美精品刺激| 午夜日韩欧美国产| 最近最新中文字幕大全免费视频| 中文亚洲av片在线观看爽| 别揉我奶头~嗯~啊~动态视频| 人人妻,人人澡人人爽秒播| 亚洲欧美精品综合久久99| 久久天堂一区二区三区四区| 国产精品电影一区二区三区| 免费女性裸体啪啪无遮挡网站| 老熟妇乱子伦视频在线观看| 制服诱惑二区| 99国产极品粉嫩在线观看| 人妻久久中文字幕网| 国产欧美日韩一区二区三| 国产99久久九九免费精品| 国产一区二区三区在线臀色熟女 | 欧美精品一区二区免费开放| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲欧美精品永久| 国产精品二区激情视频| 国产亚洲精品一区二区www| 51午夜福利影视在线观看| 国产精品1区2区在线观看.| 午夜精品在线福利| 村上凉子中文字幕在线| 国产亚洲欧美精品永久| 亚洲熟妇中文字幕五十中出 | 人成视频在线观看免费观看| 身体一侧抽搐| 性欧美人与动物交配| 精品国产美女av久久久久小说| 看免费av毛片| 一夜夜www| 亚洲成人免费电影在线观看| 精品国产一区二区久久| av欧美777| 黄色怎么调成土黄色| 国产精品av久久久久免费| 欧美日韩av久久| 波多野结衣一区麻豆| 国产av在哪里看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品美女久久久久99蜜臀| 看黄色毛片网站| 免费人成视频x8x8入口观看| 少妇粗大呻吟视频| 精品午夜福利视频在线观看一区| 男女床上黄色一级片免费看| 久久精品亚洲精品国产色婷小说| a在线观看视频网站| 欧美+亚洲+日韩+国产| 欧美另类亚洲清纯唯美| 日韩三级视频一区二区三区| www日本在线高清视频| 久久九九热精品免费| 最好的美女福利视频网| 国产人伦9x9x在线观看| 精品久久久久久久毛片微露脸| 亚洲色图av天堂| 真人一进一出gif抽搐免费| 久久久国产欧美日韩av| 久久热在线av| 久久人妻福利社区极品人妻图片| 亚洲成人免费电影在线观看| 桃色一区二区三区在线观看| 丰满的人妻完整版| 亚洲av成人一区二区三| 在线观看午夜福利视频| 午夜免费成人在线视频| 成人影院久久| 黄色视频,在线免费观看| 亚洲五月色婷婷综合| 美女 人体艺术 gogo| 精品人妻在线不人妻| 国产1区2区3区精品| 精品福利永久在线观看| 999久久久国产精品视频| 校园春色视频在线观看| 女生性感内裤真人,穿戴方法视频| 欧美日韩瑟瑟在线播放| 69精品国产乱码久久久| 在线观看一区二区三区| 国产激情欧美一区二区| cao死你这个sao货| 精品电影一区二区在线| 国产精品 欧美亚洲| 老汉色av国产亚洲站长工具| 久久精品亚洲熟妇少妇任你| 国产欧美日韩一区二区精品| 免费少妇av软件| 免费在线观看日本一区| 男女下面插进去视频免费观看| 女人被躁到高潮嗷嗷叫费观| 美女扒开内裤让男人捅视频| 久久久久久久午夜电影 | 99在线视频只有这里精品首页| 日本欧美视频一区| 亚洲精品av麻豆狂野| av网站免费在线观看视频| 老熟妇仑乱视频hdxx| 欧美中文综合在线视频| 成人亚洲精品一区在线观看| 一级黄色大片毛片| 国产亚洲欧美98| 美女高潮喷水抽搐中文字幕| 欧美性长视频在线观看| 中文字幕人妻熟女乱码| 男女做爰动态图高潮gif福利片 | 老熟妇乱子伦视频在线观看| 国产国语露脸激情在线看| 少妇被粗大的猛进出69影院| 热re99久久国产66热| 又大又爽又粗| 男女午夜视频在线观看| 夜夜躁狠狠躁天天躁| 久久久久久久精品吃奶| 韩国精品一区二区三区| 制服人妻中文乱码| 一夜夜www| 欧美日韩亚洲高清精品| 欧美成人免费av一区二区三区| 免费高清视频大片| 777久久人妻少妇嫩草av网站| 亚洲人成伊人成综合网2020| 黄片小视频在线播放| 亚洲va日本ⅴa欧美va伊人久久| 后天国语完整版免费观看| avwww免费| 亚洲av成人一区二区三| 91国产中文字幕| 日韩成人在线观看一区二区三区| 黄色怎么调成土黄色| 精品福利永久在线观看| 精品人妻1区二区| 亚洲在线自拍视频| 国产成人精品无人区| 亚洲 国产 在线| 天堂影院成人在线观看| 中文字幕另类日韩欧美亚洲嫩草| 亚洲成人久久性| 成人18禁在线播放| 日韩欧美免费精品| 一进一出抽搐gif免费好疼 | 国产亚洲欧美在线一区二区| 欧美日韩黄片免| 亚洲男人的天堂狠狠| 嫩草影院精品99| 亚洲精品中文字幕在线视频| www.精华液| 99国产综合亚洲精品| 亚洲国产欧美日韩在线播放| 成人18禁在线播放| 日本免费a在线| 成人av一区二区三区在线看| 黄色丝袜av网址大全| 天堂中文最新版在线下载| 18禁观看日本| 精品第一国产精品| 日韩av在线大香蕉| 男女午夜视频在线观看| 亚洲 欧美一区二区三区| 久久久久国内视频| 啦啦啦免费观看视频1|