紫外辐射(290～400 nm)占太阳总辐射的比例很小，在大气层顶也只有8%左右(Lovengreen et al.,2000；Santos et al.,2011)，但由于其光量子能级较高，所产生的光化学作用、植物光合作用等生物效应十分显著。紫外辐射对气候、生态环境及人类健康都有重要影响(邓雪娇等，2003)，特别是随着臭氧洞的发现，紫外辐射的观测估测已逐渐成为当前研究的焦点之一。

紫外辐射可以通过破坏蛋白质中的化学键来损坏生物个体细胞，通过破坏植物叶片来抑制光合作用，从而对整个生态系统产生影响。另外，平流层臭氧的减少会导致对流层紫外辐射上升，这将加快近地面层光化学反应的发生，产生更多的二次污染物，使得近地面空气质量恶化，对人类健康和生活环境产生直接影响(Diffey，1991)。澳大利亚每10万人中有800人、美国每10万人中有250人得皮肤癌，在盛行日光浴的美国和加拿大，每年新增的皮肤癌患者超过100万人。因此有必要大量开展紫外辐射的观测与估测研究。

目前紫外辐射的预报方法可分为模式和统计预报两种方法(吴兑，2000)。沈元芳等(2009)应用GRAPES(Global/Regional Assimilation and PrEdiction System)模式中的Goddard短波辐射方案[①]，建立了紫外线数值预报系统。傅炳珊等(2002)在大气辐射传输理论的基础上，利用中分辨率大气辐射传输模式，建立晴空或少云天气状况下石家庄市紫外波段辐射强度及指数预报模式。李春(2003)依据拉萨地区不同月份紫外辐射占全波段太阳总辐射的比例关系和全波段太阳总辐射的气候学计算方法，提出了一套紫外线指数统计估测方法。王繁强等(2005)、毛则剑和张立峰(2007)、丛菁等(2009)、曲晓黎等(2011)分析紫外辐射与能见度、云量、温度、相对湿度等气象要素的相关性，选出相关性较大的要素作为因子，采用多元回归方法估测紫外辐射，但是直接采用气温日较差来进行紫外辐射估测的研究较少，本文拟利用2005～2011年长时间的辐射观测资料分析总辐射与紫外辐射的相互关系，在此基础上通过气温日较差建立总辐射估测方法，从而建立紫外辐射估测模型。

禹城站(36°40′N，116°22′E)位于我国华北平原东部，该地区是典型暖温带半湿润季风气候 区，年平均气温13.1°C，降雨量582 mm，水面蒸发力952 mm，太阳辐射总量5225 MJ m⁻²，日照时数2640小时，大于0°C积温为4951°C，大于10°C积温为4441°C，无霜期200天，光热资源丰富，雨热同期，有利于农业生产。

气象要素使用Milos520(Vaisala，芬兰)自动气象站观测，其中温度精度为0.2°C，相对湿度精度为±2%。自动站辐射观测使用仪器为Kipp & Zonen(荷兰)生产的辐射表，辐射表的参数如表 1所示。气象要素和辐射数据采用DM520数据采集器采集，采集频率为每分钟采集一次。

表1(Table 1)

表 1 辐射表参数 Table 1 The parameters for pyranometers 辐射类型 表型 准确度 总辐射(0.3～3 μm) CM11 ±3% 紫外辐射(0.29～0.4 μm) CUV3 ±10%

表 1 辐射表参数 Table 1 The parameters for pyranometers

总辐射表每年在春季采用“交替法”进行标定，标定精度为±3%，紫外辐射表采用辐射标准灯和紫外—可见光光谱仪进行标定，标定精度为±10%。

(1)同一地理位置观测的总辐射和紫外辐射应分别小于大气层顶的总辐射与紫外辐射，否则直接将数据剔除(Geiger et al.,2002)。大气层顶的总辐射和紫外辐射量采取Foyo-Moreno et al.(1999)提出的计算方案。

