Next Article in Journal
Quantifying the Evapotranspiration Rate and Its Cooling Effects of Urban Hedges Based on Three-Temperature Model and Infrared Remote Sensing
Next Article in Special Issue
Correction: Yahui, Che et al. Validation of Aerosol Products from AATSR and MERIS/AATSR Synergy Algorithms—Part 1: Global Evaluation. Remote Sensing 2018, 10, 1414
Previous Article in Journal
Combination Analysis of Future Polar-Type Gravity Mission and GRACE Follow-On
Previous Article in Special Issue
Retrieval of Vertical Mass Concentration Distributions—Vipava Valley Case Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Time Scale Analysis of Regional Aerosol Optical Depth Changes in National-Level Urban Agglomerations in China Using Modis Collection 6.1 Datasets from 2001 to 2017

1
School of Geodesy and Geomatics, Wuhan University, Wuhan 430079, China
2
School of Resources and Environmental Engineering, Wuhan University of Technology, Wuhan 430070, China
*
Author to whom correspondence should be addressed.
Submission received: 13 December 2018 / Revised: 17 January 2019 / Accepted: 18 January 2019 / Published: 21 January 2019
(This article belongs to the Special Issue Remote Sensing of Air Quality)

Abstract

:
With the rapid development of China’s economy and industry, characterizing the spatial and temporal changes of aerosols in China has attracted widespread attention from researchers. The national-level urban agglomerations are the most concentrated areas of China’s economic, population and resource. Studying the spatial and temporal changes of aerosol optical depth (AOD) in these regions has practical guiding significance for effective monitoring of atmospheric particulate pollution. This paper analyzed the spatial and temporal variations of AOD in China’s urban agglomerations during 2001–2017 by using Terra Moderate resolution Imaging Spectroradiometer (MODIS) Collection 6.1 (C6.1) Level 2 aerosol products (MOD04_L2). Five national-level urban agglomerations were chosen: Yangtze River Delta (YRD), Pearl River Delta (PRD), Beijing-Tianjin-Hebei (BTH), Yangtze River Middle-Reach (YRMR) and Cheng-Yu (CY). We analyzed the change patterns of AOD in different urban agglomerations at multi-time scales and built a time series decomposition model to mine the long-term trend, seasonal variation and abnormal change information of AOD time series. The result indicated that averaged AOD values in the five urban agglomerations were basically increased first and then decreased at the annual time scale during 2001–2017. The averaged AOD showed strong seasonal differences and AOD values in spring and summer were typically higher than those in autumn and winter. At the monthly time scale, the AOD typically varied from low in cold months to high in warm months and then decreased during the rainy periods. Time series decompositions revealed that a notable transition around 2007–2008 dominated the long-term overall trend over the five selected urban agglomerations and an initial upward tendency followed by a downward tendency was observed during 2001–2017. This study can be utilized to provide decision-making basis for atmospheric environmental governance and future development of urban agglomerations.

Graphical Abstract

1. Introduction

Atmospheric aerosols are solid or liquid particles that are suspended in the air. Atmospheric aerosols have great spatial and temporal variability and their distribution characteristics have a significant influence on the formation of cloud and precipitation [1,2]. Aerosols can reduce the amount of solar radiation and sunlight reaching the ground by absorbing and scattering solar radiation, thereby responding to the Earth’s radiation balance [3]. Additionally, aerosol particles can act as cloud condensation nuclei to affect the optical properties of the cloud, the amount of clouds and the lifetime of the cloud, thereby affecting precipitation and water cycles [4,5]. Furthermore, aerosols can alter atmospheric chemical processes, which in turn affect the concentration and distribution of greenhouse gas [6,7]. Therefore, the monitoring of the spatial and temporal distribution of aerosols has been used as an important basis for accurately assessing the climatic effects of aerosols [8,9] and the impacts of aerosols on air quality [10].
Aerosol optical depth (AOD), which is defined as the total attenuation of absorption and scattering of solar radiation by aerosols along the radiation transmission path, is one of the core parameters widely used for characterizing atmospheric aerosol properties [11] and evaluating their radiation radiative forcing and climatic effects [12]. In addition, AOD has been used as a proxy for particulate pollution by the air quality community [13]. Satellite-derived AOD products provide a unique tool to monitor aerosol properties at regional and global scales. Moderate resolution Imaging Spectroradiometer (MODIS) AOD is widely used to characterize the spatiotemporal change of atmospheric aerosols because it has been extensively validated and provides continuous observations extending for nearly two decades [14,15,16]. Numerous studies have used MODIS AOD to study the long-term spatiotemporal variation of aerosols at regional and global scales [17,18,19,20,21,22,23,24,25].
As the world’s largest developing country, China is experiencing rapid economic expansion, industrial development and urbanization, accompanied by increasing air pollution and aerosol particles are the major harmful pollutants that are responsible for the poor air quality of most urban areas in China. Several studies have quantitatively analyzed the spatial and temporal changes of aerosols in China. Xie and Xia (2008) analyzed the spatial and temporal variations of AOD in north China from 1980 to 2001 by using Total Ozone Mapping Spectrometer (TOMS) monthly AOD data [26]. Guo et al. (2011) conducted a long-term trend analysis of AOD at one-degree over eight typical regions in China, both spatially and temporally, using TOMS AOD products (1980–2001), along with MODIS AOD data (2000–2008) [27]. Luo et al. (2014) analyzed the trends and seasonal variations of AOD over 10 regions in China from 2000 to 2010 by using the Level 3 MODIS AOD products [28]. Che et al. (2015) characterized the aerosol optical properties over a large area of China based on an extensive set of ground-based observations collected from the China Aerosol Remote Sensing Network (CARSNET) during 2002–2013 [29]. He et al. (2016) identified the spatio-temporal trends of AOD in China and quantified the plausible factors that influence AOD, using Aqua MODIS Collection 6 aerosol data obtained between 2002 and 2015 [30]. Zhao et al. (2017) investigated inter-annual AOD trends during 2001–2015 over Eastern and Central China using observations from Multi-angle Imaging SpectroRadiometer (MISR) and MODIS [31]. Sogacheva et al. (2018) analyzed the spatial patterns and inter-annual and seasonal variations of AOD over China based on one-degree Level 3 monthly AOD data retrieved from AlongTrack Scanning Radiometer (ATSR-2, 1995–2002) and Advanced ATSR (AATSR, 2002–2012) and Level 3 Terra MODIS merged Dark Target and Deep Blue (DTDB) AOD product (2000–2017) [24,25].
However, few studies have analyzed the distribution and dynamic characteristics of aerosols over national-level urban agglomerations in China, which are regions with the most developed economy and the highest population density. The five national-level urban agglomerations accounted for approximately half of China’s GDP in 2010 and the population in these regions has been steadily increasing [25,32]. The migration of people from rural areas to these large cities has resulted in fast urbanization and further industrialization and economic growth in these regions [25]. Along with the rapid economic growth and expanding urbanization, these regions are also the most polluted areas in China. Therefore, analysis of the spatial and temporal changes of long-term AOD time series can lay a foundation for the study of atmospheric aerosols in these areas and has practical guiding significance for the sustainable development of urban agglomerations and the construction of an environment-friendly society.
In this study, the latest MODIS Collection 6.1 (C6.1) aerosol products are used to analyze the spatiotemporal variations of AOD over five national-level urban agglomerations in China. The analysis is performed at annual, seasonal and monthly time scales and a time series analysis model is built to explore the dynamic characteristics of the monthly AOD time series. The main goal of this study is to provide a reference for understanding the spatial and temporal change patterns of atmospheric aerosols over these highly urbanized areas.

