Aerosol optical properties determined from sky-radiometer over Loess Plateau of Northwest China

Aerosol optical properties determined from sky-radiometer over Loess Plateau of Northwest China Y. Liu, J. Huang, G. Shi, T. Takamura, P. Khatri, J. Bi, J. Shi, T. Wang, X. Wang, and B. Zhang Key Laboratory for Semi-Arid Climate Change of the Ministry of Education and College of Atmospheric Sciences, Lanzhou University, Lanzhou, 730000, China LASG, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China Center for Environmental Remote Sensing, Chiba University, Chiba, Japan


Introduction
Human activities may create a lot of aerosol particles which are important variables in the Earth's energy budget.Aerosols influence the energy balance either by directly absorbing and scattering the shortwave and absorbing and emitting longwave radiation (direct effect) (Charlson et al., 1992;Li and Kou, 1998;Che et al., 2005) or by affecting the size distribution of cloud droplets (indirect effect) (Li, 1998;Ramanathan and Crutzen, 2003).Recently there are many studies on the radiative effect of aerosols, as described in the assessment report of the International Panel on Climate Change (IPCC, 2007), but the quantitative effect of aerosols on energy balance and climate Figures system is still largely uncertain (IPCC, 2007).Recently, special attention has been dedicated to cloud interactions with desert aerosol particles (Rosenfeld et al., 2001;Br éon et al., 2002;DeMott et al., 2003;Kawamoto and Nakajima, 2003;Huang et al., 2006a, b).However, the knowledge of the effects of Asian dust aerosols on arid climate is still very limited due to the lack of observations.Aerosols are generally believed to exert a cooling influence on the climate directly by scattering solar radiation and through their indirect effects on clouds.However, the semi-direct effect has the potential to offset this cooling by reducing low cloud cover and water path.Analysis of the satellite observations indicates that, on average, the water path of dusty clouds is considerably smaller than that from dust-free clouds in the same frontal systems (Huang et al., 2006a, b).The absorption or diabatic heating of Asian dust particles can cause the evaporation of cloud droplets and reduce the cloud water path.The key mechanism may be related to the dust aerosol warming effect through the absorption of solar radiation.
To understand the radiative effect of aerosols, it's critical to know the optical and physical properties and their temporal and spatial variability.Many efforts have been devoted to measure aerosol optical and physical properties (Holben et al., 2001;Takemura et al., 2002;Dubovik et al., 2002;Smirnov et al., 2002;Nakajima et al., 2003;Aoki and Fujiyoshi, 2003;Kim et al., 2004;Eck et al., 2005).Uchiyama et al. (2005) studied the characteristics of aeolian dust observed by sky-radiometer at several sites including some in source region.Li et al. (2007)  from Sun photometer and sky-radiometer at Beijing.Huang et al. (2008a) observed the vertical distribution of aerosols in the free troposphere using surface Micro-Pulse Lidar (MPL).Semi-arid region over Northwest of China has been widely affected by naturally occurring dust storms event especially during the period of early spring to summer.Such event causes not only the serious damage to agriculture, economics, environment and human health in this region but also a climatic impact on the global scale.Ge et al. (2010) retrieved dust aerosol optical properties from Multi-filter Rotating Shadowband Radiometer (MFRSR) and estimated the dust aerosol direct radiative forcing over Northwestern China during the 2008 China-US Joint Field Experiment.Bi et al. (2010) analyzed the aerosol optical properties measured at Semi-Arid Climate and Environment Observatory of Lanzhou University over Loess Plateau of Northwestern China obtained through a CIMEL Sun photometer of the ground-based Aerosol Robotic Network for the period of August 2006-October 2008.However, there are some difference between the results of Ge et al. (2010) and Bi et al. (2010).This study will use another measurement (sky-radiometer) to calculate and check the optical and radiative properties of the dust aerosol at SACOL over Loess Plateau of Northwest China.
2 Instrument and observational data SACOL, located approximately 48 km from the center of Lanzhou on the southern bank of the Yellow River in Gansu province northwest China (35 • 57 N, 104 • 08 E, 1965.8 m a.s.l.), is one of the reference sites of the international Coordinated Energy and Water Cycle Observations Project (CEOP).The parent soil material is mainly the Quaternary aeolian Loess with the main soil type of sierozem.The terrain where the measurements were carried out was flat and covered with short grass.Vegetation was usually less than 15 cm tall and covered less than 80 % of the surface in summer and autumn.SACOL is located within temperate continental climate zone, and the annual mean air temperature is 6.7 • .The annual precipitation at SACOL is about 381.8 mm, Introduction

