北京地面紫外辐射（光谱）的观测与分析

通过对北京大气物理研究所与长春光机所合作研制的地基太阳紫外辐射光谱仪观测资料的分析，和用辐射传输模式UVSS计算的结果，对影响到达地面的UVB辐照度的主要因子太阳高度角（SZA）、臭氧总量和地表反照率进行了分析研究。最后对紫外光谱仪的观测资料做了总量和谱分析。     

基于CMIP6强迫模拟分析人为气溶胶的气候效应(一)——介绍NUIST模式评估结果

为减少不同气候模式评估气溶胶气候效应的差异,第六次耦合模式比较计划(Coupled Model Intercomparison Project Phase 6,CMIP6)直接给定了人为气溶胶强迫数据。因此,有必要基于此强迫数据重新评估气溶胶气候效应。本研究首先将CMIP6给出的描述人为气溶胶强迫的模块引入南京信息工程大学(Nanjing University of Information Science and Technology,NUIST)的地球系统模式(The NUIST Earth System Model,NESM)。之后,利用NESM模式评估地球辐射收支平衡对此人为气溶胶强迫的响应,并分析模式模拟结果的不确定性。评估给出的人为气溶胶有效辐射强迫为-0. 45(±0. 28) W·m~(-2)。其中,气溶胶直接辐射效应为-0. 34(±0. 01) W·m~(-2),与第二次气溶胶比较计划(The second phase of Aerosol Comparisons between Observations and M odels,Aero ComⅡ)的评估结果基本一致;气溶胶对云辐射强迫的影响(包括半直接效应和间接效应)为-0. 10(±0. 30) W·m~(-2),明显受到模式内部变率的干扰,具有较大的不确定性。     

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.     

高云和气溶胶辐射效应对边界层的影响

通过WRF(Weather Research and Forecasting) 模式嵌套包含了高云和气溶胶辐射效应的大气边界层模式, 结合激光雷达资料, 进行数值模拟, 定量分析高云和气溶胶辐射效应对城市边界层的影响。模式能很好地模拟出在高云和气溶胶辐射效应下城市边界层的特征。夜间, 气溶胶在低层起到保温作用, 高云使得保温作用得到加强, 地表增温达1.5 K。中高层, 气溶胶降低所在气层温度, 高云使得降温幅度减少, 降温达0.2～0.7 K。白天, 高云和气溶胶减少到达地面的太阳短波辐射, 导致低层温度降低, 地表降温达1.3 K。中高层, 气溶胶加热所在的气层, 高云使得这一增温幅度减少, 在500 m处增温最大, 达0.85 K。无论白天还是夜间, 气溶胶的辐射效应都会抵消一部分形成山谷风的热力条件, 使得中低层的风速减少, 这种影响在白天显得尤为明显。高云的存在使得这种抵消得到少量的补偿。     

On the use of a cloud modification factor for solar UV (290–385 nm) spectral range

Summary  Knowledge of ultraviolet radiation is necessary in different applications, in the absence of measurements, this radiometric flux must be estimated from available parameters. To compute this flux under all sky conditions one must consider the influence of clouds. Clouds are the largest modulators of the solar radiative flux reaching the Earth’s surface. The amount and type of cloud cover prevailing at a given time and location largely determines the amount and type of solar radiation received at the Earth’s surface. This cloud radiative effect is different for the different solar spectral bands. In this work, we analyse the cloud radiative effect over ultraviolet radiation (290–385 nm). This could be done by defining a cloud modification Factor. We have developed such cloud modification Factor considering two different types of clouds. The efficiency of the cloud radiative effect scheme has been tested in combination with a cloudless sky empirical model using independent data sets. The performance of the model has been tested in relation to its predictive capability of global ultraviolet radiation. For this purpose, data recorded at two radiometric stations are used. The first one is located at the University of Almería, a seashore location (36.83° N, 2.41° W, 20 m a.m.s.l.), while the second one is located at Granada (37.18° N, 3.58° W, 660 m a.m.s.l.), an inland location. The database includes hourly values of the relevant variables that cover the years 1993–94 in Almería and 1994–95 in Granada. Cloud cover information provided by the Spanish Meteorological Service has been include to compute the clouds radiative effect. After our study, it appears that the combination of an appropriate cloudless sky model with the cloud modification Factor scheme provides estimates of ultraviolet radiation with mean bias deviation of about 5% that is close to experimental errors. Comparisons with similar formulations of the cloud radiative effect over the whole solar spectrum provides evidence for the spectral dependency of the cloud radiative effect. Received November 15, 1999 Revised September 11, 2000     

