[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Application of Intelligent Dynamic Bayesian Network with Wavelet Analysis for Probabilistic Prediction of Storm Track Intensity Index
Previous Article in Journal
Changes in Cold Surge Occurrence over East Asia in the Future: Role of Thermal Structure
Previous Article in Special Issue
The Variability of Ozone Sensitivity to Anthropogenic Emissions with Biogenic Emissions Modeled by MEGAN and BEIS3
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Atmosphere
Volume 9
Issue 6
10.3390/atmos9060223
atmosphere-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of the Impact of Ozone on Four Economically Important Crops in the City Belt of Central Mexico

by
Luis Gerardo Ruiz-Suárez
1,*,
Bertha Eugenia Mar-Morales
1,
José Agustín García-Reynoso
1,
Gema Luz Andraca-Ayala
1,
Ricardo Torres-Jardón
1,
José Santos García-Yee
2,
Hugo Alberto Barrera-Huertas
2,
Arturo Gavilán-García
3 and
Roberto Basaldud Cruz
3
1
Centro de Ciencias de la Atmósfera, Universidad Nacional Autónoma de México, Circuito de la Investigación Científica s/n, Ciudad Universitaria, Ciudad de México 04510, Mexico
2
Programa de Posgrado en Ciencias de la Tierra. Unidad de Posgrado, edificio D, 1er piso. Circuito de Posgrados, Ciudad Universitaria. Delegación Coyoacán, Ciudad de México C.P. 04510, Mexico
3
Instituto Nacional de Ecología y Cambio Climático (INECC), Blvd. Adolfo Ruíz Cortines 4209, Jardines en la Montaña, Delegación Tlalpan, Ciudad de México C.P. 14210, Mexico
*
Author to whom correspondence should be addressed.
Atmosphere 2018, 9(6), 223; https://doi.org/10.3390/atmos9060223
Submission received: 18 January 2018 / Revised: 15 April 2018 / Accepted: 20 April 2018 / Published: 11 June 2018
(This article belongs to the Special Issue Tropospheric Ozone and Its Precursors)
Browse Figures
Figure 1
<p>Study area with orography and boundaries between federal states and Mexico City. Four data point categories are shown.</p> "> Figure 2
<p>Conceptual model of the process used to obtain the AOT40 exceedance maps.</p> "> Figure 3
<p>The average passive monitoring values at any site (Co) and the Mesoscale Climate and Chemistry Model (MCCM) average (Cm), which is obtained by averaging the profile with the hourly means from MCCM calculated over the two-month passive sampling period.</p> "> Figure 4
<p>Evolution of the hybrid exceedance maps as observational data were gradually added to correct the maps obtained only with the air quality model. (<b>a</b>) Using only modeled and corrected regular grid points. (<b>b</b>) Panel (a) + RAMA + REMA. (<b>c</b>) Panel (b) + García-Yee + Barrera-Huertas. (<b>d</b>) Panel (c) + MILAGRO + CARIEM + Tula.</p> "> Figure 5
<p>Maize: (<b>a</b>) Yield (t/ha). (<b>b</b>) Relative yield loss (%). (<b>c</b>) Production (t). (<b>d</b>) Crop production loss (t). (<b>e</b>) Production value (MXN). (<b>f</b>) Economic loss (MXN), for 2011.</p> ">
Versions Notes

Abstract

:
In this work, we report the economic impact of exposure to high ozone concentrations on four important crops in the area of influence of the Mexico City Megalopolis. Estimated yield losses were as follows: maize: 3%; oats: 26%; beans: 14%; sorghum: 15%. The information needed to estimate the impact of air pollution in Mexico is decidedly deficient. Regarding ozone, the coverage provided by the monitoring networks is strongly focused on urban monitoring and its consistency over time is highly irregular. Apart from the Mexico City Metropolitan Area (MCMA) and less than a handful of other cities, the quality of the data is poor. Ozone in rural areas can be estimated with air quality models. However, these models depend on a high-resolution emissions inventory, which has only been done through validation processes in the MCMA. With these limitations, we set out to estimate the economic impact of exposure to ozone in these crops with a varying degree of sensitivity to ozone in the city belt of Central Mexico. To this end, we developed a procedure that makes optimal use of the sparse information available for construction of AOT40 (accumulated exposure over the threshold of 40 ppb) exceedance maps for the 2011 growing season. We believe that, due to the way in which we dealt with the sparse information and the uncertainty regarding the available data, our findings lie on the safe side of having little knowledge such that they may be useful to decision-makers. We believe that this procedure can be extended to the rest of the country, and that it may be useful to developing countries with similar monitoring and modeling capacities. In addition, these impacts are not evenly distributed in the region and sometimes they were greater in municipalities that have a higher index of poverty. Air pollution arriving from urban areas increases the social inequalities to which these already vulnerable populations are exposed.

1. Introduction

Ozone is a delocalized free radical with high oxidant potential. In the troposphere it is formed by the photolysis of NO2 and by lightning [1]. Subsequent titration with NO yields a steady state driven by UV irradiance. In the presence of VOCs, the peroxy radicals (RO2•), formed during their photo degradation, oxidize NO, opening a reaction path to NO2 without consuming O3, thus allowing its accumulation. In urban areas VOC and NOx emissions enhance these processes, both formation and titration of O3. There, O3 production is generally VOC-limited. In rural and periurban areas with lower NOx levels, O3 production is NOx-limited [2]. Increased ozone production occurs during times of high temperature and solar radiation, such as the stagnant high pressure systems in the spring and summer [3]. The different spatial distributions of NOx and VOC production as well as NO destruction of ozone often result in higher ozone concentrations downwind of urban areas, rather than in the urban areas themselves, where they hinder plant growth [4,5]. Ozone can also be brought into a region by long range transport or by stratosphere-troposphere exchange [6].
Plants exposed to O3 show a reduction in growth [7,8], which in commercial crops means lesser yield and revenue. In the year 2000, economic losses attributed to the worldwide adverse effect of ground-level ozone on maize, soybean, and wheat were estimated at $11–18 billion USD [9].
One of the exposure indexes most closely associated with the injury observed in vegetation is the AOT40 cumulative exposure (accumulated exposure over the threshold of 40 ppb) index. This is calculated as the sum of the differences between the hourly mean ozone concentration and 40 ppb when the concentration exceeds 40 ppb during daylight hours. This value accumulates during periods of time that are defined based on the types of vegetation: usually three months for crops [10]. The exposure–response functions for the AOT40 index were derived from studies conducted under experimental conditions in which environmental conditions favoring the development of the plant and the maximum absorption of the pollutant were ensured [11].
There are also dose–response functions, which indicate the response of the plant to the dose of ozone absorbed using the calculation of ozone flow rates from the air to the plant [12]. By analogy with the exposure–response functions, a flow rate index for ozone accumulated during the plant’s growth period is used. This flow rate is obtained by adding the hourly ozone flow rates over a given growth period that is defined depending on the type of vegetation. The calculation of the ozone flow rate takes into account the ambient ozone concentration and stomatal conductance [13]. Dose–exposure functions are easier to use and are more widely used [2,14].
Geographic information systems (GISs) can be used in order to estimate the economic impact of ozone on crops [15,16]. O3 maps consistent with the chosen function are made and overlaid on crop maps in order to apply those functions [15,16,17]. However, worldwide, there are much fewer rural ozone monitoring networks than urban ozone monitoring networks. Ozone monitoring networks are usually established in cities in order to protect the health of the population. This deficiency is even more visible in developing countries. As of 2011, in Mexico, there was not a single continuous air quality monitoring station specifically devoted to rural monitoring, and in the city belt of Central Mexico, there were three air quality monitoring systems (AQMSs) reporting ozone records from 36 monitoring stations to the “Sistema Nacional de Indicadores de Calidad del Aire” (SINAICA) (National Air Quality Information System) [18], not all of which have the same quality standards as the Mexico City Metropolitan Area (MCMA) network [19]. The data thus obtained are not sufficient to create exceedance maps for the critical value AOT40 in the region. It is well established that ozone concentration increases downwind of urban areas due to the absence of fresh of NOx emissions, from which NO can titrate O3 [20]. Creating maps by interpolating only data from urban stations implies underestimation of the concentrations over all the rural areas in the city belt. Correlations between urban and rural sites [21] cannot be obtained in the complete absence of rural sites. The orographic barriers that separate the air basins in the center of the country impose further restrictions on the interpolation methods, when only using data from the available monitoring stations. Figure 1 shows how the monitoring stations (empty white circles) cluster in three groups. The largest group at the center corresponds to the MCMA, the cluster to the west corresponds to the Toluca Metropolitan Area, and the group located towards the east corresponds to the Puebla Metropolitan Area.
Passive samplers are regarded as a practical choice for measuring ambient ozone concentrations for the purposes of atmospheric chemistry and ecological assessment [22,23]. They are inexpensive and easy to deploy in the field, requiring no power supply or supporting equipment, and they also offer flexibility of placement. These features make them attractive for monitoring in rural areas and in difficult environments, such as within the forest canopy [24]. However, they only provide continuous cumulative values without differentiating day from night. To obtain AOT40 estimates based on measurements with passive samplers, it is necessary to correct the total quantity of O3 provided by the passive sampler to the equivalent ozone levels for daylight hours above 40 ppb. The use of mathematical functions from the Weibull probability family applied to nearby continuous ozone measurements has been the tool used in order to approximate the distribution of hourly O3 concentrations [25]. However, the resulting model must be verified against monitoring stations in the study area, which in the case of Central Mexico are available only in urban areas. Nevertheless, some authors have proposed the use of the concept of relative altitude. This is an empirical statistical approach based on using the difference between the altitude at each point at which ozone levels have been measured and the lowest altitude in each circular area around each measured point [26,27], but this approach may not apply in areas that are under the influence of nearby sources of ozone precursors; it is subject to considerable differences between each particular site in the region [26]. Therefore, exceedance maps for the city belt of Central Mexico cannot be obtained using the statistical modeling nor the relative altitude approaches.
Air quality models have also been used to estimate the formation and transport of ozone over the areas covered by the vegetation of interest [15,28], but the results of the model need to be validated against observations. Additionally, if available, data from monitoring networks can be used to adjust the results of air quality models [29,30,31]. Tarrasón et al. [29] created difference fields in the simulation domain by interpolating the differences observed between the results of a regional air quality model and the measurements taken. This technique allows for correction of the results of the model at any point in the simulation domain. In our case, the absence of rural monitoring stations precluded us from creating such difference fields outside of urban areas.
Damage to vegetation by ozone has been reported for Mexico. These reports have been of a qualitative nature [32,33], or with too coarse resolution in global modeled studies [17]. The same lack of information about ozone in rural areas that hindered quantitative estimates in those early studies still prevails. In this work, we report for the first time in Mexico, the estimated cost of agricultural production loss due to the exposure of several crops of economic importance to high concentrations of ozone present in the city belt of Central Mexico, which includes the Mexico City megalopolis. We address the methodological challenges facing an emerging economy such as that of Mexico in the estimation of the economic impacts of crop ozone exposure and we describe the equity issues that such exposure exacerbates.

2. Methods

2.1. Study Area

The study region, located in the central region of Mexico, is demarcated by the coordinates 20.6200–18.2276° N and 100.2030–97.5693° W and has an area of 72,900 km2. It encompasses the states of Morelos, Mexico State, Tlaxcala, and Mexico City and includes parts of the states of Guerrero, Hidalgo, Queretaro, and Puebla. The total population of this region is 33.77 million inhabitants, which represents 31.2% of the total national population [34]. It has an estimated surface area of rain-fed crops of 28,154 km2 [35].
The following percentages show these states’ participation in the current total national GDP values (2011) by primary activity (agriculture, animal breeding and farming, forestry, fishing and hunting): Morelos: 1.2%; Mexico State: 3.5%; Tlaxcala: 0.5%; Mexico City: 0.2%; Guerrero: 2.4%; Hidalgo: 1.9%; Queretaro: 1.3%; Puebla: 4.1% [36].

2.2. Ozone Data

2.2.1. Passive Monitoring

The passive sampler (Ogawa & Co., USA, Inc., Pompano Beach, FL, USA) consists of a cylindrical polymer body (2 × 3 cm) and a plastic pin clip holder. There are two cavities on the ends of the cylinder, each of which holds one coated filter between two stainless steel screens. The passive sampler uses two nitrite-coated filters in the interior part. When exposed to ambient O3, nitrite ion is oxidized to nitrate ion. Both exposed filters are extracted and analyzed by ion chromatography to determine cumulative nitrate concentration for the two-week exposure period. This allows for a determination of average cumulative O3 exposure in ppb [37].
For analytical determination, a Shimadzu HPLC model CDD-10AVP was used. It is equipped with an isocratic pump LC-20AD and an ionic conductivity detector CD-10A-VP. The column was Shodex SI-90 4E (4.0 × 250 mm, particle 9 μm). The mobile phase was 1.8 mM Na2CO3 + a 1.7 mM NaHCO3 solution, with a flow of 2.0 mL/min and a 200 mL injection volume. Analysis time was 15 min.
Our implementation of this method was tested by co-colocation of passive samplers with a standard continuous O3 monitor. In two independent experiments, the Ogawa passive samplers always overestimated accumulated O3 concentrations by no more than 10% in samples taken between 1 day and 3 weeks. The remaining nitrite in the passive sampler followed an exponential decay curve as the sampling period increased, staying in its linear-like part for up to two weeks.
The passive sampling was conducted in a smaller area, the mountain zone south of the MCMA (Figure 1, triangles with dots).

2.2.2. Continuous Monitoring

The MCMA continuous ozone monitoring data were obtained from the databases of the environment authority of Mexico City [38]. Data from the Metropolitan Area of Puebla and from the Metropolitan Area of Toluca were obtained from SINAICA [18]. All of them are represented by empty white circles in Figure 1. Short-term (up to five weeks) continuous data were obtained in several rural and periurban sites within the region using mobile units as part of studies funded for other purposes [39,40,41] (Figure 1, stars).

2.3. Maps

Maps were created using ArcGis software v.10.2 by ESRI, Redlands, CA, USA and consisted of five layers in shapefile format, with Lambert conformal conic projection (LCC), Datum International Terrestrial Reference Frame 92 (ITRF92) epoch 1988.0 (DATUM ITRF92). These layers are as follows: state, municipal, land use administrative boundaries [42], the study area limit (own elaboration), and the final layer, in which the spatial distribution of the ozone exposure index studied in this work, AOT40 for the year 2011, is calculated and represented, with our data.
The maps were created on a 1 ha scale. This exceeds the spatial resolution of the information available to us on agricultural production. It also exceeds the spatial resolution in our Mesoscale Climate and Chemistry Model (MCCM) runs (3 km × 3 km). However, we kept this resolution because we expect that the government agencies that have information at the parcel level will eventually make it available. We are developing the capability of running WRF-Chem, instead of MCCM at 1 km2 for Central Mexico. For AOT40, the selected interpolation method was kriging, which is widely applied for interpolating air quality data [31,43,44].

2.4. Exposure–Response Functions

The exposure–response functions (ERFs) from the literature were derived from open-top chamber studies in which the effect on the growth of plants exposed to different concentrations of ozone was measured. The functions compiled by Mills et al. [12] were used; these functions give a linear exposure–response relationship as a function of the AOT40 metric for the crops studied (Table 1). For sorghum, oats and beans, in the absence of specific functions we choose to use the generic sensitive function suggested by Mills. The functions were modified to force RY = 1.0 at x = 0 and RY = 0.95 at the critical level [16,17].
Production loss per unit area was calculated using the following equation:
C P L i = R Y L i 1 − R Y L i ∗ C P i
where CPLi is the crop production loss [t] per unit area i [ha], RYLi is the relative yield loss, and CPi is the actual crop yield [t/ha] [16].
Municipal crop production loss used in this study is the mean of the crop production loss in all the crop cells within each municipality.
Finally, the economic damage was calculated using the following equation:
E L i = C P L i ∗   A R P i
where ELi is the value of the economic loss [$/ha], and ARPi was the average rural price per municipality [$/t] in 2011. In Mexico, this is the highest resolution crop price available. Monetary units are current Mexican pesos in 2011. Final aggregate values are given in USD using the average 2011 exchange rate.

2.5. Statistical Information

The first tests were carried out using information from the database of “Oficina Estatal de Información para el Desarrollo Rural Sustentable” (OEIDRUS) [45], which provided the spatial distribution of the plots of the crops of interest, with a resolution of 1 ha. Subsequently, we realized that the information available to us regarding other states was very irregular. In order to standardize procedures, we decided to use statistical information available from “Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación” (SAGARPA) [46], which provides the following municipal scale information by crop: hectares planted, hectares harvested, average yield per hectare, production value, and average production value.
In order to include a measure of environmental equity associated with externalities arising from air pollution we used the marginalization index, which is defined as a measurement summary that allows for differentiation between locations in the country according to the overall impact of the deficiencies suffered by the population as a result of their lack of access to education and assets, and the fact that many reside in inadequate housing, all of which is associated with a low income [47].

2.6. Air Quality Model

The air quality model (AQM) used was the MCCM [48] in its non-hydrostatic mode. The chemical model is consistent with the meteorological model, using both the same transport schemes and the same horizontal and vertical grid. The chemical mechanism was RADM2 [49]. The meteorological information was obtained from the global forecast (GFS). The modeling domain has 73 × 73 cells of 3 × 3 km and 22 vertical levels. The parametrizations for the cloud radiation scheme were as follows: non-hydrostatic option, Grell’s cumulus parameterization, the Burk–Thompson PBL scheme, warm-rain scheme, Reisner scheme for explicit moisture, and a five-layer soil model.
A comprehensive assessment of model performance may be done using a set of metrics. Among them, the index of agreement Ic [50] allows for a quick and intuitive verification of model performance.
I c = 1 − ∑ i = 1 N ( p i − o i ) 2 ∑ i = 1 N ( | p i − o ¯ | + | o i − o ¯ | ) 2
where p i is the model-predicted concentration at time interval i. o i is the hourly average of the observed concentration at time interval i. I c varies from 0 to 1. When I c = 1 , there is a perfect match between modeled and observed concentration values.

3. Results and Discussion

3.1. Creation of the Exceedance Maps

The limiting factor for estimating crop loss due to ozone exposure is the availability of cumulative AT40 exceedance maps. The method developed for arriving at a reliable estimate, on the safe side of having sparse information, was the first result of this work. Figure 2 shows a conceptual model of the procedure followed in order to obtain the exceedance maps with the information available in Mexico. The following is a description of the contents of each of the model’s nine boxes.
● Box 1
The first layer of information necessary to obtain the AOT40 exceedance maps was the results from MCCM for January 2009 to July 2011 for the regular grid shown in Figure 1 (white circles with dots). The high spatial and temporal resolution emissions model used in the simulations (INEM-AR 2008, Spanish initials), was obtained from the National Emissions Inventory disaggregated to a high spatial resolution (3 km × 3 km) using by default temporal emissions profiles or local profiles if available. As is commonly necessary in the MCMA, NOx and VOC emissions were adjusted in order to optimize the adjustment of the O3 profiles modeled using the observations of the local monitoring network [51,52]. In any case, the first task when using the model in each specific application is to evaluate its performance by comparing the results of the model against observations of the MCMA monitoring network during the episode or scenario under study [53] (Figure 1, white empty circles). The O3 comparison between the model and data from several MCMA monitoring stations has an index of agreement greater than 0.81 [54]. The index value decreased to 0.5 < I c < 0.66 when compared with the rural and periurban sites monitored during two intensive campaigns [40,41]. The aim of some tasks in the following boxes is to improve on this metric.
● Box 2
With a network of seven passive monitors, O3 samples were collected in the Sierra de Chichinautzin (Figure 1, triangles with dots) between 2011 and 2012. In the central region of Mexico, the seasons are divided into cold dry (November–February), warm dry (March–May), and wet (June–October) [55]. Passive samplers were placed in four periods of two weeks each, for two months in each season. Thus, eight samples were collected for each season during the 24 months of the study.
● Box 3
To prepare the results from the model for comparison with the measurements taken by the passive samplers, which are of lower temporal resolution, the hourly ozone concentrations reported by MCCM were averaged over the two months of sampling for each season, resulting in the MCCM profile shown in Figure 3. The averages over each point of the grid, 12 sites, in the part of the simulation domain that contains the passive monitor sampling region were interpolated in the GIS in order to create a field of modeled averages and thus extract the averages over the passive monitoring sites for comparison.
● Box 4
The average value of the modeled field extracted at the passive sampling sites (Cmi) was compared with the average of the results observed with the passive samplers during the two months at the site ( C o i ) as follows:
C m = ∑ i = 1 n C p i n
where Cm is the mean of the modeled concentrations from MCCM for the period of time similar to the passive monitoring time. Cp is the concentration predicted by the simulation model throughout the simulation period. n is the number of time intervals for the entire simulation period. Simultaneously,
C o = ∑ j = 1 l C j l
where Co is the mean of the observed concentrations obtained from the passive monitoring operations for the sampling period. Cj is the concentration obtained by a passive monitoring operation for a period of 15 days. l is the number of samples collected during a campaign, be it during the cold dry, warm dry, or wet seasons.
For each passive monitoring site,
x = C o C m .
At this point, a field of adjustment quotients such as that used by Tarrasón et al. [29] can be created, but a review of the exceedance map shows that the region where an interpolation can be applied is very small. Therefore, at this stage, we opted for an average adjustment factor for the MCCM results. An average adjustment factor for the results of the model in the entire region was obtained as follows:
X ¯ =   1 n ∑ i = 1 n x i .
Table 2 shows the correction factors obtained for the passive sampling sites. The wet season data could not be used because, in that season, at all those sites, saturated relative humidity conditions are persistent; this causes condensation in the passive sampler. There is no specific pattern. Tres Marías is a small community in the State of Morelos that is under the influence of the Mexico City-Cuernavaca toll and non-toll highways. However, Parres, in Mexico City, is also under the influence of the same highways. Only Desierto de los Leones site showed the opposite ratio.
● Boxes 5 and 6
In 2011 and 2012, our group carried out three intensive monitoring campaigns at several sites in two regions of the city belt of Central Mexico. The sites selected were rural or periurban. From the 2011 campaign [40], the Amecameca site was chosen. This site, located in the State of Mexico at the Tenango del Aire mountain pass, connects the Valley of Mexico with the valleys of Cuautla and Cuernavaca. From the 2012 campaign [41], the Chipilo site was chosen, and they were used as described below.
With the adjustment factor obtained in Box 4, the results of the model were corrected at all grid points and at the coordinates of the mobile units in Amecameca and Chipilo for the periods of time during which the campaigns were carried out; these corrected results were compared with the hourly profiles observed at those sites. At the urban sites, the comparison showed that the model results corrected with the average adjustment value obtained from passive monitoring considerably overestimated the reported values; however, for Amecameca and Chipilo, the comparison showed that the daytime profiles of the corrected model are more similar to the observations, with an overestimate between 10 and 30%. For completeness, a small additional correction was made. The nighttime profiles are also overestimated, a possible sign of the absence of enough nocturnal NO emissions in the emissions inventory. During nighttime, excess modeled boundary layer height and/or wind intensity could also reduce the NO concentrations available to titrate O3. These results are also valuable in top-down validations of emissions inventories, but that is the object of another work to be reported elsewhere. With two sites, a nighttime adjustment factor of 0.47 was obtained, from 19:00 to 07:00 h. Its impact on the AOT40 indicator is small, producing an effect only at dawn and at dusk. Nevertheless, it was applied so as not to overestimate the calculation of the AOT40 indicator for those two hours. At this stage, comparison of modeled and observed concentrations in the rural and periurban available pairs of data yielded much better values of the index of agreement 0.62 < I c < 0.90 .
● Box 7
With the results from MCCM optimized for the MCMA, the hourly ozone profiles were obtained for the 2011 crop season. These profiles, for all sites in the regular grid, were corrected with the daytime and nighttime adjustment factors, and the AOT40 indicator was obtained for all points of the regular grid, using the kriging interpolating method (Map a, Figure 4).
● Box 8
The values of the AOT40 indicator were obtained for all monitoring sites with data from the local MCMA network (RAMA) and that of the Metropolitan Area of Puebla (REMA) and added to the interpolating data set yielding Map b in Figure 4. These allowed for the correction, in urban areas, of the overestimation in the interpolation introduced by the indiscriminate use of the adjustment factor applied to the results of the model. In the absence of sufficient data from a station, the exceedances for the growing period were estimated using a ratio proportional to the number of observation days [26].
● Box 9
The AOT40 indicator calculated for the rural and periurban sites in the 2011 [40] and 2012 [56] campaigns were then appended to the interpolating data set, yielding Figure 4c.
At this point, we still have some data from a few rural sites not yet used and must decide if it is safe to use them. Some were from the campaign carried out in 2009 in the State of Morelos (CARIEM) [39]. We also had data from some sites, obtained by the mobile unit of the air quality monitoring system of the State of Mexico between 2006 and 2009, two sites, Tula [57] and ININ [58] plus data for periurban sites during MILAGRO [59]. The AOT40 indicator was estimated for all of these sites. We decided it was safe to use them after the following empirical analysis.
In total, we had 16 data sets from short time campaigns. Two (AME and CHI) were used to confirm the use of the single average adjustment factor. Two from the 2011 campaign, OZU and TEN, and three from the 2012 campaign, AMO, CA, and HUA, had already been used to obtain the map in Figure 4c. To produce the final map (Figure 4d), we added the estimated AOT40 exceedance of the remaining nine rural and periurban sites to the interpolating data and analyzed the resulting map. Altogether, when making Figure 4c,f, of the 16 sites, 11 of them warmed the colors to the already adjusted map. Five of them modulated an excess correction by the average adjustment factor
In the absence of a region wide of several years of AOT40 exceedance data, Figure 4d offers the best representation possible of the AOT40 exceedance using the available information for July–September months. We call it a hybrid map, as it uses results of a mesocale air quality model fed with the 2008 INEM-AR emissions inventory and driven by 2011 meteorology. These estimates were corrected with passive sampling data collected during 2011 and 2012 and with short-term continuous monitoring data in rural and periurban sites, 16 sites, in four field campaigns in 2006, 2009, 2011, and 2012.
Figure 4d is consistent with some known facts that help to explain the shapes and shades on the maps. The parts of the map with the greatest accumulation of exceedances, the hottest parts, are the high-altitude O3 in the higher parts of the Sierra Nevada to the east of the MCMA, the Sierra Chichinautzin to the south, and the Sierra de las Cruces to the west of the MCMA; on the other hand, the red patch to the north of the Toluca Valley Metropolitan Area (TVMA) represents a net receptor area of air masses, rich in ozone and still photochemically active, arriving, depending on the dominant direction of surface winds, from the TVMA [60] to the south-southwest, from the MCMA [61,62] to the east, or from the Tula-Tepeji industrial zone [63] to the north. The impact of the Tula-Tepeji industrial zone on the MCMA is well documented [54,63]; however, due to the presence of the Sierra de las Cruces, a difference of only a few degrees can divert emissions from Tula-Tepeji towards the TVMA instead of MCMA.
The maps in Figure 4 show that the parts with the lower AOT40 exceedances systematically correspond to areas where there is no measurement of ozone, despite the upward correction with a constant factor of the MCCM results. Additionally, in Figure 4a–d, in most cases, as the scarce observational data from rural and periurban sites not previously used were added to the interpolation database, the map showed warmer colors meaning higher AOT40 exceedances. This supports the following hypothesis: It is highly likely that the cumulative estimate of injury to vegetation and economic costs due to ozone exposure in the city belt of Central Mexico will be an underestimate of the actual costs and damages.

3.2. Estimation of Economic Losses

Figure 5a–f shows the results of the application of Equations (1)–(3) to the corresponding layers of information in the GIS. Table 3 shows input values and results from the application of Equations (1)–(3) in the 20 municipalities most affected by the exposure of maize to ozone throughout the simulation domain.
The impacts are not homogeneous in the region; they depend on the production of the crop, on the planted area, the cumulative AOT40, and the price paid to the farmer, which is not the same in all municipalities. For example, in the case of maize, Table 3 shows that the farmers in the municipality of Huamantla in the State of Tlaxcala lost 3% of their production due to ozone exposure. The estimated value of their production is 156 million pesos MXN and the economic loss is 5 million pesos MXN. In the municipality of Villa de Allende in the State of Mexico farmers lost 7% of their production due to ozone exposure. The estimated value of their production is 152 million pesos MXN, and the economic loss is 11 million pesos MXN. In terms of equity, it can be observed that the degree of urban marginalization, used here as a euphemism for poverty, in Huamantla is low, and in Villa de Allende it is high; surely such a potential loss of income represents a high cost to already impoverished communities.
The total production value for maize in these 20 municipalities was 1.590 billion pesos MXN, and, due to ozone exposure, 85 million pesos MXN (5%) were lost. In the entire simulation domain, the value of maize production in 2011 was 6.598 billion pesos MXN; 46,641 tons were lost to yield reduction, and the value of that unfulfilled production was 191 million pesos MXN. This represents 3% of the potential production.
The same procedure was performed to estimate the effect of ozone exposure in oats, beans, and sorghum, which are more sensitive to ozone than maize. Table 4 shows the economic loss for each of the four crops throughout the simulation domain. The total economic loss was 671 million pesos MXN, which at the 2011 exchange rate equals 50 million USD. For the other crops, the relative impact of the degree of marginalization was not uniform in the municipalities that make up the simulation domain. It is also clear that these costs may be underestimated due to the likely underestimation of the AOT40 exceedances in the areas of the map with colder colors.

4. Conclusions

This work is the first attempt to estimate the cost of air pollution in Mexico, in terms of its impact on agriculture. Maize, the most important cereal in the Mexican diet, showed low sensitivity to ozone and economic losses of up to 3%. However, if the other crops are as sensitive as their proxies, then farmers who grow beans may lose up to 14% in potential revenue. Those farming sorghum may lose 15%, and for oats the loss may be up to 26% of the yield under optimal growth conditions. The impact on other species such as vegetables, which are highly sensitive to ozone and are cultivated near the cities, was not estimated, nor was the impact on species that are native to Mexico and economically important in the region, such as nopal and agave, due to the lack of exposure–response functions, even by proxy.
This attempt comes more than 30 years after continuous ozone monitoring was first performed in Mexico. It has been done with almost non-existent information for the affected rural areas. Although it is incomplete, the evidence accumulated by the incorporation of data from rural and periurban sites seems to show that, in terms of AOT40 exceedance maps, the costs reported here may be an underestimation.
The other source of uncertainty or bias is the use of proxy species which points to the need for national or local determination of exposure–response functions for crops of economic importance not reported previously in the literature.
It is evident that the agricultural areas affected by the metropolitan areas of the city belt of Central Mexico are strongly influenced by the cumulative AOT40 exceedances. Some of the affected municipalities are ranked high on the marginality index, tightening the screw on already-impoverished farming communities.
It falls upon the decision-makers and the affected groups to decide if the information provided here is sufficient to initiate actions for the mitigation of the impacts that negatively affect the availability of food and reduce the sustainability of the agricultural activities in the region. The first of such actions should be the long-term rural monitoring of air pollutants nationwide.

Author Contributions

Luis Gerardo Ruiz-Suárez wrote the paper and conceived and designed the project and the field campaigns; Bertha Eugenia Mar-Morales was in charge of the GIS calculations and maps; José Agustin García-Reynoso was responsible for the AQM work; Gema Luz Andraca-Ayala was in charge of the passive sampling, lab and data analysis; Ricardo Torres-Jardón was in charge of the mobile unit observations in all the periurban and rural sites; José Santos García-Yee carried out the day-to-day work and data analysis of the 2011 field campaign; Hugo Alberto Barrera-Huertas carried out the day-to-day work and data analysis of the 2012 field campaign; Arturo Gavilán-García and Roberto Basaldud-Cruz were liaisons with INECC and assisted with information collection.

Acknowledgments

This study was conducted with support from the following projects: Conacyt/Semarnat 2006 23081, Conacyt/Fomix PUE-2007-C01-77293, and UNAM/INECC INECC/A1-004/2014 (Ecaim). The air quality model was run using high performance computing provided be DGTIC-UNAM. We also wish to recognize Patrick Weill for his translation of this work into English, with support from Conacyt 280403 (RedCAM).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thompson, A.M.; Yorks, J.E.; Miller, S.K.; Witte, J.C.; Dougherty, K.M.; Morris, G.A.; Baumgardner, D.; Ladino, L.; Rappenglück, B. Tropospheric ozone sources and wave activity over Mexico City and Houston during MILAGRO/Intercontinental Transport Experiment (INTEX-B) Ozonesonde Network Study, 2006 (IONS-06). Atmos. Chem. Phys. 2008, 8, 5113–5125. [Google Scholar] [CrossRef] [Green Version]
  2. Felzer, B.S.; Cronin, T.; Reilly, J.M.; Melillo, J.M.; Wang, X. Impacts of ozone on trees and crops. C. R. Geosci. 2007, 339, 784–798. [Google Scholar] [CrossRef] [Green Version]
  3. Logan, J.A. Tropospheric ozone: Seasonal behavior, trends, and anthropogenic influence. J. Geophys. Res. Atmos. 1985, 90, 10463–10482. [Google Scholar] [CrossRef]
  4. Gregg, J.W.; Jones, C.G.; Dawson, T.E. Urbanization effects on tree growth in the vicinity of New York City. Nature 2003, 424, 183–187. [Google Scholar] [CrossRef] [PubMed]
  5. Mauzerall, D.L.; Wang, X. Protecting agricultural crops from the effects of tropospheric ozone exposure: Reconciling science and standard setting in the United States, Europe, and Asia. Annu. Rev. Energy Environ. 2001, 26, 237–268. [Google Scholar] [CrossRef]
  6. Ganguly, N.D.; Tzanis, C. Study of stratosphere-troposphere exchange events of ozone in India and Greece using ozonesonde ascents. Meteorol. Appl. 2011, 18, 467–474. [Google Scholar] [CrossRef]
  7. De Bauer, M.D.; Hernandez-Tejeda, T. A review of ozone-induced effects on the forests of central Mexico. Environ. Pollut. 2007, 147, 446–453. [Google Scholar] [CrossRef] [PubMed]
  8. Hayes, F.; Mills, G.; Harmens, H.; Norris, D. Evidence of Widespread Ozone Damage to Vegetation in Europe (1990–2006); ICP Vegetation Programme Coordination Centre, CEH: Bangor, UK, 2007. [Google Scholar]
  9. Avnery, S.; Mauzerall, D.L.; Liu, J.; Horowitz, L.W. Global crop yield reductions due to surface ozone exposure: 1. Year 2000 crop production losses and economic damage. Atmos. Environ. 2011, 45, 2284–2296. [Google Scholar] [CrossRef]
  10. Fuhrer, J.; Skärby, L.; Ashmore, M.R. Critical levels for ozone effects on vegetation in Europe. Environ. Pollut. 1997, 97, 91–106. [Google Scholar] [CrossRef]
  11. Fuhrer, J. Critical level for ozone to protect agricultural crops: Interaction with water availability. Water Air Soil Pollut. 1995, 85, 1355–1360. [Google Scholar] [CrossRef]
  12. Mills, G.; Buse, A.; Gimeno, B.; Bermejo, V.; Holland, M.; Emberson, L.; Pleijel, H. A synthesis of AOT40-based response functions and critical levels of ozone for agricultural and horticultural crops. Atmos. Environ. 2007, 41, 2630–2643. [Google Scholar] [CrossRef]
  13. Musselman, R.C.; Lefohn, A.S.; Massman, W.J.; Heath, R.L. A critical review and analysis of the use of exposure- and flux-based ozone indices for predicting vegetation effects. Atmos. Environ. 2006, 40, 1869–1888. [Google Scholar] [CrossRef]
  14. Krupa, S.V.; Manning, W.J. Atmospheric ozone: Formation and effects on vegetation. Environ. Pollut. 1988, 50, 101–137. [Google Scholar] [CrossRef]
  15. Vlachokostas, C.; Nastis, S.A.; Achillas, C.; Kalogeropoulos, K.; Karmiris, I.; Moussiopoulos, N.; Chourdakis, E.; Banias, G.; Limperi, N. Economic damages of ozone air pollution to crops using combined air quality and GIS modelling. Atmos. Environ. 2010, 44, 3352–3361. [Google Scholar] [CrossRef]
  16. Van Dingenen, R.; Dentener, F.J.; Raes, F.; Krol, M.C.; Emberson, L.; Cofala, J. The global impact of ozone on agricultural crop yields under current and future air quality legislation. Atmos. Environ. 2009, 43, 604–618. [Google Scholar] [CrossRef]
  17. Chuwah, C.; van Noije, T.; van Vuuren, D.P.; Stehfest, E.; Hazeleger, W. Global impacts of surface ozone changes on crop yields and land use. Atmos. Environ. 2015, 106, 11–23. [Google Scholar] [CrossRef]
  18. Instituto Nacional de Ecología y Cambio Climático (INECC). Sistema Nacional de Indicadores de Calidad del Aire. Available online: http://sinaica.inecc.gob.mx/ (accessed on 9 March 2018).
  19. Semarnat. Programa de Gestión Federal para Mejorar la Calidad del Aire de la Megalópolis 2017–2030; Semarnat: Ciudad de México, México, 2017; p. 330. [Google Scholar]
  20. Goldan, P.D.; Trainer, M.; Kuster, W.C.; Parrish, D.D.; Carpenter, J.; Roberts, J.M.; Yee, J.E.; Fehsenfeld, F.C. Measurements of hydrocarbons, oxygenated hydrocarbons, carbon monoxide, and nitrogen oxides in an urban basin in Colorado: Implications for emission inventories. J. Geophys. Res.-Atmos. 1995, 100, 22771–22783. [Google Scholar] [CrossRef]
  21. Coyle, M.; Smith, R.I.; Stedman, J.R.; Weston, K.J.; Fowler, D. Quantifying the spatial distribution of surface ozone concentration in the UK. Atmos. Environ. 2002, 36, 1013–1024. [Google Scholar] [CrossRef]
  22. Krupa, S.V.; Legge, A.H. Passive sampling of ambient, gaseous air pollutants: An assessment from an ecological perspective. Environ. Pollut. 2000, 107, 31–45. [Google Scholar] [CrossRef]
  23. Carmichael, G.R.; Ferm, M.; Thongboonchoo, N.; Woo, J.-H.; Chan, L.Y.; Murano, K.; Viet, P.H.; Mossberg, C.; Bala, R.; Boonjawat, J. Measurements of sulfur dioxide, ozone and ammonia concentrations in Asia, Africa, and south America using passive samplers. Atmos. Environ. 2003, 37, 1293–1308. [Google Scholar] [CrossRef]
  24. Cox, R.M.; Malcolm, J.W. Passive ozone monitoring for forest health assessment. Water Air Soil Pollut. 1999, 116, 339–344. [Google Scholar] [CrossRef]
  25. Krupa, S.; Nosal, M.; Peterson, D.L. Use of passive ambient ozone (O3) samplers in vegetation effects assessment. Environ. Pollut. 2001, 112, 303–309. [Google Scholar] [CrossRef]
  26. Gerosa, G.; Ferretti, M.; Bussotti, F.; Rocchini, D. Estimates of ozone aot40 from passive sampling in forest sites in south-western Europe. Environ. Pollut. 2007, 145, 629–635. [Google Scholar] [CrossRef] [PubMed]
  27. Loibi, W.; Winiwarter, W.; Kopsca, A.; Zufger, J.; Baumann, R. Estimating the spatial distribution of ozone concentrations in complex terrain. Atmos. Environ. 1994, 28, 2557–2566. [Google Scholar] [CrossRef]
  28. Simpson, D. Photochemical model calculations over Europe for two extended summer periods: 1985 and 1989. Model results and comparison with observations. Atmos. Environ. Part A Gen. Top. 1993, 27, 921–943. [Google Scholar] [CrossRef]
  29. Tarrasón, L.; Semb, A.; Hjellbrekke, A.G.; Tsyro, S.; Schaug, J.; Bartnicki, J.; Solberg, S. Geographical Distribution of Sulphur and Nitrogen Compounds in Europe Derived both from Modelled and Observed Concentrations; 71; Norwegian Meteorological Institute: Kjeller, Norway, 1998. [Google Scholar]
  30. Bush, T.; Targa, J.; Stedman, J. UK Air Quality Modelling for Annual Reporting 2004 on Ambient Air Quality Assessment under Council Directives 96/62/EC and 2002/3/EC Relating to Ozone in Ambient Air; Department for Environment, Food and Rural Affairs, Welsh Assembly Government, the Scottish Executive and the Department of the Environment for Northern Ireland: Didcot, UK, 2007.
  31. Denby, B.; Horálek, O.J.; Walker, S.E.; Eben, K.; Fiala, J. Interpolation and Assimilation Methods for European Scale Air Quality Assessment and Mapping Part I: Review and Recommendations; European Topic Centre on Air and Climate Change: Bilthoven, The Netherlands, 2005; p. 51. [Google Scholar]
  32. De Bauer, L.I.; Krupa, S.V. The valley of Mexico: Summary of observational studies on its air quality and effects on vegetation. Environ. Pollut. 1990, 65, 109–118. [Google Scholar] [CrossRef]
  33. Ortíz-García, C.F.; Laguette-Rey, H.D.; de Bauer, L.I. Effects of oxidants in ambient air on annual crops in the basin of Mexico. In Urban Air Pollution and Forests. Ecological Studies (Analysis and Synthesis); Fenn, M.E., de Bauer, L.I., Hernández-Tejeda, T., Eds.; Springer: New York, NY, USA, 2002; Volume 156, pp. 320–333. [Google Scholar]
  34. Instituto Nacional de Estadística y Geografía (INEGI). Censo de Población y Vivienda 2010. Principales Resultados por Área Geoestadística Básica (AGEB); Instituto Nacional de Estadística y Geografía: Aguascalientes, México, 2011. [Google Scholar]
  35. Instituto Nacional de Estadística y Geografía (INEGI). Conjunto de Datos Vectoriales de Uso de Suelo y Vegetación Escala 1:250,000; Serie IV (2007–2010); Instituto Nacional de Estadística y Geografía: Aguascalientes, México, 2010. [Google Scholar]
  36. Instituto Nacional de Estadística y Geografía (INEGI). Banco de Información Económica (BIE); Instituto Nacional de Estadística y Geografía: Aguascalientes, México, 2011. [Google Scholar]
  37. Koutrakis, P.; Wolfson, J.M.; Bunyaviroch, A.; Froehlich, S.E.; Hirano, K.; Mulik, J.D. Measurement of ambient ozone using a nitrite-coated filter. Anal. Chem. 1993, 65, 209–214. [Google Scholar] [CrossRef]
  38. Sedema. Estaciones de Monitoreo. Available online: http://www.aire.df.gob.mx/default.php?opc=‘ZaBhnmI=&dc=‘ZA== (accessed on 8 March 2018).
  39. Salcedo, D.; Castro, T.; Ruiz-Suárez, L.G.; García-Reynoso, A.; Torres-Jardón, R.; Torres-Jaramillo, A.; Mar-Morales, B.E.; Salcido, A.; Celada, A.T.; Carreón-Sierra, S.; et al. Study of the regional air quality south of Mexico City (Morelos state). Sci. Total Environ. 2012, 414, 417–432. [Google Scholar] [CrossRef] [PubMed]
  40. García-Yee, J.S.; Torres-Jardón, R.; Barrera-Huertas, H.; Castro, T.; Peralta, O.; García, M.; Gutiérrez, W.; Robles, M.; Torres-Jaramillo, A.; Ortínez, A.; et al. Characterization of NOx-Ox relationships during daytime interchange of air masses over a mountain pass in the Mexico City megalopolis. Atmos. Environ. 2018, 177, 100–110. [Google Scholar] [CrossRef]
  41. Barrera-Huertas, H.; Torres, R.; Ruiz-Suárez, L.; Garcia, J.; Gutierrez, W.; Torres, A. Analysis of ozone transportation in Tlaxcala-Puebla Mexico air basin. In AGU Fall Meeting Abstracts; American Geophysical Union: Washington, DC, USA, 2014; p. 3213. [Google Scholar]
  42. Instituto Nacional de Estadística y Geografía (INEGI). Marco Geoestadístico Nacional (MGN) V.5.0.; Instituto Nacional de Estadística y Geografía: Aguascalientes, México, 2010. [Google Scholar]
  43. Hengl, T. A Practical Guide to Geostatistical Mapping of Environmental Variables; JRC Scientific and Technical Reports; Publications Office of the European Union: Luxembourg, 2007; p. 270. [Google Scholar]
  44. Li, J.; Heap, A.D. A Review of Spatial Interpolation Methods for Environmental Scientists; Geoscience Australia: Canberra, Australia, 2008; p. 137.
  45. Oficina Estatal de Información para el Desarrollo Rural Sustentable (OEIDRUS). Sistema de Información Geográfica de Morelos; Oficina Estatal de Información para el Desarrollo Rural Sustentable: Cuernavaca, Mexico, 2014. [Google Scholar]
  46. Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación (SAGARPA). Servicio de Información Agroalimentaria y Pesquera (SIAP); Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación: Ciudad de México, México, 2011. [Google Scholar]
  47. Consejo Nacional de Población (CONAPO). Índice de Marginación por Localidad; Consejo Nacional de Población: Ciudad de México, México, 2012. [Google Scholar]
  48. Grell, G.A.; Emeis, S.; Stockwell, W.R.; Schoenemeyer, T.; Forkel, R.; Michalakes, J.; Knoche, R.; Seidl, W. Application of a multiscale, coupled MM5/chemistry model to the complex terrain of the VOTALP valley campaign. Atmos. Environ. 2000, 34, 1435–1453. [Google Scholar] [CrossRef]
  49. Stockwell, W.R.; Middleton, P.; Chang, J.S.; Tang, X. The second generation regional acid deposition model chemical mechanism for regional air quality modeling. J. Geophys. Res. Atmos. 1990, 95, 16343–16367. [Google Scholar] [CrossRef]
  50. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar]
  51. Song, J.; Lei, W.; Bei, N.; Zavala, M.; De Foy, B.; Volkamer, R.; Cardenas, B.; Zheng, J.; Zhang, R.; Molina, L.T. Ozone response to emission changes: A modeling study during the MCMA-2006/MILAGRO Campaign. Atmos. Chem. Phys. 2010, 10, 3827–3846. [Google Scholar] [CrossRef] [Green Version]
  52. Lei, W.; de Foy, B.; Zavala, M.; Volkamer, R.; Molina, L.T. Characterizing ozone production in the Mexico City metropolitan area: A case study using a chemical transport model. Atmos. Chem. Phys. 2007, 7, 1347–1366. [Google Scholar] [CrossRef]
  53. Jazcilevich, A.; Garcia, A.; Ruiz-Suarez, L. A modeling study of air pollution modulation through land-use change in the valley of Mexico. Atmos. Environ. 2002, 36, 2297–2307. [Google Scholar] [CrossRef]
  54. Escalante García, J.S.; García Reynoso, J.A.; Jazcilevich Diamant, A.; Ruiz-Suárez, L.G. The influence of the Tula, Hidalgo complex on the air quality of the Mexico City Metropolitan Area. Atmósfera 2014, 27, 215–225. [Google Scholar] [CrossRef] [Green Version]
  55. Ruiz Suárez, L.; Longoria, R.; Hernandez, F.; Segura, E.; Trujillo, A.; Conde, C. Emisiones biogénicas de hidrocarburos no-metano y de oxido nítrico en la cuenca del valle de México. Atmosfera 1999, 12, 89–100. [Google Scholar]
  56. Barrera-Huertas, T.-J.R.; Ruíz-Suárez, L.G.; García-Yee, J.S.; Torres-Jaramillo, A.; Martínez, A.P.; Gutierrez, W.; García, L.M.; Robles, M.; Retama, A.; García-Reynoso, J.A. Análisis del transporte de ozono en la cuenca atmosférica de puebla-tlaxcala en el centro de México. Atmosfera 2017. submitted. [Google Scholar]
  57. Mendoza Flores, A. Análisis de la Calidad del Aire Alrededor del Complejo Industrial en Tula Durante el 2008; Universidad Nacional Autónoma de México (UNAM): Ciudad de México, México, 2011. [Google Scholar]
  58. Granada-Macias, L.M.; Rosas-Perez, R.; Torres-Jardon, R. Ozone levels at a high altitude suburban abies religiosa forest on the western slope of a mountain range under the influence of the Mexico City and Toluca Valley urban areas. Environ. Monit. Assess. 2018. submited. [Google Scholar]
  59. Molina, L.T.; Madronich, S.; Gaffney, J.S.; Apel, E.; de Foy, B.; Fast, J.; Ferrare, R.; Herndon, S.; Jimenez, J.L.; Lamb, B.; et al. An overview of the milagro 2006 campaign: Mexico city emissions and their transport and transformation. Atmos. Chem. Phys. 2010, 10, 8697–8760. [Google Scholar] [CrossRef] [Green Version]
  60. Garcia-Reynoso, A.; Jazcilevich, A.; Ruiz-Suarez, L.G.; Torres-Jardon, R.; Lastra, M.S.; Juarez, N.A.R. Ozone weekend effect analysis in Mexico City. Atmosfera 2009, 22, 281–297. [Google Scholar]
  61. Instituto Nacional de Ecología y Cambio Climático (INECC). Estudios de Calidad del Aire y su Impacto en la Región Centro de México (ECAIM); Instituto Nacional de Ecología y Cambio Climático: Ciudad de México, México, 2014; p. 744. [Google Scholar]
  62. Instituto Nacional de Ecología y Cambio Climático (INECC). Evaluación de las Capacidades Técnicas y de Infraestructura de los Laboratorios Institucionales; Instituto Nacional de Ecología y Cambio Climático: Ciudad de México, Mexico, 2015; p. 25. [Google Scholar]
  63. De Foy, B.; Krotkov, N.A.; Bei, N.; Herndon, S.C.; Huey, L.G.; Martínez, A.-P.; Ruiz-Suárez, L.G.; Word, E.C.; Zavala, M.; Molina, L.T. Hit from both sides: Tracking industrial and volcanic plumes in Mexico City with surface measurements and OMI SO2 retrievals during the MILAGRO field campaign. Atmos. Chem. Phys. 2009, 9, 9599–9617. [Google Scholar] [CrossRef]
