Simulating the Effects of Urban Parameterizations on the Passage of a Cold Front During a Pollution Episode in Megacity Shanghai

Urbanization significantly influences meteorological conditions and air quality. Statistically, air pollution in the megacity of Shanghai usually occurs with cold weather fronts. An air pollution episode during a cold front was simulated using weather research and forecasting and the Community Multi-scale Air Quality model system. In this study, we used two urban schemes, a simple bulk scheme (denoted BULK) and the single-layer urban canopy model (SLUCM), to check the effects of urban parameterization. Our results showed that SLUCM better predicted the arrival time and cooling process of the cold front and more realistically simulated the moving process of the cold front. The improvement in the temperature and relative humidity simulation achieved using SLUCM was more effective under higher urbanization levels, whereas the wind speed simulation was better in rural areas. The simulated concentrations at sites with high urbanization were obviously improved by urban parameterization. The barrier role of the urban canopy during a cold front was better represented and was shown to cause a wider polluted area and higher pollutant concentration using SLUCM than with BULK. Overall, accurate meteorological simulations in the atmospheric boundary layer using SLUCM are expected to provide good prediction of urban air quality.


Introduction
Shanghai is located in the eastern section of the Yangtze River Delta, which is one of the three largest city-cluster areas in China, with a population of 20 million.Shanghai City has experienced rapid economic development and urban expansion.Since the 1980s and according to the national new urbanization report in 2015, Shanghai's urbanization rate ranked first, at 87.6% in 2015 [1,2].Shanghai has experienced a significant change in population and land use and land cover change, which has resulted in buildings displacing cultivated land [3].Land use types and building density have dramatically changed in Shanghai since the 1980s.Both observations and simulations proved that urban surface change affects the meteorological conditions [4,5].Urbanization converts the city surface to urban covers including tall buildings, thereby increasing the roughness of the surface over the urban area.The friction and drag of the buildings decrease the near-surface wind speed in the urban area.This aerodynamic effect has been observed in many cities [6,7].The influence of urbanization on meteorological fields has been observed in the Yangtze River Delta [8].Zhang et al. [9] documented that urbanization may cause a wind speed loss of more than 50% in the Yangtze River Delta over urbanized areas, and that this wind speed decrease may not only occur over urban areas but also on a regional scale.
Applying an urban parameterization scheme in meteorological models is a key requirement to accurately simulate the atmospheric boundary layer (ABL) [10][11][12].Kusaka et al. [13] used the weather research and forecasting (WRF) model coupled with a single-layer canopy scheme to simulate the Tokyo heat island status.Teddy and Pullen [14] performed a 23-day temperature analysis using the WRF-urban canopy model (UCM) in the New York region.Zhang et al. [15] performed a heat island analysis using the WRF-UCM in the Suzhou area in China.Miao et al. [16] simulated the influence of the morphological characteristics of buildings on the urban boundary layer structure of Beijing using an urban boundary layer model.These efforts have shown that a mesoscale model with a canopy module has some advantages for simulating weather and fluxes in urban areas, and particularly the ability to better simulate the characteristics of an urban heat island.Although these studies identified the meteorological mechanisms related to urban schemes, the advantages and disadvantages of the urban schemes have still not been clearly identified.Additionally, the performance of urban schemes over megalopolises such as Shanghai is still uncertain.
Air quality is simultaneously affected by meteorological conditions and emission sources.The effects of meteorological conditions are known to closely interact with chemical compositions, chemical reaction processes, and physical-optical characteristics of air pollutants [17][18][19][20][21][22][23].Emissions from the same polluter could cause the concentration of air pollutant to differ by dozens or even a hundred times in different weather conditions.Therefore, meteorological field simulation limits the precision of the air quality simulation.The use of wind, temperature, and ABL height (ABLH) is necessary to understand the transport and diffusion of air pollutants within the boundary layer.Fan et al. [24] simulated the structure and characteristics of the ABL over the Pearl River Delta and noted that ABL characteristics are important for identifying air pollution problems.Therefore, accurately simulating the ABL characteristics will help to better predict the urban air quality in Shanghai.
Cold front system are important synoptic systems involved in heavy air pollution in China.These systems produce a difference in the wind ahead of and behind a cold front, and changes in ABL height (ABLH) and vertical convection, which have significant impacts on regional accumulation and transportation of air pollutants [25].Shahgedanova et al. [26] determined that the passage of a cold front causes high concentrations of pollutants in areas where the elevated sources were determined by analyzing the NO 2 concentration of multiple land use types and pollutant sources.Cannichael et al. [27] simulated the impact of a low-pressure system on the transfer process of ozone in spring time in East Asia and showed that cold front processes greatly influence the spatial distribution of ozone.Studies are lacking on the impact of cold front processes on air pollution focusing on China, especially for the megacity of Shanghai.Wang et al. [28] researched the influence of cold fronts on the air quality of Lanzhou by combining historical synoptic maps with data from observation sites.Cheng et al. [25] used the Community Multiscale Air Quality (CMAQ) and HYbrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) models to simulate the air pollution processes of a typical strong cold front moving across Northeastern China.
We aimed to evaluate the ability of different urban schemes to simulate the characteristics of the ABL in the Shanghai metropolitan area.Section 2 describes the model settings.The observational data and synoptic situations are presented in Section 3. Comparisons of the measured and simulated data are presented in Section 4. Finally, the conclusions are provided in Section 5.All the time used in this study was Beijing time.

Model Simulation and Scheme Settings
The WRF mesoscale model was employed in this study.The WRF model consists of fully compressible nonhydrostatic equations on a staggered Arakawa C grid.Its vertical coordinate is a terrain-following hydrostatic pressure coordinate [29].A three-dimensional (3D) Eulerian atmospheric chemistry and transport model, CMAQ model version 4.6, was applied to simulate the concentrations of the major pollutants.Designed as a "one-atmosphere" model, CMAQ can simultaneously address the complex couplings among several air quality issues across spatial scales ranging from local to hemispheric [30,31].The emission inventory of Shanghai was developed with a 1 × 1 km resolution by the Shanghai Environmental Monitoring Center (SEMC) [32], in which the emission in the Baoshan district was updated by Tan et al. [33].For areas outside Shanghai, the emission inventory was combined with the Regional Emission Inventory in Asia with Transport Chemical evolution over Pacific (TRACE-P).Detailed emission information was provided by Tan et al. [33].

