2.1 Spatial and temporal domain

We focus on pollution over Eastern China (25–41^∘ N, 110–123^∘ E). Guided by an EOF analysis, we contrast pollution over the southern (SEC, south of 35^∘ N) and northern (NEC) parts to address the regional differences in day-to-day pollution variability. Such latitudinal separation coincides with the Huai River climate transitional zone (Ye and Li, 2017). The orange lines in Fig. 1 separate the two regions.

Our study period is from 25 October to 25 December 2013, with a total of 1488 h in 62 days. Most air pollution data are missing in January and February 2014 because of instrumental failure or data retrieval failure, and data before 25 October are not available.

2.2 NO₂ and PM_2.5 observations

We retrieve hourly measurements of NO₂ and PM_2.5 from 193 air quality monitoring stations of the MEP. Most stations are located in urban areas, and only six stations are suburban. As almost every station has missing values in more than one day, we exclude stations that have missing values at ≥30 % of the 1488 h or during a consecutive 72 h period. We thus select 163 stations for NO₂ and 159 stations for PM_2.5 in the same 42 cities. The dots in Fig. 1a and b depict the stations and cities, respectively. The blue dots show stations and cities with both valid NO₂ and PM_2.5, the green dots with NO₂ only, and the purple dots with PM_2.5 only. The slight difference between NO₂ and PM_2.5 stations does not affect our analysis of the regional pattern of pollutants.

2.2.1 Correction of raw NO₂ measurements

At the monitoring sites, NO₂ is measured via molybdenum-catalyzed conversion to nitric oxide (NO) and a subsequent chemiluminescence measurement. The measurement technique suffers from interference by more oxidized nitrogen species, since the heated molybdenum surface exhibits low chemical selectivity (Boersma et al., 2009; Lamsal et al., 2008; L. Zhang et al., 2016)

Here we follow Lamsal et al. (2008) to correct for the interference by introducing a correction factor (CF) based on GEOS-Chem-simulated nitrogen species (NO₂, HNO₃, PAN, and all alkyl nitrates (∑AN)).

$$\begin{array}{}\text{(1)}& \mathrm{CF}={\displaystyle \frac{{\mathrm{NO}}_{\mathrm{2}}}{{\mathrm{NO}}_{\mathrm{2}}+\sum \mathrm{AN}+\mathrm{0.95}\mathrm{PAN}+\mathrm{0.35}{\mathrm{HNO}}_{\mathrm{3}}}}\end{array}$$

We multiply CF with the raw NO₂ data to obtain “corrected” NO₂ concentrations. Our sensitivity test suggests that assuming PAN and HNO₃ to be fully converted to NO₂ (i.e., assuming the coefficients to be unity for both PAN and HNO₃ in Eq. 1) does not affect our spatiotemporal analysis of NO₂. Hereafter the NO₂ corrected by Eq. (1) is discussed, unless stated otherwise.

Figure 2Regional mean hourly time series of raw and “corrected” NO₂ from the observations. The gray shading indicates 1 standard deviation across all stations.

Download

Figure 2 compares the regional mean hourly time series of raw and corrected NO₂. The correction reduces NO₂ concentrations by about 2–30 µg m⁻³ over the whole period and is higher at times when nitrogen is more oxidized. It slightly reduces the relative contribution of day-to-day variability to the total variance of NO₂ under the EOF–EEMD analysis (not shown) because excluding the more oxidized species shortens the lifetime of NO₂.

2.3 Meteorological observations

We use 3-hourly measurements of 2 m air temperature, 2 m relative humidity, and 10 m wind speed from meteorological stations recorded at the National Oceanic and Atmospheric Administration National Centers for Environment Information (NOAA NCEI). We do not use surface pressure additionally because it is highly correlated with air temperature and relative humidity on the day-to-day scale. The locations of these stations do not always coincide with air pollution stations. Thus, we select 36 meteorological stations within 10 km of air pollution stations (red hollow dots in Fig. 1). Despite the difference (in number and location) between pollution and meteorological stations, an analysis of the regional temporal patterns of pollutants and meteorology is still informative (see Sect. 3.2).