2. Materials and Methods

2.1. Study Areas

Five national-level urban agglomerations in China are selected for this study: the Yangtze River Delta urban agglomeration (YRD), the Pearl River Delta urban agglomeration (PRD), the Beijing-Tianjin-Hebei urban agglomeration (BTH), the Yangtze River Middle-Reach urban agglomeration (YRMR) and the Cheng-Yu urban agglomeration (CY). The spatial pattern of these urban agglomerations was proposed in the National New Urbanization Plan (2014–2020) issued on March 16, 2013 by the Chinese Central Government [32]. YRD sits along the central-eastern coastline of China. YRMR is located in central China and includes 3 medium-scale urban agglomerations: the Chang-Zhu-Tan urban agglomeration, the Wuhan metropolitan area and the Poyang Lake city group. YRD and YRMR have a subtropical monsoon climate characterized by hot, humid summers. PRD is located in southeastern China and has a humid subtropical climate with frequent cloudy and rainy weather. BTH is in the northern part of the North China Plain, bordering on the Bohai Sea. It has a monsoon-influenced humid continental climate with severe, dry winters. CY lies in the upstream of Yangtze River Basin in southwest China. It has a subtropical monsoon climate characterized by high humidity and many overcast days. The five selected urban agglomerations have stimulated the development of the national economy and have been the current foci of China’s New Urbanization strategy [32,33]. Their boundaries are indicated by different colors in Figure 1.

2.2. MODIS Collection 6.1 Aerosol Products

The MODIS instruments onboard polar orbiting and sun-synchronous Terra and Aqua satellites have 36 spectral channels and provide nearly daily global coverage at spatial resolutions of 250 m, 500 m and 1000 m [34]. The MODIS Adaptive Processing System (MODAPS) generates global AOD products over land at 3 km or 10 km spatial resolutions based on the Dark Target (DT) algorithm and a 10 km spatial resolution based on the Deep Blue (DB) algorithm. It has recently released an improved C6.1 for all MODIS Level 1 (L1) and higher-level Level 2 (L2) & Level 3 (L3) Atmosphere Team products.
In this study, Terra MODIS C6.1 Level 2 aerosol products (MOD04_L2) at 10 km spatial resolution from 1 January 2001 to 31 December 2017 were used. All data were downloaded from the Level 1 and Atmosphere Archive and Distribution System (LAADS) Distributed Active Archive Center (DAAC) at NASA’s Goddard Space Flight Center in Greenbelt, Maryland, United States. The merged DTDB AOD data were extracted and only land pixels were used. We validated the DTDB AOD product over China with Aerosol Robotic Network (AERONET) data for the sites with an observation period longer than six months. A total of 28 sites were selected and AERONET Version 3 Level 2.0 AOD data were used as reference. Satellite data within a 25 km radius around the AERONET site were averaged to matchup with a one-hour AERONET data segment centered on the satellite overpass time [35]. Figure 2 depicts a good agreement between the MODIS and AERONET AOD values. The coefficient of determination (R2) is 0.85 and the root mean squared error (RMSE) is 0.21. The validation results indicate that MODIS C6.1 DTDB AOD is suitable for characterizing long-term AOD variations over China. Statistics of the site-level MODIS/AERONET collocations are listed in Table S1 in the Supplement.

2.3. Multi-Time Scale Analysis of AOD Changes

2.3.1. Calculation of Mean AOD Distributions and Regional Time Series