(2)夜间辐射的观测值应该小于0，或者由于仪器零点漂移造成不大于5 W m⁻²瞬时输出值，处理时全部用0代替。

(3)白天总辐射观测的极小值应该大于连续阴天时总辐射的观测值，即剔除白天平均无云指数(地基总辐射与大气层顶总辐射的比值)小于0.03的数据(Omran，2000)。

(4)由于辐射仪器的余弦效应，在太阳高度角较低时会带来很大的观测误差，剔除太阳高度角小于10°的数据(Huang et al.,2011)。

对禹城近7年的瞬时辐射数据进行严格的数据质量控制，直接删除其中不合格的数据，删除比例为1.1%左右。

通过2005年1月1日至2011年12月31日的辐射观测数据对禹城到达地面的紫外辐射变化规律，紫外辐射与总辐射及其相关要素的关系进行分析。

从图 1可以看出，总辐射和紫外辐射日累计值的季节变化规律一致，都是冬季小，春末夏初大，总辐射日累计值的变化范围在5～30 MJ m⁻² d⁻¹之间，年均值13.525 MJ m⁻² d⁻¹，紫外辐射日累计值的变化范围在0.1～1.20 MJ m⁻² d⁻¹之间，年均值为0.468 MJ m⁻² d⁻¹。造成这种变化的原因主要是地球公转导致一年中的太阳位置季节变化。夏季太阳直射点靠近北回归线，此时华北地区的太阳高度角较大，禹城地区接收的天文辐射较大，即该区域的总辐射和紫外辐射都达到极大值；冬季则反之。Hu et al.(2010)利用北京地区2005年1月到2008年6月的辐射资料，得到紫外辐射平均日累计值为0.37±0.17 MJ m⁻² d⁻¹。张兴华等(2012)利用拉萨2005～2010年的辐射资料，得到紫外辐射平均日累计值为0.87 MJ m⁻² d⁻¹，比禹城高，主要是因为拉萨海拔高，空气比较稀薄，对太阳辐射的削弱小，到达地面的太阳辐射较强。

图 1

Fig.1

图 1 禹城2005～2011年总辐射和紫外辐射日累积值的逐日变化 Fig.1 Temporal variation of the daily total values of solar radiation and ultraviolet(UV)radiation in Yucheng during 2005−2011

2005～2011年禹城总辐射和紫外辐射的月平均日累计值变化规律如图 2所示，总辐射和紫外辐射呈现“夏季大冬季小”的季节变化特征。总辐射和紫外辐射月平均日累计值最小值出现在2006年1月，其数值分别为4.864 MJ m⁻² d⁻¹和0.139 MJ m⁻² d⁻¹；最大值出现在2009年6月，分别为22.693 MJ m⁻² d⁻¹和0.834 MJ m⁻² d⁻¹。表 2给出了禹城2005～2011年月降水量的季节变化，降水主要集中在7、8月份，云量和大气中的水汽增加，更多的太阳辐射被吸收，总辐射与紫外辐射月平均日累计最大值都出现在5、6月份，而不是相对多雨的7、8月份。

表2(Table 2)

表 2 禹城2005～2011年平均月降水量的季节变化 Table 2 Mean of the monthly precipitation in Yucheng during 2005−2011 月份 1 2 3 4 5 6 7 8 9 10 11 12 降水量/mm 1.0 8.0 14.4 27.9 50.3 57.9 132.1 116.1 45.0 13.6 11.4 2.7

表 2 禹城2005～2011年平均月降水量的季节变化 Table 2 Mean of the monthly precipitation in Yucheng during 2005−2011

图 2

Fig.2

图 2 2005～2011年紫外辐射的逐月变化Fig. 2 The monthly variation of UV radiation during 2005−2011

