Measurement and Modeling Air Quality Impacts of Dust Emissions from Unpaved Roads in Tuxtla Gutierrez, Chiapas

Dust emissions from unpaved roads are one of the main pollutants affecting air quality around the world. As part of initial air quality studies in Tuxtla Gutiérrez (TGZ), Chiapas, Mexico, urban aeolian emission events from unpaved roads and simple meteorological inputs were measured in February 2014 at two different sites located within the city to characterize emissions for representative road conditions and to produce Industrial Source Complex (ISC3) model inputs. Emissions of particulate matter of aerodynamic diameter less than 10 μm (PM10) were determined for eight wind erosion events. PM10 concentrations were measured downwind from sites using a Minivol sampler during February and March 2014. Three high PM10 concentration scenarios, associated with unstable conditions generated by cold fronts (CF) were selected to simulate dust plume dispersion to identify impacted areas. Results show that unpaved roads represent a potential source of dust that affect air quality of urban regions; in this study generating emissions ≥ 1.92 × 10−3 g·m−2·s−1 when winds ≥6 m·s−1 were present. Air pollution events that exceed the Mexico national standard for 24-h average PM10 concentration (≥75 μg·m−3) were observed, impacting different areas in the city, representing a risk to human health. This demonstrates the influence of CF over southern Mexico, generating high PM10 concentrations in urban regions.


Introduction
Mineral dust emissions are one of the main pollution sources that affect air quality around the world [1] and also have important implications on human health and visibility [2][3][4][5][6][7]. Inhalable particulates have been found to cause adverse effects on respiratory health and contribute to excess mortality [8][9][10]. Dust emissions are generated when vehicles transit over unprotected soils [11][12][13] causing air quality impacts affecting human health. Jazcilevich et al. [14] developed a stochastic model to estimate human exposure to emissions from the resuspension of road dust due to moving vehicles. They found that children experience a high risk when they are exposed to this activity. Other studies have shown the contribution of bare soils to air pollution where wind erosion is the main natural phenomenon affecting the emissions of dust particles. Vega et al. [15] have shown that the dry Lake of Texcoco is still a dust source that affects air quality in the north-east of Mexico City.

Study Site
TGZ (16°45′ N, 93°06′ W) is located in the south-east of Mexico and in the central region of Chiapas comprising an area of 334.61 km 2 at an elevation of 522 m above sea level (masl) (see Figure 1). This is a tropical zone with mean annual precipitation oscillating from 900 to 1200 mm·year −1 [34] and surrounding areas consisting of forests. From November to the middle of March, the region is influenced by CF entrances over southern Mexico, which generate high wind intensities (up to 10 m·s −1 ) coming from the west and north-west. The presence of two natural barriers, Mactumatzá Hill at the south and geological formations of the Cañón del Sumidero National Park at the north, both with elevations above 1000 masl, influence the prevailing wind direction and subsequent dispersion of air pollutants affecting the TGZ air quality. Agricultural areas outside of the city have been subjected to strong social pressure due to urban expansion generating unpaved street layouts for new human settlements. When soil aggregates and crusts are crushed by vehicle traffic, the content of loose soil material increases, causing an increase in soil susceptibility to wind erosion [11,35]. The presence of high wind speeds yields fugitive dust emissions impacting the region. . The blue lines represent the presence of rivers and streams. The brown slight lines represent elevation in meters above sea level.

Field Campaign
In order to quantify the dust emissions, a tower with five Big Spring Number Eight (BSNE) sediment catchers [36] was installed in 2014 for approximately four hours (i.e., 14:00 to 18:00 h) each day over unpaved roads during 12 February in the north at Lajas Maciel and during 21 February in the west at Plan de Ayala areas of TGZ (see Figures 1 and 2, Table 1). Although unpaved roads in the city are common, these sites were selected considering their orientation relative to prevailing wind direction and soil surface conditions of roads that are representative of places located around the city. Road emission measurements were only taken when the observed winds were within approximately ±10 degrees to the road direction. Table 2 lists some soil parameters obtained from soil sample analysis as well as some site characteristics. The roads in this study were considered source unlimited for fine loose particles given the frequent traffic that grind surfaces into new source material.