To analyze the spatial and temporal variation of AOD at annual, seasonal and monthly time scales, the annual, multi-year seasonal and multi-year monthly mean AOD distributions over China are calculated using a direct averaging of all of the valid retrievals. A minimum sampling frequency of 3% for each pixel is applied to exclude poorly sampled pixels.
The regional mean AOD values for the five selected urban agglomerations are aggregated based on their boundaries. Since MODIS AOD values are retrieved only in cloud-free conditions and the derived AOD suffers from non-uniform and incomplete sampling issues, it is expected that a day with more valid pixels will enable a better estimate of the mean than a poorly sampled day. As a result, the area-weighted mean is calculated to reduce the impact of poorly sampled days. That is, each day’s contribution is weighted by the number of valid pixels. The generic regional mean AOD (AODmean) is:
AOD mean = i AOD i × N i i N i
where AODi is the daily regional mean AOD for the day i, Ni is the number of valid pixels on day i.

2.3.2. Time Series Decomposition

In order to find the intrinsic statistical characteristics of AOD variations over the five selected national-level urban agglomerations, time series decomposition is implemented to decompose monthly AOD time series. We build an additive decomposition model to separate monthly AOD time series into three distinct components:
At = Tt + St + It
where At is the AOD value at time t, Tt is the trend component, St is the seasonal component and It represents the irregular component. The trend component is the trend of long-term changes in the AOD time series that gradually increased or decreased with time. It is a trend curve that reflects the general direction of change of the monthly AOD time series. The seasonal component is the fixed variation exhibited by the AOD time series within one year and the main reason for the seasonal changes is the impact of the natural seasonal alternation. Irregular component reflects random changes and it is unpredictable. Irregular changes may occur due to natural disasters, human activities and sudden changes in climate.
Tt is estimated by using a quadratic trend model:
Tt = α0 + α1t + α2t2
where α0, α1 and α2 are regression coefficients. Then Tt is subtracted from the original data At to obtain the de-trended series Dt. St is modeled as a function of seasonal indicator variables:
St = γ1δ1t + γ2δ2t + … + γ12δ12t
where γi is the coefficient associated with the seasonal indicator variable and δit is an indicator variable taking the value 1 in the ith month of the year and 0 in all other months. The de-trended series Dt is regressed against the seasonal indicators to estimate the coefficients of the indicator variables and then the seasonal component St can be calculated. The irregular component It is estimated by subtracting Tt and St from At.

3. Results

3.1. AOD Changes at Different Time Scales

3.1.1. Annual AOD Changes in Urban Agglomerations

The spatial distributions of annual averaged AOD in China during 2001–2017 are presented in Figure 3. It reveals that the spatial pattern of annual averaged AOD has obvious regional characteristics and is strongly influenced by topography. Since the topography of China is high in the west and low in the east (Figure 1), it disrupts the transmission of high concentrations of aerosol particles from the east to the west. In addition, the eastern part of China has many large industrial cities which are gathering centers for energy consumption and anthropogenic emission. These areas are highly prone to form areas with high AOD values. In YRMR, large industrial cities such as Wuhan, Nanchang and Changsha are areas with high AOD values. This is because the diffusion of pollutants is hindered by the surrounding mountains and the hygroscopic growth of aerosols is enhanced by excessive humidity caused by intensive lakes and rivers. As an important platform for the great development of western China, CY is the strategic support of the Yangtze River Economic Belt. However, the relatively occluded geographical environment of the Sichuan Basin makes it difficult for the aerosol particles generated in CY to be transported outside the basin and aerosol particles aggregate to some extent due to the foggy weather. As a result, the observed AOD is high in this region.
Figure 4 indicates that the temporal variations of annual averaged AOD values in five urban agglomerations are quite different. The annual averaged AOD values of YRD vary from 0.49 to 0.76 during 2001–2017. The change of AOD values in YRD shows a growth trend from 2001 to 2007, while a downward trend is observed during 2011–2017. A similar pattern is observed in YRMR, except the initial upward trend ends at 2006 and AOD values in this region are lower than that of YRD between 2007 and 2017. The annual averaged AOD values of CY vary from 0.39 to 0.83. The peak AOD value is observed in 2006 and a continuous downward trend is observed from 2011 onward. The annual averaged AOD of PRD is observed of the order of 0.5–0.7 during 2001–2017. The AOD in PRD gradually decreases after 2012 and the minimum value of 0.35 is observed in 2016. The annual averaged AOD values of BTH are lower than other regions. It is because the western and northern parts of this region are the Xishan and Yanshan mountain ranges, which are low AOD areas. Therefore, regional mean AOD values are biased toward low values despite the AOD is high over the southern part of this region.
In addition, Figure 3 and Figure 4 show that AOD typically increases before 2006 and fluctuates during 2006–2011 and then decreases after 2011 over the selected urban agglomerations. The upward trend observed before 2006 is mainly due to emission increases caused by rapid economic and industrial development. In the 11th five-year plan (2006–2010), China set targets to improve energy efficiency and reduce emissions of sulfur dioxide and primary aerosols. During the 12th five-year plan (2011–2015), the released environmental protection policies extended the previous policies and aimed to improve the air quality further [36]. Due to the emission reductions, the trend reversal observed in 2011 might be attributed to the implementation of these environmental protection policies [37,38]. In September 2013, China issued the Air Pollution Prevention and Control Action Plan, which aimed to reduce particulate matter concentrations across the country [36]. The Action Plan strengthened the targets of air quality improvement and raised stringent air pollution control policies. It seems that the Action Plan has led to a significant drop of AOD since 2014.