总辐射与紫外辐射的日变化规律一致(图 3)，都呈现为早晚小中午大，即06:00(地方时，下同)总辐射和紫外辐射值最低，分别为0.298 MJ m⁻²和0.01 MJ m⁻²；随着太阳高度角的增大总辐射和紫外辐射逐渐增大，到12:00达到最大，其值分别为1.98 MJ m⁻²和0.07 MJ m⁻²，之后随着太阳高度角的减小，辐射量逐渐减小，到18:00出现极小值；总辐射和紫外辐射的极小值分别为0.299 MJ m⁻²和0.01 MJ m⁻²。导致这种变化的主要原因是地球自转，由于太阳高度角日变化规律造成辐射的日变化特征。张兴华等(2012)利用拉萨地区2005～2010年数据得到总辐射小时累计值变化范围为0.227～2.970 MJ m⁻²，紫外辐射为0.009～0.126 MJ m⁻²，比禹城地区稍大。

图 3

Fig.3

图 3 近7年禹城地区总辐射和紫外辐射的日变化图Fig. 3 Diurnal variation of the 7-year mean hourly values of solar radiation and UV radiation in Yucheng

由于紫外辐射与总辐射有相同的季节变化和日变化规律，因此通过研究紫外辐射与总辐射比值的变化规律有助于我们更好地了解紫外辐射，利用该比值的变化特征结合总辐射的历史观测数据获得大量的紫外辐射用于相关的研究。图 4为2005～2011年紫外辐射与总辐射比值的月变化规律图，其变化范围为0.023～0.046，这与Robaa(2004)利用4年的埃及开罗数据得到的比值范围0.027～0.049一致。江灏等(1994)利用HEIFE绿洲区1991年的观测数据得到比值范围为0.026～0.034，比禹城地区略小。HEIFE绿洲区1991年平均相对湿度为49.64%，禹城2005～2011年年平均相对湿度为66.94%，可能由于禹城的水汽含量相对较高，导致比值较大。紫外辐射与总辐射比值的季节变化规律也呈现“夏季大冬季小”的特征，且夏季变化幅度大于冬季。紫外辐射与总辐射的比值最大值在2008年7月，为0.046；最小值在2005年1月，为0.023。紫外辐射与总辐射比值的季节变化主要取决于水汽的季节变化(Hu et al.,2010)，用相对湿度来表 征水汽含量，由图 4可以看出，比值F_uv与相对湿度RH的变化一致，水汽对近红外辐射的吸收较强，而对短波辐射的吸收相对较少，水汽含量的增多使得比值变大，反之则变小。

图 4

Fig.4

图 4 2005～2011年紫外辐射与总辐射比值和相对湿度(RH)的逐月变化Fig. 4 Temporal variation of the monthly mean ratio of ultraviolet radiation to solar radiation and relative humidity(RH)during 2005−2011

影响到达地面的紫外辐射量的因子除了太阳高度角、臭氧外，云和气溶胶的影响也较大。然而，由于云和气溶胶的观测较少，定量化研究也比较困难，通常通过平均晴空指数(地面观测到的太阳总辐射与大气层顶的总辐射之比)作为云和气溶胶的识别因子(张兴华等，2012)来开展气溶胶、云对紫外辐射影响的研究。图 5为紫外辐射与总辐射的比值随晴空指数K_s的变化趋势图，在K_s小于0.4时，比值的变化范围为0.02～0.08，波动较大，原因可能是K_s较小时，对应于太阳高度角较小或者云量较多的情况。在太阳高度角较小时，大气光化学路径相对较长，大气分子的瑞利散射对波长较短的紫外辐射有较强的削弱作用，致使比值较小；在云量较多时，可能由于水汽大量吸收波长较长的红外辐射致使比值较大，也可能由于水汽中的细微粒对紫外辐射的削弱较强致使比值较小。在K_s大于0.5时比值相对比较稳定，基本维持在0.03左右，与Feister and Grasnick(1992)的结论一致，可能是因为K_s较大时，总辐射较高，多在夏季晴天，天气条件比较一致，使得比值相对稳定。