Field Campaign
In order to quantify the dust emissions, a tower with five Big Spring Number Eight (BSNE) sediment catchers [36] was installed in 2014 for approximately four hours (i.e., 14:00 to 18:00 h) each day over unpaved roads during 12 February in the north at Lajas Maciel and during 21 February in the west at Plan de Ayala areas of TGZ (see Figures 1 and 2, Table 1). Although unpaved roads in the city are common, these sites were selected considering their orientation relative to prevailing wind direction and soil surface conditions of roads that are representative of places located around the city. Road emission measurements were only taken when the observed winds were within approximately ±10 degrees to the road direction. Table 2 lists some soil parameters obtained from soil sample analysis as well as some site characteristics. The roads in this study were considered source unlimited for fine loose particles given the frequent traffic that grind surfaces into new source material. Primary particle size at each site was determined following the specifications established in NOM-021-SEMARNAT-2000 [37] which is based on the Bouyoucos method. Bulk density of loose surface material was determined using a 10 mL test tube [38]. Granulometry analysis of undispersed samples was made using sieves with openings of 9.5, 3.35, 2.0, 1.4, 1.0, 0.71, 0.5, 0.35, 0.25, 0.18, 0.125, 0.075 and 0.037 mm according to ASTM D422/63 [39]. PM10 source available for site emissions was approximated using a polynomial extrapolation from the granulometry analysis according to Wachecka-Kotkowska and Kotkowski [40]. Primary particle size at each site was determined following the specifications established in NOM-021-SEMARNAT-2000 [37] which is based on the Bouyoucos method. Bulk density of loose surface material was determined using a 10 mL test tube [38]. Granulometry analysis of undispersed samples was made using sieves with openings of 9.5, 3.35, 2.0, 1.4, 1.0, 0.71, 0.5, 0.35, 0.25, 0.18, 0.125, 0.075 and 0.037 mm according to ASTM D422/63 [39]. PM 10 source available for site emissions was approximated using a polynomial extrapolation from the granulometry analysis according to Wachecka-Kotkowska and Kotkowski [40]. One Davis meteorological station (Davis Instruments, Hayward, CA, USA), model Vantage Vue (VV) with anemometer at 1 m height was installed together with BSNE catchers. Thus, every time the suspension of material was observed, the instant maximum wind speed was recorded in the station console. The equipment array was able to generate data about wind speed over the soil surface during different events. We defined a wind erosion event similar to Feng and Sharratt [41] where wind velocity exceeding the threshold for >20 consecutive seconds and terminated when the velocity was lower than the threshold for greater than 10 min. Based on this definition, a total of eight events were observed, and their associated sediment fluxes were calculated according to the expression presented in van Donk and Skidmore [42,43] and van Donk et al. [44]:

