基于0级数据的星载大气痕量气体差分吸收光谱仪在轨衰变监测

苏静明，赵敏杰，周海金，杨东上,2，洪 炎，司福祺*

1.中国科学院合肥物质科学研究院安徽光学精密机械研究所环境光学与技术重点实验室，安徽 合肥 230031 2.中国科学技术大学，安徽 合肥 230026 3.安徽理工大学电气与信息工程学院，安徽 淮南 232001

引 言

卫星光学仪器在轨运行期间，受空间环境中太阳紫外辐射和原子氧剥蚀、空间辐射等影响，元器件在轨性能可能会发生衰变，这种衰变会改变原始光谱的形状、大小，进而对光谱定标、辐射定标精度产生影响，随着时间的推移这种影响会越来越大，最终影响L1b数据产品和反演数据产品的质量。2018年5月，我国自主研发的大气痕量气体差分吸收光谱仪(environmental trace gases monitoring instrument,EMI)载荷搭载于高分五号卫星(GF-5)发射升空。为监测EMI运行状况，提高载荷数据有效性、可靠性，需要不间断地监测载荷在轨衰变并及时进行校正。

搭载于欧空局RSR-2卫星的全球臭氧监测仪(global ozone monitoring experiment，GOME)[1]，利用每日太阳光谱测量监测光学仪器性能在轨衰变；中分辨率成像光谱仪(moderate-resolution imaging spectroradiometer，MODIS)安装了比辐射计监测太阳反射板BRDF(bi-directional reflectance distribution function)的在轨衰变，通过对太阳漫反射板定标衰变系数矩阵进行查表来完成衰变校正[2]。臭氧监测仪(ozone monitoring instrument，OMI)趋势监测和定标系统TMCF通过分析仪器各项衰变参数，完成在轨性能衰变校正[3]。针对EMI载荷在轨运行参数，Zhao等对EMI发射前定标和第一年在轨定标进行了分析[4]。

本文利用EMI在轨运行期间获取的对地光谱、经漫反射板测量的太阳光谱、内置白光源监测数据等0级数据，结合发射前实验室测量的EMI器件性能及定标参数，对载荷中CCD探测器及漫反射板等关键光学器件在轨衰变进行监测，这为EMI载荷能够长期、高质量运行提供保障。

1 实验部分

高分五号卫星(GF-5)是高光谱光学遥感卫星，设计为太阳同步轨道，轨道高度约705 km。其载荷EMI是4通道背照式帧转移面阵CCD探测器成像光谱仪，探测波长范围240～710 nm，光谱分辨率0.3～0.5 nm，载荷在轨对地观测中，以天底推扫方式实现穿轨方向上的114°和飞行方向上的0.5°的瞬时视场，可以实现一天全球覆盖。

EMI采用4台带有二维电荷耦合器件探测器的光栅成像光谱仪，紫外通道选用紫外增强型帧转移CCD47-20探测器，可见通道选用深势阱型帧转移CCD55-30探测器，EMI光学路径如图1所示。

图1 EMI光学路径示意图Fig.1 EMI optical layout

测量太阳光谱时，EMI在星下点出阴影时打开太阳挡板，太阳光通过太阳网格(透射10%)进入仪器，分别经过备用漫反射板(reference solar diffuser，RSD，40 mm×16 mm)或者石英漫反射板(quartz volume diffuser，QVD，40 mm×16 mm×6 mm)反射到望远镜的折叠镜到达光谱仪。在这个位置的折叠镜阻挡了来自主镜的地球辐射。测量对地辐射光谱时，地球散射光经主副反射镜和偏振扰频器反射到入口的狭缝中，再经中继反射镜将入射光反射到分色滤光片上，最后通过紫外透镜和可见光透镜在光谱仪上成像。暗背景噪声是在星下点黑夜时光谱仪测量的地球散射光。载荷在星下点黑夜区进行内部白光源自检，白光源通过透射漫反射板照亮折叠镜，到达光谱仪。

2 结果与讨论

