孫詩雅黃富祥夏學(xué)齊

青藏高原衛(wèi)星紫外后向散射輻射的正演仿真

孫詩雅1，2，黃富祥2，3*，夏學(xué)齊1

（1.中國地質(zhì)大學(xué)（北京） 地球科學(xué)與資源學(xué)院，北京 100083；2.中國氣象局 國家衛(wèi)星氣象中心，北京 100081；3.中國氣象局 中國遙感衛(wèi)星輻射測量和定標(biāo)重點(diǎn)開放實(shí)驗(yàn)室，北京 100081）

為了定量評估像元下墊面反射率等要素對青藏高原衛(wèi)星臭氧遙感數(shù)據(jù)精度的影響，設(shè)計(jì)了青藏高原衛(wèi)星太陽紫外后向散射輻射的正演仿真計(jì)算。首先提取NOAA衛(wèi)星SBUV/2青藏高原像元下墊面的觀測數(shù)據(jù)，結(jié)合提取的臭氧垂直廓線數(shù)據(jù)，設(shè)計(jì)青藏高原正演仿真計(jì)算的輸入?yún)?shù)，最后輸入風(fēng)云三號衛(wèi)星紫外臭氧垂直探測儀正演仿真計(jì)算系統(tǒng)開展正演仿真計(jì)算，利用計(jì)算結(jié)果，考察衛(wèi)星紫外臭氧垂直探測儀通道輻射觀測值隨下墊面特征參數(shù)的變化特征。結(jié)果表明：當(dāng)太陽天頂角在30°～88°內(nèi)變化時(shí)，儀器各個(gè)通道的輻亮度觀測值都隨太陽天頂角增大而減小，變率在0.002%～46.24%，變率隨著波長的增大而增大，其中第12通道的最大達(dá)到46.24%；當(dāng)像元下墊面反射率在0～0.8內(nèi)變化時(shí)，儀器通道輻亮度觀測值的變化特征隨波長不同而有顯著差異，其中短波長的第1～6通道，對觀測值的影響幾乎可以忽略不計(jì)，而長波長的第7～12通道，觀測值隨著下墊面反射率的增大而增大，而且這種變率隨著波長的增大而快速增加，第12通道的最大變率達(dá)到2 903.2%；像元下墊面有效氣壓的改變，對儀器通道輻亮度觀測值的影響可以忽略不計(jì)。

衛(wèi)星遙感；青藏高原；紫外臭氧垂直探測儀；正演計(jì)算；輻射傳輸

1 引　言

青藏高原是全球最大的高原，被稱作世界第三極［1］，是我國“兩屏三帶”生態(tài)格局中重要的生態(tài)屏障區(qū)，也是北半球地帶臭氧損耗的重點(diǎn)地區(qū)。獨(dú)特的地理位置和海拔高度使青藏高原對全球氣候變化極其敏感，因此青藏高原也是全球臭氧和氣候變化研究的關(guān)鍵地區(qū)［2］。

高原臭氧變化監(jiān)測和研究主要依賴衛(wèi)星臭氧遙感技術(shù)。衛(wèi)星紫外臭氧遙感技術(shù)研究始于1970年美國首次在Nimbus-4衛(wèi)星上搭載紫外后向散射儀（Backscatter Ultraviolet， BUV）［3］。此后，在BUV 基礎(chǔ)上涌現(xiàn)出臭氧垂直探測儀（Solar Backscatter Ultraviolet， SBUV）和SBUV/2、臭氧總量探測儀（Total Ozone Mapping Spectrometer， TOMS）和（Ozone Monitoring Instrument， OMI）等一系列臭氧探測儀器，至今已連續(xù)獲得近40年的全球大氣臭氧資料［4］。我國首次自主研制的衛(wèi)星紫外臭氧遙感器（Solar Backscatter Ultraviolet Sounder， SBUS）和（Total Ozone Uni， TOU）于2008年由風(fēng)云三號A星攜載升空，此后FY-3B/C/D等后續(xù)星攜載同類儀器，共同組成臭氧連續(xù)觀測，在全球和極地臭氧監(jiān)測中發(fā)揮了重要作用［5］。

利用衛(wèi)星紫外臭氧探測儀監(jiān)測青藏高原大氣臭氧變化面臨的一個(gè)關(guān)鍵問題是，高原下墊面類型多樣、地形破碎化嚴(yán)重，像元下墊面反射率復(fù)雜多變，勢必引起衛(wèi)星臭氧反演數(shù)據(jù)的不確定性。然而，至今尚無人就高原獨(dú)特的地形地貌對衛(wèi)星臭氧反演造成的影響開展定量評估。