Opening
where q(z) is sediment flux (kg·m −2 ), z is the height of sampler opening above the soil surface (m), and a and b are fitting parameters considering the height and the mass collected by each BSNE catcher and vary for each storm event.
The discharge of sediments passing the BSNE catchers were then determined integrating Equation (1) from 0 to 2.5 m considering the highest BSNE catcher position (see Table 1): where Q is sediment discharge (kg·m −1 ). It should be noted that we were interested in quantifying emissions generated only by natural wind action so when traffic was present, the catchers were covered to avoid mechanical contribution. Similar to Hagen [45], the sediment discharge is divided by the distance to upwind limit to a non-erodible area to provide an approximation of soil loss per unit area of the unpaved road (see Table 2). Considering that an efficiency between 10 and 25 percent for PM10 dust emissions has been observed in BSNE sediment catchers [46,47], a mean value (17.5%) was used in this study. Based on the sediment discharge, the soil loss and the mean efficiency value, the PM10 emissions were estimated according to Equation (3): where EPM10 is PM10 emissions (g·m −2 ·s −1 ), ε is the catcher efficiency, d is the distance to upwind limit (m) and t (s) is the duration of the event (i.e., mean value of 10 s). One Davis meteorological station (Davis Instruments, Hayward, CA, USA), model Vantage Vue (VV) with anemometer at 1 m height was installed together with BSNE catchers. Thus, every time the suspension of material was observed, the instant maximum wind speed was recorded in the station console. The equipment array was able to generate data about wind speed over the soil surface during different events. We defined a wind erosion event similar to Feng and Sharratt [41] where wind velocity exceeding the threshold for >20 consecutive seconds and terminated when the velocity was lower than the threshold for greater than 10 min. Based on this definition, a total of eight events were observed, and their associated sediment fluxes were calculated according to the expression presented in van Donk and Skidmore [42,43] and van Donk et al. [44]: where q(z) is sediment flux (kg·m −2 ), z is the height of sampler opening above the soil surface (m), and a and b are fitting parameters considering the height and the mass collected by each BSNE catcher and vary for each storm event.
The discharge of sediments passing the BSNE catchers were then determined integrating Equation (1) from 0 to 2.5 m considering the highest BSNE catcher position (see Table 1): where Q is sediment discharge (kg·m −1 ). It should be noted that we were interested in quantifying emissions generated only by natural wind action so when traffic was present, the catchers were covered to avoid mechanical contribution. Similar to Hagen [45], the sediment discharge is divided by the distance to upwind limit to a non-erodible area to provide an approximation of soil loss per unit area of the unpaved road (see Table 2). Considering that an efficiency between 10 and 25 percent for PM 10 dust emissions has been observed in BSNE sediment catchers [46,47], a mean value (17.5%) was used in this study.
Based on the sediment discharge, the soil loss and the mean efficiency value, the PM 10 emissions were estimated according to Equation (3): where E PM10 is PM10 emissions (g·m −2 ·s −1 ), ε is the catcher efficiency, d is the distance to upwind limit (m) and t (s) is the duration of the event (i.e., mean value of 10 s). PM 10 emissions were related to the high wind speed at which suspended material was observed (see Table 3). In this way linear regression equations were obtained correlating PM 10 emissions with maximum wind speed from VV. These equations were used to approximate PM 10  PM10 emissions were related to the high wind speed at which suspended material was observed (see Table 3). In this way linear regression equations were obtained correlating PM10 emissions with maximum wind speed from VV. These equations were used to approximate PM10 emissions under different scenarios (Figure 4).

Model Application for Impacted Areas
The ISC3 model was used to simulate the dust plume dispersion to have an approximation of the impacted areas by PM 10 along TGZ. The model is based on the straight-line, steady-state Gaussian plume equation and includes a wide range of options for modeling air quality impacts of pollution sources. Emission sources are categorized into four basic types of sources, i.e., point sources, volume sources, area sources, and open pit sources [26]. Considering representative roads in TGZ as area sources of PM 10 discharge, in this work the ISC3 was used to calculate the ground-level concentration at a receptor located downwind of all or a portion of the source.
In order to identify the areas impacted by dust particles during days presenting high PM 10 concentrations, meteorological and emission input files were created from the field campaign for the ISC3 model execution. The meteorological input file contained wind speed, wind direction, and temperature mean hourly data which were obtained from VP2. Hourly high wind speeds from VP2 were used to create the dust emissions input file using linear regression equations previously obtained from field work data to approximate PM 10 emissions from unpaved roads. Once the ISC3 model had been executed, the PM 10 concentrations generated by the model were plotted to create the dust plume using the ordinary Kriging interpolation method in ArcGIS 10.0 [48].

