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

    LEO-LEO微波掩星探測溫度和水汽廓線研究進展*

    2022-02-11 07:56:10張志華呂達仁
    遙測遙控 2022年1期
    關鍵詞:掩星廓線水汽

    張志華,王 鑫,呂達仁

    LEO-LEO微波掩星探測溫度和水汽廓線研究進展*

    張志華1,2,王 鑫*1,呂達仁1,2

    (1中國科學院大氣物理研究所 北京 100029 2中國科學院大學 北京 100049)

    低地球軌道衛(wèi)星間(LEO-LEO)微波掩星探測技術利用水汽吸收線附近信號,通過測量信號的折射和吸收作用,能夠在沒有輔助背景信息的情況下獨立反演溫度和水汽廓線,從而克服GNSS-LEO無線電掩星技術在大氣中低對流層存在的“溫度–水汽模糊”問題。目前,國內(nèi)外還沒有實現(xiàn)在軌的LEO-LEO掩星探測。為了促進我國LEO-LEO掩星探測技術發(fā)展,本文對LEO-LEO微波掩星探測溫度和水汽廓線研究前沿及進展進行了回顧和總結。首先介紹了LEO-LEO掩星探測技術的發(fā)展歷史,闡述了LEO-LEO掩星探測的基本原理和探測頻點選擇,重點介紹了LEO-LEO微波掩星探測溫度和水汽廓線的反演技術研究動態(tài),并介紹了LEO-LEO掩星探測任務國際進展,最后進行了總結和展望。本研究可為實現(xiàn)LEO-LEO掩星空間探測技術提供基本科學理論參考。

    LEO-LEO;微波掩星探測;反演;溫度廓線;水汽廓線

    引 言

    近些年來,臺風、海嘯、暴雨等極端天氣的頻發(fā)帶來巨大災害,已使全球氣候變化與應對策略成為關乎國家發(fā)展與人民生活的極端重要科學問題。溫度和水汽是描述地球大氣狀態(tài)的兩個重要參數(shù),具有明顯的四維時空變化特性,直接影響太陽短波輻射與地–氣系統(tǒng)長波輻射的相互作用,進而影響全球的輻射能量收支平衡。政府間氣候變化專門委員會(IPCC)第六次評估報告顯示,相對于1850–1900年,2010–2019年全球地表平均溫度已上升約1 ℃,并指出從未來20年的平均溫度變化來看,全球溫升預計將達到或超過1.5 ℃[1]。報告預估,在未來幾十年里,所有地區(qū)的氣候變化都將加劇。水是地球大氣中唯一能以固、液、氣三相態(tài)存在的物質(zhì),且存在三相態(tài)互相變化,在云的演變、降水的發(fā)生及變化中都起著重要的作用。水汽在地球大氣中所占比例很?。▋H0.1%~3%),卻是大氣中最活躍的部分,對地球生態(tài)系統(tǒng)有重大影響[2]。水汽還是地球大氣的主要溫室氣體,它對于太陽短波是透明的,但能夠吸收地球的長波輻射,從而使大氣溫度升高[3]。認識和了解溫度和水汽的時空分布與變化規(guī)律,特別是高垂直分辨率、高精度、長期穩(wěn)定的垂直分布特征,是研究認識大氣能量物質(zhì)交換、全球變化等重大問題的科學基礎,對氣候變化研究、天氣現(xiàn)象分析及數(shù)值天氣預報等都具有重要意義。

    1960年4月,美國國家航空航天局(NASA)成功發(fā)射第一顆氣象衛(wèi)星Tiros-1,開創(chuàng)了空間探測科學研究的新紀元,衛(wèi)星遙感技術具有視野開闊、觀測范圍廣、觀測時次多等優(yōu)點,成為探測地球大氣的有力手段[4]。利用衛(wèi)星探測地球大氣垂直結構,經(jīng)歷了多個發(fā)展階段,從初期僅有少數(shù)幾個紅外寬波段通道、低空間分辨率觀測,逐步發(fā)展到紅外儀器、微波儀器同平臺跨衛(wèi)星軌道掃描和較高空間分辨率觀測。然而,紅外探測設備不能穿透云層,地球上大約有三分之二的區(qū)域覆蓋了云,這給紅外遙感系統(tǒng)造成了基本的采樣困擾,全球對流層的紅外采樣本質(zhì)上是不完整的和偏干的,例如AIRS像素的90%~95%受到云污染[5]。另外,強對流天氣系統(tǒng)通常形成在多云區(qū)域,而紅外儀器無法穿透,無法滿足盡可能早地捕獲強對流天氣系統(tǒng)的愿望。被動微波探測器能夠穿透云層,例如衛(wèi)星微波輻射計,通過被動接收大氣輻射的微波信號,探測輻射亮溫信息,從而獲得大氣溫度和濕度廓線。但由于大陸地表溫度的多變性和復雜性,此技術用于海洋表面的遙感比陸地更為適合,且衛(wèi)星微波輻射計的垂直分辨率有限,例如先進微波探測儀(AMSU-B),其垂直分辨率約為2 km。此外,被動微波臨邊探測儀(MLS)可以提供對流層中上層的大氣廓線,垂直分辨率約為2 km~3 km,但容易受冰云散射的限制,反演得出的對流層采樣不完整且有偏[6,7]。

    二十世紀九十年代初,美國的全球定位系統(tǒng)GPS(Global Positioning System)星座的建成為無線電掩星探測地球大氣創(chuàng)造了條件[8]。全球導航定位系統(tǒng)GNSS(Global Navigation and Positioning System)無線電掩星探測可以實現(xiàn)大氣折射率、溫度和水汽等大氣參數(shù)的高垂直分辨率、高精度探測,GNSS導航信號穿越大氣層時,不受云、雨和氣溶膠的影響,彌補了衛(wèi)星可見光、紅外遙感在這方面的缺陷,具有真正全天候的優(yōu)勢[9]。此外,GNSS-LEO無線電掩星GRO(GNSS-LEO Radio Occultation)探測是自定標系統(tǒng),具有自校準能力,衛(wèi)星成功運行后,不再需要進行校正和調(diào)整,具有長期穩(wěn)定的特點[10]。GNSS-LEO無線電掩星技術作為目前最先進的空間探測技術之一,探測資料已被各大數(shù)值天氣預報中心同化到數(shù)值天氣預報模式中,對天氣預報和氣候研究產(chǎn)生了積極的影響。然而,在對流層低層較暖區(qū)域,水汽比較豐富,僅使用GNSS導航信號進行掩星探測,無法準確區(qū)分溫度和水汽對折射率的貢獻,存在“溫度–水汽模糊”問題[11,12]。針對這一問題,二十世紀九十年代末,低地球軌道衛(wèi)星間(LEO-LEO)微波掩星探測的概念被提出。LEO-LEO微波掩星LMO(LEO-LEO Microwave Occultation)采用水汽吸收線附近的微波信號[12],基于電磁波在大氣中傳播時的折射和吸收效應,獲取電磁波信號的相位延遲和振幅衰減信息,能夠在沒有輔助大氣背景場的情況下獨立反演溫度和水汽廓線。

    LEO-LEO微波掩星探測具有高精度、高垂直分辨率、高穩(wěn)定性、全天候、全球分布等特點[13,14],它獨立于現(xiàn)有的探測方法,是現(xiàn)有地球大氣探測手段的重要補充。發(fā)展LEO-LEO掩星探測技術,對提高我國在大氣領域的探測能力和相關研究水平具有重要的科學意義和應用價值。到目前為止,國內(nèi)外還沒有實現(xiàn)在軌的LEO-LEO掩星探測衛(wèi)星,回顧LEO-LEO微波掩星探測技術的發(fā)展過程和研究其未來發(fā)展趨勢,可以更好地促進我國LEO-LEO微波掩星探測技術的發(fā)展。

    本文重點對LEO-LEO微波掩星探測溫度和水汽廓線的研究進展進行了介紹和總結。首先介紹了LEO-LEO掩星探測技術的發(fā)展歷史,闡述了LEO-LEO掩星探測的基本原理和微波掩星探測溫度、水汽廓線的技術途徑與反演方法,其次介紹了國際LEO-LEO掩星探測計劃的前沿進展,最后進行了總結和展望。

    1 LEO-LEO掩星探測技術發(fā)展歷史

    掩星技術最早可以追溯到十八世紀[15],天文學家利用月球和行星掩星技術進行時間測定,在行星掩星觀測過程中發(fā)現(xiàn),當恒星光線進入行星大氣層時,強度會逐漸變?nèi)?,這個變?nèi)醯倪^程與行星大氣標高有關,由此提出利用自然光源如太陽或星光進行掩星行星大氣探測。二十世紀六十年代,星際飛行器問世,利用無線電掩星技術探測行星大氣廓線的設想得以實現(xiàn)。在地震波反演地層結構技術的基礎上,無線電掩星技術的數(shù)學方法被提出,用于行星大氣參數(shù)廓線的探測反演研究,Abel變換法被引入到從彎曲角到折射率的反演過程中[16]。隨后,NASA進行了一系列太陽系行星探測計劃,通過水手號、海盜號和先驅者號等探測器,成功進行了行星大氣的掩星觀測,獲得了幾乎所有太陽系行星及其衛(wèi)星的大氣廓線[17]。然而,利用無線電掩星技術探測地球大氣的設想一直沒有實現(xiàn),主要原因是掩星探測要求信號的發(fā)射端和接收端都必須在被探測體之外,而在地球周圍一直沒有找到合適的觀測點,缺乏廉價信號源使得它的性價比遠低于當時其他的大氣探測手段,導致沒有相關的探測項目得以實施[8]。二十世紀七十年代,美國和蘇聯(lián)出于軍事考慮,開始籌建各自的全球定位系統(tǒng)——GPS系統(tǒng)和GLONASS系統(tǒng)。二十世紀九十年代初,美國GPS導航星座建成,為地球大氣掩星觀測提供了無線電發(fā)射端,以相對較低的成本為掩星技術探測地球大氣創(chuàng)造了條件。1995年4月3日,搭載GPS接收機的Microlab-1低軌實驗衛(wèi)星被成功發(fā)射,首次從工程方面驗證了GPS掩星探測地球大氣的可行性,為無線電掩星技術探測地球大氣拉開了序幕[8]。

    GPS/MET掩星觀測試驗儀器在軌運行了約兩年,成功進行了相當數(shù)量的掩星觀測,獲得了大量氣壓、溫度和水汽廓線等寶貴大氣資料。Kursinski等研究認為[9,18],在低對流層較暖區(qū)域(>240 K),水汽比較豐富,溫度和水汽對總折射率的貢獻無法準確區(qū)分,造成“溫度–水汽模糊性”問題。如果忽略水汽,利用實部折射率方程、流體靜力方程和理想氣體狀態(tài)方程,可求解得到干空氣假設下氣壓、密度和溫度廓線,干空氣的溫度通常在7 km以下比實際溫度偏低;如果不忽略水汽,三個方程約束密度、壓力、溫度和水汽四個未知變量,在不引入輔助先驗溫度信息時,無法反演得到水汽廓線[19]。Poli等(2002)提出將一維變分(1D Var)同化反演方法應用于濕大氣廓線的反演[20],以數(shù)值模式數(shù)據(jù)作為輔助背景場,可同時反演得到溫度和水汽廓線,但其反演精度非常依賴模式初始場和模式誤差,不同的模式數(shù)據(jù)有其固有的系統(tǒng)偏差,采用不同的模式數(shù)據(jù)作為背景場會導致不同的結果[21,22]。

    為了克服“溫度–水汽模糊性”問題,獨立反演獲得溫度和水汽廓線,在GNSS-LEO無線電掩星探測技術的基礎上,1995年,美國亞利桑那大學的研究人員開始考慮利用高于GPS的信號頻率來測量因22 GHz水汽共振線引起的水汽吸收,探索其進行掩星探測的可能性[14,23]。1997年,美國噴氣推進實驗室和亞利桑那大學的研究人員對22 GHz附近頻率進行掩星探測的性能進行了初步評估。1998年,前期相關研究成果被總結成兩大提案提交給NASA,即AMORE任務和ATOMS任務。其中,ATOMS提案獲得資助,研究人員對22 GHz和183 GHz水汽吸收線附近的微波掩星探測能力進行了評估研究[23]。Feng等(2001)和Kursinski等(2002)在GNSS-LEO無線電掩星實部折射率反演算法的基礎上[11,12],發(fā)展了LEO-LEO微波掩星虛部折射率反演算法,利用相鄰頻率的信號強度比消除散焦、閃爍和衍射等非吸收效應,得到僅由吸收引起的大氣差分透過率和大氣吸收系數(shù)。2005年,在前期研究基礎上,ATOMMS計劃被提出,用以實現(xiàn)近地面到中間層大氣溫度、水汽和氣壓的高精度、高垂直分辨率探測[9]。2007年,Kursinski等提出研制ATOMMS原型機,并進行飛機對飛機掩星演示試驗,但由于資金限制,飛機演示試驗未得到實施,僅成功研制了ATOMMS原理樣機。2010年3月,Kursinski等在亞利桑那大學校園內(nèi)進行了ATOMMS原型機“高頻段”182 GHz~205 GHz的屋頂試驗,得到了水汽吸收光譜的時間序列,由于無法測量沒有水汽吸收的“大氣頂”的信號強度,沒有得到絕對水汽含量,僅得到沿路徑上水汽含量隨時間的變化,試驗證明了ATOMMS儀器對吸收光譜線形狀的敏感性[24]。2010年8月,Kursinski等在亞利桑那州圖森市北部相距約5.4 km的兩山之間進行了ATOMMS原型機“低頻段”22 GHz通道的山頂試驗。2011年8月,又在兩山之間進行了第二次山頂試驗,在晴朗、多云和雷暴天氣條件下,測量了ATOMMS儀器在信號頻率24.4 GHz和198.5 GHz的水汽吸收光譜,得到了高達17的光學厚度條件下的精確水汽結果,在該光學厚度以上,信號已不可檢測[25]。

    在歐洲,早在1996年就提出通過添加Ku波段信號來分離溫度和水汽。2001年,WATS計劃被提出,用以探測對流層和平流層的水汽和溫度垂直分布,并對22 GHz附近頻率掩星探測的性能進行了初步評估[26]。2002年,ACE+計劃被提出,它由WATS計劃的LEO-LEO微波掩星任務與ACE計劃的GRO任務合并而成,兩顆小衛(wèi)星攜帶X和K波段的發(fā)射機,另外兩顆小衛(wèi)星攜帶X和K波段的掩星接收機,測量LEO-LEO鏈路之間的電波相位和振幅變化[6]。2004年,ACCURATE計劃被奧地利格拉茨大學魏格納氣候研究中心提出,目標是對上對流層-下平流層(UTLS)區(qū)域大氣氣候化學和氣候趨勢進行探測,它不僅繼承了ACE+計劃的LEO-LEO微波掩星探測技術,還進行低軌道衛(wèi)星間紅外激光掩星探測(LIO),構成LEO-LEO微波掩星和紅外激光掩星(LMIO)概念[27,28]。2010年,ACCURATE計劃得到進一步資金支持,開展了紅外激光掩星原型機的研制,并于2013年進行了地面驗證實驗[29]。

    國內(nèi)對LEO-LEO掩星探測技術的研究起步較晚。2007年,杜曉勇等分析了發(fā)射衛(wèi)星和接收衛(wèi)星的軌道高度、傾角、近地點角距、升交點赤經(jīng)等軌道參數(shù)對LEO-LEO掩星事件數(shù)量和分布的影響[30]。呂華平等仿真研究了一天及一個月LEO-LEO掩星事件的數(shù)量及分布,并對LEO-LEO掩星探測頻點選擇進行了初步研究[31,32]。孫立剛等對LEO-LEO掩星事件持續(xù)時間及衛(wèi)星間相對角速度進行了仿真分析[33]。2017年,中國科學院大氣物理所呂達仁等提出了我國首個LEO-LEO掩星探測任務——CACES(Climate and Atmospheric Components Exploring Satellites)計劃,旨在通過低地球軌道衛(wèi)星間微波與紅外激光掩星探測手段,獲取特定頻率與譜段的信號相位和振幅變化,反演得到氣壓、溫度和水汽等大氣參數(shù)廓線,以及H2O、CO2、CH4和O3等多種氣體成分廓線。中國科學院國家空間科學中心研制了LEO-LEO微波掩星原理樣機,探測頻點在22.23 GHz水汽吸收線附近[34]。2018年,馮梓宸等進行了LEO-LEO掩星信號接收技術研究[35],設計研制了一套針對60 km~90 km大氣探測的LEO-LEO掩星接收機原理樣機。

    自LEO-LEO掩星探測概念提出以來,依托上述提出的掩星探測計劃,研究人員開展了相關載荷研制和原理樣機地面試驗驗證等工作,并在LEO-LEO掩星軌道設計、頻率選擇評估、反演算法改進、探測性能分析等方面開展了一系列仿真模擬研究,為LEO-LEO掩星探測技術的實現(xiàn)打下了堅實的基礎,積累了寶貴的經(jīng)驗。但由于技術條件、風險控制和預算調(diào)整等多種原因,上述提出的計劃后續(xù)均未得到在軌實施,目前國際上仍沒有實現(xiàn)的LEO-LEO掩星空間探測任務。發(fā)展LEO-LEO掩星探測技術,使我國相關研究能力達到國際前沿水平,還需從硬件研制、反演技術等方面進行大量且深入的研究。

    2 LEO-LEO微波掩星探測原理

    2.1 LEO-LEO掩星探測方式

    LEO-LEO掩星探測是基于衛(wèi)星的雙雷達系統(tǒng),以主動臨邊掃描方式探測地球大氣層。LEO-LEO掩星事件是指兩顆低軌衛(wèi)星LEO(Low Earth Orbit)分別運動到地球的兩側,其中一顆作為發(fā)射衛(wèi)星,另一顆作為接收衛(wèi)星。由發(fā)射衛(wèi)星發(fā)出電波信號,穿過地球大氣到達接收衛(wèi)星,由于兩顆LEO衛(wèi)星之間的相對運動,發(fā)射信號從大氣層頂逐漸切過整層地球大氣直到地球表面。或者相反,從切過地球表面直到切過大氣層頂,稱為發(fā)生了一次LEO-LEO掩星事件[8]。LEO-LEO掩星探測技術的發(fā)射和接收衛(wèi)星軌道高度一般都在300 km~900 km高度范圍,收發(fā)距離約6 000 km,相比GNSS-LEO無線電掩星探測距離相對較短,并且發(fā)射信號源的特性可以控制,能夠提供較高的信噪比[36]。

    GNSS-LEO無線電掩星探測技術使用導航衛(wèi)星信號探測地球大氣,GNSS信號頻點最初的選擇是為了最大限度地減少大氣信息的影響,用于氣象觀測會受到特定波長(L波段,1 GHz~2 GHz)的限制。與GNSS-LEO無線電掩星探測不同,LEO-LEO微波掩星探測技術可采用22.23 GHz和183.31 GHz水汽吸收中心線附近的頻率作為發(fā)射頻率探測地球大氣,如圖1所示,微波信號從發(fā)射衛(wèi)星穿過地球大氣層時被折射和吸收,其傳播速度和路徑發(fā)生變化,還會引起電波信號衰減,這些效應的綜合導致電波信號的相位和振幅發(fā)生變化。通過測量電波信號的相位延遲和振幅衰減信息,可以在沒有輔助大氣背景信息的條件下同時反演得到密度、氣壓、溫度、水汽等大氣參數(shù)廓線。

    圖1 LEO-LEO微波掩星探測示意圖

    2.2 LEO-LEO微波掩星探測頻點選擇

    GNSS-LEO掩星探測技術是利用LEO上搭載的接收機接收GNSS導航信號,發(fā)射頻點固定,不需要考慮發(fā)射機的信號頻率。但LEO-LEO微波掩星不同,它除了需要在一顆LEO上搭載接收機,還需要在另外一顆LEO上搭載發(fā)射機,發(fā)射機頻率通道的選擇是LEO-LEO掩星探測技術需要考慮的重要問題之一。

    LEO-LEO微波掩星的信號穿過地球大氣時,會造成信號相位延遲和大氣損耗,總的大氣損耗由散焦和吸收兩部分引起,一般分開考慮。在各種大氣成分中,水汽和氧氣對電波信號的吸收衰減起主要作用。在200 GHz以下,存在22.23 GHz和183.31 GHz兩條水汽吸收中心線。通常,在頻率0處靠近共振吸收頻率的氣體(水汽)吸收系數(shù)()可以表示為[37]

    大氣中的水汽含量在整個對流層和平流層的變化超過幾個數(shù)量級,不同吸收線在不同高度的吸收強度明顯不同,因此必須將大氣分成幾個高度區(qū)域,為每個高度區(qū)域選擇適當?shù)奶綔y頻率。Yunck等(2000)評估分析了22.23 GHz水汽吸收線探測大氣濕度的能力[10],結果表明,23 GHz頻率通道對探測相對干燥的上對流層大氣很重要,而10 GHz、14 GHz和18 GHz頻率通道對探測潮濕的下對流層大氣很重要。Eriksson等(2003)評估分析了利用10 GHz~35 GHz頻率通道探測對流層和平流層水汽的性能[26],研究表明,10.3 GHz、17.2 GHz、22.6 GHz足夠用于水汽信息探測,27.4 GHz和32.9 GHz對熱噪聲更敏感,水汽反演精度更差。Herman等(2003)利用毫米波傳播模型(MPM93)計算并分析了200 GHz以下頻率通道在不同濕度條件下的吸收光譜,包括水汽和氧氣在近地面、約15 km、30 km和50 km高處的吸收衰減,研究結果表明,183.31 GHz水汽吸收線的吸收系數(shù)更大,被用于探測水汽含量低的高空,而22.23 GHz水汽吸收線被用于探測水汽含量豐富的中低對流層。由于散焦、衍射和閃爍等誤差影響,一般選擇兩個或更多頻率同時測量信號衰減,靠近水汽吸收線中心的頻率作為主頻率,測量強信號吸收,遠離吸收線中心的頻率作為參考頻率,測量弱信號吸收用于校正。兩個頻率的大氣透過率相減得到差分透過率,計算得到差分吸收系數(shù),能夠消除背景吸收和散焦等不必要影響[11-13]。LEO-LEO微波掩星探測利用22.23 GHz和183.31 GHz水汽吸收線附近的頻率,通過在不同高度范圍內(nèi)使用不同信號頻率對之間的差分吸收系數(shù)進行計算,能夠探測整個對流層和平流層的溫度和水汽廓線。

    從科學目標的角度來看,LEO-LEO微波掩星選擇頻率通道時,希望盡可能選擇接近水汽吸收中心線的頻點,以便更靈敏地探測水汽,并達到更大的高度范圍。但是,頻率通道的選擇涉及科學目標與國際無線電頻率分配之間的權衡[36]。例如,國際無線電頻率規(guī)定禁止直接使用吸收線中心的頻點進行主動探測,水汽吸收線中心頻率及其相鄰頻率已經(jīng)被保護用于其他專門的科學和商業(yè)應用,不能用于LEO-LEO微波掩星空間任務的實際探測[13,14]。因此,在LEO-LEO微波掩星相關載荷研制設計中,選擇載荷頻點時,必須在滿足國際無線電頻率規(guī)則的情況下,根據(jù)科學探測目標,在允許的頻率段范圍內(nèi),選擇能夠提供最佳水汽信息的頻點。

    3 LEO-LEO微波掩星探測反演研究

    3.1 溫度和水汽廓線反演

    LEO-LEO微波掩星探測的原始觀測量是電波信號的相位和振幅,在已知各種大氣參數(shù)條件下,正向模擬出掩星觀測的信號相位和振幅,是正演過程。由LEO-LEO掩星探測的電波信號相位和振幅數(shù)據(jù),獲得傳播路徑上的大氣介質(zhì)信息,是信號傳播的一個反演過程。

    在幾何光學近似條件下,利用LEO-LEO微波掩星反演溫度和水汽等大氣參數(shù)廓線,可以總結為三大部分[8,36,38],即:

    ①實部折射率反演。首先,對模擬或觀測到的相位延遲數(shù)據(jù)進行平滑處理,去除異常值,并進行電離層校正,剔除電離層引起的電波彎曲和相位延遲;然后,計算多普勒頻移,結合發(fā)射衛(wèi)星和接收衛(wèi)星的精密軌道信息,導出彎曲角廓線,并對高空彎曲角進行統(tǒng)計優(yōu)化處理;最后,利用經(jīng)典Abel積分變換公式,由彎曲角廓線反演得到實部大氣折射率廓線。

    ②虛部折射率反演。LEO-LEO微波掩星的電波信號因吸收作用導致振幅衰減,反映了射線路徑上的水汽含量。首先,對LEO-LEO微波掩星測得的信號振幅進行散焦和球面信號擴展校正以及歸一化處理,根據(jù)Bouguer-Lambert-Beer定律,由振幅計算得到大氣透過率廓線;然后,結合反演得到的彎曲角廓線和實部折射率廓線,通過另外的Abel積分變換可以反演得到大氣吸收系數(shù)廓線,進而求得虛部折射率廓線。

    ③大氣參數(shù)廓線反演。利用以上由相位延遲和振幅求得的實部折射率和虛部折射率,結合大氣靜力學方程、理想氣體狀態(tài)方程、光譜方程,建立非線性方程組。通過對靜力學方程向下積分,依次求解出每個高度層的氣壓、水汽、溫度等大氣參數(shù)。

    LEO-LEO微波掩星探測的實部折射率反演部分與GNSS-LEO無線電掩星探測的基本相同,虛部折射率反演是LEO-LEO微波掩星探測特有的反演過程,每個頻率對應一個虛部折射率廓線。在濕大氣條件下,在水汽的貢獻不能忽略時,有密度、氣壓、溫度和水汽四個未知變量,LEO-LEO微波掩星探測將多個頻率通道對應的虛部折射率方程與實部折射率方程、流體靜力方程和濕空氣狀態(tài)方程聯(lián)立,至少有四個方程可用于約束四個未知變量,非線性方程組是超定方程組,不需要引入任何輔助背景場,即可獨立求解得到溫度和水汽廓線,克服了GNSS-LEO無線電掩星探測存在的“溫度–濕度模糊”問題。

    Kursinski等(2002)研究了LEO-LEO微波掩星反演大氣透過率和大氣吸收系數(shù)的基本方法,通過在不同高度范圍內(nèi)使用不同信號頻率對之間的差分透過率進行計算,水汽反演誤差在約1 km~75 km內(nèi)可控制在約1%~3%之內(nèi),溫度反演誤差從幾千米到80千米可在亞開爾文量級[12]。Syndergaard等(2004)研究表明,在吸收系數(shù)水平梯度極端的情況下,由于不滿足球對稱假設,導致在5 km以下對流層低層水汽反演的不確定性為百分之幾十[39]。Schweitzer等(2011)研究發(fā)現(xiàn)氣壓、溫度和水汽廓線的誤差分別小于0.2%、0.5 K和10%,準確度基本滿足探測目標要求,在三個頻率均位于K波段時,水汽在約11 km以下的相對誤差小于10%,另外增加兩個183 GHz水汽吸收中心線附近的頻率時,準確反演能達到約18 km[36]。

    在對流層低層,水汽含量往往較為豐富并伴隨小尺度的劇烈變化,導致大氣折射率急劇變化,從而使電波信號在傳播時產(chǎn)生多路徑效應。在大氣多路徑條件下,通過幾何光學方法反演得到的彎曲角是碰撞參數(shù)的多值函數(shù)。針對這個問題,Gorbunov和Kirchengast等進行了波動光學方法反演水汽廓線的研究[40-42],利用正則變換方法(CT)由相位和振幅計算得到彎曲角和大氣透過率廓線,然后反演獲得水汽廓線,結果表明,在約3 km處水汽反演準確度可達到0.5 g/kg。

    3.2 云的影響

    LEO-LEO微波掩星信號可以穿透云層探測地球大氣信息,但電波信號在穿過云層時,會因不同云層結構的存在被不同程度地吸收衰減。云的微物理特性依賴于水粒子的尺寸、形狀和相態(tài)。通常情況下,水滴粒子是球形的,且小于0.1 mm。冰晶可以組成各種各樣的尺寸和形狀,在自然界中,除了簡單的多面體外,也有很多不規(guī)則冰晶和簡單形狀的混合體。當有云層存在時,22.23 GHz水汽吸收線的波長約為13.5 mm,明顯大于云中液態(tài)水的顆粒尺寸,散射相比吸收作用可以忽略不計[8]。一般來說,云中液態(tài)水的吸收大約比冰晶的吸收大兩個數(shù)量級。當有云中冰晶時,冰晶相對于液態(tài)水粒子直徑更大,對微波信號的散射作用更強[26],且散射作用的影響隨著頻率的增加而增加,故冰晶對頻率為183.31 GHz的散射影響比22.23 GHz的更大。

    在對流層下部,除了水汽對LEO-LEO微波掩星信號的吸收作用外,如果有云層存在,云中液態(tài)水也會對總吸收做出貢獻,為了反演得到高精度的溫度和水汽廓線,需要去除云中液態(tài)水對總吸收的貢獻。對于LEO-LEO微波掩星探測技術,如果提供足夠的頻率通道,可以區(qū)分開水汽和云中液態(tài)水引起的吸收。在反演溫度和水汽廓線時,云中液態(tài)水廓線可以作為副產(chǎn)品被同時反演出來。Kursinski等(2004)初步討論了云層對水汽和溫度反演精度的影響[43],結果表明,在對流層最低層存在液態(tài)水云的情況下,水汽反演精度比晴空條件下約低2倍。ATOMMS任務擬采用22 GHz水汽吸收線附近的至少五個頻率,以分離液態(tài)水與水汽的影響,擬采用183 GHz水汽吸收線附近頻率,探測對流層上層的水汽信息。

    在理論上,云中冰晶也可以利用LEO-LEO微波掩星技術被探測并反演出來。183 GHz附近的頻率對冰粒子的散射比較敏感[43],利用吸收通道和參考通道的比率可以消除幾乎大部分冰晶對溫度和水汽探測的影響。Kursinski等(2004)提出,通過測量單個頻率的衰減曲線來獲得冰晶信息,通過測量至少2個頻率間隔較寬的信號來獲得有關冰量和其他屬性的信息,具體探測和反演技術未來仍需要開展深入研究。

    3.3 湍流的影響

    電波信號在湍流介質(zhì)中傳播時,通過衍射產(chǎn)生干擾,從而產(chǎn)生相位和振幅閃爍,在測量水汽吸收特征時,湍流引起的振幅閃爍是一種噪聲源[9]。Kursinski等研究發(fā)現(xiàn),通過使用吸收通道和參考通道之間的差分大氣透過率,可以在很大程度上消除絕對大氣透過率中常見的不必要的寬帶效應,特別是大氣湍流中的振幅閃爍影響[11,12]。

    Kursinski和Otarola等估計分析了對流層中層22 GHz水汽吸收線附近的濕湍流貢獻和對流層上層183 GHz水汽吸收線附近的干湍流貢獻[9],研究結果表明,對于22 GHz和183 GHz,單頻率原始振幅變化的標準偏差大約為7%~8%,兩頻率振幅比的標準偏差約為1%~2%。較短的波長更容易受到湍流的影響,盡管原始振幅閃爍在22 GHz和183 GHz吸收線相當,但對于183 GHz,參考通道離吸收線中心越近,在比值后產(chǎn)生的湍流效應的抵消效果越好。Gorbunov和Kirchengast(2007)進行了2-D湍流仿真模擬研究,評估分析了22 GHz水汽吸收線附近湍流的影響,研究結果表明,在低緯度的5 km~8 km高度范圍,9.7 GHz和17.25 GHz頻率的差分大氣透過率誤差約為0.05 dB或1%。值得注意的是,Kursinski等估計的低緯度中對流層條件下的大氣透過率誤差大約是Gorbunov等估計結果的兩倍。盡管二者用于模擬湍流效應的方法非常不同,但近似的湍流噪聲水平在某種程度上是真實的,如何更好地理解和協(xié)調(diào)兩種誤差估計方法之間的差異還需要更進一步的研究。

    在反演溫度和水汽廓線時,湍流強度可作為副產(chǎn)品被同時反演出來,即對于湍流研究來說,利用LEO-LEO掩星探測技術是一個測量獲取湍流信息的潛在途徑[6]。

    4 LEO-LEO微波掩星探測計劃進展

    目前,國際上主要有三個LEO-LEO掩星探測計劃,即歐洲的ACCURATE計劃、美國的ATOMMS計劃和中國的CACES計劃。

    ① ACCURATE計劃

    ACCURATE(Atmospheric Climate and Chemistry in the UTLS Region And climate Trends Explorer)計劃由奧地利格拉茨大學魏格納氣候研究中心在2004年提出,旨在對上對流層–下平流層(UTLS)區(qū)域大氣氣候化學和氣候趨勢進行探測。ACCURATE計劃的探測目標是5 km~35 km范圍內(nèi)的氣壓、溫度和水汽等大氣參數(shù)廓線,同時探測6種主要的溫室氣體,水汽、二氧化碳和它們的同位素,以及探測氣溶膠、云層和大氣閃爍等。該計劃繼承了ACE+計劃的LEO-LEO微波掩星探測任務,同時進行LEO-LEO紅外激光掩星探測,是首個被提出的LEO-LEO微波掩星和紅外激光掩星(LMIO)概念。LEO-LEO微波掩星探測的頻點均在K波段,分別為17.25 GHz、20.2 GHz和22.6 GHz,3個頻點可以采用1個天線,從而減少了儀器體積、重量和成本。ACCURATE星座設計為4顆LEO衛(wèi)星,分布在2個軌道平面上,軌道高度分別為650 km和850 km,每個軌道面分布2顆LEO衛(wèi)星[27]。

    依托ACCURATE計劃,Kirchengast等對EGOPS軟件包進行了升級改進,可以進行端到端的LEO-LEO微波掩星仿真和反演計算,取得了一些有意義的成果,但到目前為止,該計劃沒有開展過相關LEO-LEO微波掩星載荷的硬件研制工作。Kirchengast等研制了紅外激光掩星原型機,并在2011年進行了地面驗證實驗,研究人員對地面實驗結果進行了分析研究,之后關于該計劃的研究基本處于停滯狀態(tài)。

    ② ATOMMS計劃

    ATOMMS(Active Temperature, Ozone and Moisture Microwave Spectrometer)計劃由美國亞利桑那大學在2005年提出,即主動式溫度、臭氧和水分微波波譜儀,旨在探測從近地面到中間層的密度、氣壓、溫度和水汽廓線,以及從對流層頂?shù)酵瑴貙拥某粞趵€。該計劃組合了GPS-LEO掩星探測系統(tǒng)和LEO-LEO微波掩星探測系統(tǒng),只進行微波掩星探測。LEO-LEO微波掩星的探測頻點分別在X波段、K波段和M波段,在22.23 GHz水汽吸收線附近有5個頻點,分別是9.7 GHz、13.5 GHz、17.25 GHz、20.2 GHz和22.6 GHz;在183.31 GHz水汽吸收線附近有3個頻點,分別是179.0 GHz、181.95 GHz和183.6 GHz。ATOMMS星座設計為12顆LEO衛(wèi)星,分布在2個軌道平面,軌道高度分別為650 km和800 km,每個軌道面分布6顆LEO衛(wèi)星,一天可以發(fā)生1 000次以上LEO-LEO掩星事件和6 000~8 000次GPS-LEO掩星事件[6]。

    亞利桑那大學的研究人員研制了低頻段22.23 GHz和高頻段183.31 GHz的原理樣機,分別在樓頂和山頂進行了水汽吸收線衰減測量試驗,但后續(xù)沒有對儀器再進行改進和研制,之后關于該計劃的研究也基本處于停滯狀態(tài)。

    ③ CACES計劃

    CACES(Climate and Atmospheric Components Exploring Satellites)計劃由中國科學院大氣物理研究所在2017年提出,旨在通過低地球軌道衛(wèi)星間微波與紅外激光掩星探測手段,獲取特定頻率與譜段的信號相位和振幅變化,反演大氣溫度、濕度和壓強等大氣參數(shù),獲得H2O、CO2、CH4和O3等多種氣體成分廓線,并結合紅外成像實現(xiàn)三維立體探測。LEO-LEO微波掩星的探測頻點分別在X波段和K波段,分別為9.7 GHz、13.5 GHz、17.25 GHz和22.6 GHz。CACES星座設計為4顆LEO衛(wèi)星,分布在2個軌道平面上,軌道高度分別為500 km和550 km,每個軌道面分布2顆LEO衛(wèi)星。

    CACES計劃是目前國際上在研的LEO-LEO掩星探測計劃,前期經(jīng)過科學概念研究、深化論證和預先研究三個階段,目前順利得到中國科學院“空間科學(二期)”戰(zhàn)略性先導科技專項的后續(xù)支持,進入到背景型號研制階段,關鍵載荷紅外激光掩星與微波掩星原理樣機正在研制當中,目前計劃開展外場實驗,進行山基探測驗證。

    5 LEO-LEO微波掩星探測技術展望

    基于國際研究前沿與當前進展,認為未來LEO-LEO微波掩星探測技術研究需要在目前研究基礎上,開展以下工作:

    ①反演技術

    Kursinski等地面實驗表明[24],AM6.2光譜模型比MPM93模型更適合測量光譜線的形狀。目前的反演算法中,微波信號的模擬實部折射率和虛部折射率是由MPM93模型計算的,未來仍需要通過大量仿真實驗,對比分析MPM93、AM6.2等微波傳輸模型對溫度和水汽廓線反演精度的影響,選擇最優(yōu)的微波傳輸模型建立LEO-LEO微波掩星反演算法。此外,微波信號的實部折射率和虛部折射率的權重在低對流層的設置不是最佳的,不同信號頻率的權重分配是固定的,為了進一步提高溫度和水汽廓線的反演精度,仍需要對權重函數(shù)進行定量計算分析與改進,也可根據(jù)信號頻率的強度建立動態(tài)可變的權重因子,從而提高算法的反演精度。

    對流層低層水汽含量豐富,當存在較厚云層和湍流時,溫度和水汽廓線反演精度會受到顯著不利影響。針對此問題,Kirchengast等提出兩個解決方法[36],一是引入先驗溫度場信息,利用最佳擬合溫度外推方法反演3 km~6 km以下溫度和水汽廓線,但此方法存在的缺點是需要引入輔助背景場信息,這與LEO-LEO獨立反演溫度和水汽廓線的初衷不符;二是波動光學方法,在反演大氣透過率的過程中可以減少閃爍波動,能夠在沒有先驗信息的情況下反演對流層下部的溫度和水汽廓線。目前,主要有第二類正則變換方法(CT2)、全譜法(FSI)等波動光學反演方法。Gorbunov等利用波動光學方法初步研究了22.23 GHz吸收線附近頻點反演溫度和水汽廓線的性能,仍需要進一步開展研究,比較不同波動光學方法對反演性能的影響。此外,還需對波動光學方法在183.31 GHz吸收線附近頻點反演溫度和水汽廓線的性能進行評估分析。

    LEO-LEO微波掩星理論上可以實現(xiàn)云中冰晶的探測,但目前關于這方面的研究非常少,在未來工作中,需要從基本原理出發(fā),研究云中冰晶的探測和反演技術,從而實現(xiàn)LEO-LEO微波掩星技術對水汽、云中液態(tài)水和云中冰晶三相態(tài)的探測和反演。

    ②硬件載荷研制

    國外已研制出低頻段22.23 GHz和高頻段183.31 GHz的原理樣機,并進行了地面實驗驗證,有待進行星載探測驗證[24,25];國內(nèi)硬件研制領域,低頻器件及相關技術較高頻器件成熟,目前已經(jīng)突破X/K波段高穩(wěn)定功率發(fā)射與高精度接收技術、X/K波段天線幅度穩(wěn)定度控制技術等關鍵技術研究,研制出了低頻段22.23 GHz的原理樣機,未來將進行外場、山基等實驗驗證。但尚缺乏183.31 GHz原理樣機相關的研制,主要原因是高頻振蕩器等方面的關鍵技術有待攻關突破。

    ③ LEO-LEO微波掩星技術與LEO-LEO紅外激光掩星技術協(xié)同研究

    LEO-LEO微波掩星技術利用微波信號的折射作用反演大氣的溫度、壓強等大氣熱力學參數(shù)廓線,結合水汽吸收頻點的吸收作用反演水汽廓線;再以微波掩星探測的大氣熱力學參數(shù)作為輸入量,利用各種氣體成分對相應吸收譜線頻點的吸收作用反演氣體濃度廓線,該過程既是LEO-LEO微波掩星和LEO-LEO紅外激光掩星的同步探測過程,也是協(xié)同反演的過程,從而實現(xiàn)大氣溫度、濕度和壓強等熱力學參數(shù)以及H2O、CO2、CH4、N2O、O3和CO等多種氣體成分的高精度、高垂直分辨率探測。因此,在開展上述總結的LEO-LEO微波掩星技術研究工作的同時,需要開展LEO-LEO紅外激光掩星的研究工作,主要包括多通道差分吸收激光雷達技術、短波紅外波段高精度高靈敏度探測技術等關鍵技術的攻關研究。

    6 結束語

    LEO-LEO掩星探測技術是目前無線電掩星技術最有前途的發(fā)展方向之一。LEO-LEO微波掩星探測概念自提出以來,已經(jīng)開展了一系列有意義的研究,建立了LEO-LEO微波掩星反演算法,開展了系列端到端的仿真計算和反演實驗,進行了載荷原理樣機的地面驗證試驗,初步證明了LEO-LEO掩星探測地球大氣的可行性和探測能力[5,6,8,12,14,36]。本文重點對LEO-LEO微波掩星探測溫度和水汽廓線的研究進展進行了總結,并對未來LEO-LEO微波掩星探測技術研究進行了展望,為發(fā)展LEO-LEO掩星空間探測技術提供了基本科學理論參考。

    [1] 孫穎. 人類活動對氣候系統(tǒng)的影響——解讀IPCC第六次評估報告第一工作組報告第三章[J]. 大氣科學學報, 2021, 44(5): 654–657.

    SUN Ying. Impact of humanactivities on climate system: an interpretation of Chapter Ⅲof WGⅠreport of IPCC AR6[J]. Transactions of Atmospheric Sciences, 2021, 44(5): 654–657.

    [2] 王鑫, 薛震剛, 杜曉勇, 等. 利用無線電掩星觀測數(shù)據(jù)反演大氣水汽剖面[J]. 電波科學學報, 2003, 18(4): 462–465.

    WANG Xin, XUE Zhenggan, DU Xiaoyong, et al. Inversion of atmospheric water vapour using radio occultation data[J]. Chinese Journal of Radio Science, 2003, 18(4): 462–465.

    [3] WATS. Water vapor and temperature in the troposphere and stratosphere[R]. Noordwijk: ESA/ESTEC, 2001.

    [4] 盧乃錳, 谷松巖. 氣象衛(wèi)星發(fā)展回顧與展望[J]. 遙感學報, 2016, 20(5): 832–841.

    LU Naimeng, GU Songyan. Review and prospect on the development of meteorological satellites[J]. Journal of Remote Sensing, 2016, 20(5): 832–841.

    [5] H?EG P, KIRCHENGAST G. ACE+ atmosphere and climate explorer based on GPS, GALILEO, and LEO-LEO radio occultation[R]. Austria: Institute for Geophysics, Astrophysics, Meteorology, University of Graz, 2002.

    [6] KURSINSKI E R, WARD D, OTAROLA A, et al. The active temperature, ozone and moisture microwave spectrometer (ATOMMS)[M]// STEINER A K, et al. New Horizons in Occultation Research. Berlin: Springer, 2009: 295–313.

    [7] WONG S, FETZER E J, SCHREIER M, et al. Cloud induced uncertainties in AIRS and ECMWF temperature and specific humidity[J]. Journal of Geophysical Research-Atmosphere, 2015, 120(5): 1880–1901.

    [8] SCHWEITZER S. Atmosphere and climate explorer mission ACE+: Humidity and temperature retrieval performance analysis[D]. Austria: Institute for Geophysics, Astrophysics, Meteorology, University of Graz, 2004.

    [9] KURSINSKI E R, HAJJ G A, SCHOFIELD J T, et al. Observing earth's atmosphere with radio occultation measurements using the global positioning system[J]. Journal of Geophysical Research, 1997, 102(D19): 23429–23465.

    [10] YUNCK T P, HAJJ G A, KURSINSKI E R, et al. AMORE: An autonomous constellation concept for atmospheric and ocean observation[J]. Acta Astronautica, 2000, 46(2–6): 355–364.

    [11] FENG D D, SYNDERGAARD S. HERMAN B M, et al. Deriving atmospheric water vapor and ozone profiles from active microwave occultation measurements[C]// FUJISADA H, LURIE J B, ROPERTZ A, et al. Sensors, Systems, and Next-Generation Satellites IV, Proc. SPIE, 2001, 4169: 299–308.

    [12] KURSINSKI E R, SYNDERGAARD S, FLITTNER D, et al. A microwave occultation observing system optimized to characterize atmospheric water, temperature, and geopotential via absorption[J]. Journal of Atmospheric and Oceanic Technology, 2002, 19(12): 1897–1914.

    [13] HERMAN B M, KURSINSKI E R, FENG D, et al. Active Tropospheric Ozone and Moisture Sounder (ATOMS)[R]. Tucson: Institute of Atmospheric Physics, University of Arizona, 2003.

    [14] HERMAN B M, FENG D, FLITTNER D, et al. An overview of the Univ. of Arizona ATOMS project[M]// KIRCHENGAST G, FOELSCHE U, STEINER A K. Occultations for Probing Atmosphere and Climate. Berlin: Springer, 2004: 189–200.

    [15] 樂新安, 郭英華, 曾楨, 等. 近地空間環(huán)境的GNSS無線電掩星探測技術[J]. 地球物理學報, 2016, 59(4): 1161–1188.

    YUE Xin’an, GUO Yinghua, ZENG Zhen, et al. GNSS radio occultation technique for near-Earth space environment detection[J]. Chinese Journal of Geophysics, 2016, 59(4): 1161–1188.

    [16] FJELDBO G, ESHLEMAN V R. The bistatic radar-occultation method for the study of planetary atmospheres[J]. Journal of Geophysical Research, 1965, 70(13): 3217–3225.

    [17] FJELDBO G, KLIORE A J, ESHLEMAN V R. The neutral atmosphere of Venus as studied with the Mariner V radio occultation experiments[J]. The Astronomical Journal, 1971, 76(2): 123.

    [18] KURSINSKI E R, HAJJ G A, HARDY, et al. Observing tropospheric water vapor by radio occultation using the global positioning system[J]. Geophysical Research Letters, 1995, 22(17): 2365–2368.

    [19] 王鑫, 薛震剛, 杜曉勇, 等. GPS-LEO掩星探測地球大氣的方法[J]. 全球定位系統(tǒng), 2002, 27(5): 14–20.

    WANG Xin, XUE Zhengang, DU Xiaoyong, et al. Method of GPS-LEO occultation measurements of the atmosphere[J]. GPS World of China, 2002, 27(5): 14–20.

    [20] POLI P, JOINER J, KURSINSKI E R. 1DVAR analysis of temperature and humidity using GPS radio occultation refractivity data[J]. Journal of Geophysical Research, 2002, 107(D20): 4448.

    [21] 王鑫, 呂達仁, 薛震剛. GNSS掩星中大氣水汽的非線性反演[J]. 地球物理學報, 2005, 48(1): 32–38.

    WANG Xin, LYU Daren, XUE Zhengang. A non-linear inversion method for retrieval of water vapour from radio occultation measurements[J]. Chinese Journal of Geophysics, 2005, 48(1): 32–38.

    [22] LI Y, YUAN Y, WANG, X. Assessments of the retrieval of atmospheric profiles from GNSS radio occultation data in moist tropospheric conditions using radiosonde data[J]. Remote Sensing, 2020, 12(17): 2717.

    [23] LIU C L, KIRCHENGAST G, SYNDERGAARD S, et al. A review of LEO-LEO occultation techniques using microwave and infrared-laser signals[J]. Advances in Space Research, 2017, 60(12): 2776–2811.

    [24] KURSINSKI E R, WARD D, STOVERN M, et al. Development and testing of the active temperature, ozone and moisture microwave spectrometer (ATOMMS) cm and mm wavelength occultation instrument[J]. Atmospheric Measurement Techniques, 2012, 5(2): 439–456.

    [25] KURSINSKI E R, WARD D, OTAROLA A, et al. Atmospheric profiling via satellite to satellite occultations near water and ozone absorption lines for weather and climate[C]// In: Proc SPIE 9881, Earth Observing Missions and Sensors: Development, Implementation, and Characterization IV, 2016: 98810Z.

    [26] ERIKSSON P, JIMENEZ C, MURTAGH D, et al. Measurement of tropospheric/stratospheric transmission at 10-35GHz for H2O retrieval in low Earth orbiting satellite links[J]. Radio Science, 2003, 38(4): 8069.

    [27] KIRCHENGAST G, BERNATH P, BUEHLER S, et al. ACCURATE-climate benchmark profiling of greenhouse gases and thermodynamic variables and wind from space[R]. Austria: Wegener Center, University of Graz, 2010.

    [28] KIRCHENGAST G, SCHWEITZER, S. Climate benchmark profiling of greenhouse gases and thermodynamic structure and wind from space[J]. Geophysical Research Letters, 2011, 38(13): L13701.

    [29] PROSCHEK V, KIRCHENGAST G, SCHWEITZER S, et al. Retrieval and validation of carbon dioxide, methane and water vapor for the Canary Islands IR-laser occultation experiment[J]. Atmospheric Measurement Techniques, 2015, 8(8): 3315–3336.

    [30] 杜曉勇, 符養(yǎng), 薛震剛, 等. 衛(wèi)星軌道參數(shù)對LEO-LEO掩星事件數(shù)量及分布影響的模擬研究[J]. 地球物理學報, 2007, 50(5): 1289–1297.

    DU Xiaoyong, FU Yang, XUE Zhengang, et al. Simulation of the impacts of satellite orbit parameters on the number and distribution of LEO - LEO occultation events[J]. Chinese Journal of Geophysics, 2007, 50(5): 1289–1297.

    [31] 呂華平, 嚴衛(wèi), 王迎強, 等. LEO-LEO掩星事件仿真研究[J]. 遙感信息, 2015, 30(2): 65–70.

    LYU Huaping, YAN Wei, WANG Yingqiang, et al. Investigation of LEO-LEO occultation simulation[J]. Remote Sensing Information, 2015, 30(2): 65–70.

    [32] 呂華平, 嚴衛(wèi), 柳聰亮, 等. LEO-LEO掩星探測頻點選擇分析與仿真研究[J]. 氣象科學, 2016, 36(2): 212–217.

    LYU Huaping, YAN Wei, LIU Congliang, et al. Analysis and simulation of frequency selection of LEO-LEO occultation atmosphere detection[J]. Journal of tshe Meteorological Sciences, 2016, 36(2): 212–217.

    [33] 孫立剛, 吳小成, 胡雄. LEO-LEO掩星事件持續(xù)時間及衛(wèi)星間相對角速度仿真分析[J]. 空間科學學報, 2017, 37(5): 574–584.

    SUN Ligang, WU Xiaocheng, HU Xiong. Simulation of duration of LEO-LEO occultation events and relative angular velocity between satellites[J]. Chinese Journal of Space Science, 2017, 37(5): 574–584.

    [34] WU C, SUN Y, WANG X, et al. A LEO-LEO occultation system using microwave signals[C]//IGARSS 2018: 3147–3148.

    [35] 馮梓宸. LEO-LEO掩星信號接收技術研究[D]. 北京: 中國科學院大學(中國科學院國家空間科學中心), 2018.

    [36] SCHWEITZER S, KIRCHENGAST G, SCHWAERZ M, et al. Thermodynamic state retrieval from microwave occultation data and performance analysis based on end-to-end simulations[J]. Journal of Geophysical Research, 2011, 116(D10): D10301.

    [37] KOLBE W F, BUSCHER H, LESKOVAR B. Microwave absorption coefficients of atmospheric pollutants and constituents[J]. Quantitative Spectroscopy and Radiative Transfer, 1977, 18(1): 47–64.

    [38] KIRCHENGAST G, SCHWEITZER S, SCHWAERZ M, et al. Advanced retrieval processing chain for derivation of atmospheric profiles from LEO-LEO radio occultation data[R]. Austria: Wegener Center, University of Graz, 2006.

    [39] SYNDERGAARD S, FLITTNER D E, KURSINSKI E R, et al. Simulating the influence of horizontal gradients on retrieved profiles from ATOMS occultation measurements–a promising approach for data assimilation[M]// KIRCHENGAST G, FOELSCHE U, STEINER A K. Occultations for Probing Atmosphere and Climate. Berlin: Springer, 2004: 221–232.

    [40] GORBUNOV M E, KIRCHENGAST G. Advanced wave-optics processing of LEO-LEO radio occultation data in presence of turbulence[R]. Austria: Institute for Geophysics, Astrophysics, and Meteorology, University of Graz, 2005a.

    [41] GORBUNOV M E, KIRCHENGAST G. Processing X/K band radio occultation data in the presence of turbulence[J]. Radio Science, 2005b, 40(6): RS6001.

    [42] GORBUNOV M E, KIRCHENGAST G. Fluctuations of radio occultation signals in X/K band in the presence of anisotropic turbulence and differential transmission retrieval performance[J]. Radio Science, 2007, 42(4): RS4025.

    [43] KURSINSKI E R, FENG D, FLITTNER D, et al. An active microwave limb sounder for profiling water vapor, ozone, temperature, geopotential, clouds, isotopes and stratospheric winds[M]// KIRCHENGAST G, FOELSCHE U, STEINER A K. Occultations for Probing Atmosphere and Climate. Berlin: Springer, 2004: 173–187.

    Progress of temperature and water vapor profiles detected by LEO-LEO microwave occultation

    ZHANG Zhihua1,2, WANG Xin1, LYU Daren1,2

    (1. Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China;2. University of Chinese Academy of Sciences, Beijing 100049, China)

    LEO-LEO microwave occultation measurement technique uses the frequencies near the water vapor absorption line to measure both refraction and absorption of the signals, enabling simultaneous retrieval of the temperature and water vapor profiles without an auxiliary background field, which overcomes the "temperature-humidity ambiguity" problem of GNSS-LEO radio occultation technique in the middle and low troposphere. At present, there is no any LEO-LEO occultation exploration mission in the space. In order to promote of LEO-LEO occultation measurement technique, the research frontiers and progress of temperature and water vapor profile detected by LEO-LEO microwave occultation are reviewed and summarized in this paper. Firstly, the development history of LEO-LEO occultation measurement technique is introduced, and the basic principle of LEO-LEO occultation detection and the frequencies selection are explained. Secondly, the research trends of retrieval technique of LEO-LEO microwave occultation detection temperature and water vapor profiles is emphasized, and the progress of LEO-LEO occultation missions is introduced. Finally the paper summarizes and offers a prospect. This study can provide basic scientific and theoretical reference for the realization of LEO-LEO occultation space exploration mission.

    LEO-LEO; Microwave occultation measurement; Retrieval; Temperature profiles; Water vapor profiles

    P412.27

    A

    CN11-1780(2022)01-0001-12

    10.12347/j.ycyk.20211109001

    張志華, 王鑫, 呂達仁.LEO-LEO微波掩星探測溫度和水汽廓線研究進展[J]. 遙測遙控, 2022, 43(1): 1–12.

    DOI:10.12347/j.ycyk.20211109001

    : ZHANG Zhihua, WANG Xin, LYU Daren. Progress of temperature and water vapor profiles detected by LEO-LEO microwave occultation[J]. Journal of Telemetry, Tracking and Command, 2022, 43(1): 1–12.

    張志華 1990年生,在讀博士研究生,主要研究方向為LEO-LEO微波掩星探測技術。

    王 鑫 1977年生,博士,高級工程師,主要研究方向為大氣物理、大氣與海洋遙感。

    呂達仁 1940年生,院士,研究員,大氣物理學家。

    Website: ycyk.brit.com.cn Email: ycyk704@163.com

    基金項目:中科院戰(zhàn)略性先導科技專項空間科學背景型號項目(XDA15021000)

    通訊作者:王鑫(wangx2003@mail.iap.ac.cn)

    2021-11-09

    2021-12-16

    (本文編輯:傅 杰)

    猜你喜歡
    掩星廓線水汽
    青藏高原上空平流層水汽的時空演變特征
    FY-3D 衛(wèi)星的北斗掩星分布特征與誤差特性*
    基于COSMIC掩星精密定軌數(shù)據(jù)的等離子體層電子含量研究
    利用掩星溫度數(shù)據(jù)推算大氣月平均緯向風場
    不同降水強度下風廓線雷達譜矩特征與測風準確性分析
    同化Argo海洋廓線觀測對CAS-ESM-C的上層海洋溫度模擬的改進
    1979~2011年間平流層溫度及平流層水汽的演變趨勢
    基于快速局域線性回歸的IRAS/FY-3B大氣溫濕廓線反演
    GRAPES模式頂外部背景廓線構造方法初步研究
    深圳“5·11”特大暴雨過程的水汽輸送特征分析
    日韩精品有码人妻一区| 精品99又大又爽又粗少妇毛片| 欧美最黄视频在线播放免费| 久久精品91蜜桃| 日韩欧美精品v在线| 哪里可以看免费的av片| 1024手机看黄色片| 亚洲av中文字字幕乱码综合| 精品熟女少妇av免费看| 亚洲国产欧美在线一区| 春色校园在线视频观看| 十八禁国产超污无遮挡网站| 亚洲中文字幕日韩| 精品午夜福利在线看| 最后的刺客免费高清国语| 国产69精品久久久久777片| 午夜视频国产福利| 国产亚洲5aaaaa淫片| 精品久久国产蜜桃| 在线免费十八禁| 亚洲人成网站高清观看| 午夜福利在线观看免费完整高清在 | 麻豆久久精品国产亚洲av| 男人舔女人下体高潮全视频| a级毛片免费高清观看在线播放| 亚洲精品粉嫩美女一区| 亚洲美女视频黄频| 日本av手机在线免费观看| 美女 人体艺术 gogo| 99热精品在线国产| 日韩精品青青久久久久久| 国产成人福利小说| 国产日本99.免费观看| 久久这里有精品视频免费| 身体一侧抽搐| 国国产精品蜜臀av免费| 麻豆乱淫一区二区| 亚洲,欧美,日韩| 久久99精品国语久久久| 成人漫画全彩无遮挡| 秋霞在线观看毛片| 一个人看视频在线观看www免费| 日产精品乱码卡一卡2卡三| 欧美一区二区国产精品久久精品| 亚洲成a人片在线一区二区| 看黄色毛片网站| 在线观看av片永久免费下载| 美女xxoo啪啪120秒动态图| 在线免费观看不下载黄p国产| 美女被艹到高潮喷水动态| 99久久久亚洲精品蜜臀av| 国产私拍福利视频在线观看| 国产老妇伦熟女老妇高清| 联通29元200g的流量卡| 欧美又色又爽又黄视频| 91在线精品国自产拍蜜月| 嘟嘟电影网在线观看| 精品少妇黑人巨大在线播放 | 一边亲一边摸免费视频| 最近的中文字幕免费完整| 看片在线看免费视频| 欧美一区二区精品小视频在线| 一边亲一边摸免费视频| 麻豆av噜噜一区二区三区| 中国美白少妇内射xxxbb| 国产av在哪里看| 听说在线观看完整版免费高清| 美女 人体艺术 gogo| 久久久欧美国产精品| 此物有八面人人有两片| 天美传媒精品一区二区| 国内久久婷婷六月综合欲色啪| 亚洲经典国产精华液单| 五月玫瑰六月丁香| 韩国av在线不卡| 国产亚洲精品av在线| 人人妻人人澡欧美一区二区| 免费搜索国产男女视频| 别揉我奶头 嗯啊视频| av福利片在线观看| 亚洲国产欧洲综合997久久,| 人妻系列 视频| 尾随美女入室| 亚洲美女视频黄频| 国产精品嫩草影院av在线观看| 精品久久久噜噜| 1000部很黄的大片| 我要搜黄色片| 中文在线观看免费www的网站| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 久久99热6这里只有精品| 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 亚洲精品久久国产高清桃花| 亚洲在线观看片| 99久久久亚洲精品蜜臀av| 欧美性猛交黑人性爽| 秋霞在线观看毛片| 国语自产精品视频在线第100页| 色综合亚洲欧美另类图片| 一级毛片我不卡| 日本一本二区三区精品| 国产av在哪里看| 精品免费久久久久久久清纯| 又爽又黄a免费视频| 亚洲欧美精品自产自拍| 嫩草影院精品99| 久久精品国产99精品国产亚洲性色| 少妇熟女aⅴ在线视频| 日韩精品有码人妻一区| 一本久久精品| 国产精品不卡视频一区二区| 久久这里只有精品中国| 日本一二三区视频观看| 中文资源天堂在线| 99久久精品热视频| 中文字幕精品亚洲无线码一区| 欧美日本视频| 久久这里只有精品中国| 成人性生交大片免费视频hd| 国产成人freesex在线| 亚洲美女视频黄频| 99久国产av精品国产电影| 久久这里有精品视频免费| 久久99精品国语久久久| 99国产精品一区二区蜜桃av| 99精品在免费线老司机午夜| www日本黄色视频网| 最后的刺客免费高清国语| 99精品在免费线老司机午夜| 26uuu在线亚洲综合色| av女优亚洲男人天堂| 国内精品久久久久精免费| 22中文网久久字幕| 日韩成人av中文字幕在线观看| 中文字幕制服av| 成人特级黄色片久久久久久久| 国产国拍精品亚洲av在线观看| 亚洲人与动物交配视频| 搞女人的毛片| 亚洲国产欧洲综合997久久,| 国产日本99.免费观看| 国产精品.久久久| 久久久成人免费电影| 久久午夜亚洲精品久久| 国产精品99久久久久久久久| 国产精品免费一区二区三区在线| 精品无人区乱码1区二区| 国产成人freesex在线| 免费av毛片视频| av女优亚洲男人天堂| 欧美日韩国产亚洲二区| 深夜精品福利| 精品日产1卡2卡| 精品不卡国产一区二区三区| 亚洲无线在线观看| 国产乱人偷精品视频| 99久久久亚洲精品蜜臀av| av在线天堂中文字幕| 免费黄网站久久成人精品| 久久精品91蜜桃| 噜噜噜噜噜久久久久久91| 成人亚洲精品av一区二区| 熟女电影av网| 搡老妇女老女人老熟妇| av黄色大香蕉| 亚洲成人精品中文字幕电影| 精品欧美国产一区二区三| 日韩在线高清观看一区二区三区| 高清日韩中文字幕在线| 免费看av在线观看网站| 国产精品一区二区性色av| 欧美性感艳星| 国产色婷婷99| 少妇丰满av| 禁无遮挡网站| 国产毛片a区久久久久| 免费电影在线观看免费观看| 极品教师在线视频| 高清午夜精品一区二区三区 | 色综合色国产| 国产成人福利小说| 国产在视频线在精品| 91午夜精品亚洲一区二区三区| 最近中文字幕高清免费大全6| 日韩欧美国产在线观看| 国产极品天堂在线| 国产亚洲av嫩草精品影院| 日韩欧美在线乱码| 国产成年人精品一区二区| 亚洲成a人片在线一区二区| 在线a可以看的网站| 少妇的逼好多水| 男人狂女人下面高潮的视频| 久久热精品热| 久久久久九九精品影院| 国内精品一区二区在线观看| a级毛片免费高清观看在线播放| 精品久久久久久成人av| 国产亚洲5aaaaa淫片| 性色avwww在线观看| 日韩欧美在线乱码| 久久精品国产清高在天天线| 干丝袜人妻中文字幕| 女同久久另类99精品国产91| 乱人视频在线观看| 国产成人精品一,二区 | 中文字幕久久专区| 美女大奶头视频| 在线观看美女被高潮喷水网站| 桃色一区二区三区在线观看| 精品久久久久久久久久免费视频| 国产精品国产高清国产av| 自拍偷自拍亚洲精品老妇| 欧美一区二区国产精品久久精品| 给我免费播放毛片高清在线观看| 国产91av在线免费观看| www.av在线官网国产| 欧美日韩在线观看h| 一级二级三级毛片免费看| 麻豆乱淫一区二区| 久久久欧美国产精品| 国产三级中文精品| 色综合站精品国产| 国产人妻一区二区三区在| 国产探花在线观看一区二区| 不卡一级毛片| 成人特级av手机在线观看| 秋霞在线观看毛片| 看免费成人av毛片| a级毛片a级免费在线| 国产乱人偷精品视频| 日本与韩国留学比较| 国内精品宾馆在线| 成人亚洲欧美一区二区av| 美女xxoo啪啪120秒动态图| 国产成人a区在线观看| 欧美日韩精品成人综合77777| 国产人妻一区二区三区在| 国产黄片美女视频| 狠狠狠狠99中文字幕| 亚洲丝袜综合中文字幕| 亚洲av免费在线观看| 中出人妻视频一区二区| av天堂在线播放| 成人无遮挡网站| 日韩成人伦理影院| 一级毛片电影观看 | 中文字幕久久专区| 久久婷婷人人爽人人干人人爱| 日韩一区二区三区影片| 日韩av不卡免费在线播放| 日韩欧美精品免费久久| 少妇高潮的动态图| 观看美女的网站| 国产探花在线观看一区二区| 亚洲精品色激情综合| 一卡2卡三卡四卡精品乱码亚洲| 男女啪啪激烈高潮av片| 在现免费观看毛片| 亚洲人成网站在线播放欧美日韩| 嫩草影院入口| 欧美激情久久久久久爽电影| 欧美一区二区亚洲| 午夜免费男女啪啪视频观看| 国产精品一区二区三区四区久久| 亚洲人成网站在线播| 久久欧美精品欧美久久欧美| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 欧美zozozo另类| 在线播放无遮挡| 国产成年人精品一区二区| 成年av动漫网址| 国产极品天堂在线| 午夜激情福利司机影院| 中文字幕av成人在线电影| 亚洲精品国产av成人精品| 久久久久国产网址| 一级毛片久久久久久久久女| 日日摸夜夜添夜夜添av毛片| 成人特级av手机在线观看| 亚洲人成网站高清观看| 日韩欧美一区二区三区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品成人久久久久久| 久久精品夜色国产| 国内少妇人妻偷人精品xxx网站| 国模一区二区三区四区视频| 欧美最新免费一区二区三区| 免费人成视频x8x8入口观看| 日韩欧美在线乱码| 丰满乱子伦码专区| av专区在线播放| 美女黄网站色视频| 直男gayav资源| 婷婷色综合大香蕉| 亚洲欧美精品综合久久99| 蜜臀久久99精品久久宅男| 麻豆国产av国片精品| 女同久久另类99精品国产91| 国产精品久久视频播放| 亚洲色图av天堂| 青青草视频在线视频观看| 少妇的逼水好多| 深夜精品福利| 天堂网av新在线| 国产伦精品一区二区三区视频9| 成人二区视频| 亚洲自拍偷在线| 午夜激情欧美在线| 国产女主播在线喷水免费视频网站 | 成人国产麻豆网| 国产单亲对白刺激| 国产一区二区三区av在线 | 九九久久精品国产亚洲av麻豆| 国产v大片淫在线免费观看| av在线观看视频网站免费| 能在线免费观看的黄片| 12—13女人毛片做爰片一| 久久久久九九精品影院| 亚洲国产高清在线一区二区三| 美女内射精品一级片tv| 国产老妇女一区| av卡一久久| 天天一区二区日本电影三级| 国产高清有码在线观看视频| 久久久久国产网址| 少妇猛男粗大的猛烈进出视频 | 禁无遮挡网站| 不卡一级毛片| 亚洲国产精品成人久久小说 | 成人鲁丝片一二三区免费| 91午夜精品亚洲一区二区三区| av女优亚洲男人天堂| 不卡一级毛片| 国产精品久久久久久精品电影小说 | 日韩av不卡免费在线播放| 国产淫片久久久久久久久| 夜夜看夜夜爽夜夜摸| 日韩欧美精品免费久久| 亚洲电影在线观看av| 欧美xxxx黑人xx丫x性爽| 大香蕉久久网| 美女内射精品一级片tv| 国产老妇女一区| 久久午夜福利片| 老熟妇乱子伦视频在线观看| 91在线精品国自产拍蜜月| 免费看光身美女| 国产黄色小视频在线观看| 欧美+日韩+精品| 中国国产av一级| 看黄色毛片网站| 九九热线精品视视频播放| 精品久久久久久久人妻蜜臀av| 精品久久国产蜜桃| 免费av观看视频| 久久久精品94久久精品| 欧美高清成人免费视频www| 男女边吃奶边做爰视频| 小蜜桃在线观看免费完整版高清| 国产一区二区三区在线臀色熟女| 看黄色毛片网站| 国产精品爽爽va在线观看网站| 麻豆一二三区av精品| 久久久久久久午夜电影| 国产精品人妻久久久久久| 国产视频首页在线观看| 欧美人与善性xxx| 观看美女的网站| 国产成人freesex在线| 亚洲av第一区精品v没综合| 国产精品99久久久久久久久| 欧美zozozo另类| 网址你懂的国产日韩在线| 亚洲国产精品成人综合色| 亚洲精华国产精华液的使用体验 | 亚洲不卡免费看| 精品一区二区三区人妻视频| 免费大片18禁| 国产高清三级在线| 日产精品乱码卡一卡2卡三| 午夜老司机福利剧场| 国产伦理片在线播放av一区 | 国产淫片久久久久久久久| 婷婷亚洲欧美| 噜噜噜噜噜久久久久久91| 中文在线观看免费www的网站| 久久99蜜桃精品久久| 禁无遮挡网站| 亚洲国产精品合色在线| 久久久久久国产a免费观看| 欧美高清性xxxxhd video| 亚洲激情五月婷婷啪啪| 久久午夜福利片| 欧美+日韩+精品| 国产精品日韩av在线免费观看| av黄色大香蕉| 少妇被粗大猛烈的视频| 日本撒尿小便嘘嘘汇集6| 在线天堂最新版资源| 日本熟妇午夜| 最好的美女福利视频网| 亚洲va在线va天堂va国产| 欧美不卡视频在线免费观看| 亚洲av一区综合| 精品国内亚洲2022精品成人| 欧美日韩在线观看h| 在现免费观看毛片| 小说图片视频综合网站| 超碰av人人做人人爽久久| 午夜精品一区二区三区免费看| 国产 一区精品| 国产亚洲精品久久久com| 91狼人影院| 我的女老师完整版在线观看| 国产一级毛片在线| 午夜福利高清视频| 一区二区三区免费毛片| 尾随美女入室| 亚洲内射少妇av| 亚洲天堂国产精品一区在线| 欧美日韩在线观看h| 成人av在线播放网站| 黄色日韩在线| 最近中文字幕高清免费大全6| 波多野结衣巨乳人妻| 国产一区亚洲一区在线观看| 九九爱精品视频在线观看| 日本与韩国留学比较| 狂野欧美白嫩少妇大欣赏| 1000部很黄的大片| 一本久久精品| 亚洲av二区三区四区| 免费人成视频x8x8入口观看| 天堂中文最新版在线下载 | 国产成人午夜福利电影在线观看| 免费看美女性在线毛片视频| 亚洲经典国产精华液单| 欧美色视频一区免费| 亚洲精品自拍成人| 国产v大片淫在线免费观看| 日韩制服骚丝袜av| 亚洲三级黄色毛片| 免费大片18禁| 小说图片视频综合网站| 国产人妻一区二区三区在| 久久99精品国语久久久| 欧洲精品卡2卡3卡4卡5卡区| 国产成人精品婷婷| 中国美白少妇内射xxxbb| 欧美在线一区亚洲| 国产成人精品久久久久久| 国产成人aa在线观看| 国产一区二区三区av在线 | 亚洲欧洲日产国产| 日本熟妇午夜| 美女cb高潮喷水在线观看| 日韩av不卡免费在线播放| 欧美一级a爱片免费观看看| 爱豆传媒免费全集在线观看| 中文字幕av在线有码专区| 精品免费久久久久久久清纯| 九九在线视频观看精品| 日韩 亚洲 欧美在线| 美女大奶头视频| av黄色大香蕉| 亚洲精品国产成人久久av| 日本熟妇午夜| 毛片女人毛片| 2022亚洲国产成人精品| 中文字幕免费在线视频6| 深夜精品福利| 欧美日韩国产亚洲二区| 97热精品久久久久久| 男的添女的下面高潮视频| 成人鲁丝片一二三区免费| 丰满的人妻完整版| 精品久久久久久久人妻蜜臀av| 国产黄片视频在线免费观看| 国产亚洲精品久久久久久毛片| 日韩 亚洲 欧美在线| 18禁黄网站禁片免费观看直播| 亚洲国产精品国产精品| 国产精品三级大全| 免费看光身美女| 男人的好看免费观看在线视频| 男人舔奶头视频| 欧美成人免费av一区二区三区| 久久精品国产清高在天天线| 免费av观看视频| 性插视频无遮挡在线免费观看| 国产成人精品久久久久久| 夜夜夜夜夜久久久久| 亚洲美女视频黄频| 天堂av国产一区二区熟女人妻| 岛国在线免费视频观看| 国产免费一级a男人的天堂| 国产69精品久久久久777片| 久久久久久久久大av| 国产午夜精品一二区理论片| 六月丁香七月| 亚洲精品亚洲一区二区| 国产精品永久免费网站| 久久久久久大精品| 波野结衣二区三区在线| 51国产日韩欧美| 久久精品夜夜夜夜夜久久蜜豆| 啦啦啦观看免费观看视频高清| 观看美女的网站| 综合色av麻豆| 国产精品国产高清国产av| 麻豆国产av国片精品| 97人妻精品一区二区三区麻豆| 国产亚洲av片在线观看秒播厂 | 久久久午夜欧美精品| 日本三级黄在线观看| 国产一区二区亚洲精品在线观看| 九草在线视频观看| 人妻少妇偷人精品九色| 国内揄拍国产精品人妻在线| 国产精品电影一区二区三区| 成人午夜高清在线视频| 欧美一区二区精品小视频在线| 老司机福利观看| 国产精品一区二区三区四区久久| 免费看光身美女| 1024手机看黄色片| 欧美bdsm另类| 欧美潮喷喷水| 免费人成视频x8x8入口观看| 日韩人妻高清精品专区| 99久久人妻综合| 久久久久久国产a免费观看| 国产白丝娇喘喷水9色精品| 成人性生交大片免费视频hd| 成人国产麻豆网| 波多野结衣巨乳人妻| 亚洲五月天丁香| 日韩av不卡免费在线播放| 日本在线视频免费播放| 丰满乱子伦码专区| 精品不卡国产一区二区三区| 夜夜夜夜夜久久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产黄色视频一区二区在线观看 | 国产成人影院久久av| 欧美日本亚洲视频在线播放| 亚洲av成人精品一区久久| 亚洲av第一区精品v没综合| 麻豆久久精品国产亚洲av| 欧美极品一区二区三区四区| 国产精品日韩av在线免费观看| 亚洲av二区三区四区| 免费在线观看成人毛片| 久久精品国产清高在天天线| 久久久久久久久久成人| 精品一区二区三区视频在线| 精品国产三级普通话版| 看黄色毛片网站| 性色avwww在线观看| 亚洲自拍偷在线| 成人综合一区亚洲| 男女视频在线观看网站免费| 亚洲精品久久久久久婷婷小说 | 国产精品不卡视频一区二区| 禁无遮挡网站| 黄片wwwwww| 赤兔流量卡办理| 成人无遮挡网站| 91aial.com中文字幕在线观看| 国内精品久久久久精免费| 少妇猛男粗大的猛烈进出视频 | 亚洲国产精品国产精品| 深爱激情五月婷婷| 国产极品精品免费视频能看的| 亚洲av不卡在线观看| 亚洲欧洲日产国产| 九九在线视频观看精品| 国产精品久久久久久久久免| 亚洲精品粉嫩美女一区| 亚洲成人av在线免费| 又爽又黄无遮挡网站| 国产亚洲91精品色在线| 国语自产精品视频在线第100页| 免费在线观看成人毛片| 五月玫瑰六月丁香| 看十八女毛片水多多多| 国产精品一二三区在线看| 亚洲av成人av| 免费无遮挡裸体视频| a级毛色黄片| 国产精品一区二区三区四区久久| 免费无遮挡裸体视频| 国产精品一二三区在线看| 亚洲熟妇中文字幕五十中出| 一个人看的www免费观看视频| 夜夜夜夜夜久久久久| 哪里可以看免费的av片| av天堂中文字幕网| 免费看av在线观看网站| 亚洲国产欧美人成| 国产精品久久久久久久电影| 欧美极品一区二区三区四区| 激情 狠狠 欧美| 亚洲一区二区三区色噜噜| 直男gayav资源| 日本av手机在线免费观看| 亚洲激情五月婷婷啪啪| 亚洲av成人精品一区久久| 国产精品久久久久久av不卡| 村上凉子中文字幕在线| 免费一级毛片在线播放高清视频|