黄土高原半干旱区典型日吸收性气溶胶综合观测分析

利用兰州大学半干旱气候与环境观测站的太阳光度计、激光雷达、微波辐射计综合观测资料,结合辐射传输模式分析了该地区秋季典型日2012年9月3~4日、21日和28日气溶胶物理特性、垂直分布特征,及其与气象条件的关系。研究时期的气溶胶主要为局地沙尘与人为污染混合气溶胶,吸收性明显,尺度较小。其中,4日西北风增强,远距离传输沙尘气溶胶,气溶胶光学厚度最大,粒子尺度明显增大。尝试利用灰色关联度法确定参考高度,分别为7.41 km、8.47 km、7.13 km和7.66 km,反演气溶胶消光系数,由此积分得到的光学厚度与太阳光度计观测值相关性可达0.975,反演效果较好。研究时期气溶胶的抬升主要受白天热力湍流作用,边界层发展,气溶胶向上传输,每日12时（当地时间,下同）至14时传输至最大高度,气溶胶抬升的高度对应大气加热率的高值区,低层加热率可达1 K d-1。气溶胶在大气层顶和地面造成负辐射强迫,分别为-12.707 W m-2、-25.398 W m-2,大气中表现为正辐射强迫,为12.692 W m-2,大气层顶的辐射强迫对气溶胶的物理特性最为敏感,当气溶胶吸收性明显时,大气层顶的瞬时辐射强迫会出现正值。     

中国地区夏季平均加热率的时空分布特征

The latitude-altitude distributions of radiative fluxes and heating rates are investigated by utilizing CloudSat satellite data over China during summer. The Tibetan Plateau causes the downward shortwave fluxes of the lower atmosphere over central China to be smaller than the fluxes over southern and northern China by generating more clouds. The existence of a larger quantity of clouds over central China reflects a greater amount of solar radiation back into space. The vertical gradients of upward shortwave radiative fluxes in the atmosphere below 8 km are greater than those above 8 km. The latitudinal-altitude distributions of downward longwave radiative fluxes show a slantwise decreasing trend from low latitudes to high latitudes that gradually weaken in the downward direction. The upward longwave radiative fluxes also weaken in the upward direction but with larger gradients. The maximum heating rates by solar radiation and cooling rates by longwave infrared radiation are located over 28-40°N at 7-8 km mean sea level (MSL), and they are larger than the rates in the northern and southern regions. The heating and cooling rates match well both vertically and geographically.     

从Himawari08卫星估算晴空地表长波辐射及其日变化特征初探

通过446183条全球晴空大气廓线的红外辐射传输模拟和统计回归,建立了由Himawari08成像仪通道遥测数据估算晴空地表上行、下行长波辐射通量的反演模式,模式应用于成像仪观测资料,处理出晴空地表上行、下行长波辐射通量实时产品,2016年2~6月的产品精度验证试验结果为:与相同时刻的AQUA卫星CERES仪器同类产品相比,地表上行通量均方根误差R_e=7.9 W/m2,相关系数R=0.9399,地表下行通量R_e=14.5 W/m2,R=0.9586;与由中国地面气象站地面气温和相对湿度观测经Brunt、Brutsaert经验公式计算的实时地表下行长波辐射通量相比,R_e=15.34 W/m2,R=0.8786;与用陆表温度计算的地表上行长波辐射通量相比,R_e=12.6 W/m2,R=0.9977。研究了2016年2、6月的晴空地表长波辐射产品,发现陆地晴空上、下行通量有着与太阳加热地表增温相应的明显日变化特征,峰值出现在12:00（当地时间,下同）至14:00,低谷出现在04:00至07:00,下行通量与上行通量几乎同步变化或约有延时,陆地上2个通量归一化的日变化指数类似一个半正弦曲线,而海面长波辐射通量则没有明显的日变化规律。     

不同形状冰晶权重假定对冰云光学和辐射特性的影响

