Seasonal Variation of Nitrate Concentration and Its Direct Radiative Forcing over East Asia

This study investigated the seasonal variation of nitrate concentration and its radiative forcing over East Asia by using an online-coupled regional climate model. Comparison with a series of in-situ observations from Acid Deposition Monitoring Network in East Asia (EANET) and China demonstrated a good skill of the model in reproducing the magnitude and seasonality of nitrate concentration across East Asia. It was found that nitrate concentration in Beijing and Tianjin exhibited the maximum in summer and the minimum in winter possibly due to stronger chemical oxidation and regional transport effect of larger emissions from the north China Plain in summer, whereas in most areas of East Asia, nitrate concentration was higher in winter and lower in summer, consistent with the seasonality of NOx emission. Surface nitrate concentration was higher over the lower reaches of the Yellow River, followed by the middle to lower reaches of the Yangtze River and portions of south China, and lower in Korean Peninsula and Japan. The annual mean surface nitrate concentration was predicted to be 2.9 μg ̈m ́3 for East Asia and 8.5 μg ̈m ́3 for east China. All-sky direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) exhibited the largest forcing up to ́7 W ̈m ́2 over the lower reaches of the Yellow River, and lower forcing of ~ ́2 W ̈m ́2 in the Korean Peninsula and Japan. Clear-sky DRF by nitrate reached the maximum in spring and the minimum in summer over both East Asia and east China, whereas DRF under all-sky condition showed its maximum in autumn, associated with seasonalities of nitrate column burden, relative humidity, and cloud effect. Annual mean all-sky DRFs at TOA were estimated to be  ́1.7 W ̈m ́2 and  ́3.7 W ̈m ́2 over East Asia and east China, respectively, significantly larger than global annual mean, suggesting the important role of nitrate aerosol in environment and climate change over East Asia.


Introduction
Aerosol pollution in East Asia has become a serious problem in the past decades mainly due to continuous economic growth and ineffective emission control in this region.As a result, more and more attention has been paid to aerosol pollution problems in East Asia.Sulfate, nitrate, black carbon, organic carbon, and mineral dust are major aerosol types in this region.Among these aerosols, many studies have been conducted concerning the origin, transport and transformation, removal, and environmental and climatic effects of sulfate, black carbon, and dust; however, nitrate, one of the major anthropogenic aerosol components, dose not attracted as many concerns.
In recent years, nitrate aerosol increases along with increasing emission of NO x , the gaseous precursor of nitrate, which is mainly from automobile and industry.Nitrate account for about 4%-30% in PM 2.5 /PM 10 mass concentration in urban areas of China in terms of seasonal or annual mean [1][2][3][4][5][6].
Wang et al. [5] reported that nitrate accounted for 25% of PM 2.1 mass concentration, just less than the mass fraction of sulfate (30%), at five sites in the Yangtze River Delta of east China during 2012-2013, demonstrating a comparable contribution of sulfate and nitrate to fine aerosol in east China at present time.Nitrate concentration was consistently increased in the past decade in China.For example, in Beijing, Yang et al. [3] found that although annual sulfate and ammonium concentrations decreased by 5% and 8% during 2000-2008, respectively, nitrate concentration experienced a solid growth by 20%.The increase of nitrate precursor emission led to this change.Xie and Han [7] has analyzed emission inventories from previous studies [8][9][10] and found that the annual emission of NO x in China had persistently increased from 15.4 Tg/year in 2000 to 21.0 Tg/year in 2006 and further to 28.5 Tg/year in 2010, 85% increment between 2010 and 2000.On the contrary, the annual emission of sulfur dioxide (SO 2 ), the gaseous precursor of sulfate, was estimated to increase from 22.8 Tg/year in 2000 to 31.4 Tg/year in 2006 but decrease to 28.4 Tg/year in 2010.Meanwhile, the overall increment of anthropogenic PM 2.5 emission and black carbon emission was 30% and 20%, respectively, between 2010 and 2000, while anthropogenic PM 10 emission decreased sharply by 74%.It was clear that NO x emission had the most rapid increase among pollutants relating with aerosols.This implies that nitrate aerosol plays a more and more important role in China in recent years.
Nitrate is formed through a series of multi-phase chemical transformation including thermodynamic equilibrium and in-cloud aqueous oxidation.NO x is first oxidized by OH radical to form gaseous nitric acid (HNO 3 ), and then HNO 3 is neutralized by NH 3 via thermodynamic equilibrium or/and in-cloud processes, forming nitrate aerosol.Heterogeneous reaction is another important formation pathway of nitrate.Gaseous species, such as HNO 3 , diffuses onto surface of preexisting aerosol (e.g., dust), and then reacts with alkaline components of the aerosol, finally forming nitrate.
Optical properties of nitrate, such as refractive index, are similar to those of sulfate [11,12], so the increase of nitrate concentration can also lead to increase of extinction coefficient, resulting in visibility reduction and a negative aerosol direct radiative forcing (DRF) at the surface and the top of the atmosphere.As the increase of NO x emission in recent years, the strength of DRF exerted by nitrate aerosol increases accordingly [12][13][14][15].Based on observations of aerosol optical properties (such as aerosol optical depth), total aerosol DRF can be estimated (e.g., [16,17]); however, the observation-based estimation is just available at limited sites and not able to identify DRF due to each aerosol component.Numerical simulation can overcome this shortcoming and can simulate aerosol DRF in past/present/future times.Using a global atmospheric chemistry-climate model, Liao et al. [12] suggested an increasing trend of DRF due to anthropogenic nitrate and a decreasing trend of DRF due to anthropogenic sulfate, and nitrate DRF was expected to exceed sulfate DRF in the future, implying a potentially more important role of nitrate in radiation transfer and climate.By running a global model, Bauer et al. [13] predicted that the nitrate burden and nitrate DRF in 2030 would be about 1.3 times those in 2000.Li et al. [14] projected an increasing global importance of nitrate DRF in the future especially over East Asia based on a global atmospheric general circulation model.With the Hadley Center global climate model, Bellouin et al. [15] suggested that nitrate will partially replace sulfate to become an important anthropogenic pollutant in the remainder of 21st century in terms of radiative effects.
Although nitrate plays important roles in environment and climate change, concerns paid to nitrate are much less than other aerosol types such as sulfate, black carbon, and mineral dust, especially over East Asia/China where intensive emissions of nitrate precursors (NO x ) exist.Previous studies on nitrate aerosol over East Asia/China were limited although nitrate accounts for a large fraction of anthropogenic aerosols.Most of those studies only investigated the spatial distribution of nitrate concentration in a short time period and model comparison with observation was usually insufficient (e.g., [18][19][20][21]).Thus far, only a few studies were conducted regarding DRF due to nitrate in East Asia/China [22][23][24] in which nitrate was found to be important in estimating aerosol radiative and climatic effects.Given the above limitation, this study aims to investigate seasonal variations of nitrate concentration and its direct radiative forcing over East Asia by using an online-couple regional climate model (RIEMS-Chem: Regional Integrated Environmental Model System with atmospheric chemistry).
This paper is organized as follows: First, RIEMS-Chem and relevant inputs and observational data are described.Then, model comparison with in-situ observations of surface nitrate concentration will be conducted.Finally, the regional distribution and seasonal variation of nitrate concentration and induced DRF are analyzed and the estimation of DRF over East Asia and east China are presented.This study provides valuable information for further understanding the aerosol's roles in radiation and climate change in this region and for global estimation of nitrate aerosol.