Urban Schemes
The initial urban parameterization scheme of the WRF model is a bulk scheme (denoted BULK) [34] that uses a roughness length of 0.8 m to represent the effects of urban surfaces, a surface albedo of 0.15 to represent the radiation trapped in the urban canyons, a volumetric heat capacity of 3.0 MJ m −3 K −1 , and a thermal conductivity of 3.24 W m −1 K −1 to represent large heat storage by the urban buildings and roads.This approach was employed in real-time forecasts by Liu et al. [35].Therefore, detailed land use data may affect the results of the BULK scheme predictions of the ABL characteristics.
The single-layer urban canopy model (SLUCM) is another urban parameterization in WRF, developed by Kusaka et al. [36] and Kusaka and Kimura [37].The urban geometry is represented through infinitely long street canyons, and three different urban surfaces (roofs, walls, and roads) are recognized.Shadowing, reflections, and radiation trapping in the street canyon are considered, and an exponential wind profile is prescribed to deduce the wind speed in the canyon from the wind speed above the canyon, where the lowest grid point is located in the model.All the urban surfaces are categorized as roofs, walls, or roads.The sensible heat fluxes from the roofs, walls, and roads are introduced in the lowest atmospheric layer [38].
Two simulation schemes were carried out to drive WRF: (1) the BULK urban option and (2) the SLUCM urban option.Both of the simulations were verified by measurements at four observational sites (Figure 1).concentrations of the major pollutants.Designed as a "one-atmosphere" model, CMAQ can simultaneously address the complex couplings among several air quality issues across spatial scales ranging from local to hemispheric [30,31].The emission inventory of Shanghai was developed with a 1 × 1 km resolution by the Shanghai Environmental Monitoring Center (SEMC) [32], in which the emission in the Baoshan district was updated by Tan et al. [33].For areas outside Shanghai, the emission inventory was combined with the Regional Emission Inventory in Asia with Transport Chemical evolution over Pacific (TRACE-P).Detailed emission information was provided by Tan et al. [33].

Urban Schemes
The initial urban parameterization scheme of the WRF model is a bulk scheme (denoted BULK) [34] that uses a roughness length of 0.8 m to represent the effects of urban surfaces, a surface albedo of 0.15 to represent the radiation trapped in the urban canyons, a volumetric heat capacity of 3.0 MJ m −3 K −1 , and a thermal conductivity of 3.24 W m −1 K −1 to represent large heat storage by the urban buildings and roads.This approach was employed in real-time forecasts by Liu et al. [35].Therefore, detailed land use data may affect the results of the BULK scheme predictions of the ABL characteristics.
The single-layer urban canopy model (SLUCM) is another urban parameterization in WRF, developed by Kusaka et al. [36] and Kusaka and Kimura [37].The urban geometry is represented through infinitely long street canyons, and three different urban surfaces (roofs, walls, and roads) are recognized.Shadowing, reflections, and radiation trapping in the street canyon are considered, and an exponential wind profile is prescribed to deduce the wind speed in the canyon from the wind speed above the canyon, where the lowest grid point is located in the model.All the urban surfaces are categorized as roofs, walls, or roads.The sensible heat fluxes from the roofs, walls, and roads are introduced in the lowest atmospheric layer [38].
Two simulation schemes were carried out to drive WRF: (1) the BULK urban option and (2) the SLUCM urban option.Both of the simulations were verified by measurements at four observational sites (Figure 1).

Model Settings
In this study, the coarse domain covered the area between 20.4 and 41.3° N and 110.1 and 134.2°E, with a horizontal resolution of 27 × 27 km, whereas the inner most domain covered the area between 30.35 and 31.52°N and 120.52 and 122.07°E, and it had a horizontal resolution of 1 × 1 km.

Model Settings
In this study, the coarse domain covered the area between 20.4 and 41.3 • N and 110.1 and 134.2 • E, with a horizontal resolution of 27 × 27 km, whereas the inner most domain covered the area between 30.35 and 31.52 • N and 120.52 and 122.07 • E, and it had a horizontal resolution of 1 × 1 km.The boundary layer processes were resolved using 28 vertical sigma levels from the top to the surface level, in which 17 layers were less than 2 km.The primary physics options included the WRF single-moment 6-class microphysics scheme [39], Betts-Miller-Janjic cumulus parameterization, the Rapid Radiative Transfer Model (RRTM) long-wave radiation scheme, the Dudhia shortwave radiation scheme, the Eta similarity surface layer scheme, the Noah land surface model, and the Mellor-Yamada-Janjic boundary layer scheme.For the two inner domains (D03 and D04), convection was assumed to be well-resolved by the explicit microphysical parameterization scheme, and no cumulus parameterization scheme was used.
A summer pollution period, occurring during 1-8 July 2011, was simulated in this study.The primary focus was the ABL characteristics during a cold front period on 4-5 July 2011.The lateral and initial conditions for the WRF simulations in this study were obtained from a daily NCEP (National Centers for Environmental Prediction)/NCAR (national center for atmospheric research) reanalysis with 1 × 1 grid data.In the present study, the spin-up time was 12 h.
The land use required for the WRF Preprocessing System (WPS) in WRF was updated using the interpretation data from Shanghai aerial remote sensing images from 2007.According to the national land use classification (for trial use), the data were divided into 71 categories, and the file of geo-grid sub-module was modified by matching the 71 categories with U.S. Geological Survey (USGS)-33 [40].Figure 2 shows the land use from 1980 USGS that was default packed in the WRF model and the updated land use for Shanghai City in 2007 used in this study.A comparison of the data showed that the urban area has considerable increased over the three decades from 1980 until 2007.In this simulation, three different urban classes (commercial or industrial (COI), high-intensity residential (HIR), and low-intensity residential (LIR) areas) and USGS 33-category land cover, derived from land use databases, were used (e.g., the National Land Cover Data for the United States (NLCD), developed by the USGS).Shanghai urban land was linked to three different urban classes by the impervious surface percentage: 20-49% defines LIR, 50-79% defines HIR, and 80-100% defines COI.Two different urban schemes were set in parallel simulations: (1) WRF with BULK and (2) WRF with SLUCM.Table 1 lists the morphological parameters considered for the three urban classes in the SLUCM option.Table 2 lists the thermal parameters for the SLUCM option.level, in which 17 layers were less than 2 km.The primary physics options included the WRF singlemoment 6-class microphysics scheme [39], Betts-Miller-Janjic cumulus parameterization, the Rapid Radiative Transfer Model (RRTM) long-wave radiation scheme, the Dudhia shortwave radiation scheme, the Eta similarity surface layer scheme, the Noah land surface model, and the Mellor-Yamada-Janjic boundary layer scheme.For the two inner domains (D03 and D04), convection was assumed to be well-resolved by the explicit microphysical parameterization scheme, and no cumulus parameterization scheme was used.
A summer pollution period, occurring during 1-8 July 2011, was simulated in this study.The primary focus was the ABL characteristics during a cold front period on 4-5 July 2011.The lateral and initial conditions for the WRF simulations in this study were obtained from a daily NCEP (National Centers for Environmental Prediction)/NCAR (national center for atmospheric research) reanalysis with 1 × 1 grid data.In the present study, the spin-up time was 12 h.
The land use required for the WRF Preprocessing System (WPS) in WRF was updated using the interpretation data from Shanghai aerial remote sensing images from 2007.According to the national land use classification (for trial use), the data were divided into 71 categories, and the file of geo-grid sub-module was modified by matching the 71 categories with U.S. Geological Survey (USGS)-33 [40].Figure 2 shows the land use from 1980 USGS that was default packed in the WRF model and the updated land use for Shanghai City in 2007 used in this study.A comparison of the data showed that the urban area has considerable increased over the three decades from 1980 until 2007.In this simulation, three different urban classes (commercial or industrial (COI), high-intensity residential (HIR), and low-intensity residential (LIR) areas) and USGS 33-category land cover, derived from land use databases, were used (e.g., the National Land Cover Data for the United States (NLCD), developed by the USGS).Shanghai urban land was linked to three different urban classes by the impervious surface percentage: 20-49% defines LIR, 50-79% defines HIR, and 80-100% defines COI.Two different urban schemes were set in parallel simulations: (1) WRF with BULK and (2) WRF with SLUCM.Table 1 lists the morphological parameters considered for the three urban classes in the SLUCM option.Table 2 lists the thermal parameters for the SLUCM option.λ is the thermal conductivity of the material, C is the specific heat of the material, T(int) is the initial temperature of the material and also the temperature of the deepest layer, ε is the emissivity of the surface, α is the albedo of the surface, and z 0 is the roughness length for momentum over the surface.

