Contributions of mobile , stationary and biogenic sources to air pollution in the Amazon rainforest : a numerical study with WRF-Chem model

This paper evaluates the contributions of the emissions from mobile, stationary and biogenic sources on air pollution in the Amazon rainforest by using the Weather Research and Forecasting with Chemistry (WRFChem) model. The analyzed air pollutants were CO, NOx, SO2, O3, PM2.5, PM10 and VOCs. Five scenarios were defined in order to evaluate the emissions by biogenic, mobile and stationary sources, as well as a future scenario 15 to assess the potential air quality impact of doubled anthropogenic emissions. The stationary sources explain the highest concentrations for all air pollutants evaluated, except for CO, for which the mobile sources are predominant. The anthropogenic sources considered resulted an increasing in the spatial peak-temporal average concentrations of pollutants in 3 to 2,780 times in relation to those with only biogenic sources. The future scenario showed an increase in the range of 3 to 62% in average concentrations and 45 to 109% in peak concentrations 20 depending on the pollutant. In addition, the spatial distributions of the scenarios has shown that the air pollution plume from the city of Manaus is predominantly transported west and southwest, and it can reach hundreds of kilometers in length.


Introduction
The impacts of anthropogenic emissions from urban areas on the environment and human health have been covered in several scientific studies through different methodologies (Kampa and Castanas, 2008;Martins et al., 2010).These analyses have shown that the consequences of such emissions go beyond their initial source and can reach even farther during long periods.In general, most recent studies focus mainly on the aftermath of the anthropogenic emissions from the megacities in the developing world (e.g., Zhu et al., 2013).This has happened because such cities have increased their energy demand, which has been high throughout recent decades (Lutz et al., 2001).For instance, energy consumption in China has increased more than 300 million tons of coal equivalent (MtCE), compared to the previously amount of 1349 MtCE.Most of this energy is produced from burning fossil fuels, mainly coal (Crompton and Wu, 2005).As a consequence, the atmosphere in these regions has experienced a great increase not only in particles but also in gases such as carbon monoxide (CO), sulfur dioxide (SO 2 ), nitrogen oxides (NO x ), and other secondary compounds such as ozone (O 3 ).Understanding the predominant sources of each pollutant is key to designing successful regulatory policies to improve air quality and bring benefits to public health.
Published by Copernicus Publications on behalf of the European Geosciences Union.Determining what is the real contribution of the anthropogenic air pollutants to the environment and to human health is not a straightforward task.Its difficulty lies in separately studying the emissions from natural and human sources and their complex combination.Although urban areas can have an identity associated with their economic profile, the transport and interactions of atmospheric pollutants effects hamper the understanding of the role of each emission species in the studied region.In this context, atmospheric models with an explicit treatment of the physical and chemical processes become an indispensable tool to the studies of the urban pollution impact.Several numerical studies can be found with special focus on different aspects of the pollution impact, such as pollution episodes (Jiang et al., 2012;Michael et al., 2013;Kuik et al., 2015), regional and longdistance transport (Tie et al., 2007;Guo et al., 2009;Lin et al., 2010), secondary formation of gases and particles (Yerramile et al., 2010;Jiang et al., 2012;Lowe et al., 2015), and the effects of land use and land cover changes (Capucim et al., 2015;Rafee et al., 2015).However, such studies did not investigate the impact of isolated urban plumes.
The pronounced growth of urban areas in recent decades, mainly due to the accelerated population growth rate and the migratory flow to cities, resulted in a few areas in the world where an isolated urban area would interact with an environment of natural and homogeneous characteristics.The city of Manaus is one of the unique scenarios in the world where an isolated urban area is in the center of a vast tropical forest area, the Amazon rainforest (Martin et al., 2016).Thus, the region is a valuable laboratory, for example for studying the impact of anthropogenic emissions of air pollutants on atmospheric chemical composition.Although there are a number of studies involving the measurements of atmospheric pollutants in the Amazon region (e.g., Kuhn et al., 2010;Trebs et al., 2012;Baars et al., 2012;Artaxo et al., 2013), most of the contribution was published very recently with the results of the GoAmazon2014/5 project (Martin et al., 2016(Martin et al., , 2017;;Alves et al., 2016;Pöhlker et al., 2016;de Sá et al., 2017;Kourtchev et al., 2016;Hu et al., 2016;Bateman et al., 2017;Cecchini et al., 2016;Jardine et al., 2016).In terms of numerical studies, only a few studies focusing the role of the anthropogenic emission sources in urban areas are available (Freitas et al., 2006;Andreae et al., 2012;Beck et al., 2013;Bela et al., 2015).These studies do not address the participation of each type of emission source on air pollutants, which is important for the formulation of regulatory public policies for air quality management.This type of evaluation can only be performed with the use of atmospheric modeling tools, which require the preparation of an inventory for mobile and stationary sources, since no official and public emissions inventory is available for Manaus.
Given the importance of anthropogenic emissions in the Amazon region, this paper reports a numerical study evaluating the impact of the urban pollution plume on pollutant concentrations above the preserved forest region considering individual contributions from the main mobile and stationary sources in Manaus.In addition, a companion study by Medeiros et al. (2017) focuses on how the changing energy matrix for power production affects air quality in the Manaus region.Both studies were conducted using the Weather Research and Forecasting with Chemistry (WRF-Chem) model, where the simulations for diverse scenarios were performed according to the current conditions of the region.The numerical experiments performed in this work addressed the following issues: participation of mobile and stationary anthropogenic sources and biogenic sources in the concentration of trace gases and particles for the city of Manaus and adjacent areas directly influenced by it; preferential direction and the distance of the urban plume impact from Manaus in the Amazon region; impact on air quality due to the likely future increase in anthropogenic emissions from mobile and stationary sources.

Area and period of study
The region of the study comprises the urban area of Manaus and its surroundings, with a total area of 230 560 km 2 (Fig. 1).Manaus is located in northern Brazil, at latitude 3 • 06 07 and longitude 60 • 01 30 with an urban area of 377 km 2 , presented in the center of the simulated grid.Currently, Manaus has an estimated population of 2 million people, representing more than 50 % of the population in the state of Amazonas, with 99.49 % of its population living in urban areas (IBGE, 2014).The simulation period encompasses the dates of 16-18 March 2014, and this period is part of the rainy season in the region (Fisch et al., 1998).Due to the scarce availability of air quality data, the choice of the days for simulations was made according to that availability, which is necessary to evaluate the performance of the model with the observed ground-based data.Four measurement sites of meteorological variables and air pollutants concentrations were available in different locations in the study area, as presented in Fig. 1: T1, located within the city of Manaus, at the National Institute for Amazonian Research (INPA); T3, located in the north of Manacapuru, approximately 100 km from Manaus, at stations of the project Green Ocean Amazon (GoAma-zon2014/5); Colégio Militar and Federal Institute of Amazonas (IFAM), located in the city of Manaus; and a fourth site associated with the Project REMCLAM Network of Climate Change Amazon from the University of the State of Amazonas (UEA).According to the data from stations located in the city of Manaus, there was no record of rainfall within the given period, whereas the Manacapuru site recorded 7.2 mm rainfall.Regarding the predominant wind direction, the ground-based data observed a prevalence of winds from the north, northeast and east.Figure S1 of the Supplement shows a time evolution of streamlines and the potential temperature advection, where these atmospheric conditions can be seen.In terms of wind direction climatology, the prevailing northeasterly/westerly winds blow all year round, 10-60 • azimuth (November to March) and 90-130 • azimuth (May to September) (Araújo et al., 2002).Another important factor when choosing the period of study was the low incidence of biomass burning wildfires in the rainy season.Fire outbreaks in the Amazon region are predominant during the dry season and can therefore affect the concentration of particulate matter and trace gases (Andreae and Merlet, 2001;Martins and Silva Dias, 2009).Since the purpose is only to analyze the impact of anthropogenic emissions attributed to fossil fuels, periods with the presence of regional wildfires would not be ideal.

Model description
In this study, the coupled WRF-Chem (Weather Research and Forecasting with Chemistry; Grell et al., 2005) model, version 3.2, was applied.The model is an online system that simultaneously predicts meteorological and chemical states.WRF/Chem or its previous versions have been applied in a large number of studies worldwide (e.g., Grell et al., 2005;Wang et al., 2009Wang et al., , 2010;;Tuccella et al., 2012;Vara-Vela et al., 2015).The physical options (Table 3) used were defined according to the recent options inserted in the model, as well as those already tested in simulations by other works (e.g., Vara -Vela et al., 2016).Furthermore, some combinations of physical parametrizations were tested and the combination that best represented the meteorological parameters observed (temperature and relative humidity) was used.
For the biogenic emission, the module used was based on descriptions by Guenther et al. (1993Guenther et al. ( , 1994)), Simpson et al. (1995) andSchoenemeyer et al. (2007).This module deals with isoprene, monoterpenes, volatile organic compound (VOC) emissions and emissions of nitrogen from the soil.Biogenic emissions were calculated using the categories of land use available in the model in which emission rates are estimated from the temperature and photosynthetic active radiation, which is the fraction of solar radiation comprised in the range of the visible spectrum available (0.4 to 0.7 µm) for the photosynthesis process.The aerosol parameterization used was based on the Modal Aerosol Dynamics model for Europe (MADE;Ackermann et al., 1998), developed from the Regional Particulate Model (RPM; Binkowski and Shandar, 1995), embedded with the representation of or-ganic aerosol from the Secondary Organic Aerosol Model (SORGAM, developed by Schell et al., 2001).
For the gas-phase chemistry, the chemical mechanism used was the Regional Acid Deposition Model version 2 (RADM2; Chang et al., 1989), originally developed by Stockwell et al. (1990).RADM2 is widely used in atmospheric models to predict concentrations of oxidants and other pollutants and contains 158 chemical reactions, of which 21 are related to photolysis.
The simulation domain was configured with a 3 km horizontal grid spacing, with 190 grid points in x and 136 grid points in y directions, centered in the city of Manaus (3 • 4 12 S, 59 • 59 24 W).In the vertical grid, 35 levels were defined, with the top of the model at 50 hPa, corresponding to about 20 km in height above the mean sea level.Analysis data from the Global Model Data Assimilation System (GDAS), with a horizontal grid spacing of 1 • and 26 vertical levels were used for initial and boundary conditions of the meteorological variables.For chemical variables, the initial and boundary conditions of simulations consist of idealized, northern hemispheric, mid-latitude, clean environmental conditions as described in Liu et al. (1996) and applied in studies such as those conducted by Grell et al. (2005), Wang et al. (2009Wang et al. ( , 2010)), Tuccella et al. (2012) and Vara-Vela et al. (2016).On the other hand, there are studies where the chemical compound profiles used as initial and boundary conditions were based on anthropogenic emission inventories, including CO, NO x , SO 2 , speciated VOC, black carbon (BC), organic carbon (OC), and particulate material (PM), at 0.5 • × 0.5 • (e.g., Zhang et al., 2009) or higher resolutions (e.g., Tie et al., 2010).Initial and boundary conditions for chemical species have also been extracted from the output of global chemical transport models, such as the Model for OZone And Related chemical Tracers (MOZART; Hu et al., 2010).However, it is important to evaluate the global model results before deciding which global model to use in order to provide the chemical initial and boundary conditions, considering that there are many uncertainties among the global results.
Different times to avoid a spin-up effect are suggested in the literature, with values of 12 (Tuccella et al., 2012;Carvalho et al., 2015), 24 (Wang et al., 2009), 36 (Hu et al., 2010), and 48 h (Tie et al., 2010) of simulation as a model spin-up.There are also studies considering a few days as a model spin-up (e.g., Wang et al., 2010).In this study, a 24 h spin-up time was considered.For land use, spatial data from the Moderate-resolution Imaging Spectroradiometer (MODIS) 2005 file were utilized, with 500 m grid spacing (Schneider et al., 2009).

Vehicle emissions
Emissions of all classes of light-duty and heavy-duty vehicles in the study region were taken into consideration for mobile sources.In order to calculate the emissions of air pollutants, information was collected based on the estimate of the number and type of vehicles, emission factors and average vehicle-use intensity.The calculations included individual contributions of five types/fuel combinations of vehicles (Table 2), considering the data estimated from the Brazilian National Department of Traffic (DENATRAN, 2014) for the city of Manaus.The urban area of Manaus contains 83.22 % of the vehicles in the state of Amazonas, corresponding to over 600 000 vehicles in the current fleet.Therefore, the fraction of the vehicle fleet according to type and fuel consumption for the entire study grid was considered based on data from Manaus.The fractions corresponding to light-duty vehicles using gasohol (a mixture of gasoline and ethanol ranging between 20 and 25 % of anhydrous ethanol), ethanol and flex fuel represent 22.81, 2.49 and 42.29 %, respectively.Heavyduty diesel-powered vehicles represent 8.8 %, and motorcycles using gasohol represent 23.61 %.
The emission factors for different vehicle types/fuel assumed in the calculations were based on experiments conducted inside the road traffic tunnels in the city of São Paulo (Martins et al., 2006;Sánchez-Ccoyllo et al., 2009;Brito et al., 2013), which are the only vehicle emission factor measurements available in Brazil.The values were adopted for CO, NO x , SO 2 , and particulate matter (PM) as well as the distribution of VOC emission fractions related to evaporative, liquid and exhaust emissions for vehicles burning gasohol, ethanol and diesel.These values correspond to those that have been recently used in the study by Andrade et al. (2015) and are listed in Table 2. Based on these emission factors, diurnal profiles using hourly variations for emissions of trace gases and PM were used in the model (Martins et al., 2006;Andrade et al., 2015).The distribution of the fractionation of emission of fine particulate matter, as well as the chemical characterization, was obtained based on several studies in the city of São Paulo for the measurements of the concentration mass and number of particulate matter (Ynoue and Andrade, 2004;Miranda and Andrade, 2005;Albuquerque et al., 2011), and are used in the investigations by Andrade et al. (2015) and Vara-Vela et al. (2015).
In order to calculate the vehicle-use intensity, reference estimates from the first Brazilian National Emissions Inventory of Road Motor Vehicles (MMA, 2011), DENATRAN and the Brazilian National Agency of Petroleum (ANP, 2014a) were used.The intensities for light-duty vehicles, heavy-duty vehicles and motorcycles are 21.3, 128.4 and 70.6 km per day, respectively.
The spatial distribution of the number of vehicles in each grid point was based on nighttime lights from

Stationary source emissions
The emission of thermal power plants (TPPs) of different types of fuel burning have been considered for the inventory of stationary sources in the studied region and the refinery located within the urban area of Manaus.It is noteworthy that the contributions of the industries located within the urban area have little significance in the region, due to the production concentrated mainly in transport and communication areas such as electronics, metal mechanical sectors and the production of motorcycles (MANAUS, 2002).In this case, the contributions on anthropogenic emissions occur indirectly by the high consumption of electricity supplied by the TPPs.The emission of the pollutants per type of TPP in each grid point of study has been calculated according to the estimates of emission factors, fuel consumption and power generation, using Eq. ( 1): where E p(i,j ) represents the emission of pollutant P at each grid point (i, j ), in grams per hour (g h −1 ); EF P is the emission factor of pollutant P , in grams per liter (g L −1 ); FC is the fuel consumption, in liters per kilowatt hour (L kWh −1 ); and PG (i,j,k) is the power generation of TPP (k) at each grid point (i, j ) m in kilowatt (kW).
According to the Generation Database (BIG) of the National Electric Energy Agency (ANEEL, 2014), Brazil has 1890 TPPs in operation, with a total installed capacity of about 37.8 GW.Based on the information contained in BIG, a spatial distribution of TPPs has been performed at each grid point of the study area, with 59 diesel, 6 fuel oil and 8 natural gas power plants (Fig. 2) found in the grid, generating a power of 1851 MW, corresponding to 96.76 % of all gen-eration in the state of Amazonas.The other TPPs located in state are very small.
Due to several types of technologies used in the burning fuels of TPPs, the emission factor admitted in the calculation was the intermediate values between the lower and upper limits adopted by the US Environmental Protection Agency (EPA, 1998(EPA, , 2010)), which has a complete estimation of emission factor for all air pollutants simulated.The emission factors are described in Table 3, and for comparative purposes, the values adopted in the inventory of São Paulo Environmental Protection Agency (CETESB, 2009) have also been listed.The distribution of PM has been done based on the fractionation designed for vehicle emissions.In addition, the speciation of VOC emissions have been performed for each TPP fuel type estimated by EPA (1998EPA ( , 2010)).Regarding the average fuel consumption for each type of TPP, the average value of fuel to produce 1 kWh of power has been adopted as described in the annual report of isolated operation systems for the North region (ELETRO-BRAS, 2013).For fuel oil, diesel and natural gas, the values of 0.27 L kWh −1 , 0.29 kg kWh −1 and 0.28 m 3 kWh −1 , respectively, have been considered.Some approximations that were assumed to obtain the emission factor values of Table 3 were as follows: -Fuel with 1 % sulfur content was assumed to calculate the emission factor of SO 2 .
-For NO x emissions by fuel oil power plants, the minimum value was attributed to the lowest value available for vapor generation above 50 t h −1 and the maximum value was assigned to the greatest value available for plants below 50 t h −1 .
-For the emission factors attributed to SO 2 and PM, the minimum and maximum values were defined based on the grade of fuel burned.In this case, the combustion of lighter distillate oils results in significantly lower PM formation than those from the combustion of heavier residual oils.According to the values proposed by EPA, where E p(i,j ) represents the emission of the pollutant P , in grams per hour (g h −1 ), EP p is the emission factor of pollutant P in grams per liter (g L −1 ) and V is the volume of refined petroleum in liters per hour (L h −1 ).
According to the National Petroleum Agency (ANP, 2014b), REMAN has as average production of approximately 1.7 × 10 6 L of petroleum per hour.The emission factors of the pollutants have been admitted from Presidente Bernardes Refinary (RPBC), located in Cubatão, São Paulo, Brazil.The VOC speciation profiles used have also been obtained from the EPA (1990).
The total emissions for CO, NO x , SO 2, PM and VOCs for the stationary and mobile sources that represent the C0 scenario are shown in Table S1 in the Supplement.These results demonstrate that global inventories underestimate mobile and stationary sources in Manaus.For instance, in the MAC-City anthropogenic emissions inventory (MACCity, 2014), the approximate sums of emissions in the study grid were 0.011 Tg yr −1 (versus 0.068 in this study), 0.003 Tg yr −1 (versus 0.087 in this study), 0.002 Tg yr −1 (versus 0.073 in this study) for CO, NO x and SO 2 , respectively.For example, Fig. 3 shows the contribution, in percent, of NO x and SO 2 emissions from all anthropogenic sectors considered.

Numerical scenarios
The investigations were defined by five numerical experiments (Table 4) to the study of the influence of mobile and stationary emissions in the region, which are as follows: -Scenario C0: considering the main sources in the region (biogenic natural emission, vehicle, TPPs and REMAN refinery) representing the emission inventory of current conditions of the region.
-Scenario C1: numerical experiment with only biogenic natural emission, which simulates an atmospheric condition of an environment without the interference of human activities.
-Scenario C2: simulation with natural and vehicular emissions, intended to evaluate how the anthropogenic emissions characterized on solely vehicular sources affect the atmospheric chemistry.
-Scenario C3: considering natural and TPPs and refinery anthropogenic emissions, aimed at assessing the impact of stationary sources on air pollutant emissions.-Scenario C4: simulation including natural and twice the mobile and stationary emission sources, and twice the urban area of Manaus, which aimed at simulating an environment presenting a growth of the regional population and urban area, as well as increased anthropogenic emissions.The urban expansion has been concentrated in the agricultural and forest areas, preserving the water resources, as well as forest reserves and parks around the city of Manaus.
The meteorological fields and air pollutants obtained from the simulation of scenario C0 have been compared to the ground-based data observed (see Fig. 1).The analysis of meteorological variables has been conducted based on Pielke's skill score (S pielke ) (Pielke, 2002;Hallak and Perreira Filho, 2011), Pearson's correlation coefficient (r), and mean bias (MB).For the air pollutants, in addition to r and MB, two statistical indexes were used, the mean normalized bias error (MNBE) and the mean normalized gross error (MNGE).Such parameters have often been used in several studies to assess the performance of atmospheric models (e.g., Tie et al., 2007;Han et al., 2009;Tuccella et al., 2012).Table 5 shows a summary of the statistical indexes used to evaluate the model performance.

Meteorology
Figure 4 shows the comparison between the observed and the simulated values for temperature and relative humidity at Colégio Militar, IFAM, T1 and T3, corresponding to the periods of 17-18 March 2014.Overall, there is a good representation of the temporal evolution of the temperature (e.g., r = 0.91 for T1) both for the average value and the daily minimum.The T1 station has the highest peak of temperature: about 36 • C, which is the highest among the monitoring points compared to simulations.However, the model has proven to present difficulties in simulating the maximum temperature, mainly for T1 station, which is located in the central part of the urban area of Manaus.For the IFAM station, which is also a central site, diurnal peak (e.g., 33.5 Pielke skill Regarding the statistical indexes presented in Table 6, the correlation coefficient (r) has provided satisfactory results for all stations, with the lowest values for temperature (0.87) and humidity (0.71) associated with T3 station and the highest values of 0.91 and 0.89 for Colégio Militar, respectively.According to Pielke's parameter skill defined in Hallak and Pereira Filho (2011), a good model performance occurs for index values of less than 2. In this case, the performance of the model is satisfactory for most of the stations, except for temperature at T1 (S pielke = 2.9) and relative humidity at IFAM (S pielke = 3.7) and T3 (S pielke = 3.9).According to the mean bias, it has been observed that the simulation, in general, underestimates the majority of observed values, mainly for temperature (T1, MB = −1.5)and relative humidity (IFAM, and T3).

Air pollutants
T1 is the only station that has air quality data and only for NO x , CO and PM 2.5 concentrations.Figure 5 and Table 7 show the comparisons between the simulation (scenario C0) and the observation of such pollutants.In terms of the PM 2.5 , it has been observed, in general, that the model tends to over-estimate the values observed (MNBE = 50.8),whereas with CO and NO x the tendency is the opposite.They underestimated the values in most of the simulation period (MNBE equal a −29.5 and −35.4 for NO x and CO, respectively).The most significant differences have been observed during the peak hours and mainly for CO.The weakest performance of the model happens mainly during the nocturnal peak, situated generally between 18:00 and midnight (local time).A possible explanation for the poor performance of the model during nocturnal peaks of CO is the fact that the fire outbreaks were not considered (four spots of fire occurred during the period inside the grid, although a low incidence of fire outbreaks was the criterion used for chosen the simulation period, Fig. S2), which could have influenced the simulated values.The MB indicates that the three analyzed pollutants present significant differences, with the value of PM 2.5 MB of 1.30 µg m −3 compared to a 0.66 µg m −3 average observed, NO x with values of −26.2 ppb compared to an 88.7 ppb average observed, and CO presenting −135 ppb of MB compared to a 382.6 ppb average observed.Similar to the Pearson coefficient, the model shows good correlations for PM 2.5 , NO x and CO with 0.72, 0.53 and 0.53, respectively.In addition, in order to assess the sensitivity of the model to respond to different emission scenarios, numerical experiments were  as expected.It is not possible to say that the model is overpredicting the actual value, since it is associated with different periods.In addition, the simulated concentrations are much lower than the CO concentrations observed in previous measurements for the dry season in Amazonas, under strong effects of biomass burning emission (e.g., Crutzen et al., 1985).O 3 mixing ratios reported by Kuhn et al. (2010) varied from 21 ppb in the adjacent background to a peak value of 63 ppb for O 3 within the PBL at 100 km distance from Manaus.In the simulated scenario C0, the peak values were in the range 26-29 ppb within the PBL and 40 ppb above it.O 3 mixing ratios downwind of Manaus under the influence of anthropogenic pollution were also reported by Trebs et al. (2012) and were on average 31 ± 14 ppb, with peak values of 60 ppb at a distance of 19 km downwind of Manaus.Such results are in good agreement with the results in this work: 10-43 ppb for the adjacent levels to 400 m at 10 km west of Manaus, considering the current conditions given by scenario C0.In terms of previous numerical studies in the Manaus urban-influenced area, Kuhn et al. (2010) applied a singlecolumn chemistry and meteorological model (Ganzeveld et al., 2002).Two numerical experiments were performed to assess the westward-moving plume, the first one performed by using anthropogenic emission fluxes from the EDGAR 3.2 emission database, but with the pollutant flux based on the 1 • × 1 • grid point of the location of the city of São Paulo, increased by a factor of 7, due to the absence of a local emission inventory for Manaus.In this scenario, the simulated CO, NO x and O 3 mixing ratios at approximately 400 m altitude and 10 km downwind were approximately 140, 4, and 50 ppb, respectively.The impact of stationary sources was evaluated in a second scenario by the inclusion of four ther-mal power plants, adding a total of 1.5 kg NO x day −1 .As a result, the simulated NO x and O 3 mixing ratios changed to 30 and 35 ppb, respectively.The significant effects of thermal power plants are also evidenced in the present study by the low concentration of NO x in the scenarios C1 and C2 (which do not include the stationary emission sources), compared to C0 (current conditions).

Scenario emissions
The discussions on the simulation results related to scenario emissions were analyzed in the lowest model level and divided in three topics -that is, the analysis of the mean and peak behavior of pollutants, spatial analysis, and temporal evaluation described below.

Mean and peak behavior of pollutants
Identifying the differences between scenarios in atmosphere models with elevated complexity is not a straightforward task, due to the large number of components involved, with each of them partially contributing to the evolution of the concentrations in both time and space.In this sense, their impacts on peak average values as in the spatial average should be evaluated.For this purpose, the methodology for microphysical variable proposed by Martins et al. (2009) has been adapted to evaluate the air pollutants in this study.There are two average properties.The first is the spatial averagetemporal average (SATA), which represents a mean value average both spatially and over time, and the second one is the spatial peak-temporal average (SPTA), which corresponds to the spatial peak average over time.The calculation of the values was based on 87 % of the study grid, removing five points of the grid in each extremity, which were the values recommended by Skamarock et al. (2008) to reduce the effects of lateral boundary conditions.It is important to emphasize that, due to the large number of chemicals involved in the VOC, the evaluation has been conducted by summing all output model compounds.
The impact of the different scenarios evaluated according to the spatial average (SATA) and peak (SPTA) concentrations is shown in Table 8.As expected, the lower concentrations of SATA and SPTA of all pollutants are those obtained from natural conditions (scenario C1).Similarly, the future scenario (C4) has shown the highest concentrations simulated for the parameters considered.Considering the current conditions of the region (C0), scenario C4 represents an increase of approximately 35,3,62,4,16 42 and 41 % (SATA) and 45,63,88,45,109,60 and 56 % (SPTA) for the pollutants NO x , CO, SO 2 , O 3 , VOCs, PM 2.5 and PM 10 , respectively.In terms of the individual contribution by sources, it was observed that the emission contributions from stationary sources, predominantly TPPs (scenario C3), were greater than the mobile sources (C2) for all analyzed pollutants, except for CO.For instance, the average concentration is more than double for NO x and about 1 order of magnitude greater for peak values of PM 10 .This indicates that the TPPs are mainly responsible for the high concentrations of most chemical species in the grid.The influence of matrix change is dis-cussed in the companion paper of Medeiros et al. (2017), and more details can found in this paper.

Spatial analysis
In order to analyze the spatial extent and location of impact of the plume of Manaus, the spatial distributions of pollutants have been evaluated.The scenario distributions have indicated that the pollution plume from Manaus could have a great impact on the surrounding area, with predominance to the west and southwest directions.The impact is both on the spatial average and spatial peak values.It has been observed that the high values of emission rates from TPPs significantly contribute to the increase in the air pollution plume area of Manaus during its spread.The influence is predominant in adjacent regions, but it can be extended over more than a hundred kilometers to the west and southwest of the city, to areas that were dominated by the original forest.As an example, Fig. 6 illustrates the spatial distributions of PM 2.5 concentration for 22:00 local time, on 18 March 2014.
For PM 2.5 , considering the contour lines greater than or equal to 5 µg m −3 , which is the value of 10 % of the SPTA (peak value) of the future scenario, it can be observed that the largest fractions of plume covering the area as shown in Fig. 6 were 3024 km 2 (C0), 1386 km 2 (C3), and 6102 km 2  (C4).Based on the plume covering the area and on the contour lines established, for the future scenario, an area 16-fold greater than the Manaus urban area (MUA) has been observed, which is approximately 377 km 2 .In addition, there were no concentrations on the simulation by vehicular sources (C2) within the established contour.The values of the influence in an area and their contour lines for other pollutants are summarized in Table 9 as a function of MUA.

Temporal evaluation
Considering the importance of assessing the temporal development of the concentrations for a locality, two distinct points have been selected.The first one was located within the city of Manaus near the T1 station, and the second in a predominant direction of the plume propagation at approximately 84 km southwest (T3), which can be identified in Fig. 1.It is important to emphasize that the points were chosen in accordance with the highest concentrations of special distribution observed.From this evaluation, the important role of the stationary sources in the results of the concentrations of air pollutants has become evident.For example, based on Fig. 7, the highest NO x concentrations are found in scenario C4, followed by C0 and C3, following the trends presented by SATA and SPTA.Even in locations that are distant from the city of Manaus, high concentrations when compared to concentrations in conditions not influenced by human activities (scenario C1) have been covered.In addi- tion, it has been observed that all scenarios with the participation of TPPs have a significant impact on the direction of the plume spread, while scenarios C1 and C2 remained with their concentrations near zero.This observation reinforces the fact that the purely vehicular source plume could not produce significant impact over long distances.

Discussions and conclusions
The main purpose of this study was to represent the daily cycle of anthropogenic and biogenic emissions that are not available in global inventories which promote the development of air pollution plumes in the city of Manaus.In this sense, to support new design of regulations, it is important to investigate the relative contribution of the mobile and stationary sources as well the area impacted by urban pollu- tion plume.However, the conclusions presented here should be viewed with caution since the simulations involve only two specific days and they are not representative of the entire year.
According to the statistical criteria used for the type of model considered, comparing the observed data with the simulation results to the current state of emissions defined by the baseline scenario (C0), we concluded that the WRF-Chem performance can be considered satisfactory due to the representative criteria for most meteorological fields and air pollutants.
Based on the analyses of observed average and peak concentrations of pollutants, as well as the spatial and temporal distribution of numerical studies, it is clear that stationary sources have an important role in the contribution of human activity in Manaus.The exception is given for carbon monoxide, which has shown few significant contributions.For a comparative evaluation of the results presented herein and those in the literature, it has been observed that there is no parallelism with studies of other parts of the world that have an urban area in a vast tropical forest cover.Several studies have been identified in which the contribution of thermal power plants is higher than the vehicle sources, presenting preserved natural environment conditions, as occurs in Manaus.Studies such as Reddy and Venkataramn (2002) from India show that the main emission of a particulate matter and sulfur dioxide originates from fossil fuel power plants.In this example, SO 2 emissions by power plants are responsible for over 60 % of the country's emissions.In another example, in a study performed in China, Zhang et al. (2012) conclude that most PM 2.5 (nitrate and sulfate) sources are from the energy sector, mainly power plants, exceeding the combined contribution of industrial, residential and transport sources.Also from China, results by Huang et al. (2011) estimate that stationary sources contribute in approximately 97, 86, 89, 91 and 69 % of the total emissions of SO 2 , NO x , PM 10 , PM 2.5 and VOC, respectively.
In the specific case of Manaus, it is important to emphasize that the only hydroelectric plant located in the region that is able to bring electricity to Manaus is the Balbina hydroelectric power station.It has an average capacity of only 250 MW and it contributes to less than 15 % of the electricity used in Manaus (ANEEL, 2014).Therefore, the largest source of electricity for the region comes from power plants burning petroleum fuels, particularly fuel oil and diesel, with higher emission factors if compared, for example, to natural gas.
Another relevant issue is the high electricity consumption of Manaus.According to data from the Brazilian Energy Research Company (EPE) (www.epe.gov.br/), the electricity consumption per capita in Brazil is 2.5 MWh per year, whereas in the city of Manaus the consumption is 7.2 MWh, about 3 times higher.Two aspects make the region highly dependent on electricity.The first is the incentive policies for the consolidation of the industrial area in Manaus in order to foster the economic development of the region.This resulted in a region with a variety of industrial centers, with a high demand for electricity.The second aspect is the fact that the city is located in a humid equatorial climate region, dominated by high temperatures and humidity, as well as little ventilation and torrential rains throughout the year.Such environmental characteristics induce high electricity consumption by both residential and commercial sectors, resulting in the greater burning of fossil fuel power plants that intensify the concentrations of air pollutants.
The spatial distributions of these scenarios indicated that the pollution plume from Manaus could have a great impact on the surrounding areas, mainly in west and southwest directions, reaching hundreds of kilometers from the city.Since most thermal power plants and the REMAN refinery are lo-cated near the banks of the Negro and Solimões rivers, the plume transportation could be influenced by the circulation of river breezes, which defines the trajectory of pollutants.Although the breeze effect is not focused on this paper, this evaluation has been confirmed by other studies for the large rivers in the Amazon, for instance Santos et al. (2014).
Finally, in order to evaluate the potential urban growth in that region, a future scenario has been designed.In summary, this scenario has shown the impact of a possible increase in mobile and stationary emissions in the study region, including the expansion of the urban area.The approach of future scenarios has been studied using several methodologies, mainly in Asia (Zhou et al., 2003;Ohara et al., 2007).In all cases, the increase in air pollution concentrations (e.g., 2020 scenario showed an increase of 12 % in CO concentration) could be observed if the current conditions of the energy matrix were maintained.However, there is the possibility of reductions and a decrease in emission rates by using sustainable energies.
S. A. Abou Rafee et al.: Contributions of mobile, stationary and biogenic sources to air pollution

Figure 1 .
Figure 1.Geographic location of the study area, meteorology stations (Colégio Militar, IFAM and T3) and air quality station (T1).The white contour line shows the delimitation of the city of Manaus.

Figure 2 .
Figure 2. Spatial distribution of diesel (yellow), fuel oil (red) and natural gas (blue) TPPs on the study grid, and the border of Manaus (white).

Figure 3 .
Figure 3. NO x and SO 2 emission contribution from all anthropogenic sectors considered.

×
100 * p k and o k are the predicted and observed value at time k, respectively.among the results, with MB values equal to −8.3 and −10.1, respectively.

Figure 4 .
Figure 4. Temporal evaluation of temperature and relative humidity simulated and observed for Colégio Militar, IFAM, T1 and T3.

Figure 5 .
Figure 5. Temporal evaluation of simulated and observed values for the concentration of pollutants at T1: (a) PM 2.5 , (b) NO x and (c) CO.

Figure 6 .
Figure 6.Spatial distribution of the scenarios studied for PM 2.5 concentration, calculated at 22:00, local time, on 18 March 2014.The black contours illustrate the delimitations of the municipalities on the grid, and the red one represents Manaus.

Figure 7 .
Figure 7. Temporal evaluation of NO x .Panel (a) represents the concentration within Manaus, while panel (b) refers to the point towards the pollution plume about 84 km southwest.

Table 1 .
WRF-Chem physics parametrizations used for this study.

Table 2 .
Emission factors for CO, NO x , PM, SO 2 , and VOCs, in grams per kilometer, for different vehicle type-fuel combinations.

Table 4 .
Summary of simulation scenarios.

Table 5 .
Statistical indexes used to evaluate the model performance.

Table 6 .
Evaluation indexes for model performance of meteorological variables during the simulation period (baseline scenario C0).

Table 7 .
Evaluation indexes for model performance of air pollutants during the simulation period (baseline scenario C0).

Table 8 .
Spatial average-temporal average (SATA) and spatial peak-temporal average (SPTA) of the concentration of trace gases and particles, in the various scenarios * .Baseline scenario (emissions by biogenic, mobile and stationary sources); C1, only natural emissions; C2, scenario referring to natural and mobile emissions; C3, scenario representing natural and stationary emissions; C4, future scenario that includes natural emissions and doubled emissions from mobile and stationary sources.

Table 9 .
Area (km 2 ) influenced by the pollution plume, represented as a function of Manaus urban area, for all simulated scenarios.