Model Description
RIEMS-Chem is an online-coupled regional climate-chemistry model [25].It couples the host model RIEMS with full atmospheric chemistry processes.The host model was developed based on the nonhydrostatic dynamic structure of the fifth-generation NCAR/Penn State Mesoscale Model (MM5) [26].It includes several sub-models and parameterizations to describe a series of important physical processes, such as the Biosphere-Atmosphere Transfer Scheme (BATS) for land surface process [27], the Medium-Range Forecast (MRF) scheme for planetary boundary layer process [28], the Grell cumulus convective parameterization scheme [29], and a modified radiation sub-model based on the radiation package of the NCAR Community Climate Model (CCM3) [30] for radiation transfer calculation.RIEMS has been successfully used to simulate East Asian monsoon climate and to explore the interaction among physical, biological, and chemical processes [31][32][33].RIEMS is also one of the key models participating in the Regional Climate Model Inter-comparison Project (RMIP) for Asia [34].
Chemistries in RIEMS-Chem include CB-IV gas phase mechanism [35], ISORROPIA thermodynamics model [36], aqueous phase chemical reactions [37], and heterogeneous reactions [38,39].Major aerosol components, namely sulfate, nitrate, ammonium, black carbon (BC), organic carbon (OC), other anthropogenic particulate matters (AnthPM), dust, and sea salt are considered in this model.Nitrate and ammonium aerosols are formed via thermodynamic equilibrium process represented by ISORROPIA.Nitrate formation through aqueous phase reaction is considered as well using a treatment similar to that in Chang et al. [37].Nine heterogeneous reactions of gaseous species (SO 2 , NO 2 , HNO 3 , O 3 , N 2 O 5 , NO 2 , NO 3 , OH, and HO 2 ) on dust surface are parameterized and incorporated into RIEMS by Li and Han [39].Based on the OPAC database (Optical Properties of Aerosols and Clouds), which was introduced in detail by Hess et al. [40], a lognormal size distribution is assumed for inorganic aerosols (sulfate, nitrate, ammonium), with median radius of 0.07 µm and geometric standard deviations of 2.0, respectively.Dust deflation and sea salt generation are calculated using the same schemes as those in Han et al. [41], but five size bins (0.1-1.0, 1.0-2.0,2.0-4.0,4.0-8.0,8.0-20.0µm) are used to save computation time.
Dry deposition velocities of gases are parameterized by the scheme of Walmsley and Wesely [42].Dry deposition velocities of aerosols are calculated as the inverse of total resistances plus the gravitational settling term for different aerosol types and sizes.Below-cloud scavenging of aerosols is parameterized based on an expression scavenging rate which is a function of precipitation rate and collision efficiency of particle by hydrometeor [41].
A Mie-theory based parameterization proposed by Ghan and Zaveri [43] is used to calculate aerosol optical properties including extinction coefficient, single scattering albedo, and asymmetry factor.The advantage of this parameterization is that aerosol optical properties are pre-calculated by a Mie code and then fitted by Chebyshev polynomials to create a table of polynomial coefficients.By using this table, the calculation of aerosol optical properties becomes much faster than the traditional Mie solution.The kappa (κ) parameterization [44][45][46] with the form of (1) is adopted to represent aerosol hygroscopic growth.g(RH), RH, and κ denote aerosol diameter hygroscopic growth factor, fractional relative humidity, and aerosol hygroscopicity parameter, respectively.Values of κ for inorganic components (sulfate, nitrate, and ammonium) are set to 0.65 [47].DRF is estimated as the differences in the net radiative fluxes at the top of the atmosphere (TOA) or at the earth's surface.At one model time step, radiative fluxes are calculated twice, with and without aerosols.The model with above treatment performs well for aerosol optical depth compared with AERONET and MODIS measurements in east China in previous studies [48,49].More details on RIEMS-Chem also refer to Han et al. [50].
The domain in this study covers most of East Asia.The model applies a Lambert conformal projection with 60 km horizontal resolution, and 16 sigma layers distribute unevenly in the vertical direction with eight layers in the boundary layer.Monthly varied anthropogenic emissions (SO 2 , NO x , CO, NMVOC, BC, primary OC, etc.) with 0.25 ˝ˆ0.25 ˝resolution for the base year 2010 are from the MIX emission inventory prepared for the Model Inter-Comparison Study for Asia (MICS-Asia) [51] and the biomass burning emissions with 0.5 ˝ˆ0.5 ˝resolution are derived from the IPCC AR5 (the Intergovernmental Panel on Climate Change Fifth Assessment Report) inventory [52].The above emissions are bilinearly interpolated to the model projection.NO x emission (figure not shown) is mainly concentrated in north China and east China, with intensive emissions in megacities such as Beijing, Tianjin, Nanjing, Shanghai, Wuhan, Guangzhou, etc.The NO x emission over East Asia shows the maximum in winter (Table S1).The seasonality of NO x emission over east China (20 ˝N-42 ˝N and 100 ˝E-123 ˝E) resembles that over East Asia but the mean intensity is about 15-16 times larger (Table S1).
The simulation period is from 1 December 2009 to 31 December 2010, with the first month as the model spin up.Six hourly NCEP (National Centers for Environmental Prediction) final reanalysis data [53] with 1.0 ˝ˆ1.0 ˝resolution is used to provide initial and boundary conditions for meteorology fields.Boundary conditions of chemical species are derived from MOZART-4 (Model for Ozone and Related chemical Tracers, version 4 [54] simulations and updated every six hour.

Observations
Monthly mean surface nitrate concentration observations for 2010 are obtained from Acid Deposition Monitoring Network in East Asia (EANET) [55].The EANET is a regional cooperation project aimed at sharing acid deposition data and at reaching a common understanding on the state of acid deposition problems among participating countries.This network began its regular monitoring in January 2001.Until now, 13 countries have participated in this network, namely China, Indonesia, Japan, Malaysia, Mongolia, Philippines, Republic of Korea, Russia, Thailand, Viet Nam, Cambodia, Laos PDR, and Myanmar.Nitrate is measured using a filter pack method at EANET sites.Data quality is assured using a rigid and common quality assurance/quality control (QA/QC) program.Monitoring sites used in this study (Figure 1) are classified into three types: urban, rural, and remote.There are only one urban site (Hongwen) in China, two rural sites (Kanghwa and Imsil) and one remote site (Cheju) in Republic of Korea, and two urban sites and four remote sites in Japan.
Because there is only one EANET site in China providing nitrate measurement, observations of nitrate around the year 2010 in China are also collected from previous publications for model validation. is adopted to represent aerosol hygroscopic growth.g(RH), RH, and κ denote aerosol diameter hygroscopic growth factor, fractional relative humidity, and aerosol hygroscopicity parameter, respectively.Values of κ for inorganic components (sulfate, nitrate, and ammonium) are set to 0.65 [47].DRF is estimated as the differences in the net radiative fluxes at the top of the atmosphere (TOA) or at the earth's surface.At one model time step, radiative fluxes are calculated twice, with and without aerosols.The model with above treatment performs well for aerosol optical depth compared with AERONET and MODIS measurements in east China in previous studies [48,49].More details on RIEMS-Chem also refer to Han et al. [50].
The domain in this study covers most of East Asia.The model applies a Lambert conformal projection with 60 km horizontal resolution, and 16 sigma layers distribute unevenly in the vertical direction with eight layers in the boundary layer.Monthly varied anthropogenic emissions (SO2, NOx, CO, NMVOC, BC, primary OC, etc.) with 0.25° × 0.25° resolution for the base year 2010 are from the MIX emission inventory prepared for the Model Inter-Comparison Study for Asia (MICS-Asia) [51] and the biomass burning emissions with 0.5° × 0.5° resolution are derived from the IPCC AR5 (the Intergovernmental Panel on Climate Change Fifth Assessment Report) inventory [52].The above emissions are bilinearly interpolated to the model projection.NOx emission (figure not shown) is mainly concentrated in north China and east China, with intensive emissions in megacities such as Beijing, Tianjin, Nanjing, Shanghai, Wuhan, Guangzhou, etc.The NOx emission over East Asia shows the maximum in winter (Table S1).The seasonality of NOx emission over east China (20°N-42°N and 100°E -123°E) resembles that over East Asia but the mean intensity is about 15-16 times larger (Table S1).
The simulation period is from 1 December 2009 to 31 December 2010, with the first month as the model spin up.Six hourly NCEP (National Centers for Environmental Prediction) final reanalysis data [53] with 1.0° × 1.0° resolution is used to provide initial and boundary conditions for meteorology fields.Boundary conditions of chemical species are derived from MOZART-4 (Model for Ozone and Related chemical Tracers, version 4 [54] simulations and updated every six hour.

Observations
Monthly mean surface nitrate concentration observations for 2010 are obtained from Acid Deposition Monitoring Network in East Asia (EANET) [55].The EANET is a regional cooperation project aimed at sharing acid deposition data and at reaching a common understanding on the state of acid deposition problems among participating countries.This network began its regular monitoring in January 2001.Until now, 13 countries have participated in this network, namely China, Indonesia, Japan, Malaysia, Mongolia, Philippines, Republic of Korea, Russia, Thailand, Viet Nam, Cambodia, Laos PDR, and Myanmar.Nitrate is measured using a filter pack method at EANET sites.Data quality is assured using a rigid and common quality assurance/quality control (QA/QC) program.Monitoring sites used in this study (Figure .1) are classified into three types: urban, rural, and remote.There are only one urban site (Hongwen) in China, two rural sites (Kanghwa and Imsil) and one remote site (Cheju) in Republic of Korea, and two urban sites and four remote sites in Japan.
Because there is only one EANET site in China providing nitrate measurement, observations of nitrate around the year 2010 in China are also collected from previous publications for model validation.Comparisons of monthly mean surface nitrate concentration between observation and simulation at ten EANET (Acid Deposition Monitoring Network in East Asia) sites (a-j).The siteaverage monthly mean nitrate concentration across all the ten sites is also shown (k).Note that y-axis in each panel is different.Site locations see Figure 2a and Table 1.

Comparison with EANET Data
First, the modeled nitrate concentration was compared with EANET observations.Figure 1 shows the monthly mean nitrate concentrations in 2010 at ten EANET sites and the statistics for model performance are listed in Table 1.The Hongwen site (Figure 1a) located in southeast China showed higher The site-average monthly mean nitrate concentration across all the ten sites is also shown (k).Note that y-axis in each panel is different.Site locations see Figure 2a and Table 1.

Comparison with EANET Data
First, the modeled nitrate concentration was compared with EANET observations.Figure 1 shows the monthly mean nitrate concentrations in 2010 at ten EANET sites and the statistics for model performance are listed in Table 1.The Hongwen site (Figure 1a) located in southeast China showed higher nitrate concentration observed in winter and spring (January to April and December, 13-17 µg¨m ´3) than that in summer and autumn, with annual mean of 11.0 µg¨m ´3.The model reproduced the seasonality and magnitude quite well, with the annual mean of 11.1 µg¨m ´3 and the correlation coefficient (R) of 0.87 for monthly mean values.
The Kanghwa and Imsil sites in Korea (Figure 1b,c) exhibited higher levels in winter (~9 µg¨m ´3) and late autumn (November, 6-10 µg¨m ´3), followed by spring (1-5 µg¨m ´3), and the lowest levels in summer and early autumn (July to September, <1 µg¨m ´3).The model generally reproduced the seasonal variation but tended to predict lower values in winter and higher values in autumn at both sites.The remote site Cheju showed similar variation pattern to the above two sites but with much lower levels (Figure 1d).The observed annual mean nitrate concentrations in Kanghwa, Imsil, and Cheju were 3.0 µg¨m ´3, 3.7 µg¨m ´3, and 0.7 µg¨m ´3, while the model simulations were 3.7 µg¨m ´3, 4.1 µg¨m ´3, and 1.4 µg¨m ´3, with R of 0.63, 0.81, and 0.22, respectively.The relatively large model bias at Cheju site could be associated with the coarse model resolution, which could not well represent the large gradients in terrain and land use.The Cheju site is located on a sharp sea cliff; specific local circulation and sea salt process could affect evolution of nitrate aerosol, but cannot be reflected well in the model.
The monthly variation was not evident in the remote site Rishiri in north Japan (Figure 1e) due to low emissions there and frequent influence by clean air from the sea.The model reproduced the seasonal variation but tended to generally overpredict observations.The observed and simulated nitrate concentrations in Rishiri were 0.8 µg¨m ´3 and 1.1 µg¨m ´3, respectively, with R of 0.5.Differently, the other site Tappi in north Japan displayed apparent higher nitrate levels than in Rishiri (Figure 1f), with maximum reaching 3.4 µg¨m ´3 in May, which was mainly due to local sources because long-range transport from continent was not strong at this time.The model overpredicted observation in June mainly due to an underprediction of precipitation and thus wet scavenging.Oki (Figure 1g) and Banryu (Figure 1h) are close in location and nitrate concentrations showed similar levels and monthly variations, lower in summer than those in the other seasons, which were also well captured by the model except overprediction in summer in Banryu.Nitrate level in Tokyo (Figure 1i) was the highest among the Japanese sites, with the monthly mean concentrations in a range of 2.4-6.6 µg¨m ´3, and the annual mean value of 3.8 µg¨m ´3.The seasonality of nitrate in Tokyo resembled those in Oki and Banryu, showing maximums in winter and minimums in summer.The model reproduced the monthly variation quite well except for the overprediction in June.The annual mean nitrate concentration was accurately predicted (3.8 µg¨m ´3) in Tokyo, with a moderate R of 0.56.Hedo is a remote site located in the downwind ocean area of the East Asian continent (Figure 1j).When northwesterly winds prevailed in winter and spring, aerosols can be transported to Hedo, resulting in an elevation of nitrate concentration (1.2-2.6 µg¨m ´3), whereas in summer when southeasterly winds prevailed, clean air from remote ocean dominated this site, leading to a much lower nitrate level (0.1-0.5 µg¨m ´3).The simulated seasonality was consistent with observation reasonably well, but the model predicted a lower nitrate level in summer.The simulated nitrate concentration in terms of annual mean in Hedo was 1.1 µg¨m ´3 with R of 0.78, agreeing well with observation (1.1 µg¨m ´3).
The monthly mean nitrate concentrations averaged over all the ten EANET sites are also presented to illustrate the typical seasonal variation of surface nitrate concentration and the overall model performance over East Asia (Figure 1k).It is encouraging that the model simulation agreed with observation quite well in terms of both seasonal variation pattern and magnitude, with the maximum of about 4 µg¨m ´3 in winter and the minimum of 1.5 µg¨m ´3 in summer.The annual mean nitrate concentrations for all sites were 2.7 µg¨m ´3 and 2.9 µg¨m ´3 from observation and simulation, respectively, with an overall correlation coefficient of 0.92.Both EANET observation and simulation also demonstrated that nitrate concentration was higher over the continent than in the western Pacific region (Korea, Japan).

Comparison with Observations across China
As there was only one EANET site in China reported nitrate concentration in 2010, observations conducted in China from previous publications around the year 2010 were also collected and used for model validation.Observations were mainly conducted in major cities (Table 2) locating in north China (Shangdianzi, Beijing, Tianjin, Shijiazhuang, and Chengde), east China (Nanjing), central China (Wuhan), south China (Guangzhou), and northwest China (Xi'an).Filter pack methods were applied in the above observations.It was noticed that in north China, nitrate concentrations were apparently higher in urban areas (Beijing, Tianjin, and Shijiazhuang) than those in rural and remote areas (Shangdianzi and Chengde).Shijiazhuang exhibited the highest annual mean concentration (30.4 µg¨m ´3) among all sites mainly due to the intensive emission there (Table S1).Beijing and Tianjin had nitrate levels of about 20 µg¨m ´3, whereas Chengde experienced the lowest nitrate level (5.8 µg¨m ´3).Both Beijing and Tianjin showed the highest nitrate concentrations in summer, followed by autumn and spring, and the lowest one in winter.It was noticed that although NO x emissions were low in Beijing and Tianjin in summer, NO x emission in the north China Plain was very high (in a similar level to the winter maximum) during this season (Table S1).The strong northward (under prevailing southerly winds) transport effect of emissions from the north China Plain could lead to elevated nitrate concentrations in the two cities.Besides, stronger chemical oxidation in summer in this region may also enhance nitrate formation as proposed by [56].The minimum concentration occurred in winter in these cities mainly due to frequent passage of strong westerly winds and weak photochemical activity in the cold season.The model reasonably captured the major features in seasonal variation in Beijing and Tianjin but predicted lower nitrate levels in spring and winter due to uncertainties in prediction of meteorological variables and emission inventories.Differently, the observed maximum nitrate in Shijiazhuang, another megacity in north China, appeared in autumn.The model failed to predict the autumn maximum possibly due to uncertainties in the emission inventory (peak in summer and winter) (Table S1).Chengde (a rural site) and Shangdianzi (a background site) experienced lower nitrate levels than those in Beijing, Tianjin, and Shijiazhuang.It is well known that there is little anthropogenic emission north and northwest of Chengde and Shangdianzi, and in seasons except summer, northerlies and northwesterlies prevailed in that region.At the two sites, nitrate levels were observed higher in spring than in other seasons.Heterogeneous reactions on the surface of mineral dust, which is often transported from northwest China in springtime, could be a reason [39].The model generally reproduced the springtime nitrate levels but predicted higher nitrate levels in summer possibly due to an overprediction of the influence from the north China Plain.
The seasonality of nitrate concentration in Nanjing, Wuhan, Guangzhou, and Xi'an were similar to each other and resembled to those at EANET sites (maximum in winter, followed by spring and autumn, and minimum in summer), which was mainly attributed to the seasonality of NO x emission (Table S1).The model well captured such seasonal variation in the above cities except that it predicted lower values by 30%-50% than observations in winter in Nanjing, Wuhan, and Xi'an.
In summary, the seasonality of nitrate concentration in north China generally exhibited the highest nitrate level in summer and the lowest level in winter, whereas that in the other parts of east China showed the maximum nitrate concentration in winter, then in spring and autumn, and the minimum in summer, same as those at EANET sites.In general, the annual mean nitrate concentrations averaged over all the sites in China were 15.6 µg¨m ´3 and 12.3 µg¨m ´3 from observation and simulation, respectively, and exhibited a decreasing level from north China to south China.
The above comparisons between model simulation and observation from EANET dataset and field experiments in east China demonstrated that RIEMS-Chem was capable in reasonably reproducing the major features in spatial distribution and seasonal variation of nitrate concentration over East Asia.

Nitrate Distribution in East Asia
Distributions of seasonal mean surface nitrate concentration are presented in Figure 2. Nitrate concentration in spring (Figure 2a) was general higher in the areas including south of the north China Plain and Henan province with mean values in a range of 15-25 µg¨m ´3.Hunan province and the Sichuan basin also experienced relatively high nitrate levels (15-20 µg¨m ´3).It was noted that nitrate concentration over the Bohai Sea was as high as ~14 µg¨m ´3 because of the eastward transport from the north China Plain under prevailing westerly.In summer (Figure 2b), nitrate aerosols were distributed in the north China Plain and the south portions of northeast China, with the maximum exceeding 25 µg¨m ´3 in the north China Plain.In other regions, such as the Sichuan basin, south of the Yangtze River, the Korean Peninsula, and Japan, nitrate concentrations became much lower due to larger precipitation and efficient wet scavenging as well as relatively clean marine air mass from the western Pacific under southeasterlies.Nitrate distribution in autumn (Figure 2c) resembled that in spring, but had higher level of 15-25 µg¨m ´3 in most areas of east China.It was noteworthy that the distribution pattern in winter (Figure 2d) was different from those in other seasons: nitrate concentrations were simulated to be higher in southern China and the Sichuan basin (20-25 µg¨m ´3), and decreased with increasing latitude, ranged from 15-20 µg¨m ´3 in areas between the Yangtze and the Yellow Rivers, and further to about 5-10 µg¨m ´3 in the north China Plain.Low temperature both favored formation of nitrate aerosol and reduction of vertical mixing, resulting in relatively higher nitrate level in winter, but the relatively weak photochemical activity and more effective wet scavenging by snow could reduce nitrate concentration at high latitudes.

Nitrate Distribution in East Asia
Distributions of seasonal mean surface nitrate concentration are presented in Figure 2. Nitrate concentration in spring (Figure 2a) was general higher in the areas including south of the north China Plain and Henan province with mean values in a range of 15-25 μg•m −3 .Hunan province and the Sichuan basin also experienced relatively high nitrate levels (15-20 μg•m −3 ).It was noted that nitrate concentration over the Bohai Sea was as high as ~14 μg•m −3 because of the eastward transport from the north China Plain under prevailing westerly.In summer (Figure 2b), nitrate aerosols were distributed in the north China Plain and the south portions of northeast China, with the maximum exceeding 25 μg•m −3 in the north China Plain.In other regions, such as the Sichuan basin, south of the Yangtze River, the Korean Peninsula, and Japan, nitrate concentrations became much lower due to larger precipitation and efficient wet scavenging as well as relatively clean marine air mass from the western Pacific under southeasterlies.Nitrate distribution in autumn (Figure 2c) resembled that in spring, but had higher level of 15-25 μg•m −3 in most areas of east China.It was noteworthy that the distribution pattern in winter (Figure 2d) was different from those in other seasons: nitrate concentrations were simulated to be higher in southern China and the Sichuan basin (20-25 μg•m −3 ), and decreased with increasing latitude, ranged from 15-20 μg•m −3 in areas between the Yangtze and the Yellow Rivers, and further to about 5-10 μg•m −3 in the north China Plain.Low temperature both favored formation of nitrate aerosol and reduction of vertical mixing, resulting in relatively higher nitrate level in winter, but the relatively weak photochemical activity and more effective wet scavenging by snow could reduce nitrate concentration at high latitudes.Table 3. Domain averages of annual and seasonal mean surface nitrate concentration, column nitrate concentration, and direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) and at the surface under all-sky and clear-sky conditions for East Asia and east China.

Direct Radiative Forcing Due to Nitrate
As already mentioned above, DRF by nitrate was estimated as the differences in the net radiative fluxes at the top of the atmosphere (TOA) or at the earth's surface with and without nitrate aerosol.3. Seasonality of nitrate concentrations in East Asia and east China both showed the maximum in winter (4.2 µg¨m ´3 and 12.9 µg¨m ´3), followed by spring and autumn, and the minimum in summer (1.6 µg¨m ´3 and 4.5 µg¨m ´3), with annual mean value in east China (8.5 µg¨m ´3) about 3 times that in East Asia (2.9 µg¨m ´3).The nitrate column concentration reached maximum in autumn (9.3 mg¨m ´2 and 22.7 mg¨m ´2) and minimum in summer (6.3 mg¨m ´2 and 12.9 mg¨m ´2), with annual mean burden of 8.0 mg¨m ´2 and 18.6 mg¨m ´2 for East Asia and east China, respectively.

Direct Radiative Forcing Due to Nitrate
As already mentioned above, DRF by nitrate was estimated as the differences in the net radiative fluxes at the top of the atmosphere (TOA) or at the earth's surface with and without nitrate aerosol.

The Distribution of Seasonal Mean DRF Due to Nitrate
Figure 4 shows the spatial distributions of seasonal mean DRF due to nitrate at TOA under all-sky condition.Nitrate exerted a consistent negative forcing in the studying domain in all seasons both at TOA and at the surface (not shown) due to its scattering property.In spring, DRF were mainly distributed over the lower reaches of the Yellow River including the north China Plain, the Bohai Sea, and the East China Sea, with the maximum approaching ´6 W¨m ´2 (Figure 4a).The distribution of DRF in summer was similar to that in spring, with the maximum negative forcing reaching ´12 W¨m ´2 in north China, and DRF over the East China Sea was weakened due to change in wind direction (Figure 4b).In autumn (Figure 4c), DRF in the lower reaches of the Yellow River reached ´7 W¨m ´2, and DRF in southern China and the Sichuan basin (´4-´6 W¨m ´2) were stronger than those in summer (0-´4 W¨m ´2).DRF in winter exhibited the maximum as high as ´10 W¨m ´2 in the Sichuan basin, and ´5-´7 W¨m ´2 in central China between the Yellow River and the Yangtze River (Figure 4d).It was noticed that all-sky DRF distribution (Figure 4) did not follow that of surface concentration (Figure 2) exactly because DRF was also associated with distribution of nitrate column concentration and affected by cloud.Figure 4 shows the spatial distributions of seasonal mean DRF due to nitrate at TOA under allsky condition.Nitrate exerted a consistent negative forcing in the studying domain in all seasons both at TOA and at the surface (not shown) due to its scattering property.In spring, DRF were mainly distributed over the lower reaches of the Yellow River including the north China Plain, the Bohai Sea, and the East China Sea, with the maximum approaching −6 W•m −2 (Figure 4a).The distribution of DRF in summer was similar to that in spring, with the maximum negative forcing reaching −12 W•m −2 in north China, and DRF over the East China Sea was weakened due to change in wind direction (Figure 4b).In autumn (Figure 4c 4d).It was noticed that all-sky DRF distribution (Figure 4) did not follow that of surface concentration (Figure 2) exactly because DRF was also associated with distribution of nitrate column concentration and affected by cloud.Figure 5b shows DRF under clear-sky condition, ranging from 0 to −15 W•m −2 , which was apparently stronger than that under all-sky condition.Compared with clear-sky condition (Figure 5b), DRF under all-sky condition (Figure 5a) was obviously reduced due to cloud effect, with maximum reduction occurring in south China where cloud occurrence and amount were significantly higher than in north China.Figure 5 shows the spatial distributions of annual mean DRF due to nitrate at TOA under both all-sky and clear-sky conditions.The all-sky DRF (Figure 5a) exhibited maximum (up to ´7 W¨m ´2) over the lower reaches of the Yellow River, about ´6 W¨m ´2 over the Sichuan basin, and ´2-´6 W¨m ´2 over the rest areas of east China.It is noteworthy that the nitrate DRF in the western Pacific regions was at the same level as those in south China, indicating the strong outflow from the continent.In Korean peninsula and Japan, DRF were around ´2 W¨m ´2, apparently lower than that in east China.Figure 5b shows DRF under clear-sky condition, ranging from 0 to ´15 W¨m ´2, which was apparently stronger than that under all-sky condition.Compared with clear-sky condition (Figure 5b), DRF under all-sky condition (Figure 5a) was obviously reduced due to cloud effect, with maximum reduction occurring in south China where cloud occurrence and amount were significantly higher than in north China.Wang et al. [22] and Zhang et al. [23] simulated the nitrate-induced all-sky DRF at TOA in China by using regional climate models with emission inventories from Streets et al. [8] for the base year 2000.The spatial distribution of the annual mean DRF from this study was generally consistent with their results; however, the magnitude of DRF at TOA was higher in this work, with the maximum value up to −7 W•m −2 over east China, about 2.3 times those in their studies (−3 W•m −2 ).Considering the fact that nitrate precursor (NOx) emission in China increased by about 85% from 2000 to 2010 [7], the larger DRF was expected and the estimation in this study was realistic.

Seasonal Variation of Domain-Averaged DRF Due to Nitrate
The seasonal mean DRF at TOA and at the surface by nitrate aerosol averaged over East Asia and east China under clear-sky and all-sky conditions are presented in Table 3.
Under clear-sky condition, DRF reached its maximum in spring and its minimum in summer over East Asia and east China both at TOA and at the surface, with stronger forcing in east China than that in East Asia.Annual mean DRF was up to −9 W•m −2 at TOA over east China, more than twice that over East Asia (−3.8 W•m −2 ).DRFs at TOA and at the surface are almost the same due to strong scattering property of nitrate.It was noteworthy that the nitrate column concentration reached the maximum in autumn, not consistent with the DRF spring maximum.Model results showed higher relative humidity (figure not shown) in the middle and upper troposphere in spring than that in autumn in East Asia.Given the similar level of column burdens between spring (9.0 mg•m −2 ) and autumn (9.3 mg•m −2 ), the higher RH could lead to stronger aerosol hygroscopic growth and thus larger DRF in spring than in autumn.
Under all-sky condition, the maximum DRF occurred in autumn, which was consistent with the seasonal variation of column burden.Although clear-sky DRF peaked in spring, it was noted that cloud amount in spring was more than that in autumn as well, thus cloud effect, which tended to reduce DRF by reflecting the incoming solar energy, was larger in spring than in autumn, resulting in larger DRF in autumn than in spring under cloudy condition.In terms of annual average, all-sky DRFs due to nitrate at TOA were estimated to be −1.7 W•m −2 and −3.7 W•m −2 over East Asia and east China, respectively.The global annual mean DRF due to nitrate aerosol reported in the 5th IPCC report was about −0.11 W•m −2 [61], considerably smaller than that over East Asia, especially over east China.This indicates the increasing important role of nitrate aerosol in climate change.

Conclusions
The spatial distribution and seasonal variation of nitrate concentration and its radiative forcing in East Asia were investigated by using an online-coupled region climate model RIEMS-Chem.Model validation with a series of in-situ measurements of nitrate concentration from EANET and from China demonstrated that RIEMS-Chem was able to simulate reasonably well the magnitude and seasonal variation of nitrate concentration in this region.Annual mean surface nitrate concentration was higher over the lower reaches of the Yellow River, the middle to lower reaches of the Yangtze River and portions of south China, and was lower in the Korean Peninsula and Japan.Nitrate concentration in Beijing and Tianjin exhibited the maximum in summer and the minimum in winter Wang et al. [22] and Zhang et al. [23] simulated the nitrate-induced all-sky DRF at TOA in China by using regional climate models with emission inventories from Streets et al. [8] for the base year 2000.The spatial distribution of the annual mean DRF from this study was generally consistent with their results; however, the magnitude of DRF at TOA was higher in this work, with the maximum value up to ´7 W¨m ´2 over east China, about 2.3 times those in their studies (´3 W¨m ´2).Considering the fact that nitrate precursor (NO x ) emission in China increased by about 85% from 2000 to 2010 [7], the larger DRF was expected and the estimation in this study was realistic.

Seasonal Variation of Domain-Averaged DRF Due to Nitrate
The seasonal mean DRF at TOA and at the surface by nitrate aerosol averaged over East Asia and east China under clear-sky and all-sky conditions are presented in Table 3.
Under clear-sky condition, DRF reached its maximum in spring and its minimum in summer over East Asia and east China both at TOA and at the surface, with stronger forcing in east China than that in East Asia.Annual mean DRF was up to ´9 W¨m ´2 at TOA over east China, more than twice that over East Asia (´3.8 W¨m ´2).DRFs at TOA and at the surface are almost the same due to strong scattering property of nitrate.It was noteworthy that the nitrate column concentration reached the maximum in autumn, not consistent with the DRF spring maximum.Model results showed higher relative humidity (figure not shown) in the middle and upper troposphere in spring than that in autumn in East Asia.Given the similar level of column burdens between spring (9.0 mg¨m ´2) and autumn (9.3 mg¨m ´2), the higher RH could lead to stronger aerosol hygroscopic growth and thus larger DRF in spring than in autumn.
Under all-sky condition, the maximum DRF occurred in autumn, which was consistent with the seasonal variation of column burden.Although clear-sky DRF peaked in spring, it was noted that cloud amount in spring was more than that in autumn as well, thus cloud effect, which tended to reduce DRF by reflecting the incoming solar energy, was larger in spring than in autumn, resulting in larger DRF in autumn than in spring under cloudy condition.In terms of annual average, all-sky DRFs due to nitrate at TOA were estimated to be ´1.7 W¨m ´2 and ´3.7 W¨m ´2 over East Asia and east China, respectively.The global annual mean DRF due to nitrate aerosol reported in the 5th IPCC report was about ´0.11 W¨m ´2 [61], considerably smaller than that over East Asia, especially over east China.This indicates the increasing important role of nitrate aerosol in climate change.

Conclusions
The spatial distribution and seasonal variation of nitrate concentration and its radiative forcing in East Asia were investigated by using an online-coupled region climate model RIEMS-Chem.Model validation with a series of in-situ measurements of nitrate concentration from EANET and from China demonstrated that RIEMS-Chem was able to simulate reasonably well the magnitude and seasonal variation of nitrate concentration in this region.Annual mean surface nitrate concentration was higher over the lower reaches of the Yellow River, the middle to lower reaches of the Yangtze River and portions of south China, and was lower in the Korean Peninsula and Japan.Nitrate concentration in Beijing and Tianjin exhibited the maximum in summer and the minimum in winter mainly due to stronger chemical oxidation and regional transport effect of larger emissions from the north China Plain in summer.In most areas of East Asia (including other parts of east China, the Korean Peninsula, and Japan), nitrate concentration generally peaked in winter and lower in summer, which was consistent with the seasonality of NO x emission.The domain averages over East Asia and east China both showed the maximum in winter (4.2 µg¨m ´3 and 12.9 µg¨m ´3) and the minimum in summer (1.6 µg¨m ´3 and 4.5 µg¨m ´3), with annual means of 2.9 µg¨m ´3 and 8.5 µg¨m ´3, respectively.The distribution of annual mean nitrate column concentration was similar to that of the surface.The maximum seasonal mean nitrate burden appeared in autumn (9.3 mg¨m ´2 and 22.7 mg¨m ´2) while the minimum in summer (6.3 mg¨m ´2 and 12.9 mg¨m ´2) for both East Asia and east China.
Nitrate exerted negative forcing in all seasons both at TOA and at the surface due to its scattering property.DRF due to nitrate at TOA under all-sky condition generally exhibited stronger negative forcing over the lower reaches of the Yellow River, with the maximum annual mean forcing up to ´7 W¨m ´2 there.DRF level was strong over the Sichuan basin (´6 W¨m ´2), moderate over the rest parts of east China (´2-´6 W¨m ´2), and relatively low in the Korean Peninsula and Japan (~´2 W¨m ´2).Compared with clear-sky condition, all-sky DRF was reduced by cloud effect, with maximum reduction occurring in south China.In terms of domain and seasonal average, DRF by nitrate under clear-sky condition reached the maximum in spring and the minimum in summer over East Asia and east China both at TOA and at the surface, associated with the seasonal variations of both column burden and relative humidity.On the contrary, DRF under all-sky condition showed its maximum in autumn, which was consistent with the seasonality of nitrate column burden.The domain and annual mean all-sky DRFs due to nitrate at TOA were estimated to be ´1.7 W¨m  S1: Monthly anthropogenic NO x emission used in the simulation (unit: Gg/month/grid).Emission data are extracted from the 0.25 ˝ˆ0.25 ˝inventory.Regional mean data are presented for East Asia (the whole domain), east China (20

Figure 1 .
Figure 1.Comparisons of monthly mean surface nitrate concentration between observation and simulation at ten EANET (Acid Deposition Monitoring Network in East Asia) sites (a-j).The siteaverage monthly mean nitrate concentration across all the ten sites is also shown (k).Note that y-axis in each panel is different.Site locations see Figure2aand Table1.

Table 1 .
Statistics of surface nitrate concentrations between EANET observation and simulation.Obs (μg•m −3 ) and Sim (μg•m −3 ) are annual average of monthly nitrate concentrations from observation and simulation, respectively, R is correlation coefficient, MBE (μg•m −3 ) is mean bias error, and NMB (%) is normalized mean bias.

Figure 1 .
Figure 1.Comparisons of monthly mean surface nitrate concentration between observation and simulation at ten EANET (Acid Deposition Monitoring Network in East Asia) sites (a-j).The site-average monthly mean nitrate concentration across all the ten sites is also shown (k).Note that y-axis in each panel is different.Site locations see Figure2aand Table1.

Figure 2 .
Figure 2. Seasonal mean surface nitrate concentrations (μg•m −3 ) from simulation for (a) spring; (b) summer; (c) autumn; and (d) winter in 2010.Observation sites including ten EANET sites (solid blue squares) and eight Chinese sites (open blue triangles) are marked and their locations refer to Tables 1 and 2, respectively.The annual mean distribution of nitrate concentration near the surface (Figure 3a) was characterized by high levels of 15-25 μg•m −3 in the lower reaches of the Yellow River (including south of the north China Plain and Henan province) and in the Sichuan basin.In the middle to lower reaches of the Yangtze River and south China, nitrate concentrations ranged from 5 μg•m −3 to 15 μg•m −3 .Relatively low nitrate levels (5-10 μg•m −3 ) occurred in southeast China and northeast China.The Korean Peninsula and Japan experienced lower nitrate concentration than those in east China, in a range of ~5 μg•m −3 .Figure 3b shows the distribution of annual mean column concentration of nitrate,

Figure 2 .
Figure 2. Seasonal mean surface nitrate concentrations (µg¨m ´3) from simulation for (a) spring; (b) summer; (c) autumn; and (d) winter in 2010.Observation sites including ten EANET sites (solid blue squares) and eight Chinese sites (open blue triangles) are marked and their locations refer to Tables1 and 2, respectively.
which was generally similar to that of the surface concentration, with values of 25-35 mg•m −2 over broad areas of east China, and 10-15 mg•m −2 in the Korean Peninsula and Japan.

Figure 3 .
Figure 3. Annual mean (a) surface concentration (μg•m −3 ) and (b) column burden (mg•m −2 ) of nitrate aerosol from simulation in 2010.The domain averages of annual and seasonal mean nitrate concentration over East Asia and east China (20°N-42°N and 100°E-123°E) are listed in Table 3. Seasonality of nitrate concentrations in East Asia and east China both showed the maximum in winter (4.2 μg•m −3 and 12.9 μg•m −3 ), followed by spring and autumn, and the minimum in summer (1.6 μg•m −3 and 4.5 μg•m −3 ), with annual mean value in east China (8.5 μg•m −3 ) about 3 times that in East Asia (2.9 μg•m −3 ).The nitrate column concentration reached maximum in autumn (9.3 mg•m −2 and 22.7 mg•m −2 ) and minimum in summer (6.3 mg•m −2 and 12.9 mg•m −2 ), with annual mean burden of 8.0 mg•m −2 and 18.6 mg•m −2 for East Asia and east China, respectively.
Figure4shows the spatial distributions of seasonal mean DRF due to nitrate at TOA under allsky condition.Nitrate exerted a consistent negative forcing in the studying domain in all seasons both at TOA and at the surface (not shown) due to its scattering property.In spring, DRF were mainly distributed over the lower reaches of the Yellow River including the north China Plain, the Bohai Sea, and the East China Sea, with the maximum approaching −6 W•m −2 (Figure4a).The distribution of DRF in summer was similar to that in spring, with the maximum negative forcing reaching −12 W•m −2 in north China, and DRF over the East China Sea was weakened due to change in wind direction (Figure4b).In autumn (Figure4c), DRF in the lower reaches of the Yellow River reached −7 W•m −2 , and DRF in southern China and the Sichuan basin (−4-−6 W•m −2 ) were stronger than those in summer (0-−4 W•m −2 ).DRF in winter exhibited the maximum as high as −10 W•m −2 in the Sichuan basin, and −5-−7 W•m −2 in central China between the Yellow River and the Yangtze River (Figure4d).It was noticed that all-sky DRF distribution (Figure4) did not follow that of surface concentration (Figure2) exactly because DRF was also associated with distribution of nitrate column concentration and affected by cloud.

Figure 4 .
Figure 4. Seasonal mean direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) under all-sky condition for (a) spring; (b) summer; (c) autumn; and (d) winter in 2010.

Figure 5
Figure5shows the spatial distributions of annual mean DRF due to nitrate at TOA under both all-sky and clear-sky conditions.The all-sky DRF (Figure5a) exhibited maximum (up to −7 W•m −2 ) over the lower reaches of the Yellow River, about −6 W•m −2 over the Sichuan basin, and −2-−6 W•m −2 over the rest areas of east China.It is noteworthy that the nitrate DRF in the western Pacific regions was at the same level as those in south China, indicating the strong outflow from the continent.In Korean peninsula and Japan, DRF were around −2 W•m −2 , apparently lower than that in east China.Figure5bshows DRF under clear-sky condition, ranging from 0 to −15 W•m −2 , which was apparently stronger than that under all-sky condition.Compared with clear-sky condition (Figure5b), DRF under all-sky condition (Figure5a) was obviously reduced due to cloud effect, with maximum reduction occurring in south China where cloud occurrence and amount were significantly higher than in north China.

Figure 4 .
Figure 4. Seasonal mean direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) under all-sky condition for (a) spring; (b) summer; (c) autumn; and (d) winter in 2010.

Figure 5 .
Figure 5. Annual mean direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) under (a) all-sky and (b) clear-sky conditions in 2010.

