• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Seasonal Inhomogeneity of Soot Particles over the Central Indo-Gangetic Plains,India:Influence of Meteorology

    2015-12-12 12:56:54SINGHTIWARIHOPKEPhilipSINGHBISHTSRIVASTAVASINGHDUMKASINGHRAIandSRIVASTAVAManoj
    Journal of Meteorological Research 2015年6期

    SINGH B.P.,TIWARI S.,HOPKE Philip K.,SINGH R.S.,BISHT D.S.,SRIVASTAVA A.K., SINGH R.K.,DUMKA U.C.,SINGH A.K.,RAI B.N.,and SRIVASTAVA Manoj K.?

    1 Department of Geophysics,Banaras Hindu University,Varanasi 221005,India

    2 Indian Institute of Tropical Meteorology(New Delhi Branch),Delhi 110060,India

    3 Clarkson University,Box 5708,Potsdam,NY 13699-5708,USA

    4 Department of Chemical Engineering,Banaras Hindu University,Varanasi 221005,India

    5 Aryabhatta Research Institute of Observational Sciences,Manora Peak,Nainital 263002,India

    6 Department of Physics,Banaras Hindu University,Varanasi 221005,India

    Seasonal Inhomogeneity of Soot Particles over the Central Indo-Gangetic Plains,India:Influence of Meteorology

    SINGH B.P.1,TIWARI S.2,HOPKE Philip K.3,SINGH R.S.4,BISHT D.S.2,SRIVASTAVA A.K.2, SINGH R.K.1,DUMKA U.C.5,SINGH A.K.6,RAI B.N.4,and SRIVASTAVA Manoj K.1?

    1 Department of Geophysics,Banaras Hindu University,Varanasi 221005,India

    2 Indian Institute of Tropical Meteorology(New Delhi Branch),Delhi 110060,India

    3 Clarkson University,Box 5708,Potsdam,NY 13699-5708,USA

    4 Department of Chemical Engineering,Banaras Hindu University,Varanasi 221005,India

    5 Aryabhatta Research Institute of Observational Sciences,Manora Peak,Nainital 263002,India

    6 Department of Physics,Banaras Hindu University,Varanasi 221005,India

    Black carbon(BC)particles play a unique and important role in earth's climate system.BC was measured (in-situ)in the central part of the Indo-Gangetic Plains(IGP)at Varanasi,which is a highly populated and polluted region due to its topography and extensive emission sources.The annual mean BC mass concentration was 8.92±7.0μg m?3,with 34%of samples exceeding the average value.Seasonally,BC was highest during the post-monsoon and winter periods(approximately 18μg m?3)and lower in the premonsoon/monsoon seasons(approximately 6μg m?3).The highest frequency(approximately 46%)observed for BC mass was in the interval from 5 to 10μg m?3.However,during the post-monsoon season,the most common values(approximately 23%)were between 20 and 25μg m?3.The nighttime concentrations of BC were approximately twice as much as the daytime values because of lower boundary layer heights at nighttime. The?Angstr¨om exponent was significantly positively correlated(0.55)with ground-level BC concentrations, indicating the impact of BC on the columnar aerosol properties.The estimated mean absorption?Angstr¨om exponent was 1.02±0.08μg m?3,indicating that the major source of BC was from fossil fuel combustion. Significant negative correlations between BC mass and meteorological parameters indicate a pronounced effect of atmospheric dynamics on the BC mass in this region.The highest mean BC mass concentration (18.1±6.9μg m?3)as a function of wind speed was under calm wind conditions(38%of the time).

    black carbon,Indo-Gangetic Plains,absorption?Angstr¨om exponent,biomass,fossil fuel

    1.Introduction

    Black carbon(BC)is produced by incomplete combustion of natural as well as anthropogenic substances,including the burning of biofuels and fossil fuels(Bond and Bergstrom,2006;Hyv¨arinen et al., 2010;Praveen et al.,2012;Srivastava et al.,2012a; Safai et al.,2013;Tiwari et al.,2013).The presence of BC in the atmosphere causes absorption of solar radiation in the visible and near infrared wavelengths.It exerts positive radiative forcing(i.e.,heating)at the top of the atmosphere,negative forcing(i.e.,cooling) at the surface,and heating in the lower troposphere. Because of these effects,it is an important constituent of the atmosphere that can cause global and regional warming(Satheesh and Ramanathan,2000;Jacobson,2001;IPCC,2007;Ramanathan and Carmichael, 2008;Moosm¨uller et al.,2009;Srivastava et al.,2012b;

    Bond et al.,2013).According to the recent IPCC report(Stocker et al.,2014),the contribution of BC to atmospheric radiative forcing(direct impact)is 0.4 (0.05–0.8)W m?2.However,Bond et al.(2013)reported in an extensive review that the BC radiative forcing is 1.1(0.17–2.1)W m?2. It is the second strongest contributor to current global warming,after carbon dioxide(Ramanathan and Carmichael,2008; Bond et al.,2013).

    Apart from its direct impact on the atmosphere, indirect impacts of BC are also apparent;for example, its influence on cloud droplet number concentrations and related alterations of cloud properties(Kaufman et al.,2002).The presence of higher concentrations of BC in the atmosphere makes it an important ingredient in the analysis of regional climate(Wang, 2004;Liao and Shang,2015),and instrumental in the alteration of the heat regime of regional air as well as regional atmospheric stability/vertical motion.It can also cause alterations in regional atmospheric hydrodynamics(Menon et al.,2002;Ramanathan and Carmichael,2008;Bond et al.,2013;Wang et al., 2014).

    It has been hypothesized that BC may also impact upon the precipitation patterns over the Indian subcontinent.The precipitation over the past several decades,with increased rainfall in South Asia and increased drought in North Asia,is attributed to regional,post-industrial era emissions(Menon et al.,2002;Hyv¨arinen et al.,2010).Specifically for the Indian climate,increased particulate loadings,which can absorb radiation,are believed to have caused the advancement of rainy periods as well as intensification of the Indian monsoon system(Lau and Kim, 2006).Existence of BC in the atmosphere is,however, a complicated issue,because it can be generated locally as well as transported regionally.Although BC exhibits a short atmospheric lifetime of around one week to 10 days,compared with hundreds of years for carbon dioxide(Jacobson,2002),the continuous impingement of BC due to anthropogenic emissions is of concern for climate change estimations.Furthermore,the BC mass concentration has been increasing in the recent past because of the changing socioeconomic status and rapid industrialization in many parts of the globe,especially in developing countries.High temporal resolution and real-time BC data are useful to estimate the impact of BC on climate,environment,and human health.However,BC measurements at most global locations,including the Indo-Gangetic Plains(IGP),are limited.The sources and mechanisms that affect their concentrations in the air are also mostly unknown.The scarcity of real-time BC data over the IGP is especially important since it is one of the most populated and polluted regions of the world.The common finding of studies carried out in the IGP region is one of high variability in BC concentrations(Tripathi et al.,2005;Safai et al.,2008; Beegum et al.,2009;Tiwari et al.,2009;Hyv¨arinen et al.,2010;Ramachandran and Kedia,2010;Bano et al., 2011;Moorthy and Satheesh,2011;Raju et al.,2011; Praveen et al.,2012).These findings suggest the need for continuous long-term monitoring of BC at multiple locations around the IGP,so as to estimate the impacts of BC effectively and assist with mitigation plans.To the best of our knowledge,this paper is the first report of long-term,in-situ measurements of BC over Varanasi(25.3°N,83.0°E;76 m above mean sea level),which is located in the central part of the IGP region.The major objectives of the present study were to:(1)assess the variability of BC concentrations on different timescales;(2)understand the impact of BC concentrations on columnar aerosol optical properties; and(3)examine the effect of surface meteorology on the distribution of BC over this study region.

    2.Site description and experimental setup

    2.1Sampling location and meteorological conditions

    In-situ measurements of BC mass concentrations were made on the campus of Banaras Hindu University,Varanasi,during the period 1 January to 31 December 2009.The data were gathered as part of a national program called the“Aerosol Radiative Forcing Initiative Network”funded by the Indian Space Research Organization,over the IGP region.Varanasi (25.28°N,82.95°E)is located in the central Ganges

    Valley of North India,in the eastern part of the state of Uttar Pradesh,along the left crescent-shaped bank of the Ganges River.The population of Varanasi City was approximately 1.2 million in 2011.The land surrounding the city is very fertile and supports substantial agricultural activity.The city also possesses many small-scale industries.Air quality in Varanasi is very poor due to emissions from badly maintained automobiles and heavily loaded transport vehicles,as well as from small-scale industries,domestic heating,and large-scale construction activities(Kumar et al.,2015; Tiwari et al.,2015a,b).

    There are four major seasons:winter(December to next February),the pre-monsoon season(March to June),the monsoon season(July to September), and the post-monsoon season(October and November)(Sharma et al.,2010;Mishra and Shibata,2012). During winter,the wind generally moves eastward with the passage of extratropical low pressure systems (called westerly disturbance),bringing occasional rain or fog over the IGP and snow over the Himalayan region.During the pre-monsoon season,the temperature rises rapidly(up to 49℃)from March to June, when the winds blow from the west.The region receives heavy rainfall during the monsoon season,when depressions and moderate low pressure systems form over the northern Bay of Bengal and move northwestward.

    Meteorological variables,such as temperature (TEM),relative humidity(RH),wind direction,wind speed(WS),and rainfall were obtained from an automatic weather station installed in an open field at a distance of about 250 m from the measurement site.Mixing height(MH)and visibility(VIS)data were obtained from the National Oceanic and Atmospheric Administration(NOAA)Air Resources Laboratory Hybrid Single-Particle Lagrangian Integrated Trajectory(HYSPLIT)model using NCEP FNL(final)analysis data.The daily variations of the abovementioned meteorological variables are shown in Fig. 1.During the study period,the mean temperature was approximately 27±2.1℃,varying from 17.6℃(January)to 35.8℃(June).WS varied from 1.72(December)to 3.61 m s?1(May),with a mean of about 2.5±0.8 m s?1.In addition,the means(±standard deviation)of the MH,RH,rainfall,and VIS were 750 (±172.5)m,60(±10)%,864.3(±198.7)mm,and 4.1(±0.6)km,respectively.

    Fig.1.Daily variations of meteorological parameters over Varanasi.(a)Temperature(TEM),(b)relative humidity (RH),(c)visibility(VIS),and(d)wind speed(WS).

    2.2BC monitoring instrument

    BC mass concentrations were measured with an Aethalometer(Model:AE-31;Magee Scientific Co., USA,http://www.mageesci.com)at a flow rate of 4.0 L min?1.This instrument measures the optical absorption at seven wavelengths(370,470,520,590,660, 880,and 950 nm)during a discrete time interval(2 min).The BC mass concentration is calculated assuming that the attenuation is linearly proportional to the amount of BC on the filter,and calculated in accordance with earlier studies(Ramachandran and Rajesh,2007;Tiwari et al.,2009,2013;Srivastava et al.,2012b,2014;Safai et al.,2013)to convert the filter transmittance in inverse mega meters(Mm?1;1 Mm?1=10?6m?1)to BC concentrations.The min-

    imum detection limit of the Aethalometer,defined as twice the standard deviation of the noise,is 5 ng m?3(Virkkula et al.,2007).The attenuation coefficientis calculated by multiplying the measured BC mass concentration(μg m?3)at 880 nm and constant cross-sectional value of attenuation(i.e.,

    The attenuation coefficient obtained from Eq.(1)is corrected for the particle loading effect,also known as the shadowing effect(R),and the multiple scattering effect(C),in order to obtain the aerosol absorption coefficient(babs;Mm?1):

    As a filter-based measurement,the BC measured by the Aethalometer may possess some artifacts due to the loading/shadowing effect,matrix effect and scattering effect(Weingartner et al.,2003; Arnott et al.,2005;Virkkula et al.,2007;Collaud et al.,2010).Weingartner et al.(2003)found that the loading/shadowing effect is more prominent for freshly emitted BC,while it is almost negligible for aged particles.The uncertainties in the estimates of BC mass concentration were around 7%.In another study,Dumka et al.(2010b)reported uncertainties of approximately 5%over the high altitude station at Manora Peak in the central Himalaya.In addition, literature also exists on the uncertainties of BC measurements by the Aethalometer(Weingartner et al., 2003;Sheridan et al.,2005;Corrigan et al.,2006). To obtain information on the sources of the BC particles,the wavelength exponent called the absorption ?Angstr¨om exponent(AAE,α)was estimated from the absorption at the seven wavelengths using the power law relationship(Kirchstetter et al.,2004)as previously used by Dumka et al.(2010b,2013):

    where K is a constant and λ is the wavelength(inμm). The wavelength exponent(α)is estimated by the linear regression of lnβabsand lnλ,which is a measure of the spectral dependence of aerosol absorption.An α value of approximately 1 denotes that the major BC source is from fossil fuel combustion.For biomass burning,it is between 1.5 and 3,and for dust it may be 2 to 3(Bergstrom et al.,2004;Kirchstetter et al., 2004;Sandradewi et al.,2008a,b;Russell et al.,2010).

    3.Results and discussion

    3.1BC aerosol characteristics over Varanasi

    Day-to-day and monthly variability in BC mass concentrations over Varanasi during the entire study period are shown in Fig.2.The daily mean hourly average BC concentration is 8.9±7.0μg m?3,varying from 34.6(7 January 2009)to 2.1(2 July 2009) μg m?3during January–December 2009.BC was relatively higher during the winter months and lower values were found in monsoonal months.The BC concentrations exceeded the average value on approximately 34%of the days.Monthly averages were highest in January(24.7μg m?3),followed by:November(24.6 μg m?3),December(20.4μg m?3),October(11.6μg m?3),February(10.0μg m?3),September(9.2μg m?3),March(8.6μg m?3),June(6.5μg m?3),July (4.8μg m?3),April(4.6μg m?3),August(4.5μg m?3),and May(4.2μg m?3).

    Fig.2.Daily,monthly,running mean,and annual mean BC mass concentrations over Varanasi during January–December 2009.

    Over the past decade,BC monitoring studies have been conducted in urban,rural,coastal,marine,and high altitude environments over the Indian subcontinent(Table 1)(Babu and Moorthy,2002;Tripathi et al.,2005;Pant et al.,2006;Safai et al.,2007,2008,

    2013;Awasthy et al.,2010;Dumka et al.,2010a;Pani and Verma,2010;Ramachandran and Kedia,2010; Tiwari et al.,2013;Bisht et al.,2015;Singh et al., 2015)and other global locations(Table 2),representing different environments and reflecting the variation in magnitude of BC concentrations from different regions.In India,significant spatial variability in BC was reported for higher concentrations(>10μg m?3)over 33 locations in the IGP region(Moorthy and Satheesh,2011).Lower BC concentrations(approximately 5μg m?3)were found in southern India (Hyderabad,Bangalore,Pune,Trivandrum,etc.).In mountainous regions(Nainital,Mukteshwar,and Sinhagad),concentrations have been found to be around 2μg m?3(Pant et al.,2006;Hyv¨arinen et al.,2009; Dumka et al.,2010b;Raju et al.,2011),except Kullu, where the BC mass concentration was approximately 4.6μg m?3(Kuniyal,2010).The much higher BC mass concentration at Kullu is likely due to vehicular emissions because of the area's popular tourist spots.Meanwhile,BC mass concentrations have been observed to be less than 2μg m?3at several locations, such as Godavari(Nepal),Uto(Finland),Granada (Spain),′Evora(Portugal),and Millan(Italy)–levels that are around six times lower than in Varanasi.However,in Lahore(Pakistan),BC has been measured to be approximately 2.5 times higher than the values reported in the present study(Husain et al.,2007).

    Seasonally,BC was higher during the postmonsoon and winter seasons(approximately 18μg m?3)and lower in the pre-monsoon/monsoon seasons (approximately 6μg m?3).Srivastava et al.(2012a) reported that open burning of crop residues is common practice in northern India during the pre-monsoon season,resulting in the transport of BC across the region. During most of the winter,the wind is calm,with lower temperatures,thick fog,and a low MH(Tiwari et al., 2013).Under such conditions,pollutants are not well dispersed,resulting in poor VIS and high local pollutant concentrations(Mohan and Bhati,2009).The frequency distribution of BC was classified into eight different concentration bins(in intervals of 5μg m?3from<5μg m?3to>40μg m?3)during the different seasons(Fig.3).The distribution was positively skewed,except for the bin ranging from 5 to 10μg m?3,which contributed approximately 46%of the observed concentrations.During the pre-monsoon and monsoon seasons,about 94%and 90%of the levels of BC mass were<10μg m?3,whereas only 30% and 15%were<10μg m?3in the post-monsoon and winter seasons,respectively.During the winter period, which often features deep smog/foggy conditions over the IGP,BC concentrations were found to be largely (26%)within 15–20μg m?3,while only 6%of the con-centrations was within 30–35μg m?3.A very high contribution of BC mass(39%)in the winter was observed within 20–30μg m?3.Similar characteristics were apparent during the post-monsoon season,withthe highest frequency(23%)within 20–25μg m?3. Fewer BC values(1%)were within 35–40μg m?3.

    Table 1.The mean mass concentration of BC over Varanasi in the present study compared with other earlier reported concentrations elsewhere in India

    Table 2.The mean mass concentration of BC over Varanasi in the present study compared with other earlier reported concentrations from various locations worldwide

    Fig.3.Seasonal and annual frequency distribution of BC mass concentrations over Varanasi.

    3.2 Diurnal variability of BC mass over Varanasi

    The BC concentrations in Varanasi were sorted into four different seasons above mentioned.Their seasonal mean diurnal variations are shown in Fig.4, along with their standard deviations.Different diurnal patterns were apparent during the different seasons.The variations were more pronounced during the winter and post-monsoon seasons compared with the pre-monsoon and monsoon seasons.BC showed a gradual rise at around 0700 LT(local time),peaking at around 0800–0900 LT.Low values were observed in the afternoon(1200–1900 LT),and there was a secondary maximum between 2100 and midnight(0100 LT).This large diurnal variability in BC mass was likely due to the combined effect of the diurnal variation in emissions and surface meteorology,including associated boundary layer dynamics.BC concentrations(maximum)during the pre-monsoon,monsoon,post-monsoon,and winter seasons in the morning hours(0800–0900 LT)were approximately 9,8, 25,and 21μg m?3,respectively,more than two to three times higher than the afternoon hours(1500–1700 LT).For the post-monsoon season,the rise in BC started at 0500 LT and continued until 0800 LT (24.1μg m?3),followed by a sudden drop in concentration until there incremental rises in values from 2000 LT until the next morning.In the morning,the surface inversion coupled to layers above a few hours after sunrise,resulting in the vertical mixing of ground-level pollutants with the free troposphere and causing the surface BC concentrations to suddenly decrease.Increased wind speeds after sunrise also increased mixing,thereby facilitating dispersion and dilution.After 0800 LT,BC emissions-primarily from vehicular traffic-caused increasing BC concentrations.During winter mornings,BC concentrations rose after 0800 LT(21.4μg m?3)and continued until 1000 LT.After that time,they dropped until 1700 LT.The mean BC mass concentrations at night(1800–0600 LT)were approximately twice the daytime values(0700–1700 LT).

    At night,the highest mean BC concentration was during the post-monsoon season(26.4μg m?3),followed by the winter(24.3μg m?3),monsoon(8.6μg m?3) and pre-monsoon(8.5μg m?3)seasons.The corresponding night/day ratios were 2.4,1.8,2.4,and 1.3, respectively.The BC concentration reached its minimum between 1400 and 1600 LT in every season due to the increased MH.After 1800 LT,the evening rush hour commences,as well as an increase in cooking activity and domestic heating use.Also,the surface inversion begins to form trapping pollutants.Subsequently,the BC mass concentration starts increasing from the evening onwards.The earlier onset of the morning rise and the later onset of the evening accumulation during the pre-monsoon season,compared with winter,correspond to the earlier sunrises and later sunsets in the pre-monsoon season.

    Fig.4.Seasonal diurnal variability of BC mass concentrations over Varanasi.

    Simpson and McGee(1996)described the marked effect of the local climate on the diurnal variation of pollutants due to fumigation effects,which greatly increase daily averages,particularly during winter.The diurnal BC variations in Varanasi have been observed at a variety of other locations,including a suburban site in Maryland,USA(Chen et al.,2001),the tropical coastal station of Trivandrum(Babu and Moorthy, 2002),and other inland sites in India(Sreekanth et al., 2007;Tiwari et al.,2013;Bisht et al.,2015).

    3.3 BC mass concentration in relation to aerosol optical properties

    Being highly light absorbent,BC particles affect various optical properties,such as aerosol optical depth(AOD),the?Angstr¨om exponent(AE),and single scattering albedo,which are responsible for aerosol radiative forcing(Pandithurai et al.,2008).AE is a good indicator of aerosol particle size and largely depends on the aerosol size distribution.It measures the comparative contributions of coarse-and finemode particles,where higher values of AE represent a higher contribution of fine-mode particles and lower values represent a higher contribution of coarse-mode particles(Pandithurai et al.,2008;Srivastava et al., 2011a;Singh et al.,2014).

    Fig.5.Seasonal values of(a)BC and(b)aerosol optical depth(AOD 550 nm)and?Angstr¨om exponent(AE)during January–December 2009.

    During the study period,AOD and AE were retrieved from MODIS(http://gdata1.sci.gsfc.nasa.gov/ daac-bin/G3/gui.cgi?instance?id=MODIS?DAILY?L3).The mean AOD(550 nm)was 0.65±0.3,varying from 0.39±0.2(February)to 1.01±0.4(August). The mean AE was 0.83±0.3,varying from 1.19±0.4 (January)to 0.56±0.1(July).Its seasonal variability

    is depicted in Fig.5.Prasad et al.(2005)reported an AOD over Varanasi during the monsoon period that was similar(>0.6)to the present study.The monthly AOD increased from April(0.43±0.2)to June(1.01 ±0.4),and decreased in July(0.74±0.3),followed by a sudden increase in August(0.98±0.2).The increase may be due to a long dry period/break phase of the monsoon during this month.AOD increased further in October(0.44±0.2)and continued increasing until December(0.77±0.3).

    AOD represents the attenuation of direct solar irradiance.It may be higher in situations where either coarse-or fine-mode particles,or both,increase.Along with AE,delineation of a higher contribution of fine-or coarse-mode particles may be deduced.The AOD values were negatively correlated(–0.27)with AE and,in general,the higher AOD values were associated with lower AE during the pre-monsoon season.This situation suggests that the dominance of coarse-mode particles was likely due to transported dust(Singh et al.,2005;Srivastava et al.,2011a).Conversely,high AODs were associated with high AE values during the winter and post-monsoon periods,indicating the dominance of fine-mode particles.Fine particles are mostly produced by enhanced anthropogenic activities(Srivastava et al.,2011b).Regression analysis between the AE and BC daily data revealed a significant positive correlation(0.55),indicating that fine particles, including BC,dominated during the study period.

    Reddy and Venkataraman(2002)reported that fuel wood and crop wastes were the primary contributors to biomass-based BC emissions in northern India.These estimates were,however,based on annual average emissions.It has been shown that,during November and next March,anthropogenic source contributions exceed 70%of the measured AODs in and around India(Ramachandran,2004).This high anthropogenic influence on AOD is due to the man-made submicron aerosols from local and regional sources. Regional transport is greater during the winter season.AE showed a decreasing trend from January(1.24 ±0.30)to May(0.55±0.10),and then increased from June onwards,up to a value of 1.20±0.40 in December.

    Seasonally averaged AE(Fig.5)was observed to be higher during winter(1.16±0.20),followed by the post-monsoon(0.97±0.20),monsoon(0.76±0.10), and pre-monsoon(0.58±0.10)seasons.Higher BC and AE during the post-monsoon and winter seasons indicate the dominance of fine particles in this location.A long-term analyses in Delhi by Lodhi et al. (2013),and episodic observations by Pandithurai et al. (2008),showed the lowest AE values in June(approximately 0.4)and highest values(approximately 1.0) in winter,gradually decreasing from winter months to a peak during the pre-monsoon months because of the increasing importance of soil-derived coarse-mode particles(dust transport from the nearby Thar desert) (Singh et al.,2005)along with the possible mixing of boundary layer particles with free tropospheric material caused by surface heating during the pre-monsoon period.The observation in Varanasi confirms the dominance of fine-mode particles,largely BC,during the winter and post-monsoon seasons.

    3.4 Spectral aerosol absorption characteristics

    BC particles are largely produced by the burning of fossil fuels from transport,industrial activity,the power sector,etc.,and the burning of biomass/biofuels (home heating/cooking,agricultural burning,and wild fires).As noted above,the AAE can help to identify the possible sources of BC.The AAE values were estimated from BC measurements at seven discrete wavelengths,and its variability is depicted in Fig.6.The mean value of AAE in Varanasi was 1.02±0.08,varying from 0.73 to 1.32 over the observation period,indicating a strong influence of fossil fuel combustion. Seasonal AAE was generally higher during the postmonsoon and winter seasons,and lower in the premonsoon/monsoon season,following the order:postmonsoon season(1.11±0.10),winter(1.02±0.10), monsoon season(1.01±0.10),and pre-monsoon season(0.98±0.10).These results suggest fossil fuel combustion as the major source of BC(Soni et al., 2010)during the pre-monsoon and winter periods in Varanasi. In Delhi,a similar AAE value(1.03±0.09)was observed with relatively similar seasonal values(Ganguly et al.,2005).However,in central India,

    Ganguly et al.(2005)reported the AAE value as 1.52, suggesting biofuel/fossil fuel emissions.Bond(2001) reported values of AAE between 1 and 2.9 and suggested that the BC source is residential biofuel burning.Aruna et al.(2013)reported a value of AAE around 1.1(varying from 0.9 to 1.1)for a tropical location near Chennai.The diurnal variability of AAE for different seasons was examined to separate the possible BC sources in Varanasi(Fig.7).Two major peaks were observed:morning(0800–1100 LT)and evening (2000–2300 LT),confirming the combined effect as the burning of fossil fuels and biomass during the morning and evening period.In addition,as a rural environment is very close to the city,the people burn anything (e.g.,wood,cow dung,leaves,crop residue,etc.)for heating purposes and for cooking food(during morning and evening),which produces large amounts of BC.

    Fig.6.Daily variation of absorption?Angstr¨om exponent, along with the running and annual mean,over Varanasi.

    Fig.7.Diurnal variations of absorption?Angstr¨om exponent over Varanai during January–December 2009.

    3.5Influence of meteorology on BC

    The impact of changes in surface meteorological conditions on the BC concentrations was investigated.Regression analysis between daily BC concentrations and surface meteorological parameters in Varanasi during the study period is presented in Fig.8.WS yielded a negative correlation coefficient (–0.50),indicating that increasing WS increased the dilution of local emissions.WS was found to be highest during the pre-monsoon season(11.1 km h?1)and lowest during winter(7.0 km h?1),and the corresponding mean BC values were 6.64 and 18.61μg m?3, respectively.A similar relationship was observed in Ahmadabad(Ramachandran and Rajesh,2007)and Delhi(Srivastava et al.,2012a).Cao et al.(2009) observed a strong inverse relationship(–0.66)between BC and WS at Xi'an,China,from September 2003 to August 2005,implying a local origin for BC.These studies suggested important contributions from locally generated BC that accumulate at low WS.Sharma et al.(2002)reported that low WS led to poorer dispersion such that BC remained concentrated around the emissions source.They also found that higher WS was the driving force for dilution of the BC concentration near the surface due to constant traffic density.These dilution conditions are a general feature over urban and suburban sites,and largely show negative correlation between BC and WS.The strength of the correlation between WS and BC is an indicator of the proximity of BC sources to the measurement site. Low correlation coefficients suggest that BC originates from distant sources.The mean monthly temperature in Varanasi was approximately 27℃.It was lower(approximately 19℃)during winter and higher(approximately 32℃)during the pre-monsoon period.

    Low correlation coefficients suggest that BC originates from distant sources.The mean monthly temperature in Varanasi was approximately 27℃.It was lower(approximately 19℃)during winter and higher (approximately 32℃)during the pre-monsoon season. An inverse relationship(–0.32)was observed between temperature and BC(Fig.8),similar to central and southern parts of India at Ahmedabad(Ramachandran and Rajesh,2007)and Trivendrum(Babu and

    Moorthy,2002).The high BC emissions during the late post-monsoon and winter seasons,from agricultural burning and biomass burning for heating and cooking,along with the lower MH at that time,led to the negative correlation between temperature and BC. During the study period,the seasonal MHs were 1277, 628,556,and 370 m during the pre-monsoon,monsoon,winter,and post-monsoon seasons,respectively. The MH was shallower during the post-monsoon and winter seasons over Varanasi,resulting in the trapping of locally emitted pollutants and higher BC concentrations.The increase in surface temperature with associated convective activity during the pre-monsoon season improved the dispersion,causing lower surface BC concentrations.Given the volume into which groundlevel emissions are dispersed,a negative relationship between BC and MH was observed(–0.66).Sloane and White(1986)suggested that the loss of VIS is an easily measured manifestation of air pollution,arising from the loss of contrast between the object and the background and attenuation of the light signal from the object due to scattering and absorption of light by fine particles and other atmospheric pollutants.Regression analysis between VIS and BC yielded a large negative correlation(–0.84).Horvath(1995)suggested the absorption of sunlight by BC contributes to reduce VIS in polluted regions.Xu et al.(2012)reported a negative effect of BC on VIS,with a significant correlation of–0.79 between atmospheric VIS and optical properties(σscatand σabs),at an urban site in Shanghai,China.

    Fig.8.Correlation coefficients between BC and the meterological parameters of MH,TEM,RH,VIS,and WS during the study period over Varanasi.

    A positive correlation between BC and RH was 0.41.RH is high during the monsoon season,and low during the post-monsoon and winter seasons when there are higher emissions;thus,this correlation likely reflects the seasonal variation as discussed previously.

    Regression analysis between BC and annual precipitation was studied and a significant negative correlation(–0.35)was found between them(annual rainfall of 876 mm in 2009).Rainfall is high during the monsoon season,so part of this negative correlation may again reflect the winter/monsoon seasonal differences already presented.In addition,the precipitation will remove some particles from the air through washout. A similar correlation(–0.35)between BC and rainfall was also reported in Ahmadabad(western India) (Ramachandran and Rajesh,2007).A strong correlation(–0.74)between these variables at the Trivandrum coastal station(southern India)was reported because of the heavy rainfall(Babu and Moorthy,2002). To analyze the transport of BC,the surface wind direction data along with the BC mass concentration over Varanasi was separated and five major wind categories were identified:east(comprising northeasterly, east-northeasterly,easterly,and north-northeasterly wind);west(comprising westerly,west-southwesterly, south-southwesterly,and southwesterly wind);north (comprising northerly,north-northwesterly,northwesterly,and west-northwesterly wind);south(comprising southerly,southeasterly,east-southeasterly, and south-southeasterly wind);and calm(meaning no wind).The observed percentage frequencies of wind during the study period resulted in the order:calm (38%),west(31%),east(18%),south(10%),and north(3%);however,the BC mass concentrations were highest(38%;18.1±6.9μg m?3)during calm conditions,confirming that the major source of BC was localized.

    4.Summary

    In-situ measurements of BC mass concentrations were made in Varanasi(an urban environment)in the central IGP region of India.For the year 2009,the results suggest the importance of BC at this location

    and the need for extensive mitigation activities to reduce the concentrations.The annual mean mass concentration of BC was 8.92±6.98 g m?3,with 34% of the days exceeding this average value.The highest (approximately 46%)mass BC frequency was in the range 5–10μg m?3;however,during the post-monsoon season,the higher contribution(23%)was between 20 and 25 g m?3.The BC mass showed gradual buildup at around 0700 LT,peaking at around 0800–0900 LT,with low values in the afternoon(1200–1900 LT) and secondary maxima between 2100 LT and midnight (0100 LT).Overall,the nighttime concentrations of BC were approximately two times higher than that during the day,due to the impact of boundary layer conditions.AE was significantly positively correlated (0.55)with surface BC,indicating the impact of BC on columnar aerosols.The estimated mean AAE was 1.02±0.08,indicating the combustion of fossil fuel sources.The significant negative correlation(>0.4) between BC mass and meteorological parameters indicated a pronounced effect on atmospheric dynamics of the enhancement of BC mass over this region.The highest concentration(38%;18.1±6.9μg m?3)of BC mass was found under calm wind conditions.

    In view of the above,urgent action is needed in reducing the current level of BC in Varanasi.Concentrations of BC in the atmosphere over Varanasi are much higher than that in other parts of India and elsewhere around the globe.

    Acknowledgments.The authors are grateful for the financial support of the ISRO-ARFI program. They also gratefully acknowledge the NOAA Air Resources Laboratory for the provision of the HYSPLIT transport and dispersion model and ready website(http://www.arl.noaa.gov/ready.html)used in the current study.Also acknowledged are the data from the Giovanni online data system,and the comments from the anonymous reviewers.

    Andrea,D.A.C.,and P.Artaxo,2001:Wintertime and summertime S?ao Paulo aerosol source apportionment study.Atmos.Environ.,35,4889–4902.

    Aruna,K.,T.V.L.Kumar,D.N.Rao,et al.,2013:Black carbon aerosols in a tropical semi-urban coastal environment:Effects of boundary layer dynamics and long range transport. J.Atmos. Solar-Terrest. Phys.,104,116–125.

    Arnott,W.P.,K.Hamasha,H.Moosm¨uller,et al.,2005: Towards aerosol light absorption measurements with a 7-wavelength aethalometer:Evaluation with a photoacoustic instrument and 3-wavelength nephelometer.Aerosol Sci.Technol.,39,17–29.

    Awasthy,V.,M.K.Kumar,and S.K.Satheesh,2010: Measurements of aerosol black carbon at an urban site in southern India.Aerosol and Clouds:Climate Change Perspectives,IASTA conference,19,463–466.

    Babu,S.S.,and K.K.Moorthy,2002:Aerosol black carbon over a tropical coastal station in India.Geophys. Res.Lett.,29,2098,doi:10.1029/2002GL015662.

    Bano,T.,S.Singh,N.C.Gupta,et al.,2011:Variation in aerosol black carbon concentration and its emission estimates at the mega-city Delhi.Int.J.Remote Sens.,32,6749–6764.

    Beegum,S.N.,K.K.Moorthy,S.S.Babu,et al.,2009: Spatial distribution of aerosol black carbon over India during pre-monsoon season.Atmos.Environ., 43,1071–1078.

    Bergstrom,R.W.,P.Pilewskie,J.Pommier,et al.,2004: Spectral absorption of solar radiation by aerosols during ACE-Asia.J.Geophys.Res.,109,D19S15, doi:10.1029/2003JD004467.

    Bisht,D.S.,U.C.Dumka,D.G.Kaskaoutis,et al.,2015: Carbonaceous aerosols and pollutants over Delhi urban environment: Temporal evolution,source apportionment,and radiative forcing.Sci.Total Environ.,521–522,431–445.

    Bond,T.C.,2001:Spectral dependence of visible light absorption by carbonaceous particles emitted from coal combustion.Geophys.Res.Lett.,28,4075–4078.

    Bond,T.C.,and R.W.Bergstrom,2006:Light absorption by carbonaceous particles: An investigative review.Aerosol Sci.Tech.,40,27–67.

    Bond,T.C.,S.J.Doherty,D.W.Fahey,et al.,2013: Bounding the role of black carbon in the climate system:A scientific assessment.J.Geophys.Res., 118,5380–5552.

    Cao,J.J.,C.S.Zhu,J.C.Chow,et al.,2009:Black carbon relationships with emissions and meteorology in Xi'an,China.Atmos.Res.,94,194–202.

    Chen,L.W.A.,B.G.Doddridge,R.R.Dickerson,et al.,2001:Seasonal variations in elemental carbon

    aerosol,carbon monoxide and sulfur dioxide:Implications for sources. Geophys. Res. Lett.,28, 1711–1714.

    Collaud,C.M.,E.Weingartner,A.Apituley,et al.,2010: Minimizing light absorption measurement artifacts of the Aethalometer:Evaluation of five correction algorithms.Atmos.Meas.Tech.3,457–474.

    Corrigan,C.E.,V.Ramanathan,and J.J.Schauer,2006: Impact of monsoon transitions on the physical and optical properties of aerosols. J.Geophys. Res., 111,D18208,doi:10.1029/2005JD006370.

    Dumka,U.C.,P.R.Sinha,R.K.Manchanda,et al., 2010a:Seasonal and diurnal variation of black carbon aerosols over tropical urban site Hyderabad India.Aerosol and Clouds:Climate Change Perspectives,IASTA Conference,19(1&2),440–441.

    Dumka,U.C.,K.K.Moorthy,R.Kumar,et al.,2010b: Characteristics of aerosol black carbon mass concentration over a high altitude location in the Central Himalayas from multi-year measurements.Atmos. Res.,96,510–521.

    Dumka,U.C.,R.K.Manchanda,P.R.Sinha,et al., 2013:Temporal variability and radiative impact of black carbon aerosol over tropical urban station Hyderabad.J.Atmos.Solar-Terrestr.Phys.,105–106,81–90.

    Ganguly,D.,A.Jayaraman,H.Gadhavi,et al.,2005: Features in wavelength dependence of aerosol absorption observed over central India.Geophys.Res. Lett.,32,L13821,doi:10.1029/2005GL023023.

    Horvath,H.,1995:Estimation of the average visibility in central Europe.Atmos.Environ.,29,241–246.

    Huang,X.F.,T.L.Sun,L.W.Zeng,et al.,2012: Black carbon aerosol characterization in a coastal city in South China using a single particle soot photometer. Atmos. Environ.,51,21–28,doi: 10.1016/j.atmosenv.2012.01.056.

    Husain,L.,V.A.Dutkiewics,A.J.Khan,et al.,2007: Characterization of carbonaceous aerosols in urban air.Atmos.Environ.,41,6872–6883.

    Hyv¨arinen,A.P.,H.Lihavainen,M.Komppula,et al., 2009:Continuous measurements of optical properties of atmospheric aerosols in Mukteshwar,northern India. J.Geophys.Res.,114,D08207,doi: 10.1029/2008JD011489.

    Hyv¨arinen,A.P.,H.Lihavainen,M.Komppula,et al., 2010:Aerosol measurements at the Gual Pahari EUCAARI station:Preliminary results from in-situ measurements.Atmos.Chem.Phys.,10,7241–7252.

    Hyv¨arinen,A.P.,P.Kolmonen,V.-M.Kerminen,et al., 2011:Aerosol black carbon at five background measurement sites over Finland,a gateway to the Arctic. Atmos.Environ.,45,4042–4050.

    Invernizzi,G.,A.Ruprecht,R.Mazza,et al.,2011:Measurement of black carbon concentration as an indicator of air quality benefits of traffic restriction policies within the ecopass zone in Milan,Italy.Atmos.Environ.,45,3522–3527.

    IPCC,2007:Changes in atmospheric constituents and in radiative forcing.Climate Change 2007:The Physical Science Basis.Contribution of Working Group to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change.S.Solomon,et al.,Eds.,Cambridge University Press,Cambridge, UK and New York,NY,USA,996 pp.

    Jacobson,M.Z.,2001:Strong radiative heating due to the mixing state of black carbon in atmospheric aerosols.Nature,409,695–697.

    Jacobson,M.Z.,2002:Control of fossil-fuel particulate black carbon and organic matter,possibly the most effective method of slowing global warming.J.Geophys.Res.,107,4410.doi:10.1029/2001JD001376.

    Kaufman,Y.J.,D.Tanr′e,and O.Boucher,2002:A satellite view of aerosols in the climate system.Nature,419,215–223.

    Kirchstetter,T.W.,T.Novakov,and P.V.Hobbs, 2004: Evidence that the spectral dependence of light absorption by aerosols is affected by organic carbon. J.Geophys. Res.,109,D21208,doi: 10.1029/2004JD004999.

    Kumar,M.,S.Tiwari,V.Murari,et al.,2015:Wintertime characteristics of aerosols at middle Indo-Gangetic Plain:Impacts of regional meteorology and long range transport.Atmos.Environ.,104, 162–175.

    Kuniyal,J.C.,2010:Aerosols climatology over the northwestern Indian Himalayan region.ARFI&ICARB Scientific Progress Report ISRO-GBP,India,93–99.

    Lau,K.M.,and K.M.Kim,2006:Observational relationship between aerosol and Asian monsoon rainfall, and circulation.Geophys.Res.Lett.,33,L21810, doi:10.1029/2006GL027546.

    Liao Hong and Shang Jingjing,2015:Regional warming by black carbon and tropospheric ozone:A review of progress and research challenges in China. J. Meteor.Res.,29,525–545.

    Lodhi,N.K.,S.N.Beegum,S.Singh,et al.,2013: Aerosol climatology at Delhi in the western Indo-Gangetic Plain: Microphysics,long-term trends and source strengths.J.Geophys.Res.,118,doi: 10.1002/jgrd.50165.

    Lyamani,H.,F.J.Olmo,I.Foyo,et al.,2011:Black carbon aerosols over an urban area in southeastern Spain:Changes detected after the 2008 economic crisis.Atmos.Environ.,45,6423–6432.

    Menon,S.,J.Hansen,L.Nazarenko,et al.,2002:Climate effects of black carbon aerosols in China and India. Science,297,2250–2253.

    Mishra,A.K.,T.Shibata,2012:Synergistic analyses of optical and microphysical properties of agricultural crop residue burning aerosols over the Indo-Gangetic Basin(IGB).Atoms.Environ.,57,205–218.

    Mohan,M.,and S.Bhati,2009: Why is megacity Delhi prone to high atmospheric pollution potential? TFMM-TF HTAP Joint Workshop,Paris, June 2009.

    Moorthy,K.K.,and S.K.Satheesh,2011:Black carbon aerosols over India.Black Carbon e-Bulletin,ENEP, 3.

    Moosm¨uller,H.,R.K.Chakrabarty,and W.P.Arnott, 2009:Aerosol light absorption and its measurement: A review.J.Quant.Spectrosc.Rad.Trans.,110, 844–878.

    Pandithurai,G.,S.Dipu,K.K.Dani,et al.,2008: Aerosol radiative forcing during dust events over New Delhi,India.J.Geophys.Res.,113,D13209, doi:10.1029/2008JD009804.

    Pant,P.,P.Hegde,U.C.Dumka,et al.,2006:Aerosol characteristics at a high-altitude location during ISRO-GBP land campaign-II.Curr.Sci.,91,1053–1061.

    Pani,S.K.,and S.Verma,2010:Black carbon and its contribution to aerosol optical depth over Kolkata on the eastern IGP.IASTA Conference Bulletin,19, 458–460.

    Pereira,S.N.,F.Wagner,and A.M.Silva,2012:Long term black carbon measurements in the southwestern Iberia Peninsula.Atmos.Environ.,57,63–71.

    Prasad,A.K.,R.P.Singh,M.Kafatos,et al.,2005:Proceedings of symposium S6 held during the Seventh IAHS Scientific Assembly at Foz do Igua?cu,Brazil, April 2005.IAHS Publ.,296.

    Praveen,P.S.,T.Ahmed,A.Kar,et al.,2012:Link between local scale BC emissions in the indo-gangetic plains and large scale atmospheric solar absorption.Atmos.Chem.Phys.,12,1173–1187,doi: 10.5194/acp-12-1173-2012.

    Raju,M.P.,P.D.Safai,P.S.P.Rao,et al.,2011:Seasonal characteristics of black carbon aerosols over a high altitude station in Southwest India.Atmos. Res.,100,103–110.

    Ramachandran,S.,2004:Spectral aerosol optical characteristicsduring thenortheastmonsoon over the Arabian Sea and the tropical Indian Ocean. 2: ?Angstr¨om parameters and anthropogenic influence. J.Geophys. Res.,109,D19208,doi: 10.1029/2003JD004483.

    Ramachandran,S.,and T.A.Rajesh,2007:Black carbon aerosol mass concentrations over Ahmedabad, an urban location in western India:Comparison with urban sites in Asia,Europe,Canada,and the United States.J.Geophys.Res.,112,D06211.

    Ramanathan,V.,and G.Carmichael,2008:Global and regional climate changes due to black carbon.Nature Geosci.,1,221–227.

    Ramachandran,S.,and S.Kedia,2010:Black carbon aerosols over an urban region:Radiative forcing and climate impact.J.Geophys.Res.,115,D10202, doi:10.1029/2009JD013560.

    Raysoni,A.U.,J.A.Sarnat,S.E.Sarnat,et al.,2011: Binational school-based monitoring of traffic-related air pollutants in El Paso,Texas(USA)and Ciudad Ju′arez,Chihuahua(M′exico).Environ.Pollut.,159, 2476–2486.

    Reddy,M.S.,and C.Venkataraman,2002:Inventory of aerosol and sulphur dioxide emissions from India. Part II:Biomass combustion.Atmos.Environ.,36, 699–712.

    Ruellan,S.,and H.Cachier,2001:Characterisation of fresh particulate vehicular exhausts near a paris high flow road.Atmos.Environ.,35,453–468.

    Russell,P.B.,R.W.Bergstrom,Y.Shinozuka,et al., 2010:Absorption Angstrom exponent in AERONET and related data as an indicator of aerosol composition.Atmos.Chem.Phys.,10,1155–1169.

    Safai,P.D.,S.Kewat,G.Pandithurai,et al.,2008: Aerosol characteristics during winter fog at Agra, North India.J.Atmos.Chem.,61,101–118.

    Safai,P.D.,S.Kewat,P.S.Praveen,et al.,2007:Seasonal variation of black carbon aerosols over tropical urban city of Pune,India.Atmos.Environ.,41, 2699–2709.

    Safai,P.D.,M.P.Raju,K.B.Budhavant,et al.,2013: Long term studies on characteristics of black carbon aerosols over a tropical urban station Pune,India. Atmos.Res.,132–133,173–184.

    Sandradewi,J.,A.S.Pr′ev?ot,S.Szidat,et al.,2008a: Using aerosol light absorption measurements for the quantitative determination of wood burning and traffic emission contributions to particulate matter. Environ.Sci.Technol.,42,3316–3323.

    Sandradewi,J.,A.S.Pr′ev?ot,E.Weingartner,et al., 2008b: A study of wood burning and traffic aerosols in an Alpine valley using a multi-wavelength Aethalometer.Atmos.Environ.,42,101–112.

    Satheesh,S.K.,and V.Ramanathan,2000:Large differences in tropical aerosol forcing at the top of the atmosphere and earth's surface.Nature,405,60–63, doi:10.1038/35011039.

    Sharma,S.,J.R.Brook,H.Cachier,et al.,2002:Light absorption and thermal measurements of black carbon in different regions of Canada. J.Geophys. Res.,107,4771,doi:10.1029/2002JD002496.

    Sharma,A.R.,S.K.Kharol,K.V.S.Badarinath,et al., 2010:Impact of agriculture crop residue burning on atmospheric aerosol loading-a study over Punjab State,India.Ann.Geophys.,28,367–379.

    Sheridan,P.J.,W.P.Arnott,J.Ogren,et al.,2005: The Reno aerosol optics study:An evaluation of aerosol absorption measurement methods.Aerosol Sci.Technol.,39,1–16.

    Simpson,A.J.,and O.S.McGee,1996:Analysis of the fumigation effect on pollutants over Pietermaritzburg.South African Geograph.J.,78,41–46.

    Singh,S.,S.Nath,R.Kohli,et al.,2005:Aerosols over Delhi during pre-monsoon months:Characteristics and effects on surface radiation forcing.Geophys. Res.Lett.,32,L13808,doi:10.1029/2005GL023062.

    Singh,B.P.,A.K.Srivastava,S.Tiwari,et al.,2014: Radiative impact of fireworks at a tropical Indian location:A case study.Adv.Meteor.,2014,Article ID 197072,8.

    Singh,S.,S.Tiwari,D.P.Gond,et al.,2015:Intraseasonal variability of black carbon aerosols over a coal field area at Dhanbad,India.Atmos.Res., 161–162,25–35.

    Sloane,C.S.,and W.H.White,1986:Visibility:An evolving issue.Environ.Sci.Technol.,20,760–766.

    Soni,K.,S.Singh,T.Bano,et al.,2010:Variations in single scattering albedo and Angstrom absorption exponent during different seasons at Delhi,India. Atmos.Environ.,44,4355–4363.

    Sreekanth,V.,K.Niranjan,and B.L.Madhavan,2007: Radiative forcing of black carbon over eastern India. Geophys. Res. Lett.,34,L17818,doi: 10.1029/2007GL030377.

    Srivastava,A.K.,K.Ram,P.Pant,et al.,2012a:Black carbon aerosols over Manora Peak in the Indian Himalayan foothills:Implications for climate forcing. Environ.Res.Lett.,7,014002,doi:10.1088/1748-9326/7/1/014002.

    Srivastava,A.K.,S.Tiwari,P.C.S.Devara,et al.,2011b: Pre-monsoon aerosol characteristics over the Indo-Gangetic Basin:Implications to climatic impact. Annal.Geophys.,29,789–804,doi:10.5194/angeo-29-789-2011.

    Srivastava,A.K.,S.Sachchidanand,P.Pant,et al., 2012b:Characteristics of black carbon over Delhi and Manora Peak-a comparative study.Atmos.Sci. Lett.,13,223–230.

    Srivastava,A.K.,D.S.Bisht,K.Ram,et al.,2014: Characterization of carbonaceous aerosols over Delhi in Ganga basin:Seasonal variability and possible sources.Environ.Sci.Pollut.Res.,21,8610–8619.

    Srivastava,M.K.,S.K.Srivastava,A.Saha,et al., 2011a:Aerosol optical properties over Delhi and Manora Peak during a rare dust event in early April 2005.Internat.J.Remote Sens.,32,7939–7954, doi:10.1080/01431161.2010.523732.

    Stocker,T.F.,D.Qin,G.-K.Plattner,et al.,2014: Climate Change 2013:The Physical Science Basis. Working Group I Contribution to the Fifth Assessment Report of the Intergovernmental Panel On Climate Change.Cambridge University Press,New York,N.Y.115 pp.

    Stone,E.A.,J.J.Schauer,B.B.Pradhan,et al.,2010: Characterization of emissions from South Asian biofuels and application to source apportionment of carbonaceous aerosol in the Himalayas.J.Geophys. Res.,115,D06301,doi:10.1029/2009JD011881.

    Tripathi,S.N.,S.Dey,V.Tare,et al.,2005:Aerosol black carbon radiative forcing at an industrial city in northern India.Geophys.Res.Lett.,32,L08802.

    Tiwari,S.,A.K.Srivastava,D.S.Bisht,et al.,2009: Black carbon and chemical characteristics of PM10and PM2.5at an urban site of North India.J.Atmos.Chem.,62,193–209.

    Tiwari,S.,A.K.Srivastava,D.S.Bisht,et al.,2013: Diurnal and seasonal variations of black carbon and PM2.5over New Delhi,India:Influence of meteorology.Atmos.Res.,125–126,50–62,doi:10.1016/ j.atmosres.2013.01.011.

    Tiwari,S.,U.C.Dumka,D.G.Kaskaoutis,et al., 2015a:Aerosol chemical characterization and role of carbonaceous aerosol on radiative effect over Varanasi in central Indo-Gangetic Plain. Atmos. Environ.,(in press),http://dx.doi.org/10.1016/j.atmosenv.2015.07.031.

    Tiwari,S.,A.K.Srivastava,A.K.Singh,et al.,2015b: Identification of aerosol types over Indo-Gangetic Basin:Implications to optical properties and associated radiative forcing.Environ.Sci.Pollut.Res., 22,12246–12260,doi:10.1007/s11356-015-4495-6.

    Virkkula,A.,T.M¨akel¨a,R.Hillamo,et al.,2007:A simple procedure for correcting loading effects of aethalometer data.J.Air Waste Manage.Assoc., 57,1214–1222.

    Vyas,B.M.,2010:Studies of regional features of atmospheric aerosol,total carbonaceous aerosols and their role in the atmospheric radiative forcing effect over the tropical semi-arid location,i.e.,Udaipur, western region part of India.ARFI&ICARB Scientific Progress Report ISRO-GBP,India,67–70.

    Wang,C.,2004:A modeling study on the climate impacts of black carbon aerosols.J.Geophys.Res., 109,D03106,doi:10.1029/2003JD004084.

    Wang Xin,Xu Baiqing,and Ming Jing,2014:An overview of the studies on black carbon and mineral dust deposition in snow and ice cores in East Asia. J.Meteor.Res.,28,354–370.

    Weingartner,E.,H.Saathoff,M.Schnaiter,et al.,2003: Absorption of light by soot particles: Determination of the absorption coefficient by means of AETHALOMETERS.J.Aero.Sci.,34,1445–1463.

    Xu,J.W.,J.Tao,R.J.Zhang,et al.,2012:Measurements of surface aerosol optical properties in winter of Shanghai.Atmos.Res.,109–110,25–35.

    :Singh,B.P.,S.Tiwari,K.Hopke Philip,et al.,2015:Seasonal inhomogeneity of soot particles over the central Indo-Gangetic Plains,India:Influence of meteorology.J.Meteor.Res.,29(6), 935–949,

    10.1007/s13351-015-5041-7.

    Supported by the Indian Space Research Organization(ISRO)“Aerosol Radiative Forcing over India”Program at Bananas Hindu University(P32/15 and P32/14).

    ?Corresponding author:mksriv@gmail.com.

    ?The Chinese Meteorological Society and Springer-Verlag Berlin Heidelberg 2015

    (Received April 30,2015;in final form September 22,2015)

    欧美性猛交╳xxx乱大交人| 极品教师在线免费播放| 天天躁夜夜躁狠狠躁躁| 亚洲一区中文字幕在线| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片午夜丰满| 亚洲国产精品成人综合色| 久久久久久国产a免费观看| 日韩欧美精品v在线| 变态另类成人亚洲欧美熟女| 一级毛片精品| 亚洲av成人精品一区久久| √禁漫天堂资源中文www| 国产视频内射| 99久久国产精品久久久| 在线观看免费日韩欧美大片| 久久久精品国产亚洲av高清涩受| 日韩av在线大香蕉| 最近视频中文字幕2019在线8| 久久婷婷人人爽人人干人人爱| 久久人人精品亚洲av| 在线观看66精品国产| 18禁美女被吸乳视频| 亚洲国产欧洲综合997久久,| 天堂√8在线中文| 国产精品av久久久久免费| 国产精品av久久久久免费| 亚洲九九香蕉| 久久香蕉国产精品| 欧美性长视频在线观看| 国产激情久久老熟女| 在线观看免费午夜福利视频| 听说在线观看完整版免费高清| 天堂动漫精品| 国产免费av片在线观看野外av| 悠悠久久av| 麻豆成人av在线观看| 韩国av一区二区三区四区| 成年女人毛片免费观看观看9| 国产成人aa在线观看| 久久亚洲真实| 免费看十八禁软件| 欧美日韩福利视频一区二区| 精品久久久久久久久久免费视频| 日韩中文字幕欧美一区二区| 国产蜜桃级精品一区二区三区| 国产精品久久久久久人妻精品电影| 一进一出抽搐动态| 一级毛片女人18水好多| av福利片在线观看| av在线播放免费不卡| 性色av乱码一区二区三区2| av福利片在线| 国产精品电影一区二区三区| 亚洲精品在线美女| 人妻夜夜爽99麻豆av| 全区人妻精品视频| 色哟哟哟哟哟哟| 妹子高潮喷水视频| 久久精品人妻少妇| 欧美+亚洲+日韩+国产| 国产一级毛片七仙女欲春2| 亚洲中文字幕日韩| 欧美乱色亚洲激情| √禁漫天堂资源中文www| 久久久国产精品麻豆| 国产一区二区在线观看日韩 | 999久久久国产精品视频| 日日爽夜夜爽网站| 亚洲第一电影网av| 人妻久久中文字幕网| 女警被强在线播放| 久99久视频精品免费| 国产aⅴ精品一区二区三区波| 不卡av一区二区三区| 好男人电影高清在线观看| 日韩欧美精品v在线| 久久久久久国产a免费观看| a级毛片在线看网站| 精品午夜福利视频在线观看一区| 欧美日韩国产亚洲二区| 老汉色∧v一级毛片| 天堂√8在线中文| 男插女下体视频免费在线播放| 亚洲人成伊人成综合网2020| 一卡2卡三卡四卡精品乱码亚洲| 亚洲在线自拍视频| 搡老熟女国产l中国老女人| 老熟妇仑乱视频hdxx| 婷婷六月久久综合丁香| 99riav亚洲国产免费| 亚洲avbb在线观看| 欧美乱色亚洲激情| 精品久久久久久,| 女生性感内裤真人,穿戴方法视频| 人成视频在线观看免费观看| 国产又黄又爽又无遮挡在线| 精品国产超薄肉色丝袜足j| 久久天堂一区二区三区四区| 99riav亚洲国产免费| 色综合欧美亚洲国产小说| 亚洲最大成人中文| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲欧美日韩高清专用| 精品一区二区三区av网在线观看| 国产高清视频在线播放一区| a级毛片a级免费在线| 这个男人来自地球电影免费观看| 一本精品99久久精品77| 亚洲一区中文字幕在线| 五月玫瑰六月丁香| 狂野欧美白嫩少妇大欣赏| 午夜精品一区二区三区免费看| 全区人妻精品视频| 在线免费观看的www视频| 亚洲一卡2卡3卡4卡5卡精品中文| 看免费av毛片| 久久欧美精品欧美久久欧美| 90打野战视频偷拍视频| 亚洲人与动物交配视频| 亚洲专区字幕在线| 美女大奶头视频| 国产伦人伦偷精品视频| 日本a在线网址| 老司机福利观看| 99国产精品一区二区三区| 亚洲自拍偷在线| 露出奶头的视频| 成人国语在线视频| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| 亚洲av熟女| 首页视频小说图片口味搜索| 香蕉久久夜色| 国产又色又爽无遮挡免费看| 中国美女看黄片| 黄色丝袜av网址大全| 日韩欧美在线乱码| avwww免费| av超薄肉色丝袜交足视频| 亚洲男人的天堂狠狠| 又黄又爽又免费观看的视频| 久久久国产欧美日韩av| 欧美日韩亚洲综合一区二区三区_| 中出人妻视频一区二区| 白带黄色成豆腐渣| 久久中文看片网| 18禁裸乳无遮挡免费网站照片| 国产av一区在线观看免费| 国产熟女xx| 午夜视频精品福利| 亚洲av熟女| 一卡2卡三卡四卡精品乱码亚洲| 色老头精品视频在线观看| 在线视频色国产色| 日本免费a在线| 国产成人欧美在线观看| 亚洲色图av天堂| 一区福利在线观看| 丝袜人妻中文字幕| 脱女人内裤的视频| 国产野战对白在线观看| 欧美成人一区二区免费高清观看 | 中文字幕av在线有码专区| 日韩欧美 国产精品| 麻豆国产97在线/欧美 | 国产激情欧美一区二区| 人人妻人人澡欧美一区二区| 黄色视频,在线免费观看| 非洲黑人性xxxx精品又粗又长| 看片在线看免费视频| 免费在线观看影片大全网站| 床上黄色一级片| 精品欧美国产一区二区三| 看免费av毛片| 国产亚洲精品第一综合不卡| 国产亚洲精品久久久久5区| 性欧美人与动物交配| 少妇熟女aⅴ在线视频| 亚洲成人中文字幕在线播放| 嫁个100分男人电影在线观看| 久久性视频一级片| 国产欧美日韩一区二区三| 国产精品 国内视频| 高清在线国产一区| 国产精品综合久久久久久久免费| 久久中文字幕一级| 人妻久久中文字幕网| 国语自产精品视频在线第100页| 婷婷精品国产亚洲av在线| 精品国产美女av久久久久小说| 午夜免费激情av| 国产成人av激情在线播放| 在线观看舔阴道视频| 国产成人系列免费观看| 男人的好看免费观看在线视频 | netflix在线观看网站| 国产成人系列免费观看| 毛片女人毛片| 国产激情偷乱视频一区二区| 久久人人精品亚洲av| 欧美zozozo另类| 色尼玛亚洲综合影院| 成人特级黄色片久久久久久久| 久久这里只有精品中国| av超薄肉色丝袜交足视频| 精品日产1卡2卡| 老熟妇乱子伦视频在线观看| 免费高清视频大片| 国产97色在线日韩免费| 日韩中文字幕欧美一区二区| 国产精品日韩av在线免费观看| 亚洲中文日韩欧美视频| www日本黄色视频网| 给我免费播放毛片高清在线观看| 一级毛片精品| 美女免费视频网站| 久久亚洲真实| 精品久久久久久久人妻蜜臀av| 一区二区三区激情视频| 极品教师在线免费播放| 直男gayav资源| 国产成人精品久久久久久| 91久久精品国产一区二区三区| 亚洲一区高清亚洲精品| 在线a可以看的网站| 亚洲18禁久久av| 51国产日韩欧美| 一个人看视频在线观看www免费| 一区福利在线观看| 一进一出抽搐动态| 人妻制服诱惑在线中文字幕| 国产精品日韩av在线免费观看| 国产一级毛片在线| 欧美区成人在线视频| 中文精品一卡2卡3卡4更新| 三级毛片av免费| 免费人成在线观看视频色| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 人人妻人人澡人人爽人人夜夜 | 波多野结衣高清无吗| 在线观看午夜福利视频| 看免费成人av毛片| 狂野欧美激情性xxxx在线观看| 中文字幕人妻熟人妻熟丝袜美| av在线观看视频网站免费| 毛片女人毛片| 国语自产精品视频在线第100页| 日韩强制内射视频| 国内精品久久久久精免费| 3wmmmm亚洲av在线观看| 国产69精品久久久久777片| 精品一区二区三区人妻视频| 美女内射精品一级片tv| 尾随美女入室| 国产一区二区在线av高清观看| 亚洲精品乱码久久久v下载方式| 日本五十路高清| 亚洲欧美清纯卡通| 久久鲁丝午夜福利片| 国内揄拍国产精品人妻在线| 最近手机中文字幕大全| 99久久九九国产精品国产免费| 我的女老师完整版在线观看| 欧美3d第一页| 天堂中文最新版在线下载 | 能在线免费观看的黄片| 日本熟妇午夜| 欧美一区二区亚洲| 亚洲欧美日韩高清在线视频| 精华霜和精华液先用哪个| 插阴视频在线观看视频| 可以在线观看毛片的网站| 黄片wwwwww| 亚洲av免费高清在线观看| av在线亚洲专区| 高清日韩中文字幕在线| 51国产日韩欧美| 国产午夜福利久久久久久| 听说在线观看完整版免费高清| 亚洲国产精品国产精品| 又粗又爽又猛毛片免费看| 欧美成人一区二区免费高清观看| 日韩欧美在线乱码| 国产精品久久久久久精品电影| 一进一出抽搐动态| 色综合色国产| 成人特级av手机在线观看| 十八禁国产超污无遮挡网站| 婷婷六月久久综合丁香| 国产精品久久电影中文字幕| 亚洲成人中文字幕在线播放| 色视频www国产| 国产亚洲精品久久久久久毛片| 国产精品精品国产色婷婷| 老女人水多毛片| 亚洲精品乱码久久久久久按摩| 我要看日韩黄色一级片| 高清在线视频一区二区三区 | 国产精品国产高清国产av| 亚洲av.av天堂| 秋霞在线观看毛片| 亚洲第一区二区三区不卡| 免费大片18禁| 丝袜喷水一区| 亚洲综合色惰| 性欧美人与动物交配| 在线天堂最新版资源| av在线天堂中文字幕| 99久久人妻综合| 91av网一区二区| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 亚洲欧洲日产国产| 国内久久婷婷六月综合欲色啪| 91av网一区二区| 性色avwww在线观看| 欧美最黄视频在线播放免费| 久久这里只有精品中国| 春色校园在线视频观看| 色综合色国产| 精品久久久久久久久av| 又爽又黄无遮挡网站| 国产亚洲精品久久久久久毛片| 亚洲精品国产成人久久av| 亚洲不卡免费看| 亚洲欧美精品自产自拍| 男女视频在线观看网站免费| 日本熟妇午夜| 九九热线精品视视频播放| 久久草成人影院| 麻豆久久精品国产亚洲av| 国产一区二区三区在线臀色熟女| 成人欧美大片| 我的女老师完整版在线观看| 老司机福利观看| 男女啪啪激烈高潮av片| 久久午夜福利片| 少妇熟女aⅴ在线视频| 欧美一区二区国产精品久久精品| 真实男女啪啪啪动态图| 亚洲美女视频黄频| 国产视频内射| 日本黄色视频三级网站网址| 亚洲欧美清纯卡通| 国产白丝娇喘喷水9色精品| 久久久精品欧美日韩精品| 国产在线精品亚洲第一网站| 日本色播在线视频| 欧美丝袜亚洲另类| 亚洲欧美精品自产自拍| 毛片一级片免费看久久久久| 在线播放无遮挡| 国产一区二区三区在线臀色熟女| 免费观看精品视频网站| 美女黄网站色视频| 午夜精品国产一区二区电影 | 99久久精品热视频| 一级av片app| 人妻夜夜爽99麻豆av| 日日摸夜夜添夜夜爱| 波多野结衣巨乳人妻| 国产精品,欧美在线| 亚洲成a人片在线一区二区| 欧美精品一区二区大全| 午夜激情福利司机影院| 最好的美女福利视频网| 日韩 亚洲 欧美在线| 夜夜看夜夜爽夜夜摸| 69av精品久久久久久| 成人欧美大片| 狂野欧美激情性xxxx在线观看| 欧美一级a爱片免费观看看| 欧美bdsm另类| 青春草国产在线视频 | 国产黄片视频在线免费观看| 成人欧美大片| 国产亚洲91精品色在线| 午夜久久久久精精品| 黄片无遮挡物在线观看| 日韩亚洲欧美综合| 精品不卡国产一区二区三区| 欧美激情在线99| 国产伦在线观看视频一区| 国产av不卡久久| 国产久久久一区二区三区| 人人妻人人澡欧美一区二区| 三级毛片av免费| 99视频精品全部免费 在线| 99riav亚洲国产免费| 少妇高潮的动态图| 嫩草影院新地址| 亚洲欧美成人精品一区二区| 免费电影在线观看免费观看| 国产在视频线在精品| 国产三级在线视频| 毛片一级片免费看久久久久| 日本熟妇午夜| 日本免费a在线| 免费av观看视频| 国产亚洲91精品色在线| 一本精品99久久精品77| 国产成人精品婷婷| 国产色爽女视频免费观看| 久久久久久九九精品二区国产| 97热精品久久久久久| 成人综合一区亚洲| 亚洲av.av天堂| 国产精品久久久久久久久免| 你懂的网址亚洲精品在线观看 | 亚洲av.av天堂| 国产在线精品亚洲第一网站| 欧美日韩精品成人综合77777| 久久久色成人| 国产极品精品免费视频能看的| 国产老妇伦熟女老妇高清| 变态另类成人亚洲欧美熟女| 超碰av人人做人人爽久久| 亚洲人成网站高清观看| 亚洲内射少妇av| 又爽又黄a免费视频| 少妇丰满av| 一本一本综合久久| 国产精品.久久久| 一级毛片我不卡| 国产精品1区2区在线观看.| 久久6这里有精品| 精品少妇黑人巨大在线播放 | 国产一区二区三区在线臀色熟女| 国产成人freesex在线| 深爱激情五月婷婷| 在线播放无遮挡| 综合色av麻豆| 成人美女网站在线观看视频| 久久久久久大精品| 精品久久久久久久人妻蜜臀av| 九色成人免费人妻av| 国产又黄又爽又无遮挡在线| 亚洲国产欧美人成| 亚洲人成网站在线观看播放| 成年女人永久免费观看视频| 国产色婷婷99| 99九九线精品视频在线观看视频| 国产午夜精品论理片| 久久久久久大精品| 久久久久久久久中文| 成人三级黄色视频| 国产真实乱freesex| 欧美日本亚洲视频在线播放| 国产真实乱freesex| 1000部很黄的大片| 简卡轻食公司| 亚洲国产精品sss在线观看| 国产精品伦人一区二区| 亚洲精品乱码久久久久久按摩| 人妻少妇偷人精品九色| 国产乱人视频| 非洲黑人性xxxx精品又粗又长| 午夜视频国产福利| 嫩草影院精品99| 免费一级毛片在线播放高清视频| 久久久久久国产a免费观看| 日韩一区二区视频免费看| 一边亲一边摸免费视频| 天堂影院成人在线观看| 免费不卡的大黄色大毛片视频在线观看 | 午夜精品国产一区二区电影 | 精品一区二区免费观看| 成人亚洲精品av一区二区| 亚洲精品日韩av片在线观看| 麻豆av噜噜一区二区三区| 男女做爰动态图高潮gif福利片| 女人被狂操c到高潮| 欧美激情国产日韩精品一区| av视频在线观看入口| 欧美三级亚洲精品| 久久精品国产亚洲网站| 日韩精品有码人妻一区| 国产精品电影一区二区三区| а√天堂www在线а√下载| 99视频精品全部免费 在线| 国产精品美女特级片免费视频播放器| 深夜a级毛片| 国产av一区在线观看免费| 色哟哟·www| 日韩,欧美,国产一区二区三区 | 麻豆乱淫一区二区| 精品国产三级普通话版| 一级毛片aaaaaa免费看小| 天美传媒精品一区二区| 亚洲人成网站在线播放欧美日韩| 少妇熟女欧美另类| 男女做爰动态图高潮gif福利片| 欧美最黄视频在线播放免费| а√天堂www在线а√下载| 亚洲电影在线观看av| 青春草国产在线视频 | 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 亚洲,欧美,日韩| 一区二区三区免费毛片| 波多野结衣高清无吗| 99热这里只有精品一区| 国产一区二区激情短视频| 国产精品人妻久久久久久| 日日摸夜夜添夜夜爱| av在线蜜桃| 国产精品伦人一区二区| 亚洲自偷自拍三级| 免费看光身美女| 日本欧美国产在线视频| 91av网一区二区| 成人美女网站在线观看视频| 搡老妇女老女人老熟妇| 三级国产精品欧美在线观看| 干丝袜人妻中文字幕| 美女黄网站色视频| 伦精品一区二区三区| 亚洲欧美清纯卡通| 91aial.com中文字幕在线观看| 中文字幕人妻熟人妻熟丝袜美| 狂野欧美激情性xxxx在线观看| 国产精品一区二区在线观看99 | 伊人久久精品亚洲午夜| 给我免费播放毛片高清在线观看| 国产 一区精品| 最近最新中文字幕大全电影3| 听说在线观看完整版免费高清| 精品不卡国产一区二区三区| 麻豆精品久久久久久蜜桃| 色综合亚洲欧美另类图片| 天堂av国产一区二区熟女人妻| 卡戴珊不雅视频在线播放| 国产爱豆传媒在线观看| 国产久久久一区二区三区| 成人特级av手机在线观看| 黄色日韩在线| 韩国av在线不卡| 国产成人91sexporn| 日韩欧美在线乱码| 高清毛片免费观看视频网站| 久久99精品国语久久久| 日本一本二区三区精品| 国产黄片视频在线免费观看| 99久久久亚洲精品蜜臀av| 国产v大片淫在线免费观看| АⅤ资源中文在线天堂| 大香蕉久久网| 在线免费观看不下载黄p国产| av天堂在线播放| 级片在线观看| 综合色丁香网| 12—13女人毛片做爰片一| 成年免费大片在线观看| 五月玫瑰六月丁香| 国产毛片a区久久久久| 最近2019中文字幕mv第一页| 丰满的人妻完整版| 精品久久久久久久人妻蜜臀av| 亚洲精品影视一区二区三区av| 日日啪夜夜撸| 日本一本二区三区精品| av在线播放精品| 边亲边吃奶的免费视频| 亚洲不卡免费看| 我要看日韩黄色一级片| 午夜福利在线观看免费完整高清在 | 成人亚洲欧美一区二区av| 国产精品一二三区在线看| 国产精品免费一区二区三区在线| 精品少妇黑人巨大在线播放 | 小蜜桃在线观看免费完整版高清| 最后的刺客免费高清国语| 日韩欧美精品v在线| 啦啦啦韩国在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 成人毛片a级毛片在线播放| 国产久久久一区二区三区| 黄色视频,在线免费观看| 寂寞人妻少妇视频99o| 午夜福利成人在线免费观看| 色综合站精品国产| 中文亚洲av片在线观看爽| 一级av片app| 精品久久国产蜜桃| 大香蕉久久网| 日韩一区二区三区影片| 精品午夜福利在线看| www.av在线官网国产| 国产黄色视频一区二区在线观看 | 在现免费观看毛片| 亚洲欧洲国产日韩| 国产午夜精品论理片| 亚洲国产色片| 久久这里只有精品中国| 欧美日韩国产亚洲二区| 亚洲国产精品成人久久小说 | 看黄色毛片网站| 99久久精品国产国产毛片| 午夜a级毛片| 91狼人影院| 国产老妇伦熟女老妇高清| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲电影在线观看av| 国产91av在线免费观看| 九九热线精品视视频播放| 午夜福利在线观看免费完整高清在 | 欧美又色又爽又黄视频| 成人美女网站在线观看视频| 天天躁日日操中文字幕| 人人妻人人澡欧美一区二区| 一级毛片aaaaaa免费看小| 18禁裸乳无遮挡免费网站照片| 精品久久久久久久久亚洲|