衛(wèi)星紫外臭氧反演算法由正演仿真計(jì)算和反演計(jì)算兩部分組成，其中正演計(jì)算是算法開發(fā)的第一步和基礎(chǔ)條件。正演計(jì)算是指在假設(shè)已經(jīng)獲取了像元臭氧垂直廓線和臭氧總量反演結(jié)果的條件下，構(gòu)造像元各種觀測條件作為輸入?yún)?shù)，利用紫外輻射傳輸模式，模擬計(jì)算衛(wèi)星臭氧垂直探測儀觀測的輻亮度觀測值，并考察通道觀測值隨著不同觀測參數(shù)變化而改變的特征和規(guī)律［6-7］。我國在研發(fā)FY-3 SBUS反演算法的過程中，利用SBUS的儀器通道參數(shù)，構(gòu)建了衛(wèi)星SBUS正演仿真系統(tǒng)（Solar Backscattered Ultraviolet Radiance Forward Calculation Simulating System， SBURFCSS），為本文開展青藏高原衛(wèi)星紫外后向散射輻射的正演仿真計(jì)算創(chuàng)造了基礎(chǔ)條件。

本文利用SBUV/2 L2級產(chǎn)品中提取的青藏高原衛(wèi)星SBUS的像元下墊面參數(shù)、云參數(shù)和臭氧廓線等信息，設(shè)計(jì)了衛(wèi)星紫外后向散射輻射正演模擬仿真的計(jì)算條件，利用風(fēng)云三號衛(wèi)星SBURFCSS正演仿真系統(tǒng)，開展正演仿真計(jì)算，定量評估高原像元下墊面特征對衛(wèi)星SBUS通道輻射觀測值的影響，為定量評估高原衛(wèi)星臭氧遙感產(chǎn)品精度創(chuàng)造有利條件。

2 數(shù)據(jù)和方法

2.1　衛(wèi)星紫外臭氧垂直探測儀

FY-3衛(wèi)星SBUS是我國第一次自主研制的衛(wèi)星臭氧探測儀器，在儀器結(jié)構(gòu)、通道設(shè)置和技術(shù)參數(shù)等方面，都與美國NOAA系列衛(wèi)星SBUV/2相同或相近。SBUS由探測器頭部和電箱組成，是一臺能在160～400 nm波段測量光譜輻亮度和太陽輻照度的掃描式光譜儀。其主要工作模式和功能有［8］：

（1）太陽模式：160～400 nm波段太陽連續(xù)光譜輻照度測量及250～340 nm間12個(gè)波長處太陽分立光譜輻照度測量；

（2）大氣模式：250～340 nm間12個(gè)特征波長處大氣太陽后向散射紫外光譜輻照度測量；

（3）標(biāo)準(zhǔn)燈模式：汞燈253.7 nm光譜線測量，用于儀器自身波長定標(biāo)，兩塊漫反射器分別進(jìn)行。

表1給出FY-3 SUBS與NOAA-17和18 SBUV/2的通道中心波長對比情況。對比表明，F(xiàn)Y-3 SBUS與NOAA SBUV/2觀測數(shù)據(jù)和反演產(chǎn)品具有良好的一致性［9-11］。FY-3 SBUS臭氧垂直廓線反演產(chǎn)品為反演的11層廓線，經(jīng)過插值處理生成從地表到大氣頂分為21層臭氧垂直廓線，量綱為DU，空間分辨率為200 km。在青藏高原地區(qū)，經(jīng)過大量臭氧探空和地基觀測數(shù)據(jù)檢驗(yàn)和評估，認(rèn)為FY-3 SBUS反演的臭氧總量和垂直廓線在青藏高原與地基觀測具有較好的一致性和適用性［12-13］。

表1FY-3 SUBS與NOAA-17和18衛(wèi)星SBUV/2通道中心波長比較

Tab.1　Comparison of central wavelengths between FY-3 SBUS and NOAA-17 and 18 SBUV/2 （nm）

本文利用提取的NOAA SBUV/2青藏高原像元觀測參數(shù)，構(gòu)造正演仿真計(jì)算條件，模擬計(jì)算FY-3 SBUS通道輻亮度觀測數(shù)據(jù)，考察像元下墊面參數(shù)變化對儀器觀測值的影響。

