AERONET简介、数据下载与绘图

AERONET

全球自动观测网AERONET (AErosol RObotic NETwork)是由NASA 和 LOA-PHOTONS (CNRS) 联合建立的地基气溶胶遥感观测网。该网络现已经覆盖了全球主要的区域，目前全球共有500多个站点，利用CIMEL 自动太阳光度计（SPAM）作为基本观测仪器。AERONET对于研究全球气溶胶的传输、气溶胶辐射效应、验证辐射传输模式以及校验卫星遥感气溶胶的结果，具有重要作用。

AERONET与卫星遥感相辅相成，形成一个“天地一体化”系统，对于包括气溶胶在内的大气污染物的监测具有重要意义。通常，将AERONET所测得的气溶胶光学厚度作为真值，来检验遥感反演的气溶胶光学厚度的准确性。

AERONET为气溶胶研究鉴定、卫星反演的验证、与其他数据库的协同，提供了长期连的、易获取的气溶胶光学、微观物理学以及辐射特性。

系统组成

AERONET由三个部分组成，分别是：CIMEL CE318多波段太阳光度计、Sutron 发射器和电源。
CIMEL Electronique CE318多波段太阳光度计，可测量光谱太阳辐照度和天空辐射度。 最新的CE318-T型号还可以对月球光谱辐射进行夜间测量。

部件：配备25厘米准直仪的传感器头安装在40厘米机器人底座上，根据预先设定的程序系统地将传感器头指向太阳，天空和月亮。 Cimel控制单元，电池和Sutron卫星传输设备通常部署在30厘米×62厘米×46厘米尺寸的防风雨塑料外壳中。 Cimel太阳光度计本身的总重量约为15公斤，控制箱，电池，Sutron卫星发射器，防风雨外壳和其他材料提供了额外的22公斤重量。

卫星数据传输需要Sutron SATLINK2数据记录器/发射器。

CE318太阳光度计

CE318型太阳光度计，是法国CIMEL公司制造的一种自动跟踪扫描的太阳光度计， 分为极化式和普通式太阳光度计两种。该仪器在可见和近红外波段均有 8 个光谱通道（极化式通道有：1020、870、675、440、936、870P1、870P2、870P3，P1、P2和P3为三个极化通道；普通式的通道有：1020、870、675、500、440、380、340、936），仪器的视场角1.2度，太阳跟踪精度0.1度。340nm通道带宽2nm，380nm通道带宽4nm，其余通道带宽10nm。AERONET不仅能自动跟踪太阳作太阳直射辐射测量，而且可以进行太阳等高度角天空扫 描、太阳主平面扫描和极化通道天空扫描。CE318 能自动存储测量数据，并在测 量完成后传输到计算机保存，它还可以通过卫星 DCP 平台远程传输数据，实现 无人管理自动测量、采集和远程数据传输。 CE318 测得的直射太阳辐射数据和 天空扫描数据，主要用来计算大气透过率，反演气溶胶光学和其它特性，如粒度 谱、相函数等。CE318 太阳光度计不仅是一种大气气溶胶环境监测仪器，也可在 遥感卫星传感器辐射定标时进行大气光学参数的测量。

CE318仪器组成

CE318 型太阳光度计主要由以下部分组成:(1)仪器主体:传感器头部、扫描步进马达和机械臂;(2)控制箱:提供软件控制预定的扫描和采样指令，获取数据，内置电池。

工作中的CE318

仪器工作原理

大气光学厚度

地面测得的直射太阳辐射 $E(W/m^2)$在特定波长上根据 Bouguer $(布格)^{注释1}$定律，有公式1:

$$E = E_0R^{-2}\cdot \exp(-m\tau)T_g$$
其中 $E_0$是在一个天文单位(AU)距离上的大气外界的太阳 $辐照度^{注释2}$， $R$ 是测量时刻的日地距离(AU)，$m$ 是大气质量数 $^{注释3}$， $\tau$ 为大气总的垂直光学厚度， $T_g$ 为吸收气体透过率。若仪器输出电压 $V$ 与 $E$ 成正比，则式1可写成公式2:
$$V = V_0R^{-2}\cdot \exp(-m\tau)T_g$$
其中 $V_0$ 是定标常数，在大气相对稳定条件下，进行不同太阳天顶角情况下 的太阳直射辐射测量，仪器输出电压 $V$ 是 $m$的函数， $V_0$ 是从一系列观测值外插 到 $m$为 0 时的电压值 $V$。由$\ln{V}+\ln{R^2}$与 $m$画直线，直线的斜率就是垂直光学厚度 $-\tau$，截距就是太阳光度计在大气外界测得的电压信号 $V_0$，这就是常说的 Langley 法。