To fill in missing values, we apply an interpolation process that accounts for diurnal variability using information for an adjacent day. For example, if the temperature on 26 October at 12:00 is missing, we calculate the temperature difference between 09:00 and 12:00 on the 25th as well as the difference between 15:00 and 12:00 on the 25th. We then use these differences to adjust the temperatures at 09:00 and 15:00 on the 26th, and finally use the mean of the two adjusted temperatures as the temperature on the 26th at 12:00.

For consistency with the hourly pollution data, we linearly interpolate the 3-hourly meteorological measurements to each hour. This interpolation does not distort the EOF–EEMD analysis, as confirmed by comparing the statistical analysis on 1-hourly GEOS-FP meteorological parameters versus an analysis on 3-hourly GEOS-FP data. Note that the GEOS-FP meteorology is used to drive GEOS-Chem.

2.4 Model simulations

2.5 GEOS-Chem

We use the nested GEOS-Chem CTM version 9-02 (L. Zhang et al., 2016) to simulate NO₂, PM_2.5, and other pollutants over China in October–December 2013. The model resolution is a 0.3125^∘ long. × 0.25^∘ lat. grid with 47 vertical layers, and the lowest 10 layers are of ∼130 m thickness each. The model is driven by the GEOS-FP assimilated meteorology from the National Aeronautics and Space Administration (NASA) Global Modeling and Assimilation Office, with the full O_x–NO_x–VOC–CO–HO_x gaseous chemistry (Mao et al., 2013) and online aerosol calculations. Vertical mixing in the PBL adopts a nonlocal scheme (Holtslag and Boville, 1993; Lin et al., 2010). Model convection is simulated with the relaxed Arakawa–Schubert scheme (Rienecker et al., 2008).

