2011 Vol. 28, No. 4

Display Method:
2011, 28(4): 725-734. doi: 10.1007/s00376-010-9240-y
Abstract:
The present study examined the diurnal variations of summer precipitation in the Beijing area by using subdaily precipitation and wind observations. A combined effect of topography and urbanization on the characteristics of diurnal variations was suggested. It was shown that stations located in the plain area exhibited typical night rain peaks, whereas those in the mountainous area exhibited clear afternoon peaks of precipitation diurnal variations. The precipitation peaks were associated with wind fields around the Beijing area, which were found to be highly modulated by mountain-valley circulation and urban-country circulation. The lower-tropospheric wind exhibited a clear diurnal shift in its direction from north at 0800 LST to south at 2000 LST, which reflected mountain-valley circulation. The transitions from valley to mountain wind and the opposite generally happened after sunset and sunrise, respectively, and both occurred earlier for the stations located closer to mountains. By comparing the diurnal variations of precipitation at stations in a northeast suburb, an urban area, and a southwest suburb, it was revealed that the northeast suburb group had the highest normalized rainfall frequency, but the southwest group had the lowest from late afternoon to late evening. On the contrary, in the early morning from about 0200 to 1000 LST, the southwest group and urban group had the highest normalized rainfall frequency. This pattern might originate from the combined effects of mountain-valley topography and urbanization.
2011, 28(4): 735-742. doi: 10.1007/s00376-010-0061-9
Abstract:
Accurate forecast of rainstorms associated with the mei-yu front has been an important issue for the Chinese economy and society. In July 1998 a heavy rainstorm hit the Yangzi River valley and received widespread attention from the public because it caused catastrophic damage in China. Several numerical studies have shown that many forecast models, including Pennsylvania State University National Center for Atmospheric Research's fifth-generation mesoscale model (MM5), failed to simulate the heavy precipitation over the Yangzi River valley. This study demonstrates that with the optimal initial conditions from the dimension-reduced projection four-dimensional variational data assimilation (DRP-4DVar) system, MM5 can successfully reproduce these observed rainfall amounts and can capture many important mesoscale features, including the southwestward shear line and the low-level jet stream. The study also indicates that the failure of previous forecasts can be mainly attributed to the lack of mesoscale details in the initial conditions of the models.
2011, 28(4): 743-757. doi: 10.1007/s00376-010-0021-4
Abstract:
By sampling perturbed state vectors from each ensemble prediction run at properly selected time levels in the vicinity of the analysis time, the recently proposed time-expanded sampling approach can enlarge the ensemble size without increasing the number of prediction runs and, hence, can reduce the computational cost of an ensemble-based filter. In this study, this approach is tested for the first time with real radar data from a tornadic thunderstorm. In particular, four assimilation experiments were performed to test the time-expanded sampling method against the conventional ensemble sampling method used by ensemble-based filters. In these experiments, the ensemble square-root filter (EnSRF) was used with 45 ensemble members generated by the time-expanded sampling and conventional sampling from 15 and 45 prediction runs, respectively, and quality-controlled radar data were compressed into super-observations with properly reduced spatial resolutions to improve the EnSRF performances. The results show that the time-expanded sampling approach not only can reduce the computational cost but also can improve the accuracy of the analysis, especially when the ensemble size is severely limited due to computational constraints for real-radar data assimilation. These potential merits are consistent with those previously demonstrated by assimilation experiments with simulated data.
2011, 28(4): 758-774. doi: 10.1007/s00376-010-0145-6
Abstract:
Satellite data obtained over synoptic data-sparse regions such as an ocean contribute toward improving the quality of the initial state of limited-area models. Background error covariances are crucial to the proper distribution of satellite-observed information in variational data assimilation. In the NMC (National Meteorological Center) method, background error covariances are underestimated over data-sparse regions such as an ocean because of small differences between different forecast times. Thus, it is necessary to reconstruct and tune the background error covariances so as to maximize the usefulness of the satellite data for the initial state of limited-area models, especially over an ocean where there is a lack of conventional data. In this study, we attempted to estimate background error covariances so as to provide adequate error statistics for data-sparse regions by using ensemble forecasts of optimal perturbations using bred vectors. The background error covariances estimated by the ensemble method reduced the overestimation of error amplitude obtained by the NMC method. By employing an appropriate horizontal length scale to exclude spurious correlations, the ensemble method produced better results than the NMC method in the assimilation of retrieved satellite data. Because the ensemble method distributes observed information over a limited local area, it would be more useful in the analysis of high-resolution satellite data. Accordingly, the performance of forecast models can be improved over the area where the satellite data are assimilated.
2011, 28(4): 775-788. doi: 10.1007/s00376-010-0009-0
Abstract:
A dynamic global vegetation model (DGVM) coupled with a land surface model (LSM) is generally initialized using a spin-up process to derive a physically-consistent initial condition. Spin-up forcing, which is the atmospheric forcing used to drive the coupled model to equilibrium solutions in the spin-up process, varies across earlier studies. In the present study, the impact of the spin-up forcing in the initialization stage on the fractional coverages (FCs) of plant functional type (PFT) in the subsequent simulation stage are assessed in seven classic climate regions by a modified Community Land Model's Dynamic Global Vegetation Model (CLM-DGVM). Results show that the impact of spin-up forcing is considerable in all regions except the tropical rainforest climate region (TR) and the wet temperate climate region (WM). In the tropical monsoon climate region (TM), the TR and TM transition region (TR-TM), the dry temperate climate region (DM), the highland climate region (H), and the boreal forest climate region (BF), where FCs are affected by climate non-negligibly, the discrepancies in initial FCs, which represent long-term cumulative response of vegetation to different climate anomalies, are large. Moreover, the large discrepancies in initial FCs usually decay slowly because there are trees or shrubs in the five regions. The intrinsic growth timescales of FCs for tree PFTs and shrub PFTs are long, and the variation of FCs of tree PFTs or shrub PFTs can affect that of grass PFTs.
2011, 28(4): 789-796. doi: 10.1007/s00376-010-9222-0
Abstract:
Nepal lies on the southern slope of Himalaya in Asia. In a width ranging between 150 and 250 km, the altitude varies greatly from about 100 m at its southern border to a maximum of 8848 min the northern part. Like the variation in altitude, climatic condition varies quite a lot. Long-term monthly mean erythemal UV daily dose values for Nepal are evaluated using Total Ozone Mapping Spectrometer (TOMS) estimation from the time of its overpass between 1996 and 2003. The results are presented as summer and winter maps of mean UV levels in each satellite grid. The mean winter erythemal UV daily dose ranges between 2.1 and 3.6 kJ m-2 whereas summer values are found to lie between 4.6 and 9.7 kJ m-2. The altitude variation increases the UV levels by about 0.2 kJ km-1 in winter months, and 0.9 kJ km-1 in summer. A multiyear monthly average erythemal daily dose in most of the areas shows that the summer value is about three times higher than that in winter. Although year-to-year variation is not pronounced in high- and mid-elevation regions, UV levels seemed to decrease from 1997 to 2002 in the southern part of the country in the low elevation region by about 5.35%. Due to the combined effects of the altitude, low ozone concentration in the troposphere, and thin air, surface UV radiation at higher altitudes is found to be higher than in the surrounding regions.
2011, 28(4): 797-808. doi: 10.1007/s00376-010-0027-y
Abstract:
Based on the final analyses data (FNL) of the Global Forecasting System of the NCEP and the observational radiosonde data, the evolution mechanism of an eastward-moving low-level vortex over the Tibetan Plateau in June 2008 was analyzed. The results show that the formation of the vortex was related to the convergence between the northwesterly over the central Tibetan Plateau from the westerly zone and the southerly from the Bay of Bengal at 500 hPa, and also to the divergence associated with the entrance region of the upper westerly jet at 200 hPa. Their dynamic effects were favorable for ascending motion and forming the vortex over the Tibetan Plateau. Furthermore, the effect of the atmospheric heat source (Q1) is discussed based on a transformed potential vorticity (PV) tendency equation. By calculating the PV budgets, we showed that Q1 had a great influence on the intensity and moving direction of the vortex. In t Q1 he developing stage of the vortex, the heating of the vertically integrated was centered to the east of the vortex center at 500 hPa, increasing PV tendency to the east of the vortex. As a result, the vortex strengthened and moved eastward through the vertically uneven distribution of Q1. In the decaying stage, the horizontally uneven heating of Q1 at 500 hPa weakened the vortex through causing the vortex tubes around the vortex to slant and redistributing the vertical vorticity field.
2011, 28(4): 809-820. doi: 10.1007/s00376-010-9180-6
[Abstract](635) [PDF 352KB](1583)
Abstract:
Stochastic weather generators are statistical models that produce random numbers that resemble the observed weather data on which they have been fitted; they are widely used in meteorological and hydrological simulations. For modeling daily precipitation in weather generators, first-order Markov chain--dependent exponential, gamma, mixed-exponential, and lognormal distributions can be used. To examine the performance of these four distributions for precipitation simulation, they were fitted to observed data collected at 10 stations in the watershed of Yishu River. The parameters of these models were estimated using a maximum-likelihood technique performed using genetic algorithms. Parameters for each calendar month and the Fourier series describing parameters for the whole year were estimated separately. Bayesian information criterion, simulated monthly mean, maximum daily value, and variance were tested and compared to evaluate the fitness and performance of these models. The results indicate that the lognormal and mixed-exponential distributions give smaller BICs, but their stochastic simulations have overestimation and underestimation respectively, while the gamma and exponential distributions give larger BICs, but their stochastic simulations produced monthly mean precipitation very well. When these distributions were fitted using Fourier series, they all underestimated the above statistics for the months of June, July and August.
2011, 28(4): 821-831. doi: 10.1007/s00376-010-0007-2
Abstract:
Human-induced land use changes and the resulting alterations in vegetation features are major but poorly recognized drivers of regional climatic patterns. In order to investigate the impacts of anthropogenically-induced seasonal vegetation cover changes on regional climate in China, harmonic analysis is applied to 1982--2000 National Oceanic and Atmospheric Administration (NOAA) Advanced Very High Resolution Radiometer (AVVHRR)-derived normalized difference vegetation index (NDVI) time series (ten day interval data). For two climatic divisions of South China, it is shown that the first harmonic term is in phase with air temperature, while the second and third harmonics are in phase with agricultural cultivation. The Penman-Monteith Equation and the Complementary Relationship Areal Evapotranspiration (CRAE) model suggest that monthly mean evapotranspiration is out of phase with temperature and precipitation in regions with significant second or third harmonics. Finally, seasonal vegetation cover changes associated with agricultural cultivation are identified: for cropped areas, the temperature and precipitation time series have a single maximum value, while the monthly evapotranspiration time series has a bimodal distribution. It is hypothesized that multi-cropping causes the land surface albedo to sharply increase during harvesting, thereby altering the energy distribution ratio and contributing to observed seasonal vegetation cover changes.
2011, 28(4): 832-842. doi: 10.1007/s00376-010-0047-7
Abstract:
To analyze the mechanism by which water vapor increase leads to cooling in the stratosphere, the effects of water-vapor increases on temperature in the stratosphere were simulated using the two-dimensional, interactive chemical dynamical radiative model (SOCRATES) of NCAR. The results indicate that increases in stratospheric water vapor lead to stratospheric cooling, with the extent of cooling increasing with height, and that cooling in the middle stratosphere is stronger in Arctic regions. Analysis of the radiation process showed that infrared radiative cooling by water vapor is a pivotal factor in middle-lower stratospheric cooling. However, in the upper stratosphere (above 45 km), infrared radiation is not a factor in cooling; there, cooling is caused by the decreased solar radiative heating rate resulting from ozone decrease due to increased stratospheric water vapor. Dynamical cooling is important in the middle-upper stratosphere, and dynamical feedback to temperature change is more distinct in the Northern Hemisphere middle-high latitudes than in other regions and significantly affects temperature and ozone in winter over Arctic regions. Increasing stratospheric water vapor will strengthen ozone depletion through the chemical process. However, ozone will increase in the middle stratosphere. The change in ozone due to increasing water vapor has an important effect on the stratospheric temperature change.
2011, 28(4): 843-854. doi: 10.1007/s00376-010-0071-7
Abstract:
An increase in atmospheric nitrogen (N) deposition can promote soil acidification, which may increase the release of ethylene (C2H4) under forest floors. Unfortunately, knowledge of whether increasing N deposition and C2H4 releases have synergistic effects on soil methane (CH4) uptake is limited and certainly deserves to be examined. We conducted some field measurements and laboratory experiments to examine this issue. The addition of (NH4)2SO4 or NH4Cl at a rate of 45 kg N ha-1 yr-1 reduced the soil CH4 uptake under a temperate old-growth forest in northeast China, and there were synergistic effects of N amendments in the presence of C2H4 concentrations equal to atmospheric CH4 concentration on the soil CH4 uptake, particularly in the NH4Cl-treated plots. Effective concentrations of added C2H4 on the soil CH4 uptake were smaller in NH4+-treated plots than in KNO3-treated plots. The concentration of ca 0.3 mu l C2H4 L-1 in the headspace gases reduced by 20% soil atmospheric CH4 uptake in the NH4Cl-treated plots, and this concentration was easily produced in temperate forest topsoils under short-term anoxic conditions. Together with short-term stimulating effects of N amendments and soil acidification on C2H4 production from forest soils, our observations suggest that knowledge of synergistic effects of NH4+, rather than NO3-, amendments and C2H4 on the in situ soil CH4 uptake is critical for understanding the role of atmospheric N deposition and cycling of C2H4 under forest floors in reducing global atmospheric CH4 uptake by forests. Synergistic functions of NH4+-N deposition and C2H4 release due to soil acidification in reducing atmospheric CH4 uptake by forests are discussed.
2011, 28(4): 855-865. doi: 10.1007/s00376-010-0040-1
Abstract:
A strong urban heat island (UHI) appeared in a hot weather episode in Suzhou City during the period from 25 July to 1 August 2007. This paper analyzes the urban heat island characteristics of Suzhou City under this hot weather episode. Both meteorological station observations and MODIS satellite observations show a strong urban heat island in this area. The maximum UHI intensity in this hot weather episode is 2.2oC, which is much greater than the summer average of 1.0oC in this year and the 37--year (from 1970 to 2006) average of 0.35oC. The Weather Research and Forecasting (WRF) model simulation results demonstrate that the rapid urbanization processes in this area will enhance the UHI in intensity, horizontal distribution, and vertical extension. The UHI spatial distribution expands as the urban size increases. The vertical extension of UHI in the afternoon increases about 50 m higher under the year 2006 urban land cover than that under the 1986 urban land cover. The conversion from rural land use to urban land type also strengthens the local lake-land breeze circulations in this area and modifies the vertical wind speed field.
2011, 28(4): 866-878. doi: 10.1007/s00376-010-0001-8
Abstract:
Data from the Beijing SAFIR 3000 lightning detection system and Doppler radar provided some insights into the three-dimensional lightning structure and evolution of a leading-line and trailing-stratiform (LLTS) mesoscale convective system (MCS) over Beijing on 31 July 2007. Most of the lightning in the LLTS-MCS was intracloud (IC) lightning, while the mean ratio of positive cloud-to-ground (+CG) lightning to --CG lightning was 1:4, which was higher than the average value from previous studies. The majority of CG lightning occurred in the convective region of the radar echo, particularly at the leading edge of the front. Little IC lightning and little +CG lightning occurred in the stratiform region. The distribution of the CG lightning indicated that the storm had a tilted dipole structure given the wind shear or the tripole charge structure. During the storm's development, most of the IC lightning occurred at an altitude of ～9.5 km; the lightning rate reached its maximum at 10.5 km, the altitude of IC lightning in the mature stage of the storm. When the thunderstorm began to dissipate, the altitude of the IC lightning decreased gradually. The spatial distribution of lightning was well correlated with the rainfall on the ground, although the peak value of rainfall appeared 75 min later than the peak lightning rate.
2011, 28(4): 879-886. doi: 10.1007/s00376-010-0079-z
Abstract:
By analyzing the pollutant concentrations over the urban area and over the rural area of the city of Lanzhou, Gansu Province, China, the relationships between the daytime inversion intensity and the pollutant concentration in the atmospheric boundary layer (ABL) are studied with the consideration of wind speed and direction, potential temperature, specific humidity profiles, pollutant concentration in the ABL, the surface temperature, and global radiation on the ground. It was shown that the daytime inversion is a key factor in controlling air pollution concentration. A clear and positive feedback process between the daytime inversion intensity and the air pollutants over the city was found through the analysis of influences of climatic and environmental factors. The mechanisms by which the terrain and air pollutants affect the formation of the daytime inversion are discussed. The solar radiation as the essential energy source to maintain the inversion is analyzed, as are various out-forcing factors affecting the inversion and air pollutants. At last, a physical frame of relationships of air pollution with daytime inversion and the local and out-forcing factors over Lanzhou is built.
2011, 28(4): 887-906. doi: 10.1007/s00376-010-0041-0
Abstract:
Using land-use types derived from satellite remote sensing data collected by the EOS Moderate Resolution Imaging Spectroradiometer (EOS/MODIS), the mesoscale and turbulent fluxes generated by inhomogeneities of the underlying surface over the Jinta Oasis, northwestern China, were simulated using the Regional Atmospheric Modeling System (RAMS4.4). The results indicate that mesoscale circulation generated by land-surface inhomogeneities over the Jinta Oasis is more important than turbulence. Vertical heat fluxes and water vapor are transported to higher levels by mesoscale circulation. Mesoscale circulation also produces mesoscale synoptic systems and prevents water vapor over the oasis from running off. Mesoscale circulation transports moisture to higher atmospheric levels as the land-surface moisture over the oasis increases, favoring the formation of clouds, which sometimes leads to rainfall. Large-scale wind speed has a significant impact on mesoscale heat fluxes. During the active phase of mesoscale circulation, the stronger large-scale winds are associated with small mesoscale fluxes; however, background wind seems to intensify the turbulent sensible heat flux and turbulent latent heat flux. If the area of oasis is enlarged properly, mesoscale circulation will be able to transport moisture to higher levels, favoring the formation of rainfall in the oasis and protecting its cold island'' effect. The impact of irrigation on rainfall is important, and increasing irrigation across the oasis is necessary to protect the oasis.
2011, 28(4): 907-912. doi: 10.1007/s00376-010-0081-5
Abstract:
Atmospheric response to SST variability was estimated using generalized equilibrium feedback analysis (GEFA) in the SST EOF space with synthesis data from an idealized climate model. Results show that the GEFA atmospheric response to the leading SST EOF modes is much more accurate and robust than the GEFA feedback matrix in physical space. Therefore, GEFA provides a practical method for assessing atmospheric response to large-scale SST anomalies in terms of the leading EOFs.
2011, 28(4): 913-926. doi: 10.1007/s00376-010-0100-6
Abstract:
Atmospheric circulation cells associated with anomalous East Asian Winter Monsoon (EAWM) were studied using the 1948/49 to 2002/03 NCEP/NCAR reanalysis and NCAR CAM3 AGCM simulations with monthly global sea surface temperatures from 1950 to 2000. Several atmospheric cells in the Pacific [i.e., the zonal Walker cell (ZWC) in the tropic, the Hadley cell in the western Pacific (WPHC), the midlatitude zonal cell (MZC) over the central North Pacific, and the Hadley cell in the eastern Pacific (EPHC)] are associated with anomalous EAWM. When the EAWM is strong, ZWC, WPHC, and MZC are enhanced, as opposed to EPHC. The anomalous enhanced ZWC is characterized by air parcels rising in the western tropical Pacific, flowing eastward in the upper troposphere, and descending in the tropical central Pacific before returning to the tropical western Pacific. The enhanced MZC has characteristics opposite those of the enhanced ZWC in the central North Pacific. The anomalous WPHC shows air parcels rising in the western Pacific, as in the case of ZWC, followed by flowing northward in the upper troposphere and descending in the west North Pacific, as in the case of the enhanced MZC before returning to the western tropical Pacific. The anomalous EPHC is opposite in properties to the anomalous WPHC. Opposite characteristics are found during the weak EAWM period. The model simulations and the observations show similar characteristics and indicate the important role of sea surface temperature. A possible mechanism is proposed to link interannual variation of EAWM with the central-eastern tropical Pacific sea surface temperature anomaly (SSTA).
2011, 28(4): 927-944. doi: 10.1007/s00376-010-0105-1
Abstract:
In this study, the crop growth model (CERES3.0) was coupled into the Biosphere-Atmosphere Transfer Scheme (BATS), which is called BATS_CERES, to represent interactions between the land surface and crop growth processes. The effects of crop growth and development on land surface processes were then studied based on numerical simulations using the land surface models. Six sensitivity experiments by BATS show that the land surface fluxes had substantial changes when LAI was changed from 0 to 6 m2 m-2. Numerical experiments for Yucheng and Taoyuan stations reveal that the coupled model could capture not only the responses of crop growth and development to environmental conditions, but also the feedbacks to land surface processes. For quantitative evaluating the effects of crop growth and development on surface fluxes in China, two numerical experiments by BATS_CERES (CSM-run) and the original one BATS (CTL-run) conducted over the continental China. Comparison of the two runs shows decreases of leaf area index and fractional vegetation cover when incorporating dynamic crop in land surface simulation, which lead to less canopy interception, vegetation transpiration, total evapotranspiration, top soil moisture, and more soil evaporation, surface runoff, root zone soil moisture. These changes are accompanied by decreasing latent heat flux and increasing sensible heat flux in the cropland region. In addition, the comparison between the simulation and observation proved that incorporating the crop growth and development process into the land surface model could reduce the systematic biases of the simulated leaf area index and top soil moisture, hence improve the simulation of land surface fluxes.
2011, 28(4): 945-956. doi: 10.1007/s00376-010-0073-5
[Abstract](546) [PDF 808KB](1082)
Abstract:
Millimeter-wavelength radar has proved to be an effective instrument for cloud observation and research. In this study, 8-mm-wavelength cloud radar (MMCR) with Doppler and polarization capabilities was used to investigate cloud dynamics in China for the first time. Its design, system specifications, calibration, and application in measuring clouds associated with typhoon are discussed in this article. The cloud radar measurements of radar reflectivity (Z), Doppler velocity (Vr), velocity spectrum width (Sw) and the depolarization ratio (LDR) at vertical incidence were used to analyze the microphysical and dynamic processes of the cloud system and precipitation associated with Typhoon Nuri, which occurred in southern China in August 2008. The results show the reflectivity observed using MMCR to be consistent with the echo height and the melting-layer location data obtained by the nearby China S-band new-generation weather radar (SA), but the Ka-band MMCR provided more detailed structural information about clouds and weak precipitation data than did the SA radar. The variation of radar reflectivity and LDR in vertical structure reveals the transformation of particle phase from ice to water. The vertical velocity and velocity spectrum width of MMCR observations indicate an updraft and strong turbulence in the stratiform cloud layer. MMCR provides a valuable new technology for meteorological research in China.
2011, 28(4): 957-972. doi: 10.1007/s00376-009-9057-8
Abstract:
A moist thermodynamic advection parameter, defined as an absolute value of the dot product of horizontal gradients of three-dimensional potential temperature advection and general potential temperature, is introduced to diagnose frontal heavy rainfall events in the north of China. It is shown that the parameter is closely related to observed 6-h accumulative surface rainfall and simulated cloud hydrometeors. Since the parameter is capable of describing the typical vertical structural characteristics of dynamic, thermodynamic and water vapor fields above a strong precipitation region near the front surface, it may serve as a physical tracker to detect precipitable weather systems near to a front. A tendency equation of the parameter was derived in Cartesian coordinates and calculated with the simulation output data of a heavy rainfall event. Results revealed that the advection of the parameter by the three-dimensional velocity vector, the covariance of potential temperature advection by local change of the velocity vector and general potential temperature, and the interaction between potential temperature advection and the source or sink of general potential temperature, accounted for local change in the parameter. This indicated that the parameter was determined by a combination of dynamic processes and cloud microphysical processes.
2011, 28(4): 973-984. doi: 10.1007/s00376-010-9210-4
Abstract:
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 numerical 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.