Next Article in Journal
Reproductive Effects of S. boulardii on Sub-Chronic Acetamiprid and Imidacloprid Toxicity in Male Rats
Next Article in Special Issue
Modeling Impacts of Urbanization on Winter Boundary Layer Meteorology and Aerosol Pollution in the Central Liaoning City Cluster, China
Previous Article in Journal
Editorial for the Special Issue on “Human Biomonitoring in Health Risk Assessment: Current Practices and Recommendations for the Future”
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mitigation Effect of Dense “Water Network” on Heavy PM2.5 Pollution: A Case Model of the Twain-Hu Basin, Central China

1
Hubei Meteorological Service Center, Wuhan 430205, China
2
China Meteorological Administration Basin Heavy Rainfall Key Laboratory/Hubei Key Laboratory for Heavy Rain Monitoring and Warning Research, Institute of Heavy Rain, China Meteorological Administration, Wuhan 430205, China
3
Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing 210044, China
4
Jiangsu Climate Center, Nanjing 210044, China
5
Key Laboratory of Meteorology and Ecological Environment of Hebei Province, Hebei Provincial Institute of Meteorological Sciences, Shijiazhuang 050021, China
6
State Key Laboratory of Pollution Control and Resource Reuse and School of the Environment, Nanjing University, Nanjing 210023, China
*
Authors to whom correspondence should be addressed.
Toxics 2023, 11(2), 169; https://doi.org/10.3390/toxics11020169
Submission received: 2 January 2023 / Revised: 7 February 2023 / Accepted: 8 February 2023 / Published: 10 February 2023

Abstract

:
The influence of the underlying surface on the atmospheric environment over rivers and lakes is not fully understood. To improve our understanding, this study targeted the Twain-Hu Basin (THB) in central China, with a unique underlying surface comprising a dense “water network” over rivers and lakes. In this study, the Weather Research and Forecasting Model with Chemistry (WRF-Chem) was used to simulate the impact of this dense “water network” on a wintertime heavy PM2.5 pollution event in the THB. On this basis, the regulating effects of density and area of the lake groups, with centralized big lakes (CBLs) and discrete small lakes (DSLs), on PM2.5 concentrations over the underlying surface of the dense “water network” in the THB were clarified, and the relative contributions of thermal factors and water vapor factors in the atmospheric boundary layer to the variation of PM2.5 concentrations were evaluated. The results show that the underlying surface of dense “water networks” in the THB generally decreases the PM2.5 concentrations, but the influences of different lake-group types are not uniform in spatial distribution. The CBLs can reduce the PM2.5 concentrations over the lake and its surroundings by 4.90–17.68% during the day and night. The ability of DSLs in reducing PM2.5 pollution is relatively weak, with the reversed contribution between −5.63% and 1.56%. Thermal factors and water vapor–related factors are the key meteorological drivers affecting the variation of PM2.5 concentrations over the underlying surface of dense “water networks”. The warming and humidification effects of such underlying surfaces contribute positively and negatively to the “purification” of air pollution, respectively. The relative contributions of thermal factors and water vapor–related factors are 52.48% and 43.91% for CBLs and 65.96% and 27.31% for DSLs, respectively. The “purification” effect of the underlying surface with a dense “water network” in the THB on regional air pollution highlights the importance of environmental protection of inland rivers and lakes in regional environmental governance. In further studies on the atmospheric environment, long-term studies are necessary, including fine measurements in terms of meteorology and the environment and more comprehensive simulations under different scenarios.

1. Introduction

The underlying surface of river and lake systems is an ecosystem dominated by hydrological processes, and the water–heat exchange between water and atmosphere is an important mechanism affecting the atmospheric boundary layer and physicochemical processes of aerosols and has important implications for regional weather, climate, and the environment [1,2,3,4,5,6].
The underlying surface of “water networks” has remarkable influences on the vertical structure of temperature and humidity in the atmosphere and on the physical processes in the atmospheric boundary layer [7,8]. Compared with other types of underlying surfaces, water bodies have greater heat capacity and thermal inertia, which causes slower warming during the day and less sensible heat exchange between water and atmosphere [9]. At the same time, more latent heat is released over the water surface due to evaporation. The inhomogeneous properties of land and water surface lead to thermal differences between the two types of surfaces and thus change the thermal radiation. The changes of thermal radiation and heat exchange intensity directly affect the process of heat transfer near the surface, which further change the structure of temperature and humidity in the surrounding areas and cause changes in the structure of the local atmospheric boundary layer [10]. The temperature and humidity structures in the atmospheric boundary layer are important factors for the accumulation of air particulate matter and affect the diffusion and deposition of air pollutants [11,12,13,14,15]. The local circulation of lake–land breeze driven by heat factors will change the temperature and wind fields and other meteorological conditions in the surrounding areas, which also affects the diffusion and distribution of air pollutants [5,16,17,18,19]. During the day, under the circulation of lake breeze, air pollutant particles are transported to the land and diffused aloft by the updraft of lake breeze. The lake breeze circulation can also cause air pollutant deposition on the lakeshore, where the pollutants maintain a high concentration along the lakeshore for a long time [20]. At night, the concentration of air particulate matter over the land decreases. Under the joint influence of environmental wind and land breeze, the air particulate matter accumulates over the lake region [21]. In addition, urban rivers are conducive to the generation of secondary aerosols, which are potential sources of atmospheric organic matter [22]. At the same time, the water surface is favorable for the hygroscopic growth of fine particles in the atmosphere and hence affects the distribution of air pollutants, which thus plays a crucial role in the distribution and diffusion of regional atmospheric particulate matter [23,24,25].
The evolution process of air pollution can be fully understood by using life cycle assessment (LCA) [26] and the Weather Research and Forecasting Model with Chemistry (WRF-Chem) [27]. The WRF-Chem model (version 3.9.1) developed by the National Oceanic and Atmospheric Administration (NOAA) and National Center was used to simulate a heavy pollution event in the Twain-Hu Basin (THB). This model can couple complex physical and chemical processes such as transport, sedimentation, emissions, chemical conversions, aerosol effect, photolysis, and radiation of chemical substances and has been widely used in the simulation of air pollution [28,29,30,31].
The Twain-Hu Basin (THB), located in the middle reaches of the Yangtze River Basin of central China, has typical characteristics of a monsoon climate. It is a region with abundant water resources in river and lake systems. The area of water bodies such as rivers and lakes account for more than 20% of the plain area in the THB, which is the region with the largest share of water area in China [32]. The dense river network and numerous lakes make the THB region a typical inland humid region in China and form a humid monsoon climate.
In recent years, the THB has become a region of frequent incidence of air particulate pollution in China [33,34,35]. The underlying surface of dense “water networks” in the THB makes the spatiotemporal variations of air particulate matter and the mechanisms involved more complex, requiring further exploration. However, the studies on the influences of lake water-driven thermal circulation on the local atmospheric environment have been mainly focused on coastal cities and large lakes in plateau regions [36,37,38,39,40]; there are few studies in inland areas. Given the important regulating effect of the unique underlying surface of “water networks” on the variation of regional PM2.5 concentrations, clarifying the influence on the regional atmospheric environment can deepen the scientific understanding of the mechanism for air pollution variation in the inland humid area of China.
This study focused on a wintertime PM2.5 pollution event in the THB in central China and explored the influence of dense “water networks” in the THB on the PM2.5 variation and the mechanism involved. The results are expected to provide a reference for the research and prevention of air pollution in inland humid areas as well as a decision-making basis for the research and protection of urban underlying surfaces. The remainder of this paper is organized as follows. Section 2 introduces the datasets and model setup used in the simulation. Section 3 presents the main results, and the conclusions are presented in Section 4.

2. Materials and Methods

2.1. Materials

