1
|
Ramachandran S, Rupakheti M, Cherian R, Lawrence MG. Aerosols heat up the Himalayan climate. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 894:164733. [PMID: 37327904 DOI: 10.1016/j.scitotenv.2023.164733] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 06/04/2023] [Accepted: 06/06/2023] [Indexed: 06/18/2023]
Abstract
The impact of aerosols, especially the absorbing aerosols, in the Himalayan region is important for climate. We closely examine ground-based high-quality observations of aerosol characteristics including radiative forcing from several locations in the Indo-Gangetic Plain (IGP), the Himalayan foothills and the Tibetan Plateau, relatively poorly studied regions with several sensitive ecosystems of global importance, as well as highly vulnerable large populations. This paper presents a state-of-the-art treatment of the warming that arises from these particles, using a combination of new measurements and modeling techniques. This is a first-time analysis of its kind, including ground-based observations, satellite data, and model simulations, which reveals that the aerosol radiative forcing efficiency (ARFE) in the atmosphere is clearly high over the IGP and the Himalayan foothills (80-135 Wm-2 per unit aerosol optical depth (AOD)), with values being greater at higher elevations. AOD is >0.30 and single scattering albedo (SSA) is ∼0.90 throughout the year over this region. The mean ARFE is 2-4 times higher here than over other polluted sites in South and East Asia, owing to higher AOD and aerosol absorption (i.e., lower SSA). Further, the observed annual mean aerosol-induced atmospheric heating rates (0.5-0.8 Kelvin/day), which are significantly higher than previously reported values for the region, imply that the aerosols alone could account for >50 % of the total warming (aerosols + greenhouse gases) of the lower atmosphere and surface over this region. We demonstrate that the current state-of-the-art models used in climate assessments significantly underestimate aerosol-induced heating, efficiency and warming over the Hindu Kush - Himalaya - Tibetan Plateau (HKHTP) region, indicating a need for a more realistic representation of aerosol properties, especially of black carbon and other aerosols. The significant, regionally coherent aerosol-induced warming that we observe in the high altitudes of the region, is a significant factor contributing to increasing air temperature, observed accelerated retreat of the glaciers, and changes in the hydrological cycle and precipitation patterns over this region. Thus, aerosols are heating up the Himalayan climate, and will remain a key factor driving climate change over the region.
Collapse
Affiliation(s)
- S Ramachandran
- Physical Research Laboratory, Ahmedabad, India; Research Institute for Sustainability - Helmholtz Centre Potsdam, Potsdam, Germany.
| | - Maheswar Rupakheti
- Research Institute for Sustainability - Helmholtz Centre Potsdam, Potsdam, Germany
| | - Ribu Cherian
- Leipzig Institute for Meteorology, University of Leipzig, Leipzig, Germany
| | - Mark G Lawrence
- Research Institute for Sustainability - Helmholtz Centre Potsdam, Potsdam, Germany; Institute for Environmental Sciences and Geography, University of Potsdam, Potsdam, Germany
| |
Collapse
|
2
|
Khamala GW, Makokha JW, Boiyo R, Kumar KR. Spatiotemporal analysis of absorbing aerosols and radiative forcing over environmentally distinct stations in East Africa during 2001-2018. THE SCIENCE OF THE TOTAL ENVIRONMENT 2023; 864:161041. [PMID: 36563754 DOI: 10.1016/j.scitotenv.2022.161041] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 12/11/2022] [Accepted: 12/15/2022] [Indexed: 06/17/2023]
Abstract
East Africa (EA) suffers from the inadequate characterization of atmospheric aerosols, with far-reaching consequences of its inability to quantify precisely the impacts of these particles on regional climate. The current study aimed at characterizing absorption and radiative properties of aerosols using the long-term (2001-2018) AErosol RObotic NETwork (AERONET) and Modern-Era Retrospective analysis for Research and Applications (MERRA-2) data over three environmentally specific sites in EA. The annual mean absorption aerosol optical depth (AAOD440 nm), absorption Angstrom Exponent (AAE440-870 nm), total effective radius (REff), and total volume concentration (μm3/μm2) revealed significant spatial heterogeneity over the domain. The study domain exhibited a significant contribution of fine-mode aerosols compared to the coarse-mode particles. The monthly variation in SSA440 nm over EA explains the strength in absorption aerosols that range from moderate to strong absorbing aerosols. The aerosols exhibited significant variability over the study domain, with the dominance of absorbing fine-mode aerosols over Mbita accounting for ∼40 to ∼50 %, while weakly absorbing coarse-mode particles accounted for ∼8.2 % over Malindi. The study conclusively determined that Mbita was dominated by AAOD mainly from biomass burning in most of the months, whereas Malindi was coated with black carbon. The direct aerosol radiative forcing (DARF) retrieved from both the AERONET and MERRA-2 models showed strong cooling at the top of the atmosphere (TOA; -6 to -27 Wm-2) and the bottom of the atmosphere (BOA, -7 to -66 Wm-2). However, significant warming was noticed within the atmosphere (ATM; +14 to +76 Wm-2), an indication of the role of aerosols in regional climate change. The study contributed to understanding aerosol absorption and radiative characteristics over EA and can form the basis of other related studies over the domain and beyond.
Collapse
Affiliation(s)
- Geoffrey W Khamala
- Department of Science Technology and Engineering, Kibabii University, P.O. Box 1699-50200, Bungoma, Kenya.
| | - John W Makokha
- Department of Science Technology and Engineering, Kibabii University, P.O. Box 1699-50200, Bungoma, Kenya
| | - Richard Boiyo
- Department of Physical Sciences, Meru University of Science and Technology, P.O. Box 972-60200, Meru, Kenya; Department of Environment, Water, Energy and Natural Resources, County Government of Vihiga, Maragoli, Kenya
| | - Kanike Raghavendra Kumar
- Department of Engineering Physics, College of Engineering, Koneru Lakshmaiah Education Foundation, Vaddeswaram, 522 302 Guntur, Andhra Pradesh, India
| |
Collapse
|
3
|
Ramachandran S, Rupakheti M, Lawrence MG. Aerosol-induced atmospheric heating rate decreases over South and East Asia as a result of changing content and composition. Sci Rep 2020; 10:20091. [PMID: 33208825 PMCID: PMC7676243 DOI: 10.1038/s41598-020-76936-z] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2020] [Accepted: 10/28/2020] [Indexed: 11/22/2022] Open
Abstract
Aerosol emissions from human activities are extensive and changing rapidly over Asia. Model simulations and satellite observations indicate a dipole pattern in aerosol emissions and loading between South Asia and East Asia, two of the most heavily polluted regions of the world. We examine the previously unexplored diverging trends in the existing dipole pattern of aerosols between East and South Asia using the high quality, two-decade long ground-based time series of observations of aerosol properties from the Aerosol Robotic Network (AERONET), from satellites (Moderate Resolution Imaging Spectroradiometer (MODIS) and Ozone Monitoring Instrument (OMI)), and from model simulations (Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). The data cover the period since 2001 for Kanpur (South Asia) and Beijing (East Asia), two locations taken as being broadly representative of the respective regions. Since 2010 a dipole in aerosol optical depth (AOD) is maintained, but the trend is reversed—the decrease in AOD over Beijing (East Asia) is rapid since 2010, being 17% less in current decade compared to first decade of twenty-first century, while the AOD over South Asia increased by 12% during the same period. Furthermore, we find that the aerosol composition is also changing over time. The single scattering albedo (SSA), a measure of aerosol’s absorption capacity and related to aerosol composition, is slightly higher over Beijing than Kanpur, and has increased from 0.91 in 2002 to 0.93 in 2017 over Beijing and from 0.89 to 0.92 during the same period over Kanpur, confirming that aerosols in this region have on an average become more scattering in nature. These changes have led to a notable decrease in aerosol-induced atmospheric heating rate (HR) over both regions between the two decades, decreasing considerably more over East Asia (− 31%) than over South Asia (− 9%). The annual mean HR is lower now, it is still large (≥ 0.6 K per day), which has significant climate implications. The seasonal trends in AOD, SSA and HR are more pronounced than their respective annual trends over both regions. The seasonal trends are caused mainly by the increase/decrease in anthropogenic aerosol emissions (sulfate, black carbon and organic carbon) while the natural aerosols (dust and sea salt) did not change significantly over South and East Asia during the last two decades. The MERRA-2 model is able to simulate the observed trends in AODs well but not the magnitude, while it also did not simulate the SSA values or trends well. These robust findings based on observations of key aerosol parameters and previously unrecognized diverging trends over South and East Asia need to be accounted for in current state-of-the-art climate models to ensure accurate quantification of the complex and evolving impact of aerosols on the regional climate over Asia.
Collapse
Affiliation(s)
- S Ramachandran
- Physical Research Laboratory, Ahmedabad, India. .,Institute for Advanced Sustainability Studies, Potsdam, Germany.
| | | | - Mark G Lawrence
- Institute for Advanced Sustainability Studies, Potsdam, Germany.,Institute for Environmental Sciences and Geography, University of Potsdam, Potsdam, Germany
| |
Collapse
|
4
|
Just AC, Liu Y, Sorek-Hamer M, Rush J, Dorman M, Chatfield R, Wang Y, Lyapustin A, Kloog I. Gradient boosting machine learning to improve satellite-derived column water vapor measurement error. ATMOSPHERIC MEASUREMENT TECHNIQUES 2020; 13:4669-4681. [PMID: 33193906 PMCID: PMC7665162 DOI: 10.5194/amt-13-4669-2020] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
The atmospheric products of the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm include column water vapor (CWV) at a 1 km resolution, derived from daily overpasses of NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) instruments aboard the Aqua and Terra satellites. We have recently shown that machine learning using extreme gradient boosting (XGBoost) can improve the estimation of MAIAC aerosol optical depth (AOD). Although MAIAC CWV is generally well validated (Pearson's R >0.97 versus CWV from AERONET sun photometers), it has not yet been assessed whether machine-learning approaches can further improve CWV. Using a novel spatiotemporal cross-validation approach to avoid overfitting, our XGBoost model, with nine features derived from land use terms, date, and ancillary variables from the MAIAC retrieval, quantifies and can correct a substantial portion of measurement error relative to collocated measurements at AERONET sites (26.9% and 16.5% decrease in root mean square error (RMSE) for Terra and Aqua datasets, respectively) in the Northeastern USA, 2000-2015. We use machine-learning interpretation tools to illustrate complex patterns of measurement error and describe a positive bias in MAIAC Terra CWV worsening in recent summertime conditions. We validate our predictive model on MAIAC CWV estimates at independent stations from the SuomiNet GPS network where our corrections decrease the RMSE by 19.7% and 9.5% for Terra and Aqua MAIAC CWV. Empirically correcting for measurement error with machine-learning algorithms is a postprocessing opportunity to improve satellite-derived CWV data for Earth science and remote sensing applications.
Collapse
Affiliation(s)
- Allan C. Just
- Department of Environmental Medicine and Public Health, Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Yang Liu
- Department of Environmental Medicine and Public Health, Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Meytar Sorek-Hamer
- Universities Space Research Association (USRA), Mountain View, California, USA
- NASA Ames Research Center, Mountain View, California, USA
| | - Johnathan Rush
- Department of Environmental Medicine and Public Health, Icahn School of Medicine at Mount Sinai, New York, New York, USA
| | - Michael Dorman
- Department of Geography and Environmental Development, Ben-Gurion University of the Negev, Beersheba, Israel
| | | | - Yujie Wang
- Joint Center for Earth Systems Technology, University of Maryland, Baltimore County, Baltimore, Maryland, USA
- NASA Goddard Space Flight Center, Greenbelt, Maryland, USA
| | | | - Itai Kloog
- Department of Geography and Environmental Development, Ben-Gurion University of the Negev, Beersheba, Israel
| |
Collapse
|
5
|
Franzin BT, Guizellini FC, de Babos DV, Hojo O, Pastre IA, Marchi MRR, Fertonani FL, Oliveira CMRR. Characterization of atmospheric aerosol (PM 10 and PM 2.5) from a medium sized city in São Paulo state, Brazil. J Environ Sci (China) 2020; 89:238-251. [PMID: 31892396 DOI: 10.1016/j.jes.2019.09.014] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2019] [Revised: 09/07/2019] [Accepted: 09/11/2019] [Indexed: 06/10/2023]
Abstract
Air pollution causes deleterious effects on human health with aerosols being among the most polluting agents. The objective of this work is the characterization of the PM2.5 and PM10 aerosol mass in the atmosphere. The methods of analysis include WD-XRF and EDS. Data were correlated with meteorological information and air mass trajectories (model HYSPLIT) by multivariate analysis. A morphological structural analysis was also carried out to identify the probable sources of atmospheric aerosols in the city of São José do Rio Preto, Brazil. The mean mass concentration values obtained were 24.54 μg/m3 for PM10, above the WHO annual standard value of 20 μg/m3 and 10.88 μg/m3 for PM2.5 whose WHO recommended limit is 10 μg/m3. WD-XRF analysis of the samples revealed Si and Al as major components of the coarse fraction. In the fine fraction, the major elements were Al and S. The SEM-FEG characterization allowed identifying the morphology of the particles in agglomerates, ellipsoids and filaments in the PM10, besides spherical in the PM2.5. The analysis by EDS corroborated WD-XRF results, identifying the crustal elements, aluminosilicates and elements of anthropogenic origin in the coarse fraction. For the fine fraction crustal elements were also identified; aluminosilicates, black carbon and spherical particles (C and O) originating from combustion processes were predominant. The use of multivariate analysis to correlate air mass trajectories with the results of the morpho-structural characterization of the particulate matter allowed confirmation of the complex composition of the particles resulting from the combination of both local and long-distance sources.
Collapse
Affiliation(s)
- Bruno T Franzin
- São Paulo State University (Unesp), Chemistry Institute, Analytical Chemistry Department, Araraquara, SP, 14800-060, Brazil; Centro de Química Estrutural - Faculdade de Ciências da Universidade de Lisboa, Edif. C8, Campo Grande, 1749-016 Lisboa, Portugal.
| | - Filipe C Guizellini
- São Paulo State University (Unesp), Chemistry Institute, Analytical Chemistry Department, Araraquara, SP, 14800-060, Brazil
| | - Diego V de Babos
- Federal University of São Carlos - UFSCAr, Chemistry Department, São Carlos, SP, 13565-905, Brazil
| | - Ossamu Hojo
- São Paulo State University (Unesp), Chemistry Institute, Analytical Chemistry Department, Araraquara, SP, 14800-060, Brazil
| | - Iêda Ap Pastre
- Chemistry and Environmental Sciences Department, Biosciences, Languages and Exact Sciences Institute, Ibilce, São Paulo State University (Unesp), São José do Rio Preto, SP, 15054-000, Brazil
| | - Mary R R Marchi
- São Paulo State University (Unesp), Chemistry Institute, Analytical Chemistry Department, Araraquara, SP, 14800-060, Brazil
| | - Fernando L Fertonani
- São Paulo State University (Unesp), Chemistry Institute, Analytical Chemistry Department, Araraquara, SP, 14800-060, Brazil; Chemistry and Environmental Sciences Department, Biosciences, Languages and Exact Sciences Institute, Ibilce, São Paulo State University (Unesp), São José do Rio Preto, SP, 15054-000, Brazil
| | - Cristina M R R Oliveira
- Centro de Química Estrutural - Faculdade de Ciências da Universidade de Lisboa, Edif. C8, Campo Grande, 1749-016 Lisboa, Portugal
| |
Collapse
|
6
|
Wang H, Li Z, Lv Y, Xu H, Li K, Li D, Hou W, Zheng F, Wei Y, Ge B. Observational study of aerosol-induced impact on planetary boundary layer based on lidar and sunphotometer in Beijing. ENVIRONMENTAL POLLUTION (BARKING, ESSEX : 1987) 2019; 252:897-906. [PMID: 31212251 DOI: 10.1016/j.envpol.2019.05.070] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/28/2018] [Revised: 05/13/2019] [Accepted: 05/13/2019] [Indexed: 06/09/2023]
Abstract
Atmospheric aerosols have been found to influence the development of planetary boundary layer (PBL) and hence to aggravate haze pollution in megacities. PBL height (PBLH) determines the vertical extent to which the most pollutant effectively disperses and is a key argument in pollution study. In this study, we quantitatively evaluate aerosol radiation effect on PBL, as well as assessment of surface cooling effect and atmosphere heating effect. All the data are measured at a site of Beijing from 2014 to 2017, of which PBLH is retrieved from micro pulse lidar and aerosol optical depth (AOD) from sunphotometer. Case study shows qualitatively that relative high aerosol load reduces PBLH, and in turn causes a high surface PM2.5 concentration. We preliminarily reveal the influential mechanism of aerosol on PBL. The influence of aerosol on the radiation flux of PBL is analyzed, with the correlation coefficient (R) of 0.938 between AOD and radiative forcing of BOA (RFBOA) and R = 0.43 between RFBOA and PBLH. Also, AOD is found to negatively correlate with PBLH (R = -0.41). With the increase of AOD, the cooling effect of surface is enhanced, and further impede the development of PBL. Due to aerosol-induced reduction of PBLH, near surface PM2.5 concentration surges and presents an exponential growth following AOD. Then, it is speculated and testified that the relationship between SSA (single scatting albedo) and PBLH would be determined by the location of absorbing aerosol within PBL. The upper PBL absorbing aerosol may decrease PBLH, while the lower absorbing aerosol appear to enhance PBLH. The study probably can provide effective observational evidence for understanding the effect of aerosol on PBL and be a reference of air pollution mitigation in Beijing and its surrounding areas.
Collapse
Affiliation(s)
- Haofei Wang
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China; University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Zhengqiang Li
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China.
| | - Yang Lv
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China; University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Hua Xu
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China
| | - Kaitao Li
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China
| | - Donghui Li
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China
| | - Weizhen Hou
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China
| | - Fengxun Zheng
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China; University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Yuanyuan Wei
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China; University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Bangyu Ge
- State Environment Protection Key Laboratory of Satellite Remote Sensing, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing, 100101, China; University of Chinese Academy of Sciences, Beijing, 100101, China
| |
Collapse
|
7
|
Rupakheti D, Kang S, Rupakheti M, Cong Z, Tripathee L, Panday AK, Holben BN. Observation of optical properties and sources of aerosols at Buddha's birthplace, Lumbini, Nepal: environmental implications. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2018; 25:14868-14881. [PMID: 29546514 DOI: 10.1007/s11356-018-1713-z] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Accepted: 03/08/2018] [Indexed: 06/08/2023]
Abstract
For the first time, aerosol optical properties are measured over Lumbini, Nepal, with CIMEL sunphotometer of the Aerosol Robotic Network (AERONET) program. Lumbini is a sacred place as the birthplace of Lord Buddha, and thus a UNESCO world heritage site, located near the northern edge of the central Indo-Gangetic Plains (IGP) and before the Himalayan foothills (and Himalayas) to its north. Average aerosol optical depth (AOD) is found to be 0.64 ± 0.38 (0.06-3.28) over the sampling period (January 2013-December 2014), with the highest seasonal AOD during the post-monsoon season (0.72 ± 0.44). More than 80% of the daily averaged AOD values, during the monitoring period, are above 0.3, indicating polluted conditions in the region. The levels of aerosol load observed over Lumbini are comparable to those observed at several heavily polluted sites in the IGP. Based on the relationship between AOD and Ångstrom exponent (α), anthropogenic, biomass burning, and mixed aerosols are found to be the most prevalent aerosol types. The aerosol volume-size distribution is bi-modal during all four seasons with modes centered at 0.1-0.3 and 3-4 μm. For both fine and coarse modes, the highest volumetric concentration of ~ 0.08 μm-3 μm-2 is observed during the post-monsoon and pre-monsoon seasons. As revealed by the single-scattering albedo (SSA), asymmetry parameter (AP), and refractive index (RI) analyses, aerosol loading over Lumbini is dominated by absorbing, urban-industrial, and biomass burning aerosols.
Collapse
Affiliation(s)
- Dipesh Rupakheti
- Key Laboratory of Tibetan Environment Changes and Land Surface Processes, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing, 100101, China.
- University of Chinese Academy of Sciences, Beijing, 100049, China.
- State Key Laboratory of Cryospheric Science, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, 730000, China.
| | - Shichang Kang
- University of Chinese Academy of Sciences, Beijing, 100049, China.
- State Key Laboratory of Cryospheric Science, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, 730000, China.
- Center for Excellence in Tibetan Plateau Earth Sciences, Chinese Academy of Sciences, Beijing, 100101, China.
| | - Maheswar Rupakheti
- Institute for Advanced Sustainability Studies (IASS), Potsdam, Germany
- Himalayan Sustainability Institute (HIMSI), Kathmandu, Nepal
| | - Zhiyuan Cong
- Key Laboratory of Tibetan Environment Changes and Land Surface Processes, Institute of Tibetan Plateau Research, Chinese Academy of Sciences, Beijing, 100101, China
- Center for Excellence in Tibetan Plateau Earth Sciences, Chinese Academy of Sciences, Beijing, 100101, China
| | - Lekhendra Tripathee
- State Key Laboratory of Cryospheric Science, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences, Lanzhou, 730000, China
| | - Arnico K Panday
- International Centre for Integrated Mountain Development (ICIMOD), Lalitpur, Nepal
| | | |
Collapse
|
8
|
Haze Optical Properties from Long-Term Ground-Based Remote Sensing over Beijing and Xuzhou, China. REMOTE SENSING 2018. [DOI: 10.3390/rs10040518] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Aerosol haze pollution has had a significant impact on both global climate and the regional air quality of Eastern China, which has a high proportion of high level pollution days. Statistical analyses of aerosol optical properties and direct radiative forcing at two AERONET sites (Beijing and Xuzhou) were conducted from 2013 to 2016. Results indicate: (1) Haze pollution days accounted for 26% and 20% of days from 2013 to 2016 in Beijing and Xuzhou, respectively, with the highest proportions in winter; (2) The averaged aerosol optical depth (AOD) at 550 nm on haze days were about 3.7 and 1.6 times greater than those on clean days in Beijing and Xuzhou, respectively. At both sites, the maximum AOD occurred in summer; (3) Hazes were dominated by fine particles at both sites. However, as compared to Xuzhou, Beijing had larger coarse mode AOD and higher percentage of small α. This data, together with an analysis of size distribution, suggests that the hazes in Beijing were more susceptible to coarse dust particles than Xuzhou; (4) During hazes in Beijing, the single scattering albedo (SSA) is significantly higher when compared to clean conditions (0.874 vs. 0.843 in SSA440 nm), an increase much less evident in Xuzhou. The most noticeable differences in both SSA and the imaginary part of the complex refractive index between Beijing and Xuzhou were found in winter; (5) In Beijing, the haze radiative forcing produced an averaged cooling effect of −113.6 ± 63.7 W/m2 at the surface, whereas the averaged heating effect of 77.5 ± 49.7 W/m2 within the atmosphere was at least twice as strong as clean days. In Xuzhou, such a radiative forcing effect appeared to be much smaller and the difference between haze and clean days was insignificant. Derived from long-term observation, these findings are more significant for the improvement of our understanding of haze formation in China and the assessment of its impacts on radiative forcing of climate change than previous short-term case studies.
Collapse
|
9
|
Song SK, Shon ZH, Park YH. Diurnal and seasonal characteristics of the optical properties and direct radiative forcing of different aerosol components in Seoul megacity. THE SCIENCE OF THE TOTAL ENVIRONMENT 2017; 599-600:400-412. [PMID: 28482298 DOI: 10.1016/j.scitotenv.2017.04.195] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Revised: 04/26/2017] [Accepted: 04/26/2017] [Indexed: 05/25/2023]
Abstract
The temporal variations (diurnal and seasonal) of the optical properties and direct aerosol radiative forcing (DARF) of different aerosol components (water-soluble, insoluble, black carbon (BC), and sea-salt) were analyzed using the hourly resolution data (PM2.5) measured at an urban site in Seoul, Korea during 2010, based on a modeling approach. In general, the water-soluble component was predominant over all other components (with a higher concentration) in terms of its impact on the optical properties (except for absorbing BC) and DARF. The annual mean aerosol optical depth (AOD, τ) at 500nm for the water-soluble component was 0.38±0.07 (0.06±0.01 for BC). The forcing at the surface (DARFSFC) and top of the atmosphere (DARFTOA), and in the atmosphere (DARFATM) for most aerosol components (except for BC) during the daytime were highest in spring and lowest in late fall or early winter. The maximum DARFSFC occurred in the morning during most seasons (except for the water-soluble components showing peaks in the afternoon or noon in summer, fall, or winter), while the maximum DARFTOA occurred in the morning during spring and/or winter and in the afternoon during summer and/or fall. The estimated DARFSFC and DARFATM of the water-soluble component were in the range of -49 to -84Wm-2 and +10 to +22Wm-2, respectively. The DARFSFC and DARFATM of BC were -26 to -39Wm-2 and +32 to +51Wm-2, respectively, showing highest in summer and lowest in spring, with morning peaks regardless of the season. This positive DARFATM of BC in this study area accounted for approximately 64% of the total atmospheric aerosol forcing due to strong radiative absorption, thus increasing atmospheric heating by 2.9±1.2Kday-1 (heating rate efficiency of 39K day-1τ-1) and then causing further atmospheric warming.
Collapse
Affiliation(s)
- Sang-Keun Song
- Department of Earth and Marine Sciences, Jeju National University, Jeju 63243, Republic of Korea
| | - Zang-Ho Shon
- Department of Environmental Engineering, Dong-Eui University, Busan 47340, Republic of Korea.
| | - Yeon-Hee Park
- Department of Earth and Marine Sciences, Jeju National University, Jeju 63243, Republic of Korea
| |
Collapse
|
10
|
Adesina AJ, Piketh S, Kanike RK, Venkataraman S. Characteristics of columnar aerosol optical and microphysical properties retrieved from the sun photometer and its impact on radiative forcing over Skukuza (South Africa) during 1999-2010. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2017; 24:16160-16171. [PMID: 28537035 DOI: 10.1007/s11356-017-9211-2] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/03/2017] [Accepted: 05/07/2017] [Indexed: 06/07/2023]
Abstract
The detailed analysis of columnar optical and microphysical properties of aerosols obtained from the AErosol RObotic NETwork (AERONET) Cimel sun photometer operated at Skukuza (24.98° S, 31.60° E, 150 m above sea level), South Africa was carried out using the level 2.0 direct sun and inversion products measured during 1999-2010. The observed aerosol optical depth (AOD) was generally low over the region, with high values noted in late winter (August) and mid-spring (September and October) seasons. The major aerosol types found during the study period were made of 3.74, 69.63, 9.34, 8.83, and 8.41% for polluted dust (PD), polluted continental (PC), non-absorbing (NA), slightly absorbing (SA), and moderately absorbing (MA) aerosols, respectively. Much attention was given to the aerosol fine- and coarse-modes deduced from the particle volume concentration, effective radius, and fine-mode volume fraction. The aerosol volume size distribution pattern was found to be bimodal with the fine-mode showing predominance relative to coarse-mode during the winter and spring seasons, owing to the onset of the biomass burning season. The mean values of total, fine-, and coarse-mode volume particle concentrations were 0.07 ± 0.04, 0.03 ± 0.03, and 0.04 ± 0.02 μm3 μm-2, respectively, whereas the mean respective effective radii observed at Skukuza for the abovementioned modes were 0.35 ± 0.17, 0.14 ± 0.02, and 2.08 ± 0.02 μm. The averaged shortwave direct aerosol radiative forcing (ARF) observed within the atmosphere was found to be positive (absorption or heating effect), whereas the negative forcing in the surface and TOA depicted significant cooling effect due to more scattering type particles.
Collapse
Affiliation(s)
- Ayodele Joseph Adesina
- School of Geo- and Spatial Science, Unit for Environmental Sciences and Management, North-West University, Potchefstroom, 2520, South Africa
| | - Stuart Piketh
- School of Geo- and Spatial Science, Unit for Environmental Sciences and Management, North-West University, Potchefstroom, 2520, South Africa
| | - Raghavendra Kumar Kanike
- Key Laboratory of Meteorological Disasters, Ministry of Education (KLME), Collaborative Innovation Centre on Forecast and Evaluation of Meteorological Disasters, Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing, Jiangsu, 210044, China.
| | - Sivakumar Venkataraman
- Discipline of Physics, School of Chemistry and Physics, University of KwaZulu-Natal, Durban, Kwazulu-Natal, 4000, South Africa
| |
Collapse
|
11
|
Mor V, Dhankhar R, Attri SD, Soni VK, Sateesh M, Taneja K. Assessment of aerosols optical properties and radiative forcing over an Urban site in North-Western India. ENVIRONMENTAL TECHNOLOGY 2017; 38:1232-1244. [PMID: 27564392 DOI: 10.1080/09593330.2016.1221473] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/04/2016] [Accepted: 08/02/2016] [Indexed: 06/06/2023]
Abstract
The present work is aimed to analyze aerosols optical properties and to estimate aerosol radiative forcing (ARF) from January to December 2013, using sky radiometer data over Rohtak, an urban site in North-Western India. The results reveal strong wavelength dependency of aerosol optical depth (AOD), with high values of AOD at shorter wavelengths and lower values at longer wavelength during the study period. The highest AOD values of 1.07 ± 0.45 at 500 nm were observed during July. A significant decline in Ångström exponent was observed during April-May, which represents the dominance of coarse mode particles due to dust-raising convective activities. Aerosols' size distribution exhibits a bimodal structure with fine mode particles around 0.17 µm and coarse mode particles with a radius around 5.28 µm. Single scattering albedo values were lowest during November-December at all wavelengths, ranging from 0.87 to 0.76, which corresponds to the higher absorption during this period. Aerosols optical properties retrieved during observation period are used as input for SBDART (Santa Barbara DISORT Atmospheric Radiative Transfer) to estimate the direct ARF at the surface, in the atmosphere and at the top of the atmosphere (TOA). The ARF at the TOA, surface and in the atmosphere are found to be in the range of -4.98 to -19.35 W m-2, -8.01 to -57.66 W m-2 and +3.02 to +41.64 W m-2, respectively. The averaged forcing for the whole period of observations at the TOA is -11.26 W m-2, while at the surface it is -38.64 W m-2, leading to atmospheric forcing of 27.38 W m-2. The highest (1.168 K day-1) values of heating rate was estimated during November, whereas the lowest value (0.084 K day-1) was estimated for the February.
Collapse
Affiliation(s)
- Vikram Mor
- a Department of Environment Science , Maharshi Dayanand University , Rohtak , India
| | - Rajesh Dhankhar
- a Department of Environment Science , Maharshi Dayanand University , Rohtak , India
| | - S D Attri
- b India Meteorological Department , Ministry of Earth Sciences , New Delhi , India
| | - V K Soni
- b India Meteorological Department , Ministry of Earth Sciences , New Delhi , India
| | - M Sateesh
- b India Meteorological Department , Ministry of Earth Sciences , New Delhi , India
| | - Kanika Taneja
- b India Meteorological Department , Ministry of Earth Sciences , New Delhi , India
| |
Collapse
|
12
|
Kang N, Kumar KR, Yu X, Yin Y. Column-integrated aerosol optical properties and direct radiative forcing over the urban-industrial megacity Nanjing in the Yangtze River Delta, China. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2016; 23:17532-17552. [PMID: 27234827 DOI: 10.1007/s11356-016-6953-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2016] [Accepted: 05/20/2016] [Indexed: 06/05/2023]
Abstract
Aerosol optical properties were measured and analyzed through the ground-based remote sensing Aerosol Robotic Network (AERONET) over an urban-industrial site, Nanjing (32.21° N, 118.72° E, and 62 m above sea level), in the Yangtze River Delta, China, during September 2007-August 2008. The annual averaged values of aerosol optical depth (AOD500) and the Ångström exponent (AE440-870) were measured to be 0.94 ± 0.52 and 1.10 ± 0.21, respectively. The seasonal averaged values of AOD500 (AE440-870) were noticed to be high in summer (autumn) and low in autumn (spring). The characterization of aerosol types showed the dominance of mixed type followed by the biomass burning and urban-industrial type of aerosol at Nanjing. Subsequently, the curvature (a 2) obtained from the second-order polynomial fit and the second derivative of AE (α') were also analyzed to understand the dominant aerosol type. The single scattering albedo at 440 nm (SSA440) varied from 0.88 to 0.93 with relatively lower (higher) values during the summer (spring), suggesting an increase in black carbon and mineral dust (desert dust) aerosols of absorbing (scattering) nature. The averaged monthly and seasonal evolutions of shortwave (0.3-4.0 μm) direct aerosol radiative forcing (DARF) values were computed from the Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART) model both at the top of atmosphere (TOA) and bottom of atmosphere (SUR) during the study period. Further, the aerosol forcing efficiency (AFE) and the corresponding atmospheric heating rates (AHR) were also estimated from the forcing within the atmosphere (ATM). The derived DARF values, therefore, produced a warming effect within the atmosphere due to strong absorption of solar radiation.
Collapse
Affiliation(s)
- Na Kang
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/International Joint Laboratory on Climate and Environmental Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing, 210044, Jiangsu, China
| | - K Raghavendra Kumar
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/International Joint Laboratory on Climate and Environmental Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing, 210044, Jiangsu, China.
| | - Xingna Yu
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/International Joint Laboratory on Climate and Environmental Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing, 210044, Jiangsu, China
| | - Yan Yin
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/International Joint Laboratory on Climate and Environmental Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing, 210044, Jiangsu, China
| |
Collapse
|
13
|
Yu X, Lü R, Kumar KR, Ma J, Zhang Q, Jiang Y, Kang N, Yang S, Wang J, Li M. Dust aerosol properties and radiative forcing observed in spring during 2001-2014 over urban Beijing, China. ENVIRONMENTAL SCIENCE AND POLLUTION RESEARCH INTERNATIONAL 2016; 23:15432-15442. [PMID: 27117151 DOI: 10.1007/s11356-016-6727-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/19/2015] [Accepted: 04/19/2016] [Indexed: 06/05/2023]
Abstract
The ground-based characteristics (optical and radiative properties) of dust aerosols measured during the springtime between 2001 and 2014 were investigated over urban Beijing, China. The seasonal averaged aerosol optical depth (AOD) during spring of 2001-2014 was about 0.78 at 440 nm. During dust days, higher AOD occurred associated with lower Ångström exponent (AE). The mean AE440-870 in the springtime was about 1.0, indicating dominance of fine particles over the region. The back-trajectory analysis revealed that the dust was transported from the deserts of Inner Mongolia and Mongolia arid regions to Beijing. The aerosol volume size distribution showed a bimodal distribution pattern, with its highest peak observed in coarse mode for all episodes (especially for dust days with increased volume concentration). The single scattering albedo (SSA) increased with wavelength on dust days, indicating the presence of more scattering particles. Furthermore, the complex parts (real and imaginary) of refractive index showed distinct characteristics with lower imaginary values (also scattering) on dust days. The shortwave (SW; 0.2-4.0 μm) and longwave (LW; 4-100 μm) aerosol radiative forcing (ARF) values were computed from the Santa Barbara DISORT Atmospheric Radiative Transfer (SBDART) model both at the top of atmosphere (TOA) and the bottom of atmosphere (BOA) during dust and non-dust (dust free) days, and the corresponding heating rates and forcing efficiencies were also estimated. The SW (LW) ARF, therefore, produced significant cooling (warming) effects at both the TOA and the BOA over Beijing.
Collapse
Affiliation(s)
- Xingna Yu
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China
| | - Rui Lü
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China
| | - K Raghavendra Kumar
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China.
| | - Jia Ma
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China
- Guangzhou Hexin Analytical Instrument Company Limited, Guangzhou, 510530, China
| | - Qiuju Zhang
- Baoshan District Center for Disease Control and Prevention, Shanghai, 201901, China
| | - Yilun Jiang
- Guangzhou Hexin Analytical Instrument Company Limited, Guangzhou, 510530, China
| | - Na Kang
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China
| | - Suying Yang
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China
| | - Jing Wang
- Key Laboratory of Meteorological Disaster, Ministry of Education (KLME)/Joint International Research Laboratory of Climate and Environment Change (ILCEC)/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD)/Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing, 210044, China
| | - Mei Li
- Institute of Atmospheric Environment Safety and Pollution Control, Jinan University, Guangzhou, 510632, China.
| |
Collapse
|
14
|
Yu X, Kumar KR, Lü R, Ma J. Changes in column aerosol optical properties during extreme haze-fog episodes in January 2013 over urban Beijing. ENVIRONMENTAL POLLUTION (BARKING, ESSEX : 1987) 2016; 210:217-226. [PMID: 26735167 DOI: 10.1016/j.envpol.2015.12.021] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/15/2015] [Revised: 12/09/2015] [Accepted: 12/10/2015] [Indexed: 06/05/2023]
Abstract
Several dense haze-fog (HF) episodes were occurred in the North China Plain (NCP), especially over Beijing in January 2013 characterized by a long duration, a large influential region, and an extremely high PM2.5 values (>500 μg m(-3)). In this study, we present the characteristics of aerosol optical properties and radiative forcing using Cimel sun-sky radiometer measurements during HF and no haze-fog (NHF) episodes occurred over Beijing during 1-31 January, 2013. The respective maximum values of daily mean aerosol optical depth at 440 nm (AOD440) were observed to be 1.21, 1.43, 1.52, and 2.21 occurred on 12, 14 19, and 28 January. It was found that the Ångström exponent (AE) values were almost higher than 1.0 during all the days with its maximum on 26 January (1.53), suggests the dominance of fine-mode particles. The maximum (minimum) aerosol volume size distributions occurred during dense HF (NHF) days with larger particle volumes of fine-mode. The single scattering albedo, asymmetry parameter, and complex refractive index values during HF events suggest the abundance of fine-mode particles from anthropogenic (absorbing) activities mixed with scattering dust particles. The average shortwave direct aerosol radiative forcing (DARF) values at the bottom-of-atmosphere (BOA) during HF and NHF days were estimated to be 112.29 ± 42.18 W m(-2) and -58.61 ± 13.09 W m(-2), while at the top-of-atmosphere (TOA) the forcing values were -45.78 ± 22.17 W m(-2) and -18.64 ± 5.84 W m(-2), with the corresponding heating rate of 1.61 ± 0.48 K day(-1) and 1.12 ± 0.31 K day(-1), respectively. The DARF values retrieved from the AERONET were in good agreement with the SBDART computed both at the TOA (r = 0.95) and the BOA (r = 0.97) over Beijing in January 2013.
Collapse
Affiliation(s)
- Xingna Yu
- Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing 210044, Jiangsu province, China
| | - K Raghavendra Kumar
- Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing 210044, Jiangsu province, China.
| | - Rui Lü
- Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing 210044, Jiangsu province, China
| | - Jia Ma
- Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, School of Atmospheric Physics, Nanjing University of Information Science and Technology, Nanjing 210044, Jiangsu province, China
| |
Collapse
|