潘張宇萱

摘要：針對(duì)星載差分吸收激光雷達(dá)，構(gòu)建了一套包含6個(gè)子模塊的正演模型，并圍繞卷云對(duì)探測(cè)能力的影響，開展了典型場(chǎng)景下的仿真模擬研究。通過對(duì)比薄冰云和厚冰云兩個(gè)典型場(chǎng)景的模擬結(jié)果，得出如下結(jié)論：當(dāng)場(chǎng)景中僅存在薄冰云時(shí)，可以有效探測(cè)云下CO₂；當(dāng)場(chǎng)景中存在較強(qiáng)的云層和氣溶膠時(shí)，難以有效探測(cè)云下CO₂，云和氣溶膠的衰減效應(yīng)對(duì)探測(cè)效果有顯著的影響。這些結(jié)果也說明，星載差分吸收激光雷達(dá)雖然空間分辨率極高，可以獲取垂直方向上的廓線信息，但與此同時(shí)，探測(cè)能力的精細(xì)化評(píng)估和數(shù)據(jù)的高精度處理技術(shù)仍是難點(diǎn)。

關(guān)鍵詞：星載激光雷達(dá)；差分吸收技術(shù)；CO₂；仿真模擬

一、前言

溫室氣體含量的逐年增加破壞了原來在天然溫室氣體成分下所維持的源和匯的自然平衡，其中，二氧化碳（CO₂）含量相對(duì)較大，也是全球碳循環(huán)的主要載體。被動(dòng)遙感中，地面碳柱觀測(cè)網(wǎng)絡(luò)TCCON是基于地基傅立葉紅外高光譜儀按照標(biāo)準(zhǔn)規(guī)范開展太陽紅外光譜觀測(cè)，反演CO₂等多種大氣溫室氣體成分[1]。2006年10月歐洲發(fā)射的METOP-A、2009年1月日本發(fā)射的GOSAT、2014年7月美國發(fā)射的OCO-2都是用于探測(cè)CO₂等溫室氣體含量的衛(wèi)星[2-4]。近年來，星載主動(dòng)遙感技術(shù)逐漸發(fā)展起來，如差分吸收激光雷達(dá)技術(shù)探測(cè)大氣主要?dú)怏w成分，美國NASA和歐洲ESA分別提出了A-SCOPE和ASENDS衛(wèi)星計(jì)劃，均以主動(dòng)式路徑積分差分吸收IPDA激光雷達(dá)為設(shè)計(jì)載荷[5-6]。2022年4月，我國成功發(fā)射大氣環(huán)境監(jiān)測(cè)衛(wèi)星，攜帶了全球第一臺(tái)用于探測(cè)二氧化碳的激光雷達(dá)，以期實(shí)現(xiàn)對(duì)地球大氣二氧化碳的全天時(shí)、高精度探測(cè)[7]。

星載差分吸收激光雷達(dá)在探測(cè)二氧化碳時(shí)會(huì)受到云層吸收和散射的影響，從而帶來探測(cè)結(jié)果的不確定性。本文將針對(duì)這一問題，建立正演模型，基于兩個(gè)典型冰云場(chǎng)景開展仿真模擬分析，探討冰云對(duì)于探測(cè)能力的影響。

二、探測(cè)原理

星載二氧化碳激光雷達(dá)采用差分吸收技術(shù)來實(shí)現(xiàn)有效探測(cè)，其基本原理如下：

三、正演模型

正演模型是指在儀器設(shè)計(jì)和探測(cè)模式等論證過程中，通過構(gòu)造數(shù)學(xué)模型，計(jì)算得到其理論值，然后用仿真模擬驗(yàn)證其設(shè)計(jì)的正確性和可行性，同時(shí)為設(shè)計(jì)方案提供優(yōu)化方法的一種有效手段。

