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

    短周期彗星表面水冰升華分布研究*

    2023-02-01 08:39:34劉燦趙玉暉季江徽
    天文學(xué)報 2023年1期
    關(guān)鍵詞:水冰進(jìn)動面元

    劉燦 趙玉暉 季江徽

    (1 中國科學(xué)院紫金山天文臺 南京 210023)

    (2 中國科學(xué)技術(shù)大學(xué)天文與空間科學(xué)學(xué)院 合肥 230026)

    (3 中國科學(xué)院行星科學(xué)重點實驗室 南京 210023)

    1 引言

    太陽系內(nèi)的彗星大多起源于柯伊伯帶或更遠(yuǎn)的奧爾特云,這些區(qū)域的冰質(zhì)小天體受到引力擾動進(jìn)而偏離原軌道進(jìn)入太陽系內(nèi)部成為彗星[1].彗核是構(gòu)成彗星的固體部分,當(dāng)靠近太陽時,強(qiáng)烈的太陽輻射會加熱彗核表面以及亞表層的物質(zhì),使其升華產(chǎn)生氣體并攜帶著塵埃向外噴發(fā),彗核周圍的氣體和塵埃會散射太陽光,形成人們可觀測到的彗發(fā)和彗尾[2].彗星形成于太陽系邊緣,相比于內(nèi)太陽系天體,彗星經(jīng)歷的太陽輻射侵蝕更少、成分更為原始,這使得彗星相關(guān)的研究成為了行星科學(xué)前沿課題.由于彗星的彗核質(zhì)量較小,沒有如大行星那樣復(fù)雜的地質(zhì)活動,因此組成彗核的物質(zhì)非常接近原始太陽星云的物質(zhì)[3],對彗核的研究能夠幫助我們理解太陽系其他天體的演化歷史.地面與空間望遠(yuǎn)鏡受限于探測距離的限制,只能以較低的分辨率觀測彗星,而深空探測器能夠抵達(dá)彗星并進(jìn)行近距離觀測,這為我們細(xì)致地研究彗星的彗核提供了條件.

    彗核的形狀通常極不規(guī)則,其中最受關(guān)注的是67P/C-G (Churyumov-Gerasimenko)的“雙頭形(Bi-lobed)”彗核,歐空局的Rosetta任務(wù)獲取了大量的67P/C-G觀測數(shù)據(jù),這使其成為被天文學(xué)家研究得最為深入的彗星.有關(guān)“雙頭形”彗核的形成和演化機(jī)制,許多工作都開展了相關(guān)研究.Massironi等[4]對比了67P/C-G彗核重力場方向與表面分層地貌結(jié)構(gòu)的法向指向之間的差異,認(rèn)為67P/C-G的雙頭形結(jié)構(gòu)可能是由兩個天體緩慢碰撞形成的.Jutzi等[5]進(jìn)行了相對低速的撞擊模擬,原始天體被撞擊后碎裂為兩個主要部分,經(jīng)過緩慢合并形成了雙頭形結(jié)構(gòu).Massironi等[4]與Jutzi等[5]研究的低速碰撞過程發(fā)生在太陽系早期,而Schwartz等[6]做了在太陽系形成任何時期都可以存在的高速撞擊的研究,他們模擬了一個天體受到高速撞擊后完全瓦解為小碎片群,并重新吸積形成雙頭形結(jié)構(gòu)的過程.彗核表面物質(zhì)噴發(fā)對彗核本身有反沖作用,由于形狀不規(guī)則和表面物質(zhì)升華的不均勻性,彗核受到的合外力矩不為零,從而處于加速或減速自轉(zhuǎn)的狀態(tài).對于雙頭形結(jié)構(gòu)的彗星,以67P/C-G為例,Hirabayashi等[7]認(rèn)為由于加速自轉(zhuǎn),它的雙頭連接部分會受到更大的張力,當(dāng)自轉(zhuǎn)足夠快時,可能發(fā)生坍塌或斷裂的現(xiàn)象.除了這種機(jī)械侵蝕效應(yīng),Matonti等[8]考慮了物質(zhì)升華造成的侵蝕效應(yīng),并論證,當(dāng)彗核受到的太陽輻射非常強(qiáng)時,物質(zhì)升華侵蝕效應(yīng)將在彗核形貌的演化中占主導(dǎo)地位.

    Vavilov等[9]提出了由物質(zhì)升華造成的不規(guī)則形狀彗核的演化機(jī)制.彗核表面的光照分布取決于其軌道和自轉(zhuǎn)狀態(tài).光照越強(qiáng),物質(zhì)升華越快,在一個公轉(zhuǎn)周期內(nèi),彗核表面不同區(qū)域由于受光照程度不同,物質(zhì)升華速度不同,彗核的形狀就會變得復(fù)雜.即使是形狀規(guī)則(如橢球形狀)的彗核,在自轉(zhuǎn)軸不穩(wěn)定的情況下,其表面的光照分布會發(fā)生變化,進(jìn)而每個位置的物質(zhì)升華速度也會相應(yīng)地改變.若彗核物質(zhì)分布不均勻,形狀會更加復(fù)雜,最終能夠形成像“雙頭形”這種不規(guī)則的形狀.

    基于在太陽輻射作用下彗核物質(zhì)升華的機(jī)制,Zhao等[10]在之前的工作中提出了小天體形貌演化模型MONET (Mass-loss-driven shape evolution model),認(rèn)為軌道偏心率、自轉(zhuǎn)軸存在傾角以及彗星在近日點和遠(yuǎn)日點處的太陽輻射通量在升華曲線的拐點兩側(cè)是造成彗核南北半球升華分布不對稱的原因.Zhao等[11]把MONET演化模型應(yīng)用在柯伊伯帶小天體Arrokoth的研究中,通過模擬一氧化碳冰的升華過程,很好地解釋了Arrokoth扁平形狀的成因,并且揭示了“初始自轉(zhuǎn)軸指向”這一關(guān)鍵因素,它將極大地影響小天體形狀演化的結(jié)果.

    本文利用MONET演化模型,通過數(shù)值模擬的方法對1顆哈雷型彗星和5顆木星族彗星的水冰升華分布進(jìn)行研究.在第2節(jié)將介紹MONET演化模型,對于本文所研究的彗核的形狀數(shù)據(jù)以及動力學(xué)參數(shù)在這一部分會給出準(zhǔn)確描述.第3節(jié)給出計算結(jié)果,并且討論水冰升華的南北不對稱分布的原因以及彗星自轉(zhuǎn)、進(jìn)動及公轉(zhuǎn)對太陽直射點造成的影響.第4節(jié)做出總結(jié),對比本文所研究的6顆短周期彗星.

    2 模型與方法

    通常,彗星的形狀模型數(shù)據(jù)可以通過3個方法獲得: (1)光變曲線反演法[12],使用地面望遠(yuǎn)鏡對彗星進(jìn)行長時間的測光觀測,通過光變曲線能夠反演出彗核的幾何形狀,但除近日點附近之外,彗星在大部分時間距離地球都較遠(yuǎn),因此精度較差; (2)雷達(dá)反演法[13],使用雷達(dá)發(fā)射信號經(jīng)過彗核反射后接收,同樣能夠得到彗核的形狀,但由于雷達(dá)的靈敏度有限以及大部分時間地球與彗星間的距離較遠(yuǎn),因此分辨率較低; (3)就位探測法[14],人類的深空探測器直接飛行至彗星周圍,環(huán)繞彗核進(jìn)行圖像或高程的數(shù)據(jù)采集,并使用這些數(shù)據(jù)推演出形狀模型,此方法得到的彗核形狀精度較高.

    NASA (National Aeronautics and Space Administration)的PDS(Planetary Data System)小天體節(jié)點1NASA PDS SBN.https://pdssbn.astro.umd.edu/data_sb/target_comets.shtml.保存了已被探測過的部分彗核的幾何形狀,這些模型數(shù)據(jù)是以“三角面元多面體”的形式存在的,構(gòu)成多面體的三角形被稱為“面元”.形狀模型數(shù)據(jù)是在彗核體固坐標(biāo)系下描述的,包括所有頂點的坐標(biāo)以及所有面元所包含的頂點索引值,本文研究的彗核形狀模型均來自該網(wǎng)站.

    2.1 MONET模型

    MONET模型[10]用于研究物質(zhì)升華引起的天體形貌的改變,它包括了氣體升華模型和形狀演化模型兩個部分.氣體升華模型基于熱物理方程和能量守恒,給出了太陽入射能量與水分子生產(chǎn)率之間的關(guān)系,而形狀演化模型是考慮了彗核的自轉(zhuǎn)、公轉(zhuǎn)后,研究物質(zhì)損失導(dǎo)致其形狀變化的物理模型,它決定了天體數(shù)字形狀模型如何發(fā)生形變.在本工作中,我們將MONET模型應(yīng)用在了短周期彗星水冰升華分布的研究中,下面簡要介紹MONET模型,并對該模型中算法的實現(xiàn)給出詳細(xì)的描述.

    2.1.1 氣體升華模型

    彗核表面的能量輸入來自太陽輻射,這些能量的一部分會以熱輻射的形式再次向外輻射,另一部分將會被吸收,使表面揮發(fā)性物質(zhì)發(fā)生升華[10].假設(shè)物質(zhì)為均勻的純水冰,不考慮熱傳導(dǎo)過程,太陽入射能量完全用于熱輻射與水冰升華,能量守恒方程為[10]

    其中,T為彗核表面局部區(qū)域(面元)的溫度,Esol為太陽輻射,Eir(T)為熱輻射的能量,Esub(T)是由水冰揮發(fā)吸收掉的能量.

    我們在之前工作中[10]通過求解上述方程,計算出了水分子生產(chǎn)率與入射太陽輻射通量之間的關(guān)系,如圖1所示,水分子生產(chǎn)率與入射太陽輻射通量之間存在非線性的正相關(guān)關(guān)系,本文基于此關(guān)系開展了短周期彗星物質(zhì)損失的研究.

    圖1 單位面積內(nèi)入射輻射通量與水分子生產(chǎn)率的對應(yīng)關(guān)系Fig.1 The corresponding relation between incident radiation flux and H2O molecular productivity per unit area

    2.1.2 形狀演化模型

    彗核表面單位面積內(nèi)太陽入射的能量與面元的指向有關(guān),因此每個面元的能量輸入是太陽輻射通量在面元法向方向的投影值,當(dāng)點積為負(fù)時,說明該面元背對太陽,值應(yīng)為零.根據(jù)能量守恒定律可得,第i個面元在日心距r處單位時間內(nèi)單位面積內(nèi)接收到的來自太陽的能量Fi(r)為

    單位為W·m-2,其中F⊙=1367 W·m-2為太陽常數(shù),為第i個面元的法向量,為太陽在彗核體固坐標(biāo)系中的單位方向向量.

    由于一些彗核的形狀極為復(fù)雜,在計算表面太陽輻射能量時,我們考慮了由不規(guī)則形狀引起的陰影和自加熱效應(yīng)(Self-heating effects)[10],最終可以得到面元j在單位時間內(nèi)的總能量輸入Ej為

    式中的ESj為太陽的入射能量,若面元被遮擋則該項為0,否則該項等于Fj(r);Nface為彗核形狀包含的面元總數(shù);Ai為面元i的面積,Fij為面元i和j間的可視因子(View Factor)[15];σSB為斯特藩-玻爾茲曼常數(shù);Ti為面元i的熱力學(xué)溫度.

    計算一個公轉(zhuǎn)周期內(nèi)累積的能量本質(zhì)上是一個對時間積分的過程,開普勒第二定律決定了彗星的公轉(zhuǎn)角速度在近日點附近比遠(yuǎn)日點附近更快,所以不能把整個公轉(zhuǎn)按照時間均勻分割,否則彗星在近日點附近會采樣不足,無法準(zhǔn)確地計算光照情況.考慮到這一點,我們選取的積分變量為偏近點角而非平近點角,間接地計算角度增量對應(yīng)的時間增量[10].開普勒方程給出了偏近點角(與位置有關(guān))和平近點角(與時間有關(guān))的關(guān)系[16],因此我們把整個軌道按照偏近點角在空間上均勻分割,然后計算每一小段的能量并積分.

    把周期為Porbit的公轉(zhuǎn)按照偏近點角平均分割為Norbit份,用E'(為了和能量E區(qū)分開)表示偏近點角,由開普勒方程以及平角速度的定義可知[10,16]:

    在上式中,k ∈[0,Norbit),E'0為初始位置偏近點角,E'k和Mk分別為第k個軌道片段起始處的偏近點角與平近點角,e為彗星軌道的偏心率,ΔMk為第k個軌道片段對應(yīng)的平近點角增量,Δtk為第k個軌道片段對應(yīng)的持續(xù)時間(用軌道片段所占平近點角增量除以平角速度求得),n為彗星的平角速度.

    我們做了這樣一個近似[10]: 在劃分后的每個軌道片段上,彗核的每個自轉(zhuǎn)周期都近似為完全相同,即這個軌道片段內(nèi)的運動包含了多次完全相同的自轉(zhuǎn)過程.因此面元i在整個公轉(zhuǎn)周期內(nèi)累積的太陽輻射能量Eiorbit為[10]:

    為其在第k個軌道片段內(nèi)的一個自轉(zhuǎn)周期內(nèi)的太陽輻射能量,Pspin為自轉(zhuǎn)周期.

    在彗核體固坐標(biāo)系下,太陽做周日視運動,其位置隨時間變化,周期等于彗核自轉(zhuǎn)周期,面元i在單位時間內(nèi)的能量輸入Ei(t)是時間的函數(shù).把彗核的一個自轉(zhuǎn)周期分割為Nspin個片段,分別計算每個片段起始位置處的能量分布,對時間進(jìn)行積分,便可得到面元i在一整個自轉(zhuǎn)周期內(nèi)的太陽輻射能量Eispin(單位面積內(nèi))[10]:

    其中,l ∈[1,Nspin],tl為第l個片段所對應(yīng)的起始時刻,積分時間微元為Δtl=Pspin/Nspin.

    以上對太陽輻射能量的推導(dǎo),同樣適用于侵蝕深度的計算.我們使用MONET氣體升華模型中的數(shù)據(jù),用插值的方法來得到不同太陽輻射通量所對應(yīng)的水分子數(shù)量生產(chǎn)率.取水分子質(zhì)量為2.99×10-26kg,水冰密度為ρ=920 kg·m-3,計算出面元i在任意時間段內(nèi)單位面積上損失掉的質(zhì)量Qi,單位為kg·m-2.如圖2所示,為第i個面元的單位法向量,把三棱柱體積看作被侵蝕掉的水冰所占據(jù)的體積,定義單位時間內(nèi)的侵蝕深度為[10]

    圖2 單位時間內(nèi)侵蝕深度的定義Fig.2 The definition of erosion depth in unit time

    在計算出太陽的瞬時入射能量Ei(t)時,可一并計算出瞬時的水冰侵蝕深度Diunit(t),前文所述的積分方法皆適用于Diunit(t),于是有[10]

    其中,為第i個面元在一個公轉(zhuǎn)周期內(nèi)累積的侵蝕深度,為第i個面元在第k個軌道片段內(nèi)每個自轉(zhuǎn)周期內(nèi)累積的侵蝕深度,為第i個面元在任意一個自轉(zhuǎn)周期內(nèi)累積的侵蝕深度.

    2.2 研究對象

    2.2.1 形狀模型

    本文選取了6顆短周期彗星作為研究對象,包括了1P/Halley、9P/Tempel 1、19P/Borrelly、67P/C-G、81P/Wild 2和103P/Hartley 2,后文正文及圖表中將使用彗星簡稱(1P、9P、19P、67P、81P、103P).這幾顆彗星的近日點均位于“冰”線以內(nèi),強(qiáng)烈的太陽輻射使得水冰升華的侵蝕效應(yīng)顯著,因此,相比于長周期彗星,在同樣的時間尺度下,短周期彗星形狀受到水冰升華侵蝕影響更加顯著.

    如圖3所示,(a)至(e)分別表示1P/Halley、9P/Tempel 1、67P/C-G、81P/Wild 2和103P/Hartley 2的形狀模型[17-21],圖中的坐標(biāo)系為每個彗核各自的體固坐標(biāo)系XY Z,從左到右的4列圖片,分別是每個彗核的正視圖、右視圖、俯視圖與仰視圖.由于Deep Space任務(wù)沒有獲取到19P/Borrelly的形狀數(shù)據(jù),從僅有的觀測圖像[22]來看,它與103P/Hartley 2的紡錘形狀非常相似,所以本文使用經(jīng)過縮放后的103P/Hartley 2的形狀數(shù)據(jù)來研究19P/Borrelly.我們統(tǒng)計了103P/Hartley 2形狀數(shù)據(jù)所有頂點坐標(biāo)后,得到其長軸約為2.34 km,Rayman[23]指出19P/Borrelly的長軸為8 km,因此本文把103P/Hartley 2形狀數(shù)據(jù)放大為3.42倍后,作為19P/Borrelly的形狀開展研究.

    圖3 彗星的形狀模型[17-21]Fig.3 Shape model of comets[17-21]

    2.2.2 動力學(xué)參數(shù)

    本文研究的短周期彗星的軌道根數(shù)是從IAU(International Astronomical Union) Minor Planet Center2IAU Minor Planet Center: https://minorplanetcenter.net/.獲取的,其軌道根數(shù)等信息已在表1中給出,各列數(shù)據(jù)依次為軌道根數(shù)所對應(yīng)的歷元(儒略日)、近點角距ω、升交點赤經(jīng)Ω、軌道傾角inc、軌道偏心率e、軌道半長軸a、平近點角M以及軌道周期Porbit.每顆彗星的彗核自轉(zhuǎn)狀態(tài)參數(shù)及形狀參數(shù)如表2所示,其中Ra和Dec分別表示彗核角動量方向指向的赤經(jīng)、赤緯,Pspin為彗核的自轉(zhuǎn)周期,Pφ為彗核的進(jìn)動周期,θ為彗核自轉(zhuǎn)軸與角動量軸(進(jìn)動軸)之間的夾角(后文簡稱為進(jìn)動角),“-”表示沒有此類數(shù)據(jù).表2中最后4列是本文計算所得數(shù)據(jù): 第8-9列分別為彗核角動量在其軌道面坐標(biāo)系(以軌道面法向為ZorbitPlane軸,以太陽-近日點方向為XorbitPlane軸)下的指向,其中的α(Obliquity)是指彗核角動量方向與軌道面法向間的夾角,β(Longitude)是指彗核角動量方向在軌道面的投影與太陽-近日點方向之間的夾角; 最后兩列是每顆彗星的彗核形狀參數(shù)統(tǒng)計,包括組成彗核形狀的面元數(shù)量和面元的平均面積.

    表1 彗星的軌道根數(shù)Table 1 The orbital elements of comets

    表2 彗核的旋轉(zhuǎn)以及形狀參數(shù)Table 2 The rotation and shape parameters of cometary nucleus

    對于自轉(zhuǎn)軸發(fā)生進(jìn)動的彗核,我們?nèi)∑溥M(jìn)動周期與自轉(zhuǎn)周期的最小公倍數(shù)的近似值(自轉(zhuǎn)周期、進(jìn)動周期各自的整數(shù)倍與該近似值的相對誤差均小于5%)作為“等效自轉(zhuǎn)周期Pspin”,在這個等效的自轉(zhuǎn)周期內(nèi),彗核完成了整數(shù)次的自轉(zhuǎn)與進(jìn)動.本文取1P/Halley、67P/C-G和103P/Hartley 2的等效自轉(zhuǎn)周期分別為7.2 d、10.5 d 和55.2 h.彗星19P/Borrelly與81P/Wild 2的自轉(zhuǎn)周期未知,由于它們不存在進(jìn)動,所以自轉(zhuǎn)周期的大小并不會影響到整個公轉(zhuǎn)周期內(nèi)的計算,本文假設(shè)其與67P/C-G的自轉(zhuǎn)周期相同,并進(jìn)行計算.

    3 結(jié)果及討論

    本文計算了表1所列出的短周期彗星彗核在一個公轉(zhuǎn)周期內(nèi)累積的太陽輻射能量以及侵蝕深度的空間分布,其中兩個數(shù)值模擬參數(shù)分別為Norbit=360,Nspin=100.圖4給出了面元侵蝕深度的統(tǒng)計直方圖,橫坐標(biāo)為侵蝕深度,單位為m,縱坐標(biāo)為面元數(shù).圖5給出了本文計算出的每個彗核的太陽輻射能量與侵蝕深度的表面分布情況,坐標(biāo)系為每個彗核各自的體固坐標(biāo)系XY Z,圖左半部分是太陽輻射能量在彗核表面的分布(單位為J),從藍(lán)色到黃色代表能量由低到高,右半部分為侵蝕深度的空間分布(單位為m),從藍(lán)色到紅色代表侵蝕深度由小到大.每幅小圖從左至右分別為彗核的北半球(沿Z軸負(fù)方向看)、赤道區(qū)域(沿Y軸負(fù)方向看)以及南半球(沿Z軸正方向看)的視圖.圖5中67P/C-G的侵蝕深度最大值與Keller等[27]的工作中給出的結(jié)果處于同一數(shù)量級,且相差較小,驗證了本文計算結(jié)果的可靠性.

    圖5 6顆彗星彗核的太陽輻射能量與侵蝕深度的分布Fig.5 Distribution of solar radiation energy and erosion depth of cometary nucleus of 6 comets

    3.1 侵蝕深度的頻數(shù)統(tǒng)計

    1P/Halley、9P/Tempel 1和19P/Borrelly的侵蝕深度直方圖表明,在一個公轉(zhuǎn)周期內(nèi),這些彗核的大多數(shù)面元都具有較高的侵蝕深度,半數(shù)以上的面元都集中在直方圖右側(cè),這說明彗核的大部分區(qū)域?qū)S著時間推移而向內(nèi)部快速收縮.67P/C-G的直方圖相對平均一些,侵蝕量由小到大均有分布,快速侵蝕的區(qū)域只占一小部分,這會導(dǎo)致67P/C-G的形狀并不像前述彗核一樣大范圍收縮,而是只有少數(shù)局部收縮較快,67P/C-G的形狀將會變得更加復(fù)雜.81P/Wild 2和103P/Hartley 2的直方圖呈現(xiàn)出接近高斯分布的情況,雖然有很多面元處于中間位置,但圖中81P/Wild 2的右半部分面元數(shù)量多于左半部分,彗核的形狀仍然會在局部快速收縮.103P/Hartley 2與前述彗核都不同,它的絕大部分面元都處在中間的位置,兩側(cè)面元的數(shù)量遠(yuǎn)小于中間面元的數(shù)量,另外,從坐標(biāo)軸刻度上還能發(fā)現(xiàn),103P/Hartley 2的侵蝕深度最大值與最小值差別很小,處于同一個數(shù)量級.103P/Hartley 2的這些特點表明它的侵蝕深度在空間上的分布非常均勻,形狀具有均勻縮小的演化趨勢.

    3.2 侵蝕深度的南北半球非對稱性分布

    Zhao等[10]在前期的工作中還計算了在不同軌道半長軸與偏心率的情況下,彗星在遠(yuǎn)日點與近日點處的水分子生產(chǎn)率之比,如圖6所示,比值越小,彗星在遠(yuǎn)日點和近日點處的水分子生產(chǎn)率差異越大.若彗星在遠(yuǎn)日點處的太陽入射能量低于圖1中的拐點,則彗星遠(yuǎn)日點和近日點處水分子生產(chǎn)率差距顯著,反之若遠(yuǎn)日點處的太陽入射能量高于拐點,則遠(yuǎn)日點與近日點處的水分子生產(chǎn)率差距不明顯.該圖中深藍(lán)色點線代表此位置軌道上的彗星在遠(yuǎn)日點處的水分子生產(chǎn)率剛好位于圖1所示的拐點處,意味著太陽入射能量-水分子生產(chǎn)率熱物理關(guān)系的非線性特點,加劇了該軌道上近日點和遠(yuǎn)日點之間的水分子生產(chǎn)率差異,該軌道上天體的水冰升華會存在南北半球非對稱分布的現(xiàn)象.紅色圖標(biāo)為幾個木星族彗星按照半長軸和偏心率經(jīng)過理論計算后在圖中所處的位置,包含了本文的研究對象

    圖6 遠(yuǎn)日點與近日點處的水分子生產(chǎn)率之比Fig.6 The ratio of H2O molecular productivity at aphelion to perihelion

    9P/Tempel 1、19P/Borrelly、67P/C-G、81P/Wild 2和103P/Hartley 2,并且本文所有研究對象都處在白色分界線的右側(cè)(彗星1P/Halley具有足夠大的半長軸與偏心率,顯然它會在該圖的右上方位置).因此,本文的所有研究對象的軌道都滿足產(chǎn)生南北半球非對稱侵蝕現(xiàn)象的條件.

    彗星1P/Halley的能量分布特點是赤道地區(qū)入射能量較多,兩極地區(qū)入射能量較少,入射能量最大值和最小值相差了一個數(shù)量級.對于侵蝕深度的分布情況,同樣也是赤道地區(qū)被侵蝕得較多兩極地區(qū)較少,但是和入射能量分布不同,侵蝕深度的南北半球分布呈現(xiàn)出非對稱現(xiàn)象,南半球被侵蝕得比北半球更多.

    9P/Tempel 1的太陽輻射能量分布與彗星1P/Halley類似,呈現(xiàn)出赤道地區(qū)入射能量較多而兩極地區(qū)入射能量較少的特點.9P/Tempel 1侵蝕深度的分布同樣具有南北半球非對稱的現(xiàn)象,然而和1P/Halley不同的是,其北半球被侵蝕得比南半球更多.

    19P/Borrelly的太陽輻射能量分布情況與前述彗星相同,赤道地區(qū)接收到的太陽輻射能量多于兩極地區(qū).侵蝕深度的南北半球不對稱性非常明顯,面朝南的面元被侵蝕的平均深度大于面朝北的面元.

    67P/C-G的太陽輻射能量與侵蝕深度分布的特點和19P/Borrelly相似.赤道地區(qū)接收到的太陽輻射能量高于兩極地區(qū),南半球被侵蝕得比北半球更多.67P/C-G凹陷處的特殊地形,即雙頭連接的脖子部分,無論是太陽輻射能量還是侵蝕深度,都比周圍鄰近區(qū)域(同一半球)的少,這是因為受到陰影效應(yīng)的影響,凹陷處地形只在較短的時間內(nèi)才能受太陽照射.

    彗星81P/Wild 2沒有完整的形狀數(shù)據(jù).Stardust探測器飛掠這顆彗星時,太陽只照亮了它一部分表面,探測器無法獲取處于黑暗部分的形狀數(shù)據(jù),NASA的小天體節(jié)點用不完整的形狀數(shù)據(jù)混合了三軸橢球體,給出了一個閉合的形狀模型,本文使用了該形狀模型進(jìn)行計算.這顆彗星的太陽輻射能量與侵蝕深度的分布特點,和1P/Halley、19P/Borrelly、67P/C-G的特點相同,此處不再贅述.

    彗星103P/Hartley 2的情況與前面所述的彗星有所不同,它最顯著的特點就是分布非常均勻.對于太陽輻射能量,赤道和兩極地區(qū)并沒有明顯差別,都處于同一個數(shù)量級,兩極地區(qū)的值比赤道地區(qū)略高.對于侵蝕深度同樣如此,分布較為均勻,南極區(qū)域比其他區(qū)域侵蝕得略多,相差不大.這些特點可能與非主軸自轉(zhuǎn)有關(guān),我們將在3.3節(jié)詳細(xì)討論自轉(zhuǎn)、進(jìn)動和公轉(zhuǎn)對彗核侵蝕深度分布的影響.

    3.3 太陽直射點位置的定量計算結(jié)果

    彗核的自轉(zhuǎn)和進(jìn)動使得太陽直射點掃過的軌跡覆蓋了彗核表面由兩個緯線圍成的帶狀區(qū)域,當(dāng)彗核自轉(zhuǎn)軸不發(fā)生進(jìn)動時,則太陽直射點軌跡覆蓋區(qū)域會退化成一條緯線.彗星公轉(zhuǎn)會影響彗核角動量軸與太陽之間的幾何相對關(guān)系,進(jìn)而使太陽直射點軌跡所覆蓋的帶狀區(qū)域的兩條邊界緯線隨時間改變位置.

    我們計算了本文所研究的6顆彗星彗核表面太陽直射點的覆蓋情況,包括在軌道上真近點角分別為β和(β±π)這兩個關(guān)鍵位置處(原因是,在這兩個位置處太陽-彗核連線與彗核角動量的夾角會分別取到最大值和最小值,對應(yīng)著彗星公轉(zhuǎn)所引起的太陽直射點運動軌跡的緯度邊界)太陽直射點軌跡覆蓋區(qū)的緯度范圍以及在整個公轉(zhuǎn)周期內(nèi)太陽直射點軌跡覆蓋的緯度范圍,如表3所示.彗核角動量的方向在其軌道面坐標(biāo)系內(nèi)表示為(α,β).Θ0是太陽-彗核連線與彗核角動量之間的夾角(在彗核體固坐標(biāo)系下該角度Θ0即為太陽的天頂距,Θ0∈[0,π]),彗星公轉(zhuǎn)會使Θ0隨著真近點角f的變化而改變,當(dāng)f等于β時,Θ0會取到最大值(或最小值); 當(dāng)f等于(β±π)時,Θ0會取到最小值(或最大值).Lat是太陽直射點的緯度范圍.r⊙是彗星的日心距,單位為au.

    表3 彗核表面太陽直射點緯度Table 3 Latitude of solar point on surface of comets

    對于本文所研究的彗星,我們在其真近點角為β和β±π的位置分別計算了Θ0、直射點緯度范圍Lat和日心距r⊙,表3很好地解釋了前文數(shù)值模擬的計算結(jié)果.例如,1P/Halley、19P/Borrelly、67P/C-G和81P/Wild 2都有一個共同的特點: 在f=β時,太陽直射點在南半球,且日心距較小,而相隔π相位的位置太陽直射點在北半球,且日心距更大.由于輻射通量與日心距成平方反比以及太陽輻射-水分子生產(chǎn)率的非線性關(guān)系,彗核在一個公轉(zhuǎn)周期內(nèi),南北半球的侵蝕深度會表現(xiàn)出明顯的差異,即南半球被侵蝕得比北半球更多.9P/Tempel 1與前面4顆彗星相反,它的日心距在太陽直射北半球時比直射南半球時更小,導(dǎo)致北半球被侵蝕得比南半球更多,這從圖5也能看出.另外,在整個公轉(zhuǎn)周期內(nèi),9P/Tempel 1的太陽直射點幾乎都在赤道附近移動,因此赤道區(qū)域受到的侵蝕效應(yīng)比南北半球都要更明顯.103P/Hartley 2彗星則有點特殊,如表3所示,它擁有非常大的進(jìn)動角θ,這導(dǎo)致無論彗星位于f=β還是f=β±π處,在自轉(zhuǎn)和進(jìn)動的作用下,第一,太陽直射點總能穿行在南北兩個半球之間; 第二,直射點軌跡在南北兩個半球的覆蓋范圍相差不大(直射點軌跡在南北半球的邊界緯度其絕對值相差較小).這兩個特點最終使得103P/Hartley 2這顆彗星的彗核在一整個公轉(zhuǎn)周期內(nèi)累積的侵蝕深度并沒有明顯的南北非對稱性,反而十分均勻.1P/Halley彗星也擁有較大的進(jìn)動角θ,其太陽直射點也能夠穿行在南北兩個半球,但是由于其軌道偏心率非常大,這使彗星在f=β和f=β±π時的日心距相差非常巨大,已經(jīng)有了數(shù)量級差距,所以微小的南北半球非對稱性被放大,體現(xiàn)在圖5中就是南北半球不對稱的侵蝕深度分布.由于1P/Halley在表格中f取最值的兩個位置的太陽直射點都同時覆蓋了較大的緯度范圍,形成了緯度重疊區(qū),重疊區(qū)內(nèi)的表面永遠(yuǎn)被太陽直射,所以1P/Halley除兩極區(qū)域以外的大部分區(qū)域均受到明顯侵蝕.

    如前文模型與方法部分所述,本文的工作建立在3個假設(shè)之上: 第一,組成彗核的物質(zhì)是各向同性的純水冰,物質(zhì)分布是均勻的.若彗核含有其他非水冰成分,則只會影響在相同條件下的水分子生產(chǎn)率的大小,不會對本文所研究的水冰升華的空間分布有影響.本文假設(shè)物質(zhì)分布是均勻的,若彗核的物質(zhì)分布不均勻,則彗核表面不同區(qū)域的水分子生產(chǎn)率會有差異,最終影響水冰升華的空間分布,所以在未來的工作中根據(jù)彗核探測數(shù)據(jù),引入更復(fù)雜的非均勻物質(zhì)分布模型,將會得到更為真實的結(jié)果.第二,忽略彗核表面物質(zhì)的熱慣量影響,不考慮熱傳導(dǎo)過程,彗核表面的太陽入射能量等于熱輻射消耗的能量與水冰升華消耗的能量之和.第三,彗星的軌道穩(wěn)定、彗核的自轉(zhuǎn)及進(jìn)動狀態(tài)在一個軌道周期內(nèi)不發(fā)生改變.

    4 結(jié)論

    在MONET模型的基礎(chǔ)上,本文開展了對彗核進(jìn)動效應(yīng)的模擬,研究短周期彗星1P/Halley、9P/Tempel 1、19P/Borrelly、67P/C-G、81P/Wild 2和103P/Hartley 2水冰升華的分布情況.短時間內(nèi)(彗核自轉(zhuǎn)、進(jìn)動周期的時間尺度),彗核自轉(zhuǎn)會使太陽直射點在相同緯度上移動,進(jìn)動會使太陽直射點運動的覆蓋范圍擴(kuò)展為兩個緯線包圍的帶狀區(qū)域,區(qū)域大小受Θ0(彗星-太陽連線與彗核角動量間的夾角)和進(jìn)動角θ(彗核自轉(zhuǎn)軸與彗核角動量間的夾角)影響.彗星公轉(zhuǎn)會使Θ0發(fā)生變化,造成太陽直射點運動的覆蓋區(qū)域改變大小并發(fā)生移動,彗核自轉(zhuǎn)角動量與軌道角動量的夾角α決定了Θ0的變化上下限,彗核自轉(zhuǎn)角動量在軌道面的投影與太陽-近日點方向之間的夾角β決定了Θ0取得最大和最小值時彗星的真近點角.

    橢圓軌道偏心率的存在使得彗核表面的同一個面元在近日點附近和遠(yuǎn)日點附近單位時間內(nèi)接收到的太陽輻射能量存在差異,太陽輻射與水分子生產(chǎn)率間的非線性關(guān)系會使這兩個位置處的水分子生產(chǎn)率呈現(xiàn)出巨大差異,而水分子生產(chǎn)率的差異會隨時間持續(xù)累積,這是水冰侵蝕深度分布出現(xiàn)南北半球不對稱現(xiàn)象的根本原因.角度β決定了被侵蝕更多的是南半球還是北半球(若β=則侵蝕深度應(yīng)是南北半球?qū)ΨQ分布的).進(jìn)動使得彗核表面的太陽直射點在短時間內(nèi)(自轉(zhuǎn)、進(jìn)動周期的時間尺度)覆蓋的范圍擴(kuò)大,進(jìn)動角越接近太陽直射點在南北半球的覆蓋范圍就會越對稱,最終會減弱甚至抵消侵蝕深度分布的南北半球不對稱性.

    彗星1P/Halley是哈雷型彗星,其公轉(zhuǎn)周期和軌道半長軸均明顯高于其他5顆木星族彗星,但計算結(jié)果表明它們的侵蝕深度在同一量級,相差較小.彗星1P/Halley的大偏心率軌道使其在近日點處的日心距比其他5顆木星族彗星都更小,由于太陽輻射-水分子生產(chǎn)率的非線性關(guān)系,近日點附近更強(qiáng)的太陽輻射導(dǎo)致了更高的水分子生產(chǎn)率.雖然1P/Halley在軌道上的大部分時間距離太陽都很遠(yuǎn),但其在近日點附近累積的水冰侵蝕深度足以達(dá)到和另外5顆木星族彗星相同的數(shù)量級.1P/Halley和5顆木星族彗星之間的對比表明,影響水冰侵蝕深度數(shù)量級的主導(dǎo)因素是彗星的軌道.

    木星族彗星9P/Tempel 1、19P/Borrelly、67P/C-G、81P/Wild 2和103P/Hartley 2的水冰侵蝕深度最大值相差較小,但是它們在侵蝕深度的分布上有很大差異.自轉(zhuǎn)軸的初始指向以及進(jìn)動的存在,使彗核表面的水冰侵蝕深度出現(xiàn)了南北半球非對稱分布的現(xiàn)象,當(dāng)進(jìn)動角足夠接近南北半球分布的非對稱性會被減弱甚至抵消.彗核的自轉(zhuǎn)和進(jìn)動決定了水冰侵蝕深度在空間上的分布情況.

    猜你喜歡
    水冰進(jìn)動面元
    隨機(jī)粗糙面散射中遮蔽效應(yīng)算法的改進(jìn)
    嚴(yán)水冰
    月球兩極確有水冰
    月球極區(qū)
    基于窄帶雷達(dá)網(wǎng)的彈道目標(biāo)三維進(jìn)動特征提取
    基于改進(jìn)Gordon方程的RCS快速算法
    進(jìn)動錐體目標(biāo)平動補(bǔ)償及微多普勒提取
    面元細(xì)分觀測系統(tǒng)應(yīng)用分析
    化工管理(2014年14期)2014-08-15 00:51:32
    基于雷達(dá)距離像的錐體目標(biāo)進(jìn)動參數(shù)估計方法
    時間距離像消隱時彈頭進(jìn)動與結(jié)構(gòu)參數(shù)估計
    精品午夜福利视频在线观看一区| 久久中文看片网| 亚洲精品美女久久av网站| 母亲3免费完整高清在线观看| 久久这里只有精品19| 九色国产91popny在线| 国产精品乱码一区二三区的特点| 无遮挡黄片免费观看| 久久香蕉国产精品| 亚洲欧美日韩高清专用| 精品久久久久久久久久免费视频| 亚洲 欧美 日韩 在线 免费| 国产精品一及| 久久精品成人免费网站| 欧美一级毛片孕妇| 国产精品国产高清国产av| 嫩草影视91久久| 男女床上黄色一级片免费看| 日本一二三区视频观看| 久久热在线av| 亚洲中文字幕一区二区三区有码在线看 | 亚洲午夜精品一区,二区,三区| 亚洲av成人不卡在线观看播放网| 成在线人永久免费视频| 久久久国产成人精品二区| 人人妻人人看人人澡| 亚洲成人久久性| 美女高潮喷水抽搐中文字幕| 亚洲一卡2卡3卡4卡5卡精品中文| 精品国产乱子伦一区二区三区| 一区二区三区激情视频| 一二三四社区在线视频社区8| 久久久久九九精品影院| 人妻丰满熟妇av一区二区三区| 日本黄色视频三级网站网址| 久久婷婷人人爽人人干人人爱| 久久精品国产亚洲av高清一级| 国产97色在线日韩免费| 国内精品一区二区在线观看| 日本五十路高清| 亚洲av电影在线进入| 中亚洲国语对白在线视频| 精品久久久久久成人av| 在线观看66精品国产| 日本 av在线| 国产亚洲精品久久久久5区| 精品无人区乱码1区二区| 波多野结衣巨乳人妻| 久久久久久九九精品二区国产 | 久久久久久久久中文| 国产高清视频在线播放一区| 午夜福利免费观看在线| 欧美日韩福利视频一区二区| 99久久99久久久精品蜜桃| 身体一侧抽搐| a级毛片a级免费在线| bbb黄色大片| 麻豆一二三区av精品| 国产精品综合久久久久久久免费| 亚洲中文av在线| 可以在线观看的亚洲视频| 国产激情欧美一区二区| 亚洲精品在线观看二区| 欧美性猛交黑人性爽| 午夜老司机福利片| 欧美 亚洲 国产 日韩一| 男女午夜视频在线观看| 在线观看美女被高潮喷水网站 | 老熟妇仑乱视频hdxx| 床上黄色一级片| 日本一区二区免费在线视频| 国产亚洲精品一区二区www| 好男人电影高清在线观看| 亚洲中文av在线| 午夜激情福利司机影院| netflix在线观看网站| av天堂在线播放| 99国产综合亚洲精品| 久久久久性生活片| 欧美绝顶高潮抽搐喷水| 国产高清激情床上av| 日本 欧美在线| 国内久久婷婷六月综合欲色啪| 变态另类成人亚洲欧美熟女| 老汉色av国产亚洲站长工具| 两人在一起打扑克的视频| 国产成+人综合+亚洲专区| 日韩成人在线观看一区二区三区| 日韩av在线大香蕉| 淫妇啪啪啪对白视频| 黄色 视频免费看| 亚洲精品粉嫩美女一区| 免费在线观看成人毛片| 99国产精品一区二区三区| 美女免费视频网站| 久久久久九九精品影院| 一区二区三区国产精品乱码| 色在线成人网| 一个人观看的视频www高清免费观看 | 国产亚洲欧美98| 一本大道久久a久久精品| 亚洲成人精品中文字幕电影| 免费看十八禁软件| 国产精品野战在线观看| 18禁黄网站禁片免费观看直播| 精品免费久久久久久久清纯| 一本精品99久久精品77| 每晚都被弄得嗷嗷叫到高潮| 久久中文字幕人妻熟女| 黄色丝袜av网址大全| 国产精品亚洲一级av第二区| videosex国产| 亚洲avbb在线观看| 色哟哟哟哟哟哟| 色综合站精品国产| 首页视频小说图片口味搜索| 老司机在亚洲福利影院| 国产精品九九99| xxxwww97欧美| 在线观看66精品国产| 国产成年人精品一区二区| 国产欧美日韩一区二区精品| 好看av亚洲va欧美ⅴa在| 一级a爱片免费观看的视频| 国产精品 国内视频| 成人18禁高潮啪啪吃奶动态图| 久久婷婷成人综合色麻豆| 久久久久久国产a免费观看| 97超级碰碰碰精品色视频在线观看| 特级一级黄色大片| 午夜免费观看网址| 老司机午夜十八禁免费视频| 久热爱精品视频在线9| 欧美黑人欧美精品刺激| 首页视频小说图片口味搜索| 亚洲国产高清在线一区二区三| 久久久水蜜桃国产精品网| 精品一区二区三区四区五区乱码| 一本大道久久a久久精品| 亚洲片人在线观看| 久久性视频一级片| 两性午夜刺激爽爽歪歪视频在线观看 | 禁无遮挡网站| 人人妻人人看人人澡| 丁香欧美五月| 国产成人影院久久av| 国产99久久九九免费精品| www.自偷自拍.com| 中国美女看黄片| 法律面前人人平等表现在哪些方面| 欧美 亚洲 国产 日韩一| 人人妻人人澡欧美一区二区| 亚洲熟妇中文字幕五十中出| 高清毛片免费观看视频网站| 日本三级黄在线观看| 国产精品乱码一区二三区的特点| 亚洲自拍偷在线| 天天添夜夜摸| 欧美人与性动交α欧美精品济南到| 成人18禁高潮啪啪吃奶动态图| 曰老女人黄片| 啦啦啦免费观看视频1| 亚洲av第一区精品v没综合| 天堂影院成人在线观看| 亚洲精品在线观看二区| 欧美在线一区亚洲| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 制服诱惑二区| 白带黄色成豆腐渣| 男女床上黄色一级片免费看| 97超级碰碰碰精品色视频在线观看| 激情在线观看视频在线高清| 岛国视频午夜一区免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品永久免费网站| 777久久人妻少妇嫩草av网站| 美女高潮喷水抽搐中文字幕| 午夜福利成人在线免费观看| 小说图片视频综合网站| 黄色视频,在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 国产午夜福利久久久久久| 嫩草影院精品99| 首页视频小说图片口味搜索| 麻豆成人av在线观看| 麻豆国产97在线/欧美 | 国产又黄又爽又无遮挡在线| 久久久久久免费高清国产稀缺| 亚洲国产欧美人成| 中文在线观看免费www的网站 | 亚洲 国产 在线| 久久久久久大精品| 成人国语在线视频| 成人国产综合亚洲| 五月玫瑰六月丁香| 久久午夜亚洲精品久久| 熟女少妇亚洲综合色aaa.| 在线观看www视频免费| 国产精品自产拍在线观看55亚洲| 久久久久国产一级毛片高清牌| 9191精品国产免费久久| 一进一出好大好爽视频| 国产成人av教育| 国产精品一区二区三区四区免费观看 | cao死你这个sao货| 9191精品国产免费久久| 国产不卡一卡二| 99精品久久久久人妻精品| 中文字幕av在线有码专区| 少妇粗大呻吟视频| 久久天躁狠狠躁夜夜2o2o| 一级毛片高清免费大全| 久久精品夜夜夜夜夜久久蜜豆 | 最新美女视频免费是黄的| av天堂在线播放| 日韩精品免费视频一区二区三区| 国产精华一区二区三区| 不卡av一区二区三区| 2021天堂中文幕一二区在线观| 久久久久国内视频| 欧美高清成人免费视频www| 亚洲成人精品中文字幕电影| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品免费视频内射| 在线观看美女被高潮喷水网站 | 国产精品久久久久久久电影 | 动漫黄色视频在线观看| 亚洲国产精品成人综合色| 欧美中文日本在线观看视频| 日本免费一区二区三区高清不卡| av超薄肉色丝袜交足视频| 中文字幕av在线有码专区| 国产精品久久久人人做人人爽| 99riav亚洲国产免费| 最近最新中文字幕大全电影3| 一边摸一边做爽爽视频免费| 久久久久久久精品吃奶| 亚洲在线自拍视频| 国产精品久久久人人做人人爽| 在线观看日韩欧美| 国产免费男女视频| 深夜精品福利| svipshipincom国产片| 神马国产精品三级电影在线观看 | 成人永久免费在线观看视频| 国产欧美日韩一区二区三| 国产亚洲欧美98| 在线观看www视频免费| 18禁美女被吸乳视频| 国产激情欧美一区二区| 2021天堂中文幕一二区在线观| 最近视频中文字幕2019在线8| 久久热在线av| 午夜福利18| 床上黄色一级片| 欧美中文日本在线观看视频| 亚洲av成人不卡在线观看播放网| 亚洲精华国产精华精| 搡老熟女国产l中国老女人| 日韩欧美一区二区三区在线观看| 真人一进一出gif抽搐免费| 国产伦在线观看视频一区| 亚洲av成人av| xxx96com| 在线十欧美十亚洲十日本专区| av有码第一页| 高清在线国产一区| 啦啦啦免费观看视频1| 国产成人av教育| 夜夜看夜夜爽夜夜摸| 国产精品av久久久久免费| 国产av一区在线观看免费| 久久人人精品亚洲av| 成人午夜高清在线视频| 香蕉久久夜色| 琪琪午夜伦伦电影理论片6080| 国内久久婷婷六月综合欲色啪| 日韩高清综合在线| 亚洲av五月六月丁香网| 黄色视频不卡| 久久草成人影院| 男人舔女人的私密视频| 男女视频在线观看网站免费 | av有码第一页| 国产成人欧美在线观看| 精品国产亚洲在线| 亚洲七黄色美女视频| 在线国产一区二区在线| 草草在线视频免费看| 久久香蕉激情| 国产精品一区二区精品视频观看| 岛国在线免费视频观看| 可以在线观看毛片的网站| 亚洲成人精品中文字幕电影| 欧美黑人巨大hd| 国产99久久九九免费精品| 午夜免费观看网址| bbb黄色大片| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 国产激情久久老熟女| 精品一区二区三区视频在线观看免费| 欧美性猛交黑人性爽| 日本在线视频免费播放| 国产三级在线视频| 三级毛片av免费| 麻豆久久精品国产亚洲av| 国产亚洲精品久久久久5区| 可以在线观看毛片的网站| 美女午夜性视频免费| tocl精华| 人妻夜夜爽99麻豆av| 久久人妻福利社区极品人妻图片| 精品高清国产在线一区| 丰满人妻一区二区三区视频av | 在线视频色国产色| 两性夫妻黄色片| 麻豆成人av在线观看| 99久久精品热视频| 丰满人妻熟妇乱又伦精品不卡| 18禁黄网站禁片免费观看直播| 久久香蕉激情| 亚洲一区二区三区色噜噜| 一个人免费在线观看电影 | 久久精品亚洲精品国产色婷小说| 精品乱码久久久久久99久播| 国产成人aa在线观看| 色哟哟哟哟哟哟| 日日干狠狠操夜夜爽| 国产成+人综合+亚洲专区| 草草在线视频免费看| 两性夫妻黄色片| 熟女电影av网| 国产成人精品无人区| 国产伦在线观看视频一区| 岛国视频午夜一区免费看| 黄色丝袜av网址大全| 2021天堂中文幕一二区在线观| 国产成人影院久久av| av免费在线观看网站| 露出奶头的视频| 国产精品永久免费网站| 激情在线观看视频在线高清| 91国产中文字幕| 日日摸夜夜添夜夜添小说| 99久久精品热视频| 久99久视频精品免费| 成人av一区二区三区在线看| 又大又爽又粗| 国产亚洲精品久久久久久毛片| 国产精品av久久久久免费| 亚洲av成人一区二区三| 免费高清视频大片| 精品久久久久久久久久久久久| 丁香欧美五月| 女警被强在线播放| 亚洲av成人不卡在线观看播放网| 午夜精品在线福利| 麻豆久久精品国产亚洲av| 69av精品久久久久久| 欧美激情久久久久久爽电影| 黄色视频,在线免费观看| 欧美日本视频| 亚洲精品国产精品久久久不卡| 成年免费大片在线观看| 亚洲真实伦在线观看| 亚洲成av人片免费观看| 成人特级黄色片久久久久久久| 一夜夜www| 搞女人的毛片| 亚洲av电影在线进入| 日韩欧美在线乱码| 日韩欧美三级三区| 欧美一级毛片孕妇| 日韩精品免费视频一区二区三区| 亚洲人成伊人成综合网2020| 日韩精品青青久久久久久| 在线观看午夜福利视频| 麻豆成人av在线观看| 狠狠狠狠99中文字幕| 亚洲一区中文字幕在线| 最近视频中文字幕2019在线8| 韩国av一区二区三区四区| 欧美日本亚洲视频在线播放| 女警被强在线播放| 午夜福利18| 可以在线观看的亚洲视频| 嫩草影视91久久| 国产亚洲av嫩草精品影院| 不卡一级毛片| 又大又爽又粗| 亚洲 欧美 日韩 在线 免费| 1024视频免费在线观看| 夜夜看夜夜爽夜夜摸| 人成视频在线观看免费观看| 一边摸一边做爽爽视频免费| 免费搜索国产男女视频| 亚洲熟女毛片儿| 夜夜爽天天搞| 中文资源天堂在线| 香蕉丝袜av| 亚洲人成网站高清观看| 免费电影在线观看免费观看| 美女 人体艺术 gogo| 黄色 视频免费看| 俺也久久电影网| 白带黄色成豆腐渣| 日本免费a在线| 人人妻人人看人人澡| 亚洲激情在线av| 精品一区二区三区av网在线观看| 91麻豆精品激情在线观看国产| 91老司机精品| 国产免费av片在线观看野外av| 欧美成人免费av一区二区三区| 免费看a级黄色片| 一边摸一边抽搐一进一小说| 手机成人av网站| 国模一区二区三区四区视频 | 国产视频一区二区在线看| 18美女黄网站色大片免费观看| 欧美 亚洲 国产 日韩一| 亚洲黑人精品在线| 午夜激情av网站| 午夜免费激情av| 国产精品国产高清国产av| 午夜福利在线观看吧| 岛国在线观看网站| 男插女下体视频免费在线播放| 精品福利观看| 后天国语完整版免费观看| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 国产精品国产高清国产av| 国产欧美日韩一区二区精品| 黄色成人免费大全| 夜夜爽天天搞| 亚洲人成网站在线播放欧美日韩| 久久亚洲真实| 成人特级黄色片久久久久久久| 亚洲国产精品成人综合色| 亚洲欧美精品综合一区二区三区| cao死你这个sao货| 精品不卡国产一区二区三区| 免费人成视频x8x8入口观看| 黑人巨大精品欧美一区二区mp4| 国产成人精品久久二区二区91| 黄色成人免费大全| 国产精品爽爽va在线观看网站| 久久草成人影院| 久久99热这里只有精品18| 热99re8久久精品国产| 麻豆久久精品国产亚洲av| 中亚洲国语对白在线视频| 中文亚洲av片在线观看爽| 最近最新免费中文字幕在线| 日本黄大片高清| 黄片大片在线免费观看| 又粗又爽又猛毛片免费看| 51午夜福利影视在线观看| 亚洲 欧美一区二区三区| 麻豆国产97在线/欧美 | 搡老岳熟女国产| 两人在一起打扑克的视频| 国产午夜福利久久久久久| 色播亚洲综合网| 亚洲中文字幕日韩| 一级片免费观看大全| 一边摸一边做爽爽视频免费| 啪啪无遮挡十八禁网站| 欧美日韩黄片免| 一边摸一边抽搐一进一小说| 少妇粗大呻吟视频| www.www免费av| 少妇粗大呻吟视频| 亚洲欧美日韩高清专用| 欧美又色又爽又黄视频| 一进一出抽搐动态| www日本黄色视频网| 久久久久国内视频| 亚洲成人久久性| 99久久国产精品久久久| www.自偷自拍.com| 色综合欧美亚洲国产小说| 久久国产乱子伦精品免费另类| 91麻豆av在线| 日韩成人在线观看一区二区三区| 久久久久精品国产欧美久久久| 久久婷婷人人爽人人干人人爱| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品美女久久久久99蜜臀| netflix在线观看网站| 毛片女人毛片| 亚洲精品一区av在线观看| 99热这里只有精品一区 | 国产av不卡久久| 精品一区二区三区av网在线观看| 中文字幕熟女人妻在线| 中文资源天堂在线| 亚洲av熟女| 国产精品乱码一区二三区的特点| 免费搜索国产男女视频| 亚洲午夜精品一区,二区,三区| 91国产中文字幕| 精品国产超薄肉色丝袜足j| 非洲黑人性xxxx精品又粗又长| 国产三级中文精品| 婷婷丁香在线五月| avwww免费| 怎么达到女性高潮| 亚洲一区中文字幕在线| 51午夜福利影视在线观看| 亚洲一区中文字幕在线| 51午夜福利影视在线观看| 不卡av一区二区三区| cao死你这个sao货| 又黄又爽又免费观看的视频| 国产黄a三级三级三级人| 丁香六月欧美| 亚洲精品美女久久久久99蜜臀| 国产精品野战在线观看| 人成视频在线观看免费观看| 在线观看66精品国产| 老鸭窝网址在线观看| 亚洲一区二区三区色噜噜| av福利片在线| 久久亚洲精品不卡| 日韩欧美在线二视频| 国产av不卡久久| 法律面前人人平等表现在哪些方面| 麻豆国产97在线/欧美 | 美女大奶头视频| 狂野欧美激情性xxxx| 精品国产美女av久久久久小说| 亚洲第一电影网av| 欧美一区二区国产精品久久精品 | 日韩欧美一区二区三区在线观看| 黄色丝袜av网址大全| 国产精品爽爽va在线观看网站| 精品久久久久久久人妻蜜臀av| 久久中文字幕一级| 老司机深夜福利视频在线观看| 香蕉国产在线看| 曰老女人黄片| 高清毛片免费观看视频网站| 午夜影院日韩av| 人妻久久中文字幕网| 国产午夜精品久久久久久| 9191精品国产免费久久| 狠狠狠狠99中文字幕| 日韩高清综合在线| www日本在线高清视频| 麻豆国产av国片精品| 成人国产一区最新在线观看| 国产片内射在线| 丁香欧美五月| svipshipincom国产片| 亚洲人成77777在线视频| 熟妇人妻久久中文字幕3abv| 好男人在线观看高清免费视频| 一边摸一边做爽爽视频免费| 国产69精品久久久久777片 | 老汉色∧v一级毛片| 特大巨黑吊av在线直播| 色噜噜av男人的天堂激情| 久久 成人 亚洲| 老鸭窝网址在线观看| 日本黄色视频三级网站网址| 欧美一级毛片孕妇| 久久婷婷人人爽人人干人人爱| 黑人操中国人逼视频| 免费观看人在逋| 两个人的视频大全免费| 亚洲男人的天堂狠狠| 亚洲男人天堂网一区| 国内精品一区二区在线观看| 这个男人来自地球电影免费观看| 亚洲成人精品中文字幕电影| 成人三级做爰电影| 九九热线精品视视频播放| 亚洲成人免费电影在线观看| 操出白浆在线播放| 久久精品国产清高在天天线| 国产激情欧美一区二区| 久久天躁狠狠躁夜夜2o2o| 亚洲中文字幕一区二区三区有码在线看 | 欧美成狂野欧美在线观看| 国产三级中文精品| 两个人看的免费小视频| 欧美日韩亚洲综合一区二区三区_| 亚洲色图av天堂| av片东京热男人的天堂| 一区福利在线观看| 久久 成人 亚洲| 国产精品九九99| 亚洲人成网站在线播放欧美日韩| 欧美中文日本在线观看视频| 欧美午夜高清在线| 天天一区二区日本电影三级| 十八禁人妻一区二区| 搡老岳熟女国产| 午夜福利成人在线免费观看| or卡值多少钱| 日本黄大片高清| 夜夜夜夜夜久久久久| 变态另类丝袜制服| 亚洲国产日韩欧美精品在线观看 | 亚洲精品粉嫩美女一区| 国产午夜精品论理片| 精品熟女少妇八av免费久了| 免费av毛片视频| 亚洲av中文字字幕乱码综合| 在线观看免费日韩欧美大片| 波多野结衣高清无吗|