Conclusions References
Tables Figures

Back Close
Full while the annual mean evaporation is about 1528.5 mm (Huang et al., 2008b).POM-01 sky-radiometer has been installed on 15 March 2009 at SACOL to measure the optical properties of aerosols in this region.It has been continuously running since then in 2009.A set of Kipp & Zonen CM21 pyranometer was also installed to measure the global solar irradiance every 10 s automatically at SACOL in March 2009, which is a high precision pyranometer with strictly selected domes.An automatic sun and sky scanning radiometer (CIMEL Sun photometer) was run at SACOL from August 2006, which will be useful to perform inter-comparison of aerosol optical parameters obtained from the two different instruments.The characteristics of the POM-01 sky-radiometer and CIMEL are shown in Table 1.On the calibration of these two instruments, skyradiometer can be calibrated by using not only traditional Langley method but also Improved Langley method.
The sky-radiometer at SACOL was calibrated for the sky radiance using an integrating sphere at Tsukuba Space Center and for the direct solar irradiance using the Improved Langley plot method.The precision of the in situ method has been estimated to be within 1-2.5 %, depending on the wavelength (Campanelli, 2004).The CIMEL Sun photometer measures the direct sun and diffuse sky radiances within the spectral range 340-1020 nm (Holben et al., 1998).The CIMEL Sun photometer located at SACOL is calibrated at the pristine Mauna Loa Observatory (MLO; latitude: 19 • 53 N, longitude: 155 • 57 W, elevation: 3400 m) once per year.
In view of the advantage on the calibration of sky-radiometer, this study will analyze the optical properties observed by sky-radiometer and calculate the radiative effect of the dust aerosols.But comparison of the aerosol optical properties from the skyradiometer and those from CIMEL will first be performed.
In this study, the latest version of Skyrad package (Version 4.2), which is developed by Nakajima et al. (1996), was used to analyze the data measured by sky-radiometer.
During the analysis of sky-radiometer, we need the total amount of ozone and the surface pressure to consider the ozone absorption and the Raleigh scattering of the air molecules.According to TOMS data and the meteorological observation data at Introduction

Conclusions References
Tables Figures

Back Close
Full SACOL, we set ozone as 0.3 cm and surface pressure as 1 atm in this study.The pyranometer data was used to remove cloudy data observed by sky-radiometer (Pradeep and Takamura, 2009).
The comparison between the aerosol optical depths observed by CIMEL Sun photometer and POM-01 sky-radiometer is first carried out, which is based on the 703 clear-sky measurements from 27 days in MAM.The detailed information of observational days is shown in Table 2.The data can represent the daily and monthly situation according to the measurement days in Table 2. Figure 1 shows the plots of AOD values derived from POM-01 and CIMEL at the two common wavelengths of 675 and 870 nm.Since there is no 500 nm channel for this CIMEL in the study, we compared the AOD values derived from POM-01 sky-radiometer and CIMEL at the two common wavelengths of 675 and 870 nm.High correlation coefficient of about 0.99 between the AOD from the POM-01 and CIMEL confirms the high consistency of the two instruments.The relative difference in the AOD between the POM-01 and CIMEL is about 4 %.Note that the observations shown in Fig. 1 exclude cloudy-sky measurements and the following analysis will only consider the clear-sky observation.

Optical depth
The daily mean AODs at 500 nm in March, April and May are shown in Fig. 2a, b, and  c, respectively.In March the averaged AOD is the largest (0.86) on 21 March and smallest (0.25) on 23 March 2009 (Fig. 2a).In April (Fig. 2b), the averaged AOD is the largest (0.58) on 8 April and smallest (0.23) is on 7 April 2009.The AOD is the largest (0.68) on 28 May and smallest (0.14) on 22 May (Fig. 2c).
By examining the global solar radiative flux observed by CM21 pyranometer and the weather record at SACOL along with the derived AOD, we found that 21 March and 8 Introduction