Langley原理示意图

大气总的消光光学厚度$\tau$由分子散射(Rayleigh)、气体吸收消光(如臭氧， 水汽)和气溶胶散射三部分组成式3，

$$\tau =\tau_r +\tau_a +\tau_g$$
其中 Rayleigh 光学厚度 $\tau_r$ 由地面气压测值计算出来，在可见近红外波段气体吸收主要是臭氧和水汽的吸收。在没有气体吸收的通道，上式右边的第三项可以忽略，那么从总的光学厚度减去 Rayleigh 光学厚度，气溶胶的光学厚度就 计算出来。

气溶胶参数

对于气溶胶光学厚度，假定气溶胶粒子谱分布遵循容格(Junge)分布，垂 直大气柱气溶胶粒子尺度谱分布如式4：

$$n(r)=\frac{dN(r)}{dr}=c(z)\cdot r^{-(v+1)}$$
其中 $r$是球形粒子的半径， $N(r)为$单位面积上气溶胶粒子总数， $ν$ 是 Junge 参数，因子$C(z)$与高度 $z$有关，正比于气溶胶的浓度。在 Junge 气溶胶谱类型和 气溶胶复折射指数与波长无关条件下，气溶胶光学厚度与波长的关系满足公式5:
$$\tau(\lambda)=k\cdot \lambda^{-v+2}$$
式中 $k$ 为 Angstrom 大气浑浊度系数，是波长 $1μm$ 处大气气溶胶光学 厚度。由式5可知，我们可以通过测量气溶胶光学厚度的谱分布就能求出 $ν 和 k，利用 k 和 ν$继而可以求出其它波长上的气溶胶光学厚度。

改进 Langley 法

在地面测得的直射太阳辐射信号在 $940nm$ 附近水汽吸收带不符合 Bouguer 定律，Bouguer 指数消光定律是对单色辐射而言。依照 Bruegge 和 Halthore 在 1992 年的研究成果，水汽透过率这时用两个参数表达式来模拟，式6:

$$T_w=exp(-aw^b)$$
其中 $T_w$是通道上的水汽吸收透过率， $w$ 是大气路径水汽总量，$a 和 b$是常数，在给定的大气条件下，它们与太阳光度计 $940nm$ 通道滤光片的波长位置、宽度和形状有关，还与大气中的温压递减率和水汽的垂直分布有关。 $a 和 b$ 由辐 射传输方程模拟来确定。为了在各种大气条件下能有效利用太阳光度计反演水汽 量，有必要研究 $a 和 b$对这些条件的灵敏度。
在 940 $nm$ 水汽吸收带，太阳光度计对太阳直射辐照度的响应可表示为，式7:
$$V=V_0R^{-2}\cdot \exp(-m\tau)\cdot T_w$$
其中 $\tau$ 是 Rayleigh 散射和气溶胶散射光学厚度，它们相互独立，气溶胶光学 厚度通过其它通道(如 870 $nm$和 1020 $nm$)内插得到。斜程水汽量 $w=m*PW$， $PW$ 为垂直水汽柱总量。将(6)式代入并两边取对数，得
$$\ln V +m\tau = \ln{V_0R^{-2}}-am^b\cdot PW^b$$
在稳定和无云大气条件下，以 $mb 值为 X$轴，以上式左边为 $Y$轴画直线， 直线的斜率为 $-a PW^b$， $Y 截距为 \ln (V_0R^{-2})$。这就是通常所说的改进的 Langley 法。

数据产品

AERONET数据通过网络工具和FTP通过互联网分发，使用ASCII文本格式进行数据传播，这些数据压缩在zip文件内。

AERONET可提供全球分布的气溶胶光学厚度（AOD）反演产品。版本3 的AOD数据是具有三个数据质量水平：Level 1.0（未筛选），Level 1.5（云筛选和质量控制）和Level 2.0（保证质量）。处理算法已经从版本1.0发展到版本2.0，现在发展到版本3.0。

