<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">GEP</journal-id><journal-title-group><journal-title>Journal of Geoscience and Environment Protection</journal-title></journal-title-group><issn pub-type="epub">2327-4336</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/gep.2020.85020</article-id><article-id pub-id-type="publisher-id">GEP-100556</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Phytotoxic Effects of Surface Ozone Exposure on Rice Crop—A Case Study of Tropical Megacity of India
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Pallavi</surname><given-names>Saxena</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Monojit</surname><given-names>Chakraborty</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Saurabh</surname><given-names>Sonwani</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Division of Climate Change and Environment, IORA Ecological Solutions, New Delhi, India</addr-line></aff><aff id="aff3"><addr-line>Department of Environmental Studies, Zakir Husain Delhi College, University of Delhi, New Delhi, India</addr-line></aff><aff id="aff1"><addr-line>Department of Environmental Sciences, Hindu College, University of Delhi, New Delhi, India</addr-line></aff><pub-date pub-type="epub"><day>06</day><month>05</month><year>2020</year></pub-date><volume>08</volume><issue>05</issue><fpage>322</fpage><lpage>334</lpage><history><date date-type="received"><day>22,</day>	<month>April</month>	<year>2020</year></date><date date-type="rev-recd"><day>26,</day>	<month>May</month>	<year>2020</year>	</date><date date-type="accepted"><day>29,</day>	<month>May</month>	<year>2020</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Increasing tropospheric ozone concentration is a big threat to food security due to its phytotoxicity. It causes a huge damage to crop production across the globe, especially in the C
  <sub>3</sub> plants (paddy (
  <em>Oryza sativa</em>)). The present study focuses on exposure-plant response index over different O
  <sub>3</sub> concentration. In this study, two metrics 
  <em>viz</em>. the average ozone for 7 h during daytime (M7) and accumulated exposure above a threshold of X ppb (AOTX) have been used in examining crop yield decline in Delhi, India. Eight AOTX indices (AOT0, AOT5, AOT10, AOT15, AOT20, AOT25, AOT30 and AOT40) were analysed and potential crop reduction was predicted. The regular monitoring of O
  <sub>3</sub> was done for 24 hours in year 2013. As per the European benchmark, a 5% yield loss was expected when AOT40 values crosses 3000 
  ppb&amp;#183h, however this study revealed that AOT40 threshold value ranged between 695 ppb to 17645 ppb which had exceeded the European benchmark in most of the months. The crop reduction was found to be ~6.3% as evaluated by AOT40 index, whereas, total AOTX contributed up to 23% of rice yield reduction in Delhi NCR. On the other side, only 2% of rice yield loss has been predicted using M7 index, which is not comparable with AOTX indices. The M7 index was also found incomparable to the calculated net yield loss (13%) for year 2013 to 2016. Hence, AOT40 may be a better index to predict the vulnerable impact of O
  <sub>3</sub> into the crop production. The total vulnerability of O
  <sub>3</sub> calculated as 57% in the crops reduction, while impacts of O
  <sub>3</sub> was calculated and summed up for both the significant and non-significant paddy growing seasons. Hence, this study highlights an alarming situation in crop yield reduction due to O
  <sub>3</sub> exposure in Delhi NCR which further threatens food security.
 