在BCC_RAD辐射传输模式和包含多形状冰晶粒子的冰云光学性质参数化方案的基础上,详细分析了不同冰晶粒子权重选取对冰云光学和辐射特性的影响。结果显示,不同形状冰晶粒子权重的选取对长波带平均消光系数、单次散射比、不对称因子和短波带平均不对称因子均有较大的影响。冰晶粒子权重选取对长波辐射通量有很大影响:对长波向下辐射通量,权重选择不同可在云底处造成高达10.50 W/m2的差别;对长波向上辐射通量,权重选择不同可在云顶处造成高达15.05 W/m2的差别。冰晶粒子权重选择对短波辐射通量也存在较大影响:对短波向下辐射通量,权重选择不同可在云底处造成高达12.48 W/m2的差别;对短波向上辐射通量,权重选择不同可在云顶处造成高达10.23 W/m2的差别。冰晶粒子权重选择对长波加热率影响较大,在云顶处和云底处分别可达1.31和-2.06 K/d。研究表明,不同形状冰晶粒子权重的选取对冰云光学性质和辐射计算均有较大的影响,在长波区间尤其明显。      

长波区间太阳辐射对气候模拟的影响

长波区间的太阳辐射在气候模式中往往被忽略。利用国家气候中心BCC_AGCM2.0.1大气环流模式,采用矩阵算子辐射传输算法,研究了长波区间太阳辐射对气候模式辐射通量和温度模拟结果的影响。结果表明,以ISCCP和CERES辐射资料为标准,考虑长波区间太阳辐射后,长波区间晴空大气地表向下辐射通量平均误差减小2.05 W/m2,均方根误差减少1.29 W/m2;长波区间晴空大气模式顶向上辐射通量平均误差减小0.70 W/m2,均方根误差减小0.21 W/m2;长波区间有云大气地表向下辐射通量平均误差减小1.38 W/m2,均方根误差减小1.03 W/m2;长波区间有云大气模式顶向上辐射通量平均误差减小0.99 W/m2,均方根误差减小0.30 W/m2。以ECMWF再分析资料为标准,考虑长波区间太阳辐射后,赤道地区上对流层—下平流层区域温度的冷偏差得到改善,对流层顶温度平均误差减小0.27 K,均方根误差减小0.25 K。     

Impact of 3-D topography on surface radiation budget over the Tibetan Plateau

The 3-D complex topography effect on the surface solar radiative budget over the Tibetan Plateau is investigated by means of a parameterization approach on the basis of “exact” 3-D Monte Carlo photon tracing simulations, which use 90 m topography data as building blocks. Using a demonstrative grid size of 10?×?10 km², we show that differences in downward surface solar fluxes for a clear sky without aerosols between the 3-D model and the conventional plane-parallel radiative transfer scheme are substantial, on the order of 200 W/m² at shaded or sunward slopes. Deviations in the reflected fluxes of the direct solar beam amount to about +100 W/m² over snow-covered areas, which would lead to an enhanced snowmelt if the 3-D topography effects had been accounted for in current climate models. We further demonstrate that the entire Tibetan Plateau would receive more solar flux by about 14 W/m², if its 3-D mountain structure was included in the calculations, which would result in larger sensible and latent heat transfer from the surface to the atmosphere.     

Simulation of Mt. Pinatubo Volcanic Aerosol with the Hamburg Climate Model ECHAM4

Summary  We study the three-dimensional transport of Mt. Pinatubo volcanic cloud with the climate model ECHAM4. In order to obtain model results comparable with observations a Newtonian relaxation technique was applied, which forces prognostic model variables towards the observations. A comparison of the simulated aerosol distribution with satellite data reveals good agreement for the first months after the eruption. The model, however, is unable to simulate the tropical aerosol maximum in 1992 and also overestimates the vertical downward and northward transport of aerosols. Substantial improvement was achieved with the introduction of reduced advective vertical transport through the 380 K isentropic layer. Heating rates and top of the atmosphere fluxes, which were calculated online for the first half year after the eruption, are in the observed range. A comparison of Pinatubo simulations between three different vertical ECHAM4 versions (ECHAM4 L19, ECHAM4 L39, MA/ECHAM4) indicates that a vertical resolution of ≈ 700 m in the tropopause region is sufficient to realistically reduce the vertical transport through the tropopause. Consideration of the upper branch of the Brewer Dobson circulation in the MA/ECHAM4 model improves the geographical distribution of the volcanic cloud. The application of a relaxation technique can further reduce major shortcomings of stratospheric simulations with the standard climate model. There remain, however some critical points in the global transport characteristics in all three models which are not fully understood. Received December 19, 1997 Revised July 22, 1998     

Direct Radiative Forcing of Anthropogenic Aerosols over Oceans from Satellite Observations