Statistical Synoptic Situations During Air-Pollution Days in 2011-2013
There are no uniform standards to judge pollution days from 2011 to 2013 since the Air Quality Index (AQI) was introduced in 2012.In this study, we defined days with a daily average concentration of PM 2.5 of 115 µg/m 3 or more as pollution days and those with more than 150 µg/m 3 as heavy pollution days.Based on these definitions, 83 pollution days were selected from the three years and classified into four kinds of weather pattern: equalizing pressure, saddle shaped field pattern, westerly in rear of high-pressure pattern, and cold front pattern.Detailed pollution days for different weather patterns are shown in Table 3.This shows that the probability of heavy pollution caused by cold front weather is higher than the other weather types.The cold front event that occurred in July 2011, a typical summer cold front that is difficult to predict in summer seasons, is important for simultaneously forecasting the meteorological field and pollution conditions.

Observational Data for the Simulated Case
We chose a pollution episode caused by a cold front to simulate the effects of urbanization on the meteorological conditions and concentrations of air pollutants.The test period occurred from Beijing time 8:00 a.m. on 1 July 2011 to 8:00 a.m.local time (UTC + 8) on 8 July 2011, including a pollution period.The four observational sites for meteorological conditions are shown in Figure 1b, including Baoshan (s1), Fudan (s2), Hongqiao (s3), and Pudong (s4).The observational data included hourly wind speed, wind direction, temperature, and relative humidity for the four sites.Site s1 was located in the Baoshan Meteorological Bureau, a national benchmark weather station in Northern Shanghai.At s1, radio soundings were performed to determine the mean velocity, wind direction, and temperature.Sounding balloons were launched two times per day (at 8:00 a.m. and 8:00 p.m. UTC + 8).Station s1 represents LIR areas.Station s2 was located on the roof of a 10-m building at Fudan University, a mixed-pollution urban site with residential areas and trafficked roads in Northeastern Shanghai.At s2, vertical measurements were conducted with light detection and ranging (LIDAR) technology.Station s2 represents HIR areas.Sites s3 and s4 were surface automatic observation stations.Station s3 is located at Hongqiao Airport, approximately in the middle of Shanghai.Station s4 was located at Pudong Airport, in Eastern Shanghai.Station s3 represents LIR areas, and s4 represents a rural environment.To effectively represent the spatiality of geography, we chose measurement sites that could widely cover the Shanghai domain.The temperature in other sites at Minghan, Chomgming, Nanhui, Qingpu, Chongming, Jinshan, and Fengxian were used to show the time evolution of the temperature field of the cold front.The ambient concentrations of air pollutants were obtained from monitoring sites named Baoshan, Xuhui, Minhang, Songjiang, Qingpu, Pudongbinhai, Pudong liuli, and Fengxian.

Weather and Pollution Situations for the Simulated Case
During the period of 1-8 July, the average temperature in Shanghai was approximately 31 • C, and the overall weather conditions were hot and rainy on some days (Figure 3).A cooling trend occurred on 4-5 July, caused by a cold front arriving in Shanghai, as shown in Figure 4. Severe convective weather appears to have promoted the spread of pollutants.Figure 5 shows the hourly concentrations of SO 2 , NO X , and PM 2.5 in Shanghai on 1-8 July.The normal diurnal variation of the three pollutants reached a peak at 8:00 a.m. and 8:00 p.m. each day.However, a special diurnal variation occurred on 4-5 July, caused by the cold front.This variation was characterized by an extremely low concentration during 2:00 p.m. July 4-12:00 a.m.July 5 followed by an extremely high concentration during 2:00-8:00 p.m. on July 5 (except for SO 2 ).In our study, during 10:00 a.m. to 3:00 p.m. on 4 July, the cold front was passing over Shanghai from the north of Shanghai, where the Baoshan site was located.The SO 2 concentration in Baoshan site had a particularly high concentration during 2:00-5:00 p.m. on 5 July, and a similar high concentration was not recorded at other sites.To test the ability of the WRF model to predict the characteristics of the ABL using different urban schemes, this period, 1-8 July, was used as a case study to analyze the passing of a cold front.

Weather and Pollution Situations for the Simulated Case
During the period of 1-8 July, the average temperature in Shanghai was approximately 31 °C, and the overall weather conditions were hot and rainy on some days (Figure 3).A cooling trend occurred on 4-5 July, caused by a cold front arriving in Shanghai, as shown in Figure 4. Severe convective weather appears to have promoted the spread of pollutants.Figure 5 shows the hourly concentrations of SO2, NOX, and PM2.5 in Shanghai on 1-8 July.The normal diurnal variation of the three pollutants reached a peak at 8:00 a.m. and 8:00 p.m. each day.However, a special diurnal variation occurred on 4-5 July, caused by the cold front.This variation was characterized by an extremely low concentration during 2:00 p.m. July 4-12:00 a.m.July 5 followed by an extremely high concentration during 2:00-8:00 p.m. on July 5 (except for SO2).In our study, during 10:00 a.m. to 3:00 p.m. on 4 July, the cold front was passing over Shanghai from the north of Shanghai, where the Baoshan site was located.The SO2 concentration in Baoshan site had a particularly high concentration during 2:00-5:00 p.m. on 5 July, and a similar high concentration was not recorded at other sites.To test the ability of the WRF model to predict the characteristics of the ABL using different urban schemes, this period, 1-8 July, was used as a case study to analyze the passing of a cold front.

Comparison of Statistical Parameters
To validate the WRF simulation results, the model results were compared with the available surface observations furnished by the meteorological stations in Domain 4. Some common statistical variables were used to evaluate the performance of the model simulations, including the mean bias (MB), mean absolute gross error (MAGE), root mean squared error (RMSE), fractional absolute error (FAE), and the correlation coefficient (r).MB, MAGE, RMSE, FAE and r are defined in Table 4 [41], where M is the number of stations, N is the number of numerical hours used in the calculation, excluding the spin-up time, and C m and C o represent the modeled and observed values, respectively.From the definitions given above, as MB, MAGE, RMSE, and FAE approach zero, the model simulations approach the observed values.These four statistics primarily reflect the accuracy of the simulated values, and the correlation coefficient (r) primarily reflects the quality of the simulated trend.The statistics for the wind direction were calculated by finding the actual angle between the simulated and observed wind directions.The rule is if

