地表净辐射是驱动大气运动的主要能量来源，在地气系统的相互作用中起着关键作用.其中，地表短波净辐射(net surface shortwave radiation，NSSR)是地表净辐射的一个重要组成部分，地表短波净辐射的研究，对地表辐射平衡、蒸散以及气候变化等众多领域具有重要研究意义.

传统获取地表短波净辐射通量的信息只能通过仪器实地测量而得.虽然这种需要利用地面台站气象数据方法计算的精度较高，但是计算麻烦，需要输入很多难以准确获取的气象参数，并且需要投入大量的人力物力财力，很难完成区域尺度的地表短波净辐射数据的获取.随着遥感技术的发展，基于航空和卫星数据来获取地表短波净辐射的方案被提出.自1984年起，许多学者针对大气顶部短波净辐射和地表短波净辐射之间的线性关系进行了研究^[1-4].Li等^[5-6]发现对于一个特定的太阳天顶角、辐射通量和地表反照率之间的关系是线性的，给定一个大气顶部的反射通量和太阳天顶角，仅存在唯一的地表短波净辐射通量值，基于此提出一个参数化模型.Masuda等^[7]考虑到云顶高度与类型、地表压强、臭氧数量和气溶胶类型与含量的影响，从而改进了Li等的参数化模型.随后，很多研究者们基于Li等和 Masuda等的参数化模型，计算了不同地区地表短波辐射通量^[8-12].由于他们的模型输入参数中，大多采用的是宽通道卫星观测数据，2008年Tang 等^[11-12]利用极轨卫星MODIS(moderate-resolution imaging spectroradiometer)基于Li等的方法，通过中等波谱分辨率的大气辐射传输模型MODTRAN模拟，建立窄通道观测数据与大气层顶表观反照率之间的关系，提出利用MODIS 7个窄通道大气顶部观测值估算地表短波净辐射的新方法，解决了传感器大气顶部窄通道表观反射率向宽通道反照率转换的问题.

虽然当前已有众多学者利用遥感数据进行了地表短波净辐射的反演研究，但是这些方法大多是基于极轨卫星数据.对于一个给定的研究区来说，当前极轨卫星在白天最多只能提供1至2个时相的数据，如遇云雨等不利天气条件，可用的数据将难以得到保障，从而严重地制约了对地表短波净辐射的研究.此外，虽然当前全球能量与水循环试验地表辐射收支产品数据集(GEWEX-SRB)能够提供时间分辨率为3h的辐射产品，但其空间分辨率较低(1°)，导致其在局地尺度的研究上难以应用.如何有效兼顾地表短波净辐射数据的时间和空间分辨率，发展更为合适的卫星地表短波净辐射遥感估算方法，使其更有效地在区域气象、水文和生态环境等领域发挥其应有的作用，具有十分重要的科学意义和实用价值.考虑到静止气象卫星对给定的像元具有固定的观测角度，一天中能够提供 48~96个时相的多光谱数据，并且其相对较高的空间分辨率(3~5km)通常也能够满足各研究和应用领域的需求.因此，本文针对风云静止气象卫星(FY-2E)数据，发展了地表短波净辐射估算方法，并利用实测数据对反演结果进行验证.

1 研究区与数据 1.1 研究区概况

黄河源区位于青藏高原东部，是黄河上游的重要水源补给区.该地区地下特征显著且复杂，海拔落差较大，其地表辐射通量在量级与分布格局上存在显著的差异.分析该区域的净辐射变化情况，如此独特环境的辐射平衡，对了解黄河源区地气间相互作用过程、区域能量收支等方面具有重要意义.但由于地面观测资料的限制，对该地区不同时空尺度净辐射的研究还不多.

本文以黄河源区中的玛曲县为研究区(图 1)，该区地理坐标为33°30′~34°15′N，101°38′~102°45′E，位于甘肃省甘南藏族自治州西南部，青藏高原东端，甘、青、川三省交界处，黄河第一弯曲部.地势西高东低，由西北向东南倾斜，海拔在3300~4806m，属于高寒湿润型气候，年平均气温1.1℃，年平均降水量615.5mm^[13].

1.2 FY-2E数据

风云二号气象卫星(FY-2)是中国自行研制的第一代地球静止轨道气象卫星，E星(FY-2E)于2008年12月23日发射成功，并定点于东经123.5°赤道上空，接替FY-2C，与极地轨道气象卫星相辅相成，构成中国气象卫星应用体系.FY-2E可实现每1h获取经纬度从N60°~S60°，E45°~165°的一幅完整影像；汛期每间隔30min可获取一幅影像.FY-2E卫星主要有效载荷为红外和可见光自旋扫描辐射器VISSR^[14]，其技术指标如表 1所示.采用FY-2E卫星，一天中最多能够获取48个时相的数据，并且其有相对较高的空间分辨率(5km)，能够满足研究的需求.