2.2　使用的數(shù)據(jù)

本文使用的衛(wèi)星SBUS青藏高原像元下墊面數(shù)據(jù)，從NOAA-16/17/18/19四顆衛(wèi)星SBUV/2 2009年11月至12月的L2產(chǎn)品中提取。NOAA在設(shè)計(jì)SBUV/2 L2產(chǎn)品格式時(shí)，保存了L1B像元的各種觀測參數(shù)，不僅有反演的臭氧垂直廓線信息，還有像元下墊面反射率、表面氣壓、儀器觀測幾何、先驗(yàn)和初估臭氧廓線等。因此，這里從SBUV/2 L2產(chǎn)品提取的高原像元信息包括：像元下墊面反射率、云蓋百分率、云頂氣壓、下墊面氣壓、反演臭氧垂直廓線、先驗(yàn)臭氧垂直廓線、先驗(yàn)溫度廓線、衛(wèi)星方位角和太陽天頂角等。利用這些參數(shù)構(gòu)建正演仿真計(jì)算條件。

青藏高原的地理位置選取27.5°N～37.5°N，75.5°E～104.5°E。在數(shù)據(jù)提取分析過程中，根據(jù)高原臭氧垂直廓線第一層數(shù)值為0的特征篩查剔除了少數(shù)非高原像元數(shù)據(jù)。

2.3　正演計(jì)算方法和FY-3 SBUS正演仿真系統(tǒng)SBURFCSS

衛(wèi)星SBUS正演計(jì)算，根據(jù)波長不同分為單次散射、單次散射和多次散射與反射兩種類型進(jìn)行計(jì)算［14］。對于波長小于290 nm的通道，由于臭氧的強(qiáng)烈吸收，光子無法到達(dá)30 km以下的大氣層［15］，這些通道的紫外后向散射輻射計(jì)算問題通過求解Rayleigh單次散射（Single Scattering，SS）得以解決［14］。而對于波長大于290 nm的通道，由于臭氧的吸收能力較弱，光子能夠進(jìn)入大氣深層到達(dá)對流層甚至地表，衛(wèi)星從大氣頂觀測到的太陽紫外后向散射輻射除了單次散射外，還包括進(jìn)入對流層深層受到大氣分子的多次散射輻射，以及從下墊面、氣溶膠和云等的表面反射輻射，因此，這些長波長通道的輻射計(jì)算需要處理多次散射和反射（Multiple Scattering and Reflection， MSR）問題。在產(chǎn)品生成算法中，SS部分可用數(shù)值積分精確計(jì)算，MSR部分則利用輻射查算表進(jìn)行快速插值計(jì)算［16-17］。

大氣層頂太陽紫外后向散射輻射中的SS部分可表示為［18］：

其中：F為波長處的太陽輻照度，β*為單位大氣壓的有效Rayleigh散射系數(shù)，（）是散射角為時(shí)的Rayleigh散射相函數(shù)，α*為單位質(zhì)量臭氧的有效吸收系數(shù)，S（）是氣壓層的散射輻射穿過的斜程臭氧質(zhì)量，S（）是斜程大氣質(zhì)量，（）是氣壓層以上臭氧氣柱含量，s是表面氣壓。

MSR部分的計(jì)算建立在有效朗伯面概念的基礎(chǔ)之上［19］。借鑒有效朗伯面概念，假設(shè)大氣只包含臭氧和Rayleigh散射粒子，沒有Mie散射，并限制在下墊面不透明的朗伯面上，這個(gè)假想面的有效反照率*和有效氣壓*通過一個(gè)臭氧吸收很弱的波長確定，忽略反照率隨波長的變化，認(rèn)為對于所有通道的下墊面反照率都等于有效反照率。在有效反照率為*的朗伯面上，包括單次散射和多次散射與反射輻射在內(nèi)的總輻射M可表示為：

其中：0，和S都是*，*，*，*以及太陽天頂角0，臭氧廓線（）和溫度廓線（）的函數(shù)［20］。

在獲得總輻射M、單次散射輻射SS計(jì)算結(jié)果的基礎(chǔ)上，通過二者相減得到多次散射和反射輻射，從而建立了總輻射、多次散射和反射輻射之間的對應(yīng)關(guān)系，利用衛(wèi)星觀測向量，通過線性插值計(jì)算，可以得到總輻射觀測值對應(yīng)的多次散射和反射輻射［21］。我們在研究開發(fā)FY-3 SBUS臭氧垂直廓線反演算法的過程中，通過輻射查算表方法來快速計(jì)算多次散射和反射輻射［22］。