EMI光机头部至卫星数传系统的图像数据采用LVDS传输格式，其在轨衰变分析是以0级数据为基础，挖掘元件性能衰变特征并建模，在0级到1级数据处理中插入衰变校正算法，实现在轨光谱和辐射准确定标。

2.1 载荷温度监测

EMI面阵CCD探测器温控单元通过PID控制内部制冷器来稳定光学平台的温度，降低温度变化对波长定标和CCD暗背景噪声的影响。为了去除载荷内部水汽和挥发物，2018年5月对载荷进行加热去污处理，自2018年6月后，载荷光学平台温度稳定性较好，UV1温度维持在-18.62 ℃，UV2温度维持在-17.97 ℃，VIS1和VIS2温度稳定在10.63 ℃。仪器温度稳定，对载荷测量数据影响可忽略。各通道CCD探测器温度变化如图2所示。

图2 CCD探测器温度变化Fig.2 Temperature variation of CCD detector

2.2 CCD探测器暗背景噪声监测

EMI载荷虽然在外部设计了铝板屏蔽，来自空间环境的质子仍然会对CCD探测器造成损伤。通过对每轨星下点黑暗区进行暗背景测量，可实现CCD暗背景噪声监测。图3所示为EMI在2019年6月3日5682轨4binning模式测量得到的UV1，UV2，VIS1和VIS2通道暗背景图像，积分时间分别为0.5，2，0.5和2 s。各通道暗背景DN(digital number)值分布均匀，极少量像素点受到太空粒子的辐射，成为响应值大于均值的热点。

图3 CCD探测器暗背景噪声Fig.3 Dark background noise of CCD detector

(1)

(2)

式中，m和n分别为CCD成像面的行数和列数(4binning模式下，紫外通道CCD47-20,258行×1 072列，可见通道CCD55-30,144行×1 286列)，pi,j是每个像素的DN值。

图4(a)所示为各通道暗背景均值变化，从第600轨到第10 060轨，暗背景均值基本成线性递增趋势，各通道暗背景均值年增加率约0.25%～1%。图4(b)所示各通道暗背景标准差整体稳定，震荡幅度在1.5%以内。由于EMI在每个轨道的暗背景信号可在0-1级数据处理中扣除，有效消除暗背景噪声的增大对2级反演的负面影响。

图4 暗背景均值及标准差变化趋势(a)：暗背景均值变化趋势；(b)：暗背景标准差变化趋势Fig. 4 Trend of mean and std. Deviation of dark background(a)：Trend of mean of dark background；(b):Trend of Std.Deviation of dark background

2.3 在轨狭缝函数监测

仪器响应函数(instrument spectral response function，ISRF)是波长定标和大气痕量气体反演中一个关键量[5]。使用EMI在轨测量的0级太阳谱定标数据，结合地面实验室测试得到的ISRF参数[6-7]，完成ISRF在轨参数的反演运算和实时更新。图5所示为EMI测量的中间行0级太阳光谱定标数据。

图5 EMI中间行太阳光谱定标数据Fig.5 Solar spectrum in center field of view of EMI

由发射前实验室分析，其狭缝函数可用二阶高斯函数模型表示，如式(3)所示。

(3)

式(3)中，λ是像素波长,λ0为峰值对应中心波长,σ为仪器狭缝函数半高宽，A为仪器狭缝函数拟合系数。将高分辨的太阳参考光谱SAO2010[8]Srefer(λ)和EMI狭缝函数ISRFEMI(λi)在设定监测光谱窗口内卷积，如式(4)所示

Sfit(λi)=Srefer(λ)⊗ISRFEMI(λi)

(4)

式(4)中，λi指光谱仪第i个像素点波长。Sfit(λi)是在EMI各像素点波长抽样值。再计算实际测量太阳谱Smeas(λi)和理论拟合参考太阳谱Sfit(λi)的残差χ，如式(5)所示。

(5)

式(5)中，N表示CCD每行测量太阳光谱的像素点数。在χ最低时，得到最佳拟合参数ISRFEMI。EMI在轨ISRF参数变化如表1所示。

表1 EMI在轨ISRF参数Table 1 On-orbit ISRF parameters