3.1.2. Seasonal AOD Changes in Urban Agglomerations

The spatial distributions of seasonal averaged AOD in China during 2001–2017 are presented in Figure 5. It reveals that the averaged AOD values in most parts of the five selected urban agglomerations in spring and summer are significantly higher than autumn and winter. The spatial distribution of areas with high and low AOD values of different urban agglomerations is basically unchanged in the four seasons but the intensity of these areas reflects their respective regional seasonal changes.
Figure 6 compares the seasonal averaged AOD values in five urban agglomerations during 2001–2017. For the five selected regions, AOD is highest in spring, except for BTH where the maximum AOD is observed in summer. In PRD and YRD, AOD values are similar in spring and summer. The high summer AOD over the two regions can be partially caused by the hygroscopic growth of aerosols under high relative humidity conditions. For PRD, YRD and YRMR, AOD is lowest in winter, whereas the minimum AOD is observed in autumn for BTH and CY. In BTH, intensive wintertime warming activities enhance the anthropogenic emission, thus the AOD increases in winter. For CY, the high winter AOD is mainly caused by temperature inversions [39]. Since the Sichuan Basin is surrounded by mountains, the strong inversion layer above the atmospheric boundary layer will trap any pollutant that enters the air and therefore haze gradually builds up.

3.1.3. Monthly AOD Changes in Urban Agglomerations

Figure 7 presents the spatial distributions of monthly averaged AOD in China during 2001–2017. The spatial distribution of areas with high and low AOD values shows similar patterns as the seasonal averaged AOD. However, monthly maps can provide useful information for qualifying the temporal variations of AOD due to different weather conditions. AOD values show an increasing trend during spring and early summer months due to the continuous increase of temperature and humidity. The monthly averaged AOD is prone to yielding large amounts of data gaps in June and July over PRD due to its frequent cloudy and rainy days. In addition, significant amounts of missing AOD areas are observed in December and January over CY because of frequent fog events in the Sichuan basin.
Figure 8 compares the monthly averaged AOD values in five urban agglomerations. Typically, AOD varies seasonally from low in cold months to high in warm months and then decreases during the rainy periods. In BTH, YRD and YRMD, the highest AOD is observed in June. This is because the progressive increase of air temperature from January to June will accelerate photochemical reactions and increase aerosol concentrations [40]. In addition, the dust aerosol transport from the arid and semi-arid regions also contributes to the increase in aerosol loading [24]. Although high relative humidity during summer months strengthens the local aerosol accumulation, East Asian summer monsoon rainfalls wash out aerosol particles and result in low AOD values over these regions. Variations of monthly averaged AOD are much more irregular over PRD and CY than over the other regions, the zigzag pattern may be caused by the sampling issues of satellite measurements due to the frequent foggy, cloudy or rainy days over the two regions.

3.2. AOD Changes Using Time Series Decomposition

Time series decompositions of monthly AOD time series in the five selected urban agglomeration are shown in Figure 9, Figure 10, Figure 11, Figure 12 and Figure 13. In these figures, SC represents the seasonal component and IC represents the irregular component.
In YRD, the trend component of AOD time series fluctuates between 0.51 and 0.70. AOD is in the rising phase during 2001–2008 and then it decreases. The maximum value of the seasonal component is 0.35 occurred in June and the minimum seasonal component is −0.19 occurred in December.
In PRD, the trend component of AOD time series fluctuates between 0.41 and 0.67. AOD is in the rising stage between 2001 and 2007 and then a decline phase is observed. The maximum value of the seasonal component is 0.19 occurred in April and the minimum seasonal component is −0.16 occurred in November.
In BTH, the trend component of AOD time series fluctuates between 0.39 and 0.52. AOD is in the rising stage during 2001–2008 and then it decreases. The maximum value of the seasonal component is 0.24 occurred in June and the minimum seasonal component is −0.21 occurred in November.
In YRMR, the trend component of AOD time series fluctuates between 0.41 and 0.68. AOD is in a rising phase between 2001 and 2007 and then decreases in the period 2008–2017. The maximum value of the seasonal component is 0.21 occurred in June and the minimum seasonal component is −0.17 occurred in December.
In CY, the trend component of AOD time series fluctuates between 0.37 and 0.69. AOD is in a rising phase during 2001–2007 and then a decline phase is observed. The maximum value of the seasonal component is 0.18 occurred in March and the minimum seasonal component is −0.20 occurred in December.
The trend component of monthly averaged AOD time series in the five selected urban agglomerations indicates that AOD increases at first and then decreases between 2001 and 2017 and the transition occurs during 2007–2008. In November 2007, the State Council of China approved the 11th five-year plan for environmental protection, which aimed to improve energy efficiency and reduce emissions of sulfur dioxide and primary aerosols [41]. As a result, the observed trend reversal may be attributed to the implementation of environmental measures. The seasonal component typically has positive values in spring and summer and negative values in autumn and winter. This is in accordance with the fact that high AOD is observed in spring and summer, as is shown in Figure 5. Extreme values of the irregular component are typically observed during rainy months and this may be caused by the sampling issues of satellite measurements. That is, the regional mean AOD is biased because a large number of pixels are missed due to frequent cloudy and rainy days.

4. Discussion