本文構(gòu)建的正演模型由6個(gè)子模塊構(gòu)成，如圖1所示：典型場(chǎng)景模擬、光學(xué)特性計(jì)算、輻射傳輸計(jì)算、衰減后向散射系數(shù)計(jì)算、回波信號(hào)仿真、參數(shù)敏感性分析。將NCEP再分析數(shù)據(jù)輸入至WRF和WRF-CHEM中進(jìn)行三維場(chǎng)景模擬，并耦合初始設(shè)置的CO₂含量廓線；再利用OPAC數(shù)據(jù)庫對(duì)水凝物、氣溶膠的光學(xué)特性進(jìn)行計(jì)算，得到后向散射特性；結(jié)合LBLTRM模式，建立大氣輻射傳輸模式，計(jì)算衰減后向散射系數(shù)；根據(jù)衛(wèi)星平臺(tái)參數(shù)和激光雷達(dá)儀器參數(shù)，模擬探測(cè)器上獲得的仿真信號(hào)；根據(jù)仿真結(jié)果可對(duì)關(guān)鍵參數(shù)進(jìn)行敏感性分析。典型場(chǎng)景一旦確定，就需要以儀器、衛(wèi)星平臺(tái)參數(shù)作為輸入值，如表1所示。

四、典型模擬結(jié)果分析

（一）典型場(chǎng)景與CO₂廓線

作為正演模型的起點(diǎn)，我們選取了兩個(gè)典型場(chǎng)景。圖2為第一個(gè)場(chǎng)景的三維示例，模擬一次熱帶地區(qū)斯里蘭卡島上空的冰云。在該場(chǎng)景中，海洋（水體）的反照率低于0.1，陸地的反照率高于0.15，這對(duì)于靠地表反射的差分吸收技術(shù)來說很重要。在場(chǎng)景中耦合CO₂濃度時(shí)，陸地和海洋采用了不同的廓線。從全球統(tǒng)計(jì)來看，海洋上空CO₂含量略高于陸地。

（二）薄冰云場(chǎng)景模擬結(jié)果分析

第一個(gè)模擬場(chǎng)景緯度橫跨約4°，為典型熱帶地區(qū)薄冰云，云頂高度位于15km左右，冰晶含水量最大區(qū)域在7.0～7.5°N，均值達(dá)到約0.66mg/m³，垂直光學(xué)厚度基本在0.3以內(nèi)。該場(chǎng)景主要考察激光雷達(dá)穿透薄冰云層后對(duì)于CO₂的探測(cè)能力。

將場(chǎng)景中的各項(xiàng)大氣條件、水凝物、CO₂等背景參數(shù)輸入正演模型，得到的仿真結(jié)果如圖3所示。圖3（a）中15km處的明顯回波可以看到1572nm波長仍然能夠用于探測(cè)云層信息，且由于云層很薄，光學(xué)厚度小，云層下方的信號(hào)較為平滑，地表和水面處偶爾有強(qiáng)反射信號(hào)需要做過飽和處理。圖3（b）給出了λ_on與λ_off兩個(gè)通道探測(cè)信號(hào)的差值，平均值在14.3μV左右，且云測(cè)信息會(huì)被削弱，其差值越明顯，對(duì)后期CO₂濃度的反演越有利。取場(chǎng)景中每條廓線λ_on與λ_off通道信號(hào)的比值P_on/P_off，并在高度上對(duì)每條廓線進(jìn)行累加，從而考察CO₂吸收作用導(dǎo)致的信號(hào)變化。通過對(duì)比發(fā)現(xiàn)，盡管圖3（c）中P_on/P_off累加值在一定區(qū)間內(nèi)出現(xiàn)不穩(wěn)定，但其基本趨勢(shì)與圖3（d）中輸入正演模型的CO₂柱濃度一致，且相關(guān)性達(dá)到0.83，證明在薄冰云中，能夠?qū)O₂柱濃度實(shí)施有效探測(cè)，但信噪比需要通過參數(shù)優(yōu)化來提升。

（三）厚冰云場(chǎng)景模擬結(jié)果分析

第二個(gè)模擬選取了中緯度地區(qū)一次典型的陸地冰云與氣溶膠同時(shí)出現(xiàn)的場(chǎng)景。冰云位于8～10km之間，分布范圍廣闊，2km下方有氣溶膠分布，設(shè)置為“大陸污染型”氣溶膠。由于均為陸地，反射率均在0.2以上。該場(chǎng)景主要考察厚冰云對(duì)于CO₂探測(cè)能力的影響。