数据下载（以AOD数据下载为例）

打开官网：

准备下载版本3的气溶胶光学厚度数据，在AERONET DATA ACCESS中选择版本3中的Data Display。
根据需要设置筛选条件

数据下载成功后将文件后缀修改为csv后使用excel打开

绘图与分析

在同一图上绘制每个波长的总AOD（全部为JD）。 可以注意到较小的波长会看到较高的光学深度。这是因为气溶胶光学深度通常随着波长的增加而减小。

从上面两图可以看到这个月有什么趋势吗？图中反映的是AOD由于Rayleigh散射对于相同波长与JD的关系。

使用Lev20文件，绘制JD和Solar Zenith Angle。太阳光度计运行的时间大部分是在中午（90度）左右开始，一直持续到当天晚些时候，这一点通过仔细观察JD得到证实。

现在我们将使用500 nm波长的ONeill_15文件，制作JD与总体，细和粗粒气溶胶外径，JD与总气溶胶OD的标准偏差细粒与粗粒气溶胶外径的标准。
大多数时间这两种类型的数量之间没有太大的差别，但有6种情况，粗气溶胶外径占总AOD的比例，这些可能来自内陆地区盐雾气候的强烈海风，也可能在中部 这段时间，细微的气溶胶OD值有轻微的增加。

结语

数据下载很容易，重要的是拿到数据后的科学分析，加油。

注释：
[1]Bouguer（布格）定律，又称比尔-朗伯定律、朗伯-比尔定律、布格-朗伯-比尔定律（Bouguer–Lambert–Beer law），是光吸收的基本定律，适用于所有的电磁辐射和所有的吸光物质，包括气体、固体、液体、分子、原子和离子。
概述：
一束单色光照射于一吸收介质表面，在通过一定厚度的介质后，由于介质吸收了一部分光能，透射光的强度就要减弱。
设入射光强度为$I_0$，吸收光强度为$I_a$，投射光强度为$I_t$，反射光强度为$I_r$，则：

$$I_0=I_a+I_t+I_r$$
投射光的强度 $I_t$与入射光强度 $I_0$之比称为透射率，用 $T$表示：
$$T=\frac{I_t}{I_0}$$
透光率越大，溶液对光的吸收越少；反之，透光率越小，溶液对光的吸收越多。
透光率的负对数称为吸光度，用A表示。A越大，溶液对光的吸收越多。
$$A=-\log T = \log \frac{I_0}{I_t}$$

Lambert定律

$$A=k_1 b$$
其中 $b$溶液厚度，$k_1$比例系数。
Beer定律
$$A=k_2 c$$
其中 $c$是物质的量浓度，$k_2$与被测物质的性质、入射光的波长、溶剂、液层厚度及温度有关的常数。

Lambert-Beer定律
当一适当波长的单色光通过溶液时，若溶液厚度一定，则吸光度与溶液浓度成正比，且与吸光物质种类、溶剂、入射光波长、液层厚度和溶液温度有关，且对所有的均匀介质和单色光都适用。

吸收介质的浓度$(c)$愈大，介质的厚度$(l)$愈大，则光强度的减弱愈显著，其关系为：

$$A=K\cdot l\cdot c$$
其中， $A$ 为吸光度，$K$是系数，可以是吸收系数或摩尔吸收系数。
当介质厚度 以 $cm$ 为单位，吸光物质浓度 以 $g/L$ 为单位时,称为吸收系数，其单位为 。当吸光物质浓度 以 $mol/L$ 为单位时,称为摩尔吸收系数。

[2]辐照度（入射度）：当电磁波由体外穿入体内时所构成的辐射通量密度，通称为入射度(Irradiance)，单位$W/m^2$ 。
辐射通量：单位时间内穿过某一面积的电磁辐射能量 单位（$w$）。
通量密度：单位时间内，穿过单位面积的电磁辐射能量($W/m^2$)。
辐射强度：描述点辐射源强度的物理量，描述了辐射通量在2π空间中的变化状况。
辐射亮度：表述面光源强弱的物理量，单位投影面积上单位立体角内的辐射通量。
…..遥感物理…..

[3]大气质量：大气对地球表面接收太阳光的影响程度被定义为大气质量（air mass）。

参考
https://aeronet.gsfc.nasa.gov/