1.3 实测数据

在研究区内选取若尔盖高原湿地生态系统研究站玛曲观测场 (102.14°N，33.89°E，如图 1中三角形所示)，该地海拔3423m.本文获取了2010年北京时间上午11点该站点的下行短波辐射与上行短波辐射实测数据，其地面观测上行和下行短波辐射由kipp&zonen公司的CNR-1 传感器观测得到，其光谱范围为0.3~2.8μm.

2 方法研究

地表短波净辐射可以由Pinty和Ramond^[15] 及Pinker et al.^[16]提出的下式求得：

${{S}_{n}}=R_{s}^{\downarrow }-R_{s}^{\uparrow }=\left( 1-r \right)R_{s}^{\downarrow },$

(1)

其中，S_n表示地表短波净辐射(W/m²)；R_s^↓和R_s^↑分别表示下行短波辐射(W/m²)和上行短波辐射(W/m²)；r表示地表反照率，无量纲.

在利用FY-2E数据计算地表短波净辐射时，地表反照率是一个关键的地表参数.由于FY-2E在可见光波段为一个较宽的观测波段(0.510~0.905μm)，因此，本文不进行窄波段到宽波段的反演，只将 FY-2E 观测的每小时的表观反射率通过简化黑暗像元法进行大气校正得到地表反射率后，直接作为短波波段地表反照率参与短波净辐射的计算^[17-18].

2.1 太阳总辐射的确定

太阳总辐射R_s^↓可以表示为^[17]

$R_{s}^{\downarrow }=G\times cos(SZA)\times {{d}_{r}}\times \tau ,$

(2)

其中，G为太阳常数(1367W/m²)；SZA表示太阳天顶角；d_r指日地距离因子，无量纲；τ是大气透过率，无量纲.

日地距离因子可以由下式计算得到：

$\begin{align} & {{d}_{r}}=1.00011+0.034221cos\alpha +0.00128sin\alpha + \\ & 0.000719cos2\alpha +0.000077sin2\alpha ', \\ \end{align}$

(3)

其中，α=2π(DOY－1)/365，DOY表示年积日.

2.2 大气透过率的确定

因为大气透过率的变化主要取决于大气水汽含量的动态变化^[19]，本文通过大气水汽含量来计算大气透过率，利用FY-2E的2个相邻热红外通道的地表比辐射率之比反演相应时刻的大气水汽含量^[20-21]，计算过程如下：

$WVC={{C}_{1}}+{{C}_{2}}\times {{\tau }_{IR2}}{{\tau }_{IR1}}$

(4)

$\frac{{{\tau }_{IR2}}}{{{\tau }_{IR1}}}=\frac{{{\varepsilon }_{IR1}}}{{{\varepsilon }_{IR2}}}\times \frac{\sum _{k=1}^{N}({{T}_{IR1,k}}-\overline{{{T}_{IR1}}})({{T}_{IR2,k}}-\overline{{{T}_{IR2}}})}{\sum _{k=1}^{N}{{({{T}_{IR1,k}}-\overline{{{T}_{IR1}}})}^{2}}},$

(5)

其中，WVC表示大气水汽含量(g/cm²)； τ_IR1与τ_IR2为通道IR1和IR2的大气透过率，无量纲；ε_IR1与ε_IR2表示通道IR1和IR2的比辐射率，无量纲；T_IR1与T_IR2分别表示通道IR1、IR2的N个相邻像元的大气顶层亮温的平均值(K)；C₁与C₂ 可以通过如下公式计算得到：

${{C}_{1}}=28.104-\frac{14.996}{cos(VZA)}+\frac{3.211}{co{{s}^{2}}(VZA)}$

(6)

${{C}_{2}}=-28.056+\frac{14.954}{cos(VZA)}-\frac{3.206}{co{{s}^{2}}(VZA)},$

(7)

其中，VZA表示观测天顶角.

利用大气透过率的变化主要取决于大气水汽含量的动态变化这一特征，本文基于MODTRAN模拟大气水汽含量与大气透过率之间的关系，通过得到的大气水汽含量计算大气透过率.又由于本研究拟反演中国地区的地表短波净辐射，故根据MODTRAN中标准大气模型的设置，挑选出适合中国地区的大气模型——中纬度地区夏季和冬季大气模型，选用Rural-VIS=23km的气溶胶模式.分别模拟中纬度地区夏季大气模型和中纬度地区冬季大气模型中FY-2E可见光波段大气透过率(τ)与大气水汽含量(w)的关系(表 2).