图 5

Fig.5

图 5 紫外辐射和总辐射比值与晴空指数(Ks)的关系Fig. 5 Relationship between the ratio of ultraviolet radiation to solar radiation and the clearness index(Ks)

图 6为不同K_s范围时，紫外辐射(R_uv)随太阳高度角正弦值(μ)的变化，K_s从0.03开始，按0.01增长，不同颜色代表不同的K_s值，对于特定的K_s，紫外辐射随μ以指数形式增长(Hu et al.,2010)，即

图 6

Fig.6

图 6 紫外辐射随太阳高度角正弦值的变化(不同颜色代表不同的Ks)Fig. 6 UV irradiances as function of sine of the solar elevation angle(different Ks values were represented by different colors)

式中，${R_{{\rm{uvm}}}}$表示单位K_s区间内紫外辐射最大值，随机选取2005～2011年一半的数据做其与K_s的散点图，如图 7所示，满足三次函数关系，即

图 7

Fig.7

图 7 单位晴空指数区间内紫外辐射最大值(Ruvm)与晴空指数Ks的关系Fig. 7 The relationship between Ruvm(the maximum value of UV radiation per Ks) and Ks at Yucheng station

紫外辐射按小时累计值的计算公式为

紫外辐射日累计值的计算公式为

其中，t_d为太阳高度角大于10°的小时数，$\overline {{K_s}} $为平均晴空指数(到达地面的日累计总辐射与日累计天文辐射之比)。代入参数a，b，c，d，N的拟合值，得到紫外辐射日累计值的计算公式

把随机选取后剩余的一半数据代入公式(5)计算紫外辐射日累计值，图 8是计算值与实测值得对比，平均相对误差为0.09，91%的数据相对误差小于0.2，决定系数R²为0.97，图中RMSE为均方根误差，N_s为样本数量。根据公式(5)计算紫外辐射日累计值的精度较高，只要能够估测出太阳总辐射即可实现紫外辐射的估测。

图 8

Fig.8

图 8 紫外辐射计算值与实测值的对比Fig. 8 The comparison of calculated and measured UV radiation, Ns represents the number of samples

公式(5)中太阳高度角正弦值μ和t_d可以通过纬度、经度和日期等计算得到，只要能够估测出晴空指数K_s即可得到紫外辐射的估测。

目前很多学者利用常规气象观测要素建立太阳辐射的估算模型，其中，最常用的有日照百分率模型和气温日较差模型。Angstrom(1924)最先提出太阳辐射的日照百分率模型

式中，S为日照百分率，a，b为经验参数。

图 9为晴空指数的季节变化，该变化趋势与图 10中气温日较差(Δt)的季节变化趋势一致。根据地表能量平衡，白天到达地表的太阳辐射较多时，云量相对较少，昼夜温差增大(Almorox et al.,2011)，所以晴空指数与气温日较差的变化趋势一致。Bristow and Campbell(1984)提出了一个基于气温日较差的太阳辐射估算模型，该模型认为日总辐射是天文辐射和气温日较差的函数，即

图 9

Fig.9

图 9 晴空指数的季节变化Fig. 9 The seasonal variation of clearness index

图 10

Fig.10

图 10 温度日较差开方()的季节变化Fig. 10 The seasonal variation of the square root of diurnal temperature range()

Wu et al.(2007)、叶建圣(2010)、邓艳君(2012)等人比较了两种模型，发现日照百分率模型的计算值更接近于观测值。但是，次日的日照百分率估测比较困难，因此选用气温日较差模型对太阳总辐射进行估测。

到达地面的太阳总辐射与气温日较差、水汽等相关，很多学者对气温日较差模型进行了改进，Hargreaves et al.(1985)、Hunt et al.(1998)、Wu et al.(2007)等人提出了改进的气温日较差模型，用Hunt^a和Hunt^b表示Hunt et al.(1998)提出的两种计算方法。

