Next Article in Journal
Evaluating the Performance of Pulsed and Continuous-Wave Lidar Wind Profilers with a Controlled Motion Experiment
Previous Article in Journal
An Efficient Task Implementation Modeling Framework with Multi-Stage Feature Selection and AutoML: A Case Study in Forest Fire Risk Prediction
Previous Article in Special Issue
Geometric Factor Correction Algorithm Based on Temperature and Humidity Profile Lidar
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Distribution, Sources, and Impact on Atmospheric Oxidation of Reactive Nitrogen Oxides in the North China Plain Agricultural Regions in Summer

1
Department of Precision Machinery and Precision Instrumentation, University of Science and Technology of China, Hefei 230026, China
2
Key Laboratory of Ecosystem Carbon Source and Sink, China Meteorological Administration (ECSS-CMA), Wuxi University, Wuxi 214105, China
3
Key Lab of Environmental Optics and Technology, Anhui Institute of Optics and Fine Mechanics, Hefei Institutes of Physical Science, Chinese Academy of Sciences, Hefei 230031, China
4
School of Environmental Science and Optoelectronic Technology, University of Science and Technology of China, Hefei 230026, China
5
College of Oceanography, Hohai University, Nanjing 210024, China
6
Center for Excellence in Regional Atmospheric Environment, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China
7
Key Laboratory of Precision Scientific Instrumentation of Anhui Higher Education Institutes, University of Science and Technology of China, Hefei 230026, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(17), 3192; https://doi.org/10.3390/rs16173192
Submission received: 6 July 2024 / Revised: 27 August 2024 / Accepted: 28 August 2024 / Published: 29 August 2024

Abstract

:
The lack of vertical observation of reactive nitrogen oxides in agricultural areas has posed a significant challenge in fully understanding their sources and impacts on atmospheric oxidation. Ground-based multi-axis differential optical absorption spectroscopy (MAX-DOAS) observations were conducted in the agricultural regions of the North China Plain (NCP) during the summer of 2019 to measure the vertical distributions of aerosols, nitrogen dioxide (NO2), and nitrous acid (HONO). This study aimed at revealing the spatiotemporal distribution, sources, and environmental effects of reactive nitrogen oxides in the NCP agricultural areas. Our findings indicated that the vertical profiles of aerosols and NO2 exhibited a near-Gaussian distribution, with distinct peak times occurring between 8:00–10:00 and 16:00–18:00. HONO reached its maximum concentration near the surface around 8:00 in the morning and decreased exponentially with altitude. After sunrise, the concentration of HONO rapidly decreased due to photolysis. Additionally, the potential source contribution function (PSCF) was used to evaluate the potential sources of air pollutants. The results indicated that the main potential pollution sources of aerosols were located in the southern part of the Hebei, Shanxi, Shandong, and Jiangsu provinces, while the potential pollution sources of NO2 were concentrated in the Beijing–Tianjin–Hebei region. At altitudes exceeding 500 m, the heterogeneous reactions of NO2 on aerosol surfaces were identified as one of the important contributors to the formation of HONO. Furthermore, we discussed the production rate of hydroxyl radicals (OH) from HONO photolysis. It was found that the production rate of OH from HONO photolysis decreased with altitude, with peaks occurring in the morning and late afternoon. This pattern was consistent with the variations in HONO concentration, indicating that HONO was the main contributor to OH production in the agricultural regions of the NCP. This study provides a new perspective on the sources of active nitrogen in agricultural regions and their contribution to atmospheric oxidation capacity from a vertical perspective.

1. Introduction

Nitrous acid (HONO) plays a critical role in the atmospheric nitrogen cycle by producing hydroxyl radical (OH) through photolysis, thereby influencing atmospheric oxidation capacity [1,2,3,4]. Up to 20–90% of OH formation has been attributed to HONO photolysis [5,6,7,8]. Primary pollutants in the atmosphere, such as carbon monoxide (CO), nitrogen oxides (NOx), and volatile organic compounds (VOCs), can be oxidized by OH, promoting the formation of ozone (O3) and secondary aerosols [9,10]. High concentrations of HONO may seriously threaten people’s health because long-term exposure may harm lung function [11,12]. Additionally, the reaction between HONO and ammonia (NH3) can produce carcinogenic nitrosamines [13]. Despite extensive research on HONO over the years, there is still no consensus on its sources. A substantial amount of research indicates that the sources of HONO can be classified into four categories [14]: (1) homogeneous reactions between nitric oxide (NO) and OH (R1) [15,16,17], (2) direct emissions from vehicle exhaust, soil, and biomass combustion [18,19,20], (3) heterogeneous reactions of nitrogen dioxide (NO2) on wet surfaces (R2) [16,21], and (4) photolysis of nitrate particles in the soil and atmosphere (R3) [15,22]. Reaction (R1) almost exclusively occurs during daytime, under conditions of high NO and OH concentrations [23]. Reaction (R2), particularly on aerosol surfaces, is an important source of HONO and the main source of nighttime HONO [5]. Furthermore, the surface area density and chemical nature of aerosols are the dominant factors controlling the conversion efficiency of NO2 to HONO [24]. Reaction (R3) is closely associated with agricultural activities.
N O + O H + M H O N O + M
2 N O 2 + H 2 O + s u r f a c e H O N O + H N O 3
H N O 3 / N O 3 + h v H O N O / N O 2 + O
In agricultural areas, the extensive use of nitrogen fertilizers has become a significant source of soil emissions of HONO, impacting air quality [25]. Therefore, it is essential to understand the sources, spatiotemporal distribution, and environmental effects of HONO in agricultural regions. In the agricultural regions of Wangdu in the NCP, it has been found that soil fertilization is an important source of HONO [26]. Observations conducted in the agricultural areas of the Huai River Basin revealed that the heterogeneous reactions involving surface NO2 serve as the primary contributor to the formation of nighttime HONO [25]. Moreover, after fertilization, farms released large amounts of HONO, leading to increased HONO concentrations and a higher HONO-to-NO2 ratio [27]. However, these studies mentioned above primarily focused on near-surface observations, lacking insight into the sources and environmental effects of HONO at different altitudes in agricultural regions. Research in urban areas has shown that during haze events, high-altitude HONO mainly originates from heterogeneous reactions of NO2 on aerosol surfaces rather than near-surface diffusion [28]. Therefore, by detecting the vertical distribution of HONO, we can obtain a better understanding of its sources and contribution to atmospheric oxidation capacity [29].
Currently, common vertical observation methods for HONO include tower observations [28,30,31], airborne and balloon-borne observations [32,33,34], long-path differential optical absorption spectroscopy (LP-DOAS) observations [35,36], and multi-axis differential optical absorption spectroscopy (MAX-DOAS) remote sensing observations [1,5,37]. However, tower observations can only measure the concentrations of atmospheric pollutants up to 400 m [38]. Observations conducted via airborne and balloon platforms face constraints related to flight time, posing difficulties for prolonged continuous monitoring [39]. Although the LP-DOAS technique facilitates non-contact and continuous measurements, its operational altitude remains restricted. As a promising spectral remote sensing technology, MAX-DOAS can obtain real-time and continuous vertical distributions of atmospheric pollutants [40,41,42,43]. In recent years, MAX-DOAS has been used worldwide for HONO observations, including in Madrid, Spain [1], Melbourne, Australia [37], Beijing, China [5], and Leshan, China [44].
Even though MAX-DOAS has been used to observe HONO in several locations in China, most observation sites are in urban regions, and few studies have focused on agricultural sources. Therefore, in this study, MAX-DOAS observation experiments were carried out in the agricultural regions of the NCP to investigate the spatiotemporal distributions and potential sources of aerosols, NO2, and HONO. In addition, considering the critical role of aerosols in the heterogeneous reaction of NO2 to HONO, we analyzed the relationship between aerosols and the conversion rate from NO2 to HONO. Furthermore, based on vertical observations of HONO, this study opens a new window on rethinking the sources of reactive nitrogen oxides in agricultural regions and their contribution to atmospheric oxidation capacity from a vertical perspective.
The structure of this work is outlined as follows: Section 2 details the MAX-DOAS observation campaign, the retrieval algorithm for pollutants, the Potential Source Contribution Function (PSCF), and additional supporting data. In Section 3, the potential sources and spatiotemporal distribution of pollutants are presented. The formation of HONO through the heterogeneous reactions of NO2, and OH radicals from HONO photolysis are evaluated. In Section 4, further analyses of the results are conducted, along with comparisons to results from other regions.

2. Materials and Methods

2.1. MAX-DOAS Measurement Setup

The MAX-DOAS instrument was set up at the Rural Environmental Research Station of the Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences (115.26°E, 38.68°N), located in Wangdu County, Baoding City, in the NCP, from 1 June 2019 to 31 August 2019 (Figure 1). The Wangdu station is surrounded by farmland where winter wheat and summer maize are grown in rotation, with almost no industrial emission sources nearby [27,45]. This setup aimed to study the vertical distributions, potential sources, and contribution to atmospheric oxidation of reactive nitrogen oxides in agricultural environments. The station is situated at an altitude of 35 m, and the observational azimuth angle was set to 78°. The instrument consists of a telescope system and a spectral analysis system. In the telescope system, a quartz right-angle prism, rotated by a stepper motor, receives scattered sunlight at different elevation angles. After passing through the prism, the sunlight gathers through a plano-convex lens to the optical fiber entrance, and eventually enters the spectral analysis system. The field of view (FOV) of the telescope system is less than 0.5°, and the elevation accuracy is better than 0.1°. The spectral analysis system primarily consists of two Avantes spectrometers, which split the collected scattered sunlight and convert the light signal into an electrical signal using a Charge-Coupled Device (CCD), thus obtaining spectral data from the scattered light. Because the CCD is highly sensitive to temperature changes, its temperature is maintained at 20 ± 0.5 °C by a temperature sensor and a semiconductor cooling device near the CCD. The effective spectral range collected by the instrument is 300–500 nm, with a spectral resolution of 0.49 nm. Further detailed descriptions of the instrument can be found in previous studies [46,47].
During the observation period, each observation series consisted of 11 elevation angles: 1, 2, 3, 4, 5, 6, 8, 10, 15, 30, and 90°. The integration time for each elevation angle was adjusted according to the solar zenith angle (SZA), with the single exposure time and scan number automatically adjusted based on the light intensity. It takes about 12 min to complete an observation series.

2.2. Spectral Analysis