Anthropogenic aerosols play an important role in the atmospheric energy balance. Anthropogenic aerosol optical depth (AOD) and its accompanying shortwave radiative forcing (RF) are usually simulated by nu- merical models. Recently, with the development of space-borne instruments and sophisticated retrieval algorithms, it has become possible to estimate aerosol radiative forcing based on satellite observations. In this study, we have estimated shortwave direct radiative forcing due to anthropogenic aerosols over oceans in all-sky conditions by combining clouds and the Single Scanner Footprint data of the Clouds and Earth’s Radiant Energy System (CERES/SSF) experiment, which provide measurements of upward shortwave fluxes at the top of atmosphere, with Moderate Resolution Imaging Spectroradiometer (MODIS) aerosol and cloud products. We found that globally averaged aerosol radiative forcing over oceans in the clear-sky conditions and all-sky conditions were -1.03±0.48 W m-2 and -0.34 ±0.16 W m-2, respectively. Direct radiative forcing by anthropogenic aerosols shows large regional and seasonal variations. In some regions and in particular seasons, the magnitude of direct forcing by anthropogenic aerosols can be comparable to the forcing of greenhouse gases. However, it shows that aerosols caused the cooling effect, rather than warming effect from global scale, which is different from greenhouse gases.     

二流—四流球谐函数谱展开累加辐射传输方案在全球气候模式中的应用

本文首先构建了二流—四流球谐函数谱展开累加辐射传输的新方案,然后将其应用于国家气候中心第二代大气环流模式BCC_AGCM2.0.1的新版本中,并与模式中原有的Eddington累加方案进行了比较。由于新方案本质上是单层Eddington近似方案在四流上的推广。因此新方案在计算精度上要优于原方案。通过在全球气候模式中的应用与比较,本文发现新方案对气候模拟会产生比较大的影响。在晴空条件下,新方案计算的在南纬30°到60°区间、北大西洋东北部以及非洲北部的撒哈拉沙漠区域的地表向下年平均短波辐射通量要小于原方案结果,最大差别可以达到3.5 W/m2;同时,新方案计算的在南纬30°到60°区间和北大西洋东北部的大气顶向上年平均短波辐射通量要大于原方案结果,最大差别达到3 W/m2。在有云大气情况下,新方案计算的地表向下年平均短波辐射通量要小于原方案结果,并随着纬度的增加,新旧两种方案的差别逐渐变大,在南北极时达到最大5.5 W/m2;同时,新方案计算的在赤道区域的大气顶的年平均短波向上辐射通量要小于原方案结果,最大差别为2.5 W/m2,而在南北纬30°到60°区间,新方案计算的在大气顶的年平均短波向上辐射通量则要大于原方案结果,最大差别为1.5 W/m2。新方案计算的年平均短波加热率普遍高于原方案结果,特别是在800 hPa到地表之间的低层大气以及50 hPa到100 hPa的高层大气,最大差别可达0.03 K/d。因此,新方案有助于改善全球气候模式中普遍存在的赤道平流层中下层的温度冷偏差现象。     

The Arctic surface energy budget as simulated with the IPCC AR4 AOGCMs

Ensembles of simulations of the twentieth- and twentyfirst-century climate, performed with 20 coupled models for the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment, provide the basis for an evaluation of the Arctic (70°–90°N) surface energy budget. While the various observational sources used for validation contain differences among themselves, some model biases and across-model differences emerge. For all energy budget components in the twentieth-century simulations (the 20C3M simulation), the across-model variance and the differences from observational estimates are largest in the marginal ice zone (Barents, Kara, Chukchi Seas). Both downward and upward longwave radiation at the surface are underestimated in winter by many models, and the ensenmble mean annual net surface energy loss by longwave radiation is 35 W/m², which is less than for the NCEP and ERA40 reanalyses but in line with some of the satellite estimates. Incoming solar radiation is overestimated by the models in spring and underestimated in summer and autumn. The ensemble mean annual net surface energy gain by shortwave radiation is 39 W/m², which is slightly less than for the observational based estimates, In the twentyfirst-century simulations driven by the SRES A2 scenario, increased concentrations of greenhouse gasses increase (average for 2080–2100 minus average for 1980–2000 averages) the annual average ensemble mean downward longwave radiation by 30.1 W/m². This was partly counteracted by a 10.7 W/m² reduction in downward shortwave radiation. Enhanced sea ice melt and increased surface temperatures increase the annual surface upward longwave radiation by 27.1 W/m² and reduce the upward shortwave radiation by 13.2 W/m², giving an annual net (shortwave plus longwave) surface radiation increase of 5.8 W/m² , with the maximum changes in summer. The increase in net surface radiation is largely offset by an increased energy loss of 4.4 W/m² by the turbulent fluxes.     

