卷云對星載CO2激光雷達(dá)探測能力影響的仿真模擬研究

潘張宇萱

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

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

一、前言

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

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

二、探測原理

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

三、正演模型

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

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

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

（一）典型場景與CO₂廓線

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

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

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

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

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

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

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

五、結(jié)語

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

參考文獻(xiàn)

[1]周敏強(qiáng)，張興贏，王普才，等.二氧化碳柱濃度的衛(wèi)星反演試驗及地基驗證[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)測衛(wèi)星升空[N].中國科學(xué)報，2022-04-18（1）.

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