积云对流和云物理过程调整对气候模拟的影响

本文利用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室 (LASG/IAP) 的大气环流谱模式SAMIL, 结合观测与政府间气候变化专门委员会第四次评估报告 (IPCC AR4) 大气模式集合平均结果, 以大气辐射通量为例, 诊断分析了物理过程调整前后模式对气候模拟的影响。旧版本SAMIL对大气辐射通量的模拟存在较大偏差, 经过大气辐射过程、 积云对流和诊断云等物理过程的调整后, 新版本SAMIL模拟的全球辐射通量的年平均结果与观测的偏差大幅减小, 其中大气顶能量收支的年变化及其平均值与观测更为接近。在积云对流方案调整基础上, 通过对诊断云物理方案的进一步调整, 新版本SAMIL对云物理量模拟更为合理, 在赤道辐合带等区域, 在很大程度上克服了单一积云对流物理过程调整引起的云宏观和微观属性不匹配问题, 能模拟出夏季气候平均辐射通量的全球分布特征, 尤其在东亚区域有较好的模拟能力。研究还表明, 在热带和副热带对流活跃区域, 当前SAMIL中积云对流过程偏差对辐射通量的模拟偏差有很大影响, 而模式中较为简单的诊断云方案也会将云宏观物理量模拟偏差带入云微观量模拟中, 也是主要偏差源之一。本文结果表明, 要继续提高SAMIL的模拟性能, 急需更新云物理参数化方案以改进云辐射过程的模拟, 同时也需要有针对性的研究积云对流和云物理过程之间相互作用, 并作进一步协同调整。     

晴空或少云状况下紫外辐射强度及指数预报模式

在大气辐射传输理论的基础上，利用中分辨率大气辐射传输系统，应用可测得的实际大气物理参数，建立睛空或少云天气状况下石家庄市紫外波段（280-400nm)辐射强度及指数预报模式，考虑了臭氧的两个吸收带。晴空条件下的散射主要考虑了分子和气溶胶的多次散射，少云时还考虑了云对紫外辐射的影响。其中多次散射的计算采用子离散坐标法。且对由于臭氧实时资料的短缺造成的误差进行了系统订正，并将订正结果与实况资料进行了对比。结果表明，由模式客观预报紫外辐射强度是可行的。     

Advances in land modeling of KIAPS based on the Noah Land Surface Model

As of 2013, the Noah Land Surface Model (LSM) version 2.7.1 was implemented in a new global model being developed at the Korea Institute of Atmospheric Prediction Systems (KIAPS). This land surface scheme is further refined in two aspects, by adding new physical processes and by updating surface input parameters. Thus, the treatment of glacier land, sea ice, and snow cover are addressed more realistically. Inconsistencies in the amount of absorbed solar flux at ground level by the land surface and radiative processes are rectified. In addition, new parameters are available by using 1-km land cover data, which had usually not been possible at a global scale. Land surface albedo/emissivity climatology is newly created using Moderate-Resolution Imaging Spectroradiometer (MODIS) satellitebased data and adjusted parameterization. These updates have been applied to the KIAPS-developed model and generally provide a positive impact on near-surface weather forecasting.     

Influences of Physical Processes and Parameters on Simulations of TOA Radiance at UV Wavelengths: Implications for Satellite UV Instrument Validation