3 结果分析

基于上述方法，随机选取2010年的一个晴天(2010年9月15日)，反演当天北京时间上午11点研究区的地表短波净辐射，如图 2所示.可以看出地表短波净辐射由南至北、由东至西递减，与太阳天顶角有关.为进一步分析地物类型对地表短波净辐射的影响，同时做出下午2点与上午8点的地表短波净辐射差值图(图 3)与土地利用图(图 4)，从图中能明显看出，越高植被覆盖地区地表短波净辐射差值越低，这与实际情况完全吻合，植被覆盖率越高，地表吸收的辐射能越低.

3.1 反演结果与实测数据对照分析

为了检验方法的可靠性，选取2010年玛曲县的所有晴天(共33d)的地面实测数据对反演结果进行对比验证.图 5给出2010年全年所有晴天地表短波净辐射地面测量值与FY-2E数据反演结果的散点图.经过计算得出：晴天条件下，均方根误差(RMSE)为 27.51W/m²，决定系数(R²)为 0.99.

本研究在不考虑冬天(1、11和12月)的晴天情况下，对剩下晴天(共22d)的反演结果与地面实测数据进行验证分析，如图 6所示，R²为0.99，RMSE为15.36W/m².

3.2 误差分析

FY-2E卫星的可见光波段范围为0.510~0.905μm，而地面实测上行和下行短波辐射数据的波段范围为0.3~2.8μm，因而反演结果对比实测数据具有一定的误差.另外，研究区海拔在3300~4806m，高程落差较大，地形崎岖不平，地表粗糙度较大，同时，较大的景观多样性指数与空间缀块性进一步造成了该地区的空间异质性较为强烈，对地表比辐射率的反演的精度影响较大，从而导致大气水汽含量的误差，进一步降低地表短波净辐射的准确性.其次，风云卫星数据本身的信噪比较低，从而降低了反演精度.虽然观测塔观测的辐射数据能够代表一定范围的辐射状况，但是用观测值与5km×5km像元值对比验证，也会导致一定的误差.另外，除去冬天后的反演结果精度更高的原因可能与冬天冰雪覆盖有关.

Huang et al.^[22]的研究中对在玛曲县GLASS的下行短波辐射产品与实测数据进行验证，其中RMSE为130.4W/m²，相关系数(R)为0.84.若需得到地表短波净辐射，还需获得上行短波辐射或地表反照率，因此，最终误差必定更大.因此，可以证明此方法比GLASS产品的方法更精确可行.另外，唐伯惠^[10]的研究中提到，基于MODIS数据反演的地表短波净辐射结果与禹城(下垫面比较均一的地区)站点2003年1、4、7、10月的实测数据进行的对比验证，结果发现RMSE小于20W/m².本文的研究结果在除去冬季的1、11、12月份后，反演结果也达到了唐伯惠^[10]文中的精度，且本文选取的研究区下垫面异质性更高，有些地区常年冰雪覆盖，对反演精度影响更大.上述结论，可以证明该方法反演地表短波净辐射精确可靠.

4 结论

本文发展了适合于中国地区的地表短波净辐射的反演方法.首先，基于MODTRAN模拟大气水汽含量与大气透过率之间的关系，然后利用FY-2E的2个相邻热红外通道的地表比辐射率之比反演得到大气水汽含量，进而反演得到地表短波净辐射.

结合地表发射率空间分布图、地表短波净辐射差值图与土地利用图能明显看出反演结果真实可信.此外，基于玛曲地区2010年的地面实测气象数据，本文将对应时刻的地表短波净辐射反演结果进行了验证.结果显示，对于晴天，均方根误差为 27.51W/m²，决定系数为 0.99；在除去冬季的晴天条件下，R²为0.99，RMSE为15.36W/m².因此，该结果证明本文的地表短波净辐射反演方法精确可行，从而为反演同时具备高时间和空间分辨率的地表短波净辐射数据提供坚实依据,对长时间尺度气象变化和陆表水文循环等方面的研究具有一定意义.

中国科学院寒区旱区环境与工程研究所西部气候环境与灾害实验室提供了若尔盖高原湿地生态系统研究站玛曲观测场的辐射数据，国家气象卫星中心(http://satellite.nsmc.org.cn/portalsite/default.aspx)提供了FY-2E数据，在此一并致以衷心的感谢！