The Effects of Different HITRAN Versions on Calculated Long-Wave Radiation and Uncertainty Evaluation

Four editions of the High Resolution Transmission (HITRAN) databases (HITRAN96, HITRAN2K, HITRAN04, and HITRAN08) are compared by using a line-by-line (LBL) radiative model in the long-wave calculation for six typical atmospheres. The results show that differences in downward radiative fluxes between HITRAN96 and HITRAN08 at the surface can reach a maximum of 1.70 W m 2 for tropical atmospheres. The largest difference in heating rate between HITRAN96 and HITRAN08 can reach 0.1 K day 1 for midlatitude summer atmosphere. Uncertainties caused by line intensity and air-broadened halfwidths are also evaluated in this work using the uncertainty codes given in HITRAN08. The uncertainty is found to be 1.92 W m 2 for upward fluxes at the top of the atmosphere (TOA) and 1.97 W m 2 for downward fluxes at the surface. The largest heating rate caused by the uncertainty of line intensity and air-broadened half-width can reach 0.5 K day 1 . The differences in optical depths between 1300 and 1700 cm 1 caused by different HITRAN versions are larger than those caused by the uncertainties in intensity and air-broadened half-width. This paper suggests that there is inaccurate representation of line parameters over some spectral ranges in HITRAN and more attention should be paid to these ranges in fields such as remote sensing.     

Influence of atmospheric relative humidity on ultraviolet flux and aerosol direct radiative forcing: Observation and simulation

The atmospheric aerosols can absorb moisture from the environment due to their hydrophilicity and thus affect atmospheric radiation fluxes. In this article, the ultraviolet radiation and relative humidity (RH) data from ground observations and a radiative transfer model were used to examine the influence of RH on ultraviolet radiation flux and aerosol direct radiative forcing under the clear-sky conditions. The results show that RH has a significant influence on ultraviolet radiation because of aerosol hygroscopicity. The relationship between attenuation rate and RH can be fitted logarithmically and all of the R² of the 4 sets of samples are high, i.e. 0.87, 0.96, 0.9, and 0.9, respectively. When the RH is 60%, 70%, 80% and 90%, the mean aerosol direct radiative forcing in ultraviolet is ?4.22W m^?2, ?4.5W m^?2, ?4.82W m^?2 and ?5.4W m^?2, respectively. For the selected polluted air samples the growth factor for computing aerosol direct radiative forcing in the ultraviolet for the RH of 80% varies from 1.19 to 1.53, with an average of 1.31.     

平流层气溶胶的辐射强迫及其气候响应的水平二维分析

利用比较先进的辐射模式计算了平流层气溶胶的辐射强迫,并对之进行了参数化。结果发现平流层气溶胶的辐射强迫的水平分布不仅与其本身的水平变化有关,而且与下垫面的反照率有很大的关系。利用近期开发的二维能量平衡模式模拟了皮纳图博火山气溶胶对地面平衡温度的影响,结果表明:皮纳图博火山至喷发后1年半左右降温达最大,至喷发后第5年降温已很小。     

The vertical distribution of climate forcings and feedbacks from the surface to top of atmosphere

The radiative forcings and feedbacks that determine Earth’s climate sensitivity are typically defined at the top-of-atmosphere (TOA) or tropopause, yet climate sensitivity itself refers to a change in temperature at the surface. In this paper, we describe how TOA radiative perturbations translate into surface temperature changes. It is shown using first principles that radiation changes at the TOA can be equated with the change in energy stored by the oceans and land surface. This ocean and land heat uptake in turn involves an adjustment of the surface radiative and non-radiative energy fluxes, with the latter being comprised of the turbulent exchange of latent and sensible heat between the surface and atmosphere. We employ the radiative kernel technique to decompose TOA radiative feedbacks in the IPCC Fourth Assessment Report climate models into components associated with changes in radiative heating of the atmosphere and of the surface. (We consider the equilibrium response of atmosphere-mixed layer ocean models subjected to an instantaneous doubling of atmospheric CO₂). It is shown that most feedbacks, i.e., the temperature, water vapor and cloud feedbacks, (as well as CO₂ forcing) affect primarily the turbulent energy exchange at the surface rather than the radiative energy exchange. Specifically, the temperature feedback increases the surface turbulent (radiative) energy loss by 2.87 W m^?2 K^?1 (0.60 W m^?2 K^?1) in the multimodel mean; the water vapor feedback decreases the surface turbulent energy loss by 1.07 W m^?2 K^?1 and increases the surface radiative heating by 0.89 W m^?2 K^?1; and the cloud feedback decreases both the turbulent energy loss and the radiative heating at the surface by 0.43 and 0.24 W m^?2 K^?1, respectively. Since changes to the surface turbulent energy exchange are dominated in the global mean sense by changes in surface evaporation, these results serve to highlight the fundamental importance of the global water cycle to Earth’s climate sensitivity.     