Numerous factors can influence the radiative transfer simulation of hyper-spectral ultraviolet satellite observation,including the radiative transfer scheme, gaseous absorption coefficients, Rayleigh scattering scheme, surface reflectance, aerosol scattering, band center wavelength shifts of sensor, and accuracy of input profiles. In this study, a Unified Linearized Vector Radiative Transfer Model(UNL-VRTM) is used to understand the influences of various factors on the top of atmosphere(TOA) normalized radiance in the ultraviolet(UV) region. A benchmark test for Rayleigh scattering is first performed to verify the UNL-VRTM accuracy, showing that the model performances agree well with earlier peer-reviewed results. Sensitivity experiments show that a scalar radiative transfer approximation considering only ozone and a constant surface reflectance within the UV region may cause significant errors to the TOA normalized radiance. A comparison of the Ozone Mapping and Profiler Suite(OMPS) radiances between simulations and observations shows that the surface reflectance strongly influences the accuracy for the wavelengths larger than 340 nm. Thus, using the surface reflectivity at 331 nm as a proxy for simulating the whole OMPS hyperspectral ultraviolet radiances is problematic. The impact of rotational Raman scattering on TOA radiance can be simulated through using SCIATRAN, which can also reduce the difference between measurements and simulations to some extent. Overall, the differences between OMPS simulations and observations can be less than 3% for the entire wavelengths. The bias is nearly constant across the cross-track direction.     

A method for obtaining pre-twentieth century initial conditions for use in climate change studies

A method is proposed to initialise coupled atmosphere-ocean general circulation models (AOGCMs) developed to study climate change on multi-century time scales. The method assumes that current generation AOGCMs are developed and evaluated using present-day radiative forcing and near present day oceanic initial conditions. To find pre-twentieth century initial conditions, we propose that the radiative forcing be run backwards in time from the present to the desired starting date. The model should then be run for 3–5 centuries with the radiative forcing held constant at the desired date. In our tests, instantaneously switching to pre-twentieth century radiative forcing did not save computational time. When a sufficiently stable pre-twentieth century condition is achieved, the coupled system can be integrated forward to the present and into the future. This method is a first step toward the standardization of AOGCM initialization and suggests a framework for AOGCM initialization for the first time. It provides an internally consistent set of pre-twentieth century initial conditions, although they will vary from model to model. Furthermore, it is likely that this method will yield a fairly realistic present-day climate in transient climate change experiments of the twentieth century, if the model biases are not too large. The main disadvantage of the method is that it is fairly computationally expensive in that it requires an additional 4–6 centuries of model integration before starting historical twentieth century integrations. However, the relative cost of this technique diminishes as more simulations are conducted using the oceanic initial condition obtained using our method.     

海气耦合模式FGOALS_gl模拟的水汽和云辐射反馈过程

本文分析了中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室(LASG)发展的快速耦合模式FGOALS_gl对低纬太平洋区域水汽温室效应(Ga)和云辐射强迫特征的模拟能力,讨论了模拟偏差的成因.结果表明,FGOALS_gl能合理再现Ga、云辐射强迫的气候态空间分布特征,但也存在明显的偏差.模式低估了冷...     

基于FY-3/IRAS利用非线性模式反演OLR

FY-3系列卫星星载IRAS仪器设有26个通道,其中20个通道用于探测地球大气在红外波段的热辐射,通道辐射率代表了地球大气系统在大气顶的向外辐射光谱信息,与总波段的射出长波辐射（OLR）通量相关性高。该文基于逐线辐射传输模式计算软件LBLRTM对全球2521条大气廓线的大气顶射出辐射率模拟数据,计算了每条廓线的OLR和FY-3B/IRAS,FY-3C/IRAS通道辐射率,用统计回归方法建立了利用IRAS的多通道辐射率计算OLR的非线性理论回归模式;应用模式和FY-3B/IRAS,FY-3C/IRAS的L1级数据,处理得到2016年4月1-30日的全球日平均、月平均OLR格点产品。与Aqua/CERES,Terra/CERES仪器宽波段观测OLR产品对比表明:对于水平分辨率为1°×1°的全球月平均OLR格点产品,均方根误差为2.22 W·m-2,相关系数为0.9982 W·m-2,平均偏差为-0.2 W·m-2,表明FY-3/IRAS仪器定标及反演模式均达到较高水平。文中还回顾了历史上不同气象卫星的多种OLR反演算法模式,并对不同模式精度进行了比较。     