</p></abstract><kwd-group><kwd>Tropospheric Ozone</kwd><kwd> AOT40</kwd><kwd> M7</kwd><kwd> Paddy Yield</kwd><kwd> Relative Yield Loss</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Tropospheric ozone (O<sub>3</sub>) is considered as the single-most deleterious phytotoxic air pollutant and has shown a tremendous negative effect on growth and production of crops in different parts of the world (Feng et al., 2008; Mills et al., 2011; Wild et al. 2011; Stevenson et al. 2013). This is also a major concern in a country like India where high urbanization and industrialization leads to increase in precursor gases which finally results in high generation of tropospheric O<sub>3</sub> (Burney &amp; Ramanathan, 2014). This ozone is the result of reactions involved in photochemical smog formation in the troposphere. Ozone is also generated from its precursors like volatile organic compounds (VOCs), nitrogen oxides (NO<sub>x</sub>) and carbon monoxide (CO) in the presence of sunlight (Gurjar et al., 2016). Tropospheric O<sub>3</sub> is recorded as the most harmful air pollutant in North America, Japan and Europe (Sharma et al., 2017). In contrast with the past, over the last decade, O<sub>3</sub> is also found to be very high in large areas of Asian developing countries ( Mills et al., 2007; Michael et al., 2013). Ozone is a powerful oxidant which results in a number of health impacts (Rai et al., 2007; Saxena &amp; Sonwani, 2020), for instance, oxidative stress (Sinha et al., 2015; Saxena et al., 2016) and irritation in the lungs (Garg et al., 2015). In context of plants, various common crop species are sensitive to O<sub>3</sub>, therefore exposure to O<sub>3</sub> leads to decline in crop yields and resulted in lower crop quality (Debaje, 2014; Klingberg et al., 2014). Several studies have been reported so far about crop yield impacts and photosynthetic activities by O<sub>3</sub> exposure (Krupa et al., 1994; Chameides et al., 1999; Rai et al., 2007; Burney &amp; Ramanathan, 2014). It also causes an imbalance in redox status which produces early senescence. That’s why, any basis that might raise the level of plant metabolites interaction with antioxidant capacity will also be responsible for increasing the adverse impact of O<sub>3</sub> on crop plants (Fagnano et al., 2009; Saxena &amp; Sonwani, 2019). Due to the global rise in tropospheric ozone, the cost of agricultural losses is projected to reach about 35 million USD by the year 2030 (Avnery et al., 2011). The National Crop Loss Assessment Network (NCLAN) program is the nodal agency which implements permissible limit for studying the impacts of O<sub>3</sub> on crop yield using open-top chamber (OTC) experiments (Avnery et al., 2011). NCLAN results depicted that decreased yearly soybean and cotton yield were 10% and 12% respectively when seasonal mean O<sub>3</sub> mixing ratios were found to be greater than 50 ppb (Klingberg et al., 2014). Fuhrer et al. (1997) found that a 10% reduction was found in spring wheat crop when O<sub>3</sub> crosses 40 ppb concentration in Southern Europe. Different ozone indices are used to estimate the threshold effects which include AOTX (0, 10, 20, 30, 40, 50, 60), SUM06, W126, M7 or M12, POD0 or POD10. The AOT40 is a cumulative exposure index, which is obtained by the sum of differences between the hourly concentration of O<sub>3</sub> (in ppb) and permissible limit i.e. 40 ppb in each hour when its value raises above 40 ppb (UNECE, 1999). The SUM06 is the sum of hourly O<sub>3</sub> mixing ratios which goes above the permissible limit of 60 ppb all around the growing season and sunlight hours. The W126 index evaluates each hourly value by a sigmoidal weighting method so that the lower O<sub>3</sub> values are also taken into consideration. Among them, AOT40 and SUM06 are generally used to measure crop loss due to ozone exposure in Europe and USA, respectively. Lastly M7 is, calculated for the whole cropping period by the average hourly O<sub>3</sub> mixing ratio (ppb) during daytime hours between 0900 and 1600 h (Adams et al. 1989; Lesser et al. 1990; Lal et al., 2017). AOT40 is considered as the best linear fit between exposure to ozone concentration and decrease in crop yield. As per the European guidelines, if the accumulated ozone concentration under three-month duration is exceeding 3000 ppb∙h, then approximately 5% reduction in crop yield is expected (Ishi et al., 2007). Various studies have been reported so far regarding solutions in terms of the correlation between AOTX and crop responses, particularly connecting this index with the level of O<sub>3</sub>-damage to vegetation. In a developing country like India, recent studies so far have suggested a net reduction in agricultural production due to O<sub>3,</sub> both using model-simulated studies (Van Dingenen et al., 2009; Avnery et al., 2011; Ghude et al., 2014) or based on some data derived from selected sites (Debaje, 2014; Sinha et al., 2015). Very few studies have been reported on O<sub>3</sub> metrics, AOT40, AOTX indices, M7 or M12 indexes in India which can estimate the reduction in crop yield and loss in the agricultural sector. Therefore, the present study focuses on crop yield reduction of one of the staple crop of India i.e. rice by AOTX and M7 indices. In the case of rice, AOTX and M7 indices are more preferential to use as compared to other crops (Debaje, 2014; Sinha et al. 2015 and Lal et al., 2017). Hence, this study designed to evaluate and predict the loss in rice yield considering the Indian climate particularly in Delhi NCR by using AOTX and M7.</p></sec><sec id="s2"><title>2. Methodology</title><sec id="s2_1"><title>2.1. Study Area Description</title><p>Delhi, India (28.38˚N, 77.10˚E), lies in between the Indo-Gangetic Plain (IGP) in the east direction and semi-arid range of land of Rajasthan to the southwest. It constitutes an administrative area of 1483 km<sup>2</sup> and is part of the bigger National Capital Region (NCR) having a total area of 54,984 km<sup>2</sup>. Delhi also faces high traffic density in part due to high influx of immigrants from other Indian cities and daily commuters from the NCR region i.e. Ghaziabad, Gurugram and Faridabad. It is one of the most populated city, with approximately 16.7 million people and is highly polluted (vehicle population = approximately 7.4 million) especially in terms of aerosols (Tiwari et al., 2013). The city experiences three main seasons viz. summer, monsoon, and winters. Winters cover the months of November till February, summers span from March to June and monsoon season is from July to September.</p></sec><sec id="s2_2"><title>2.2. Ozone Monitoring</title><p>The ozone data were collected continuously at Wazirabad, Delhi, India (28.7336˚N, 77.2190˚E). This site is located near to Yamuna Floodplain and has extensive paddy cultivation. This site is also a heavy traffic intersection area and hence traffic density is very high especially during the peak hours. Regular monitoring was done from January to December 2013 and data was recorded on an hourly basis for 24 hours every day. Ozone concentrations were monitored using a UV based Absorption Ozone Analyser (Environnment S.A. Model O<sub>3</sub> 42M, France). This instrument is based on Beer-Lambert’s Law which measures sensitive values of ozone in ambient air.</p></sec><sec id="s2_3"><title>2.3. Crop Selection</title><p>Rice (Oryza sativa) was selected as a study crop and it was used to estimate the amount of reduction recorded in its crop yield. Rice or paddy is one of the main Kharif crops and hence, chosen for the present study. The sowing of rice starts at the end of May in Kerala i.e. particularly in the start-up of the south-west monsoon season. But, sowing dates vary accordingly with the pattern of movement of rains and hence, in Northern part of India i.e. in Delhi, Kharif season normally includes sowing during July and is harvested by end of October or November.</p></sec><sec id="s2_4"><title>2.4. Mathematical Models</title><p>AOTX is an analytical tool which is evaluated as the sum of differences between concentration, X ppb and the hourly ozone concentration, exceeding X ppb, for each sunlight hour with solar radiation higher than 50 W/m<sup>2</sup> for a period of 3 months.</p><p>The formula to calculate AOTX is:</p><p>AOTX = ∑ i = 1 n [ C O 3 − X ] i   for   C O 3 &gt; X   ppb (1)</p><p>where C<sub>O</sub><sub>3</sub> is the hourly ozone concentration in ppb, i is the running index, and n is the number of hours along with C<sub>O</sub><sub>3</sub> higher than X ppb, during the time evaluation (Grunhage et al., 1999).</p><p>M 7 = 1 n ∑ i = 0 n [ C O 3 ] (2)</p><p>From n = 0900 and 1600 h.</p><p>To obtain the suitable AOTX index in Delhi, firstly the approximate value of ozone concentration is evaluated. Due to fewer number of studies on impacts of ozone on crop production in Delhi, there are no available statistics that can be applied in this study. So, one of the methods is used to divide the permissible value for each AOTX index with total of ozone concentration. This value is supposed to be of initial crop reduction probably happen due to ozone exposure. Hence, the amount of paddy reduction in terms of percentage for each AOTX index is calculated by the following formula:</p><p>Calculated   paddy   reduction   ( % ) = AOTX   value Σ O   zone   concentration &#215; Paddy   Production (3)</p><p>where, X is the index of AOT.</p><p>As India doesn’t have indicators for estimation of decline in crop yield from ozone exposure, an analysis by AOTX indices should be carried out. As per European benchmark where 5% loss in yield for any agricultural crops will be expected to occur if the accumulated ozone concentration within that period is higher than 3000 ppb∙h, (Ishii et al., 2007), the similar critical limit was used to distinguish the variation of responses between European and Indian region. Therefore, Equation (4) is used to identify the possible crop loss that likely to occur for each AOTX index in Delhi. In this analysis, 3000 ppb∙h was used as the significant level, similar to the European concentration-based critical level (Manual, 2004).</p><p>Estimated   Crop   Reduction   ( % ) = Calculated   Crop   Reduction   ( % ) Total   AOTX ( ppb ⋅ h ) &#215; 3000   ppb ⋅ h (4)</p><p>where, 3000 ppb∙h is the European standard for the AOT40 index (UNECE, 1996). Therefore, to determine the loss of the crops for all AOTX indexes, 3000 ppb∙h is used as the standard consistent with the European critical level.</p><p>Moreover, the crop response functions for relative yields (RYs) for rice can also be calculated by the following formulas (Adams et al. 1989; Mills et al. 2007):</p><p>RY = − 0.0000039 &#215; AOT 40 + 0.94 (5)</p><p>RY = e − ( M 7 / 202 ) 2.47 / e − ( 25 / 202 ) 2.47 (6)</p><p>The relative yield loss (RYL) is calculated by</p><p>RYL = 1.0 − RY (7)</p></sec></sec><sec id="s3"><title>3. Results and Discussion</title><p>The year-long study in Delhi revealed that the tropospheric 24 h averaged O<sub>3</sub> concentrations vary from 15 to 40 ppb. However, day time-averaged concentrations touched the high notch as 54 ppb during the month of April and May. Prolonged exposure to high levels of O<sub>3</sub> may lead to damage to vegetation with possible relevant reductions in agricultural crop yields (Agrawal et al. 2003; Saxena &amp; Sonwani, 2019; Mills et al., 2007). The quantitative impacts analysis of O<sub>3</sub> has been reciprocated with different O<sub>3</sub> indexes, AOTX i.e. AOT0, AOT5, AOT10, AOT15, AOT20, AOT25, AOT30 and AOT40 in 2013 in Delhi at Wazirabad site were summarized in <xref ref-type="table" rid="table1">Table 1</xref>. These are the accumulated O<sub>3</sub> concentrations over a threshold of 0, 5, 10, 15, 25, 30 and 40 ppb respectively. The</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Summary of ozone concentration’s threshold ozone indices in Delhi during the year 2013</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Month</th><th align="center" valign="middle" >AOT0</th><th align="center" valign="middle" >AOT5</th><th align="center" valign="middle" >AOT10</th><th align="center" valign="middle" >AOT15</th><th align="center" valign="middle" >AOT20</th><th align="center" valign="middle" >AOT25</th><th align="center" valign="middle" >AOT30</th><th align="center" valign="middle" >AOT40</th><th align="center" valign="middle" >M7</th></tr></thead><tr><td align="center" valign="middle" >Jan-Mar</td><td align="center" valign="middle" >57,654</td><td align="center" valign="middle" >38,765</td><td align="center" valign="middle" >28,765</td><td align="center" valign="middle" >23,176</td><td align="center" valign="middle" >18,769</td><td align="center" valign="middle" >14,547</td><td align="center" valign="middle" >10,876</td><td align="center" valign="middle" >7685</td><td align="center" valign="middle" >31.02</td></tr><tr><td align="center" valign="middle" >Apr-Jun</td><td align="center" valign="middle" >81,497</td><td align="center" valign="middle" >57,232</td><td align="center" valign="middle" >50,674</td><td align="center" valign="middle" >43,969</td><td align="center" valign="middle" >36,946</td><td align="center" valign="middle" >30,024</td><td align="center" valign="middle" >23,078</td><td align="center" valign="middle" >17,645</td><td align="center" valign="middle" >49.26</td></tr><tr><td align="center" valign="middle" >July-Sept</td><td align="center" valign="middle" >38,065</td><td align="center" valign="middle" >23,958</td><td align="center" valign="middle" >18,428</td><td align="center" valign="middle" >13,537</td><td align="center" valign="middle" >9341</td><td align="center" valign="middle" >5983</td><td align="center" valign="middle" >3386</td><td align="center" valign="middle" >695</td><td align="center" valign="middle" >28.63</td></tr><tr><td align="center" valign="middle" >Oct-Dec</td><td align="center" valign="middle" >61,366</td><td align="center" valign="middle" >35,454</td><td align="center" valign="middle" >29,485</td><td align="center" valign="middle" >23,640</td><td align="center" valign="middle" >18,079</td><td align="center" valign="middle" >13,203</td><td align="center" valign="middle" >9070</td><td align="center" valign="middle" >3369</td><td align="center" valign="middle" >35.20</td></tr></tbody></table></table-wrap><p>various threshold values which were produced in <xref ref-type="table" rid="table1">Table 1</xref> were important to depict the range of threshold values at Wazirabad, Delhi, which will consequently be used to estimate the crucial level of the AOTX index. Different AOTX index give different threshold values and hence they are important to analyze for the same. In <xref ref-type="table" rid="table1">Table 1</xref>, it has been shown that the AOT40 threshold value ranged between 695 ppb∙h. to 17645 ppb∙h. As per the European standard, permissible limit which is higher than 3000 ppb∙h covering a 3 months period is expected to result in an approximate 5% crop reduction. As per the estimated range of ozone limit concentration of 40 ppb in Delhi, it is largely affected by the ozone concentration of 5% - 30% due to the increase in the concentration of accumulated ozone exposure i.e. up to 18,000 ppb∙h.</p><p>On the contrary, another exposure plant response index, M7 has shown values as high 63 ppb in May and as low as 25 ppb in July. The M7 value during the Kharif season was calculated to be 32 ppb, almost similar value as reported by Lal et al., 2017 for the north Indian region (<xref ref-type="table" rid="table1">Table 1</xref>). It has been observed that the statistical correlation between 24-h average monthly O<sub>3</sub> and monthly average M7 as well as correlations between monthly average AOT40 and monthly average M7 are very poorly fitted (R<sup>2</sup> = 0.47 &amp; R<sup>2</sup> = 0.27, respectively). However, the relations among day time average O<sub>3</sub> concentration, monthly average M7, and monthly average AOT40 have been found slightly better (R<sup>2</sup> = 0.67 &amp; R<sup>2</sup> = 0.45, respectively). This has been observed that the M7 index value was lower i.e. 32 ppb during the paddy season/significant season (July to November) than the value i.e. 38 ppb, was found during the non-paddy season/non-significant season (December to June). Moreover, during paddy season, the estimated crop reduction calculated by M7 was found to be 2% while in non-paddy season the total crop reduction was calculated as 11.2%.</p><p>The analysis on the effect of AOTX indexes in paddy reduction was performed in both significant and non-significant seasons. <xref ref-type="fig" rid="fig1">Figure 1</xref>(a) and <xref ref-type="fig" rid="fig1">Figure 1</xref>(b) shows the correlation between AOTX indexes and paddy reduction as per the Indian scenario of climate for two periods in the year 2013. In <xref ref-type="fig" rid="fig1">Figure 1</xref>(a) and <xref ref-type="fig" rid="fig1">Figure 1</xref>(b), there are two subsequent patterns between the AOTX indexes and paddy reduction in the Delhi (Indian) climate. As per observations, during both the periods, AOT0 gave the least paddy reduction while AOT40 reported the maximum loss in paddy. The highest percentage for paddy reduction was found to be 76.45%, depicted by the AOT40 index during the non-significant</p><p>season. In India especially in Northern India (one of them in Delhi NCR) where paddy production was mostly influenced by ozone exposure during the non-significant season (<xref ref-type="fig" rid="fig1">Figure 1</xref>(b)).</p><p>During this season, Delhi NCR experiences lower boundary layer, low-temperature inversions, high relative humidity and low wind speed with less rainfall as compared to other months (Lal et al., 2017; Saxena &amp; Kulshrestha, 2016). These conditions are favorable for high ozone concentrations during winters while in early summers high temperature and high sunlight intensity favor high ozone concentrations (Bhuyan et al., 2014). The most favorable season for sowing of rice is during monsoon season. This season acts as a scavenging agent which cleans up the pollutants present in the atmosphere due to heavy rainfall, high relative humidity, and high wind velocity. This is the main reason that ozone concentrations during this season are very low (Ghosh et al., 2013).</p><p>Although this occurrence will decline the rate of ozone conversions but higher humidity and moderate temperatures can give rise to higher stomatal conductance. Consequently, the ozone uptake by plant leaves results in a loss of yield (Mohammed et al., 2013). This is a clear indication that paddy crop reduction had occurred during the significant season, though the percentage reduction was less than in the non-significant season.</p><p>Statistical correlation associations between paddy reduction and AOTX indexes in 2013 are shown in <xref ref-type="table" rid="table2">Table 2</xref>. As observed, the AOT40 index has got the highest R<sup>2</sup> value as compared to the other indexes. Mohammed et al. (2013), had compared several AOTX indexes i.e. AOT15, AOT20, AOT25, AOT30, and AOT50 in Malaysia to determine the impact of different indexes on paddy. They found that AOT50 index had provided the best fit correlation between the index and crop i.e. paddy response. Hence, as per the variation of different indexes, the crop response would vary. Apart from this study, several studies find that AOT40 index has the best fit when considering the crop response (Rai et al., 2007; Beig et al., 2008; Debaje, 2014; Ghude et al., 2014; Sinha et al., 2015). A similar observation has been obtained in our present study which compare all the indices with paddy response, with AOT40 having the best paddy reduction response. For example, 16.3% of relative yield loss of paddy was reported by Ghude et al. (2008) in Delhi region by using AOT40 index. Another study was reported by Lal et al., (2017), where 0.3% and 6.3% of paddy reduction was found in total of different parts of India by using M7 and AOT40 respectively. This present study also reported that AOT40 is having the best paddy reduction response of approximately 6.3% in Delhi region.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref> depicts the approximate percentage of paddy reduction for each AOTX index in Delhi climate (Indian tropical). This has been widely accepted if the cumulative tropospheric O<sub>3</sub> concentration will exceed the limit of 3000 ppb∙h in a particular crop growing season and only then conspicuous yield loss can be seen (K&#228;renlampi &amp; Sk&#228;rby, 1996; Beig et al., 2008). During our research, the</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> R<sup>2</sup> value for the Relationship of AOTX and M7 indices with Paddy Reduction in Wazirabad, Delhi</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >AOTX indexes</th><th align="center" valign="middle" >Coefficient of determination (R<sup>2</sup>)</th><th align="center" valign="middle" >Equation (y = paddy reduction)</th></tr></thead><tr><td align="center" valign="middle" >AOT0</td><td align="center" valign="middle" >0.78</td><td align="center" valign="middle" >y = 13,612 + 0.7633x</td></tr><tr><td align="center" valign="middle" >AOT5</td><td align="center" valign="middle" >0.80</td><td align="center" valign="middle" >y = 14,671 + 1.0212x</td></tr><tr><td align="center" valign="middle" >AOT10</td><td align="center" valign="middle" >0.81</td><td align="center" valign="middle" >y = 12,243 + 1.663x</td></tr><tr><td align="center" valign="middle" >AOT15</td><td align="center" valign="middle" >0.86</td><td align="center" valign="middle" >y = 10,133 + 1.732x</td></tr><tr><td align="center" valign="middle" >AOT20</td><td align="center" valign="middle" >0.88</td><td align="center" valign="middle" >y = 8613.3 + 1.967x</td></tr><tr><td align="center" valign="middle" >AOT25</td><td align="center" valign="middle" >0.92</td><td align="center" valign="middle" >y = 8412.4 + 2.167x</td></tr><tr><td align="center" valign="middle" >AOT30</td><td align="center" valign="middle" >0.93</td><td align="center" valign="middle" >y = 8316.6 + 3.22x</td></tr><tr><td align="center" valign="middle" >AOT40</td><td align="center" valign="middle" >0.93</td><td align="center" valign="middle" >y = 11,671 + 5.13x</td></tr><tr><td align="center" valign="middle" >M7</td><td align="center" valign="middle" >0.46</td><td align="center" valign="middle" >y = 3287.7x − 90194</td></tr></tbody></table></table-wrap><p>cumulative O<sub>3</sub> concentration was found to be more than 14,000 ppbh, which is far exceeded value than the critical level. Therefore, reduction in paddy was calculated as 6.3% by AOT40, which is much higher value than the AOT40 calculated by the critical value as 3000 ppb∙h. <xref ref-type="table" rid="table2">Table 2</xref> results also support this fact, as the AOT40 index fits best with the Indian climate as compared to other indexes. The obtained percentage reduction i.e. 6.3% is different from the European benchmark which is only 5%. Thus, in the Indian climate, especially in an urban area like Delhi NCR, there is a high percentage crop (paddy) reduction and it is different from the European benchmark.</p><p>Relative Yield Loss (RYL) of rice has also been calculated as 0.23 of RYL in the Delhi NCR region at Wazirabad site which is a quite comparable to the value of 0.29 as reported in Lal et al., (2017). Earlier studies reported in India on RYL of rice had also reported the similar values (Debaje, 2014; Sinha et al., 2015; Ghude et al., 2014; Van Dingenen et al., 2009). This value is alarming for agricultural scientists as ozone exposure in Delhi NCR is much high and is highly affecting the crop production in Delhi and its nearby areas. The total yield loss of rice has been estimated by using the AOT40 index and is found to be 0.16 Mt (6%) of 2.6 Mt rice produced in Delhi NCR (<xref ref-type="table" rid="table3">Table 3</xref>). In comparison, the total AOTX contributed to a huge rice yield loss of 0.61 Mt (23%). On the contrary, M7 index contributed to the rice yield loss of only 0.06 Mt (2%). Therefore, the study has shown a clear impact of O<sub>3</sub> phytotoxicity on rice yield. The indices used to calculate the entire range of crop yield losses, shows that 2% to 23% losses of crop yield occur in Delhi NCR. This has been validated with the rice production data from ICAR (2018) which finds a decrease of up to 13% from 2013 to 2016 in the Delhi NCR. Similar study was conducted by Ghude et al. (2008), where, 16.3% of relative yield loss of paddy was reported in Delhi region by using AOT40 index.</p><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Loss of paddy by relative yield loss with respect to AOT40 and M7 in Delhi-NCR of 2013 (Source: ICAR reporter, 2018)</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >National Capital Region (districts)</th><th align="center" valign="middle"  rowspan="2"  >Avg O<sub>3</sub> (ppb)</th><th align="center" valign="middle"  rowspan="2"  >Avg AOT40 (ppb∙h)</th><th align="center" valign="middle"  rowspan="2"  >Avg M7 (ppb)</th><th align="center" valign="middle"  rowspan="2"  >RYL AOT40</th><th align="center" valign="middle"  rowspan="2"  >RYL M7</th><th align="center" valign="middle"  rowspan="2"  >Paddy Yield (t)</th><th align="center" valign="middle"  colspan="2"  >Paddy loss (t)</th></tr></thead><tr><td align="center" valign="middle" >AOT40</td><td align="center" valign="middle" >M7</td></tr><tr><td align="center" valign="middle" >Delhi (7)</td><td align="center" valign="middle"  rowspan="4"  >19</td><td align="center" valign="middle"  rowspan="4"  >318</td><td align="center" valign="middle"  rowspan="4"  >32</td><td align="center" valign="middle"  rowspan="4"  >0.06</td><td align="center" valign="middle"  rowspan="4"  >0.02</td><td align="center" valign="middle" >29,652</td><td align="center" valign="middle" >1779</td><td align="center" valign="middle" >593</td></tr><tr><td align="center" valign="middle" >Haryana (10)</td><td align="center" valign="middle" >2,153,000</td><td align="center" valign="middle" >129,180</td><td align="center" valign="middle" >43,060</td></tr><tr><td align="center" valign="middle" >Uttar Pradesh (7)</td><td align="center" valign="middle" >452,250</td><td align="center" valign="middle" >27,135</td><td align="center" valign="middle" >9045</td></tr><tr><td align="center" valign="middle" >Rajasthan (2)</td><td align="center" valign="middle" >4259</td><td align="center" valign="middle" >256</td><td align="center" valign="middle" >85</td></tr><tr><td align="center" valign="middle" >Total</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >2,639,161</td><td align="center" valign="middle" >158,350</td><td align="center" valign="middle" >52,783</td></tr></tbody></table></table-wrap></sec><sec id="s4"><title>4. Conclusion</title><p>The present study highlighted that paddy crop reduction was higher during the non-significant season as compared to a significant season. This study also demonstrates that among all the selected AOTX and M7 indexes, AOT40 is highly responsible for paddy reduction as compared to the rest. To further support this statement, linear regression analysis confirms that AOT40 best fits in the context of Indian climate particularly Delhi. Thus, this study suggests that the AOT40 is the most appropriate index especially in comparison to M7 while evaluating potential ozone effects on crops in the tropical regions. The marked distinction in the values of AOTX indices in Europe and India is due to the variation in climatic conditions and availability of longer daylight hours in India. In addition to that, it has also been reported that AOT40 index evaluated 6.3% of paddy reduction which was higher than 3000 ppb∙h over three months of the growing season. This result coincided with the European guidelines, which indicates that at least 5% crop reduction would be expected. By taking the overall view of crop loss due to ozone exposure in India, some changes can be made in the modification of AOTX indexes. If AOT40 is used as the index, this study indicates that the approximate crop reduction must be shifted from 5% to 6.3%. Moreover, it has been expected that as per the AOT40 index the total yield loss of rice would be 0.16 Mt (6%) of 2.6 Mt rice produced in Delhi NCR. Whereas, the total AOTX contributed huge rice yield loss as 0.61 Mt i.e. 23% of the total rice produced. On the contrary, M7 index contributed to the rice yield loss as 0.06 Mt (2%), which may does not truly fit for Delhi NCR or the tropical climate and the net crop yield loss calculated as 13% during the year 2013 to 2016. Therefore, the study has shown a clear impact of O<sub>3</sub> phytotoxicity on rice yield over Delhi NCR. As paddy is the main crop of India, a decreasing paddy production in the future may negatively impact the agriculture sector, particularly in Delhi.</p></sec><sec id="s5"><title>Acknowledgements</title><p>The authors pay sincere thanks to National Environmental Engineering Research Institute (NEERI), New Delhi for providing the instrumentation facility.</p></sec><sec id="s6"><title>Conflicts of Interest</title><p>The authors declare that they have no conflict of interest.</p></sec><sec id="s7"><title>Cite this paper</title><p>Saxena, P., Chakraborty, M., &amp; Sonwani, S. (2020). Phytotoxic Effects of Surface Ozone Exposure on Rice Crop—A Case Study of Tropical Me- gacity of India. Journal of Geoscience and Environment Protection, 8, 322-334. https://doi.org/10.4236/gep.2020.85020</p></sec></body><back><ref-list><title>References</title><ref id="scirp.100556-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Adams, R. M., Glyer, J. D., Johnson, S. L., &amp; McCarl, B. A. (1989). A Reassessment of the Economic Effects of Ozone on US Agriculture. JAPCA, 39, 960-968. https://doi.org/10.1080/08940630.1989.10466583</mixed-citation></ref><ref id="scirp.100556-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Agrawal, M., Singh, B., Rajput, M., Marshall, F., &amp; Bell, J. N. B. (2003). Effect of Air Pollution on Peri-Urban Agriculture: A Case Study. Environmental Pollution, 126, 323-329. https://doi.org/10.1016/S0269-7491(03)00245-8</mixed-citation></ref><ref id="scirp.100556-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Avnery, S., Mauzerall, D. L., Liu, J., &amp; Horowitz, L. W. (2011). Global Crop Yield Reductions Due to Surface Ozone Exposure: 1. Year 2000 Crop Production Losses and Economic Damage. Atmospheric Environment, 45, 2284-2296. https://doi.org/10.1016/j.atmosenv.2010.11.045</mixed-citation></ref><ref id="scirp.100556-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Beig, G., Ghude, S. D., Polade, S. D., &amp; Tyagi, B. (2008). Threshold Exceedances and Cumulative Ozone Exposure Indices at Tropical Suburban Site. Geophysical Research Letters, 35, L02802. https://doi.org/10.1029/2007GL031434</mixed-citation></ref><ref id="scirp.100556-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Bhuyan, P. K., Bharali, C., Pathak, B., &amp; Kalita, G. (2014). The Role of Precursor Gases and Meteorology on Temporal Evolution of O3 at a Tropical Location in Northeast India. Environmental Science and Pollution Research, 21, 6696-6713. https://doi.org/10.1007/s11356-014-2587-3</mixed-citation></ref><ref id="scirp.100556-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Burney, J., &amp; Ramanathan, V. (2014). Recent Climate and Air Pollution Impacts on Indian Agriculture. Proceedings of the National Academy of Sciences, 111, 16319-16324. https://doi.org/10.1073/pnas.1317275111</mixed-citation></ref><ref id="scirp.100556-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Chameides, W. L., Xingsheng, L., Xiaoyan, T., Xiuji, Z., Luo, C., Kiang, C. S., et al. (1999). Is Ozone Pollution Affecting Crop Yields in China? Geophysical Research Letters, 26, 867-870. https://doi.org/10.1029/1999GL900068</mixed-citation></ref><ref id="scirp.100556-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Debaje, S. B. (2014). Estimated Crop Yield Losses Due to Surface Ozone Exposure and Economic Damage in India. Environmental Science and Pollution Research, 21, 7329- 7338. https://doi.org/10.1007/s11356-014-2657-6</mixed-citation></ref><ref id="scirp.100556-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Fagnano, M., Maggio, A., &amp; Fumagalli, I. (2009). Crops’ Responses to Ozone in Mediterranean Environments. Environmental Pollution, 157, 1438-1444. https://doi.org/10.1016/j.envpol.2008.09.001</mixed-citation></ref><ref id="scirp.100556-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Feng, Z., Kobayashi, K., &amp; Ainsworth, E. A. (2008). Impact of Elevated Ozone Concentration on Growth, Physiology, and Yield of Wheat (Triticum aestivum L.): A Meta-Analysis. Global Change Biology, 14, 2696-2708. https://doi.org/10.1111/j.1365-2486.2008.01673.x</mixed-citation></ref><ref id="scirp.100556-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Fuhrer, J., Sk&amp;#228rby, L., &amp; Ashmore, M. R. (1997). Critical Levels for Ozone Effects on Vegetation in Europe. Environmental Pollution, 97, 91-106. https://doi.org/10.1016/S0269-7491(97)00067-5</mixed-citation></ref><ref id="scirp.100556-ref12"><label>12</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Garg</surname><given-names> A.</given-names></name>,<name name-style="western"><surname> Saxena</surname><given-names> P.</given-names></name>,<name name-style="western"><surname> &amp; Ghosh</surname><given-names> C. </given-names></name>,<etal>et al</etal>. (<year>2015</year>)<article-title>. Evaluation of Tolerance and Sensitivity of Selected Plant Species with Special Reference to Gasoline Exhaust Pollution</article-title><source> International Journal of Scientific &amp; Technology Research</source><volume> 4</volume>,<fpage> 199</fpage>-<lpage>207</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.100556-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Ghosh, D., Lal, S., &amp; Sarkar, U. (2013). High Nocturnal Ozone Levels at a Surface Site in Kolkata, India: Trade-Off between Meteorology and Specific Nocturnal Chemistry. Urban Climate, 5, 82-103. https://doi.org/10.1016/j.uclim.2013.07.002</mixed-citation></ref><ref id="scirp.100556-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Ghude, S. D., Jain, S. L., Arya, B. C., Beig, G., Ahammed, Y. N., Kumar, A., &amp; Tyagi, B. (2008). Ozone in Ambient Air at a Tropical Megacity, Delhi: Characteristics, Trends and Cumulative Ozone Exposure Indices. Journal of Atmospheric Chemistry, 60, 237- 252. https://doi.org/10.1007/s10874-009-9119-4</mixed-citation></ref><ref id="scirp.100556-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Ghude, S. D., Jena, C., Chate, D. M., Beig, G., Pfister, G. G., Kumar, R., &amp; Ramanathan, V. (2014). Reductions in India’s Crop Yield Due to Ozone. Geophysical Research Letters, 41, 5685-5691. https://doi.org/10.1002/2014GL060930</mixed-citation></ref><ref id="scirp.100556-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Grünhage, L., J&amp;#228ger, H. J., Haenel, H. D., L&amp;#246pmeier, F. J., &amp; Hanewald, K. (1999). The European Critical Levels for Ozone: Improving Their Usage. Environmental Pollution, 105, 163-173. https://doi.org/10.1016/S0269-7491(99)00029-9</mixed-citation></ref><ref id="scirp.100556-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Gurjar, B. R., Ravindra, K., &amp; Nagpure, A. S. (2016). Air Pollution Trends over Indian Megacities and Their Local-to-Global Implications. Atmospheric Environment, 142, 475-495. https://doi.org/10.1016/j.atmosenv.2016.06.030</mixed-citation></ref><ref id="scirp.100556-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">ICAR Reporter (2018). https://icar.org.in/sites/default/files/IR-April-June-2018Web.pdf</mixed-citation></ref><ref id="scirp.100556-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Ishii, S., Bell, J. N. B., &amp; Marshall, F. M. (2007). Phytotoxic Risk Assessment of Ambient Air Pollution on Agricultural Crops in Selangor State, Malaysia. Environmental Pollution, 150, 267-279. https://doi.org/10.1016/j.envpol.2007.01.012</mixed-citation></ref><ref id="scirp.100556-ref20"><label>20</label><mixed-citation publication-type="book" xlink:type="simple">K&amp;#228renlampi, L., &amp; Sk&amp;#228rby, L. (Eds.) (1996). Critical Levels for Ozone in Europe: Testing and Finalising the Concepts. Kuopio: University of Kuopio.</mixed-citation></ref><ref id="scirp.100556-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Klingberg, J., Engardt, M., Karlsson, P. E., Langner, J., &amp; Pleijel, H. (2014). Declining Ozone Exposure of European Vegetation under Climate Change and Reduced Precursor Emissions. Biogeosciences, 11, 5269-5283. https://doi.org/10.5194/bg-11-5269-2014</mixed-citation></ref><ref id="scirp.100556-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Krupa, S. V., Nosal, M., &amp; Legge, A. H. (1994). Ambient Ozone and Crop Loss: Establishing a Cause-Effect Relationship. Environmental Pollution, 83, 269-276. https://doi.org/10.1016/0269-7491(94)90147-3</mixed-citation></ref><ref id="scirp.100556-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Lal, S., Venkataramani, S., Naja, M., Kuniyal, J. C., Mandal, T. K., Bhuyan, P. K., et al. (2017). Loss of Crop Yields in India Due to Surface Ozone: An Estimation Based on a Network of Observations. Environmental Science and Pollution Research, 24, 20972- 20981. https://doi.org/10.1007/s11356-017-9729-3</mixed-citation></ref><ref id="scirp.100556-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Lesser, V. M., Rawlings, J. O., Spruill, S. E., &amp; Somerville, M. C. (1990). Ozone Effects on Agricultural Crops: Statistical Methodologies and Estimated Dose-Response Relationships. Crop Science, 30, 148-155. https://doi.org/10.2135/cropsci1990.0011183X003000010033x</mixed-citation></ref><ref id="scirp.100556-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Manual, M. (2004). Manual on Methodologies and Criteria for Modelling and Mapping Critical Loads and Levels and Air Pollution Effects, Risks and Trends. In: Mapping Manual UNECE Convention on Long-Range Transboundary Air Pollution (LRTAP) (251 p). Geneva: UNECE.</mixed-citation></ref><ref id="scirp.100556-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Mills, G., Buse, A., Gimeno, B., Bermejo, V., Holland, M., Emberson, L., &amp; Pleijel, H. (2007). A Synthesis of AOT40-Based Response Functions and Critical Levels of Ozone for Agricultural and Horticultural Crops. Atmospheric Environment, 41, 2630-2643. https://doi.org/10.1016/j.atmosenv.2006.11.016</mixed-citation></ref><ref id="scirp.100556-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Mills, G., Hayes, F., Simpson, D., Emberson, L., Norris, D., Harmens, H., &amp; Büker, P. (2011). Evidence of Widespread Effects of Ozone on Crops and (Semi-) Natural Vegetation in Europe (1990-2006) in Relation to AOT40- and Flux-Based Risk Maps. Global Change Biology, 17, 592-613. https://doi.org/10.1111/j.1365-2486.2010.02217.x</mixed-citation></ref><ref id="scirp.100556-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Mohammed, N. I., Ramli, N. A., &amp; Yahya, A. S. (2013). Ozone Phytotoxicity Evaluation and Prediction of Crops Production in Tropical Regions. Atmospheric Environment, 68, 343-349. https://doi.org/10.1016/j.atmosenv.2012.09.010</mixed-citation></ref><ref id="scirp.100556-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Rai, R., Agrawal, M., &amp; Agrawal, S. B. (2007). Assessment of Yield Losses in Tropical Wheat Using Open Top Chambers. Atmospheric Environment, 41, 9543-9554. https://doi.org/10.1016/j.atmosenv.2007.08.038</mixed-citation></ref><ref id="scirp.100556-ref30"><label>30</label><mixed-citation publication-type="book" xlink:type="simple">Saxena, P. Sonwani, S. &amp; Kulshrestha, U.C. (2016). Impact of Tropospheric Ozone and Particulate Matter on Plant Health. In D. Pandey, &amp; A. Sarkar (Eds.), Sustaining Future Food Security in Changing Environments (pp. 19-60). UK: Nova Science Publisher.</mixed-citation></ref><ref id="scirp.100556-ref31"><label>31</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Saxena</surname><given-names> P.</given-names></name>,<name name-style="western"><surname> &amp; Kulshrestha</surname><given-names> U.C. </given-names></name>,<etal>et al</etal>. (<year>2016</year>)<article-title>. Contribution of Carbonaceous Species in SOA Formation during Fog and Non-Fog Period</article-title><source> Journal of Chemical Biological Physical Science</source><volume> 7</volume>,<fpage> 31</fpage>-<lpage>48</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.100556-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Saxena, P., &amp; Sonwani, S. (2019). Criteria Air Pollutants and Their Impact on Environmental Health. Singapore: Springer. https://doi.org/10.1007/978-981-13-9992-3</mixed-citation></ref><ref id="scirp.100556-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Saxena, P., &amp; Sonwani, S. (2020). Remediation of Ozone Pollution by Ornamental Plants in Indoor Environment. Global Journal of Environmental Science and Management, 6, No. 4. (In Press)</mixed-citation></ref><ref id="scirp.100556-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Sharma, S., &amp; Chatterjee, S. (2017). Microplastic Pollution, a Threat to Marine Ecosystem and Human Health: A Short Review. Environmental Science and Pollution Research, 24, 21530-21547. https://doi.org/10.1007/s11356-017-9910-8</mixed-citation></ref><ref id="scirp.100556-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Sinha, B., Sangwan, K. S., Maurya, Y., Kumar, V., Sarkar, C., Chandra, B. P., &amp; Sinha, V. (2015). Assessment of Crop Yield Losses in Punjab and Haryana Using Two Years of Continuous In-Situ Ozone Measurements. Atmospheric Chemistry &amp; Physics Discussions, 15. https://doi.org/10.5194/acpd-15-2355-2015</mixed-citation></ref><ref id="scirp.100556-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Stevenson, D. S., Young, P. J., Naik, V., Lamarque, J. F., Shindell, D. T., Voulgarakis, A., et al. (2013). Tropospheric Ozone Changes, Radiative Forcing and Attribution to Emissions in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP). Atmospheric Chemistry &amp; Physics, 13, 3063-3085. https://doi.org/10.5194/acp-13-3063-2013</mixed-citation></ref><ref id="scirp.100556-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Tiwari, S., Srivastava, A. K., Bisht, D. S., Parmita, P., Srivastava, M. K., &amp; Attri, S. D. (2013). Diurnal and Seasonal Variations of Black Carbon and PM2.5 over New Delhi, India: Influence of Meteorology. Atmospheric Research, 125, 50-62. https://doi.org/10.1016/j.atmosres.2013.01.011</mixed-citation></ref><ref id="scirp.100556-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">UNECE (1996). Manual on Methodologies and Criteria for Mapping Critical Levels/ Loads and Geographical Areas Where They Are Exceeded. Berlin: Umweltbundesamt.</mixed-citation></ref><ref id="scirp.100556-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">UNECE (1999). Protocol to the 1979 Convention on Long-Range Transboundary Air Pollution to Abate Acidification, Eutrophication and Ground-Level Ozone. Gothenburg, Thirtieth Day of November. New York and Geneva: United Nations Treaty Collection.</mixed-citation></ref><ref id="scirp.100556-ref40"><label>40</label><mixed-citation publication-type="other" xlink:type="simple">Van Dingenen, R., Dentener, F. J., Raes, F., Krol, M. C., Emberson, L., &amp; Cofala, J. (2009). The Global Impact of Ozone on Agricultural Crop Yields under Current and Future Air Quality Legislation. Atmospheric Environment, 43, 604-618. https://doi.org/10.1016/j.atmosenv.2008.10.033</mixed-citation></ref><ref id="scirp.100556-ref41"><label>41</label><mixed-citation publication-type="other" xlink:type="simple">Wild, O., Fiore, A. M., Shindell, D. T., Doherty, R. M., Collins, W. J., Dentener, F. J., et al. (2011). Modelling Future Changes in Surface Ozone: A Parameterized Approach. Atmospheric Chemistry &amp; Physics Discussions, 12, 2037-2054. https://doi.org/10.5194/acpd-11-27547-2011</mixed-citation></ref></ref-list></back></article>