Comparison of Statistical Parameters
To validate the WRF simulation results, the model results were compared with the available surface observations furnished by the meteorological stations in Domain 4. Some common statistical variables were used to evaluate the performance of the model simulations, including the mean bias (MB), mean absolute gross error (MAGE), root mean squared error (RMSE), fractional absolute error (FAE), and the correlation coefficient (r).MB, MAGE, RMSE, FAE and r are defined in Table 4 [41], where M is the number of stations, N is the number of numerical hours used in the calculation, excluding the spin-up time, and C m and C o represent the modeled and observed values, respectively.From the definitions given above, as MB, MAGE, RMSE, and FAE approach zero, the model simulations approach the observed values.These four statistics primarily reflect the accuracy of the simulated values, and the correlation coefficient (r) primarily reflects the quality of the simulated trend.The statistics for the wind direction were calculated by finding the actual angle between the simulated and observed wind directions.The rule is if Table 5 shows the mean statistical parameters for four sites between the measured and the simulated 2-m air temperature (T), wind speed (WS), wind direction (WD), relative humidity (RH), and ABLH on 1-8 July 2011.For the air temperature, wind speed, relative humidity, and ABLH, the statistics on the MB, MAGE, RMSE, and FAE were all similar, i.e., the four statistics (indicated by X) of the two schemes behave as |X BULK | > |X SLUCM | > 0, indicating that the SLUCM simulation performed better than the BULK simulation for all four sites' measures.The total number of observations for all sites was 768.The model mostly underestimated the air temperature, a phenomenon that may be caused by systematic errors of the WRF model during the extreme high-temperature period.The values of MB for the BULK and SLUCM schemes were −1.95 and −1.13 • C, respectively, indicating that the mean simulated T increased by 0.8 • C in SLUCM, which was closer to the observations than the BULK value.Scheme SLUCM decreased the MAGE by 0.6 • C, the RMSE by 0.6 • C, and the FAE by 0.03, showing that the simulated values of SLUCM were better than those of BULK.The correlation coefficient of the temperature for each scheme was greater than 0.7, which shows that the simulated and observed values of the temperature were well correlated.The relative humidity showed similar statistical features to those of air temperature, except that the model mostly overestimated the relative humidity.The model mostly overestimated the wind speed.The values of MB for the BULK and SLUCM schemes were 1.75 and 1.08 m s −1 , respectively, indicating that SLUCM predicted a mean simulated speed that was lower by 0.7 m s −1 and that the SLUCM simulations were closer to the observations.Scheme SLUCM decreased the MAGE by 0.51 m s −1 , the RMSE by 0.48 m s −1 , and the FAE by 0.08, showing that SLUCM performed better than BULK.The predictions of changing wind speed trends by BULK were better than those by SLUCM because the correlation coefficients of the BULK and SLUCM schemes were 0.49 and 0.33, respectively.
The calculated MB values for the wind direction using the BULK and SLUCM schemes were 10.17 • and 15.41 • .Therefore, SLUCM increased the mean simulated wind direction by 5.24 • , deviating further from the observations.Scheme SLUCM increased the MAGE by 2.31 • , the RMSE by 5.05 • and the FAE by 0.01, indicating that the BULK simulation performed better than the SLUCM simulation.Notably, the difference in FAE between the BULK and SLUCM simulations was small, at only 0.01.The correlation coefficients of the two schemes were 0.47 and 0.41, indicating that BULK better predicts the changing trends in wind direction.All statistical parameters showed that BULK is better than SLUCM for wind direction simulation; however, the difference between BULK and SLUCM was minimal.The reason for this difference may be that although SLUCM considers the urban morphological parameters in principle, the accuracy of the parameters may be inadequate.Default values were used for parameters such as ε, α, and z 0 in Table 2 because of the lack of local data.This lack may cause errors that could exceed the improvement in the wind simulation on the ground by SLUCM.
The ABLH is a critical parameter for vertical dispersion.The WRF model underestimated the ABLH in most simulations.Using the SLUCM scheme increased the average ABLH by up to 37 m, resulting in simulations closer to the observations.Scheme SLUCM decreased the MAGE by 40.78 m, the RMSE by 17.86 m, and the FAE by 0.04, indicating that the SLUCM simulation performed better than the BULK simulation.The correlation coefficients of the two schemes were 0.52 and 0.59, respectively, indicating that using the SLUCM scheme resulted in better predictions for the changing trends of the ABLH.