沙尘辐射效应对天气和沙尘输送影响的数值模拟

现有的沙尘气溶胶数值预报模式大多是单向模式, 即只考虑气象场对沙尘气溶胶的输送作用, 而忽略沙尘气溶胶的辐射效应对气象场的反馈作用。本文利用一个考虑了沙尘气溶胶短波辐射效应的天气－沙尘双向模式, 对南疆盆地一次冷空气东侵而引起的沙尘天气进行数值模拟。通过比较考虑与不考虑辐射效应的两组试验, 分析了在沙尘暴－强浮尘过程中沙尘辐射效应对气象要素和沙尘输送的影响, 结果表明: 沙尘辐射效应的考虑明显改进后期强浮尘天气过程中海平面气压变化、地面气温日变化和气温垂直变化的预报; 前期沙尘暴过程沙尘辐射效应引起中低层主导风增大, 但同时使起沙量降低及沉降量减少; 沙尘辐射效应与地形的共同作用使沙尘的沉降量比起沙量减少的更显著, 从而使盆地内浮尘天气维持。     

Research on the Application of the Radiative Transfer Model Based on Deep Neural Network in One-dimensional Variational Algorithm

As a typical physical retrieval algorithm for retrieving atmospheric parameters, one-dimensional variational(1 DVAR) algorithm is widely used in various climate and meteorological communities and enjoys an important position in the field of microwave remote sensing. Among algorithm parameters affecting the performance of the 1 DVAR algorithm, the accuracy of the microwave radiative transfer model for calculating the simulated brightness temperature is the fundamental constraint on the retrieval ...     

气候模式中云的垂直重叠及其辐射传输问题研究进展

概述了全球气候模式中云的垂直重叠的处理方法及其辐射物理过程的最新研究进展。从云垂直重叠模型的构造、模型在气候模式中的实现方式,得到与观测一致的云重叠结构所采用的数据和方法、重叠云的辐射传输等方面,给出了针对这一国际研究难点问题的最新研究进展。关于气候模式中云的垂直重叠问题的研究至今已取得了许多成果,表现在:重叠模型上有了更为科学的描述形式(如指数衰减重叠);重叠云的辐射传输也有了更快速的处理方法(如蒙特卡洛独立柱近似)并被广泛应用;连续的三维云遥感观测(如CloudSat/CALIPSO)和云分辨尺度的三维云模式的发展为在气候模式中精确描述云的垂直结构提供了丰富的观测资料和模式数据。但是,气候模式中现有的云重叠结构处理及其辐射传输方法还远不够完善,仍然存在很多没有解决的问题需要在未来进行探索。     

The influence of clouds on surface UV erythemal irradiance

The purpose of this study is to examine the effect of clouds on the ultraviolet erythemal irradiance. The study was developed at three stations in the Iberian Peninsula: Madrid and Murcia, using data recorded in the period 2000–2001, and Zaragoza, using data recorded in 2001. In order to determine the cloud effect on ultraviolet erythemal irradiance, we considered a cloud modification factor defined as the ratio between the measured values of ultraviolet erythemal irradiance and the corresponding clear-sky ultraviolet erythemal irradiance, which would be expected for the same time period and atmospheric conditions. The dependence of this cloud modification factor on total cloud amount, cloud type and solar elevation angle was investigated. The results suggest that the effect of cloud on ultraviolet erythemal irradiance can be parameterized in a simple way in terms of the cloud amount. Our results suggest that the same cloud modification factor model can be used at the three analysed locations estimating the ultraviolet erythemal irradiance with mean bias deviation (MBD) in the range of the expected experimental errors. This cloud modification factor is lower than that associated to the whole solar spectral range, indicating that the attenuation for the ultraviolet erythemal irradiance is lower than that associated to other solar spectral ranges. The cloud modification factor for ultraviolet erythemal irradiance presents dependence with solar elevation, with opposite dependencies with solar elevation for overcast and partial cloud cover conditions, a fact that can be explained in terms of the influence of reflection-enhancement of the ultraviolet irradiance in the last case. Concerning the influence of cloud type, a limited study of two cloud categories, low and medium level and high level, indicated that for overcast conditions, lower clouds presents an attenuation of ultraviolet erythemal irradiance 20% greater than that associated to high level clouds.     