Hargreaves et al.(1985)：

Hunt^a：

Hunt^b：

Wu et al.(2007)：

随机选取2005～2011年一半数据来拟合(7)～(11)式，a₀，a₁，a₂，a₃，a₄的值如表 3所示，R_s为到达地面的太阳总辐射(MJ m⁻² d⁻¹)，R₀为到达大气顶的太阳总辐射(MJ m⁻² d⁻¹)，t_max为日最高气温(°C)，t_min为日最低气温(°C)，P为日降水量(mm)，t_m为日平均气温(°C)，T_dm为日平均露点温度(°C)，P_t代表降水量对太阳辐射的影响，当降雨量大于0 mm时，P_t=1，在原公式上加一个常数，否则P_t为0。

表3(Table 3)

表 3 计算总辐射的参数 Table 3 The parameters of total radiation calculation 计算方法 a0 a1 a2 a3 a4 R2 Bristow and Campbell(1984) 0.585 0.030 1.649 — — 0.71 Hargreaves et al.(1985) 0.138 －0.568 — — — 0.68 Hunta 0.170 －0.141 －0.403 0.005 −0.734 0.73 Huntb 0.578 －0.014 0.044 0.009 — 0.71 Wu et al.(2007) －0.014 0.129 0.008 －0.068 －0.009 0.74

表 3 计算总辐射的参数 Table 3 The parameters of total radiation calculation

由于Wu et al.(2007)计算的总辐射与实测总辐射的R²最大(表 3)，且平均相对误差最小，选取该方法来估测总辐射。

把随机选取后剩余的另一半数据代入Wu et al.(2007)的计算公式，得到总辐射计算值，与实测值的比较如图 11所示，R²达0.80，平均相对误差为0.19，拟合效果较好。

图 11

Fig.11

图 11 总辐射计算值与实测值比较Fig. 11 The comparison of calculated and measured solar radiation

由计算出来的太阳辐射日总量与天文辐射日总量之比得到平均晴空指数，代入公式(5)估测紫外辐射日总量。图 12为紫外辐射估测值与实测值的对比，R²达到0.73，平均相对误差为0.22，估测效果较好。

图 12

Fig.12

图 12 紫外辐射计算值与实测值的对比Fig. 12 The comparison of calculated and measured UV radiation

以紫外辐射月平均日总量最大值的十分之一作为紫外线指数的一个标准单位(祝青林等，2005)，并按照中气预发[2000] 11号文标准，对紫外线指数进行等级划分，有95 %的紫外线等级估测与实测值相差不大于1，估测效果较好。

通过禹城站近7年观测数据获得了该区域紫外辐射的变化规律，并建立了利用温度日较差估测紫外辐射的方法，结论如下：

(1)紫外辐射与总辐射有一致的年、月、日变化规律，由于地球公转导致总辐射与紫外辐射都呈冬季小夏季大的季节分布特征，由于地球自转使得总辐射和紫外辐射都呈早晚小中午大的日变化特点。紫外辐射日累计值基本维持在0.1～1.20MJ m⁻² d⁻¹，平均值为0.468 MJ m⁻² d⁻¹。

(2)紫外辐射与总辐射的比值范围在0.023～0.046，呈夏大冬小的特征，并且夏季的变化幅度大于冬季，该比值随晴空指数的增大而减小，在晴空指数大于0.5时比较稳定。

(3)可以由温度日较差估测总辐射，R²达0.80，平均相对误差为0.19，再由总辐射计算晴空指数K_s，结合太阳高度角和可照时数估测紫外辐射日累计值，R²达0.73，相对均方根误差为0.22，估测紫外线等级与实测紫外线等级相差不大于1的数据占95%，估测效果较好。

(4)影响到达地面总辐射的因子除了气温、水汽外还有其他很多因子，本文考虑得不全面，气温日较差模型准确率需要进一步改进。