The differential slant column densities (DSCDs) of aerosols, NO2, and HONO were retrieved using the QDOAS software (version 3.4) developed by the Belgian Institute for Space Aeronomy (BIRA-IASB) (http://uv-vis.aeronomie.be/software/QDOAS/, accessed on 22 March 2024). Before spectral analysis, dark current and electronic offset were subtracted from the spectra, and wavelength calibration was performed with a standard mercury lamp [48,49]. Additionally, cloud effects are filtered out based on color index [50,51,52,53], as detailed in Section S1. Table 1 provides the detailed retrieval settings for O4, H2O, NO2, and HONO. The details about the Ring effect and I0 effect can be found in Section S2. An example of DOAS fitting results for a measured spectrum with an elevation of 4° and a solar zenith angle of 72.32° on 6 June 2019, is shown in Figure 2. To ensure the reliability of the DOAS retrieval, NO2 and HONO DSCDs were excluded when the root mean square (RMS) of fitting residuals exceeded 2 × 10−3. As a result, the corresponding qualified DSCD results of NO2 and HONO were 85.84% and 90.81%, respectively.

2.3. Vertical Profile Retrievals

In this study, the vertical profiles of aerosols, NO2, and HONO were retrieved based on the Heidelberg Profile retrieval algorithm (HEIPRO) [63,64]. This algorithm relies on the Optimal Estimation Method (OEM) and utilizes the radiative transfer model SCIATRAN (version 2.2) as the forward model [65]. The so-called maximum a posteriori (MAP) solution is obtained from the following cost function:
χ 2 x = y F x , b ] T S ε 1 y F x , b + x x a ] T S a 1 x x a
Here, x is the aerosols or trace gas profiles; y is the DSCDs of different components at different elevation angles; F x , b is the forward model simulation value of DSCDs; b is meteorological parameters (e.g., pressure and temperature vertical profiles); x a is prior information; S a is the covariance matrix of x ; and S ε is the covariance matrix of y . In the retrieved process, the optical parameters of the aerosols were selected as follows: the asymmetry parameter was 0.7, the surface albedo was 0.05, and the single scattering albedo was 0.95 [66]. The atmosphere was divided into 20 layers from 0 to 2 km, each with a thickness of 100 m. The time resolution of the vertical profiles was set to 15 min to ensure that a complete observation elevation series could be covered. The aerosol optical depth (AOD), which measures the extent to which aerosols prevent the transmission of sunlight by absorption or scattering, and the vertical column densities (VCD), which represent the total number of molecules of trace gases in a column of air extending from the ground to the top of the atmosphere, were obtained by integrating the concentrations at different heights. The detailed retrievals process can be found in Figure S1.

2.4. Calculation of Pollutant Growth Percentages during Haze Days Compared to Non-Haze Days

According to the China Meteorological Administration’s haze standard (QX/T 113-2010), weather conditions are classified into two categories based on visibility and relative humidity: non-haze days (visibility > 10 km) and haze days (visibility < 10 km and RH < 80%) [67,68,69]. Haze days accounted for 11% of the data, while non-haze days accounted for 89%.
Compared to non-haze days, the growth percentages of atmospheric pollutants on haze days can be expressed as follows:
G = P n o n h a z e P h a z e P n o n h a z e × 100 %
Here, G represents the growth percentages of atmospheric pollutants. The average concentrations of atmospheric pollutants at different altitudes on non-haze and haze days are, respectively, represented by P n o n h a z e and P h a z e .

2.5. PSCF Analysis

The Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model with the Global Data Assimilation System (GDAS) was developed by the National Oceanic and Atmospheric Administration Air Resource Laboratory to calculate the backward trajectories of air masses.
The Potential Source Contribution Function (PSCF) is an effective method for identifying pollution source distributions and transport paths based on pollutant concentrations and backward trajectories [70,71]. First, the study area is divided into several grids, and N i j is defined as the total number of trajectory endpoints within the i j th grid cell. Using the average aerosol extinction and NO2 concentrations during the observation period as the threshold values, M i j is defined as the total number of endpoints with a pollutant concentration higher than the threshold within the i j th grid cell. The P S C F within the i j th grid cell is defined as follows:
P S C F i j = M i j N i j
Because PSCF is a conditional probability, considerable uncertainty occurs when the value of N i j is small. In order to reduce the influence of uncertainty, W i j was defined as the weight factor of the i j th grid cell, determined by N i j . Therefore, the P S C F i j is multiplied by the W i j , resulting in the weighted PSCF (WPSCF). The detailed settings of W i j are as follows:
W i j = 1.0                                             N i j > 3 N a v e     0.7                 3 N a v e > N i j > 3 2 N a v e 0.4                 3 2 N a v e > N i j > N a v e 0.2                                                 N a v e > N i j
where N a v e represents the average number of endpoints per grid. The W P S C F within the i j th grid cell is defined as follows:
W P S C F i j = W i j × P S C F i j
In this study, the backward trajectories and WPSCF of the observation site were analyzed using TrajStat software (version 1.5.5) [72,73,74].

2.6. Calculation of Conversion Rate from NO2 to HONO

The data showing an increase in the ratio of HONO to NO2 over a continuous three-hour period were used to calculate the conversion rate from NO2 to HONO [62,66]:
C H O N O = H O N O t 2 N O 2 t 2 H O N O t 1 N O 2 t 1 t 2 t 1
Here, C H O N O represents the conversion rate from NO2 to HONO; H O N O and N O 2 represent the concentration of NO2 and HONO, respectively; and t 2 and t 1 represent the initial and final times, respectively.

2.7. TUV Model and OH Generation Rate Calculation

The Tropospheric Ultraviolet–Visible Model (TUV) was used to calculate the photolysis frequency of HONO [75]. The total ozone column input into the model was obtained from TROPOMI observation, the AOD at 360.8 nm was obtained from the MAX-DOAS instrument, and the single scattering albedo (SSA) was 0.95 [76]. The TUV model has been widely used in observations based on the MAX-DOAS instrument [5,37,77].
The contribution of HONO to the OH generation rate ( P O H H O N O ) can be calculated using the following formula [44,66]:
P O H H O N O = J H O N O × H O N O
Here, H O N O is the HONO concentration obtained by the MAX-DOAS instrument. J H O N O represents the photolysis frequency of HONO.

2.8. Ancillary Data