中国东部大气温湿廓线特征及其对辐射收支计算影响的分析

辐射平衡是气候决定性因素之一, 温湿廓线对辐射计算有重要影响。利用中国东部地区1960—2008年74个探空站资料, 将SBDART（Santa Barbara DISORT Atmospheric Radiative Transfer）辐射传输模式自带的大气温湿廓线和探空测得的实际大气温湿廓线进行对比, 并代入SBDART模式中, 分别计算晴空条件下不同温湿廓线对辐射通量的影响。结果表明:（1）中国东部地区的大气温湿廓线有明显的区域和季节特点;（2）模式地面向下辐射通量与实际相比有较明显差异, 尤其是在东北地区的冬季、中东部地区的夏季和华南地区的冬季, 差值达20—35 W/m2, 相对误差达2.01%—3.18%;（3）在东北地区的冬季、中东部地区的冬季以及华南地区的夏、冬季, 模式计算的大气顶向上辐射通量与实际相比差值达10—22 W/m2, 相对误差可达3.67%—8.94%;（4）模式与实际辐射加热率的差值在0.03—0.29 K/d。研究表明, 模式自带的大气温湿廓线区域和季节划分并不够细致, 不足以代表中国东部地区各个区域及季节的温湿特点。建立一套中国东部地区的大气标准廓线可以为辐射模拟提供更准确的输入量。      

2011年2月—2012年1月东南极高原辐射平衡观测研究

利用东南极高原熊猫-1自动气象站2011年2月—2012年1月观测的辐射资料和相关资料,对辐射分量和辐射平衡的季节变化进行了研究。结果表明,夏季是东南极高原获得太阳能的主要时段,总辐射通量夏季平均为365.0 W/m²,总量达到2752.1 MJ/m²,占全年总辐射量的58%。各个季节均能出现总辐射瞬时值大于大气顶水平总辐射,春季发生频率最高,冬季最小,总辐射平均日变化呈单峰型。大气长波辐射除夏季外,日变化不明显。冰雪面长波辐射除冬季外,各季节平均日变化呈明显的单峰单谷型。净辐射12月和1月为很小的正值,其他月份为负值。年平均净辐射为 -8.7 W/m²,表明地表相对于大气为冷源。该站的辐射平衡特征与其他南极内陆高原站相似,雪面具有强烈的辐射冷却效应,导致净辐射绝对值都小于下降风区。     

火山气溶胶的辐射影响

本文应用辐射传输计算研究了火山气溶胶,特别是大的火山爆发后短时期及火山周围地区内的火山气溶胶,对到达地面的向下的总的太阳辐射通量、行星反照率及大气加热率的可能影响.     

The energy budget of the Sable Island ocean region

The heat exchange between ocean and atmosphere over cold water is studied by calculating all terms in the energy balance twice each day for the year 1971 for the Sable Island region.

The atmospheric long‐wave radiation is relatively constant because of frequent overcast and low clouds. The surface long‐wave balance is markedly negative in winter but slightly positive for a short time in summer, due to strong advection of warm moist air over the cold water. In winter, the turbulent fluxes are directed upwards and are strong, the upward fluxes beginning after the middle of August and lasting until mid‐March. The maximum daily values of latent heat flux are 400 to 500 ly day^?1 (194 to 242 W m^?2), about a third or a quarter of the magnitude over the warmer Gulf Stream water. The summer fluxes are fairly constant and directed downward.

The water of the Labrador Current in the Sable Island region warms substantially from March to September and conversely cools intensely in the period November‐January.