Sogacheva et al. (2018) investigated the long-term AOD variations over China during 1995–2017 by combining AOD data retrieved from ATSR-2, AATSR and MODIS instruments. Linear regression was applied to estimate the AOD tendencies for two periods: 1995–2006 and 2011–2017 [25]. Although the choice of study regions is different from this paper, it also covers major urban/industrial regions. Their results also revealed an increasing tendency before 2006 and a reverse tendency after 2011 in BTH, YRD, PRD and CY.
To calculate the regional mean AOD at different time scales, we weighted the daily regional mean AOD by the number of valid pixels, bearing in mind that a day with more valid pixels would enable a better estimate of the mean than a poorly sampled day. However, the calculated mean AOD might still be biased when large amounts of data are missing during the selected period. Figure 14 indicates that the sampling frequencies of AOD show significant seasonal changes. Typically, the sampling frequency is lower than 10% in summer over PRD and in winter over CY. This is because AOD data are only available during clear-sky conditions. Clouds, fog and precipitation will lead to a large amount of missing AOD data. Figure 15 presents the seasonal distributions of the standard error of the mean (SEM) derived from AOD data. SEM is used as an indirect measure of the reliability of the estimated mean. It depends on both the variability of AOD and the number of observations. For example, the cloudy fog in the Sichuan Basin in winter makes it difficult to obtain adequate amounts of effective observations and large SEM values indicate that mean AOD values have a low reliability over this region. That is, the derived mean AOD might be unrepresentative and the sampling issue of satellite-derived AOD products should be taken into account.
A time series decomposition model was used to analyze the intrinsic statistical characteristics of AOD dynamics in this study, whereas there are two issues still deserving to be further discussed. One issue is about the premise of time series decomposition. The time series decomposition model assumes that the seasonal factors are the same for each year, whereas climate change can lead to different seasonal factors in different year, such as changes in clouds and water vapor, wind strength and the amount of sand and dust. This study does not address these seasonal changes and the impact of seasonal factors on aerosols. Another issue is associated with outliers. Although time series decomposition can accurately locate outliers, it does not handle outliers well. Further efforts should be made to eliminate the effects of outliers in the time series analysis model.

5. Conclusions

Using the latest MODIS Collection 6.1 (C6.1) Level 2 aerosol products, this paper analyzed the change patterns of AOD over five national-level urban agglomerations in China at multi-time scales and built a time series decomposition model to study the long-term dynamics of AOD in the five national-level urban agglomerations in China from 2001 to 2017. The spatial distribution of AOD in China had obvious regional characteristics and high AOD regions were mainly concentrated in east China where most urban agglomerations located. During the study period, changes of annual averaged AOD values in five urban agglomerations were basically increased first and then decreased. Time series decompositions also indicated that long-term AOD trends in the five selected urban agglomerations were basically rising first and then falling and the transition occurred during 2007–2008. The AOD of the five selected urban agglomerations was higher in spring and summer than that in autumn and winter. Although the seasonal variations of AOD in these urban agglomerations were different, the AOD typically varied from low in cold months to high in warm months and then decreased during the rainy periods. The detailed analysis of trends in AODs over different urban agglomerations could provide useful references for the assessments of regional climate and environmental impacts due to aerosols.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-4292/11/2/201/s1, Table S1: Statistics of the site-level comparisons between MODIS Collection 6.1 DTDB AOD products and AERONET observations in China.

Author Contributions

Conceptualization, M.Z. and L.Z.; Methodology, M.Z.; Validation, M.Z. and L.Z.; Formal analysis, M.Z., L.Z. and Y.Y.; Investigation, Y.Y.; Resources, L.Z.; Supervision, Y.Y.; Writing—original draft preparation, M.Z. and L.Z.; Writing—review and editing, L.Z. and Y.Y.

Funding

This research was funded by National Postdoctoral Program for Innovative Talents to Lei Zhang (grant number BX201700175), China Postdoctoral Science Foundation (grant number 2018M630882) and Natural Science Foundation of Hubei Province of China (grant number 2018CFB341).

Acknowledgments