Chinese anthropogenic emissions of NO_x and other pollutants adopt the monthly MEIC inventory with a base year of 2010 (http://www.meicmodel.org, last access: 1 December 2015) (Geng et al., 2017). We further use the monthly DOMINO v2 NO₂ data to scale monthly anthropogenic NO_x emissions from 2010 to the simulation year (Lin et al., 2015). The emission scaling improves the simulation of NO₂ (Cui et al., 2016). Other model setups are referred to Lin et al. (2015) and Yan et al. (2016).

GEOS-Chem modeled PM_2.5 includes secondary inorganic aerosols (sulfate, nitrate, and ammonium), black carbon, primary organic carbon, natural dust, and sea salt. Secondary organic aerosols are not included in this study, considering the severe underestimate in China due to missing precursor emissions and formation pathways (Fu et al., 2012; L. Zhang et al., 2016). Anthropogenic dust is also not included.

The nested model simulation is from 15 October to 25 December in 2013, allowing for a 10-day spin-up period. Its lateral boundary conditions of chemicals are updated every 3 h by results from a corresponding global simulation on a 2.5^∘ long. × 2^∘ lat. grid. Modeled NO₂ and PM_2.5 in the first layer are sampled at city centers and times with valid observations, unless stated otherwise.

2.6 EOF–EEMD analysis visualization package

Figure 3The flowchart of the EOF–EEMD analysis visualization package. The red boxes represent the quantities visualized.

Download

As shown in Fig. 3, our EOF–EEMD analysis visualization package consists, in order, of an EOF analysis (Lorenz, 1956), an EEMD analysis (Wu et al., 2009), a Hilbert transform (HT) with marginal spectrum analysis (MSA), and a visualization step to quantitatively depict the spatial–temporal scales of measurement or model data.

The basic purpose of our package is to quickly and simultaneously identify and visualize various spatial and temporal scales of interest in the observation or model datasets. As shown by Feng et al. (2014) and Wu et al. (2016), combining EOF with EEMD to decompose the datasets leads to a faster calculation than MEEMD by 1 or 2 orders of magnitude because here the EEMD is applied to the temporal components (i.e., PCs) out of an EOF analysis rather than to all dimensions. Also, our EOF–EEMD package conducts additional HT–MSA and provides visualization of all spatial and temporal scales of interest.

• EOF analysis to decompose a two-dimensional dataset (time series at multiple locations) into spatial and temporal components.

  Suppose there are n locations, each having a time series of length p. The associated dataset Z is an n×p matrix. An EOF analysis of Z gives

  $$\begin{array}{}\text{(2)}& \mathbf{Z}=\mathbf{U}\sum {\mathbf{W}}^{\mathrm{T}}.\end{array}$$

  Here Σ is a diagonal q×q matrix containing the first q singular values of Z, and it represents the contribution of each pattern to the total variance of Z. The diagonal values of Σ are in a descending order, and thus the first several modes are the dominant ones. U is an n×q matrix representing the spatial component, and each column of U represents a spatial mode. W is a p×q matrix representing the temporal component, and each column of W represents a principal component (PC) for temporal variation associated with the corresponding spatial mode.

• EEMD analysis of each PC time series to obtain its “intrinsic mode functions” (IMFs) of descending frequencies.

  Each PC is mixed with multiple scales, which requires further decomposition in the time domain. Unlike fast Fourier transform (FFT) or wavelet transform (WT), EEMD does not need a priori bases, and it can be appropriately applied to delineate nonlinear and nonstationary time series, as in our pollution study.

  EEMD consists of an ensemble of empirical mode decomposition (EMD) performed on each PC time series (denoted as x(t) in Eq. 3). Each EMD linearly decomposes x(t) into individual IMFs c_j (of ascending timescales and descending frequencies) and a residual r_n:

  $$\begin{array}{}\text{(3)}& x\left(t\right)=\sum _{j=\mathrm{1}}^n{c}_j\left(t\right)+r_n\left(t\right).\end{array}$$

  EMD is based on finding the local maxima and minima of the time series. A detailed decomposition process can be found in Huang et al. (1998, 1999). EMD is much less susceptible to missing values and data interpolation than approaches that are based on an analysis of the whole time series (e.g., FFT and WT).

  EMD may be sensitive to noise in the real data to encounter a “mode mixing” problem (Wu et al., 2009). EEMD solves this problem by performing an ensemble of hundreds of EMDs, each with certain white noise added to x(t). Hence, the noise in the real data is incorporated as part of the white noise, and the ensemble further minimizes the effects of noise. The white noise is assumed to follow the standard Gaussian distribution (Wu et al., 2009). Figure 4 shows an example of the EEMD analysis.

• Hilbert transform and marginal spectrum analysis of each IMF to reveal its representative frequency range.

  There are no discrete periods or frequencies in the pollution and meteorological time series. Correspondingly, an IMF also has a continuous frequency range (rather than a constant frequency) that can be determined by HT–MSA. The HT reveals the IMF energy–frequency–time distribution (Huang et al., 1999). The MSA further shows the IMF distribution of variance (energy) with respect to different frequencies. The spectral peak represents the largest contribution to total variance.

  A spurious oscillation may occur near the edges of certain IMF time series, resulting in an inaccurate calculation of variance under HT–MSA. We apply a box-car filter (Gubbins, 2004) to select the internal 60 % of an IMF time series (from 20 % to 80 % of the 1488 h) to perform HT–MSA. Figure 4b shows an example of the visualized result of HT–MSA, in which the horizontal axis is the number of occurrences within the whole period (frequency, in h⁻¹, multiplied by the time length, 1488 h) and the vertical axis is the energy contribution. IMF2–IMF5 are visualized and analyzed in this study. The higher-frequency IMF1 is noisy as the energy is distributed over a wide range of occurrence numbers. IMF6–IMF10 represent the longest temporal scales that contribute little to the total variance of the decomposed PC. Thus IMF1 and IMF5–IMF10 are not further analyzed.

  Based on HT–MSA, we determine a representative frequency range (RFR) such that the range encompasses the peak frequency and that the frequencies within the range contribute 50 % of the total variance of an IMF. The frequencies below and above the RFR bounds each contribute 25 % of the total variance of the IMF. Before calculating the RFR, we smooth the marginal spectrum by connecting all local maxima of the spectrum with a cubic spline.

• Visualization of the spatial and temporal scales in a two-dimensional plot.

  Finally, we simultaneously visualize the spatial and temporal scales as well as their contributions to the total variance of Z in a two-dimensional plot for easy observational diagnosis and model evaluation. In this plot, an IMF is represented by a vertical “error bar” and a horizontal bar. The length of the error bar stands for the representative period range (RPR, the inverse of RFR), and a shorter length means a more stationary variation mode (i.e., towards a fixed frequency or period). The length of the horizontal bar stands for the contribution to the total variance. For clearer presentation, the plot does not include IMFs that do not pass the white noise examination, that lay outside the range of scales considered here (hours to days), or that contribute little to the total variance of the original data (e.g., less than 1 %).

3.1 General characteristics

Figure 4EEMD–HT–MSA result for PC1 of observed NO₂.

Download

Figure 5Observed (filled circles) and modeled (color maps) NO₂ and PM_2.5 averaged over 25 October–25 December 2013. Here the model results are averaged over all days rather than sampled at times of valid observations.

Download

The colored dots in Fig. 5a and b show the observed spatial distributions of city mean NO₂ and PM_2.5 averaged over the time period. Both NO₂ and PM_2.5 are largest over Beijing–Tianjin–Hebei (BTH) in the north and the Yangtze River Delta (YRD) in the east. NO₂ concentrations exceed 60 µg m⁻³ at many sites. The range of PM_2.5 is larger, from below 10 µg m⁻³ in some northern and coastal cities to about 200 µg m⁻³ in several cities of BTH.

Figure 6(a) Diurnal variation of observed NO₂ averaged over 25 October–25 December 2013. The black vertical bars represent 1 standard deviation across the days. PC1 from the EOF analysis is overlaid in red. (b) Similar to (a) but for PM_2.5. (c) Day-to-day variation of daily mean NO₂ over 25 October–25 December 2013. Data are de-trended. The black vertical bars represent 1 standard deviation due to the diurnal variation. PC1 and PC2 from the EOF analysis are overlaid in red. (d) Similar to (c) but for PM_2.5.

Download

Figure 6a and b show the diurnal variations of NO₂ and PM_2.5 over Eastern China, NEC, and SEC averaged over all days. Similarly, over the three regions, NO₂ peaks around 19:00 due to evening rush hour emissions, reduced PBL mixing, and a lengthened lifetime. NO₂ reaches a minimum at 14:00 because of the shortest lifetime and strongest PBL mixing. The diurnal range (maximum minus minimum) is about 30 µg m⁻³. The PM_2.5 level also reaches a minimum in the early afternoon. It has a much smaller diurnal range at 10 µg m⁻³. The vertical error bars in Fig. 6a and b depict the standard deviation for the day-to-day variation of NO₂ and PM_2.5 at any given hour. At a given hour, the PM_2.5 level is much more variable across the days than NO₂. In particular, the day-to-day standard deviation for PM_2.5 at a given hour is as large as the diurnal range of PM_2.5.

Figure 6c and d further show the time series of daily mean NO₂ and PM_2.5. All data are de-trended (trends are at 0.01 µg m⁻³ h⁻¹ for NO₂ and 0.05 µg m⁻³ h⁻¹ for PM_2.5). Although local maxima and minima (peaks and troughs of the time series) occur every several days, there is no single period or amplitude for the variation of each species. For NO₂ over Eastern China (black line in Fig. 6c), the local maxima vary from 60 to 100 µg m⁻³, and the local minima vary from 20 to 40 µg m⁻³. For PM_2.5 over Eastern China (black line in Fig. 6d), the local maxima vary from 100 to 300 µg m⁻³, and the local minima vary from 20 to 120 µg m⁻³. Furthermore, comparing the green and blue lines reveals that pollutants over NEC and SEC synchronize on some days but are out of phase on others; this feature is quantitatively analyzed in Sect. 3.2. These day-to-day variation patterns are associated with meteorological conditions and pollutant lifetimes.

Figure 7Daily anomalies of observed meteorological parameters and pollutant concentrations averaged over NEC and SEC, as well as their correlations. All data are de-trended. Correlation coefficients with “*” and “**” are statistically significant with P values below 0.05 and 0.01, respectively.

Download

Figure 7 shows day-to-day anomalies of observed pollutant concentrations and meteorological parameters over NEC and SEC. All data are de-trended. Over NEC, wind speed is clearly anticorrelated with pollutant levels. The correlation coefficient reaches −0.73 between NO₂ and wind speed and −0.60 between PM_2.5 and wind speed. Over this region, stronger winds are often associated with lower RH and lower temperature, characteristic of a cold air passage that brings cleaner, colder, and drier air from the north to NEC and transports the NEC pollution out of the region. Correspondingly, RH is strongly positively correlated with NO₂ (R=0.62) and PM_2.5 (R=0.69). The meteorology-associated day-to-day variability is more apparent after mid-November, when the variations of the two pollutants are more synchronous.

Over SEC (Fig. 7), the relationship between pollutant levels and meteorological parameters is more complex. The correlation between daily mean PM_2.5 and wind speed is relatively weak ($R=-\mathrm{0.44}$ compared to −0.60 over NEC), and its correlation with RH is even weaker (R=0.29). This indicates that the northerly air does not reduce PM_2.5 levels over SEC as effectively as over NEC, as PM_2.5 from NEC may be transported to SEC. By comparison, NO₂ is still highly anticorrelated with wind speed ($R=-\mathrm{0.77}$) over SEC, likely a result of the short lifetime of NO₂. Compared to PM_2.5 whose lifetime is sufficiently long (several days) for transport from NEC to SEC (Hu et al., 2014), NO₂ has a much shorter lifetime (below 1 day; Lin et al., 2012) and cannot undergo effective long-distance transport. However, almost all pollution measurement sites are urban, and weaker (stronger) winds allow for rapid accumulation (removal) of urban NO₂ pollution.

3.2 EOF–EEMD analyses of pollutants and meteorological parameters

Although informative, the time series analyses of regional mean pollution in Sect. 3.1 do not provide adequate quantitative information on the spatiotemporal variability and embedded scales. In fact, the separate discussion on NEC and SEC in Sect. 3.1 is largely inspired by the following EOF–EEMD analysis that suggests distinctive features between these two subregions. In this section, we use the EOF–EEMD package to distinguish and visualize the quantitative contributions of individual spatial and temporal modes to variations in the pollutant and meteorological data.

Figure 8EOF–EEMD–HT–MSA results for the observed temperature, RH, wind speed, NO₂, and PM_2.5. The first two rows depict EOF1 and EOF2, and the third row shows the EEMD–HT–MSA result for PC1 and PC2. In each panel of the third row, the length of the vertical “error bar” shows the RPR of an IMF, while the length of the horizontal bar represents the percentage contribution of the IMF to the total variance of the original data (as such, the horizontal lengths for different IMFs across different PCs can be compared). The blue (red) color indicates diurnal (day-to-day) variation.

Download

The columns in Fig. 8 show the EOF–EEMD results for the observed temperature, RH, wind speed, NO₂, and PM_2.5. The first two rows show the first two spatial patterns (EOF1 and EOF2) from the EOF analysis. The third row visualizes the EEMD–HT–MSA results for PC1 and PC2, the temporal counterparts of EOF1 and EOF2. For all variables, the first two PCs contribute more than 50 % of the total variance of the original data. The following PCs (PC3, PC4…) contain small variances and are not discussed here.

3.2.1 EOF–EEMD analyses of pollutants

The fourth column in Fig. 8 for NO₂ shows a primary pattern (EOF1 and PC1) with synchronous variation over the entirety of Eastern China. This pattern contributes 42 % of the total variance of NO₂. The two dominant IMFs of PC1 have time periods at 24 and 12 h, respectively, and together they contribute 30.4 % of the total variance of NO₂. Thus, PC1 mainly reflects the diurnal variation of NO₂. PC1 also contains some day-to-day variability in IMFs, which contribute about 10 % of the total variance of NO₂. The second pattern (EOF2) of NO₂ reveals opposite temporal variations between NEC and SEC. This temporal contrast is mainly reflected in the day-to-day variability, with RPRs around 2–5 days contributing 10.9 % of the total variance in NO₂. The day-to-day components of PC1 and PC2 correspond to the finding in Sect. 3.1 that NO₂ over NEC and SEC is synchronous on some days but out of phase on others.

We further investigate the physical meanings of PC1 and PC2 for NO₂. The red solid and red dashed lines in Fig. 6a and c show the diurnal and day-to-day variations of PC1 and PC2 in comparison to regional mean NO₂ levels over Eastern China (black line), NEC (green line), and SEC (blue line). Table 1 shows the associated correlation coefficients. PC1 is synchronous with Eastern China mean NO₂ for both diurnal and day-to-day variations (R reaches 1.0), confirming this regionally synchronous pattern. The day-to-day variation of PC2 is correlated with NEC NO₂ (R=0.66) but anticorrelated with SEC NO₂ ($R=-\mathrm{0.45}$, Table 1), again confirming this NEC–SEC contrasting pattern.

Table 1Correlation between PCs and regional mean values in terms of diurnal and day-to-day variability for NO₂.

** The correlation coefficient is statistically significant with the P value <0.01.

Download Print Version | Download XLSX

The last column in Fig. 8 shows the EOF–EEMD result for PM_2.5. As for NO₂, EOF1 and PC1 of PM_2.5 reflect a temporally synchronous pattern over Eastern China, which contributes 44 % of the total variation of PM_2.5. Again, PC1 is synchronous with Eastern China mean PM_2.5 (red versus black lines in Fig. 6b, d) in terms of both diurnal and day-to-day variations, with correlation coefficients approaching 1.0 (Table 2). However, the IMFs of PC1 representing diurnal variation are relatively weak, consistent with the noisy diurnal cycle of PM_2.5 discussed in Sect. 3.1. The dominant IMF of PC1 shows a period of around 7 days. PC2 of PM_2.5 reflects the day-to-day contrast between NEC and SEC (Fig. 6d and Table 2) with RPRs of 2–5 days, similar to PC2 of NO₂.

Table 2Correlation between PCs and regional mean values in terms of diurnal and day-to-day variability for PM_2.5.

** The correlation coefficient is statistically significant with the P value <0.01.

Download Print Version | Download XLSX

4.1 General evaluation

The color contours in Fig. 5a–d show the horizontal distributions of NO₂ and PM_2.5 simulated by GEOS-Chem and CMAQ. The model results here are averaged from all days over the time period rather than sampled from days with valid observations. Both models capture the general spatial patterns of observed NO₂ and PM_2.5, with the heaviest pollution over the north and east.

Figure 9Observed and simulated diurnal and day-to-day variations of (a) NO₂ and (b) PM_2.5 over NEC and SEC (µg m⁻³).

Download

Figure 9 evaluates the regional mean diurnal and day-to-day variations of modeled pollutant levels over NEC and SEC. Here model data are sampled from days and locations with valid observations. All trends are negligible and have been removed, consistent with the observational analysis. GEOS-Chem underestimates the observations by about 17 µg m⁻³ (21 µg m⁻³ over NEC and 13 µg m⁻³over SEC) for NO₂ and by 35 µg m⁻³ over Eastern China (31 µg m⁻³ over NEC and 41 µg m⁻³ over SEC) for PM_2.5 averaged over the whole period. The model bias is relatively consistent across individual hours. GEOS-Chem captures the observed diurnal variability for both pollutants as well as the day-to-day variability of PM_2.5, although it greatly underestimates the day-to-day variability of NO₂. More model evaluation statistics are shown in Table 3.

Table 3Observed and simulated pollutants and their correlations.

¹ Correlation between observed and simulated variables. ** indicates the correlation coefficient is statistically significant with the P value <0.01, while * indicates it passed a statistical test with the P value <0.01. ² Revised GEOS-Chem NO₂ and PM_2.5 by multiplying the ratio of the first layer to the second layer of CMAQ values.

Download Print Version | Download XLSX

Figure 9 also shows that WRF/CMAQ overestimates the nighttime observations by about 30 µg m⁻³ for NO₂ and 60 µg m⁻³ for PM_2.5 averaged over Eastern China, although it reproduces the daytime pollutant levels. This means an overestimate of the diurnal range, as is also revealed by the EOF–EEMD analysis in Sect. 4.2. CMAQ captures the day-to-day variability of daily mean NO₂ and PM_2.5 much better than GEOS-Chem (R=0.63–0.84 versus 0.25–0.37 over NEC and SEC for NO₂; and 0.87–0.88 versus 0.55–0.75 for PM_2.5). Note that the correlations shown here mainly reflect the model capabilities to capture Eastern China synchronous day-to-day variation, and they do not imply the model performance in simulating the NEC–SEC contrast, which is revealed in Sect. 4.2. More model evaluation statistics are shown in Table 3.

4.2 Model evaluation based on the EOF–EEMD analysis

Figure 10EOF–EEMD–HT–MSA results for observed, GEOS-Chem, and CMAQ NO₂. See Sect. 4.2 for detailed descriptions.

Download

Figure 11EOF–EEMD results for observed, GEOS-Chem, and CMAQ PM_2.5. See Sect. 4.2 for detailed descriptions.

Download

Figures 10 and 11 evaluate the EOF–EEMD results for modeled NO₂ and PM_2.5, respectively. Prior to the EOF–EEMD analysis, modeled NO₂ and PM_2.5 were sampled at times and locations with valid observations and then underwent the same interpolation procedure to fill in the missing values. In these figures, the last three rows visualize the EOF–EEMD–HT–MSA results in different ways (manifested in different lengths of the horizontal bar for each IMF). In the third row, the variance of each IMF is normalized to the total variance of the original data (NO₂ or PM_2.5). In the fourth row, the variance of each IMF is normalized to the variance of its respective PC in order to better visualize the signals from PC2 (which has a much smaller variance than PC1); as such, only the IMFs from the same PC are intercomparable. The fifth row visualizes the absolute variance of each IMF without any normalization.

The first two rows in Fig. 10 show EOF1 and EOF2 of NO₂. Both GEOS-Chem and CMAQ exhibit a synchronous pattern (EOF1) and an NEC–SEC contrasting pattern (EOF2), consistent with the observation. However, the CMAQ-simulated NEC–SEC contrast in EOF2 is much weaker than the observed. Table 1 shows that for modeled NO₂, PC1 is highly correlated with Eastern China mean NO₂ for diurnal (R=1.0 for GEOS-Chem and CMAQ) and day-to-day (R=0.56–0.97) variability and that PC2 is correlated with NEC NO₂ (R=0.74–0.81) and anticorrelated with SEC NO₂ ($R=-\mathrm{0.47}$ to −0.32) in terms of day-to-day variability, in line with the observational analysis.

The last three rows in Fig. 10 show that both models underestimate the contribution of day-to-day variability to the total variance of NO₂ (with a shorter length of horizontal bar). For PC1, CMAQ captures the RPR (position of “error bar”) and variance (length of horizontal bar) of the observed IMFs fairly well. By comparison, GEOS-Chem underestimates the day-to-day variance (too-small horizontal length) and does not capture its RPR. These results are consistent with the analysis in Sect. 4.1 (Fig. 9) showing that CMAQ is correlated with the observed Eastern China synchronous NO₂ time series much better than GEOS-Chem. For PC2, which reflects the NEC–SEC contrasting pattern, GEOS-Chem outperforms CMAQ in capturing the RPR and variance of the observed day-to-day IMFs (red colored in fourth row). This model characteristic is not seen from the time series discussion in Sect. 4.1.

Figure 11 shows that both GEOS-Chem and CMAQ capture the synchronous pattern (EOF1) and the NEC–SEC contrasting pattern (EOF2) of PM_2.5. For PC1, GEOS-Chem captures the variance of each IMF but not its RPR (especially for the day-to-day IMFs). CMAQ simulates too-strong diurnal IMFs, consistent with its overestimated diurnal cycle discussed in Sect. 4.1. CMAQ outperforms GEOS-Chem in capturing the RPR of day-to-day IMFs of PC1, in line with its better correlation with the observations (Fig. 9). For PC2, GEOS-Chem captures the variance and RPR of the observed day-to-day IMFs better than CMAQ.

4.3 Discussion on model deficiencies

WRF/CMAQ overestimates the diurnal variation of NO₂ and PM_2.5. The causes are multifaceted. The ACM2 PBL mixing scheme in WRF v3.5.1 and CMAQ v5.0.1 (used here) assumes the same value for the eddy diffusivity of momentum (K_m) and heat (K_h), which implies a Prandtl number ($Pr=K_{\mathrm{m}}/K_{\mathrm{h}}$) of unity and too-weak mixing under stable atmospheric conditions (i.e., at night). This deficiency has been alleviated in WRF v3.7 and CMAQ v5.1. Also, there is inconsistency between CMAQ and WRF in the Monin–Obukhov length in the surface layer module. This error has been corrected in CMAQ v5.1. For more model update details, please refer to the online document (https://www.airqualitymodeling.org/index.php/CMAQ_version_5.1_(November_2015_release)_Technical_Documentation, last access: 3 September 2018).

GEOS-Chem (the first model layer) underestimates surface NO₂ by about 17 µg m⁻³ and PM_2.5 by 35 µg m⁻³ averaged over Eastern China. The underestimate of PM_2.5 is in part because this simulation of GEOS-Chem does not include secondary organic aerosols, which likely contribute as much as 21 % of PM_2.5 over Eastern China (Fu et al., 2012). Also, the model does not include anthropogenic dust. Furthermore, although the observation stations are close to the ground, the first layer of GEOS-Chem is too thick (130 m) to fully capture the vertical gradient of pollution concentrations. Figure 12 shows Eastern China mean vertical profiles of NO₂ in the two models. The center of the first layer of CMAQ (40 m) is closer to the ground, and the center of its second layer is located at a height similar to the center of the first layer of GEOS-Chem. CMAQ shows a strong vertical gradient of NO₂ from its first to second layer. Had we used the CMAQ-simulated ratio of the first over the second layer to extrapolate GEOS-Chem first-layer NO₂ to 40 m, this would significantly increase the model's “ground-level” NO₂ (by 24 % over NEC and 17 % over SEC) and PM_2.5 (by 45 % and 17 %). However, the extrapolation does not improve the day-to-day correlation to the observations, indicating the important roles played by other factors. See Table 3 for more evaluation statistics.

Figure 12Eastern China mean NO₂ vertical profiles simulated by GEOS-Chem and CMAQ averaged over 25 October–25 December 2013. The black and red dots denote the center of each vertical layer in the two models. The evening is from 20:00 to 23:00 LT, while the afternoon is from 12:00 to 15:00 LT.

Download

GEOS-Chem (the first model layer) also underestimates the Eastern China synchronous day-to-day variation of NO₂. When averaged over the 10 lowest layers (below 850 hPa), GEOS-Chem NO₂ captures the day-to-day variability of observed surface NO₂. This suggests that the model deficiency in day-to-day variability may be specific to the first layer. Moreover, the first layer of GEOS-Chem captures the day-to-day variation of observed NO₂ in the afternoon (12:00–15:00 LT, R=0.9 over NEC and 0.8 over SEC), but the model performance is rather poor in the evening (20:00–23:00 LT, R=0.1 over NEC and SEC), suggesting nighttime-specific model inadequacies. A further analysis of nighttime ozone and the NO:NO₂ ratio suggests that GEOS-Chem greatly underestimates the observed nighttime ozone by 49.2 % on average over NEC and 54.6 % over SEC, particularly on days when its NO:NO₂ ratio is much greater than the CMAQ-modeled ratio. The mean NO:NO₂ ratio in GEOS-Chem is 1.8 over NEC and 1.4 over SEC, greater than the ratio in CMAQ (1.0 over NEC and 0.4 over SEC) by a factor of 2–3. Overall, it appears that the nighttime chemistry is poorly represented in the first layer of GEOS-Chem, the causes of which warrant further investigations.

The magnitude of emission differences between the two models plays an insignificant role in the differences between their simulated NO₂ or PM_2.5 concentrations. The Chinese anthropogenic emissions in 2010 used in GEOS-Chem (except for NO_x) are close to the emissions in 2013 used in CMAQ (within 10 % for both gases and primary aerosols, mostly within 5 %; see Zheng et al., 2018). NO_x emissions in GEOS-Chem are scaled to 2013 using satellite NO₂ data, which further eliminates the differences from those used in CMAQ. The difference in the spatial distribution of emissions is also small (Geng et al., 2017; Zheng et al., 2018).

We further use CMAQ simulations to investigate whether the inclusion of SOA affects our analysis of the spatiotemporal patterns of PM_2.5. Supplement Fig. S1 compares the time series of CMAQ-simulated PM_2.5 with versus without including SOA. Although SOA contributes about 8–9 µg m⁻³ of PM_2.5 averaged over the days, inclusion of SOA does not affect the temporal variability. The EOF–EEMD results in Supplement Fig. S2 further confirm that the spatiotemporal scales are very consistent whether or not SOA is included.