Conclusions References
Tables Figures

Back Close
Full according to the weather record, there are even more heavy dust events which were not caught by POM-01 but occurred on 23 April 2009, for example.As described by Uchiyama et al. (2005), if it is cloudy or very heavy dust storm, we cannot make the measurements.Thus the sky-radiometer might miss some very heavy dust events during the observational period in spring at SACOL.In the following, through the weather record and aerosol optical analysis from CIMEL and sky-radiometer observation, we define the observational days without obvious dust events as the clean day as background atmosphere.
Table 2 shows that the overall (including background and dusty days) monthly averaged AOD at 500 nm is 0.46, 0.42 and 0.26 in March, April and May, respectively.One can find that the monthly mean AOD is the largest in March and smallest in May, which agrees with the high frequency of dust event in March over Loess Plateau.The overall seasonal averaged AOD at 500 nm for spring is 0.40, which is comparable to the result derived from CIMEL (0.415 ± 0.23) at SACOL (Bi et al., 2010).
For the background aerosols at SACOL, the monthly averaged AODs are significant different (shown in Table 3).The averaged AOD of background aerosols (0.37) at 500 nm in MAM is much smaller than that during dust event (0.72). Figure 3 gives the frequency distributions of AOD at 500 nm observed by skyradiometer in MAM at SACOL.As also shown in Fig. 2, the AOD at 500 nm in MAM fluctuates with range from 0.1 to 0.9 during the measurement days.We find that most of the AOD values (about 60 %) are in the range from 0.1 to 0.4, in which the highest frequency of the AOD values (about 30 %) is in the range from 0.2 to 0.3.

Single scattering albedo and asymmetry factor
During the observation days measured by POM-01, the daily mean SSA at 500 nm derived from sky-radiometer ranging from 0.85 to 1.0, 0.9 to 1.0 and 0.8 to 1.0 in March, April and May, respectively.The monthly mean SSA retrieved directly is the largest (0.97) for April and smallest (0.94) in May at SACOL and the monthly mean ASY at SACOL ranges from 0.70 to 0.72 in MAM Introduction

Conclusions References
Tables Figures

Back Close
Full We will further examine the SSA and ASY using the radiative closure experiments in 4.1 where we found that there may be smaller than these retrieved values.

Comparison between the observed and simulated radiative fluxes
Before calculating the aerosol radiative forcing (ARF), we must determine the accuracy of the retrieved aerosol optical properties from POM-01 sky-radiometer.Using SBDART model, we simulated the solar radiation fluxes (total, direct and diffuse) at the surface using the retrieved aerosol optical properties including AOD, SSA and ASY.In the simulation, ozone concentrations are taken from TOMS, precipitable water content (PWC) are taken from microwave radiometer (MWR), and the surface albedo is derived from the calculation based on the CM21 observation at SACOL (the mean value is 0.25 in this study).
Comparing the calculated broadband radiative fluxes with these observed by pyranometer and pyrheliometer during the same days at SACOL, the relative difference in the direct, total and diffuse radiation are 1.8 %, 3.13 % and 12.16 %, respectively.The good agreement between the calculated and observed direct radiations indicates that the retrieved AOD are accurate, which was also shown in the comparison with the AOD by CIMEL (Fig. 1).The disagreement between the diffuse fluxes indicates that the SSA derived from sky-radiometer may have significant errors.To obtain the accurate ARF, we must adjust SSA to match the observed radiative fluxes.Figur 4 shows an example of radiative closure experiment.In Fig. 4, result 1 is the calculated result using the retrieved aerosol optical properties directly including AOD, SSA and ASY; result 2 is the simulation after adjusting the SSA at 500 nm from 0.95 to 0.8.We found that the values of SSA observed by POM-01 are larger than the real values at SACOL, i.e., the SSA retrieved from sky-radiometer maybe underestimate the absorption of the aerosols over the Loess Plateau in Northwest of the China.Introduction