Figure 5 .
Figure 5. Annual mean direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) under (a) all-sky and (b) clear-sky conditions in 2010.

´2 and ´3. 7
W¨m ´2 over East Asia and east China, respectively, which were significantly larger than global annual mean, indicating the important role of nitrate aerosol in environment and climate change in this region.Supplementary Materials: The following are available online at www.mdpi.com/2073-4433/7/8/105/s1, Figure S1: Monthly anthropogenic NOx emission.(a) Beijing, Tianjin, and Shijiazhuang; (b) domain averages for north China (35 ˝N-42 ˝N and 100 ˝E-123 ˝E), south China (20 ˝N-35 ˝N and 100 ˝E-123 ˝E), and north China Plain (36 ˝N-42 ˝N and 113 ˝E-120 ˝E); (c) domain averages for East Asia (the whole domain) and east China (20 ˝N-42 ˝N and 100 ˝E-123 ˝E).Grid resolution is 0.25 ˝ˆ0.25 ˝, Table

Table 1 .
Statistics of surface nitrate concentrations between EANET observation and simulation.Obs (µg¨m ´3) and Sim (µg¨m ´3) are annual average of monthly nitrate concentrations from observation and simulation, respectively, R is correlation coefficient, MBE (µg¨m ´3) is mean bias error, and NMB (%) is normalized mean bias.

Table 2 .
Comparison of observed and simulated nitrate concentrations across China (units: µg¨m ´3).PM 2.5 was measured in all sites except Xi'an where PM 1 was collected. *

Table 3 .
Domain averages of annual and seasonal mean surface nitrate concentration, column nitrate concentration, and direct radiative forcing (DRF) due to nitrate at the top of the atmosphere (TOA) and at the surface under all-sky and clear-sky conditions for East Asia and east China.East Asia: the whole domain; East China: 20 ˝N-42 ˝N and 100 ˝E-123 ˝E. * ˝N-42 ˝N and 100 ˝E-123 ˝E), north China (35 ˝N-42 ˝N and 100 ˝E-123 ˝E), south China (20 ˝N-35 ˝N and 100 ˝E-123 ˝E), and north China Plain (36 ˝N-42 ˝N and 113 ˝E-120 ˝E).