Atmosphere 09 00223 g001 550
Figure 1. Study area with orography and boundaries between federal states and Mexico City. Four data point categories are shown.
Figure 1. Study area with orography and boundaries between federal states and Mexico City. Four data point categories are shown.
Atmosphere 09 00223 g001
Atmosphere 09 00223 g002 550
Figure 2. Conceptual model of the process used to obtain the AOT40 exceedance maps.
Figure 2. Conceptual model of the process used to obtain the AOT40 exceedance maps.
Atmosphere 09 00223 g002
Atmosphere 09 00223 g003 550
Figure 3. The average passive monitoring values at any site (Co) and the Mesoscale Climate and Chemistry Model (MCCM) average (Cm), which is obtained by averaging the profile with the hourly means from MCCM calculated over the two-month passive sampling period.
Figure 3. The average passive monitoring values at any site (Co) and the Mesoscale Climate and Chemistry Model (MCCM) average (Cm), which is obtained by averaging the profile with the hourly means from MCCM calculated over the two-month passive sampling period.
Atmosphere 09 00223 g003
Atmosphere 09 00223 g004 550
Figure 4. Evolution of the hybrid exceedance maps as observational data were gradually added to correct the maps obtained only with the air quality model. (a) Using only modeled and corrected regular grid points. (b) Panel (a) + RAMA + REMA. (c) Panel (b) + García-Yee + Barrera-Huertas. (d) Panel (c) + MILAGRO + CARIEM + Tula.
Figure 4. Evolution of the hybrid exceedance maps as observational data were gradually added to correct the maps obtained only with the air quality model. (a) Using only modeled and corrected regular grid points. (b) Panel (a) + RAMA + REMA. (c) Panel (b) + García-Yee + Barrera-Huertas. (d) Panel (c) + MILAGRO + CARIEM + Tula.
Atmosphere 09 00223 g004
Atmosphere 09 00223 g005 550
Figure 5. Maize: (a) Yield (t/ha). (b) Relative yield loss (%). (c) Production (t). (d) Crop production loss (t). (e) Production value (MXN). (f) Economic loss (MXN), for 2011.
Figure 5. Maize: (a) Yield (t/ha). (b) Relative yield loss (%). (c) Production (t). (d) Crop production loss (t). (e) Production value (MXN). (f) Economic loss (MXN), for 2011.
Atmosphere 09 00223 g005
Table
Table 1. Exposure–response functions for the four crops of economic interest analyzed in this work in the city belt of Central Mexico [12]. AOT40: accumulated exposure over the threshold of 40 ppb.
Table 1. Exposure–response functions for the four crops of economic interest analyzed in this work in the city belt of Central Mexico [12]. AOT40: accumulated exposure over the threshold of 40 ppb.
CropERFs Used for RYL Calculation, AOT40 [ppm·h]Critical Level * [ppm·h]
Maize−0.0035 * AOT4013.9
Sensitive crops (oats, sorghum, and beans)−0.01212 * AOT403
* For 5% yield reduction.
Table
Table 2. Empirical correction factors for x = c o c m showing that MCCM systematically underestimates O3 in the rural sites sampled using passive monitors, with X ¯ = 2.06   ±   1.0 .
Table 2. Empirical correction factors for x = c o c m showing that MCCM systematically underestimates O3 in the rural sites sampled using passive monitors, with X ¯ = 2.06   ±   1.0 .
Tres MaríasCICSDesierto de los LeonesLomas de AhuatlanParresSanta CatarinaSan Nicolas
Cold dry 20093.211.631.141.482.293.021.68
Cold dry 20104.201.641.131.502.054.371.52
Cold dry 20113.401.581.091.452.213.161.62
Cold dry average3.601.621.121.482.183.521.61
Warm dry 20093.241.370.761.532.023.261.29
Warm dry 20103.501.370.791.541.953.681.30
Warm dry 20113.391.370.761.511.923.581.25
Warm dry average3.371.370.771.521.963.511.28
Dry season average3.491.490.941.502.073.511.44
Std dev cold + warm dry0.360.140.190.030.150.490.19
Table
Table 3. Results for 2011 in the 20 municipalities most affected in terms of maize production, showing the input values and the results of the application of Equations (1)–(3) in the geographic information system (GIS).
Table 3. Results for 2011 in the 20 municipalities most affected in terms of maize production, showing the input values and the results of the application of Equations (1)–(3) in the geographic information system (GIS).
StateMunicipalityAOT40 (ppm·h)Relative Yield Loss (%)Yield (t/ha)Yield Loss (t/ha)Economic Loss/ha (MXN/ha)Sown Area (ha)Production Value (Millions of Pesos MXN)Economic Loss (Millions of Pesos MXN)Degree of Marginalization
MexicoVilla de Allende23.9872.500.1879913,50015211High
MexicoToluca18.2852.500.1249717,8221859Very low
MexicoAlmoloya de Juarez24.9571.340.1038221,0001078Medium
MexicoVilla Victoria25.8471.400.1149612,100766High
TlaxcalaHuamantla14.3332.400.0835214,4801565Low
MexicoTenango del Valle18.0542.600.125208891985Low
MexicoIxtlahuaca27.4080.800.0731312,903474Medium
MexicoLerma21.8361.500.093559989574Very low
MexicoSan Jose del Rincón23.2560.850.0624613,173483High
MexicoAmecameca18.3153.270.166994553663Low
MexicoTemoaya26.2871.940.165565600393Medium
MexicoZinacantepec18.9551.660.083209557603Low
MexicoTianguistenco19.7652.500.135795251563Low
MexicoAmanalco22.1661.600.104586500473Medium
GuerreroHuitzuco 13.8433.000.093228971943High
MexicoDonato Guerra22.8161.320.093947200433High
MichoacanZitacuaro23.7471.750.123597550393Medium
MorelosOcuituco23.6063.180.228383010363Medium
GuerreroIguala11.3225.500.1256043371142Low
MichoacanContepec14.5533.000.102718545692Medium
Total 194,932159085
Table
Table 4. Planted area, production value, and economic loss estimated for four crops of economic importance in the city belt of Central Mexico.
Table 4. Planted area, production value, and economic loss estimated for four crops of economic importance in the city belt of Central Mexico.
CropSown Area (ha)Production Value (In Millions of Pesos MXN)Economic Loss (In Millions of Pesos MXN)
Maize1,132,1506598191
Oats (grain)26,0879825
Oats (forage)153,1541133310
Beans61,12835751
Sorghum53,52462094