Conclusions References
Tables Figures

Back Close
Full The daily mean values of adjusted SSA at 500 nm are shown in Fig. 5 ranging from 0.82 to 0.97 at 500 nm in MAM at SACOL.The adjusted mean SSA at 500 nm is 0.92 in MAM, while the unadjusted (shown in bracket in Table 2) mean SSA is 0.95.The relative difference in the mean SSA between adjusted and unadjusted value is about 3.16 % in MAM.Compared with the SSA observed at the other observatory over the Loess Plateau (Ge et al., 2010), the values of the SSA (adjusted) are larger in this study.The adjusted ASY at 500 nm ranges from 0.65 to 0.85 in MAM (not shown in figure ).
As shown in Table 3, the background monthly mean SSA is the larger (0.93) in March and April and smaller (0.90) in May at SACOL.The averaged SSA of background is 0.92 in MAM at SACOL, while the mean value during dust events is 0.89 which agrees with the 0.89 of Dunhuang (Kim et al., 2004).The smaller SSA can be interpreted as the result of larger particles during dust events (Fu et al., 2009).The averaged ASY during dust period is much larger than that of background aerosols, indicating that the lager dust particles that lead to stronger forward scattering (Fu et al., 2009).
Figure 6 shows the comparison of the surface solar irradiances observed by CM21 pyranometer and those calculated by radiative transfer model using adjusted SSAs and ASYs, indicating an excellent agreement.
Based on the calculated solar radiation at the TOA and surface by using the adjusted SSAs and ASYs, we calculate the aerosol shortwave radiative heating rate at SACOL (not shown).The aerosols during dust event period have a larger impact on the radiative energy budget than background aerosols over Loess Plateau.The maximum 24-h averaged radiative heating rate reaches 2.01 K day −1 on 21 March which was a strong dusty day.The aerosols heat the atmosphere (daily averaged) by up to 2 K day −1 due to dust event, which is 67 % larger than the heating (1.2 K day −1 ) of background aerosols.Introduction

Conclusions References
Tables Figures

Back Close
Full

Radiative forcing due to dust aerosol
The ARF at the surface is estimated as follows, where R is the surface reflectance, and F abs (t) is the observed solar irradiance at the surface (SFC) and time t.F SRF,Ray (t) is the downward solar irradiance at the surface and time t for aerosol-free atmosphere, which is estimated by a radiative transfer calculation using the SBDART model.The integration term in Eq. ( 2) shows the radiative effect of cloud and/or aerosol.Using Eqs. ( 1) and ( 2 ) can be calculated as follows, where F N SFC and F N SFC,Ray are the daily mean net solar radiation at the surface with and without aerosols, respectively.
Figure 7 shows the comparison of the ARF calculated by the radiative transfer model and that derived from the CM21 observation.ARF where F N TOA and F N TOA,Ray are daily mean net solar radiative flux at the TOA with and without aerosols, respectively.Figure 8 shows the daily mean radiative forcing at the surface (green), TOA (red)  .In general, the ARF in the atmosphere, which represents absorbed solar radiation within the atmosphere, heats the atmosphere.In this study, aerosols have a heating effect in the atmosphere especially during the dust event period, which range from 12.47 to 100.24 wm −2 .
The monthly averaged values of ARF are shown in Fig. 9.The monthly mean values of ARF at the surface and TOA in MAM at SACOL are negative forcing, while that in the atmosphere is positive.The overall seasonal mean ARF is −16.79,−53.51 and 36.72 wm −2 , respectively, at the TOA, surface and in the atmosphere in MAM at SACOL.As shown in Table 3, the heating effect in the atmosphere (positive ARF, 80.08 wm −2 ) of the aerosols during dust event period is much more significant than background aerosol (ARF in the atmosphere is 33.26 wm

Conclusions and summary
In this study, we investigated the aerosol optical properties and their associated radiative forcing retrieved from the POM-01 sky-radiometer and surface solar radiation data collected at the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL) in MAM 2009.The aerosol radiative forcing was simulated by SBDART model.The aerosol optical property was retrieved from the sky-radiometer.SSA was obtained by matching the simulation using radiative transfer model with observation.
The averaged SSA of background aerosols is 0.92 in MAM, while the mean value due to dust events is 0.89.The maximum 24-h averaged radiative heating rate reaches 2.01 K day −1 on 21 March which was a strong dusty day.The aerosols heat the atmosphere (daily averaged) by up to 2 K day −1 due to dust event, which is 67 % larger than the heating (1.2 K day −1 ) of background aerosols.The heating effect in the atmosphere (positive ARF, 80.08 wm −2 ) of the aerosols during dust event period is much more significant than background aerosol (ARF in the atmosphere is 33.26 wm −2 ) in MAM 2009.Dust aerosols may heat the atmosphere strongly at SACOL.Introduction