In this study, a heavy PM2.5 pollution event induced by stagnant weather conditions in the THB during 11–24 January 2018 was selected as the research object. The Weather Research and Forecasting Model (WRF) with Chemistry (WRF-Chem) was used to conduct a series of simulations. The Final Operational Global Analysis (FNL) data provided by the National Center for Environmental Prediction (NECP)/National Center for Atmospheric Research (NCAR) was used to adopt the initial meteorological field and boundary conditions for WRF-Chem model simulations. The spatial and temporal resolutions of the FNL data are 1° × 1° and 6 h, respectively. Surface meteorological variables such as 10 m wind speed and sea level pressure were also used in this study, which are from the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis (ERA5) dataset (https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land?tab=form, accessed on 7 February 2023), with the spatiotemporal resolutions being 0.25° × 0.25° and 1 h. The hourly surface observation data obtained from the National Meteorological Center (http://data.cma.cn/, accessed on 7 February 2023) were also adopted, including temperature, humidity, wind speed and direction, atmospheric pressure, and precipitation. The hourly dataset of PM2.5 mass concentrations was obtained from the China Air Quality Online Monitoring and Analysis Platform (https://www.aqistudy.cn/, accessed on 7 February 2023), which currently collects the PM2.5 concentration data from 367 cities in China, and all the data are automatically updated every hour. The gridded dataset of emission sources adopts the 2016 Multi-resolution Emission Inventory for China (MEIC; http://meicmodel.org.cn/, accessed on 7 February 2023) [41,42], which is output by month, with the spatial resolution being 0.25° × 0.25°.
This study also updated the underlying surface data for the WRF-Chem with the MODIS Land Cover Type (MCD12Q1) product. This product adopts the storage mode of sinusoidal projection and contains multiple classification schemes, which describe land cover properties derived from observations spanning a year’s input of Terra and Aqua data. The spatial resolution is 500 m. Each image has a size of 1200 × 1200 and contains 16 layers. The primary land cover scheme product data were selected for application in this study, which identified 17 land cover classes defined by the International Geosphere Biosphere Programme (IGBP), including 11 classes of vegetation cover, 3 developed and mosaicked land classes, and 3 non-vegetated land classes. We used ArcGIS to splice and transcode Land Cover Type 1 tif data products of MCD12Q1 in 2018 (simulation period). The underlying surface data are used in the WRF-Chem model after being re-defined. In addition, the image remote sensing data of the environmental and disaster monitoring and prediction satellite from the China Resources Satellite Application Center were used to show the distribution of water and land (Figure S1c).

2.2. Model Setting

This study employed the WRF-Chem version 3.9.1 to carry out the simulations, and the simulation period was from 20:00 10 January to 20:00 24 January 2018 (local time, the same as below). The first 24 h were used as the spin-up time, and the subsequent simulation results were used for numerical diagnosis and analysis. In this study, the center of the simulation region was (31.5° N, 115° E), and the vertical direction was divided into 38 layers. The model contained three nested regions (Figure 1). The outermost domain covers central and eastern China, with the spatial resolution being 9 km × 9 km. The middle domain includes the entire Twain-Hu Basin, with the spatial resolution being 3 km × 3 km. The inner nested region includes several small lakes scattered around Wuhan and large-lake areas of the Dongting Lake and Hong Lake, with the spatial resolution being 1 km × 1 km. The lateral boundary was updated every 6 h. The model was used in nudging the four-dimensional assimilation process, with grid analysis nudging above the planetary boundary layer, and the nudging process was applied for air temperature, wind, and air moisture, with a nudging coefficient of 3.0 × 10−4 s−1 in this study. The simulation result was output every hour.
The simulations used the Lin microphysics scheme [43], RRTM long-wave radiation scheme [44], Goddard short-wave radiation scheme [45], MM5 Monin-Obukhov near-stratum scheme, Noah land surface process scheme [46], and ACM2 planetary boundary layer scheme; the Grell 3D cumulus parameterization scheme was only applied to the coarse domain d01 with the horizontal resolution of 9 km (Figure 1a). The model was also coupled with the lake module of the community land model [47]. The results showed that the lake–atmosphere coupled model with the lake module and optimized parameters is superior to the model without the lake scheme [47,48]. The biogenic emissions were calculated through an online Model of Emissions of Gases and Aerosols from Nature (MEGAN, https://www.acom.ucar.edu/wrf-chem/download.shtml (accessed on 7 February 2023). The wintertime biological VOC emissions were much weaker than anthropogenic NMVOC emissions in the THB region over central China (Figure S2), which could indicate that anthropogenic emissions were the main sources of PM2.5 for the wintertime PM2.5 pollution event.
In order to analyze the influence of dense “water networks” on the atmospheric environment in the THB, a set of simulation experiments was performed in this study. The real land use data were used in the control experiment (denoted as lake experiment), and the water bodies in the model domain were changed to farmland underlying surface in the sensitive experiment (denoted as no-lake experiment). By studying the differences between the meteorological elements and air pollutants simulated from the lake experiment and no-lake experiment, the effect of the underlying surface with dense “water networks” in the THB was investigated as follows.

2.3. Methods of Statistical Analysis

Several statistical indicators, including correlation coefficient ( R ), root mean square error ( R M S E ), mean bias ( M B ), general error ( G E ), mean fractional bias ( M F B ), and mean fractional error ( M F E ) were used to evaluate the simulation results. The indicators were defined as follows:
R = i = 1 N ( X i X ¯ ) ( Y i Y ¯ ) i = 1 N ( X i X ¯ ) 2 i = 1 N ( Y i Y ¯ ) 2 ,
R M S E = i = 1 N ( X i Y i ) 2 N ,
M B = i = 1 N ( X i Y i ) N ,
G E = i = 1 N | X i Y i | N ,
M F B = 1 N i = 1 N ( X i Y i ) ( X i / 2 + Y i ) × 100 % ,
M F E = 1 N i = 1 N | X i Y i | ( X i / 2 + Y i ) × 100 % ,
where X i and Y i represent the simulated and observed values, respectively. The indicator of R with a value closer to 1, or M B with a value closer to 0, or R M S E and G E with smaller values indicate a better simulation effect [49].
In order to quantify the effects of meteorological elements on PM2.5 variations induced by the underlying surface of “water networks”, multiple linear regression is adopted to quantify the contribution of meteorological factors to the change in air pollutants by using the software MATLAB (https://ww2.mathworks.cn/, accessed on 7 February 2023).

3. Results and Analysis

3.1. Simulation and Verification of Heavy Pollution Process

From 14 to 22 January in 2018, a local heavy pollution event induced by stagnant weather conditions occurred in the THB (Figure S3). From 14 to 17 January, the local pollutants were gradually accumulated, and the air quality gradually deteriorated in the THB. Air pollution in the THB became worse on 18 January due to the pollutants transported by weak cold air from northern China. From 19 to 22 January, as the THB was controlled by the uniform pressure field after the weak cold air passed through, the wind speed was small, with poor diffusion conditions in the ambient atmosphere. Therefore, the cumulative growth of pollutants led to continuous severe PM2.5 pollution. Figure 2 shows the spatial distributions of observed and simulated PM2.5 concentrations during 14–17 January and the 10 m wind field during 19–22 January. It can be seen that the THB in both stages was affected by weak wind and stagnant weather conditions in the uniform pressure field. The model can well reproduce the continuous accumulation process of air pollutants in the THB under stagnant weather conditions. The simulations are consistent with the observations with respect to the large-value centers of PM2.5 concentrations in North China and the THB. We also estimated the spatial distributions of simulated PM2.5 concentrations and weather conditions against observations over domain 01 from 14 January to 24 January (Figure S3); the simulated wind pressure fields were reasonably close to the observed, and the simulated spatial distribution characteristics of PM2.5 were relatively consistent with the observations.
We used the hourly observations of 2 m air temperature, surface air pressure, 2 m humidity, 10 m wind speed, 10 m wind direction, and precipitation rate at the representative station in the inner domain d03 (Figure 3 and Table 1). Figure 3 shows the simulated and observed near-surface PM2.5 concentrations, surface air pressure, 10 m wind speed, and 10 m wind direction at Wuhan national meteorological observation station from 20:00 11 January to 20:00 24 January 2018. The model well describes the evolution of this pollution process and the variations of the meteorological elements. From 14 to 17 January, the wind speed is maintained at about 2 m·s−1, and the PM2.5 concentration increases slowly. On 18 January, the wind speed increases as the weak cold air moves southward, and the pollution increases due to the pollutant transport. From 19 to 22 January, the wind speed decreases remarkably to about 1 m·s−1, which contributes to the continuous accumulation of pollutants. On 23 January, the wind speed increases, and the pollution quickly dissipates. During the early stage of the pollution process, the simulated PM2.5 concentration is close to the observations. However, in the later stage, the simulations present some deviation from the observations in terms of the spatial distribution due to the underestimated precipitation and higher model resolution than that of the emission inventory. During the pollution process, the PM2.5 concentration is overestimated in simulations, but the changing trend is generally consistent with the observations. Boylan et al. [49] believed that the MFB and MFE are more suitable for evaluating the model performance for prediction of particulate matter and proposed that “reasonable” prediction corresponds to −60% ≤ MFB ≤ 60% and MFE ≤ 75%, and “ideal” prediction corresponds to −30% ≤ MFB ≤ 30% and MFE ≤ 50%. Table 1 shows the evaluation results for each simulated variable at Wuhan and Yueyang stations in the simulation domain d03. We also used the method mentioned in Astruop’s paper [50], with the direction differences normalized to stay within ±180 degrees to calculate the statistical parameters of wind direction (Table S1). All the values of evaluation indicators are within reasonable ranges, indicating that the simulation effect is relatively ideal, which further proves that the model has good applicability over the underlying surface of “water networks”.

3.2. Influence of Dense “Water Network” on PM2.5 Concentrations

The simulation results show that the air pollution in the THB reaches the level of heavy PM2.5 pollution during 19–22 January 2018 (Figure S4). In order to clarify the regulating effect of the underlying surface of “water networks” on heavy air pollution, the influence of such underlying surfaces on the spatiotemporal distribution of PM2.5 concentrations in the THB during the heavy pollution process was first analyzed. The comparison between the lake experiment and no-lake experiment reveals that the underlying surface of “water networks” in the THB generally reduces the regional PM2.5 concentrations (Figure 4), especially over the lake and its surrounding areas. However, the influence from the underlying surface is not uniform in the spatial distribution. In areas surrounding the Dongting Lake and Honghu Lake, the PM2.5 concentration decreases remarkably. The PM2.5 concentration over the lake surface decreases by more than 20 μg/m3, reducing the regional average PM2.5 concentrations by 9.43%. However, for numerous scattered small lakes around Wuhan, the impact from the underlying surface of “water networks” on PM2.5 concentrations is rather complex, which is manifested as the local aggravation or alleviation of PM2.5 pollution, contributing negatively to the regional average PM2.5 concentrations by 0.79%. Compared with the land, water bodies have greater heat capacity and thermal inertia, which most directly affects the vertical water vapor flux [8]. Although the area proportions of water networks for centralized big lakes (CBLs) and discrete small lakes (DSLs) are similar (12.5% and 11.8% of total grid points for the water networks in the two types of lakes, respectively), the simulation results show that the upward fluxes of water vapor over these two types of underlying surface are very different (Figure S2). Different degrees of water-heat exchange cause differences in local meteorological conditions, which further regulates the spatial distribution of PM2.5 concentrations. This may be one reason for the non-uniform distribution of the impacts from the dense water network on the regional atmospheric environment in the THB.
Figure 5 shows the diurnal variations of near-surface PM2.5 concentration differences between the simulations from the lake experiment and no-lake experiment averaged over CBLs and DSLs during the air pollution process. Figure S5 shows the PM2.5 concentrations in the lake experiment simulation. It is revealed that the influences of the “water network” on the PM2.5 exhibit obvious diurnal variations, which vary with the density and area of “water network”. Around sunrise, water bodies cause the greatest decreases in PM2.5 concentrations over both CBLs and DSLs. The CBLs reduced PM2.5 by an average of 16.28 μg/m3 at 08:00 (maximum 29.28 μg/m3), while the DSLs reduced PM2.5 by an average of 13.29 μg/m3 (maximum 24.42 μg/m3). More significantly, the CBLs reduced the PM2.5 concentrations all day, with the stronger potential to “purify” PM2.5 particles at night. The DSLs aggravated local PM2.5 pollution from 10:00 to 00:00, with the peak of 17.91 μg/m3 at 00:00 and decreasing PM2.5 concentrations in other periods. In general, CBLs can result in a decrease of the local PM2.5 concentrations by 4.15–16.28 μg/m3 (4.90–17.68%), while the underlying “water network” of DSLs has weak potential to clean the air environment, with the reverse impacts on PM2.5 levels ranging from −13.29–3.76 μg/m3 in the relative change rates from a negative value of −5.63% to a positive value of 1.56%.
Figure 6 shows the influences of the density and area of “water networks” on the vertical distribution of PM2.5 concentrations below 2 km in the atmosphere. Figure S6 shows the PM2.5 concentrations in the lake experiment simulation. It can be seen that the CBLs improve the local air quality remarkably, which can reduce the PM2.5 concentrations from the ground up to 1500 m at 02:00. Moreover, the PM2.5 concentration within the height of 500 m above the ground can be reduced by about 9.60–13.11%. At 14:00, the PM2.5 concentration near the ground in the surrounding area of the CBLs decreases by 1.71–6.24%, and the influence range reaches up to 1000 m in the vertical direction. In contrast to the CBLs, the influence of DSLs on PM2.5 concentration is more complicated, which is manifested by greater uncertainty and weak effects near the surface. In addition, the DSLs tend to deteriorate the air quality near the surface during the daytime.

3.3. Influence Mechanism of Dense “Water Network” on Atmospheric Environment

The river and lake systems are usually cold in the day and warm at night. That is, the lake acts as a cold (warm) source, which has a good cooling (warming) effect in the daytime (nighttime). Meanwhile, the water vapor evaporation will cause more latent heat release. The inhomogeneity of water and land properties can affect the spatiotemporal distribution of sensible heat flux and latent heat flux [8].
Figure 7 shows the diurnal variations for the differences of simulated, spatially averaged sensible heat flux, latent heat flux, 2 m temperature, and atmospheric boundary layer height between the lake experiment and no-lake experiment for the CBLs and DSLs. Figure S7 shows the variables in the lake experiment simulation. During the day, the sensible heat flux over the underlying surface of “water networks” is reduced, which causes little evaporation and weak variation of latent heat flux over the CBLs. The decreases of sensible and latent heat fluxes over the DSLs during the day are relatively more obvious, which make the surrounding areas of the DSLs a strong cold source. At night, the sensible and latent heat fluxes over the underlying surface of “water networks” increase remarkably, especially for the latent heat release over the CBLs. The common changes of sensible and latent heat fluxes make the underlying surface of “water networks” show a warming effect, and the warm lake effect of the CBLs is more obvious. The diurnal variation of lake effect on the local atmospheric environment is consistent with the results of many previous studies [8,51]. However, different types of dense “water networks” exhibit different effects in the THB. The DSLs show a stronger cold effect during the day, while the CBLs present a stronger warm effect at night. Both are important reasons for the uneven distribution of the influences on the PM2.5 concentrations from the underlying surface of “water networks” in the THB.
The underlying surface of dense “water networks” in the THB obviously regulates the local thermal conditions, causing the changes in the structure of atmospheric boundary layer, which plays a crucial role in the distribution and diffusion of atmospheric particulate matter [8]. Figure 8 shows the spatial differences of daytime and nighttime 2 m temperature, atmospheric boundary layer height, and water vapor flux transport in the horizontal direction between the lake experiment and no-lake experiment. Figure S8 shows the variables in the lake experiment simulation. In the daytime, the cold lake effect of the DSLs is stronger (Figure 8a), which cools the atmospheric environment near the surface by 0.23 °C, and the vertical influence range reaches up to 1 km (Figure S9). As a result, the atmospheric stratification tends to be stable, and the height of the atmospheric boundary layer in the surrounding area of the lakes decreases remarkably (Figure 8c). The diurnal variation of the planetary boundary layer height is consistent with the diurnal variation of air temperature (Figure 7), especially in the DSL area; the height of local atmospheric boundary layer decreases by more than 100 m, which favors the accumulation of air pollutants. At night, the underlying surface of “water networks” shows a warm lake effect (Figure 8b), which causes an unstable and thicker atmospheric boundary layer (Figure 8d), indicating that the boundary layer changes over the water network are dominated by thermal effects, affecting the spatial and temporal variations of local PM2.5 over the “water network”. Especially in the CBL area, the stronger warm lake effect increases the temperature near the surface by 0.36 °C, and the vertical influence range reaches up to 500 m (Figure S9). The uplift of the local boundary layer even exceeds 100 m in local area (Figure 8d), which is favorable for the diffusion and alleviation of local pollution. The cold (warm) lake effect near the surface during the day (night) leads the atmospheric boundary layer to be more stabilized (unstable) and decreases (increases) the boundary layer height, which is conducive to the accumulation (diffusion) of local PM2.5 pollution.
The regulating effect of the underlying surface of “water networks” on local temperature causes the temperature differences between the lake surface and surrounding land area, which triggers local atmospheric circulation. During the day, the cold lake effect promotes the divergence of the local wind field and water vapor flux (Figure 8e). The smooth lake surface strengthens the local wind speed, which facilitates the diffusion and alleviation of local air pollution. The warm lake effect at night promotes the local convergence over the underlying surface of “water networks” (Figure 8f). Hence, the water vapor flux transported from the surrounding area is enhanced. The precursor and water vapor transport in surrounding source regions may aggravate the secondary generation and accumulation of PM2.5 particles.
In summary, the “water network” regulates the local air temperature by influencing the sensible and latent heat fluxes, which therefore changes the energy between the surface and atmosphere and the structure of the atmospheric boundary layer. Moreover, the local atmospheric circulation formed by the lake–land thermal difference drives the transport of air pollutants and water vapor from the surrounding source regions and further influences the spatiotemporal distribution of pollutants by affecting the physical processes in the atmospheric boundary layer and the generation of secondary aerosols.
Table 2 shows the correlations between the difference of PM2.5 concentration and the difference of meteorological elements between the lake experiment and no-lake experiment during the simulation period (312 samples during 11–24 January 2018). It is found that the PM2.5 changes induced by the underlying surface of “water networks” are significantly correlated with the changes in meteorological elements, proving that the atmospheric boundary layer over the underlying surface of dense “water networks” in the THB significantly affects the variation of regional PM2.5 concentration.
However, a problem to be further explored is which factor plays the dominant role under the joint action of these factors. Multiple linear regression is adopted to quantify the contribution of each factor to the changes in the dependent variables. The differences in the meteorological factors between the lake experiment and no-lake experiment are taken as independent variables, including sensible heat flux ( x 1 ), latent heat flux ( x 2 ), 2-m air temperature ( x 3 ), 2-m relative humidity ( x 4 ), atmospheric boundary layer height ( x 5 ), and 10-m wind speed ( x 6 ). With the difference of PM2.5 concentration between the lake experiment and no-lake experiment as the dependent variable ( y ), a set of standardized multiple linear regression equations are established, as follows:
y = 0.2988 x 1 + 0.6379 x 2 0.8500 x 3 + 0.3232 x 4 0.0190 x 5 + 0.0601 x 6 ,
y = 0.2702 x 1 + 0.3668 x 2 0.8317 x 3 + 0.0894 x 4 0.0364 x 5 0.0759 x 6 ,
Equation (7) is the standardized multiple linear regression equation averaged over the CBLs, and Equation (8) is for the DSLs. The correlation coefficients between the fitting results and dependent variables ( y ) in Equations (7) and (8) are 0.77 and 0.76, respectively, and both pass the t-test at the 0.01 significance level. The positive (negative) coefficients in the equation indicate that a certain independent variable exerts the same (opposite) effect on the dependent variable as the joint influence of all variables.
The interaction between thermal factors and water vapor over the underlying surface of “water networks” is an important mechanism affecting the changes of the atmospheric boundary layer and the aerosol physicochemical processes. Equations (7) and (8) reveal that the changes of sensible heat flux ( x 1 ) and 2-m air temperature ( x 3 ), as thermal factors, pose negative effects on the change of PM2.5 concentration between the presence and absence of lakes, that is, the warm (cold) lake effect promotes the decrease (increase) in PM2.5 concentration. The increases in latent heat flux ( x 2 ) and 2-m relative humidity ( x 4 ) over the underlying surface of “water networks” contribute positively to the increase in local PM2.5 concentration. The latent heat flux represents the vertical transport of water vapor in water–heat exchange between the “water network” and atmosphere. The lake acts a long-term water vapor source in the lower atmosphere, and the higher latent heat flux and relative humidity have obvious humidification effects on the atmosphere. The water vapor is an important factor affecting the hygroscopic growth of aerosols, which can promote the liquid-phase heterogeneous reaction and accelerate the gas-particle transformation, thus accelerating the generation of secondary aerosols. Hence, the PM2.5 concentration increases [52,53]. The height of the atmospheric boundary layer has a negative effect on the PM2.5 concentration, where the increase in the height is beneficial to the decrease in PM2.5 concentration. In addition, the change of wind speed over the underlying surface of “water networks” of the CBLs positively affects the change of PM2.5 concentration, which may be due to the transports of pollutants or water vapor by the wind. However, for the underlying surface of “water networks” of the DSLs, the change of wind speed poses a negative impact, showing the removal effect of enhanced wind speed on PM2.5 particles over the smooth lake surface.
The standardized regression coefficients were used to estimate the relative contributions of thermal, water vapor, and dynamic factors to the change in PM2.5 concentration over the underlying surface of “water networks” of the CBLs and DSLs (Table 3). For the CBLs, the relative contribution rates of thermal factors (sensible heat flux and air temperature) and water vapor-related factors (latent heat flux and relative humidity) are 52.48% and 43.91%, respectively. These two factors jointly dominate the change in PM2.5 concentrations, while the relative contribution rate of dynamic factors (atmospheric boundary layer height and wind speed) is only 3.62%.
For the DSLs, the thermal factors dominate the change in PM2.5 concentrations, which contributes 65.96% to the total. The relative contribution rate of water vapor-related factors is 27.31%, indicating that the water vapor-related factors contribute obviously less to the PM2.5 over the DSLs than over the CBLs. The relative contribution rate of dynamic factors is 6.72%. These results reflect the difference of lake effect over different types of lakes.
In summary, the thermal and water vapor-related factors are the key meteorological factors through which the underlying surface of “water networks” affects the atmospheric environment of the THB. The thermal factors mainly affect the structure of the atmospheric boundary layer, and the water vapor-related factors mainly contribute to the humidification effect on the atmosphere. The warming and humidification effects of the “water networks” contribute positively and negatively to the “purification” of air pollution, respectively.
As reported by previous studies, lake-land breeze triggers the distinctive wind patterns between daytime and nighttime, which change both the dispersion and transport process of PM2.5 [20,21]. Although the intensity of lake breeze in our study on heavy PM2.5 pollution is much weaker compared with a typical sea breeze in previous reports, the thermal effects of lake breeze in inland lakes on the local meteorology could significantly affect the change of PM2.5 concentrations in heavy PM2.5 pollution. This may be related with the thermal and dynamic effects as affected by the density and area of water networks under different environment and climate backgrounds.

4. Conclusions

This study conducted a set of experiments on the impact of dense “water networks” on a wintertime PM2.5 pollution event in the THB. By comparing the simulation results of a lake experiment and no-lake experiment, the influence of the underlying surface of dense “water networks” in the THB on the variations of PM2.5 concentrations were investigated. Moreover, the key thermal factors and water vapor-related factors in the changes of PM2.5 concentrations in the atmospheric boundary layer were evaluated. The main conclusions are as follows:
The underlying surface of dense “water networks” in the THB leads to the change in PM2.5 concentrations in association with the density and area of lake groups. The CBLs reduce the local PM2.5 concentrations by 4.90–17.68% during the day and night. The influence of the lake on PM2.5 concentrations reaches up to 1.5 km in the vertical direction, and the PM2.5 concentration can be reduced by 9.60–13.11% from the ground up to 500 m. On the contrary, the ability of the DSLs to alleviate PM2.5 pollution is weaker. The influences of the DSLs on the average PM2.5 concentrations over the lake and surrounding areas range from −5.63% to 1.56%. Moreover, the PM2.5 pollution may even be aggravated during the daytime.
The CBLs show a stronger warm source effect at night, while the DSLs show a stronger cold source effect during the day, which are important reasons for the unevenly-distributed influences on regional PM2.5 concentration by the underlying surface of “water networks” in the THB. The dense “water network” in the THB regulates local air temperature by influencing the changes in sensible and latent heat fluxes. The changes in the land-air energy budget change the structure of atmospheric boundary layer, and the warm (cold) lake effect promotes the decrease (increase) in PM2.5 concentration. Meanwhile, the local atmospheric circulation formed by the lake-land thermal difference also drives the transport of air pollutants and water vapor from the surrounding source areas, which further affects the spatiotemporal distribution of PM2.5 concentration by influencing the physical process in the atmospheric boundary layer and the generation of secondary aerosols.
Thermal factors and water vapor-related factors are the key meteorological factors affecting the atmospheric environment of the THB over the underlying surface of the “water network”. The warming and humidification effects of the “water network” contribute positively and negatively to the “purification” of air pollution, respectively. For the CBLs, the relative contribution rates of thermal factors and water vapor–related factors are 52.48% and 43.91%, respectively. The two factors jointly dominate the change in PM2.5 concentrations. For the DSLs, the thermal factors contribute 65.96% of the total, which dominate the change in PM2.5 concentrations, while the relative contribution rate of water vapor-related factors is 27.31%. This indicates that the contribution of water vapor-related factors over the CBLs is stronger.
In order to fully understand the potential impact of the underlying “water network” surface on air pollution, long-term studies should be carried out, with the fine meteorological and environmental measurements involving air-surface exchanges of energy and water as well as physical and chemical processes of air pollution and more comprehensive simulations of lake breeze and local air pollution under different scenarios.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxics11020169/s1, Figure S1: The distribution of (a) land and water, (b) simulated water vapor flux (QFX), and (c) raster pixel image of CCD2 Level2 on HJ1A Satellite (horizontal resolution of 30 m) in the modeling domain d03. The small box indicates the DSLs, the large box indicates the CBLs; Figure S2: The spatial distributions of (a) biogenic VOC emissions from MEGAN, and (b) NMVOC emissions from MEIC in the modeling domain d01 averaged for the simulation period; Figure S3: Spatial distributions of (left panel) observed and (right panel) simulated near–surface PM2.5 concentrations (color dots (left panel), color contours (right panel), μg/m3), 10-m wind vectors (m/s) and sea level pressure (hPa) at 00:00 (UTC) during 14–24 January, 2018. Blue lines indicate sea level air pressure. Blue area denotes the Pacific Ocean in the left panels; Figure S4: Time variations of simulated near-surface PM2.5 concentrations at Wuhan station during the simulation period. Orange line shows the threshold of heavy pollution level; Figure S5: Diurnal variations of PM2.5 concentrations averaged over CBLs and DSLs during air pollution process from 19–22 January 2018 in the lake experiment of simulation; Figure S6: Vertical distribution of daytime (14:00 in local time) and nighttime (02:00 in local time) PM2.5 concentrations below the height of 2 km in the atmosphere averaged over CBLs and DSLs during air pollution process from 19 to 22 January 2018 in the lake experiment of simulation; Figure S7: Diurnal variations (hourly Time) of (a) sensible heat flux (HFX, W/m2), (b) latent heat flux (LH, W/m2), (c) 2-m air temperature (T2, °C), and (d) planetary boundary layer height (PBLH, m) averaged over CBLs and DSLs during air pollution process over January 19–22, 2018 in the lake experiment of simulation; Figure S8: Spatial distributions of (a,b) 2-m temperature, (c, d) atmospheric boundary layer height, and (e,f) water vapor flux transport in the horizontal direction averaged in the (a,c,e) daytime (14:00 in local time) and (b,d,f) nighttime (02:00 in local time) from the simulations of the lake experiment during the pollution process from 19 to 22 January 2018; Figure S9: Vertical distribution of air temperature differences (Temperature) between the simulations of lake experiment and no-lake experiment averaged for CBLs and DSLs (Height) below 2 km in the atmosphere during the air pollution process from 19–22 January 2018; Table S1: Statistical evaluation of 10-m wind direction at Wuhan and Yueyang station from 20:00 11 January to 20:00 24 January 2018.

Author Contributions

Conceptualization, Y.Z. (Yan Zhu), Y.B., J.X. (Jie Xiong) and T.Z.; Resources, Y.B. and J.X. (Jie Xiong); Data curation, Y.B. and K.M.; Methodology, Y.Z. (Yan Zhu), Y.B., J.X. (Jie Xiong) and T.Z.; Investigation, J.X. (Jiaping Xu), Y.Z. (Yue Zhou), K.M., C.M., X.S. and W.H.; Writing—original draft, Y.Z. (Yan Zhu); Writing—review and editing, Y.B., J.X. (Jie Xiong) and T.Z.; Funding acquisition, Y.B. and T.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant Nos. 41830965, 42075186).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this paper may be provided by Yan Zhu ([email protected]) upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scott, R.W.; Huff, F.A. Impacts of the Great Lakes on Regional Climate Conditions. J. Great Lakes Res. 1996, 22, 845–863. [Google Scholar] [CrossRef]
  2. Long, Z.; Perrie, W.; Gyakum, J.; Caya, D.; Laprise, R. Northern Lake Impacts on Local Seasonal Climate. J. Hydrometeorol. 2007, 8, 881–896. [Google Scholar] [CrossRef]
  3. Thiery, W.; Davin, E.L.; Panitz, H.J.; Demuzere, M.; Lhermitte, S.; van Lipzig, N. The Impact of the African Great Lakes on the Regional Climate. J. Clim. 2015, 28, 4061–4085. [Google Scholar] [CrossRef]
  4. Wen, L.; Lv, S.; Li, Z.; Zhao, L.; Nagabhatla, N. Impacts of the Two Biggest Lakes on Local Temperature and Precipitation in the Yellow River Source Region of the Tibetan Plateau. Adv. Meteorol. 2015, 2015, 248031. [Google Scholar] [CrossRef]
  5. Wentworth, G.R.; Murphy, J.G.; Sills, D.M.L. Impact of Lake Breezes on Ozone and Nitrogen Oxides in the Greater Toronto Area. Atmos. Environ. 2015, 109, 52–60. [Google Scholar] [CrossRef]
  6. Diallo, I.; Giorgi, F.; Stordal, F. Influence of Lake Malawi on Regional Climate from a Double-Nested Regional Climate Model Experiment. Clim. Dyn. 2018, 50, 3397–3411. [Google Scholar] [CrossRef]
  7. Li, Z.; Lyu, S.; Wen, L.; Zhao, L.; Ao, Y.; Wang, S. Effect of a Cold, Dry Air Incursion on Atmospheric Boundary Layer Processes over a High-Altitude Lake in the Tibetan Plateau. Atmos. Res. 2017, 185, 32–43. [Google Scholar] [CrossRef]
  8. Su, D.; Wen, L.; Gao, X.; Leppäranta, M.; Song, X.; Shi, Q.; Kirillin, G. Effects of the Largest Lake of the Tibetan Plateau on the Regional Climate. J. Geophys. Res. Atmos. 2020, 125, e2020JD033396. [Google Scholar] [CrossRef]
  9. Fu, M. Research on the Impacts of Poyang Lake on Typical Weather Process and the Characteristics of near Surface Boundary Layer. Ph.D. Thesis, Nanjing University of Information Science and Technology, Nanjing, China, 2013. [Google Scholar]
  10. Ren, X.; Wang, Y.; Zhang, Z.; Yang, Y.; Hu, C.; Kang, H. Simulation Studies for Lake Taihu Effect on Surrounding Cities Thermal Environment. Acta Meteorol. Sin. 2017, 75, 645–660. [Google Scholar] [CrossRef]
  11. Ding, Y.; Liu, Y. Analysis of Long-Term Variations of Fog and Haze in China in Recent 50 Years and Their Relations with Atmospheric Humidity. Sci. China Earth Sci. 2014, 57, 36–46. [Google Scholar] [CrossRef]
  12. Lai, S.; Zhao, Y.; Ding, A.; Zhang, Y.; Song, T.; Zheng, J.; Ho, K.F.; Lee, S.; Zhong, L. Characterization of PM2.5 and the Major Chemical Components during a 1-Year Campaign in Rural Guangzhou, Southern China. Atmos. Res. 2016, 167, 208–215. [Google Scholar] [CrossRef]
  13. Lu, D.; Mao, W.; Yang, D.; Zhao, J.; Xu, J. Effects of Land Use and Landscape Pattern on PM2.5 in Yangtze River Delta, China. Atmos. Pollut. Res. 2018, 9, 705–713. [Google Scholar] [CrossRef]
  14. Zhu, C.; Zeng, Y. Effects of Urban Lake Wetlands on the Spatial and Temporal Distribution of Air PM10 and PM2.5 in the Spring in Wuhan. Urban For. Urban Green. 2018, 31, 142–156. [Google Scholar] [CrossRef]
  15. Sun, W.; Liu, Z.; Zhang, Y.; Xu, W.; Lv, X.; Liu, Y.; Lyu, H.; Li, X.; Xiao, J.; Ma, F. Study on Land-Use Changes and Their Impacts on Air Pollution in Chengdu. Atmosphere 2019, 11, 42. [Google Scholar] [CrossRef]
  16. Haurwitz, B. Comments on the Sea-Breeze Circulation. J. Atmos. Sci. 1947, 4, 1–8. [Google Scholar] [CrossRef]
  17. Lyons, W.A. The Climatology and Prediction of the Chicago Lake Breeze. J. Appl. Meteorol. Climatol. 1972, 11, 1259–1270. [Google Scholar] [CrossRef]
  18. Physick, W.L. Numerical Experiments on the Inland Penetration of the Sea Breeze. Q. J. R. Met. Soc. 1980, 106, 735–746. [Google Scholar] [CrossRef]
  19. Liu, J.; Zhu, L.; Wang, H.; Yang, Y.; Liu, J.; Qiu, D.; Ma, W.; Zhang, Z.; Liu, J. Dry Deposition of Particulate Matter at an Urban Forest, Wetland and Lake Surface in Beijing. Atmos. Environ. 2016, 125, 178–187. [Google Scholar] [CrossRef]
  20. Keen, C.S.; Lyons, W.A. Lake/Land Breeze Circulations on the Western Shore of Lake Michigan. J. Appl. Meteorol. Climatol. 1978, 17, 1843–1855. [Google Scholar] [CrossRef]
  21. Wang, F.; Wang, Y.; Gao, S.; Hu, N. Simulation Analysis of the Influence of Lake-Land Breeze Circulation on High Ozone Concentration Events. Acta Sci. Circumstantiae 2019, 39, 1392–1401. [Google Scholar]
  22. Zhao, F.; Zhou, Y.; Xu, H.; Zhu, G.; Zhan, X.; Zou, W.; Zhu, M.; Kang, L.; Zhao, X. Oxic Urban Rivers as a Potential Source of Atmospheric Methane. Environ. Pollut. 2022, 297, 118769. [Google Scholar] [CrossRef] [PubMed]
  23. Keeler, J.M.; Kristovich, D.A.R. Observations of Urban Heat Island Influence on Lake-Breeze Frontal Movement. J. Appl. Meteorol. Climatol. 2012, 51, 702–710. [Google Scholar] [CrossRef]
  24. Wu, J.; Li, J.; Peng, J.; Li, W.; Xu, G.; Dong, C. Applying Land Use Regression Model to Estimate Spatial Variation of PM₂.₅ in Beijing, China. Env. Sci. Pollut. Res. Int. 2015, 22, 7045–7061. [Google Scholar] [CrossRef] [PubMed]
  25. Fu, X.; Wang, X.; Hu, Q.; Li, G.; Ding, X.; Zhang, Y.; He, Q.; Liu, T.; Zhang, Z.; Yu, Q.; et al. Changes in Visibility with PM2.5 Composition and Relative Humidity at a Background Site in the Pearl River Delta Region. J. Environ. Sci. 2016, 40, 10–19. [Google Scholar] [CrossRef]
  26. Li, G.; Hao, Y.; Yang, T.; Xiao, W.; Pan, M.; Huo, S.; Lyu, T. Enhancing Bioenergy Production from the Raw and Defatted Microalgal Biomass Using Wastewater as the Cultivation Medium. Bioengineering 2022, 9, 637. [Google Scholar] [CrossRef]
  27. Grell, G.A.; Peckham, S.E.; Schmitz, R.; McKeen, S.A.; Frost, G.; Skamarock, W.C.; Edere, B. Fully coupled “online” chemistry within the WRF model. Atmos. Environ. 2005, 39, 6957–6975. [Google Scholar] [CrossRef]
  28. Sha, T.; Ma, X.; Wang, J.; Tian, R.; Zhao, J.; Cao, F.; Zhang, Y.L. Improvement of inorganic aerosol component in PM2.5 by constraining aqueous-phase formation of sulfate in cloud with satellite retrievals: WRF-Chem simulations. Sci. Total Environ. 2022, 804, 150229. [Google Scholar] [CrossRef]
  29. Gustafson, W.I.; Chapman, E.G.; Ghan, S.J.; Easter, R.C.; Fast, J.D. Impact on modeled cloud characteristics due to simplified treatment of uniform cloud condensation nuclei during NEAQS 2004. Geophys. Res. Lett. 2007, 34, 255–268. [Google Scholar] [CrossRef]
  30. Guo, J.; He, J.; Liu, H.; Miao, Y.; Liu, H.; Zhai, P. Impact of various emission control schemes on air quality using WRF-Chem during APEC China 2014. Atmos. Environ. 2016, 140, 311–319. [Google Scholar] [CrossRef]
  31. Lennartson, E.M.; Wang, J.; Gu, J.; Garcia, L.C.; Ge, C.; Gao, M.; Choi, M.; Saide, P.; Carmichael, G.R.; Kim, J.; et al. Diurnal variation of aerosol optical depth and PM2.5 in South Korea: A synthesis from AERONET, satellite (GOCI), KORUS-AQ observation, and the WRF-Chem model. Atmos. Chem. Phys. 2018, 18, 15125–15144. [Google Scholar] [CrossRef]
  32. Liang, C.; Zhang, X.; Li, Z.; Xiao, Y. Temporal and Spatial Evolution Characteristics of River and Lake Systems in the Urbanization Process of Wuhan City. J. North China Univ. Water Resour. Electr. Power 2019, 40, 61–68. [Google Scholar] [CrossRef]
  33. Shen, L.; Wang, H.; Zhao, T.; Liu, J.; Bai, Y.; Kong, S.; Shu, Z. Characterizing Regional Aerosol Pollution in Central China Based on 19 Years of MODIS Data: Spatiotemporal Variation and Aerosol Type Discrimination. Environ. Pollut. 2020, 263, 114556. [Google Scholar] [CrossRef]
  34. Yu, C.; Zhao, T.; Bai, Y.; Zhang, L.; Kong, S.; Yu, X.; He, J.; Cui, C.; Yang, J.; You, Y.; et al. Heavy Air Pollution with a Unique “Non-Stagnant” Atmospheric Boundary Layer in the Yangtze River Middle Basin Aggravated by Regional Transport of PM2.5 over China. Atmos. Chem. Phys. 2020, 20, 7217–7230. [Google Scholar] [CrossRef]
  35. Zhu, Y.; Zhao, T.; Bai, Y.; Xu, J.; Sun, X.; Hu, W.; Chang, J.; Yang, J.; Zhu, C. Characteristics of Atmospheric Particulate Matter Pollution and the Unique Wind and Underlying Surface Impact in the Twain-Hu Basin in Winter. Environ. Sci. 2021, 42, 4669–4677. [Google Scholar] [CrossRef]
  36. Lin, C.; Chen, F.; Huang, J.C.; Chen, W.C.; Liou, Y.A.; Chen, W.N.; Liu, S.C. Urban Heat Island Effect and Its Impact on Boundary Layer Development and Land–Sea Circulation over Northern Taiwan. Atmos. Environ. 2008, 42, 5635–5649. [Google Scholar] [CrossRef]
  37. Blanken, P.D.; Spence, C.; Hedstrom, N.; Lenters, J.D. Evaporation from Lake Superior: 1. Physical Controls and Processes. J. Great Lakes Res. 2011, 37, 707–716. [Google Scholar] [CrossRef]
  38. Hou, P.; Chen, Y.; Qiao, W.; Cao, G.; Jiang, W.; Li, J. Near-Surface Air Temperature Retrieval from Satellite Images and Influence by Wetlands in Urban Region. Appl. Clim. 2013, 111, 109–118. [Google Scholar] [CrossRef]
  39. Hu, X.; Xue, M. Influence of Synoptic Sea-Breeze Fronts on the Urban Heat Island Intensity in Dallas–Fort Worth, Texas. Mon. Weather. Rev. 2016, 144, 1487–1507. [Google Scholar] [CrossRef]
  40. Zhu, C.; Przybysz, A.; Chen, Y.; Guo, H.; Chen, Y.; Zeng, Y. Effect of Spatial Heterogeneity of Plant Communities on Air PM10 and PM2.5 in an Urban Forest Park in Wuhan, China. Urban For. Urban Green. 2019, 46, 126487. [Google Scholar] [CrossRef]
  41. Li, M.; Zhang, Q.; Kurokawa, J.I.; Woo, J.H.; He, K.; Lu, Z.; Ohara, T.; Song, Y.; Streets, D.G.; Carmichael, G.R.; et al. MIX: A Mosaic Asian Anthropogenic Emission Inventory under the International Collaboration Framework of the MICS-Asia and HTAP. Atmos. Chem. Phys. 2017, 17, 935–963. [Google Scholar] [CrossRef]
  42. Qi, J.; Zheng, B.; Li, M.; Yu, F.; Chen, C.; Liu, F.; Zhou, X.; Yuan, J.; Zhang, Q.; He, K. A High-Resolution Air Pollutants Emission Inventory in 2013 for the Beijing-Tianjin-Hebei Region, China. Atmos. Environ. 2017, 170, 156–168. [Google Scholar] [CrossRef]
  43. Lin, Y.L.; Farley, R.D.; Orville, H.D. Bulk Parameterization of the Snow Field in a Cloud Model. J. Appl. Meteorol. Climatol. 1983, 22, 1065–1092. [Google Scholar] [CrossRef]
  44. Mlawer, E.J.; Taubman, S.J.; Brown, P.D.; Iacono, M.J.; Clough, S.A. Radiative Transfer for Inhomogeneous Atmospheres: RRTM, a Validated Correlated-k Model for the Longwave. J. Geophys. Res. 1997, 102, 16663–16682. [Google Scholar] [CrossRef]
  45. Chou, M.D.; Suarez, M.J.; Ho, C.H.; Yan, M.M.H.; Lee, K.T. Parameterizations for Cloud Overlapping and Shortwave Single-Scattering Properties for Use in General Circulation and Cloud Ensemble Models. J. Clim. 1998, 11, 202–214. [Google Scholar] [CrossRef]
  46. Tewari, M.; Wang, W.; Dudhia, J.; LeMone, M.A.; Mitchell, K.; Ek, M.; Gayno, G.; Wegiel, J.; Cuenca, R. Implementation and Verification of the United NOAH Land Surface Model in the WRF Model. In Proceedings of the 20th Conference on Weather Analysis and Forecasting/16th Conference on Numerical Weather Prediction, Seattle, WA, USA, 10–12 January 2004; p. 15. [Google Scholar]
  47. Gu, H.; Jin, J.; Wu, Y.; Ek, M.B.; Subin, Z.M. Calibration and Validation of Lake Surface Temperature Simulations with the Coupled WRF-Lake Model. Clim. Chang. 2015, 129, 471–483. [Google Scholar] [CrossRef]
  48. Gu, H.; Shen, X.; Jin, J.; Zhao, L.; Xiao, W.; Wang, Y. Evaluation on Simulation of Coupled WRF Lake Model to Lake Surface Temperature in Taihu Lake. Meteorol. Mon. 2014, 40, 166–173. [Google Scholar] [CrossRef]
  49. Boylan, J.W.; Russell, A.G. PM and Light Extinction Model Performance Metrics, Goals, and Criteria for Three-Dimensional Air Quality Models. Atmos. Environ. 2006, 40, 4946–4959. [Google Scholar] [CrossRef]
  50. Astrup, P.; Mikkelsen, T.K. Comparison of NWP Prognosis and Local Monitoring Data from NPPs. Radioprotection 2010, 45, S97–S111. [Google Scholar] [CrossRef]
  51. Huang, A.; Lazhu; Wang, J.; Dai, Y.; Yang, K.; Wei, N.; Wen, L.; Wu, Y.; Zhu, X.; Zhang, X.; et al. Evaluating and Improving the Performance of Three 1-D Lake Models in a Large Deep Lake of the Central Tibetan Plateau. J. Geophys. Res. Atmos. 2019, 124, 3143–3167. [Google Scholar] [CrossRef]
  52. Wang, T.; Huang, X.; Wang, Z.; Liu, Y.; Zhou, D.; Ding, K.; Wang, H.; Qi, X.; Ding, A. Secondary Aerosol Formation and Its Linkage with Synoptic Conditions during Winter Haze Pollution over Eastern China. Sci. Total Environ. 2020, 730, 138888. [Google Scholar] [CrossRef]
  53. Nie, W.; Yan, C.; Huang, D.D.; Wang, Z.; Liu, Y.; Qiao, X.; Guo, Y.; Tian, L.; Zheng, P.; Xu, Z.; et al. Secondary Organic Aerosol Formed by Condensing Anthropogenic Vapours over China’s Megacities. Nat. Geosci. 2022, 15, 255–261. [Google Scholar] [CrossRef]
Figure 1. (a) Three nested domains d01, d02, and d03 (in different boxes) and topographical height above sea level (color contours, unit: m) for WRF-Chem simulation; (b) land use type in the domain d03. The big box in the lower left corner of d03 denotes the centralized big lakes (CBLs) Dongting Lake and Honghu Lake, and the small box in the upper right corner denotes the discrete small lakes (DSLs) around Wuhan.
Figure 1. (a) Three nested domains d01, d02, and d03 (in different boxes) and topographical height above sea level (color contours, unit: m) for WRF-Chem simulation; (b) land use type in the domain d03. The big box in the lower left corner of d03 denotes the centralized big lakes (CBLs) Dongting Lake and Honghu Lake, and the small box in the upper right corner denotes the discrete small lakes (DSLs) around Wuhan.
Toxics 11 00169 g001
Figure 2. Spatial distribution of (left panel) observed and (right panel) simulated near-surface PM2.5 concentrations (color dots (left panel), color contours (right panel), μg/m3), and 10 m wind vectors (m·s−1) averaged (a,b) over 14–17 January and (c,d) over 19–22 January in 2018. Blue lines indicate sea level pressure (hPa). Blue area denotes the Pacific Ocean in the left panels.
Figure 2. Spatial distribution of (left panel) observed and (right panel) simulated near-surface PM2.5 concentrations (color dots (left panel), color contours (right panel), μg/m3), and 10 m wind vectors (m·s−1) averaged (a,b) over 14–17 January and (c,d) over 19–22 January in 2018. Blue lines indicate sea level pressure (hPa). Blue area denotes the Pacific Ocean in the left panels.
Toxics 11 00169 g002
Figure 3. Comparison between simulation and observation of near-surface PM2.5 concentrations, surface air pressure (P), 10-m wind speed (10-m WS), and 10-m wind direction (10-m WD) in Wuhan from 20:00 11 January to 20:00 24 January 2018.
Figure 3. Comparison between simulation and observation of near-surface PM2.5 concentrations, surface air pressure (P), 10-m wind speed (10-m WS), and 10-m wind direction (10-m WD) in Wuhan from 20:00 11 January to 20:00 24 January 2018.
Toxics 11 00169 g003
Figure 4. Spatial distribution of averaged PM2.5 concentration difference (μg/m3) between the lake experiment and no-lake experiment during the pollution process, 19–22 January 2018.
Figure 4. Spatial distribution of averaged PM2.5 concentration difference (μg/m3) between the lake experiment and no-lake experiment during the pollution process, 19–22 January 2018.
Toxics 11 00169 g004
Figure 5. Diurnal variations of the near-surface PM2.5 concentration differences (μg/m3) between the simulations from the lake experiment and no-lake experiment, averaged over CBLs and DSLs during the air pollution process of 19–22 January 2018. The shaded regions are the extreme value intervals over CBLs and DSLs.
Figure 5. Diurnal variations of the near-surface PM2.5 concentration differences (μg/m3) between the simulations from the lake experiment and no-lake experiment, averaged over CBLs and DSLs during the air pollution process of 19–22 January 2018. The shaded regions are the extreme value intervals over CBLs and DSLs.
Toxics 11 00169 g005
Figure 6. Vertical distribution of PM2.5 contribution differences (%) between the lake experiment and no-lake experiment averaged for CBLs and DSLs below 2 km in the atmosphere during the air pollution process of 19–22 January 2018. The yellow and purple lines represent the daytime at 14:00 (local time) and the nighttime at 02:00 (local time).
Figure 6. Vertical distribution of PM2.5 contribution differences (%) between the lake experiment and no-lake experiment averaged for CBLs and DSLs below 2 km in the atmosphere during the air pollution process of 19–22 January 2018. The yellow and purple lines represent the daytime at 14:00 (local time) and the nighttime at 02:00 (local time).
Toxics 11 00169 g006
Figure 7. Diurnal variations in differences of (a) sensible heat flux (HFX, W/m2), (b) latent heat flux (LH, W/m2), (c) 2-m air temperature (T2, °C), and (d) planetary boundary layer height (PBLH, m) between the simulations of the lake experiment and no-lake experiment averaged for CBLs and DSLs during the air pollution process from 19–22 January 2018.
Figure 7. Diurnal variations in differences of (a) sensible heat flux (HFX, W/m2), (b) latent heat flux (LH, W/m2), (c) 2-m air temperature (T2, °C), and (d) planetary boundary layer height (PBLH, m) between the simulations of the lake experiment and no-lake experiment averaged for CBLs and DSLs during the air pollution process from 19–22 January 2018.
Toxics 11 00169 g007
Figure 8. Spatial differences of (a,c,e) daytime (14:00 local time) and (b,d,f) nighttime (02:00 local time) (a,b) 2 m temperature, (c,d) atmospheric boundary layer height, and (e,f) water vapor flux transport in the horizontal direction between the lake experiment and no-lake experiment during the air pollution process from 19–22 January 2018.
Figure 8. Spatial differences of (a,c,e) daytime (14:00 local time) and (b,d,f) nighttime (02:00 local time) (a,b) 2 m temperature, (c,d) atmospheric boundary layer height, and (e,f) water vapor flux transport in the horizontal direction between the lake experiment and no-lake experiment during the air pollution process from 19–22 January 2018.
Toxics 11 00169 g008aToxics 11 00169 g008b
Table 1. Statistical evaluation of simulation results at Wuhan and Yueyang stations from 20:00 11 January to 20:00 24 January 2018.
Table 1. Statistical evaluation of simulation results at Wuhan and Yueyang stations from 20:00 11 January to 20:00 24 January 2018.
StationVariablesRRMSEMBGEMFB (%)MFE (%)
Wuhan2-m temperature (°C)0.78 **3.883.113.2547.3349.42
2-m relative humidity (%)0.71 **26.43−23.3423.37−20.9420.97
10-m wind speed (m/s)0.57 **1.280.561.0344.8359.11
10-m wind direction (°)0.27**139.65−5.6883.889.2039.95
Surface air pressure (hPa)0.99 **0.88−0.460.71−0.030.05
Precipitation rate (mm/h)0.60 **0.140.000.04//
Near-surface PM25 (μg/m3)0.56 **56.3428.6838.9916.4122.03
Yueyang2-m temperature (°C)0.80 **3.583.323.3233.2633.29
2-m relative humidity (%)0.82 **16.44−12.6914.18−10.1412.14
10-m wind speed (m/s)0.64 **2.411.882.0550.4555.39
10-m wind direction (°)0.41**117.98−8.3671.319.3442.92
Surface air pressure (hPa)0.99 **1.491.301.320.090.09
Precipitation rate (mm/h)0.70 **0.21−0.020.09//
Near-surface PM25 (μg/m3)0.31 **42.701.8635.635.1526.78
Note: ** indicates that the values have passed the significance test at the 0.01 significance level.
Table 2. Correlation coefficients of PM2.5 concentration difference with influencing factor difference between lake experiment and no-lake experiment.
Table 2. Correlation coefficients of PM2.5 concentration difference with influencing factor difference between lake experiment and no-lake experiment.
VariablesCBLDSL
Sensible heat flux−0.35 **−0.47 **
Latent heat flux−0.37 **−0.57 **
2-m temperature−0.64 **−0.75 **
2-m relative humidity0.54 **−0.08
Atmospheric boundary layer height−0.40 **−0.52 **
10-m wind speed0.01−0.40 **
Note: ** indicates that the values have passed the significance test at 0.01 significance level.
Table 3. Contributions of different factors affected by the underlying surface of “water networks” to the change of PM2.5 concentrations.
Table 3. Contributions of different factors affected by the underlying surface of “water networks” to the change of PM2.5 concentrations.
VariablesCBLsDSLs
Sensible heat flux13.65%16.17%
Latent heat flux29.14%21.96%
2-m temperature38.83%49.79%
2-m relative humidity14.77%5.35%
Atmospheric boundary layer height0.87%2.18%
10-m wind speed2.75%4.54%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhu, Y.; Bai, Y.; Xiong, J.; Zhao, T.; Xu, J.; Zhou, Y.; Meng, K.; Meng, C.; Sun, X.; Hu, W. Mitigation Effect of Dense “Water Network” on Heavy PM2.5 Pollution: A Case Model of the Twain-Hu Basin, Central China. Toxics 2023, 11, 169. https://doi.org/10.3390/toxics11020169

AMA Style

Zhu Y, Bai Y, Xiong J, Zhao T, Xu J, Zhou Y, Meng K, Meng C, Sun X, Hu W. Mitigation Effect of Dense “Water Network” on Heavy PM2.5 Pollution: A Case Model of the Twain-Hu Basin, Central China. Toxics. 2023; 11(2):169. https://doi.org/10.3390/toxics11020169

Chicago/Turabian Style

Zhu, Yan, Yongqing Bai, Jie Xiong, Tianliang Zhao, Jiaping Xu, Yue Zhou, Kai Meng, Chengzhen Meng, Xiaoyun Sun, and Weiyang Hu. 2023. "Mitigation Effect of Dense “Water Network” on Heavy PM2.5 Pollution: A Case Model of the Twain-Hu Basin, Central China" Toxics 11, no. 2: 169. https://doi.org/10.3390/toxics11020169

APA Style

Zhu, Y., Bai, Y., Xiong, J., Zhao, T., Xu, J., Zhou, Y., Meng, K., Meng, C., Sun, X., & Hu, W. (2023). Mitigation Effect of Dense “Water Network” on Heavy PM2.5 Pollution: A Case Model of the Twain-Hu Basin, Central China. Toxics, 11(2), 169. https://doi.org/10.3390/toxics11020169

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

Article Metrics

Back to TopTop