Relative effects of solar and infrared radiative forcing in a mesoscale model

A new efficient parameterization scheme for solar short-wave radiative heating, as a component of the net radiative effects in the atmosphere, is tested in a three-dimensional mesoscale model. This model is designed with moist convective processes in mind, so that the radiative parameterization (solar plus thermal infrared) are interactive with the cloud field. Previous work by the authors with only an infrared scheme has demonstrated that cloud-radiation interactions are characterized by strong cloud-top cooling, leading to upper cloud-layer destabilization. The effects of including solar heating are to modulate the strength of the strong infrared cooling, thereby leading to weaker interactions between clouds, radiation, and mesoscale fields. The present study shows that even on the mesoscale and for relatively short time-spans, radiative processes in the presence of clouds are not negligible.As a further step, a simple fractional cloud cover parameterization is introduced and the model response is compared with results omitting this parameterization.     

GRAPES紫外线（UV）数值预报

应用GRAPES（Global/Regional Assi milation PrEdiction System）模式中的Goddard短波辐射方案,创建了紫外线数值预报系统（GRAPES-UV）。介绍了Goddard短波辐射方案,给出了GRAPES-UV系统的运行和个例分析。研究结果表明,紫外线指数（UVI）除了与纬度、地形和日变化有关外,还与云的分布以及天气形势密切相关,GRAPES模式中云的微物理方案对UV预报有较大的影响,UVI在晴空和对流云降水地区的强烈反差是UVI的重要特征。在个例试验中应用了臭氧总量预报模式,通过T213模式为化学输送模式提供气象背景,利用卫星资料同化技术建立臭氧的初始场,预报大气臭氧总量。应用国际上通常采用的临界成功指数（CSI）对2007年夏季北京和上海UV预报进行统计检验。北京和上海24 h紫外线强度等级为强和很强的CSI分别为0.625和0.780,接近同样方法的美国检验结果0.677。该系统从2006年3月起在中国气象局大气成分观测与服务中心的业务系统中正常运行,并在大气成分中心网站上发布UVI和紫外线强度等级预报,并提醒公众采取合适的方法保护自身免受UV的过度辐射。     

Temperature stabilization, ocean heat uptake and radiative forcing overshoot profiles

Political leaders in numerous nations argue for an upper limit of the global average surface temperature of 2 K above the pre-industrial level, in order to attempt to avoid the most serious impacts of climate change. This paper analyzes what this limit implies in terms of radiative forcing, emissions pathways and abatement costs, for a range of assumptions on rate of ocean heat uptake and climate sensitivity. The primary aim is to analyze the importance of ocean heat uptake for radiative forcing pathways that temporarily overshoot the long-run stabilization forcing, yet keep the temperature increase at or below the 2 K limit. In order to generate such pathways, an integrated climate-economy model, MiMiC, is used, in which the emissions pathways generated represent the least-cost solution of stabilizing the global average surface temperature at 2 K above the pre-industrial level. We find that the level of overshoot can be substantial. For example, the level of overshoot in radiative forcing in 2100 ranges from about 0.2 to 1 W/m2, where the value depends strongly and positively on the effective diffusivity of heat in the oceans. Measured in relative terms, the level of radiative forcing overshoot above its longrun equilibrium level in 2100 is 20% to 60% for high values of climate sensitivity (i.e., about 4.5 K) and 8% to 30% for low values of climate sensitivity (i.e., about 2 K). In addition, for cases in which the radiative forcing level can be directly stabilized at the equilibrium level associated with a specific climate sensitivity and the 2 K limit, the net present value abatement cost is roughly cut by half if overshoot pathways are considered instead of stabilization of radiative forcing at the equilibrium level without an overshoot.     

The UVic earth system climate model: Model description,climatology, and applications to past,present and future climates

Abstract