Comparison of Measurements and Simulations for Different Urban Classes
To further show the effects of urbanization parameterizations in different urban classes, the simulations were compared with the measurements in different urban-level areas.Table 6 shows the statistical values of three different urban classes for the measured and the simulated parameters.Generally, the SLUCM scheme provided a greater improvement in the temperature, relative humidity, and wind speed simulations in areas of higher urbanization than in lower urbanization and rural areas.As shown in Table 6 the MB values for temperature calculated for the BULK and SLUCM schemes in LIR areas were −2.21 and −1.39 • C, respectively, and the correlation coefficients were 0.85 and 0.81, respectively.Scheme SLUCM increased the mean simulated temperature by 0.82 • C in LIR areas.In HIR areas, the MB values of the BULK and SLUCM schemes in HIR areas were −2.29 and −1.31 • C, respectively, and the correlation coefficients were 0.56 and 0.66, respectively.Thus, SLUCM increased the mean simulated temperature by 0.98 • C and the correlation coefficient by 0.1 in HIR areas.There was a greater improvement in the temperature simulation by SLUCM in the HIR areas than in the LIR areas (Table 6).
For relative humidity in HIR areas, as shown in Table 6, the MB values of the BULK and the SLUCM schemes in HIR areas were 8.86% and 4.63%, respectively, and the correlation coefficients were 0.50 and 0.66, respectively.Thus, SLUCM decreased the mean simulated relative humidity by 4.23% and increased the correlation coefficient by 0.16 in HIR areas.There was a greater improvement in the relative humidity simulation using SLUCM in HIR areas than in LIR areas.
As shown in Table 6, the MB values for wind speeds of BULK and SLUCM the HIR areas were 1.13 and 0.54 m/s, respectively, and the correlation coefficients were 0.49 and 0.40, respectively.Scheme SLUCM decreased the mean simulated wind speed by 0.59 m/s and the correlation coefficient by 0.09 in the HIR areas.In the rural areas, the MB values of BULK and SLUCM in the rural areas were 1.43 and 0.57 m/s, respectively, and the correlation coefficients were 0.39 and 0.31, respectively.Scheme SLUCM decreased the mean simulated wind speed by 0.86 m/s and the correlation coefficients by 0.08 in the rural areas.Generally, the model performed better in the HIR areas.Scheme SLUCM improved the simulation of wind speed but did not improve the simulation of wind speed trends.The observed data revealed that the ABLH on non-cold front days exhibited the well-known diurnal variation.As shown in Figure 5, the concentration of pollutants (except for PM 2.5 ) remained at low levels between 8:00 a.m. and 8:00 p.m. and at high values between 8:00 p.m. and 8:00 a.m. the next day during the same period, indicating that the pollutant concentration was correlated well with ABLH.For the ABLH, the model performed better during low-value periods than in high-value periods of diurnal variation.The duration of the simulated high values was low, and the values fluctuated significantly during the high-value periods.Scheme SLUCM showed no superiority over BULK during the regular diurnal variations.At 8:00 a.m. on 4 July, the cold front arrived, and the ABLH began to decrease until it reached a nadir at approximately 12:00 p.m.Then, it remained stable at approximately 500 m until 2:00 p.m. on 5 July.Between 2:00 and 8:00 p.m. on 5 July, a peak occurred at approximately 1000 m, far less than the normal daily peak at 1500 m; then, the ABLH was stable again at approximately 500 m until 8:00 a.m. on July 6.Scheme SLUCM accurately predicted the times that the ABLH began to decrease on 4 July and increase on 6 July, and it approximately predicted the peak between 2:00 and 8:00 p.m. on 5 July.The BULK scheme predicted that the ABLH would decrease four hours earlier and increase two hours earlier, and it predicted the ABLH peak on 5 July as occurring six hours earlier.From 8:00 p.m. on 5 July to 2:00 a.m. on 6 July, BULK predicted a virtual trough near 0 m, but SLUCM predicted a constant at approximately 500 m, which was more consistent with the observed values.To compare the detailed vertical wind field during the cold front, the vertical wind profile obtained by radio soundings is shown with the model results at station s1 in Figure 7.The selected time is the moment of the greatest change in ABLH during the cold front.There were nulls at some heights because the data from the radio soundings below 2000 m were low-resolution.The simulated wind field from the ground to a high altitude showed a change from a northeasterly wind to a westerly wind at 8:00 a.m. on 4 July.Near the surface (0-100 m), the observed wind speed was low, and the wind direction was southwesterly.The simulated wind direction was northeasterly, and BULK overestimated the wind speed.At approximately 800 m, the observed wind direction changed to a westerly wind.The BULK-simulated wind was still in the process of changing direction because the wind speed approached 0 m/s, but the SLUCM-simulated wind direction had already changed to a westerly wind with the wind speed being closer to the measured wind speed.This result indicates The MAGEs of BULK and SLUCM were 472.73 m and 431.95 m, respectively, and the correlation coefficients were 0.52 and 0.59, respectively.Scheme SLUCM decreased the MAGE by 40.78 m and increased the correlation coefficient by 0.07.Between 2:00 p.m. on 4 July and 8:00 a.m. on 6 July, the ABLH remained low.Between 2:00 p.m. on 4 July and 2:00 a.m. on 5 July, i.e., the first 10 h after the cold front arrived, the pollutant concentration remained low.The reason for this low concentration may be that the arrival of the cold front brought strong convective weather.From 2:00 a.m. on 5 July until 8:00 a.m. on 6 July, the pollutant concentration began to increase, and the NO x and PM 2.5 levels remained high.Therefore, the pollutant concentration was correlated well with the ABLH during a cold front period, but strong wind and precipitation brought by a cold front may reduce pollution.In general, for the ABLH simulation, SLUCM performed better than BULK during the cold front period but did not show any superiority during the conventional diurnal variation period.Additionally, the pollutant concentration correlated well with the ABLH, indicating that SLUCM might perform better during high pollution periods.

Comparison of Atmospheric Boundary Layer Height and Vertical Profile of Wind
To compare the detailed vertical wind field during the cold front, the vertical wind profile obtained by radio soundings is shown with the model results at station s1 in Figure 7.The selected time is the moment of the greatest change in ABLH during the cold front.There were nulls at some heights because the data from the radio soundings below 2000 m were low-resolution.The simulated wind field from the ground to a high altitude showed a change from a northeasterly wind to a westerly wind at 8:00 a.m. on 4 July.Near the surface (0-100 m), the observed wind speed was low, and the wind direction was southwesterly.The simulated wind direction was northeasterly, and BULK overestimated the wind speed.At approximately 800 m, the observed wind direction changed to a westerly wind.The BULK-simulated wind was still in the process of changing direction because the wind speed approached 0 m/s, but the SLUCM-simulated wind direction had already changed to a westerly wind with the wind speed being closer to the measured wind speed.This result indicates that, at this moment, SLUCM better predicted the vertical change in the wind direction.At 8:00 p.m. on 4 July, both the observed and simulated wind directions changed from an easterly wind on the ground to a westerly wind at high altitudes.Between 200 and 800 m, the observed wind entirely finished the shift in direction, and SLUCM predicted this process well.At 800 m, the wind direction of SLUCM completely changed to a westerly wind, and the wind speed was close to the observed value.At this time and height, the wind direction predicted by BULK was still easterly.The SLUCM scheme was better than the BULK scheme at 8:00 p.m. on 4 July.At 8:00 a.m. on 5 July, the observed wind direction from the ground to a high altitude changed from easterly to westerly.The two schemes did not perform well at this moment in time.The predicted height of the change in wind direction was lower than the observation, and the simulated wind speed at the high altitude was larger than the observed speed.At 8:00 p.m. on 5 July, the wind changed from a southerly wind on the ground to a westerly wind at a high altitude.Neither scheme reproduced the change in wind direction well at this moment.Scheme BULK overestimated the wind speed, and SLUCM performed better than BULK.In general, for vertical wind simulations between 200 and 1000 m, SLUCM performed better than BULK.The data in Figure 7 and Table 3 show that SLUCM simulated the wind better than BULK for the 200-1000 m vertical heights in urban areas; however, on the ground, BULK performed better.The reason for this difference may be that the urban parameters are not sufficiently comprehensive.Although SLUCM considers the overall urban morphology, the available parameters were not accurate enough to represent the actual urban situation.Default values were used for parameters such as ε, α, and z 0 in Table 2 due to the lack of local data.This may have caused errors that exceeded the improvement in the wind simulation on the ground due to SLUCM.However, the improvement that SLUCM adds to the simulation is apparent as the errors from the buildings disappeared between 200 and 1000 m.Overall, SLUCM is recommended in urban areas to effectively simulate the wind field above the urban canopy, which shows the city morphology has a strong impact on the vertical wind field simulation of WRF.Given sufficiently accurate land use data, SLUCM may simulate the wind field in the entire boundary layer better than BULK.
improvement that SLUCM adds to the simulation is apparent as the errors from the buildings disappeared between 200 and 1000 m.Overall, SLUCM is recommended in urban areas to effectively simulate the wind field above the urban canopy, which shows the city morphology has a strong impact on the vertical wind field simulation of WRF.Given sufficiently accurate land use data, SLUCM may simulate the wind field in the entire boundary layer better than BULK.