The authors would like to acknowledge the NASA Goddard Earth Science Distributed Active Archive Centre for the MODIS Level 2 data and we would like to thank the principal investigators of the AERONET stations used in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rosenfeld, D.; Lohmann, U.; Raga, G.B.; O’dowd, C.D.; Kulmala, M.; Fuzzi, S.; Reissell, A.; Andreae, M.O. Flood or drought: How do aerosols affect precipitation? Science 2008, 321, 1309–1313. [Google Scholar] [CrossRef] [PubMed]
  2. Li, Z.; Niu, F.; Fan, J.; Liu, Y.; Rosenfeld, D.; Ding, Y. Long-term impacts of aerosols on the vertical development of clouds and precipitation. Nat. Geosci. 2011, 4, 888–894. [Google Scholar] [CrossRef]
  3. Ramanathan, V.; Crutzen, P.; Kiehl, J.; Rosenfeld, D. Aerosols, climate and the hydrological cycle. Science 2001, 294, 2119–2124. [Google Scholar] [CrossRef] [PubMed]
  4. Liepert, B.G. Observed reductions of surface solar radiation at sites in the United States and worldwide from 1961 to 1990. Geophys. Res. Lett. 2002, 29, 1–4. [Google Scholar] [CrossRef]
  5. Mishchenko, M.I.; Cairns, B.; Kopp, G.; Schueler, C.F.; Fafaul, B.A.; Hansen, J.E.; Hooker, R.J.; Itchkawich, T.; Maring, H.B.; Travis, L.D. Accurate monitoring of terrestrial aerosols and total solar irradiance: Introducing the Glory Mission. Bulle. Am. Meteorol. Soc. 2007, 88, 677–692. [Google Scholar] [CrossRef]
  6. Smith, K.R.; Jerrett, M.; Anderson, H.R.; Burnett, R.T.; Stone, V.; Derwent, R.; Atkinson, R.W.; Cohen, A.; Shonkoff, S.B.; Krewski, D. Public health benefits of strategies to reduce greenhouse-gas emissions: Health implications of short-lived greenhouse pollutants. Lancet 2009, 374, 2091–2103. [Google Scholar] [CrossRef]
  7. Delworth, T.L.; Dixon, K.W. Have anthropogenic aerosols delayed a greenhouse gas-induced weakening of the North Atlantic thermohaline circulation? Geophys. Res. Lett. 2006, 33, 1–4. [Google Scholar] [CrossRef]
  8. Mishchenko, M.I.; Cairns, B.; Hansen, J.E.; Travis, L.D.; Burg, R.; Kaufman, Y.J.; Martins, J.V.; Shettle, E.P. Monitoring of aerosol forcing of climate from space: Analysis of measurement requirements. J. Quant. Spectrosc. Radiat. Transf. 2004, 88, 149–161. [Google Scholar] [CrossRef]
  9. Liu, X.-Y.; Zhang, Y.; Zhang, Q.; He, K.-B. Application of online-coupled WRF/Chem-MADRID in East Asia: Model evaluation and climatic effects of anthropogenic aerosols. Atmos. Environ. 2016, 124, 321–336. [Google Scholar] [CrossRef] [Green Version]
  10. Chin, M.; Diehl, T.; Ginoux, P.; Malm, W. Intercontinental transport of pollution and dust aerosols: Implications for regional air quality. Atmos. Chem. Phys. 2007, 7, 5501–5517. [Google Scholar] [CrossRef]
  11. Li, J.; Jiang, Y.; Xia, X.; Hu, Y. Increase of surface solar irradiance across East China related to changes in aerosol properties during the past decade. Environ. Res. Lett. 2018, 13, 1–9. [Google Scholar] [CrossRef]
  12. Chylek, P.; Henderson, B.; Mishchenko, M. Aerosol radiative forcing and the accuracy of satellite aerosol optical depth retrieval. J. Geophys. Res. Atmos. 2003, 108, 1–8. [Google Scholar] [CrossRef]
  13. An, X.-Q.; Zhu, T.; Wang, Z.-F.; Li, C.-Y.; Wang, Y.-S. A modeling analysis of a heavy air pollution episode occurred in Beijing. Atmos. Chem. Phys. 2007, 7, 3103–3114. [Google Scholar] [CrossRef] [Green Version]
  14. Sayer, A.M.; Munchak, L.A.; Hsu, N.C.; Levy, R.C.; Bettenhausen, C.; Jeong, M.J. MODIS Collection 6 aerosol products: Comparison between Aqua’s e-Deep Blue, Dark Target and “merged” data sets and usage recommendations. J. Geophys. Res. Atmos. 2014, 119, 13965–13989. [Google Scholar] [CrossRef]
  15. Tian, X.; Liu, Q.; Li, X.; Wei, J. Validation and Comparison of MODIS C6.1 and C6 Aerosol Products over Beijing, China. Remote Sens. 2018, 10, 2021. [Google Scholar] [CrossRef]
  16. Zhang, M.; Huang, B.; He, Q. An Evaluation of Four MODIS Collection 6 Aerosol Products in a Humid Subtropical Region. Remote Sens. 2017, 9, 1173. [Google Scholar] [CrossRef]
  17. Ramachandran, S.; Kedia, S.; Srivastava, R. Aerosol optical depth trends over different regions of India. Atmos. Environ. 2012, 49, 338–347. [Google Scholar] [CrossRef]
  18. Kumar, K.R.; Sivakumar, V.; Yin, Y.; Reddy, R.; Kang, N.; Diao, Y.; Adesina, A.J.; Yu, X. Long-term (2003–2013) climatological trends and variations in aerosol optical parameters retrieved from MODIS over three stations in South Africa. Atmos. Environ. 2014, 95, 400–408. [Google Scholar] [CrossRef]
  19. Mehta, M.; Singh, R.; Singh, A.; Singh, N. Recent global aerosol optical depth variations and trends—A comparative study using MODIS and MISR level 3 datasets. Remote Sens. Environ. 2016, 181, 137–150. [Google Scholar] [CrossRef]
  20. Floutsi, A.A.; Korras-Carraca, M.B.; Matsoukas, C.; Hatzianastassiou, N.; Biskos, G. Climatology and trends of aerosol optical depth over the Mediterranean basin during the last 12 years (2002–2014) based on Collection 006 MODIS-Aqua data. Sci. Total Environ. 2016, 551–552, 292–303. [Google Scholar] [CrossRef]
  21. Zhang, J.L.; Reid, J.S.; Alfaro-Contreras, R.; Xian, P. Has China been exporting less particulate air pollution over the past decade? Geophys. Res. Lett. 2017, 44, 2941–2948. [Google Scholar] [CrossRef]
  22. Shokr, M.; El-Tahan, M.; Ibrahim, A.; Steiner, A.; Gad, N. Long-Term, High-Resolution Survey of Atmospheric Aerosols over Egypt with NASA’s MODIS Data. Remote Sens. 2017, 9, 1027. [Google Scholar] [CrossRef]
  23. De Leeuw, G.; Sogacheva, L.; Rodriguez, E.; Kourtidis, K.; Georgoulias, A.K.; Alexandri, G.; Amiridis, V.; Proestakis, E.; Marinou, E.; Xue, Y.; et al. Two decades of satellite observations of AOD over mainland China using ATSR-2, AATSR and MODIS/Terra: Data set evaluation and large-scale patterns. Atmos. Chem. Phys. 2018, 18, 1573–1592. [Google Scholar] [CrossRef]
  24. Sogacheva, L.; de Leeuw, G.; Rodriguez, E.; Kolmonen, P.; Georgoulias, A.K.; Alexandri, G.; Kourtidis, K.; Proestakis, E.; Marinou, E.; Amiridis, V.; et al. Spatial and seasonal variations of aerosols over China from two decades of multi-satellite observations—Part 1: ATSR (1995–2011) and MODIS C6.1 (2000–2017). Atmos. Chem. Phys. 2018, 18, 11389–11407. [Google Scholar] [CrossRef]
  25. Sogacheva, L.; Rodriguez, E.; Kolmonen, P.; Virtanen, T.H.; Saponaro, G.; de Leeuw, G.; Georgoulias, A.K.; Alexandri, G.; Kourtidis, K.; van der A, R.J. Spatial and seasonal variations of aerosols over China from two decades of multi-satellite observations—Part 2: AOD time series for 1995–2017 combined from ATSR ADV and MODIS C6.1 and AOD tendency estimations. Atmos. Chem. Phys. 2018, 18, 16631–16652. [Google Scholar] [CrossRef]
  26. Xie, J.; Xia, X. Long-term trend in aerosol optical depth from 1980 to 2001 in north China. Particuology 2008, 6, 106–111. [Google Scholar] [CrossRef]
  27. Guo, J.-P.; Zhang, X.-Y.; Wu, Y.-R.; Zhaxi, Y.; Che, H.-Z.; La, B.; Wang, W.; Li, X.-W. Spatio-temporal variation trends of satellite-based aerosol optical depth in China during 1980–2008. Atmos. Environ. 2011, 45, 6802–6811. [Google Scholar] [CrossRef]
  28. Luo, Y.; Zheng, X.; Zhao, T.; Chen, J. A climatology of aerosol optical depth over China from recent 10 years of MODIS remote sensing data. Int. J. Climatol. 2014, 34, 863–870. [Google Scholar] [CrossRef]
  29. Che, H.; Zhang, X.-Y.; Xia, X.; Goloub, P.; Holben, B.; Zhao, H.; Wang, Y.; Zhang, X.-C.; Wang, H.; Blarel, L. Ground-based aerosol climatology of China: Aerosol optical depths from the China Aerosol Remote Sensing Network (CARSNET) 2002–2013. Atmos. Chem. Phys. 2015, 15, 7619–7652. [Google Scholar] [CrossRef]
  30. He, Q.; Zhang, M.; Huang, B. Spatio-temporal variation and impact factors analysis of satellite-based aerosol optical depth over China from 2002 to 2015. Atmos. Environ. 2016, 129, 79–90. [Google Scholar] [CrossRef]
  31. Zhao, B.; Jiang, J.H.; Gu, Y.; Diner, D.; Worden, J.; Liou, K.-N.; Su, H.; Xing, J.; Garay, M.; Huang, L. Decadal-scale trends in regional aerosol particle properties and their linkage to emission changes. Environ. Res. Lett. 2017, 12, 054021. [Google Scholar] [CrossRef] [Green Version]
  32. Fang, C. Important progress and future direction of studies on China’s urban agglomerations. J. Geogr. Sci. 2015, 25, 1003–1024. [Google Scholar] [CrossRef] [Green Version]
  33. Fang, C.; Yu, D. Urban agglomeration: An evolving concept of an emerging phenomenon. Landsc. Urban Plan. 2017, 162, 126–136. [Google Scholar] [CrossRef]
  34. Turner, W.; Spector, S.; Gardiner, N.; Fladeland, M.; Sterling, E.; Steininger, M. Remote sensing for biodiversity science and conservation. Trends Ecol. Evol. 2003, 18, 306–314. [Google Scholar] [CrossRef]
  35. Petrenko, M.; Ichoku, C. Coherent uncertainty analysis of aerosol measurements from multiple satellite sensors. Atmos. Chem. Phys. 2013, 13, 6777–6805. [Google Scholar] [CrossRef] [Green Version]
  36. Cai, S.; Wang, Y.; Zhao, B.; Wang, S.; Chang, X.; Hao, J. The impact of the “Air Pollution Prevention and Control Action Plan” on PM 2.5 concentrations in Jing-Jin-Ji region during 2012–2020. Sci. Total Environ. 2017, 580, 197–209. [Google Scholar] [CrossRef]
  37. Georgoulias, A.K.; van der A, R.J.; Stammes, P.; Boersma, K.F.; Eskes, H.J. Trends and trend reversal detection in two decades of tropospheric NO2 satellite observations. Atmos. Chem. Phys. Discuss. 2018, 2018, 1–38. [Google Scholar]
  38. Ma, Z.; Liu, R.; Liu, Y.; Bi, J. Effects of air pollution control policies on PM2.5 pollution improvement in China from 2005 to 2017: A satellite based perspective. Atmos. Chem. Phys. Discuss. 2018, 2018, 1–23. [Google Scholar]
  39. Ning, G.; Wang, S.; Yim, S.H.L.; Li, J.; Hu, Y.; Shang, Z.; Wang, J.; Wang, J. Impact of low-pressure systems on winter heavy air pollution in the northwest Sichuan Basin, China. Atmos. Chem. Phys. 2018, 18, 13601–13615. [Google Scholar] [CrossRef]
  40. Wang, D.; Zhou, B.; Fu, Q.; Zhao, Q.; Zhang, Q.; Chen, J.; Yang, X.; Duan, Y.; Li, J. Intense secondary aerosol formation due to strong atmospheric photochemical reactions in summer: Observations at a rural site in eastern Yangtze River Delta of China. Sci. Total Environ. 2016, 571, 1454–1466. [Google Scholar] [CrossRef]
  41. Zhao, B.; Wang, S.; Wang, J.; Fu, J.S.; Liu, T.; Xu, J.; Fu, X.; Hao, J. Impact of national NOx and SO2 control policies on particulate matter pollution in China. Atmos. Environ. 2013, 77, 453–463. [Google Scholar] [CrossRef]