A new earth system climate model of intermediate complexity has been developed and its climatology compared to observations. The UVic Earth System Climate Model consists of a three‐dimensional ocean general circulation model coupled to a thermodynamic/dynamic sea‐ice model, an energy‐moisture balance atmospheric model with dynamical feedbacks, and a thermomechanical land‐ice model. In order to keep the model computationally efficient a reduced complexity atmosphere model is used. Atmospheric heat and freshwater transports are parametrized through Fickian diffusion, and precipitation is assumed to occur when the relative humidity is greater than 85%. Moisture transport can also be accomplished through advection if desired. Precipitation over land is assumed to return instantaneously to the ocean via one of 33 observed river drainage basins. Ice and snow albedo feedbacks are included in the coupled model by locally increasing the prescribed latitudinal profile of the planetary albedo. The atmospheric model includes a parametrization of water vapour/planetary longwave feedbacks, although the radiative forcing associated with changes in atmospheric CO₂ is prescribed as a modification of the planetary longwave radiative flux. A specified lapse rate is used to reduce the surface temperature over land where there is topography. The model uses prescribed present‐day winds in its climatology, although a dynamical wind feedback is included which exploits a latitudinally‐varying empirical relationship between atmospheric surface temperature and density. The ocean component of the coupled model is based on the Geophysical Fluid Dynamics Laboratory (GFDL) Modular Ocean Model 2.2, with a global resolution of 3.6° (zonal) by 1.8° (meridional) and 19 vertical levels, and includes an option for brine‐rejection parametrization. The sea‐ice component incorporates an elastic‐viscous‐plastic rheology to represent sea‐ice dynamics and various options for the representation of sea‐ice thermodynamics and thickness distribution. The systematic comparison of the coupled model with observations reveals good agreement, especially when moisture transport is accomplished through advection.

Global warming simulations conducted using the model to explore the role of moisture advection reveal a climate sensitivity of 3.0°C for a doubling of CO₂, in line with other more comprehensive coupled models. Moisture advection, together with the wind feedback, leads to a transient simulation in which the meridional overturning in the North Atlantic initially weakens, but is eventually re‐established to its initial strength once the radiative forcing is held fixed, as found in many coupled atmosphere General Circulation Models (GCMs). This is in contrast to experiments in which moisture transport is accomplished through diffusion whereby the overturning is reestablished to a strength that is greater than its initial condition.

When applied to the climate of the Last Glacial Maximum (LGM), the model obtains tropical cooling (30°N‐30°S), relative to the present, of about 2.1°C over the ocean and 3.6°C over the land. These are generally cooler than CLIMAP estimates, but not as cool as some other reconstructions. This moderate cooling is consistent with alkenone reconstructions and a low to medium climate sensitivity to perturbations in radiative forcing. An amplification of the cooling occurs in the North Atlantic due to the weakening of North Atlantic Deep Water formation. Concurrent with this weakening is a shallowing of, and a more northward penetration of, Antarctic Bottom Water.

Climate models are usually evaluated by spinning them up under perpetual present‐day forcing and comparing the model results with present‐day observations. Implicit in this approach is the assumption that the present‐day observations are in equilibrium with the present‐day radiative forcing. The comparison of a long transient integration (starting at 6 KBP), forced by changing radiative forcing (solar, CO₂, orbital), with an equilibrium integration reveals substantial differences. Relative to the climatology from the present‐day equilibrium integration, the global mean surface air and sea surface temperatures (SSTs) are 0.74°C and 0.55°C colder, respectively. Deep ocean temperatures are substantially cooler and southern hemisphere sea‐ice cover is 22% greater, although the North Atlantic conveyor remains remarkably stable in all cases. The differences are due to the long timescale memory of the deep ocean to climatic conditions which prevailed throughout the late Holocene. It is also demonstrated that a global warming simulation that starts from an equilibrium present‐day climate (cold start) underestimates the global temperature increase at 2100 by 13% when compared to a transient simulation, under historical solar, CO₂ and orbital forcing, that is also extended out to 2100. This is larger (13% compared to 9.8%) than the difference from an analogous transient experiment which does not include historical changes in solar forcing. These results suggest that those groups that do not account for solar forcing changes over the twentieth century may slightly underestimate (～3% in our model) the projected warming by the year 2100.