The Tropospheric Monitoring Instrument (TROPOMI) is aboard the Sentinel-5 Precursor satellite, with a transit at 13:30 local time. TROPOMI has a UV–VIS spectral range of 270–495 nm and a spatial resolution of 5.5 km × 3.5 km [78]. The retrieval algorithm for determining the horizontal distribution of NO2 was developed by the University of Science and Technology of China [79]. The Moderate Resolution Imaging Spectroradiometer (MODIS) is aboard the Earth Observing System (EOS) Terra and Aqua satellites, with transit times at 10:30 and 13:30 local time. The average MODIS AOD and TROPOMI NO2 VCD within a 25 km*25 km area centered on the observation station were selected as verification data for MAX-DOAS AOD and NO2 VCD. Visibility and humidity information was obtained from Shijiazhuang Zhengding International Airport (https://www.wunderground.com/weather/ZBSJ, accessed on 21 May 2024), located approximately 65 km from the observation station. Temperature and boundary layer height information was obtained from the simulation using the Weather Research and Forecasting model (WRF) [80]. Details regarding the WRF model configurations are provided in Section S3 and Table S1.

3. Results

3.1. Validations of NO2 VCD and AOD Measured by MAX-DOAS

The average AOD at 360.8 nm measured by MAX-DOAS between 10:30 and 13:30 local time was selected for correlation verification with AOD measured by MODIS. The AOD at 360.8 nm from MODIS is derived by calculating the Ångström exponent using the AOD measurements at 470 nm and 550 nm, and then applying this exponent to estimate the AOD at 360.8 nm [81,82]. The results indicated a strong correlation between MAX-DOAS and MODIS AOD, with a correlation coefficient of 0.87 (Figure 3a). Similarly, the average NO2 VCD measured by MAX-DOAS between 13:20 and 13:40 local time was selected for correlation verification with NO2 VCD measured by TROPOMI. The results showed a strong correlation between MAX-DOAS and TROPOMI NO2 VCD, with a correlation coefficient of 0.77 (Figure 3b). The comparison results demonstrated that the data observed by MAX-DOAS were reliable.

3.2. Measurement Overview

In this study, in order to discuss the vertical distribution characteristics of aerosols, NO2, and HONO during the summer, four typical altitudes ranges were selected: (1) 0–100 m, representing the bottom boundary layer; (2) 100–300 m, representing the lower boundary layer; (3) 500–700 m, representing the middle boundary layer; and (4) 900–1100 m, representing the upper boundary layer. Time series of trace gases (NO2 and HONO) and aerosol extinction in the bottom, lower, middle, and upper layers are shown in Figure 4. The average and maximum values (aerosol extinction, NO2 and HCHO concentration) are presented in Table 2. During the observation period, the time series of meteorological information (including boundary layer height, temperature, relative humidity, and visibility) in the Wangdu area is shown in Figure S2. The scatter plots between meteorological information (RH and temperature) and pollutants (aerosols, NO2, HONO, and NO2/HONO) are shown in Figure S3.
The results showed that NO2 was mainly distributed in the bottom layer (1.659 ppb) and lower layer (1.926 ppb), with the concentration in the lower layer slightly higher than that in the bottom layer. Compared to the bottom layer, the ratios of NO2 concentration in the lower, middle, and upper layers were 116.12%, 9.07%, and 2.44%, respectively. Aerosol extinction in the bottom layer (0.358 km−1), lower layer (0.370 km−1), and middle layer (0.342 km−1) was approximately equal and higher than that in the upper layer (0.161 km−1). Compared to the bottom layer, the ratios of aerosol extinction in the lower, middle, and upper layers were 103.47%, 95.67%, and 44.97%, respectively. Unlike aerosols and NO2, the HONO concentration varied significantly in the four layers, with the highest concentration in the bottom layer (0.079 ppb) and the lowest concentration in the upper layer (0.012 ppb). Compared to the bottom layer, the ratios of HONO concentration in the lower, middle, and upper layers were 53.39%, 26.25%, and 15.20%, respectively. Similar HONO vertical gradients were observed in Colorado using a mass spectrometer on the tower [31], in Los Angeles and Houston using LP-DOAS [35,36], and in Melbourne using MAX-DOAS [37].
Figure 5 illustrates the diurnal variations in the vertical distributions of aerosol extinction, NO2, and HONO observed during the study period, along with the corresponding boundary layer height. In the Wangdu area, the boundary layer height reached its peak at noon and was at its lowest during the morning and evening, with an average boundary layer height of 1064 m during the day. This boundary layer height plays a significant role in shaping the diurnal variations of pollutant concentration [83]. The average aerosol extinction near the surface varied in the order of July (0.377 km−1) > August (0.346 km−1) > June (0.343 km−1), slightly higher than the level observed at the nearby rural site of Raoyang (aerosol extinction: 0.33 km−1 during the summer of 2014) [84]. The average NO2 concentration near the surface varied in the order of June (1.958 ppb) > July (1.341 ppb) > August (1.317 ppb). The average HONO concentration near the surface varied in the order of June (0.094 ppb) > July (0.075 ppb) > August (0.073 ppb).
The vertical profiles of NO2 and aerosol extinction exhibited a near-Gaussian pattern, characterized by an initial increase in concentration and followed by a decrease with altitude during the summer months. The maximum aerosol extinction was observed around 0.4 km, and the aerosol extinction decreased with the increase in altitude beyond 0.4 km. The maximum NO2 concentration was observed around 0.15 km, and the NO2 concentration decreased with the increase in altitude beyond 0.4 km. These trends could be related to the transportation of pollutants. Unlike aerosols and NO2, the maximum HONO concentration consistently occurred near the surface during different summer months. In addition, HONO concentration decreased exponentially with altitude. Above an altitude of 1 km, the HONO concentration was less than 10% of that near the surface. Similar vertical distributions of HONO have been found in previous studies [1,5,31,44,77].
Similar aerosol diurnal variation characteristics were observed during different months. The maximum aerosol extinction occurred near the surface at 8:00. Subsequently, the altitude of maximum aerosol extinction rose with the increasing boundary layer height, peaking at an altitude of 0.5 km at around 10:00. Between 10:30 and 16:00, the distribution of aerosol concentration at various altitudes exhibited a stable pattern, with the maximum altitude consistently observed at 0.5 km. After 16:00, aerosol extinction decreased with the decrease in the atmospheric boundary layer, concentrating at lower altitudes. By around 18:00, the maximum aerosol extinction for the day occurred at an altitude of approximately 0.1 km.
In June, July, and August, the maximum concentrations of NO2 were found at 8:00, 9:00, and 10:00, respectively, showing significant variation characteristics. The peak of NO2 concentrations appeared at 8:00–10:00 and 16:00–18:00. The minimum NO2 concentration appeared at noon, which might be related to the photolysis of NO2 [85]. In addition, the rise of the boundary layer at noon also facilitated the vertical diffusion of NO2 [39]. After 16:00, the reduction in boundary layer height led to an accumulation of NO2 at ground level, resulting in the formation of a secondary peak. Similar diurnal variations have been previously observed in the Wangdu area [86].
The highest concentration of HONO was observed near the surface at 8:00 across all months. This phenomenon can be attributed to the decline in solar radiation after sunset, which reduced the photolytic loss of HONO. Additionally, heterogeneous reactions of NO2 led to the accumulation of HONO during the night, reaching its maximum concentration in the early morning [25,31,44,77]. After sunrise, HONO concentration decreased rapidly due to photolysis, remaining a low concentration from 10:00 to 14:00. After 16:00, HONO concentration increased with the rise in the precursor NO2 concentration and the decrease in sunlight, reaching a second peak near the surface. Similar HONO diurnal variation characteristics were also found in farmland environment studies Wangdu [27] and Xingtai [51] in the North China Plain, China, as well as in Cabauw, the Netherlands [87], and Grignon, France [88].

3.3. The Difference of Aerosol, NO2, and HONO Profiles during Non-Haze and Haze Days

To assess the impact of haze on pollutants, the profiles of aerosol, NO2, and HONO during non-haze and haze days are shown in Figure 6, with the corresponding error bars illustrated in Figure S5.
As shown in Figure 6, the average profile shapes of aerosols, NO2, and HONO were generally similar on both non-haze and haze days. Aerosols and NO2 exhibited a near-Gaussian distribution, while HONO showed an exponential decrease with altitude. During non-haze days, within the 0–0.4 km altitude range, the average aerosol extinction, NO2 concentration, and HONO concentration were 0.345 km−1, 1.457 ppb, and 0.043 ppb, respectively. During haze days, within the same altitude range, the average aerosol extinction, NO2 concentration, and HONO concentration were 0.519 km−1, 1.554 ppb, and 0.046 ppb, respectively. The aerosol extinction on haze days was significantly higher than that on non-haze days, likely due to the increased conversion of directly emitted pollutants into secondary pollutants [9,24,89]. The NO2 and HONO concentrations were slightly higher on haze days compared to non-haze days. On haze days, the higher aerosol concentration enhanced heterogeneous reactions on aerosol surfaces, leading to an increase in HONO concentration [90]. Additionally, the weaker sunlight on haze days prevented the photolysis of NO2, resulting in higher NO2 concentrations [91].
To further quantify the differences in the vertical distribution of pollutants between non-haze and haze days at different altitudes, the growth percentages of atmospheric pollutants were calculated. As shown in Figure 6, the highest growth percentage of aerosol occurred at 0.4 km, approximately 62%, and decreased with altitude above 0.4 km, reaching approximately zero above 0.9 km. This indicated that aerosol pollution on haze days mainly occurred at mid-altitude and low-altitude layers. The highest growth percentages of NO2 were observed at altitudes of 0.1 km and 0.2 km, approximately 21% and 26%, respectively. A negative growth percentage of approximately −40% was observed at 0.4 km. The highest growth percentage of HONO occurred at an altitude of 0.1 km, approximately 13%, and decreased slowly with altitude above 0.1 km. Compared to aerosols and NO2, the growth percentage of HONO is much smaller, indicating that haze exerts a lesser influence on HONO levels.

3.4. Aerosols and NO2 Potential Sources in Agricultural Regions in Summer

The observation station is located in farmland with no apparent emission sources nearby. To further elucidate the potential sources of aerosols and NO2 at the station and to assess transport influences, we calculated the backward trajectories of the air mass at the bottom boundary layer (50 m), lower boundary layer (200 m), middle boundary layer (600 m), and upper boundary layer (1000 m) arriving at Wangdu station in the summer of 2019. Due to the different lifetimes of atmospheric pollutants, 72 h and 24 h backward trajectories were calculated for aerosols and NO2, respectively, with a time interval of 1 h.
Figure 7 and Table 3, respectively, summarize the clustering pathways and clustering characteristics in the bottom layer (50 m), lower layer (200 m), middle layer (600 m), and upper layer (1000 m) at Wangdu station during the observation period. For aerosols (Figure 7a–d), the maximum aerosol extinction was reached in Cluster 3 (0.429 ± 0.081 km−1) at 50 m, Cluster 4 (0.501 ± 0.102 km−1) at 200 m, Cluster 4 (0.377 ± 0.062 km−1) at 600 m, and Cluster 3 (0.174 ± 0.021 km−1) at 1000 m. Cluster 3 at 50 m represented long-distance trajectories from the Inner Mongolia Plateau. Cluster 4 at 200 m was a short track from the border of Shanxi and Hebei. Cluster 4 at 600 m and Cluster 3 at 1000 m originated from the Yellow Sea and passed through Shandong province into Hebei province. For NO2 (Figure 7e–h), the maximum concentrations were found in Cluster 4 (1.967 ± 0.313 ppb) at 50 m, Cluster 3 (2.278 ± 0.302 ppb) at 200 m, Cluster 2 (0.192 ± 0.035 ppb) at 600 m, and Cluster 3 (0.042 ± 0.013 ppb) at 1000 m. Among them, Cluster 4 at 50 m was a short-distance transmission, originating from central Hebei province. Cluster 3 at 200 m originated from the Bohai Sea coast and passed through Tianjin province, similar to the cluster observed at the nearby Gucheng site [92]. These results indicated that short-range transport from central Hebei and mid-range transport from the Bohai Sea coast significantly impacted the NO2 pollution at Wangdu station.
In addition, the WPSCF values of aerosol and NO2 were analyzed at different altitudes to characterize potential pollutant sources. A higher WPSCF value indicated a more significant contribution of pollutants from the region to the observation station. In this study, regions with WPSCF values between 0.3 and 0.4 were defined as having low pollution contribution, between 0.4 and 0.5 as middle pollution contribution, and WPSCF values greater than 0.5 as high pollution contribution regions [93,94].
As shown in Figure 8, pollutants at different altitudes generally originated from the same source directions, but potential pollution sources were more widely distributed in mid- to high-altitude regions compared to low-altitude areas, with similar phenomena also observed in Xingtai [51]. High pollution contribution regions for aerosol were mainly distributed in the Shanxi, Shandong, and Jiangsu provinces and the southern part of Hebei province. The predominant potential aerosol source below 200 m was Shanxi province, known for its abundant coal resources, where coal combustion and biomass burning significantly contribute to PM2.5 emissions [95]. Above 500 m, aerosols primarily originated from the southeast, especially in the Shandong and Jiangsu provinces. High pollution contribution regions for NO2 were mainly distributed in the Beijing–Tianjin–Hebei region, due to heavy industrial emissions. The combustion of fossil fuels in heavy industries was an important source of NO2 emissions [96]. In addition, Beijing and Tianjin, with high population densities and heavy traffic flow, contributed to high NO2 concentrations [97,98]. In the high-altitude region above 1000 m, the southern part of Shanxi and the Bohai Sea were also high pollution contribution regions for NO2.

3.5. Heterogeneous Contribution to HONO from NO2 in Agricultural Regions in Summer

The heterogeneous reactions of NO2 on wet surfaces are also an important source of HONO. Data indicating an increase in the ratio of HONO to NO2 over a continuous three-hour period between 8:00 and 18:00 were used to calculate the conversion rate from NO2 to HONO.
Figure 9 describes the variation of C H O N O and the aerosol extinction coefficient with altitude during the summer months. Table S2 shows the monthly average temperature and humidity during the summer. We found that the mean C H O N O near the surface varied in the order of June (0.0226 h−1) > July (0.0208 h−1) > August (0.0186 h−1). C H O N O decreased with altitude between 0 and 150 m, reaching the first low point at 150 m. This may be attributed to the neglect of HONO emissions from agricultural sources, leading to an overestimation of near-surface C H O N O [66]. C H O N O increased with altitude from 150 m to 500 m. The maximum mean C H O N O occurred at 500 m, varying in the order of August (0.0395 h−1) > June (0.0333 h−1) > July (0.0324 h−1). C H O N O continued to decrease with altitude above 500 m. The minimum mean C H O N O occurred at the highest altitude of 2000 m, varying in the order of June (0.0022 h−1) > July (0.0018 h−1) > August (0.0015 h−1). The maximum aerosol extinction also occurred near 500 m and decreased continuously with altitude above 500 m, consistent with the trend of C H O N O .

3.6. OH Radicals Produced from HONO Photolysis at Different Altitudes

HONO photolysis is one of the important sources of atmospheric OH radicals, affecting the atmospheric oxidation capacity [2,10]. To evaluate the role of HONO in the production of OH radicals in agricultural regions, the photolysis frequency of HONO and its contribution to the OH generation rate were determined using the TUV model.
The TUV model performs best on cloudless and non-haze days, so data on such days were selected for analysis. The monthly average HONO photolysis frequency and OH generation rate during the observation period are shown in Figure 10. The average HONO photolysis frequency near the surface varied in the order of June (0.00124 s−1) > July (0.00122 s−1) > August (0.00115 s−1). The HONO photolysis frequency increased with altitude, typically reaching its peak around 12:00 to 13:00 daily. In contrast to the trend of photolysis frequency, the OH production rate decreased with altitude, with peaks occurring in the morning and late afternoon.

4. Discussion

In the time series of trace gases in Wangdu (Figure 4), we observed a peak in HONO concentrations on 20 August, reaching approximately 0.49 ppb, even though the NO2 levels on that day were not high. Further analysis of the diurnal variation in HONO on 20 August, as shown in Figure S5, indicated that high levels of HONO occurred in the morning (8:00–9:00), at noon (11:00–13:00), and in the evening (16:00–18:00). Similar diurnal variation patterns were observed at the same site (Wangdu) during the fertilization period [99]. Therefore, we suspect this phenomenon may be related to local agricultural activities.
The peak levels of NO2 occurred at 8:00–10:00 and 16:00–18:00 (Figure 5e–g), even though the station is located in an agricultural area. This is consistent with previous observations in the rural areas of Gucheng [92] and Raoyang [100] in the NCP. On the one hand, despite being situated in an agricultural area, the observation site experienced significant NO2 emissions from a major highway located nearby to the northwest (Figure S6). Moreover, weak photolysis of NO2 during the night led to an accumulation of NO2 after sunset, resulting in peak concentrations in the morning [85]. On the other hand, the low boundary layer in the morning and evening hindered the diffusion of NO2 [39].
Above 500 m, the trend of the C H O N O profile was consistent with aerosols (Figure 9). It can be inferred that above 500 m, heterogeneous reactions on aerosol surfaces were one of the important HONO sources. Furthermore, we analyzed the correlation between the ratio of HONO to NO2 and aerosols at different altitudes, with the results shown in Figure S7. The results indicate that the correlation between the ratio of HONO to NO2 and aerosols above 500 m was significantly stronger than that below 500 m, further supporting this inference.
The production and loss of OH may involve complex chemical reactions, such as the photolysis of O3 and HCHO leading to OH formation [101], and the homogeneous reaction between OH and NO resulting in OH loss [102]. This study aims to provide a preliminary quantitative analysis of the relationship between OH production from photolysis, HONO concentration, and photolysis rates, with a particular focus on the contribution of HONO to atmospheric OH levels. The summer average OH production rate near the surface was 0.297 ppb/h, slightly higher than that reported in the urban regions of Hefei (0.287 ppb/h) and lower than that in the rural regions of Hefei (0.36 ppb/h) [66]. The maximum OH production rate near the surface in summer was 0.872 ppb/h, slightly higher than that in Beijing (approximately 0.700 ppb/h) [5]. The OH production rate decreased with altitude, with peaks in the morning and late afternoon, consistent with the variation in HONO concentration, contrary to the HONO photolysis rate. This trend indicated that the OH production rate from HONO photolysis was mainly controlled by HONO concentration.

5. Conclusions

MAX-DOAS observations were conducted in the agricultural regions of the NCP to measure vertical distributions of aerosols, NO2, and HONO. MAX-DOAS NO2 VCD and AOD showed correlations of 0.77 and 0.87, respectively, with TROPOMI NO2 VCD and MODIS AOD. The vertical profiles of aerosols and NO2 exhibited a near-Gaussian distribution. Aerosol extinction peaked typically before 10:00 in the morning and after 17:00, influenced by diurnal variations in the boundary layer height. NO2 peaks appeared at 8:00–10:00 and 16:00–18:00. HONO reached its maximum concentration near the surface around 8:00 due to nighttime accumulation and exponentially decreased with height. After sunrise, HONO concentration rapidly reduced due to photolysis.
Potential source analysis suggested that the main potential sources of aerosol pollution were in the southern part of the Hebei, Shanxi, Shandong, and Jiangsu provinces. The primary potential pollution sources for NO2 were in the Beijing–Tianjin–Hebei region. The south of Shanxi province and the Bohai Sea area were also significant potential pollution sources for NO2 above 1000 m. Additionally, the NO2 to HONO conversion rate through heterogeneous reactions was calculated, revealing that the maximum values of C(HONO) and aerosol extinction occurred near 500 m, and decreased continuously with altitude above 500 m. This suggested that the heterogeneous reactions of NO2 on aerosol surfaces play a significant role as sources of HONO at altitudes above 500 m.
HONO serves as a crucial source of OH in the atmosphere. In this study, we observed that the maximum photolysis frequency of HONO consistently occurred between 12:00 and 13:00 each day. The production rate of OH from HONO photolysis exhibited a decrease with altitude, with peaks occurring in the morning and late afternoon, which correlated with the variation in HONO concentration. This indicated that the OH production rate was primarily controlled by HONO concentration.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/rs16173192/s1, Section S1. The Cloud effect. Section S2. The detailed introduction to the Ring effect and the I0 effect. Section S3. WRF model configurations. Table S1. Model configuration options. Table S2. The monthly average temperature and humidity during the observation period. Figure S1. Process of vertical profile retrievals. Figure S2. The time series of meteorological information (boundary layer height, temperature, relative humidity, and visibility) in the Wangdu area during the observation period. Figure S3. The scatter plots between meteorological information (RH and temperature) and pollutants (aerosols, NO2, HONO, and NO2/HONO). Figure S4. Mean profiles and error bars of aerosol extinction, NO2 and HONO concentrations, and growth percentages during haze and non-haze days. Figure S5. (a) Vertical diurnal variation of HONO on 20 August. (b) Diurnal variation of HONO within 0–100 m altitude on 20 August. Figure S6. Surrounding environment of the observation site. Figure S7. The scatter plots between the ratio of HONO to NO2 and aerosols at different altitudes [103,104,105,106,107,108,109,110,111,112].

Author Contributions

Conceptualization, S.W., Q.H., C.L. (Cheng Liu), and C.X.; Validation, S.W., X.W. and J.X.; Writing—original draft, S.W.; Writing—review and editing, S.W., C.X., Q.H. and A.T.; Software, T.L., J.F., C.L. (Chao Liu), T.L. and J.C.; Formal analysis, X.W., W.T., C.X. and Q.H.; Resources, C.L. (Cheng Liu). All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (U21A2027, 42207113, 42307135), the National Key Research and Development Program of China (2023YFC3705300), the Natural Science Fund of Jiangsu Province (BK20210485), and the Wuxi University Research Start-up Fund for Introduced Talents (2024r040).

Data Availability Statement

QDOAS spectral fitting software (http://uv-vis.aeronomie.be/software/QDOAS/, accessed on 22 March 2024). The HYSPLIT (Hybrid Single particle Lagrangian Integrated Trajectory) model (https://www.ready.noaa.gov/HYSPLIT.php, accessed on 22 March 2023). The NCAR Tropospheric Ultraviolet and Visible (TUV) radiation model (https://www2.acom.ucar.edu/modeling/tropospheric-ultraviolet-and-visible-tuv-radiation-model, accessed on 22 March 2024). Visibility and humidity information from Shijiazhuang Zhengding International Airport (https://www.wunderground.com/weather/ZBSJ, accessed on 21 May 2024). The MODIS AOD data (https://modis.gsfc.nasa.gov/data/, accessed on 10 May 2024). The Weather Research and Forecasting (WRF) model (https://www.mmm.ucar.edu/models/wrf, accessed on 10 May 2024).

Acknowledgments

We acknowledge the Belgian Institute for Space Aeronomy (BIRAIASB), Brussels, Belgium, for their freely accessible QDOAS software. Additionally, we thank the National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory (ARL) for the HYSPLIT transport and dispersion model. We appreciate the help from IUP Heidelberg University on HEIPRO. We recognize the great contributions of the Weather Research and Forecasting (WRF) model, the SCIATRAN radiative transfer model, and the NCAR Tropospheric Ultraviolet and Visible (TUV) radiation model. We also thank the Research Center of Eco-Environmental Sciences, Chinese Academy of Sciences, for maintaining the MAX-DOAS instrument.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Garcia-Nieto, D.; Benavent, N.; Saiz-Lopez, A. Measurements of Atmospheric HONO Vertical Distribution and Temporal Evolution in Madrid (Spain) Using the MAX-DOAS Technique. Sci. Total Environ. 2018, 643, 957–966. [Google Scholar] [CrossRef] [PubMed]
  2. Huang, R.-J.; Yang, L.; Cao, J.; Wang, Q.; Tie, X.; Ho, K.-F.; Shen, Z.; Zhang, R.; Li, G.; Zhu, C.; et al. Concentration and Sources of Atmospheric Nitrous Acid (HONO) at an Urban Site in Western China. Sci. Total Environ. 2017, 593–594, 165–172. [Google Scholar] [CrossRef]
  3. Michoud, V.; Kukui, A.; Camredon, M.; Colomb, A.; Borbon, A.; Miet, K.; Aumont, B.; Beekmann, M.; Durand-Jolibois, R.; Perrier, S.; et al. Radical Budget Analysis in a Suburban European Site during the MEGAPOLI Summer Field Campaign. Atmos. Chem. Phys. 2012, 12, 11951–11974. [Google Scholar] [CrossRef]
  4. Perner, D.; Platt, U. Detection of Nitrous Acid in the Atmosphere by Differential Optical Absorption. Geophys. Res. Lett. 1979, 6, 917–920. [Google Scholar] [CrossRef]
  5. Hendrick, F.; Müller, J.-F.; Clémer, K.; Wang, P.; De Mazière, M.; Fayt, C.; Gielen, C.; Hermans, C.; Ma, J.Z.; Pinardi, G.; et al. Four Years of Ground-Based MAX-DOAS Observations of HONO and NO2 in the Beijing Area. Atmos. Chem. Phys. 2014, 14, 765–781. [Google Scholar] [CrossRef]
  6. Su, H.; Cheng, Y.F.; Cheng, P.; Zhang, Y.H.; Dong, S.; Zeng, L.M.; Wang, X.; Slanina, J.; Shao, M.; Wiedensohler, A. Observation of Nighttime Nitrous Acid (HONO) Formation at a Non-Urban Site during PRIDE-PRD2004 in China. Atmos. Environ. 2008, 42, 6219–6232. [Google Scholar] [CrossRef]
  7. Acker, K.; Möller, D.; Wieprecht, W.; Meixner, F.X.; Bohn, B.; Gilge, S.; Plass-Dülmer, C.; Berresheim, H. Strong Daytime Production of OH from HNO2 at a Rural Mountain Site. Geophys. Res. Lett. 2006, 33, L02809. [Google Scholar] [CrossRef]
  8. Song, Y.; Zhang, Y.; Xue, C.; Liu, P.; He, X.; Li, X.; Mu, Y. The Seasonal Variations and Potential Sources of Nitrous Acid (HONO) in the Rural North China Plain. Environ. Pollut. 2022, 311, 119967. [Google Scholar] [CrossRef]
  9. Huang, R.-J.; Zhang, Y.; Bozzetti, C.; Ho, K.-F.; Cao, J.-J.; Han, Y.; Daellenbach, K.R.; Slowik, J.G.; Platt, S.M.; Canonaco, F.; et al. High Secondary Aerosol Contribution to Particulate Pollution during Haze Events in China. Nature 2014, 514, 218–222. [Google Scholar] [CrossRef] [PubMed]
  10. Lelieveld, J.; Dentener, F.J.; Peters, W.; Krol, M.C. On the Role of Hydroxyl Radicals in the Self-Cleansing Capacity of the Troposphere. Atmos. Chem. Phys. 2004, 4, 2337–2344. [Google Scholar] [CrossRef]
  11. Jarvis, D.L.; Leaderer, B.P.; Chinn, S.; Burney, P.G. Indoor Nitrous Acid and Respiratory Symptoms and Lung Function in Adults. Thorax 2005, 60, 474–479. [Google Scholar] [CrossRef]
  12. Ohyama, M.; Nakajima, T.; Minejima, C.; Azuma, K.; Oka, K.; Itano, Y.; Kudo, S.; Takenaka, N. Association between Indoor Nitrous Acid, Outdoor Nitrogen Dioxide, and Asthma Attacks: Results of a Pilot Study. Int. J. Environ. Health Res. 2019, 29, 632–642. [Google Scholar] [CrossRef] [PubMed]
  13. Sleiman, M.; Gundel, L.A.; Pankow, J.F.; Jacob, P.; Singer, B.C.; Destaillats, H. Formation of Carcinogens Indoors by Surface-Mediated Reactions of Nicotine with Nitrous Acid, Leading to Potential Thirdhand Smoke Hazards. Proc. Natl. Acad. Sci. USA 2010, 107, 6576–6581. [Google Scholar] [CrossRef]
  14. Ge, M.; Tong, S.; Wang, W.; Zhang, W.; Chen, M.; Peng, C.; Li, J.; Zhou, L.; Chen, Y.; Liu, M. Important Oxidants and Their Impact on the Environmental Effects of Aerosols. J. Phys. Chem. A 2021, 125, 3813–3825. [Google Scholar] [CrossRef] [PubMed]
  15. Li, X.; Brauers, T.; Häseler, R.; Bohn, B.; Fuchs, H.; Hofzumahaus, A.; Holland, F.; Lou, S.; Lu, K.D.; Rohrer, F.; et al. Exploring the Atmospheric Chemistry of Nitrous Acid (HONO) at a Rural Site in Southern China. Atmos. Chem. Phys. 2012, 12, 1497–1513. [Google Scholar] [CrossRef]
  16. Su, H.; Cheng, Y.F.; Shao, M.; Gao, D.F.; Yu, Z.Y.; Zeng, L.M.; Slanina, J.; Zhang, Y.H.; Wiedensohler, A. Nitrous Acid (HONO) and Its Daytime Sources at a Rural Site during the 2004 PRIDE-PRD Experiment in China. J. Geophys. Res. Atmos. 2008, 113, D14312. [Google Scholar] [CrossRef]
  17. Calvert, J.G.; Yarwood, G.; Dunker, A.M. An Evaluation of the Mechanism of Nitrous Acid Formation in the Urban Atmosphere. Res. Chem. Intermed. 1994, 20, 463–502. [Google Scholar] [CrossRef]
  18. Li, S.; Song, W.; Zhan, H.; Zhang, Y.; Zhang, X.; Li, W.; Tong, S.; Pei, C.; Wang, Y.; Chen, Y.; et al. Contribution of Vehicle Emission and NO2 Surface Conversion to Nitrous Acid (HONO) in Urban Environments: Implications from Tests in a Tunnel. Environ. Sci. Technol. 2021, 55, 15616–15624. [Google Scholar] [CrossRef]
  19. Oswald, R.; Behrendt, T.; Ermel, M.; Wu, D.; Su, H.; Cheng, Y.; Breuninger, C.; Moravek, A.; Mougin, E.; Delon, C.; et al. HONO Emissions from Soil Bacteria as a Major Source of Atmospheric Reactive Nitrogen. Science 2013, 341, 1233–1235. [Google Scholar] [CrossRef]
  20. Stemmler, K.; Ammann, M.; Donders, C.; Kleffmann, J.; George, C. Photosensitized Reduction of Nitrogen Dioxide on Humic Acid as a Source of Nitrous Acid. Nature 2006, 440, 195–198. [Google Scholar] [CrossRef]
  21. Han, C.; Yang, W.; Wu, Q.; Yang, H.; Xue, X. Heterogeneous Photochemical Conversion of NO2 to HONO on the Humic Acid Surface under Simulated Sunlight. Environ. Sci. Technol. 2016, 50, 5017–5023. [Google Scholar] [CrossRef]
  22. Su, H.; Cheng, Y.; Oswald, R.; Behrendt, T.; Trebs, I.; Meixner, F.X.; Andreae, M.O.; Cheng, P.; Zhang, Y.; Pöschl, U. Soil Nitrite as a Source of Atmospheric HONO and OH Radicals. Science 2011, 333, 1616–1618. [Google Scholar] [CrossRef]
  23. Wang, S.; Zhou, R.; Zhao, H.; Wang, Z.; Chen, L.; Zhou, B. Long-Term Observation of Atmospheric Nitrous Acid (HONO) and Its Implication to Local NO2 Levels in Shanghai, China. Atmos. Environ. 2013, 77, 718–724. [Google Scholar] [CrossRef]
  24. Cui, L.; Li, R.; Zhang, Y.; Meng, Y.; Fu, H.; Chen, J. An Observational Study of Nitrous Acid (HONO) in Shanghai, China: The Aerosol Impact on HONO Formation during the Haze Episodes. Sci. Total Environ. 2018, 630, 1057–1070. [Google Scholar] [CrossRef] [PubMed]
  25. Meng, F.; Qin, M.; Fang, W.; Duan, J.; Tang, K.; Zhang, H.; Shao, D.; Liao, Z.; Feng, Y.; Huang, Y.; et al. Measurement of HONO Flux Using the Aerodynamic Gradient Method over an Agricultural Field in the Huaihe River Basin, China. Atmos. Chem. Complex Air Pollut. 2022, 114, 297–307. [Google Scholar] [CrossRef]
  26. Tang, K.; Qin, M.; Duan, J.; Fang, W.; Meng, F.; Liang, S.; Xie, P.; Liu, J.; Liu, W.; Xue, C.; et al. A Dual Dynamic Chamber System Based on IBBCEAS for Measuring Fluxes of Nitrous Acid in Agricultural Fields in the North China Plain. Atmos. Environ. 2019, 196, 10–19. [Google Scholar] [CrossRef]
  27. Xue, C.; Ye, C.; Zhang, C.; Catoire, V.; Liu, P.; Gu, R.; Zhang, J.; Ma, Z.; Zhao, X.; Zhang, W.; et al. Evidence for Strong HONO Emission from Fertilized Agricultural Fields and Its Remarkable Impact on Regional O3 Pollution in the Summer North China Plain. ACS Earth Space Chem. 2021, 5, 340–347. [Google Scholar] [CrossRef]
  28. Meng, F.; Qin, M.; Tang, K.; Duan, J.; Fang, W.; Liang, S.; Ye, K.; Xie, P.; Sun, Y.; Xie, C.; et al. High-Resolution Vertical Distribution and Sources of HONO and NO2 in the Nocturnal Boundary Layer in Urban Beijing, China. Atmos. Chem. Phys. 2020, 20, 5071–5092. [Google Scholar] [CrossRef]
  29. Acker, K.; Febo, A.; Trick, S.; Perrino, C.; Bruno, P.; Wiesen, P.; Möller, D.; Wieprecht, W.; Auel, R.; Giusto, M.; et al. Nitrous Acid in the Urban Area of Rome. Atmos. Environ. 2006, 40, 3123–3133. [Google Scholar] [CrossRef]
  30. Gall, E.T.; Griffin, R.J.; Steiner, A.L.; Dibb, J.; Scheuer, E.; Gong, L.; Rutter, A.P.; Cevik, B.K.; Kim, S.; Lefer, B.; et al. Evaluation of Nitrous Acid Sources and Sinks in Urban Outflow. Atmos. Environ. 2016, 127, 272–282. [Google Scholar] [CrossRef]
  31. VandenBoer, T.; Brown, S.; Murphy, J.; Keene, W.; Young, C.; Pszenny, A.; Kim, S.; Warneke, C.; de Gouw, J.; Maben, J.; et al. Understanding the Role of the Ground Surface in HONO Vertical Structure: High Resolution Vertical Profiles during NACHTT-11. J. Geophys. Res. 2013, 118, 10155–10171. [Google Scholar] [CrossRef]
  32. Zhang, N.; Zhou, X.; Shepson, P.B.; Gao, H.; Alaghmand, M.; Stirm, B. Aircraft Measurement of HONO Vertical Profiles over a Forested Region. Geophys. Res. Lett. 2009, 36, L15820. [Google Scholar] [CrossRef]
  33. Ye, C.; Zhou, X.; Pu, D.; Stutz, J.; Festa, J.; Spolaor, M.; Tsai, C.; Cantrell, C.; Mauldin III, R.L.; Weinheimer, A.; et al. Tropospheric HONO Distribution and Chemistry in the Southeastern US. Atmos. Chem. Phys. 2018, 18, 9107–9120. [Google Scholar] [CrossRef]
  34. Li, X.; Rohrer, F.; Hofzumahaus, A.; Brauers, T.; Häseler, R.; Bohn, B.; Broch, S.; Fuchs, H.; Gomm, S.; Holland, F.; et al. Missing Gas-Phase Source of HONO Inferred from Zeppelin Measurements in the Troposphere. Science 2014, 344, 292–296. [Google Scholar] [CrossRef]
  35. Young, C.J.; Washenfelder, R.A.; Roberts, J.M.; Mielke, L.H.; Osthoff, H.D.; Tsai, C.; Pikelnaya, O.; Stutz, J.; Veres, P.R.; Cochran, A.K.; et al. Vertically Resolved Measurements of Nighttime Radical Reservoirs in Los Angeles and Their Contribution to the Urban Radical Budget. Environ. Sci. Technol. 2012, 46, 10965–10973. [Google Scholar] [CrossRef] [PubMed]
  36. Wong, K.W.; Tsai, C.; Lefer, B.; Haman, C.; Grossberg, N.; Brune, W.H.; Ren, X.; Luke, W.; Stutz, J. Daytime HONO Vertical Gradients during SHARP 2009 in Houston, TX. Atmos. Chem. Phys. 2012, 12, 635–652. [Google Scholar] [CrossRef]
  37. Ryan, R.G.; Rhodes, S.; Tully, M.; Wilson, S.; Jones, N.; Frieß, U.; Schofield, R. Daytime HONO, NO2 and Aerosol Distributions from MAX-DOAS Observations in Melbourne. Atmos. Chem. Phys. 2018, 18, 13969–13985. [Google Scholar] [CrossRef]
  38. Li, L.; Lu, C.; Chan, P.-W.; Zhang, X.; Yang, H.-L.; Lan, Z.-J.; Zhang, W.-H.; Liu, Y.-W.; Pan, L.; Zhang, L. Tower Observed Vertical Distribution of PM2.5, O3 and NOx in the Pearl River Delta. Atmos. Environ. 2020, 220, 117083. [Google Scholar] [CrossRef]
  39. Kang, Y.; Tang, G.; Li, Q.; Liu, B.; Cao, J.; Hu, Q.; Wang, Y. Evaluation and Evolution of MAX-DOAS-Observed Vertical NO2 Profiles in Urban Beijing. Adv. Atmos. Sci. 2021, 38, 1188–1196. [Google Scholar] [CrossRef]
  40. Wagner, T.; Dix, B.; Friedeburg, C.V.; Frieß, U.; Sanghavi, S.; Sinreich, R.; Platt, U. MAX-DOAS O4 Measurements: A New Technique to Derive Information on Atmospheric Aerosols—Principles and Information Content. J. Geophys. Res. Atmos. 2004, 109, D22205. [Google Scholar] [CrossRef]
  41. Xing, C.; Xu, S.; Song, Y.; Liu, C.; Liu, Y.; Lu, K.; Tan, W.; Zhang, C.; Hu, Q.; Wang, S.; et al. A New Insight into the Vertical Differences in NO2 Heterogeneous Reaction to Produce HONO over Inland and Marginal Seas. Atmos. Chem. Phys. 2023, 23, 5815–5834. [Google Scholar] [CrossRef]
  42. Bobrowski, N.; Hönninger, G.; Galle, B.; Platt, U. Detection of Bromine Monoxide in a Volcanic Plume. Nature 2003, 423, 273–276. [Google Scholar] [CrossRef] [PubMed]
  43. Vlemmix, T.; Hendrick, F.; Pinardi, G.; De Smedt, I.; Fayt, C.; Hermans, C.; Piters, A.; Wang, P.; Levelt, P.; Van Roozendael, M. MAX-DOAS Observations of Aerosols, Formaldehyde and Nitrogen Dioxide in the Beijing Area: Comparison of Two Profile Retrieval Approaches. Atmos. Meas. Tech. 2015, 8, 941–963. [Google Scholar] [CrossRef]
  44. Xing, C.; Liu, C.; Hu, Q.; Fu, Q.; Wang, S.; Lin, H.; Zhu, Y.; Wang, S.; Wang, W.; Javed, Z.; et al. Vertical Distributions of Wintertime Atmospheric Nitrogenous Compounds and the Corresponding OH Radicals Production in Leshan, Southwest China. J. Environ. Sci. 2021, 105, 44–55. [Google Scholar] [CrossRef]
  45. Xue, C.; Ye, C.; Zhang, Y.; Ma, Z.; Liu, P.; Zhang, C.; Zhao, X.; Liu, J.; Mu, Y. Development and Application of a Twin Open-Top Chambers Method to Measure Soil HONO Emission in the North China Plain. Sci. Total Environ. 2019, 659, 621–631. [Google Scholar] [CrossRef]
  46. Javed, Z.; Tanvir, A.; Bilal, M.; Su, W.; Xia, C.; Rehman, A.; Zhang, Y.; Sandhu, O.; Xing, C.; Ji, X.; et al. Recommendations for HCHO and SO2 Retrieval Settings from MAX-DOAS Observations under Different Meteorological Conditions. Remote Sens. 2021, 13, 2244. [Google Scholar] [CrossRef]
  47. Ji, X.; Liu, C.; Wang, Y.; Hu, Q.; Lin, H.; Zhao, F.; Xing, C.; Tang, G.; Zhang, J.; Wagner, T. Ozone Profiles without Blind Area Retrieved from MAX-DOAS Measurements and Comprehensive Validation with Multi-Platform Observations. Remote Sens. Environ. 2023, 284, 113339. [Google Scholar] [CrossRef]
  48. Lampel, J.; Wang, Y.; Hilboll, A.; Beirle, S.; Sihler, H.; Puķīte, J.; Platt, U.; Wagner, T. The Tilt Effect in DOAS Observations. Atmos. Meas. Tech. 2017, 10, 4819–4831. [Google Scholar] [CrossRef]
  49. Xing, C.; Liu, C.; Wang, S.; Chan, K.L.; Gao, Y.; Huang, X.; Su, W.; Zhang, C.; Dong, Y.; Fan, G.; et al. Observations of the Vertical Distributions of Summertime Atmospheric Pollutants and the Corresponding Ozone Production in Shanghai, China. Atmos. Chem. Phys. 2017, 17, 14275–14289. [Google Scholar] [CrossRef]
  50. Wagner, T.; Apituley, A.; Beirle, S.; Dörner, S.; Friess, U.; Remmers, J.; Shaiganfar, R. Cloud Detection and Classification Based on MAX-DOAS Observations. Atmos. Meas. Tech. 2014, 7, 1289–1320. [Google Scholar] [CrossRef]
  51. Wang, Y.; Dörner, S.; Donner, S.; Böhnke, S.; De Smedt, I.; Dickerson, R.R.; Dong, Z.; He, H.; Li, Z.; Li, Z.; et al. Vertical Profiles of NO2, SO2, HONO, HCHO, CHOCHO and Aerosols Derived from MAX-DOAS Measurements at a Rural Site in the Central Western North China Plain and Their Relation to Emission Sources and Effects of Regional Transport. Atmos. Chem. Phys. 2019, 19, 5417–5449. [Google Scholar] [CrossRef]
  52. Jin, J.; Zhang, Y.; Wang, Y.; Zhou, Q.; Lv, S.; Ma, J. Aerosol and Cloud Classifications Derived from MAX-DOAS Measurements in Urban North China and Their Comparisons to Multiple Remote Sensing Datasets. In Proceedings of the 2019 International Conference on Meteorology Observations (ICMO), Chengdu, China, 28–31 December 2019; pp. 1–4. [Google Scholar]
  53. Wagner, T.; Senne, T.; Erle, F.; Otten, C.; Stutz, J.; Pfeilsticker, K.; Platt, U. Determination of Cloud Properties and Cloud Type from DOAS-Measurements; Varotsos, C., Ed.; Springer: Berlin/Heidelberg, Germany, 1997; Volume 53, pp. 327–336. [Google Scholar]
  54. Vandaele, A.C.; Hermans, C.; Simon, P.C.; Carleer, M.; Colin, R.; Fally, S.; Mérienne, M.F.; Jenouvrier, A.; Coquart, B. Measurements of the NO2 Absorption Cross-Section from 42,000 cm−1 to 10,000 cm−1 (238–1000 Nm) at 220 K and 294 K. J. Quant. Spectrosc. Radiat. Transf. 1998, 59, 171–184. [Google Scholar] [CrossRef]
  55. Serdyuchenko, A.; Gorshelev, V.; Weber, M.; Chehade, W.; Burrows, J.P. High Spectral Resolution Ozone Absorption Cross-Sections—Part 2: Temperature Dependence. Atmos. Meas. Tech. 2014, 7, 625–636. [Google Scholar] [CrossRef]
  56. Thalman, R.; Volkamer, R. Temperature Dependent Absorption Cross-Sections of O2–O2 Collision Pairs between 340 and 630 Nm and at Atmospherically Relevant Pressure. Phys. Chem. Chem. Phys. 2013, 15, 15371–15381. [Google Scholar] [CrossRef] [PubMed]
  57. Meller, R.; Moortgat, G.K. Temperature Dependence of the Absorption Cross Sections of Formaldehyde between 223 and 323 K in the Wavelength Range 225–375 Nm. J. Geophys. Res. Atmos. 2000, 105, 7089–7101. [Google Scholar] [CrossRef]
  58. Rothman, L.S.; Gordon, I.E.; Barber, R.J.; Dothe, H.; Gamache, R.R.; Goldman, A.; Perevalov, V.I.; Tashkun, S.A.; Tennyson, J. HITEMP, the High-Temperature Molecular Spectroscopic Database. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 2139–2150. [Google Scholar] [CrossRef]
  59. Fleischmann, O.C.; Hartmann, M.; Burrows, J.P.; Orphal, J. New Ultraviolet Absorption Cross-Sections of BrO at Atmospheric Temperatures Measured by Time-Windowing Fourier Transform Spectroscopy. J. Photochem. Photobiol. Chem. 2004, 168, 117–132. [Google Scholar] [CrossRef]
  60. Stutz, J.; Kim, E.S.; Platt, U.; Bruno, P.; Perrino, C.; Febo, A. UV-Visible Absorption Cross Sections of Nitrous Acid. J. Geophys. Res. Atmos. 2000, 105, 14585–14592. [Google Scholar] [CrossRef]
  61. Chance, K.; Kurucz, R.L. An Improved High-Resolution Solar Reference Spectrum for Earth’s Atmosphere Measurements in the Ultraviolet, Visible, and near Infrared. J. Quant. Spectrosc. Radiat. Transf. 2010, 111, 1289–1295. [Google Scholar] [CrossRef]
  62. Alicke, B.; Platt, U.; Stutz, J. Impact of Nitrous Acid Photolysis on the Total Hydroxyl Radical Budget during the Limitation of Oxidant Production/Pianura Padana Produzione Di Ozono Study in Milan. J. Geophys. Res. Atmos. 2002, 107, LOP 9-1–LOP 9-17. [Google Scholar] [CrossRef]
  63. Frieß, U.; Monks, P.S.; Remedios, J.J.; Rozanov, A.; Sinreich, R.; Wagner, T.; Platt, U. MAX-DOAS O4 Measurements: A New Technique to Derive Information on Atmospheric Aerosols: 2. Modeling Studies. J. Geophys. Res. Atmos. 2006, 111, D14203. [Google Scholar] [CrossRef]
  64. Frieß, U.; Sihler, H.; Sander, R.; Pöhler, D.; Yilmaz, S.; Platt, U. The Vertical Distribution of BrO and Aerosols in the Arctic: Measurements by Active and Passive Differential Optical Absorption Spectroscopy. J. Geophys. Res. Atmos. 2011, 116, D00R04. [Google Scholar] [CrossRef]
  65. Rozanov, A.; Rozanov, V.; Buchwitz, M.; Kokhanovsky, A.; Burrows, J.P. SCIATRAN 2.0—A New Radiative Transfer Model for Geophysical Applications in the 175–2400 nm Spectral Region. Adv. Space Res. 2005, 36, 1015–1019. [Google Scholar] [CrossRef]
  66. Xing, C.; Liu, C.; Li, Q.; Wang, S.; Tan, W.; Zou, T.; Wang, Z.; Lu, C. Observations of HONO and Its Precursors between Urban and Its Surrounding Agricultural Fields: The Vertical Transports, Sources and Contribution to OH. Sci. Total Environ. 2024, 915, 169159. [Google Scholar] [CrossRef] [PubMed]
  67. Chong, X.; Wang, Y.; Liu, R.; Zhang, Y.; Zhang, Y.; Zheng, W. Pollution Characteristics and Source Difference of Gaseous Elemental Mercury between Haze and Non-Haze Days in Winter. Sci. Total Environ. 2019, 678, 671–680. [Google Scholar] [CrossRef]
  68. Javed, Z.; Liu, C.; Khokhar, M.F.; Xing, C.; Tan, W.; Subhani, M.A.; Rehman, A.; Tanvir, A. Investigating the Impact of Glyoxal Retrieval from MAX-DOAS Observations during Haze and Non-Haze Conditions in Beijing. J. Environ. Sci. 2019, 80, 296–305. [Google Scholar] [CrossRef]
  69. Hong, Q.; Xie, Z.; Liu, C.; Wang, F.; Xie, P.; Kang, H.; Xu, J.; Wang, J.; Wu, F.; He, P.; et al. Speciated Atmospheric Mercury on Haze and Non-Haze Days in an Inland City in China. Atmos. Chem. Phys. 2016, 16, 13807–13821. [Google Scholar] [CrossRef]
  70. Wang, Y.Q.; Zhang, X.Y.; Arimoto, R. The Contribution from Distant Dust Sources to the Atmospheric Particulate Matter Loadings at XiAn, China during Spring. Sci. Total Environ. 2006, 368, 875–883. [Google Scholar] [CrossRef]
  71. Ren, B.; Xie, P.; Xu, J.; Li, A.; Tian, X.; Hu, Z.; Huang, Y.; Li, X.; Zhang, Q.; Ren, H.; et al. Use of the PSCF Method to Analyze the Variations of Potential Sources and Transports of NO2, SO2, and HCHO Observed by MAX-DOAS in Nanjing, China during 2019. Sci. Total Environ. 2021, 782, 146865. [Google Scholar] [CrossRef]
  72. Wang, Y.Q.; Zhang, X.Y.; Draxler, R.R. TrajStat: GIS-Based Software That Uses Various Trajectory Statistical Analysis Methods to Identify Potential Sources from Long-Term Air Pollution Measurement Data. Environ. Model. Softw. 2009, 24, 938–939. [Google Scholar] [CrossRef]
  73. Wang, Y.Q. MeteoInfo: GIS Software for Meteorological Data Visualization and Analysis. Meteorol. Appl. 2014, 21, 360–368. [Google Scholar] [CrossRef]
  74. Stein, A.F.; Draxler, R.R.; Rolph, G.D.; Stunder, B.J.B.; Cohen, M.D.; Ngan, F. NOAA’s HYSPLIT Atmospheric Transport and Dispersion Modeling System. Bull. Am. Meteorol. Soc. 2015, 96, 2059–2077. [Google Scholar] [CrossRef]
  75. Madronich, S.; Flocke, S. The Role of Solar Radiation in Atmospheric Chemistry. In Environmental Photochemistry; Boule, P., Ed.; Springer: Berlin/Heidelberg, Germany, 1999; pp. 1–26. ISBN 978-3-540-69044-3. [Google Scholar]
  76. Xing, C.; Liu, C.; Wang, S.; Hu, Q.; Liu, H.; Tan, W.; Zhang, W.; Li, B.; Liu, J. A New Method to Determine the Aerosol Optical Properties from Multiple-Wavelength O4 Absorptions by MAX-DOAS Observation. Atmos. Meas. Tech. 2019, 12, 3289–3302. [Google Scholar] [CrossRef]
  77. Xu, S.; Wang, S.; Xia, M.; Lin, H.; Xing, C.; Ji, X.; Su, W.; Tan, W.; Liu, C.; Hu, Q. Observations by Ground-Based MAX-DOAS of the Vertical Characters of Winter Pollution and the Influencing Factors of HONO Generation in Shanghai, China. Remote Sens. 2021, 13, 3518. [Google Scholar] [CrossRef]
  78. Cai, K.; Li, S.; Lai, J.; Xia, Y.; Wang, Y.; Hu, X.; Li, A. Evaluation of TROPOMI and OMI Tropospheric NO2 Products Using Measurements from MAX-DOAS and State-Controlled Stations in the Jiangsu Province of China. Atmosphere 2022, 13, 886. [Google Scholar] [CrossRef]
  79. Zhu, Y.; Liu, C.; Hu, Q.; Teng, J.; You, D.; Zhang, C.; Ou, J.; Liu, T.; Lin, J.; Xu, T.; et al. Impacts of TROPOMI-Derived NOX Emissions on NO2 and O3 Simulations in the NCP during COVID-19. ACS Environ. Au 2022, 2, 441–454. [Google Scholar] [CrossRef] [PubMed]
  80. Peña, A.; Gryning, S.-E.; Hahmann, A.N. Observations of the Atmospheric Boundary Layer Height under Marine Upstream Flow Conditions at a Coastal Site. J. Geophys. Res. Atmos. 2013, 118, 1924–1940. [Google Scholar] [CrossRef]
  81. Jin, J.; Henzing, B.; Segers, A. How Aerosol Size Matters in Aerosol Optical Depth (AOD) Assimilation and the Optimization Using the Ångström Exponent. Atmos. Chem. Phys. 2023, 23, 1641–1660. [Google Scholar] [CrossRef]
  82. Mali, P.; Biswas, M.S.; Beirle, S.; Wagner, T.; Hulswar, S.; Inamdar, S.; Mahajan, A.S. Aerosol Measurements over India: Comparison of MAX-DOAS Measurements with Ground-Based (AERONET) and Satellite-Based (MODIS) Data. Aerosol Air Qual. Res. 2024, 24, 230076. [Google Scholar] [CrossRef]
  83. Liu, L.; Guo, J.; Miao, Y.; Liu, L.; Li, J.; Chen, D.; He, J.; Cui, C. Elucidating the Relationship between Aerosol Concentration and Summertime Boundary Layer Structure in Central China. Environ. Pollut. 2018, 241, 646–653. [Google Scholar] [CrossRef]
  84. Cheng, S.; Jin, J.; Ma, J.; Xu, X.; Ran, L.; Ma, Z.; Chen, J.; Guo, J.; Yang, P.; Wang, Y.; et al. Measuring the Vertical Profiles of Aerosol Extinction in the Lower Troposphere by MAX-DOAS at a Rural Site in the North China Plain. Atmosphere 2020, 11, 1037. [Google Scholar] [CrossRef]
  85. Iqbal, A.; Ahmad, N.; Mohy ud Din, H.; Van Roozendael, M.; Anjum, M.S.; Zeeshan Ali Khan, M.; Khokhar, M.F. Retrieval of NO2 Columns by Exploiting MAX-DOAS Observations and Comparison with OMI and TROPOMI Data during the Time Period of 2015–2019. Aerosol Air Qual. Res. 2022, 22, 210398. [Google Scholar] [CrossRef]
  86. Hu, X.; Yang, G.; Liu, Y.; Lu, Y.; Wang, Y.; Chen, H.; Chen, J.; Wang, L. Atmospheric Gaseous Organic Acids in Winter in a Rural Site of the North China Plain. J. Environ. Sci. 2022, 113, 190–203. [Google Scholar] [CrossRef]
  87. Wang, Y.; Apituley, A.; Bais, A.; Beirle, S.; Benavent, N.; Borovski, A.; Bruchkouski, I.; Chan, K.L.; Donner, S.; Drosoglou, T.; et al. Inter-Comparison of MAX-DOAS Measurements of Tropospheric HONO Slant Column Densities and Vertical Profiles during the CINDI-2 Campaign. Atmos. Meas. Tech. 2020, 13, 5087–5116. [Google Scholar] [CrossRef]
  88. Laufs, S.; Cazaunau, M.; Stella, P.; Kurtenbach, R.; Cellier, P.; Mellouki, A.; Loubet, B.; Kleffmann, J. Diurnal Fluxes of HONO above a Crop Rotation. Atmos. Chem. Phys. 2017, 17, 6907–6923. [Google Scholar] [CrossRef]
  89. Gil, J.; Kim, J.; Lee, M.; Lee, G.; Lee, D.; Jung, J.; An, J.; Hong, J.; Cho, S.; Lee, J.; et al. The Role of HONO in O3 Formation and Insight into Its Formation Mechanism during the KORUS-AQ Campaign. Atmos. Chem. Phys. Dis. 2019, 2019, 1–30. [Google Scholar] [CrossRef]
  90. He, S.; Wang, S.; Zhang, S.; Zhu, J.; Sun, Z.; Xue, R.; Zhou, B. Vertical Distributions of Atmospheric HONO and the Corresponding OH Radical Production by Photolysis at the Suburb Area of Shanghai, China. Sci. Total Environ. 2023, 858, 159703. [Google Scholar] [CrossRef] [PubMed]
  91. Meng, Z.Y.; Ding, G.A.; Xu, X.B.; Xu, X.D.; Yu, H.Q.; Wang, S.F. Vertical Distributions of SO2 and NO2 in the Lower Atmosphere in Beijing Urban Areas, China. Sci. Total Environ. 2008, 390, 456–465. [Google Scholar] [CrossRef]
  92. Jin, J.; Ma, J.; Lin, W.; Zhao, H.; Shaiganfar, R.; Beirle, S.; Wagner, T. MAX-DOAS Measurements and Satellite Validation of Tropospheric NO2 and SO2 Vertical Column Densities at a Rural Site of North China. Atmos. Environ. 2016, 133, 12–25. [Google Scholar] [CrossRef]
  93. Hong, Q.; Liu, C.; Hu, Q.; Xing, C.; Tan, W.; Liu, H.; Huang, Y.; Zhu, Y.; Zhang, J.; Geng, T.; et al. Evolution of the Vertical Structure of Air Pollutants during Winter Heavy Pollution Episodes: The Role of Regional Transport and Potential Sources. Atmos. Res. 2019, 228, 206–222. [Google Scholar] [CrossRef]
  94. Zong, Z.; Wang, X.; Tian, C.; Chen, Y.; Fu, S.; Qu, L.; Ji, L.; Li, J.; Zhang, G. PMF and PSCF Based Source Apportionment of PM2.5 at a Regional Background Site in North China. Atmos. Res. 2018, 203, 207–215. [Google Scholar] [CrossRef]
  95. Li, L.; Qi, H.; Li, X. Composition, Source Apportionment, and Health Risk of PM2.5-Bound Metals during Winter Haze in Yuci College Town, Shanxi, China. Toxics 2022, 10, 467. [Google Scholar] [CrossRef]
  96. Xiao, C.; Chang, M.; Guo, P.; Gu, M.; Li, Y. Analysis of Air Quality Characteristics of Beijing–Tianjin–Hebei and Its Surrounding Air Pollution Transport Channel Cities in China. J. Environ. Sci. 2020, 87, 213–227. [Google Scholar] [CrossRef]
  97. Liu, Y.; Zhang, Y.; Lian, C.; Yan, C.; Feng, Z.; Zheng, F.; Fan, X.; Chen, Y.; Wang, W.; Chu, B.; et al. The Promotion Effect of Nitrous Acid on Aerosol Formation in Wintertime in Beijing: The Possible Contribution of Traffic-Related Emissions. Atmos. Chem. Phys. 2020, 20, 13023–13040. [Google Scholar] [CrossRef]
  98. Cheng, N.; Li, Y.; Sun, F.; Chen, C.; Wang, B.; Li, Q.; Wei, P.; Cheng, B. Ground-Level NO2 in Urban Beijing: Trends, Distribution, and Effects of Emission Reduction Measures. Aerosol Air Qual. Res. 2018, 18, 343–356. [Google Scholar] [CrossRef]
  99. Ran, H.; An, J.; Zhang, J.; Huang, J.; Qu, Y.; Chen, Y.; Xue, C.; Mu, Y.; Liu, X. Impact of Soil–Atmosphere HONO Exchange on Concentrations of HONO and O3 in the North China Plain. Sci. Total Environ. 2024, 928, 172336. [Google Scholar] [CrossRef] [PubMed]
  100. Cheng, S.; Jin, J.; Ma, J.; Lv, J.; Liu, S.; Xu, X. Temporal Variation of NO2 and HCHO Vertical Profiles Derived from MAX-DOAS Observation in Summer at a Rural Site of the North China Plain and Ozone Production in Relation to HCHO/NO2 Ratio. Atmosphere 2022, 13, 860. [Google Scholar] [CrossRef]
  101. Nan, J.; Wang, S.; Guo, Y.; Xiang, Y.; Zhou, B. Study on the Daytime OH Radical and Implication for Its Relationship with Fine Particles over Megacity of Shanghai, China. Atmos. Environ. 2017, 154, 167–178. [Google Scholar] [CrossRef]
  102. Elshorbany, Y.F.; Kurtenbach, R.; Wiesen, P.; Lissi, E.; Rubio, M.; Villena, G.; Gramsch, E.; Rickard, A.R.; Pilling, M.J.; Kleffmann, J. Oxidation Capacity of the City Air of Santiago, Chile. Atmos. Chem. Phys. 2009, 9, 2257–2273. [Google Scholar] [CrossRef]
  103. Grainger, J.; Ring, J. Anomalous fraunhofer line profiles. Nature 1962, 193, 762. [Google Scholar] [CrossRef]
  104. Solomon, S.; Schmeltekopf, A.L.; Sanders, R.W. On the Interpretation of Zenith Sky Absorption Measurements. J. Geophys. Res. Atmos. 1987, 92, 8311–8319. [Google Scholar] [CrossRef]
  105. Platt, U.; Marquard, L.; Wagner, T.; Perner, D. Corrections for Zenith Scattered Light DOAS. Geophys. Res. Lett. 1997, 24, 1759–1762. [Google Scholar] [CrossRef]
  106. Aliwell, S.R.; Van Roozendael, M.; Johnston, P.V.; Richter, A.; Wagner, T.; Arlander, D.W.; Burrows, J.P.; Fish, D.J.; Jones, R.L.; Tørnkvist, K.K.; et al. Analysis for BrO in Zenith-Sky Spectra: An Intercomparison Exercise for Analysis Improvement. J. Geophys. Res. Atmos. 2002, 107, ACH 10-1. [Google Scholar] [CrossRef]
  107. Liu, H.; Liu, C.; Xie, Z.; Li, Y.; Huang, X.; Wang, S.; Xu, J.; Xie, P. A Paradox for Air Pollution Controlling in China Revealed by “APEC Blue” and “Parade Blue”. Sci. Rep. 2016, 6, 34408. [Google Scholar] [CrossRef] [PubMed]
  108. Chen, S.-H.; Sun, W.-Y. A One-Dimensional Time Dependent Cloud Model. J. Meteorol. Soc. Jpn. Ser II 2002, 80, 99–118. [Google Scholar] [CrossRef]
  109. Iacono, M.J.; Delamere, J.S.; Mlawer, E.J.; Shephard, M.W.; Clough, S.A.; Collins, W.D. Radiative Forcing by Long-Lived Greenhouse Gases: Calculations with the AER Radiative Transfer Models. J. Geophys. Res. Atmos. 2008, 113. [Google Scholar] [CrossRef]
  110. Grell, G.A.; Freitas, S.R. A Scale and Aerosol Aware Stochastic Convective Parameterization for Weather and Air Quality Modeling. Atmos. Chem. Phys. 2014, 14, 5233–5250. [Google Scholar] [CrossRef]
  111. Tewari, M.; Chen, F.; Wang, W.; Dudhia, J. Implementation and Verification of the Unified Noah Land Surface Model in the WRF Model. In Proceedings of the 16th Conference on Numerical Weather Prediction, Seattle, WA, USA, 12–16 January 2004. [Google Scholar]
  112. Hong, S.-Y.; Noh, Y.; Dudhia, J. A New Vertical Diffusion Package with an Explicit Treatment of Entrainment Processes. Mon. Weather Rev. 2006, 134, 2318–2341. [Google Scholar] [CrossRef]
Figure 1. Map of (a) Wangdu in North Chain Plain and (b) the campaign site. (c) The MAX-DOAS instrument used during this campaign.
Figure 1. Map of (a) Wangdu in North Chain Plain and (b) the campaign site. (c) The MAX-DOAS instrument used during this campaign.
Remotesensing 16 03192 g001
Figure 2. Typical DOAS fitting results for (a) O4, (b) NO2, and (c) HONO.
Figure 2. Typical DOAS fitting results for (a) O4, (b) NO2, and (c) HONO.
Remotesensing 16 03192 g002
Figure 3. Correlation analysis of AOD and tropospheric NO2 VCD. (a) AOD measured by MAX-DOAS and MODIS. (b) Tropospheric NO2 VCD measured by MAX-DOAS and TROPOMI.
Figure 3. Correlation analysis of AOD and tropospheric NO2 VCD. (a) AOD measured by MAX-DOAS and MODIS. (b) Tropospheric NO2 VCD measured by MAX-DOAS and TROPOMI.
Remotesensing 16 03192 g003
Figure 4. Time series of daily average (a) aerosol extinction coefficient, (b) NO2, and (c) HONO concentrations retrieved from MAX-DOAS in the near the surface (0–100 m), lower (100–300 m), middle (500–700 m), and upper (900–1100 m) layers during this campaign.
Figure 4. Time series of daily average (a) aerosol extinction coefficient, (b) NO2, and (c) HONO concentrations retrieved from MAX-DOAS in the near the surface (0–100 m), lower (100–300 m), middle (500–700 m), and upper (900–1100 m) layers during this campaign.
Remotesensing 16 03192 g004
Figure 5. The diurnal variations of the vertical distributions of (ad) aerosol extinction, (eh) NO2, and (il) HONO from July to August 2019. The black lines in the (ac), (eg), and (ik) represent the atmospheric boundary layer height.
Figure 5. The diurnal variations of the vertical distributions of (ad) aerosol extinction, (eh) NO2, and (il) HONO from July to August 2019. The black lines in the (ac), (eg), and (ik) represent the atmospheric boundary layer height.
Remotesensing 16 03192 g005
Figure 6. Mean profiles under haze and non-haze days, and growth percentages of (a) aerosol extinction, (b) NO2, and (c) HONO concentrations from MAX-DOAS observation.
Figure 6. Mean profiles under haze and non-haze days, and growth percentages of (a) aerosol extinction, (b) NO2, and (c) HONO concentrations from MAX-DOAS observation.
Remotesensing 16 03192 g006
Figure 7. Mean (ad) 72 h/(eh) 24 h backward trajectories of each trajectory cluster arriving at Wangdu at bottom (50 m), lower (200 m), middle (600 m), and upper (1000 m) altitudes in the summer of 2019 and the percentage of allocation to each cluster.
Figure 7. Mean (ad) 72 h/(eh) 24 h backward trajectories of each trajectory cluster arriving at Wangdu at bottom (50 m), lower (200 m), middle (600 m), and upper (1000 m) altitudes in the summer of 2019 and the percentage of allocation to each cluster.
Remotesensing 16 03192 g007
Figure 8. Spatial distributions of WPSCF values for (ad) aerosol extinction coefficient and (eh) NO2 concentration at the bottom (50 m), lower (200 m), middle (600 m), and upper (1000 m) in the summer of 2019 over Wangdu.
Figure 8. Spatial distributions of WPSCF values for (ad) aerosol extinction coefficient and (eh) NO2 concentration at the bottom (50 m), lower (200 m), middle (600 m), and upper (1000 m) in the summer of 2019 over Wangdu.
Remotesensing 16 03192 g008
Figure 9. (a) Observed NO2-to-HONO conversion rates and (b) aerosol extinction at different height layers and seasons.
Figure 9. (a) Observed NO2-to-HONO conversion rates and (b) aerosol extinction at different height layers and seasons.
Remotesensing 16 03192 g009
Figure 10. (ac) The photolysis frequency of HONO ( J H O N O ) and (df) the OH produced from HONO ( P O H H O N O ) at different seasons.
Figure 10. (ac) The photolysis frequency of HONO ( J H O N O ) and (df) the OH produced from HONO ( P O H H O N O ) at different seasons.
Remotesensing 16 03192 g010
Table 1. Summary of DOAS retrieval settings of O4, NO2, and HONO.
Table 1. Summary of DOAS retrieval settings of O4, NO2, and HONO.
ParameterData SourceTrace Gases
O4NO2HONO
Wavelength Range (nm) 338–370338–370335–373
NO2220 K, 298 K, I0-corrected * [54]
O3223 K, 243 K, I0-corrected * [55]
O4293 K [56]
HCHO298 K [57]
H2OHITEMP [58]
BrO223 K [59]
HONO296 K [60]
RingCalculated with QDOAS
Polynomial degree Order 5Order 5Order 5
Intensity offset ConstantConstantConstant
Wavelength calibrationBased on a high resolution solar reference spectrum (SAO2010 solar spectral) [61]
* Solar I0-correction [62].
Table 2. The average and maximum concentrations of NO2, HONO, and aerosol extinction coefficient.
Table 2. The average and maximum concentrations of NO2, HONO, and aerosol extinction coefficient.
PollutantsAverage Extinction Coefficient (km−1) or Concentration (ppb)Maximum Extinction Coefficient (km−1) or Concentration (ppb)
SurfaceLower LayerMiddle LayerUpper LayerSurfaceLower LayerMiddle LayerUpper Layer
Aerosol0.3580.3700.3420.1610.7220.9290.6290.357
NO21.6591.9260.1500.0403.9584.3530.5540.049
HONO0.0790.0420.0210.0120.4900.0970.0230.013
Table 3. Trajectory ratios and averages of aerosol extinction coefficient (corresponding to the 72 h backward trajectories) and NO2 concentration (corresponding to the 24 h backward trajectories) for all trajectory clusters arriving at Wangdu at bottom (50 m), lower (200 m), middle (600 m), and upper (1000 m) in the summer of 2019.
Table 3. Trajectory ratios and averages of aerosol extinction coefficient (corresponding to the 72 h backward trajectories) and NO2 concentration (corresponding to the 24 h backward trajectories) for all trajectory clusters arriving at Wangdu at bottom (50 m), lower (200 m), middle (600 m), and upper (1000 m) in the summer of 2019.
ClusterTraj_RatioAerosol Extinction (km−1)ClusterTraj_RatioNO2 (ppb)
Mean ± S.D.Mean ± S.D.
Bottom (50 m)118.71%0.262 ± 0.045122.89%1.745 ± 0.262
237.74%0.356 ± 0.062218.21%1.898 ± 0.325
37.24%0.429 ± 0.081331.20%1.859 ± 0.258
436.29%0.405 ± 0.068427.70%1.967 ± 0.313
All100%0.358 ± 0.062All100%1.873 ± 0.281
Lower (200 m)148.23%0.419 ± 0.063141.33%2.098 ± 0.304
215.40%0.247 ± 0.040213.76%1.527 ± 0.159
313.08%0.172 ± 0.032319.53%2.278 ± 0.302
423.30%0.501 ± 0.102425.39%2.236 ± 0.426
All100%0.378 ± 0.058All100%2.091 ± 0.307
Middle (600 m)127.11%0.190 ± 0.028130.70%0.117 ± 0.020
228.29%0.367 ± 0.055220.75%0.192 ± 0.035
313.35%0.263 ± 0.042336.29%0.118 ± 0.017
431.24%0.377 ± 0.062412.26%0.143 ± 0.031
All100%0.315 ± 0.045All100%0.137 ± 0.026
Upper (1000 m)120.23%0.108 ± 0.017152.41%0.040 ± 0.006
222.94%0.108 ± 0.016216.44%0.041 ± 0.018
329.51%0.174 ± 0.02139.22%0.042 ± 0.013
427.32%0.129 ± 0.020421.93%0.040 ± 0.021
All100%0.138 ± 0.019All100%0.040 ± 0.015
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

Wei, S.; Hong, Q.; Tan, W.; Chen, J.; Li, T.; Wang, X.; Xue, J.; Fang, J.; Liu, C.; Tanvir, A.; et al. Spatiotemporal Distribution, Sources, and Impact on Atmospheric Oxidation of Reactive Nitrogen Oxides in the North China Plain Agricultural Regions in Summer. Remote Sens. 2024, 16, 3192. https://doi.org/10.3390/rs16173192

AMA Style

Wei S, Hong Q, Tan W, Chen J, Li T, Wang X, Xue J, Fang J, Liu C, Tanvir A, et al. Spatiotemporal Distribution, Sources, and Impact on Atmospheric Oxidation of Reactive Nitrogen Oxides in the North China Plain Agricultural Regions in Summer. Remote Sensing. 2024; 16(17):3192. https://doi.org/10.3390/rs16173192

Chicago/Turabian Style

Wei, Shaocong, Qianqian Hong, Wei Tan, Jian Chen, Tianhao Li, Xiaohan Wang, Jingkai Xue, Jiale Fang, Chao Liu, Aimon Tanvir, and et al. 2024. "Spatiotemporal Distribution, Sources, and Impact on Atmospheric Oxidation of Reactive Nitrogen Oxides in the North China Plain Agricultural Regions in Summer" Remote Sensing 16, no. 17: 3192. https://doi.org/10.3390/rs16173192

APA Style

Wei, S., Hong, Q., Tan, W., Chen, J., Li, T., Wang, X., Xue, J., Fang, J., Liu, C., Tanvir, A., Xing, C., & Liu, C. (2024). Spatiotemporal Distribution, Sources, and Impact on Atmospheric Oxidation of Reactive Nitrogen Oxides in the North China Plain Agricultural Regions in Summer. Remote Sensing, 16(17), 3192. https://doi.org/10.3390/rs16173192

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