FY-3 SBUS業(yè)務(wù)反演算法建立的衛(wèi)星SBUS正演仿真系統(tǒng)SBURFCSS，由儀器參數(shù)配置、像元和觀測幾何配置、單次散射計(jì)算、多次散射和反射輻射計(jì)算等部分組成。其中，儀器參數(shù)配置部分主要提供儀器通道中心波長、帶寬等參數(shù)；像元和觀測幾何部分則是在軌觀測時(shí)像元下墊面云蓋百分率、有效反射率、有效氣壓、太陽天頂角、方位角和衛(wèi)星觀測角等各種參數(shù)的配置和組合。這樣，在設(shè)計(jì)了各種輸入?yún)?shù)后，就可以通過輻射傳輸計(jì)算輸出衛(wèi)星載荷各通道的輻亮度觀測值。

2.4　正演計(jì)算條件設(shè)計(jì)

本文正演仿真計(jì)算，以提取的SBUV/2在青藏高原獲得的像元下墊面參數(shù)、觀測幾何、反演臭氧垂直廓線等為參考依據(jù)，設(shè)計(jì)如下正演計(jì)算條件：

（1）表面反射率取值5個(gè)水平：0，0.2，0.4，0.6，0.8；

（2）下墊面表面氣壓取值3個(gè)水平：60.6，55.55，50.5 kPa；

（3）云頂氣壓選取6個(gè)水平：80.8，75.75，60.6，50.5，45.45，40.4 kPa；

（4）有效云蓋百分率選取5個(gè)水平：0，0.2，0.4，0.6，0.9；

（5）太陽天頂角選取9個(gè)水平：30°，45°，60°，70°，77°，81°，84°，86°，88°；

（6）臭氧垂直廓線和大氣溫度廓線選取高原反演臭氧廓線和溫度廓線。

正演計(jì)算使用的衛(wèi)星探測器通道參數(shù)是表1中的FY-3 SBUS參數(shù)。計(jì)算輸出結(jié)果是SBUS 12個(gè)通道的輻亮度觀測值。

3 正演計(jì)算結(jié)果及分析

3.1　SBUS通道輻亮度觀測數(shù)據(jù)隨太陽天頂角的變化特征

圖1給出了SBUS 12個(gè)通道輻亮度觀測值隨太陽天頂角（Solar Zenith Angle， SZA）的變化情況。從圖1可見，SBUS輻亮度觀測值的變化特征如下：

圖1　青藏高原像元SBUS 12個(gè)通道輻亮度觀測值隨SZA的變化情況

（1）當(dāng)太陽天頂角在30°～88°內(nèi)變化時(shí)，儀器各通道輻亮度觀測值隨著太陽天頂角的增大而減小；

（2）不同波長的通道，輻亮度觀測值隨太陽天頂角的變化速率存在顯著差異，變率取值在0.001 7%～46.2%；波長越長，變率越大，第12通道的最大變率達(dá)到46.2%。

3.2　SBUS通道輻亮度觀測值隨像元表面反射率的變化特征

圖2給出了SBUS 12個(gè)通道輻亮度觀測值隨像元表面反射率的變化情況。從圖2可見，當(dāng)表面反射率在0～0.8內(nèi)變化時(shí)，儀器不同波長通道輻亮度觀測值呈現(xiàn)不同的特征。對波長較短的第1～6通道觀測值，像元反射率變化的影響可以忽略不計(jì)；而對波長較長的第7～12通道，其影響則比較顯著，通道觀測值變率在19.03%～2 900%。其中，第12通道受到的影響最大，達(dá)到2 900%。

短波長通道中光子被臭氧強(qiáng)烈吸收無法進(jìn)入大氣深層，更無法到達(dá)下墊面，因此下墊面反射率變化對這些通道的輻亮度觀測值不會產(chǎn)生影響；而波長超過290 nm的通道，由于臭氧吸收能力較弱，光子可以進(jìn)入大氣深層，尤其是在青藏高原，臭氧總量較少，更多光子可以達(dá)到像元下墊面被反射進(jìn)入儀器，當(dāng)下墊面反射率發(fā)生變化時(shí)，這些長波長通道的觀測值受到的影響更大。