Field Measurements
As can be seen in Table 2, surface soil samples of two sites were mainly composed of sand. According to particle size analysis, Plan de Ayala presented higher sand content (70.8%) than Lajas Maciel (64.1%). Conversely, clay and silt content were higher in Lajas Maciel (23.9 and 12.0% respectively) than in Plan de Ayala (19.0 and 10.3%, respectively). Bulk density results show that the value in Plan de Ayala was higher (1.6 g·cm −3 ) than in Lajas Maciel (1.4 g·cm −3 ). Lajas Maciel had roads with slopes slightly more pronounced than Plan de Ayala which were determined using a clinometer. As can be seen in Table 3, these differences may contribute to the discharge of sediments in Plan de Ayala being higher than in Lajas Maciel. Granulometry showed that there was a low natural PM 10 content at the two sites, but saltation and suspension size material content were higher in Plan de Ayala than in Lajas Maciel. However, Mirzamostafa et al. [49] cited that abrasion of soil clods accounted for 14-27% additional suspension size material than available as loose erodible material in the soil. Hagen [50] also found that breakdown of saltating aggregates during erosion contributed on average, 4.9% of PM 10 emissions greater than that originally present in the source material for a range of soils from nine western U.S. states.
Field aeolian emission events from unpaved roads were measured at sites in order to characterize emissions for representative road conditions in Lajas Maciel during 12 February and in Plan de Ayala during 21 February. The discharge of sediments and PM 10 emissions were calculated from sediment collected in BSNE catchers installed over unpaved roads in two different sites in TGZ (see Table 3). Linear fits (LF) relating maximum observed wind speeds from VV and PM 10 emissions are shown in Figure 4. The LF permitted us to know the quantity of PM 10 emitted in different scenarios as a response to different values of wind speed. In this sense, the emission input file to simulate particulate dispersion during days with high PM 10 concentration using the ISC3 model was created using LF and hourly high wind speed from VP2. PM10 emissions were related to the high wind speed at which suspended material was observed (see Table 3). In this way linear regression equations were obtained correlating PM10 emissions with maximum wind speed from VV. These equations were used to approximate PM10 emissions under different scenarios (Figure 4).   In Lajas Maciel (R 2 = 0.808) and Plan de Ayala (R 2 = 0.944) a positive linear trend was observed between PM 10 emissions and high wind speeds. This is as expected due to near unlimited source of fine loose surface material which is characteristic of the many unpaved roads located around TGZ (Figure 4).
Wind erosion emissions were observed when wind speeds were ≥9 m·s −1 in Lajas Maciel. By comparison, in Plan de Ayala, the erosive process was initiated when surface wind speeds of ≥6 m·s −1 were present (Figure 4). These wind threshold differences could possibly be due to the greater content of erodible size source at Plan de Ayala (Table 2). Storm events at this site were also visually more intense. Considering the results obtained from the field work, dust emissions indicate not only the influence of high wind speeds (Figure 4) but also wind turbulence generated by soil surface conditions and street layouts (Figure 2). The Minivol measured the 24-h average PM 10 Figures 1 and 5) allowed the particulate transport to monitoring site. Note that the wind speeds on days associated with the entrances of CF exhibit a wider distribution (i.e., higher winds throughout the day) compared to adjacent days without a CF event. The wider distribution of wind speed during the CF event would result in greater total wind energy during the CF for driving the erosion and dust dispersion processes. Considering that the entrances of the CF during those days resulted in high winds in TGZ, the high measured PM 10 concentrations on these days were associated with wind erosion phenomena on unpaved roads aligning with the wind direction.
Some peaks of PM 10 concentrations were present with low measured wind speeds while low PM 10 concentrations were present with high measured wind speeds ( Figure 5). This may be attributed to minor variations in wind direction relative to the road direction, which would be expected to vary measured wind speeds and concentrations. In addition, discharge generated by vehicle traffic or by burning of garden waste near to the monitoring site during those days would affect concentrations. Some peaks of PM10 concentrations were present with low measured wind speeds while low PM10 concentrations were present with high measured wind speeds ( Figure 5). This may be attributed to minor variations in wind direction relative to the road direction, which would be expected to vary measured wind speeds and concentrations. In addition, discharge generated by vehicle traffic or by burning of garden waste near to the monitoring site during those days would affect concentrations.