A comparison of the energy exchange for a current and for water without motion shows that the surface temperatures would be similar in summer, and the temperature drop would be about equal until November. From that time on, the surface temperature would level off for a water body with no current, but in actual conditions the surface temperature continues to drop to a late winter minimum of about 1°C.

Atmospheric advection of latent heat was calculated by assuming that the daily precipitation was always caused first by condensation of all locally evaporated water with any remainder being supplied by water‐vapour advection. The main cause for atmospheric heating in the Sable Island area was found to be condensation of imported water vapour. The region is, in summer, a marked sink for atmospheric heat and water content. For water it remains a sink even in winter. For sensible heat it becomes a source from November to March. The warming of the atmosphere is caused by release of latent heat of advected water vapour in the period February‐August. During the months September‐January the heat sources are both water‐vapour advection and surface turbulent terms.     

A review of global stratospheric aerosol: Measurements, importance, life cycle, and local stratospheric aerosol

Stratospheric aerosol, noted after large volcanic eruptions since at least the late 1800s, were first measured in the late 1950s, with the modern continuous record beginning in the 1970s. Stratospheric aerosol, both volcanic and non-volcanic are sulfuric acid droplets with radii (concentrations) on the order of 0.1–0.5 µm (0.5–0.005 cm^− 3), increasing by factors of 2–4 (10–10³) after large volcanic eruptions. The source of the sulfur for the aerosol is either through direct injection from sulfur-rich volcanic eruptions, or from tropical injection of tropospheric air containing OCS, SO₂, and sulfate particles. The life cycle of non-volcanic stratospheric aerosol, consisting of photo-dissociation and oxidation of sulfur source gases, nucleation/condensation in the tropics, transport pole-ward and downward in the global planetary wave driven tropical pump, leads to a quasi steady state relative maximum in particle number concentration at around 20 km in the mid latitudes. Stratospheric aerosol have significant impacts on the Earth's radiation balance for several years following volcanic eruptions. Away from large eruptions, the direct radiation impact is small and well characterized; however, these particles also may play a role in the nucleation of near tropopause cirrus, and thus indirectly affect radiation. Stratospheric aerosol play a larger role in the chemical, particularly ozone, balance of the stratosphere. In the mid latitudes they interact with both nitrous oxides and chlorine reservoirs, thus indirectly affecting ozone. In the polar regions they provide condensation sites for polar stratospheric clouds which then provide the surfaces necessary to convert inactive to active chlorine leading to polar ozone loss. Until the mid 1990s the modern record has been dominated by three large sulfur-rich eruptions: Fuego (1974), El Chichón (1982) and Pinatubo (1991), thus definitive conclusions concerning the trend of non-volcanic stratospheric aerosol could only recently be made. Although anthropogenic emissions of SO₂ have changed somewhat over the past 30 years, the measurements during volcanically quiescent periods indicate no long term trend in non-volcanic stratospheric aerosol.     

Influences of the Internal Mixing of Anthropogenic Aerosols on Global Aridity Change

Influences of the mixing treatments of anthropogenic aerosols on their effective radiative forcing (ERF) and global aridity are evaluated by using the BCC_AGCM2.0_CUACE/Aero, an aerosol–climate online coupled model. Simulations show that the negative ERF due to external mixing (EM, a scheme in which all aerosol particles are treated as independent spheres formed by single substance) aerosols is largely reduced by the partial internal mixing (PIM, a scheme in which some of the aerosol particles are formed by one absorptive and one scattering substance) method. Compared to EM, PIM aerosols have much stronger absorptive ability and generally weaker hygroscopicity, which would lead to changes in radiative forcing, hence to climate. For the global mean values, the ERFs due to anthropogenic aerosols since the pre-industrial are–1.02 and–1.68 W m^–2 for PIM and EM schemes, respectively. The variables related to aridity such as global mean temperature, net radiation flux at the surface, and the potential evaporation capacity are all decreased by 2.18/1.61 K, 5.06/3.90 W m^–2, and 0.21/0.14 mm day^–1 since 1850 for EM and PIM schemes, respectively. According to the changes in aridity index, the anthropogenic aerosols have caused general humidification over central Asia, South America, Africa, and Australia, but great aridification over eastern China and the Tibetan Plateau since the pre-industrial in both mixing schemes. However, the aridification is considerably alleviated in China, but intensified in the Arabian Peninsula and East Africa in the PIM scheme.