圖2　青藏高原像元SBUS 12個(gè)通道輻亮度觀測值隨下墊面反射率的變化情況

3.3　SBUS通道輻亮度觀測值隨像元表面氣壓的變化特征

圖3給出了SBUS 12個(gè)通道輻亮度觀測值隨高原像元表面氣壓的變化情況。從圖3可見，像元下墊面表面氣壓變化對儀器各個(gè)通道輻亮度觀測值的影響都很小，可以忽略不計(jì)。

圖3　青藏高原像元SBUS 12個(gè)通道輻亮度觀測值隨像元表面氣壓的變化特征

3.4　像元要素對通道觀測值的綜合影響

實(shí)際中，高原像元下墊面不同要素往往是共同對儀器通道觀測值產(chǎn)生影響。上文分析表明，青藏高原像元表面反射率和太陽天頂角對儀器通道輻亮度觀測值的影響較大，而且隨著波長的改變，影響程度差異顯著，而表面氣壓對儀器觀測值的影響可以忽略不計(jì)。因此，本文考察像元表面反射率與太陽天頂角變化對儀器通道輻亮度觀測值的綜合影響。圖4給出像元表面反射率和太陽天頂角兩個(gè)因子對第1和第12通道輻亮度觀測值的綜合影響。

選取受到影響最小的第1通道和影響最大的第12通道作為例子，表現(xiàn)像元參數(shù)變化對儀器輻亮度觀測值的綜合影響。輻亮度觀測值主要隨著太陽天頂角增大而減小，表面反射率變化對觀測值幾乎不產(chǎn)生影響，因此，當(dāng)太陽天頂角為30°時(shí)，儀器輻亮度的觀測值取最大值，不隨下墊面反射率的變化而改變。在第12通道，當(dāng)太陽天頂角取30°、下墊面反射率取0.8，輻亮度觀測值最大；當(dāng)太陽天頂角固定，輻亮度觀測值隨著下墊面反射率的增大而增大；當(dāng)太陽天頂角取88°、下墊面反射率取0時(shí)，輻亮度觀測值最小。

圖4　青藏高原像元反射率和太陽天頂角變化對SBUS通道輻亮度觀測值的綜合影響

4 結(jié)　論

本文利用提取的NOAA衛(wèi)星SBUV/2像元數(shù)據(jù)，結(jié)合FY-3 SBUS儀器參數(shù)，利用在研發(fā)的衛(wèi)星SBUS正演仿真系統(tǒng)SBURFCSS，開展了青藏高原衛(wèi)星SBUS正演仿真計(jì)算。計(jì)算結(jié)果表明：太陽天頂角變化對SBUS各通道輻亮度觀測值產(chǎn)生顯著影響，這種影響隨著波長的增大顯著增大；當(dāng)太陽天頂角在30°～88°內(nèi)變化時(shí)，隨著太陽天頂角的增大，SBUS各通道輻亮度觀測值隨之減小。其中，第1通道的變率最小為0.002%，第12通道輻亮度變率最大為46.24%。像元下墊面反射率變化對SBUS通道輻亮度觀測值的影響隨波長不同存在顯著差異。短波長的第1～6通道幾乎不受影響；而長波長的第7～12通道，輻亮度觀測值隨著反射率的增加而增大，變率在19.03%～2 900%之間，而且變率隨著波長的增大快速增加，第12通道變率最大，達(dá)到2 900%。像元下墊面氣壓變化對SBUS各通道輻亮度觀測值的影響很小，可以忽略不計(jì)。

青藏高原地形高聳，臭氧總量相對較少，導(dǎo)致長波長通道光子更多地到達(dá)像元下墊面，下墊面反射率變化對這些通道輻亮度觀測值的影響遠(yuǎn)大于低海拔地區(qū)。此外，高原獨(dú)特的地貌導(dǎo)致像元下墊面破碎化，反射率多變，這勢必影響儀器長波長通道的觀測值，進(jìn)而影響高原衛(wèi)星臭氧反演結(jié)果的精度。

［1］ YAN J J， WANG G L， YANG P C. Study on the sensitivity of summer ozone density to the enhanced aerosol loading over the Tibetan Plateau［J］.， 2020， 11（2）： 138.