Simulations
The ISC3 model was applied to simulate the high PM10 concentration events and to identify areas impacted by dust particles. Mean hourly wind speed, wind direction, and temperature data from VP2 were used to define the meteorological input files of the model. PM10 emission input files were generated using hourly high wind speed from VP2 and applying the LF obtained from results of field work shown in Figure 4. The ISC3 model results were processed using the ordinary Kriging interpolation method in ArcGIS to generate dust plumes for 27 February and 8 and 13 March 2014, considering unpaved roads from Lajas Maciel and Plan de Ayala as area sources. As seen in Figure 5, these dates correspond to the entrance of cold front events, characterized by greater than normal winds associated with unstable conditions. Observed PM10 concentration and those generated by the ISC3 model are shown in Table 4. As can be seen, the magnitude of PM10 concentrations were underestimated by the model for the 27 February 2014 event, however the model could reproduce the event. If other areas with similar characteristics to Lajas Maciel and located to the north of the monitoring site had been considered, the results from the ISC3 would have been similar to those observed during that day due to cumulative concentrations from overlap of adjacent plumes. PM10 concentrations were

Simulations
The ISC3 model was applied to simulate the high PM 10 concentration events and to identify areas impacted by dust particles. Mean hourly wind speed, wind direction, and temperature data from VP2 were used to define the meteorological input files of the model. PM 10 emission input files were generated using hourly high wind speed from VP2 and applying the LF obtained from results of field work shown in Figure 4. The ISC3 model results were processed using the ordinary Kriging interpolation method in ArcGIS to generate dust plumes for 27 February and 8 and 13 March 2014, considering unpaved roads from Lajas Maciel and Plan de Ayala as area sources. As seen in Figure 5, these dates correspond to the entrance of cold front events, characterized by greater than normal winds associated with unstable conditions. Observed PM 10 concentration and those generated by the ISC3 model are shown in Table 4. As can be seen, the magnitude of PM 10 concentrations were underestimated by the model for the 27 February 2014 event, however the model could reproduce the event. If other areas with similar characteristics to Lajas Maciel and located to the north of the monitoring site had been considered, the results from the ISC3 would have been similar to those observed during that day due to cumulative concentrations from overlap of adjacent plumes. PM 10 concentrations were over-estimated by the model for 8 and 13 March, where the presence of tall vegetation (e.g., trees) and buildings around the monitoring site, represent natural barriers that affect the particles distribution over the site which is not seen by the model. This explains why the model generates a higher concentration than observed. It should be noted that the US-EPA set their National Ambient Air Quality Standards for PM 10 at 150 µg·m −3 of 24-h average concentration "to protect and enhance the quality of the Nation's air resources so as to promote the public health and welfare and the productive capacity of its population" [51]. The standard for PM 10 in Mexico however was established by NOM-025-SSA1-2014 at 75 µg·m −3 of 24-h average concentration [52]. Figure 6 shows the simulated distribution of the PM 10 dust plume coming from Lajas Maciel during 27 February 2014. The plume traveled from north-west to south-east of the city starting in Lajas Maciel. It can be seen that the monitoring site was affected by dust emissions. The high PM 10 concentration reported by the Minivol sampler during 27 February 2014 ( Figure 5) can be likely explained by the emissions coming from other bare areas that may be potential dust sources (e.g., other unpaved roads and bare land) located to the north of the monitoring site that were not considered in this study. Vehicle traffic could represent another important contribution of dust emissions affecting air quality according to Etyemezian et al. [12] and Kuhns et al. [13]. During this day, a PM 10 concentration of 1307 µg·m −3 was obtained from the sampler, while a concentration of 37 µg·m −3 is shown by the modeled plume over the monitoring site (see Figure 7).
Surface distribution of the PM 10 dust plume from Plan de Ayala for March 8th is shown in Figure 7. The plume affected mainly the north-western part of TGZ but meteorological conditions enhanced the dust transport to the central part of the city. The ISC3 model showed that PM 10 concentrations at the monitoring site could be influenced by the emissions from unpaved roads located in Plan de Ayala. During this day the model reproduced a concentration of 177 µg·m −3 while the Minivol sampler reported 114 µg·m −3 (see Figures 5 and 7). Figure 8 shows the surface distribution of the PM 10 dust plume from Plan de Ayala for 13 March 2014. Variation in wind direction during the day generated two plumes that affected a large part of TGZ. The ISC3 model showed that PM 10 concentrations at the monitoring site could be influenced by the emissions from unpaved roads located in Plan de Ayala. During this day the model reproduced a concentration of 403 µg·m −3 while the Minivol sampler reported 63 µg·m −3 (see Figures 5 and 8). Considering these results and the presence of tall vegetation and buildings around the monitoring site, it is likely that part of the population was exposed to particulate concentrations exceeding the maximum permissible values set in Mexico environmental standards [52].
The west part of the city is mainly affected by emissions coming from Plan de Ayala generating concentrations higher than the PM 10 standard for Mexico (75 µg·m −3 of 24-h average concentration) [52] while the Central part of the city could be affected considering the prevailing wind direction. Emissions coming from Lajas Maciel affected mainly the north part of TGZ, but the central and east part could be affected as well. Considering the results during simulation days, emissions from Plan de Ayala impacted a larger population area with high PM 10 concentrations than those from Lajas Maciel (Figures 6-8). This is a consequence of the wind speeds presented during those days and to the characteristic of sites affecting dust emissions as previously discussed ( Figure 5 and Table 2). Considering the population density variation of TGZ and the impacted areas, emissions from Plan de Ayala represent a greater environmental risk for population health. explained by the emissions coming from other bare areas that may be potential dust sources (e.g., other unpaved roads and bare land) located to the north of the monitoring site that were not considered in this study. Vehicle traffic could represent another important contribution of dust emissions affecting air quality according to Etyemezian et al. [12] and Kuhns et al. [13]. During this day, a PM10 concentration of 1307 µg·m −3 was obtained from the sampler, while a concentration of 37 µg·m −3 is shown by the modeled plume over the monitoring site (see Figure 7). material increases, causing an increase in soil susceptibility to wind erosion [11,35]. The presence of high wind speeds yields fugitive dust emissions impacting the region.   at an elevation of 522 m above sea level (masl) (see Figure 1). This is a tropical zone with mean annual precipitation oscillating from 900 to 1200 mm·year −1 [34] and surrounding areas consisting of forests. From November to the middle of March, the region is influenced by CF entrances over southern Mexico, which generate high wind intensities (up to 10 m·s −1 ) coming from the west and north-west. The presence of two natural barriers, Mactumatzá Hill at the south and geological formations of the Cañón del Sumidero National Park at the north, both with elevations above 1000 masl, influence the prevailing wind direction and subsequent dispersion of air pollutants affecting the TGZ air quality. Agricultural areas outside of the city have been subjected to strong social pressure due to urban expansion generating unpaved street layouts for new human settlements. When soil aggregates and crusts are crushed by vehicle traffic, the content of loose soil material increases, causing an increase in soil susceptibility to wind erosion [11,35]. The presence of high wind speeds yields fugitive dust emissions impacting the region.  In the current study, Lajas Maciel and Plan de Ayala were only two sites out of many settlements having unpaved roads within TGZ. There are many other sites with the same type of roads which means that it is very probable that PM 10 concentrations and the extent of impacted areas over the city may be significantly greater than those generated by the model. In addition, the presence of burned or other bare areas along the city edge as well as forest fires could represent additional potential sources of PM 10 that affect the air quality of TGZ and it is likely that some of these sources could have contributed to the high PM 10 concentrations reported 27 February.   This is a tropical zone with mean annual precipitation oscillating from 900 to 1200 mm·year −1 [34] and surrounding areas consisting of forests. From November to the middle of March, the region is influenced by CF entrances over southern Mexico, which generate high wind intensities (up to 10 m·s −1 ) coming from the west and north-west. The presence of two natural barriers, Mactumatzá Hill at the south and geological formations of the Cañón del Sumidero National Park at the north, both with elevations above 1000 masl, influence the prevailing wind direction and subsequent dispersion of air pollutants affecting the TGZ air quality. Agricultural areas outside of the city have been subjected to strong social pressure due to urban expansion generating unpaved street layouts for new human settlements. When soil aggregates and crusts are crushed by vehicle traffic, the content of loose soil material increases, causing an increase in soil susceptibility to wind erosion [11,35]. The presence of high wind speeds yields fugitive dust emissions impacting the region.