圖4（a）中8～10km處有明顯云層回波，35.3～35.7°N和36.4～36.7°N兩個(gè)區(qū)域內(nèi)的云層回波較強(qiáng)，光學(xué)厚度大于4.5，吸收作用顯著，其下方高度區(qū)域內(nèi)幾乎沒有信號(hào)。圖4（b）中，λ_on與λ_off兩個(gè)通道探測(cè)的信號(hào)差值較小，主要是受云層吸收作用的影響。圖4（c）中，將每條廓線的P_on/P_off在高度上進(jìn)行累加，在上述云層回波較強(qiáng)的兩個(gè)區(qū)域內(nèi)，P_on/P_off有很強(qiáng)的不確定性，幾乎無法分辨趨勢(shì)，但在其他云層較薄的區(qū)域，不確定性減弱。對(duì)比圖4（d），兩者基本沒有任何相似的特征，且相關(guān)系數(shù)也僅為0.09，表明在云層較厚、含水量較大時(shí)，無法對(duì)CO₂實(shí)施有效探測(cè)。

五、結(jié)語

本文針對(duì)星載差分吸收激光雷達(dá)，構(gòu)建了一套包含6個(gè)子模塊的正演模型，并圍繞卷云對(duì)探測(cè)能力的影響，開展了典型場(chǎng)景下的仿真模擬研究，并得出如下結(jié)論：當(dāng)場(chǎng)景中僅存在薄冰云時(shí)，可以有效探測(cè)云下CO₂；當(dāng)場(chǎng)景中存在較強(qiáng)的云層和氣溶膠時(shí)，難以有效探測(cè)云下CO₂，云和氣溶膠的衰減效應(yīng)對(duì)探測(cè)效果有顯著的影響。星載差分吸收激光雷達(dá)雖然空間分辨率極高，可以獲取垂直方向上的廓線信息，但與此同時(shí)，探測(cè)能力的精細(xì)化評(píng)估和數(shù)據(jù)的高精度處理技術(shù)仍是巨大挑戰(zhàn)。

參考文獻(xiàn)

[1]周敏強(qiáng)，張興贏，王普才，等.二氧化碳柱濃度的衛(wèi)星反演試驗(yàn)及地基驗(yàn)證[J].中國科學(xué)：地球科學(xué)，2015，45（6）：856-863.

[2]García OE，Sepúlveda E，Schneider M，et al.Consistency and quality assessment of the Metop-A/IASI and Metop-B/IASI operational trace gas products （O3， CO， N2O， CH4， and CO₂） in the subtropical North Atlantic[J].Atmospheric Measurement Techniques，2016，9（5）：2315-2333.

[3]Yokota T，Yoshida Y，Eguchi N，et al.Global concentrations of CO₂and CH4 retrieved from GOSAT： First preliminary results[J].SOLA，2009，5：160-163.

[4]Wunch D，Wennberg PO，Osterman G，et al.Comparisons of the orbiting carbon observatory-2 （OCO-2） XCO₂measurements with TCCON[J].Atmospheric Measurement Techniques，2017，10（6）：2209-2238.

[5]Caron J，Durand Y，Bezy JL，et al.Performance modeling for A-SCOPE： a space-borne lidar measuring atmospheric CO₂[C]//Lidar Technologies， Techniques， and Measurements for Atmospheric Remote Sensing V.SPIE，2009.

[6]Abshire JB，Riris H，Allan GR，et al.A lidar approach to measure CO₂concentrations from space for the ASCENDS Mission[C]//Lidar Technologies， Techniques， and Measurements for Atmospheric Remote Sensing VI.SPIE，2010.

[7]王敏，黃辛.大氣環(huán)境監(jiān)測(cè)衛(wèi)星升空[N].中國科學(xué)報(bào)，2022-04-18（1）.

作者單位：南京信息工程大學(xué)大氣物理學(xué)院