［2］ 萬凌峰，郭棟，劉仁強(qiáng)，等. WACCM3對夏季青藏高原臭氧谷雙心結(jié)構(gòu)的模擬性能評估［J］. 高原氣象， 2017， 36（1）： 57-66.

WAN L F， GUO D， LIU R Q，. Evaluation of WACCM3 performance on simulation of the double core of ozone valley over the Qinghai-Xizang in summer［J］.， 2017， 36（1）： 57-66.（in Chinese）

［3］ ZOU X L. Studies of FY-3 observations over the past 10 years： a review［J］.， 2021， 13（4）： 673.

［4］ STOLARSKI R S， KRUEGER A J， SCHOEBERL M R，. Nimbus 7 satellite measurements of the springtime Antarctic ozone decrease［J］.， 1986， 322（6082）： 808-811.

［5］ 王詠梅，王英鑒，王維和，等. FY-3衛(wèi)星紫外臭氧總量探測儀［J］. 科學(xué)通報(bào)， 2009， 54（23）： 3778-3783.

WANG Y M， WANG Y J， WANG W H，. FY-3 satellite ultraviolet total ozone unit［J］.， 2009， 54（23）： 3778-3783.（in Chinese）

［6］ WENG F Z， YU X W， DUAN Y H，. Advanced radiative transfer modeling system （ARMS）： a new-generation satellite observation operator developed for numerical weather prediction and remote sensing applications［J］.， 2020， 37（2）： 131-136.

［7］ 黃富祥，劉年慶，張鵬，等. 風(fēng)云三號A星紫外臭氧垂直探測儀反演產(chǎn)品的比較和評估［J］. 光學(xué)精密工程， 2010， 18（7）： 1568-1576.

HUANG F X， LIU N Q， ZHANG P，. Comparison and evaluation of retrieved ozone profiles from solar backscatter ultraviolet sounder on FY-3A satellite［J］.， 2010， 18（7）： 1568-1576.（in Chinese）

［8］ 王淑榮，李福田，宋克非，等. FY-3A氣象衛(wèi)星紫外臭氧垂直探測儀［J］. 光學(xué)學(xué)報(bào)， 2009， 29（9）： 2590-2593.

WANG S R， LI F T， SONG K F，. Ultrariolet ozone vertical profile probe for FY-3A meteorological satellite［J］.， 2009， 29（9）： 2590-2593.（in Chinese）

［9］ 陳樹. 臭氧廓線和紫外輻射觀測及對衛(wèi)星產(chǎn)品的驗(yàn)證分析［D］. 北京：中國氣象科學(xué)研究院， 2015.

CHEN SH.-［D］. Beijing： Chinese Academy of Meteorological Sciences， 2015. （in Chinese）

［10］ HUANG F X， HUANG Y， FLYNN L E，. Radiometric calibration of the solar backscatter ultraviolet sounder and validation of ozone profile retrievals［J］.， 2012， 50（12）： 4956-4964.

［11］ 黃富祥，黃煜， LAWRENCE E F，等. 紫外臭氧垂直探測儀（SBUS）輻射定標(biāo)和反演臭氧垂直廓線驗(yàn)證［J］. 氣象科技進(jìn)展，2013， 3（4）： 108-115.

HUANG F X， HUANG Y， FLYNN L，. Radiometric calibration of the solar backscatter ultraviolet sounder and validation of ozone profile retrievals［J］.， 2013， 3（4）： 108-115.（in Chinese）

［12］ 劉利，鄭向東，陳樹，等. 風(fēng)云-3A衛(wèi)星TOU臭氧總量地基對比驗(yàn)證分析［J］. 空間科學(xué)學(xué)報(bào)， 2015， 35（6）： 696-706.

LIU L， ZHENG X D， CHEN SH，. Validation study of FY-3A-TOU total ozone using ground-based data［J］.， 2015， 35（6）： 696-706.（in Chinese）

［13］ 黃富祥，劉年慶，趙明現(xiàn)，等. 風(fēng)云三號衛(wèi)星紫外臭氧垂直廓線產(chǎn)品反演試驗(yàn)［J］. 科學(xué)通報(bào)， 2009， 54（17）： 2556-2561.

HUANG F X， LIU N Q， ZHAO M X，. Vertical ozone profiles deduced from measurements of SBUS on FY-3 satellite［J］.， 2009， 54（17）： 2556-2561.（in Chinese）

［14］ DAVE J V. Investigation of the effect of atmospheric dust on the determination of total ozone from the earth's ultraviolet reflectivity measurements［C］.：7724692. Gaithersburg： Int Bus Mach Corp，1977.