Share and Cite

MDPI and ACS Style

Ruiz-Suárez, L.G.; Mar-Morales, B.E.; García-Reynoso, J.A.; Andraca-Ayala, G.L.; Torres-Jardón, R.; García-Yee, J.S.; Barrera-Huertas, H.A.; Gavilán-García, A.; Basaldud Cruz, R. Estimation of the Impact of Ozone on Four Economically Important Crops in the City Belt of Central Mexico. Atmosphere 2018, 9, 223. https://doi.org/10.3390/atmos9060223

AMA Style

Ruiz-Suárez LG, Mar-Morales BE, García-Reynoso JA, Andraca-Ayala GL, Torres-Jardón R, García-Yee JS, Barrera-Huertas HA, Gavilán-García A, Basaldud Cruz R. Estimation of the Impact of Ozone on Four Economically Important Crops in the City Belt of Central Mexico. Atmosphere. 2018; 9(6):223. https://doi.org/10.3390/atmos9060223

Chicago/Turabian Style

Ruiz-Suárez, Luis Gerardo, Bertha Eugenia Mar-Morales, José Agustín García-Reynoso, Gema Luz Andraca-Ayala, Ricardo Torres-Jardón, José Santos García-Yee, Hugo Alberto Barrera-Huertas, Arturo Gavilán-García, and Roberto Basaldud Cruz. 2018. "Estimation of the Impact of Ozone on Four Economically Important Crops in the City Belt of Central Mexico" Atmosphere 9, no. 6: 223. https://doi.org/10.3390/atmos9060223

APA Style

Ruiz-Suárez, L. G., Mar-Morales, B. E., García-Reynoso, J. A., Andraca-Ayala, G. L., Torres-Jardón, R., García-Yee, J. S., Barrera-Huertas, H. A., Gavilán-García, A., & Basaldud Cruz, R. (2018). Estimation of the Impact of Ozone on Four Economically Important Crops in the City Belt of Central Mexico. Atmosphere, 9(6), 223. https://doi.org/10.3390/atmos9060223

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

Article Metrics

Back to TopTop