Comparison of the Surface Temperature Fields
Figure 8 shows the observed air temperature at the 11 sites in and the simulated air temperature field over Shanghai at 10:00 and 13:00 on 4 July 2011.The seven sites (except s1, s2, s3, and s4) shown in Figure 8 were used as supplemental sites for comparison of the temperature recordings.The observed temperatures at Chongming, s1, s2, and s4 began to decrease at 10:00 a.m. on 4 July, and the cooling process continued until 3:00 p.m. when the temperature dropped to 26 °C.The observed temperatures at all sites showed that the cold front arrived north of Shanghai during the time period between 10:00 and 11:00 a.m. and that the temperatures then decreased from the north to the south of Shanghai, with the cooling continuing for five hours.The BULK scheme predicted a significant frontal surface at 10:00 and 11:00 a.m..At 1:00 p.m., the cold front covered all of Shanghai, and the temperature did not decrease significantly after 1:00 p.m.At Chongming, Jiading, s1, and s2, the temperature began to decrease at 10:00 a.m..The cold front simulated by SLUCM arrived north of Shanghai at 10:00 a.m.; however, the edge of the cold front was not clear.During the time period between 10:00 and 11:00 a.m., the temperature north of Shanghai began to decrease, and cooling persisted for more than four hours, i.e., until 3:00 p.m.At Chongming, Jiading, and s1, the temperature began to decrease at 10:00 a.m.. SLUCM more realistically simulated the time of the arrival of the cold front, thus predicting the movement of the cold front.
Figure 9 shows that the time series of temperatures during the cold front passed over Shanghai at 11 sites divided into two areas: (1) urban sites containing s1, s2, s3, and Minhang; and (b) suburban

Comparison of the Surface Temperature Fields
Figure 8 shows the observed air temperature at the 11 sites in Shanghai and the simulated air temperature field over Shanghai at 10:00 and 13:00 on 4 July 2011.The seven sites (except s1, s2, s3, and s4) shown in Figure 8 were used as supplemental sites for comparison of the temperature recordings.The observed temperatures at Chongming, s1, s2, and s4 began to decrease at 10:00 a.m. on 4 July, and the cooling process continued until 3:00 p.m. when the temperature dropped to 26 • C. The observed temperatures at all sites showed that the cold front arrived north of Shanghai during the time period between 10:00 and 11:00 a.m. and that the temperatures then decreased from the north to the south of Shanghai, with the cooling continuing for five hours.The BULK scheme predicted a significant frontal surface at 10:00 and 11:00 a.m.At 1:00 p.m., the cold front covered all of Shanghai, and the temperature did not decrease significantly after 1:00 p.m.At Chongming, Jiading, s1, and s2, the temperature began to decrease at 10:00 a.m.The cold front simulated by SLUCM arrived north of Shanghai at 10:00 a.m.; however, the edge of the cold front was not clear.During the time period between 10:00 and 11:00 a.m., the temperature north of Shanghai began to decrease, and cooling persisted for more than four hours, i.e., until 3:00 p.m.At Chongming, Jiading, and s1, the temperature began to decrease at 10:00 a.m.SLUCM more realistically simulated the time of the arrival of the cold front, thus predicting the movement of the cold front.
Figure 9 shows that the time series of temperatures during the cold front passed over Shanghai at 11 sites divided into two areas: (1) urban sites containing s1, s2, s3, and Minhang; and (b) suburban sites containing s4, Nanhui, Qingpu, Chongming, Jinshan, and Fengxian.The cold front simulated by BULK cooled significantly faster than what was observed.The spatial field of temperatures simulated by SLUCM is much closer to the observed temperatures.Additionally, the advantage of SLUCM on temperature simulation was more obvious in urban areas.The change in temperatures was closer to the observed temperatures and the absolute error was smaller.The cooling process simulated by SLUCM was closer to the observations.In general, SLUCM simulated the cooling process better than

Comparison of the Concentrations of Air Pollutants
In terms of the effect of meteorological conditions on pollutant concentrations, we used NOx as an example to show the influence of urban scheme on the simulation of air quality.The observed and simulated mean concentrations of NOx at eight stations on 1-8 July 2011 are shown in Figure 10.The seven sites (except Baoshan) shown in Figure 11 were used as supplemental sites for a comparison of the concentrations of air pollutants.This shows that different urban schemes influence NOx.At most sites, SLUCM simulations performed better than the BULK scheme.the simulated concentrations at urban sites like Xuhui and Minhang were obviously improved by WRF-SLUCM-CMAQ.The NOx concentration determined by SLUCM (59.3 μg/m 3 at Xuhui and 60.9 μg/m 3 at Minhang) was much closer to the observations (57.4 μg/m 3 at Xuhui and 68.4 μg/m 3 at Minhang) than BULK (49.4 μg/m 3 at Xuhui and 46.2 μg/m 3 at Minhang).