［15］ WU T F， MA J S， YUAN R Z，. Single-scatter model for short-range ultraviolet communication in a narrow beam case［J］.， 2019， 31（3）： 265-268.

［16］ ADAMS C N， KATTAWAR G W. Radiative transfer in spherical shell atmospheres： I. Rayleigh scattering［J］.， 1978， 35（1）： 139-151.

［17］ KORKIN S， YANG E S， SPURR R，. Revised and extended benchmark results for Rayleigh scattering of sunlight in spherical atmospheres［J］.， 2020， 254： 107181.

［18］ SOX L， WICKWAR V B， YUAN T，. Comparison of Rayleigh-scatter and sodium resonance lidar temperatures［J］.， 2018， 176： 03005.

［19］ ZHOU Y， LIU D， PEITUO X U，. Detecting atmospheric-water optical property profiles with a polarized lidar［J］.， 2019， 23（1）： 108-115.

［20］ OKAMOTO H， SATO K， NISHIZAWA T，. Development of a multiple-field-of-view multiple-scattering polarization lidar： comparison with cloud radar［J］.， 2016， 24（26）： 30053-30067.

［21］ ZHANG K， ZHANG Y P， ZHOU Y D，. Multiple scattering effects in lidar and its potential application for water cloud detection［C］. 2019（）1720，2019，，， 2019： 928-935.

［22］ SATO K， OKAMOTO H. Application of single and multiple-scattering theories to analyses of space-borne cloud radar and lidar data［J］.， 2020： 1-13.

Forward simulation of solar backscattered ultraviolet radiation for satellite over the Tibetan plateau

SUN Shiya1，2，HUANG Fuxiang2，3*，XIA Xueqi1

（1，（），100083，；2，，100081，；3，，100081，），，：

To quantitatively evaluate the influence of elements such as reflectivity of the underlying pixel surface on the precision of ozone data from satellite remote sensing over the Tibetan plateau， a forward simulation calculation of solar backscattered ultraviolet radiation for a satellite over the Tibetan plateau was conducted. First， observation data of the underlying pixel surface from NOAA satellites SBUV/2 over the Tibetan plateau were extracted and combined with the extracted ozone vertical profile. The input parameters of the forward simulation calculation for the Tibetan plateau were adjusted. Then， the parameters were input to the forward simulation system of the FY-3 Solar Backscatter Ultraviolet Sounder for forward simulation calculation. Finally， we investigated the characteristics of the observation value changes from the Solar Backscatter Ultraviolet Sounder channel radiation with different parameters through the calculated results. The results show that， when the solar zenith angle varies in the range of 30°–88°， the radiance observation values for each channel of the instrument decrease with an increase in the solar zenith angle. The rate of variation ranges from0.002% to46.24%. The variability increases with an increase in wavelength， and the maximum value of the 12-th channel was46.24%. When the reflectivity of the underlying pixel surface changes within the range of 0-0.8， the variation characteristics of the observed radiance value from the instrument channel significantly vary with the wavelength. For the short-wavelength channels from 1 to 6， the variation of the channel observations was almost negligible. For the long-wavelength channels from 7 to 12， the observed value increases as the reflectivity of the underlying surface increases. Moreover， this rate of change increases rapidly as the wavelength increases， with a maximum variation of 2 903.2% for channel 12. The variation of the effective air pressure on the underlying pixel surface has almost negligible influence on the observation of the radiance from the instrument channel.

satellite remote sensing； Tibetan plateau； solar backscatter ultraviolet sounder； forward calculation； radiation transmission

P414

A

10.37188/OPE.20223003.0256

1004-924X（2022）03-0256-08

2021-04-08；

2021-04-30.

國家自然科學(xué)基金資助項(xiàng)目（No.41675031）

孫詩雅（19-），女，湖北襄陽人，碩士研究生，主要從事衛(wèi)星遙感大氣臭氧等方面的研究。E-mail： sunsy_rs@163.com

黃富祥（1967-），男，湖北蘄春人，博士，研究員，2000年于中國科學(xué)院獲得博士學(xué)位，2002年從中國科學(xué)院大氣物理研究所博士后出站，主要從事衛(wèi)星探測閃電和衛(wèi)星遙感大氣臭氧等方面的研究。 E-mail：huangfx@cma.gov.cn