Conclusions
A preliminary field campaign was carried out to quantify dust emissions and parameterize the ISC3 model for two different sites within TGZ. Aeolian emission events from unpaved roads and simple meteorological inputs were measured at two different sites located within the city to characterize emissions for representative road conditions and to produce ISC3 model inputs. Emissions of PM 10 were measured for eight wind erosion events.
Overall, results show that unpaved roads represent a potential source of dust that affects air quality of urban regions. In this study we found unpaved roads generated emissions ≥1.92 × 10 −3 g.m −2 ·s −1 when winds ≥6 m·s −1 were observed. Air quality of this urban region was also affected as a result of high wind generated under unstable atmospheric conditions of CF, which resulted in at least one event with a measured PM 10 concentration over 1300 µg·m −3 , which exceeds the Mexico national standard. Air pollution events that exceeded the Mexico national standard for PM 10 (≥75 µg·m −3 ) were observed, impacting different areas in the city, representing an environmental risk. Considering soil parameters and street layouts, emissions from Plan de Ayala were higher than those generated in Lajas Maciel. These results indicate that prioritizing dust controls (e.g., paving or application of dust suppressants) on unpaved roads in those north-west parts of TGZ would provide greater reduction in PM 10 exposures to the general population downwind to the south-east than applying controls to the south-east part of TGZ where winds will carry PM 10 to relatively unpopulated areas. However, if paving activities are considered, it is necessary to evaluate materials that could be used to prevent other environmental problems such as urban heat island formation (e.g., low radiation adsorbing pavement) or flooding of low parts of the city during the rainy season due to surface water runoff (e.g., porous pavement).
Dust emissions from unpaved roads by wind erosion are less likely than those generated by mechanical action, but as this study indicates, they are no less important. In this sense, isolated dust emissions events from unpaved roads resulting from wind erosion can affect air quality of urban and even rural areas and can generate air pollution scenarios with PM 10 concentration higher than established standards.

Future Research
The results of this study represent the beginning of air pollution research into the urban air quality of southern Mexico and show the need for air quality monitoring networks in such urban environments. Studies are needed that partition various sources like vehicle traffic, fires, burning of wood and coal in cooking or other activities, dust emissions from agricultural lands, fugitive emissions from construction activities, among others to determine their contribution to air quality in TGZ. During the field measurement portion of this study, there were only a few events where we observed both high winds that were along the direction of the road. Therefore, in order to improve modelling of dust emissions from unpaved roads, future research should include detailed measurements of source sites, measurements of more storm events, the effect of natural barriers on dust dispersion, as well as wind speed profiles at the eroding surface for friction velocity determination. More advanced models, such as WRF-CHEM, should also be used for better simulation of urban settings to understand the problem.
Considering the accelerated growth of TGZ and other areas in the region, the forecast of model emissions and dispersion paths of PM 10 for new unpaved street layouts and other development for urban areas could represent an important tool for land planning. Impacts of PM 10 on human health, dry deposition studies, traffic contributions, effect of turbulence generated by street layouts and roughness on dust emissions, the effect of the orientation of unpaved roads relative to wind direction and particle composition analysis are complementary areas of interest that are needed.