Comparison of the Concentrations of Air Pollutants
In terms of the effect of meteorological conditions on pollutant concentrations, we used NOx as an example to show the influence of urban scheme on the simulation of air quality.The observed and simulated mean concentrations of NO x at eight stations on 1-8 July 2011 are shown in Figure 10.The seven sites (except Baoshan) shown in Figure 11 were used as supplemental sites for a comparison of the concentrations of air pollutants.This shows that different urban schemes influence NO x .At most sites, SLUCM simulations performed better than the BULK scheme.the simulated concentrations at urban sites like Xuhui and Minhang were obviously improved by WRF-SLUCM-CMAQ.The NOx concentration determined by SLUCM (59.3 µg/m 3 at Xuhui and 60.9 µg/m 3 at Minhang) was much closer to the observations (57.4 µg/m 3 at Xuhui and 68.4 µg/m 3 at Minhang) than BULK (49.4 µg/m 3 at Xuhui and 46.2 µg/m 3 at Minhang).
The simulated spatial field and observed NO x concentrations at eight sites in Shanghai with the arrival of the cold front on 4 July 2011 at 10:00 a.m. and 3:00 p.m. are shown in Figure 11 to demonstrate the response of air pollutants to WRF-SLUCM/CMAQ.WRF-CMAQ of BULK and SLUCM both simulated the delayed arrival of the cold front leading to a slow response to the change in air pollutants with the cold front.The simulated concentration of NO x in the area surrounding the Baoshan site was still high at 10:00 a.m. but the observations dropped below 50 µg/m 3 .However, the prediction is more realistic in that the large point-source pollutants in the north of Shanghai diffused more slowly toward the south with the SLUCM scheme compared to with the BULK scheme.At 3:00 to 4:00 p.m., the simulated concentration of NO x was dropping in some sites, like the Minhang site, in WRF-BULK-CMAQ; however, the simulated concentration of NO x was still in the high value area in WRF-SLUCM-CMAQ, which is closer to the observation of 150 µg/m 3 at the Minhang site.Overall, the reduction in wind speed and the boundary layer due to the barrier role of the urban area that led to a wider polluted urban area and higher concentration was demonstrated by WRF-SLUCM-CMAQ compared with WRF-BULK-CMAQ, which is closer to the observation during the cold front.The simulated spatial field and observed NOx concentrations at eight sites in Shanghai with the arrival of the cold front on 4 July 2011 at 10:00 a.m. and 3:00 p.m. are shown in Figure 11 to demonstrate the response of air pollutants to WRF-SLUCM/CMAQ.WRF-CMAQ of BULK and SLUCM both simulated the delayed arrival of the cold front leading to a slow response to the change in air pollutants with the cold front.The simulated concentration of NOx in the area surrounding the Baoshan site was still high at 10:00 a.m. but the observations dropped below 50 μg/m 3 .However, the prediction is more realistic in that the large point-source pollutants in the north of Shanghai diffused more slowly toward the south with the SLUCM scheme compared to with the BULK scheme.At 3:00 to 4:00 p.m., the simulated concentration of NOx was dropping in some sites, like the Minhang site, in WRF-BULK-CMAQ; however, the simulated concentration of NOx was still in the high value area in WRF-SLUCM-CMAQ, which is closer to the observation of 150 μg/m 3 at the Minhang site.Overall, the reduction in wind speed and the boundary layer due to the barrier role of the urban area that led to a wider polluted urban area and higher concentration was demonstrated by WRF-SLUCM-CMAQ compared with WRF-BULK-CMAQ, which is closer to the observation during the cold front.Figure 12 shows that the time series of NO x concentrations on 3-5 July.The results in Figure 12 show that the observed concentrations rose right before the arrival of the cold front and then declined continuously during the cold front.Due to the northerly wind, the cold front was coming from the north of Shanghai, where the Baoshan site is located.NO x concentration declined with the cold front arrival at 10:00 a.m.Simultaneously, the area influenced by the northerly winds expanded with the moving of the cold front.The Xuhui site was directly downstream from the cold front at this moment; the concentration of NO x there increased and continued to increase to 124 µg/m 3 at 2:00 p.m. Subsequently, NO x concentration decreased during the cold front passing the Xuhui site.Figure 13 shows vertical cross sections of simulated NO x concentrations along the longitude 121.4 • E during the process of the cold front passing.It reveals again that surface NO x tended toward high concentrations before the cold front, significantly decreasing as the cold front passed, and had higher concentrations after the cold front had passed by.The results are similar to those reported by Cheng et al. [25].At 10:00 a.m., the cold front passed the Baoshan site with the Minhang site downstream from the cold front.The high NO x concentration lifted to the upper atmosphere, which appeared to form hotspots around 500 m at the Baoshan site, while the concentration rose at the Minhang site.At 12:00-1:00 p.m., the cold front covered Shanghai, the surface NO x concentration rapidly decreased, and high NO x concentrations lifted to the upper air layer.After the cold front passed the Xuhui site, the NO x concentration increased to a higher value at around 2:00-3:00 p.m.Both WRF-SLUCM-CMAQ and WRF-BULK-CMAQ simulated similar vertical variations with the passage of a cold front, but the results using WRF-SLUCM-CMAQ were much more accurate.Figure 12 shows that the time series of NOx concentrations on 3-5 July.The results in Figure 12 show that the observed concentrations rose right before the arrival of the cold front and then declined continuously during the cold front.Due to the northerly wind, the cold front was coming from the north of Shanghai, where the Baoshan site is located.NOx concentration declined with the cold front arrival at 10:00 a.m.Simultaneously, the area influenced by the northerly winds expanded with the moving of the cold front.The Xuhui site was directly downstream from the cold front at this moment; form hotspots around 500 m at the Baoshan site, while the concentration rose at the Minhang site.At 12:00-1:00 p.m., the cold front covered Shanghai, the surface NOx concentration rapidly decreased, and high NOx concentrations lifted to the upper air layer.After the cold front passed the Xuhui site, the NOx concentration increased to a higher value at around 2:00-3:00 p.m.Both WRF-SLUCM-CMAQ and WRF-BULK-CMAQ simulated similar vertical variations with the passage of a cold front, but the results using WRF-SLUCM-CMAQ were much more accurate.

Conclusions
In this study, the BULK and SLUCM schemes were implemented in WRF-CMAQ to analyze their simulation performance of urban meteorology and pollutant changes for the entire Shanghai area.For meteorological aspects including air temperature, wind speed, relative humidity, and ABLH, the simulations performed by the SLUCM scheme were better than the BULK scheme.The SLUCM scheme simulations were closer to the actual observations than the BULK scheme simulations.The SLUCM scheme simulations decreased the mean wind speed by 0.7 m s −1 , increased the mean temperature by up to 0.8 °C, and increased the mean ABLH by up to 37 m.For both schemes, the simulated and observed values of temperature and relative humidity were well correlated.The

Conclusions
In this study, the BULK and SLUCM schemes were implemented in WRF-CMAQ to analyze their simulation performance of urban meteorology and pollutant changes for the entire Shanghai area.For meteorological aspects including air temperature, wind speed, relative humidity, and ABLH, the simulations performed by the SLUCM scheme were better than the BULK scheme.The SLUCM scheme simulations were closer to the actual observations than the BULK scheme simulations.The SLUCM scheme simulations decreased the mean wind speed by 0.7 m s −1 , increased the mean temperature by up to 0.8 • C, and increased the mean ABLH by up to 37 m.For both schemes, the simulated and observed values of temperature and relative humidity were well correlated.The SLUCM scheme better simulated changing trends to predict the ABLH but did not accurately simulate the wind speed or wind direction.The SLUCM simulation was better for the air temperature, wind speed, relative humidity, and ABLH, but did not work well for wind direction.For the concentration field of air pollutants, the simulation results using the SLUCM scheme performed comparatively better than those using the BULK scheme.The SLUCM scheme predicted a wider area of pollution than the BULK scheme, and the concentrations were closer to the values observed using the SLUCM scheme.
The SLUCM simulation performed significantly better than the BULK simulation for the daily minimum temperature predictions, especially in urban areas.However, for the daily maximum temperature simulations, SLUCM did not perform better than BULK.The improvement in the temperature simulation achieved using SLUCM was more effective at higher urbanization levels.The SLUCM scheme performed better for the mean relative humidity simulation, especially in LIR areas.However, the improvement in the relative humidity simulation by the SLUCM scheme was greater in the HIR areas than in the LIR areas.The implementation of SLUCM reduced the virtual peak that occurred in the wind speed simulation.The greatest improvement in the wind speed simulation by the implementation of SLUCM occurred in the rural areas.However, the SLUCM scheme did not show more merit than the BULK scheme for wind direction simulations in the rural and HIR areas.However, the SLUCM scheme performed better in the LIR areas.For the ABLH simulations, the SLUCM scheme performed better than the BULK scheme during a cold front episode but was not better for conventional diurnal variations.The SLUCM scheme better simulated wind than the BULK scheme for the vertical layers between 200 and 1000 m in urban areas, but on the ground, the BULK scheme performed better.The SLUCM scheme exhibited advantages when the errors from buildings disappeared between 200 and 1000 m.Therefore, in urban areas, SLUCM is the recommended scheme to apply with localized urban parameters to effectively simulate the wind field above the urban canopy.During the passage of the cold front, SLUCM predicted the arrival time of the cold front and more realistically simulated the moving cold front.The cooling process simulated by SLUCM was closer to the observation than that simulated by BULK.The barrier role of the urban canopy during a cold front was better represented by SLUCM than by BULK.
For the pollutant concentration field, the mean values predicted by the SLUCM scheme were closer to the observations at urban sites.The barrier role of the urban area led to a wider polluted urban area and a higher concentration being predicted by SLUCM than by BULK, which is closer to the observations during the cold front.
The limitations of the urban parameters in Shanghai can limit the simulation accuracy of the SLUCM scheme.Given sufficiently accurate land use data, the SLUCM scheme can simulate the wind field in the entire boundary layer better than the BULK scheme, which then affects the precision of the pollutant concentration simulation.