EMI在轨ISRF参数变化相对平稳，变化幅度约为2.3%。定期计算在轨ISRF数值，在反演痕量气体数据时使用实时更新的ISRF参数可以保障二级数据产品的质量。

2.4 白光源在轨衰变监测

EMI内部白光源采用石英卤钨素灯(QTH，6 V，10 W)监测CCD像素性能和辐射通量。QTH监测光路中，白光经漫透射板(TD)和折叠镜(FM)到达主光路。由于QTH是高度稳定的，漫透射板设计了足够的保护免受污染和宇宙辐射，QTH信号变化也可以部分反映太阳光路衰变。UV2和VIS1 对QTH的响应如图6所示。

图6为UV2和VIS1通道EMI中心视场测量得到的QTH响应。2018/05/28，2018/12/24，2019/05/21的QTH光谱辐射通量数据表明，QTH监测光路在轨响应两年衰变小于1%。由于主镜(PM)的材料与FM相同，由此可以推断在地球太阳光谱测量中辐射通量衰变较小。

2.5 漫反射板在轨衰变监测及校正因子

EMI使用石英漫反射板测量得到的太阳光谱，完成在轨光谱和辐射定标。图7是EMI通过QVD于2018/12/16测量的3221轨0级全视场太阳光谱图，水平和垂直方向分别是光谱维和空间维。

图7 EMI太阳光谱Fig.7 Raw irradiance of EMI

为了滤除太阳光路其他元件的变化，基于备用铝漫反射板(RSD)保护良好，衰变很小，可以采用石英漫反射板(QVD)与RSD的相对值估算石英漫反射板衰变情况，再对太阳光谱辐射定标系数进行校正[9]。

QVD与RSD太阳谱的相对响应计算如式(6)所示

(6)

式(6)中，ΔB(t)为t时刻太阳光谱相对响应，B(t=0)从地面测试中得到，DNQVD(t)和DNRSD(t)是分别经QVD和RSD得到的太阳光谱DN值，Isun是测量得到的太阳光谱，τ为太阳光路光学系统透过率，ηi,j是探测像元i行j列的响应系数。

由于尽量减少了RSD暴露时间，最小化了潜在的衰变，因此可假设BRDFRSD(t)不随时间t变化。如果BRDFRSD(t)在EMI任务寿命期间是恒定的，则其变化ΔB(t)仅由QVD衰变引起。对衰变曲线进行三阶多项式拟合，去除残存的太阳谱夫琅禾费线结构，得到衰变校正因子α(t)，如式(7)所示

α(t)=ΔB(t)BRDFRSD(t)

(7)

式(7)中，ΔB(t)是随着时间变化的相对衰变。α(t=0)是在发射前实验室测定的预设值，衰变校正因子α(t)应用于校正0到1级辐照度辐射定标系数。

图8中黑线ΔB(t=two years)所示为QVD在轨两年来相对衰变，利用三阶多项式拟合，得到衰变校正因子α(t=two years)如红线所示。结果表明，UV2通道的QVD年衰变率小于1.75%，VIS1通道的年衰变率小于1%，VIS2通道的年衰变率小于0.2%。将计算的QVD衰变校正因子作为L1b数据处理器的输入，可对太阳光谱辐射定标系数进行校正。

图8 漫反射板相对衰变及校正示意图Fig.8 Diagram of relative degradation and correction factors for diffuse reflector

3 结 论

以EMI在轨测量的0级数据为基础，对入轨两年来的EMI仪器各项指标进行衰变监测。结果表明，仪器温度基本保持稳定，暗背景噪声均值年增幅度在1%以内，在轨载荷ISRF函数值平稳，变化幅度在2.3%以内，实时更新ISRF函数参数，可有效降低其对二级反演结果的影响，主光路的前两年响应变化幅度小于1%，太阳光路的年衰变率分别小于1.75% (UV2)，1% (VIS1)，0.2% (VIS2)，引入衰变校正因子，可有效校正太阳光谱辐射定标系数，为Llb数据处理提供了关键参数，也为监测EMI在轨性能提供了必要的信息。