Conclusions References
Tables Figures

Back Close
Full  Climate, 15, 333-352, 2002. Uchiyama, A., Yamazaki, A., Togawa, H., and  Full  Full  Full Discussion Paper | Discussion Paper | Discussion Paper | examined continuous high-quality measurements of aerosol at Xianghe station in China during the period of September 2004 to September 2005.Retrieving from sky-radiometer measured at Anmyeon, Korea, Sohn et al. (2007) pointed out that Asian dusts can become more absorbing aerosols when mixed with pollutants.Hayasaka et al. (2007) observed the vertical distributions and optical properties of nonspherical dust and spherical aerosol particles from March to May 2005 by Lidar and sky-radiometer systems at Sapporo, Toyama, and Nagasaki, Japan.Xia et al. (2007) studied aerosol optical properties by a Sun photometer at Taihu, China.Che et al. (2008) compared the optical properties derived Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ), we can calculated the ARF (ARF CM21 SRF ) based on the global solar irradiance observed by CM21 at SACOL We can regard ARF CM21 SRF as a reference to the simulated ARF by SBDART because of the high precision of the CM21 pyranometer.Using the SBDAT model, ARF at the surface (ARF SBDART SRF

TOA
correlation with a mean bias of about 3.62 wm −2 .ARF at the top of atmosphere and in the atmosphere (ATM) could be calculated using the radiative transfer model as follows ARF TOA = F N TOA − F N Discussion Paper | Discussion Paper | Discussion Paper | ARF ATM = RF TOA − RF SFC (4) −2 ) in MAM 2009, which is determined by the larger optical depth and small SSA of the dust aerosols during dust event over Loess Plateau in Northwest of the China.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | and radiative forcing of various aerosol species with a global three-dimensional model, J.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |
and in the atmosphere (blue) in MAM 2009.Aerosols at SACOL have a strong negative radiative forcing at the surface, ranging from −22.92 to −100.41 wm −2 .The negative forcing is the largest (−100.41wm −2 ) on 21 March 2009, which corresponds with the largest value of AOD in MAM 2009.The surface ARF at other site over loess plateau of Northwest China (Ge et al., 2010) ranges from −7.9 to −35.8 wm −2 from 25 April to 25 May, while the averaged ARF at the surface in our study ranges from −22.92 to −72.26wm −2 from April to May.The critical reason of such difference is the larger optical depth in our study.Daily averaged values of radiative forcing at the surface in May are smaller than that in March and April.Most of the ARF values at the TOA are negative, and the value of ARF fluctuates from 4.35 to −42.8wm −2 Asano, J.: Characteristics of Aeolian dust observed by sky radiometer in the ADEC Intensive Observation Period I (IOP1), J.Meteorol.Soc.Jpn., 83A, 291-305, 2005.

Table 2 .
The number of clear-sky days and measurements by sky-radiometer in spring of 2009 and the retrieved averaged AOD and adjusted SSA and ASY in MAM (unadjusted values are shown in brackets).

Table 3 .
Averaged AOD, SSA, ASY and ARF of background and dusty days in MAM (AOD is retrieved from sky-radiometer directly; SSA and ASY are adjusted values; ARF is calculated by radiative transfer model with retrieved AOD and adjusted SSA and ASY).