Figure 1 .
Figure 1.(a) Domain settings in the weather research and forecasting (WRF) model and (b) the measurement sites in Shanghai.

Figure 1 .
Figure 1.(a) Domain settings in the weather research and forecasting (WRF) model and (b) the measurement sites in Shanghai.

Figure 2 .
Figure 2. The land use of Shanghai in (a) 1980 and (b) 2007.

Table 1 .
Urban morphological parameters considered for the three urban classes: commercial or industrial (COI), high-intensity residential (HIR), and low-intensity residential (LIR).

Figure 2 .
Figure 2. The land use of Shanghai in (a) 1980 and (b) 2007.

Figure 4 .
Figure 4. Weather maps of East Asia at 2:00 a.m. and 2:00 p.m. UTC + 8 on 4 July and 5 July 2011.The maps are courtesy of the Korea Meteorological Administration.

Figure 4 .
Figure 4. Weather maps of East Asia at 2:00 a.m. and 2:00 p.m. UTC + 8 on 4 July and 5 July 2011.The maps are courtesy of the Korea Meteorological Administration.

Figure 6
Figure 6 shows the observed and simulated ABLH values at station s2 on 1-8 July 2011.The observed data revealed that the ABLH on non-cold front days exhibited the well-known diurnal variation.As shown in Figure5, the concentration of pollutants (except for PM 2.5 ) remained at low levels between 8:00 a.m. and 8:00 p.m. and at high values between 8:00 p.m. and 8:00 a.m. the next day during the same period, indicating that the pollutant concentration was correlated well with ABLH.For the ABLH, the model performed better during low-value periods than in high-value periods of diurnal variation.The duration of the simulated high values was low, and the values fluctuated significantly during the high-value periods.Scheme SLUCM showed no superiority over BULK during the regular diurnal variations.At 8:00 a.m. on 4 July, the cold front arrived, and the ABLH began to decrease until it reached a nadir at approximately 12:00 p.m.Then, it remained stable at approximately 500 m until 2:00 p.m. on 5 July.Between 2:00 and 8:00 p.m. on 5 July, a peak occurred at approximately 1000 m, far less than the normal daily peak at 1500 m; then, the ABLH was stable again at approximately 500 m until 8:00 a.m. on July 6.Scheme SLUCM accurately predicted the times that the ABLH began to decrease on 4 July and increase on 6 July, and it approximately predicted the peak between 2:00 and 8:00 p.m. on 5 July.The BULK scheme predicted that the ABLH would decrease four hours earlier and increase two hours earlier, and it predicted the ABLH peak on 5 July as occurring six hours earlier.From 8:00 p.m. on 5 July to 2:00 a.m. on 6 July, BULK predicted a virtual trough near 0 m, but SLUCM predicted a constant at approximately 500 m, which was more consistent with the observed values.

Figure 7 .
Figure 7. Vertical wind profile and ABLH measured by radio soundings and model results at station s1 at 8:00 and 20:00 on 4 July and 5 July 2011.

Figure 7 .
Figure 7. Vertical wind profile and ABLH measured by radio soundings and model results at station s1 at 8:00 and 20:00 on 4 July and 5 July 2011.
passage of the cold front.The reason for this improvement might be that the urban canopy layer increases the surface roughness of the Shanghai metropolitan area, slowing the advance of the cold front.Without considering the urban canopy, BULK simulated a more rapidly cooling cold front, and SLUCM simulated the cooling closer to the actual measurements by considering the urban canopy.Atmosphere 2018, 9, x FOR PEER REVIEW 14 of 22 sites containing s4, Nanhui, Qingpu, Chongming, Jinshan, and Fengxian.The cold front simulated by BULK cooled significantly faster than what was observed.The spatial field of temperatures simulated by SLUCM is much closer to the observed temperatures.Additionally, the advantage of SLUCM on temperature simulation was more obvious in urban areas.The change in temperatures was closer to the observed temperatures and the absolute error was smaller.The cooling process simulated by SLUCM was closer to the observations.In general, SLUCM simulated the cooling process better than BULK during the passage of the cold front.The reason for this improvement might be that the urban canopy layer increases the surface roughness of the Shanghai metropolitan area, slowing the advance of the cold front.Without considering the urban canopy, BULK simulated a more rapidly cooling cold front, and SLUCM simulated the cooling closer to the actual measurements by considering the urban canopy.

Figure 12 .
Figure 12.Time series of NOx concentration at Xuhui site on 3-5 July 2011.Red dash line: Cold front passing time.

Figure 12 .
Figure 12.Time series of NO x concentration at Xuhui site on 3-5 July 2011.Red dash line: Cold front passing time.

Figure 13 .
Figure 13.Vertical distribution of simulated NOx concentrations along the longitude 121.4°E during the process of the cold front passing.

Figure 13 .
Figure 13.Vertical distribution of simulated NO x concentrations along the longitude 121.4 • E during the process of the cold front passing.

Table 1 .
Urban morphological parameters considered for the three urban classes: commercial or industrial (COI), high-intensity residential (HIR), and low-intensity residential (LIR).

Table 2 .
Thermal parameters used in the single-layer urban canopy model (SLUCM) for every urban class.

Table 3 .
Pollution days for different weather patterns.
Station s2 represents HIR areas.Sites s3 and s4 were surface automatic observation stations.Station s3 is located at Hongqiao Airport, approximately in the middle of Shanghai.Station s4 was located at Pudong Airport, in Eastern Shanghai.Station s3 represents LIR areas, and s4 represents a rural environment.To effectively represent the spatiality of geography, we chose measurement sites that could widely cover the Shanghai domain.The temperature in other sites at Minghan, Chomgming, Nanhui, Qingpu, Chongming, Jinshan, and Fengxian were used to show the time evolution of the temperature field of the cold front.The ambient concentrations of air pollutants were obtained from monitoring sites named Baoshan, Xuhui, Minhang, Songjiang, Qingpu, Pudongbinhai, Pudong liuli, and Fengxian.

Table 4 .
Calculation formulae for the primary statistical parameters.

Table 4 .
Calculation formulae for the primary statistical parameters.

Table 5 .
The mean statistical parameters of four sites between the measured and the simulated 2-m air temperature (T), wind speed (WS), wind direction (WD), relative humidity (RH), and ABL height (ABLH, station s2 measured data only) on 1-8 July 2011.

Table 6 .
Statistical parameters of three different areas between the measured and simulated 2-m air temperature (T), relative humidity (RH), and wind speed (WS) on 1-8 July 2011.