Figure 1. Locations of five selected urban agglomerations in China.
Figure 1. Locations of five selected urban agglomerations in China.
Remotesensing 11 00201 g001
Figure 2. Density scatterplot of MODIS C6.1 merged Dark Target and Deep Blue (DTDB) AOD vs. AOD from AERONET stations in China during 2001–2017. A MODIS-AERONET collocation is the spatial average of MODIS AOD within a 25 km radius around the AERONET site and the temporal average of all AERONET observations within ±30 min of satellite overpass.
Figure 2. Density scatterplot of MODIS C6.1 merged Dark Target and Deep Blue (DTDB) AOD vs. AOD from AERONET stations in China during 2001–2017. A MODIS-AERONET collocation is the spatial average of MODIS AOD within a 25 km radius around the AERONET site and the temporal average of all AERONET observations within ±30 min of satellite overpass.
Remotesensing 11 00201 g002
Figure 3. The annual AOD changes over five urban agglomerations in China.
Figure 3. The annual AOD changes over five urban agglomerations in China.
Remotesensing 11 00201 g003
Figure 4. Annual averaged AOD values in five urban agglomerations.
Figure 4. Annual averaged AOD values in five urban agglomerations.
Remotesensing 11 00201 g004
Figure 5. The seasonal AOD changes over five urban agglomerations in China.
Figure 5. The seasonal AOD changes over five urban agglomerations in China.
Remotesensing 11 00201 g005
Figure 6. Seasonal averaged AOD values in five urban agglomerations.
Figure 6. Seasonal averaged AOD values in five urban agglomerations.
Remotesensing 11 00201 g006
Figure 7. The monthly AOD changes over five urban agglomerations in China.
Figure 7. The monthly AOD changes over five urban agglomerations in China.
Remotesensing 11 00201 g007
Figure 8. Monthly averaged AOD values in five urban agglomerations.
Figure 8. Monthly averaged AOD values in five urban agglomerations.
Remotesensing 11 00201 g008
Figure 9. Time series decomposition of monthly AOD time series in the Yangtze River Delta urban agglomeration.
Figure 9. Time series decomposition of monthly AOD time series in the Yangtze River Delta urban agglomeration.
Remotesensing 11 00201 g009
Figure 10. Time series decomposition of monthly AOD time series in the Pearl River Delta urban agglomeration.
Figure 10. Time series decomposition of monthly AOD time series in the Pearl River Delta urban agglomeration.
Remotesensing 11 00201 g010
Figure 11. Time series decomposition of monthly AOD time series in the Beijing-Tianjin-Hebei urban agglomeration.
Figure 11. Time series decomposition of monthly AOD time series in the Beijing-Tianjin-Hebei urban agglomeration.
Remotesensing 11 00201 g011
Figure 12. Time series decomposition of monthly AOD time series in the Yangtze River Middle-Reach urban agglomeration.
Figure 12. Time series decomposition of monthly AOD time series in the Yangtze River Middle-Reach urban agglomeration.
Remotesensing 11 00201 g012
Figure 13. Time series decomposition of monthly AOD time series in the Cheng-Yu urban agglomeration.
Figure 13. Time series decomposition of monthly AOD time series in the Cheng-Yu urban agglomeration.
Remotesensing 11 00201 g013
Figure 14. Seasonal distributions of AOD sampling frequency over five urban agglomerations in China.
Figure 14. Seasonal distributions of AOD sampling frequency over five urban agglomerations in China.
Remotesensing 11 00201 g014
Figure 15. Seasonal distributions of the standard error of the mean (SEM) over five urban agglomerations in China.
Figure 15. Seasonal distributions of the standard error of the mean (SEM) over five urban agglomerations in China.
Remotesensing 11 00201 g015

Share and Cite

MDPI and ACS Style

Zhang, L.; Zhang, M.; Yao, Y. Multi-Time Scale Analysis of Regional Aerosol Optical Depth Changes in National-Level Urban Agglomerations in China Using Modis Collection 6.1 Datasets from 2001 to 2017. Remote Sens. 2019, 11, 201. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11020201

AMA Style

Zhang L, Zhang M, Yao Y. Multi-Time Scale Analysis of Regional Aerosol Optical Depth Changes in National-Level Urban Agglomerations in China Using Modis Collection 6.1 Datasets from 2001 to 2017. Remote Sensing. 2019; 11(2):201. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11020201

Chicago/Turabian Style

Zhang, Lei, Ming Zhang, and Yibin Yao. 2019. "Multi-Time Scale Analysis of Regional Aerosol Optical Depth Changes in National-Level Urban Agglomerations in China Using Modis Collection 6.1 Datasets from 2001 to 2017" Remote Sensing 11, no. 2: 201. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11020201

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop