Effects of Ambient Ozone on Soybean Biophysical Variables and Mineral Nutrient Accumulation

The effects of increasing ambient ozone (O3) concentrations on food security has become a major concern as the demand for agricultural productivity is projected to increase significantly over the next several decades. In this contribution, the responses of common soybean genotypes (AK-HARROW, PI88788, DWIGHT, PANA, and WILLIAMS82) to ambient O3 are characterized using hyperspectral data and foliar biophysical, mineral nutrient concentrations and soybean yield. Specifically, leaf reflectance spectra measured at different growth stages and canopy layers were used to examine the spectral indices that were most strongly correlated with leaf physiological status. The effects of elevated O3 on six important nutrients (K, Ca, Mg, Fe, Mn and Cu) were evaluated by analyzing the variations in nutrient concentrations at two critical growth stages with increasing ambient O3 concentration using Partial Least Square Regression (PLSR). Lastly, the identified best spectral indices and the robust nutrient prediction models were extrapolated to the entire growth period to explore their ability to track the effects of ambient O3 concentrations on soybean physiology and nutrient uptake. The results showed that fluorescence yield (∆F/Fm’) and photochemical quenching (qP) appear to be good indicators of soybean physiological responses to O3 stress that are echoed by the harvest index (HI). Newly identified normalized difference spectral index (NDSI) [R416, R2371] always had the highest correlation (R2 > 0.6) with ∆F/Fm’, qP and electron transport rate (ETR, μmol m−2 s−1) compared to the published indices. Additionally, there were significant and broad spectral regions in visible and near infrared region that were well-correlated with ∆F/Fm’ and selected NDSIs that were applicable to satellite observations. The results of nutrient modeling using PLSR explained 54–87% of the variance in nutrient concentrations, and the predicted mineral nutrient accumulation throughout the growing season reflected the responses of ozone tolerant and sensitive genotypes well. NDSI [R416, R2371] demonstrated great potential in regard to its sensitivity in tracking plant physiological responses to changing ambient O3 concentrations. The outcome of this research has potential implications for development of space-based observation of large-scale crop responses to O3 damage, as well as for biotechnological breeding efforts to improve ozone tolerance under future climate scenarios.


Introduction
Lack of water availability as a result of elevated temperature and increases in ambient ozone due to climate changes have been identified as major threats to global food security [1,2].This development is of great socio-economic importance given that global demand for agricultural production is projected to increase by 50% by 2050 [3].The early detection of abiotic stress is critical to avoid yield losses in crops, which can be affected by water stress and/or photosystem damage from ambient ozone [4,5].
Excluding corn, soybean is the most-planted field crop in the US and accounts for over 50% of the world's production [6].However, the yearly cost to the soybean crops due to the damage from ambient ozone concentrations likely exceeds $1 billion [1] in the US and may be as high as tens of billions of dollars globally [7].These significant costs to world agriculture are predicted to become even more costly in the next several decades with the predicted increase in tropospheric ozone [8].
As a rapid, non-destructive phenotyping technique, remote sensing has been demonstrated to be effective in monitoring plant growth inferred by leaf, canopy and ecosystem scale biophysical traits, including chlorophyll content [9], leaf area index (LAI) [10], terrestrial photosynthesis and radiation use efficiency [11], and sun-induced chlorophyll fluorescence [12].Remote sensing-derived vegetation water content has also been widely used in mapping plant health [13][14][15].Ozone damage can cause noticeable changes in one or more of these plant health indicators [16].Even though foliar damage from ambient ozone concentrations can be visually observed [17], changes in the plant biophysical variables can ultimately decrease crop productivity, especially at the early stage of exposure to ozone, and these changes are not readily visible to the human eye.Cost-effective and timely information regarding pre-visual ozone stress on crops is critical for monitoring spatio-temporal dynamics of crop growth to guide farm management practices and crop breeding efforts.
Canopy temperature has long been used to monitor crop water stress [18,19], which is based on the fact that water deficit induces a decrease in leaf stomatal conductance and transpiration rate and consequently, increases the leaf/canopy temperature.However, the relationship between the canopy temperature and stress levels is not straightforward such that the level of stomatal conductance interacts with evaporative demand as well as vegetation water content [20].Additionally, remotely sensed canopy temperature is confounded with shadow and background soil temperature variation [21], which makes plant health detection particularly challenging in heterogeneous, open tree canopy environments.In recent years, a growing body of studies has demonstrated that chlorophyll fluorescence and the photochemical reflectance index (PRI) can be used as early indicators of plant health [20,22].
Chlorophyll fluorescence is closely related to the function of photosystem II, stomatal conductance and other photosystem de-excitation pathways [23], providing a direct approach for the early detection of vegetation stress associated with photosynthetic functionality, before chlorophyll reduction occurs [24].The PRI reflects the epoxidation state of the xanthophyll cycle pigments and photosynthetic efficiency [25], indicating plant physiological status.On the other hand, narrow-band red-edge indices that are sensitive to chlorophyll concentrations have also been proven to be effective in early stress detection [26].For example, [27][28][29] found that the ratio between red-edge (694 nm) and near infrared (760 nm) was superior to thermal imagery in early stress detection, which allowed identification of stress about two weeks prior to the signs of visual damage.The reason for the 694/760 nm ratio being most effective in pre-visual stress detection is probably due to the fact that these wavelengths correspond to two important physiological measures of plant water stress-chlorophyll fluorescence and vegetation water content [30].It is worth noting, however, that translating chlorophyll fluorescence data directly to plant physiology (e.g., photosynthesis) is challenging because fluorescence, photosynthesis and non-photochemical quenching are closely connected.For example, (1) they compete for the use of absorbed light, and (2) the relationship between chlorophyll fluorescence and photosynthesis is positive under stress and negative for no-stress conditions [31].In open field environments, this can change depending on the time of the day and changes in light conditions.
Few studies of fluorescence, PRI and red-edge indices have been conducted to study the impact of ozone on crops in an open air field environment with remotely sensed data.Furthermore, only a few authors have attempted to assess ozone damage at leaf [32] and canopy [33] scales using PRI.While obvious decreases in chlorophyll, nitrogen content and therefore, crop yield were observed with elevated ozone, Ainsworth et al. [34] found no effect of ozone on photosynthetic capacity, which might be attributed to insufficient temporal sampling of spectral measurements.Thus, there is a strong need to compare the performances of fluorescence, PRI and red-edge indices on characterizing crop responses to ambient ozone concentrations in a field growth environment.Although the 531 nm wavelength has been accepted as the best PRI measurement band, there has been no consensus in literature on the use of the PRI reference band (570 nm).For example, 512 nm [35], 515 nm [20] and 586 nm [22] have all been reported as the "best" reference band in recent years, emphasizing the need for additional work to determine the ideal PRI reference band.There have been numerous studies focused on macronutrient effects of abiotic stress, including salt stress [36], heavy metal pollution [37], low temperature [38] and water stress [39,40]; no studies reported mineral nutrient accumulation in plant tissues in response to ambient ozone concentrations.
While long time exposure to water or salt stress may change the flowering time and reduce seed production, short term stress effects can be mitigated by appropriate farm management and irrigation practices before permanent damage occurs.However, foliar injury on plants from ozone can be permanent and very difficult to be mitigated.Efforts should focus on reducing ambient ozone concentration and developing ozone-tolerant crop breeds.The objectives of this study are (1) to determine the responses of soybean cultivars to ambient ozone on soybean crops in open agricultural ecosystems using hyperspectral data and plant mineral nutrient concentration, and (2) to develop remote sensing methods that improve early ozone stress detection for proactive planning with respect to food insecurity.

Experimental Design and Abiotic Conditions
The experimental site was located in Maryland Heights, Missouri, USA (38 • 42 48.10N, 90 • 30 22.81 W).The site was 4 × 4 m in dimension with five 0.76 m spaced rows, and soil type was Peers Silk Loam soil.On 23 April 2014, 100 seeds (0.015 kg) were sown in each row, corresponding to five soybean genotypes, including AK-HARROW, PI88788, DWIGHT, PANA, and WILLIAMS82.AK-HARROW and PI88788 were found to be relatively sensitive [41], while DWIGHT, PANA and WILLIAMS82 were identified as relatively tolerant to ozone [4,5].The number of plants were 55, 48, 71, 65 and 56 for AK-HARROW, PI88788, DWIGHT, PANA, and WILLIAMS82, respectively.The flower color was white for AK-HARROW and WILLIAMS82, purple-pick for DWIGHT and PANA, and purple for PI88788.All abiotic conditions were the same for the study organisms during the experiment.Monthly accumulated rainfall was 177 mm in June, 28 mm in July, 124 mm in August, and 166 mm in September.It was assumed that nutrient variations among the different soybean genotypes were most likely dominated by their tolerance to elevated ambient ozone concentration.However, the contribution of fertilizers and rainfall was not taken into account.Ozone-damaged leaves generally show tiny light-tan irregular spots, bronzing of small darkly pigmented areas, and reddening [42], which started to emerge in early August.

Hyperspectral Reflectance and Leaf Area Index Measurements
All of the measurements were made under clear-sky conditions between 10:00 a.m. to 2:00 p.m. once per week during the growing season of the soybean genotypes.Leaf reflectance spectra were collected on three leaves from six plants (n = 6) within each genotype (row).The three leaves were the third trifold from the top, middle leaf and the bottom leaf down the canopy, representing three canopy layers based on documented agronomic canopy layer standards [43,44].We used high resolution full range PSR-3500 (Portable Spectroradiometer, Spectral Revolution, Inc., Lawrence, MA, USA).The wavelength range of PSR-3500 is 350-2500 nm with a resolution of 3.5 nm in the 350-1000 nm range, 10 nm in the 1000-1900 nm range and 7 nm in the 1900-2500 nm range.A reference spectrum was taken from a 99% Spectralon calibration panel (Labsphere, Inc., North Sutton, NH, USA) before target measurement and repeated for every other five measurements to re-adjust the baseline to account for any change in illumination.Then, a leaf clip with a bifurcated fiber-optic connected to both the device and a 5-watt tungsten halogen lamp light source was used to record leaf reflectance readings with a black background.The PSR-3500 was configured to automatically average 40 scans for each sampling.The raw spectra bandwidth was interpolated to 1 nm for further analysis.
Leaf area index (LAI) was measured throughout the growing season using the LAI-2200C Plant Canopy Analyzer (LI-COR Inc., Lincoln, NE, USA).Two above canopy readings taken at the beginning and end of the transects were paired by three below canopy readings along the transect crossing two consecutive crop rows.The sensor view direction with 45 • view cap was positioned looking along the row (view direction is parallel to the row) in the first transect and looking across (perpendicular to) the row in the next transect.Extra sky radiation measurements were taken for scattering corrections.

Gas Exchange and Fluorescence Measurements
Simultaneous measurements of individual gas exchange and Chl fluorescence variables, including photosynthesis (µmol CO 2 m −2 s −1 ), stomatal conductance (mol H 2 O m −2 s −1 ) and steady-state fluorescence were made on 13 August 2014 using LI-6400XT Portable Photosynthesis System installed with a pulse amplitude-modulated leaf chamber fluorometer (Li-Cor, Inc., Lincoln, NE, USA) (Table 1).The measurements were taken on the third leaf from the top of the six plants used for leaf reflectance within each genotype.The measurement was made on this day because visual signs of foliar ozone damage occurred during this growth stage (Figure 1).The chlorophyll fluorescence yield, which measures the instantaneous photochemical efficiency of photosystem II (PSII) was determined by ∆F/Fm' = (Fm' − Fs)/Fm' [45], where Fm' is the maximum fluorescence yield under environmental conditions (light adapted), and Fs is the steady-state fluorescence (light adapted) measured before the saturation pulse.Table 1.Leaf gas exchange and Chl a fluorescence measured for the different genotypes, under conditions of light intensity = 1000 µmol m −2 s −1 , CO 2 = 400 µmol CO 2 mol −1 , light pulse for Chl fluorescence ~8000 µmol m −2 s −1 , Tblock = ambient T, air flow = 300 µmol s −1 .Values are means ± SD, n = number of measurements for each genotype.

Foliar Nutrient Concentration and Yield Determination
A total of 180 leaves, used for leaf reflectance measurements on both 17 July and 13 August 2014, were destructively collected and sent to Danforth Plant Science Center located in St. Louis, Missouri, USA for the analysis of nutrient concentrations using ion coupled plasma mass spectroscopy (ICP-MS, Perkin-Elmer SCIEX).The leaf samples for the ICP-MS were carried out following the methods described in [47].K, Ca, Mg, Fe, Mn and Cu are the important elements for plant photosynthesis and gas exchange concentration [48].Therefore, we focused our analysis on the aforementioned six nutrient elements out of the 19 elements that were measured in the lab.
For each soybean genotype, the seed yield (kg/ha) and biomass were obtained using an Almaco small bundle thresher (Almaco, Nevada, IA, USA) on 1 October 2014.Threshed biomass with the soybean pods retained over a row length of 1.21 m and a width of 0.76 m was weighted to determine above-ground biomass.Harvest index (HI) was computed using the ratio between seed yield (g) and above-ground dry biomass (g) [49].

Ozone Data
Ozone measurements were obtained from a 2B Technologies ozone sensor (Model 106-L) located ~22 km southwest of the experimental site at one of the St. Louis Ozone Gardens located in Forest Park [41,50].This instrument has operated nearly continuously since 2012 and is designed for the measurement of ambient ozone concentration down to a very low concentration.The instrument has been designated by the Environmental Protection Agency (EPA) as a Federal Equivalent Method (FEM) for monitoring compliance with the US Clean Air Act (EQOA-0914-218).The data were readily available from the instrument's data archive (http://go3project.com),which has 15-min resolution.Analyses of these data have been shown to be consistent with archived data at Missouri's Department of Natural Resources (DNR), where the archived time resolution is 1 h.None of the DNR's monitoring stations is significantly closer to the experimental site than the location of the one in the Ozone Garden used for this study.

Statistical Analysis
Among the five soybean genotypes, statistics were performed on the genotype means of nutrient concentrations to investigate the effect of early (17 July 2014) and late (13 August 2014) growth stage on nutrient variability.The day of the year was considered to be a fixed effect.A multiple comparisons test using Bonferroni adjusted t test, was carried out in order to determine which pairs of genotype had differing means.

Finding the Best Spectral Regions to Map Fluorescence Yield
One of the methods to find the link between biophysiological variables and spectral information is building indices in the form of mathematical formulae integrating spectral wavebands.In this section, we describe how we identified the wavelengths and/or normalized indices (Equation ( 1)) that are capable of estimating fluorescence yield, a sensitive indicator of crop health and ozone damage.We used the data collected on both 17 July and 13 August 2014 to develop the index.Once developed, the index has the potential to track plant biophysiological changes throughout the growth period.
We used the normalized difference spectral index (NDSI) to identify optimal wavelengths and/or indices.The NDSI is defined as: where R is the reflectance value and subscripts are wavelengths in nm.NDSI were calculated for the leaf level hyperspectral reflectance spectra using all possible combinations of available bands (i and j nm).Linear relationships between in-situ fluorescence yield measurements and NDSIs were examined, and raster maps of coefficient of determination (R 2 ) and corresponding p-values were generated.Raster maps of NDSI allowed the evaluation of the different band combinations and the selection of a sensitive NDSI for crop health and ozone damage detection.
The performance our indices was then compared with previously published indices from the relevant literature as shown in Table 2.

Partial Least Squares Regression (PLSR) Modeling
The variation in nutrient concentration within five soybean genotypes was analyzed.The variation was calculated for major nutrients by subtracting the concentration values on 17 July 2014 from the concentration values on 13 August 2014.
The development of mineral nutrient prediction models was conducted with the Partial Least Square (PLS) package in R statistical software [72].In order to increase the sample size and model robustness, leaf nutrient concentrations and reflectance data from two dates were combined to develop universal models using partial least squares regression (PLSR).To select the optimum number of latent variables (LV) and avoid overfitting, we calculated the minimum root mean squared error of cross-validation (RMSE-CV) (Equation ( 2)) by the iterative leave-one-out method and relative RMSE-CV (RMSE%) (Equation ( 3)).For each nutrient predictive model, the coefficient of determination (R 2 ) between the predicted and observed parameters was calculated.R 2 indicates the impact of each independent variable on the dependent variable-the greater the R 2 value, the greater the impact on the dependent variable of the response variable.RMSE-CV was calculated by the following equations: where n is the number of samples in test set, y p and y i are the predicted and measured nutrient concentrations of the sample i. ŷ is the mean of the measured nutrient concentration.

Leaf Gas Exchange and Yield
Figure 2 shows the genotype-average fluorescence yield (∆F/Fm'), photochemical quenching (qP) and photosynthetic rate (µmol CO 2 m −2 s −1 ) measured on 13 August 2014 along with the harvest index and crop yield.The highest and lowest fluorescence yields were recorded with AK-HARROW and PI88788, respectively, and the same trend was true for the photochemical quenching values of the genotypes.The second highest values of ∆F/Fm' and qP were found with WILLIAMS82, which was followed by PANA and DWIGHT.The highest yield was obtained from DWIGHT, and the second-highest yield from WILLIAMS82, which was closely followed by AK-HARROW and PANA.The harvested yield was the lowest for PI88788.∆F/Fm' estimates the light use efficiency of plant PSII for electron transport.Therefore, ∆F/Fm' reflects the ability of a plant to tolerate environmental stresses and the extent to which those stresses have damaged the photosynthesis [31].Chlorophyll fluorescence emission, photochemical quenching (light energy is harvested for photosynthesis) and non-photochemical quenching (excess energy is dissipated as heat) occur in competition [31].When a plant is exposed to abiotic stress, light energy used for photochemistry decreases due to protective closure of the PSII apparatus, which may increase in steady-state fluorescence and/or non-photochemical quenching.In general, increased steady-state fluorescence is an indication of decreased ∆F/Fm' and therefore, is an early sign of abiotic stress.∆F/Fm' is positively correlated with qP and photosynthetic rate.This is evident in Figure 2 with a minor exception of the photosynthetic rate such that one would expect the lowest photosynthetic rate for PI88788 instead of PANA.It is worth noting that plant tissues tend to shift toward increasing heat production at the initial and intermediate states of stress [12].This may imply that much of the excess energy absorbed by PANA is being dissipated as heat rather than steady-state fluorescence causing the higher ∆F/Fm' and lower photosynthetic rate.
Harvest index, calculated as the ratio of seed yield by the total pounds of above ground biomass, is often used as a measure of reproductive efficiency.Harvest index reflects breeding achievements and is affected by stress and/or crop management practices during crop reproductive development [73].The low-high pattern of the harvest index among the soybean genotypes was surprisingly similar to that of both ∆F/Fm' and qP, given the fact that gas exchange measurements were a snapshot of plant physiological status on a specific date of measurement, while harvest index was a reflection of crop growth throughout the entire growth cycle.This was supported by our field observations that the study organisms had experienced different degrees of ozone damage by 13 August 2014 depending on their tolerance to ambient ozone concentrations; this confirms that the measurements made on this day were reflective of overall crop health and potential harvest yield.

Mineral Nutrient Uptake Rate
A significant number of publications have reported nutrient variations and changes in mineral nutrient uptake rate in plants due to water, salt stress and cold temperature stress [38,74,75].However, to our best knowledge, very few researchers have focused on ion/nutrient effects of ambient ozone concentrations.
Mineral concentration was very different among the five soybean genotypes at the early and late growth stages.One-way ANOVA test shows that K and Fe are not statistically significant among the soybean genotypes at the early date (17 July 2014 (Figure 3)) which means there were no significant differences in K and Fe concentrations across the five soybean genotypes (α < 0.05).Data from 13 August 2014 (Figure 4) shows that there were no significant differences in Ca and Fe concentrations among the five different soybean genotypes (α < 0.05).This may imply that Fe is not a sensitive indicator of ozone stress, and the other elements may be used to characterize the ozone effects at different plant phenological stages.In contrast, Cu, Mg and Mn variations were consistently significant in both the early and late growth stages and K variation was significant at the late growth stage in regard to increased ozone damage, which may provide useful insight into characterizing plant response to ambient ozone concentration.
Figure 5 shows the percent change in mineral nutrient accumulation between the 17 July and 13 August 2014 observation dates, which was calculated by dividing the mineral nutrient accumulation difference between 17 July and 13 August over the measurement taken on 13 August.The uptake rates of all mineral nutrients were significantly affected by the increasing ozone concentration.Previous studies have demonstrated that AK-HARROW and PI88788 are relatively ozone-sensitive genotypes compared to DWIGHT, PANA and WILLIAMS82 [41].The highest increase in tissue Ca was found with PI88788, which was followed by AK-HARROW, DWIGHT and WILLIAMS82.PANA showed the smallest increase in Ca, which suggests PI88788 and PANA are the most and least ozone-sensitive genotypes among the study organisms.Ca is a crucial regulator of plant growth and development, which has been implicated in an increase in plant tissues' responses to a number of abiotic stresses, including low temperature, osmotic stress, heat, oxidative stress, anoxia, and mechanical perturbation [38].It could be possible, however, that a prolonged elevated cellular concentration of Ca is detrimental to plant cells [38].
As an essential nutrient, K plays an important role in plant biochemical and physiological processes, e.g., photosynthesis and respiration, that influence plant growth and metabolism [76].K showed very distinct decreasing and increasing trends between the two groups.The percent changes in leaf tissue K were −9.7% and −3.5% for AK-HARROW and PI88788 genotypes and 29.9%, 34.9% and 20.8% for DWIGHT, PANA and WILLIAMS82, respectively.Plant growth under stressed conditions is restricted by osmotic effects and toxic effects of ions, which results in lower nutrient uptake and inhibits the translocation of K [76].The striking decrease-increase trend of K tissue concentration between the two groups of soybean genotypes supports the argument that AK-HARROW and PI88788 are less tolerant to ozone damage than the other genotypes.It should be noted that Mg and Mn changes were opposite to the K variation among the genotypes with a minor exception that both PANA and WILLLIAMS82 showed slight decreases in Mg (−6.1%) and Mn (−8.6%).K regulates the translocation of photosynthates into sink organs, maintenance of turgor pressure, activation of enzymes, and reduces excess uptake of ions when plants are exposed to abiotic stress [77].The higher concentrations of Mg and Mn in ozone-sensitive genotypes, especially in PI88788, may be associated with decreased K concentrations in these genotypes.Although it is not possible to establish a direct one-to-one relationship, the decrease-increase pattern of nutrient concentrations among the study organisms is arguably consistent with yield (Figure 2).Overall, the increase in foliar concentrations of Ca, Mg and Mn, and the decrease of K among ozone-sensitive genotypes was similar to the observations in salt stressed plants [78], indicating possible similarities in nutrient uptake and stomatal behaviors between ozone and salt stress mechanisms.However, our work showed some unique observations of plant responses to ozone stress, which were evident from the differences between leaf tissue mineral nutrient accumulation among the ozone-sensitive and tolerant genotypes.Moreover, except for WILLLIAMS82, the difference in Cu concentration between the two measurement dates was similar to that of K, i.e., decreased in ozone-sensitive genotypes and increased in ozone-tolerant genotypes, respectively.Similarly, mixed Cu signals due to salt stress has been reported by other studies involving soybean [36] and canola [79] genotypes.However, the majority of publications have confirmed that salt stress mostly reduces the leaf tissue Cu concentration in less stress-tolerant plants [78], and the foliar Cu concentration can increase until a critical stress level is reached, before decreasing beyond the critical level [80].The decrease in foliar Cu concentration in AK-HARROW, PI88788 and WILLLIAMS82 on 13 August 2014 may be attributed to their stress tolerance limits, which might have already been reached.In contrast, the severity of ozone stress was not yet inhibitory for DWIGHT and PANA in regard to the uptake of Cu and other nutrients.

Best Spectral Regions Correlated with Fluorescence Yield
Figure 6a,b present the heat map of R 2 and p-values between the fluorescence yield (∆F/Fm') and the normalized difference spectral index (NDSI [R i ,R j ]) using complete spectrum reflectance combinations of two wavebands at R i and R j .NDSI is defined as (R i − R j )/(R i + R j ).As shown in Figure 6c, there were six significant spectral regions (A, B, C, D, E and F), where R 2 ranged between 0.36 and 0.64.The hotspot regions, C, D, E and F, within NIR, and the shortwave infrared (SWIR) bands tended to have higher R 2 (0.42 < R 2 < 0.64) compared to the regions, A and B, which only included visible and red-edge bands.∆F/Fm' had the highest R 2 of 0.64 (p < 0.001, RMSE = 0.05) with NDSI [R 416 , R 2371 ] that was found in region D. Blue and red portions of photosynthetically active radiation (400-700 nm) more efficiently excite chlorophyll than green portions [81].Compared to red light, blue light induces stronger chlorophyll fluorescence emission, which may be attributed to the fact that (1) blue light is strongly absorbed by photosynthetic pigments (chlorophyll ab and carotenoids), (2) has low penetration ability, and (3) blue light induced fluorescence is mostly from the upper leaf half, rather than from subjected reabsorption [82].This could be the reason for the selection of blue band 419 nm as a measurement band, while the SWIR band was used as a reference band as it is not affected by photosynthetic pigments.Nevertheless, due to the lower signal to noise ratio (SNR) of detectors caused by physical constraints, wavelengths in the SWIR spectral region may have low SNRs.Satellite remote sensing requires comprehensive atmospheric correction to reduce the influence of atmospheric water vapor [14].Hotspot areas were located in near infrared (NIR) bands.∆F/Fm' was most correlated with NDSI [R 523 , R 533 ] (R 2 = 0.35, p < 0.0005, RMSE = 0.08) in region A and with NDSI [R 613 , R 614 ] (R 2 = 0.36, p < 0.0005, RMSE = 0.07) in region B. The extent of these significant regions (510-590 nm and 450-550 nm for region A, 560-710 nm and 550-670 nm for region B) is spectrally wide, covering visible and red-edge regions.Visible and red-edge spectral regions are known to be sensitive to small changes in light absorbing pigments and leaf cell structure [29,83].Further, in comparison with the SWIR region, visible and infrared (VNIR) spectral regions are less effected by atmosphere and, importantly, generally available for most sensors from different platforms in a high spatial resolution [84].In our case, the wide ranges of spectral regions A and B were applicable from broad band satellite observations.In particular, regions A, B and C overlap with several VNIR bands of the Sentinel-2 sensor, which has high spatial resolution and short revisit time, indicating high potential for large scale monitoring of the effects of ambient ozone on soybean.Nevertheless, confirmation of the potential in satellite observations warrants a further, thorough investigation.The NDSI was calculated using complete spectrum reflectance combinations of two wavebands at R i and R j (n = 30).

Comparison with Existing Indices and Methods
The location of a sample leaf on a plant may have some influence on observations, because leaf biochemical concentrations and light exposition showed decreasing trends, while leaf age and LAI increased in order of the upper, middle and lower canopy layers [43,85].Our results show that for upper and middle leaf measurements, PRI, using 519 nm as the reference band wavelength, is the best for capturing plant physiological responses to ozone concentrations (Figure 7).Even though the observed relationship was not strong (0.3 < R 2 < 0.35), the selection of 519 nm aligns well with a previous study [35], which showed 512 nm as the best reference band to make PRI structurally robust.Overall, this finding is important as most optical remote sensing sensors actually measure the upper canopy.
Figure 8 compares the correlation relationship between the spectral indices and physiological parameters, including fluorescence yield (∆F/Fm'), electron transport rate (ETR, µmol m −2 s −1 ) and photochemical quenching (qP).Among the 27 indices, two indices developed in this study, i.e., PRI 519 [R 531 , R 519 ] and NDSI [R 419 , R 2371 ], showed the strongest correlations compared to those from Table 2.It is worth noting that ozone stress is best reflected by fluorescence yield, electron transport rate and photochemical quenching, while the strongest correlation with plant water stress was found with stomatal conductance [86].

of Mineral Nutrient Concentration
The intent of PLSR modeling was to evaluate the feasibility of using leaf reflectance spectra as a nondestructive means to estimate leaf mineral nutrients and apply developed models to predict the nutrient concentration of soybean cultivars throughout the growing season.The best model was achieved between averaged spectral data with respect to both growth stages and genotypes.Datasets representing different growth stages and various cultivars allowed us to build general models increasing the chance of success in developing robust models for estimating nutrient concentration, thanks to the representativeness of larger dataset [87].Table 3 presents the optimal models for mineral nutrient prediction.Coefficient of determination values (R 2 ) during cross-validation ranged between 0.54 and 0.87, indicating that 54-87% of variance in the mineral nutrients was explained by respective prediction models with RMSE% ranging between 9.68-17.90%.The models for macronutrients (Ca, K and Mg) tended to have higher R 2 values (>0.70) except for Mn, compared to the R 2 (<0.70) of the prediction models for micronutrients (Cu and Fe).The better performance of the macronutrient models could be attributed to the higher concentration of nutrients, thus stronger spectral response [88].The performances of micronutrient prediction models were comparable or superior to several existing studies.For example, [89] reported R 2 values of 0.65, 0.73, and 0.56 for Ca, K and Mg, respectively, using stepwise regression, averaged reflectance spectra of olive, heather, willow and mopane leaves.Similarly, [90] calibrated PLSR based macronutrient prediction models with continuum-removed reflectance spectra of Hyperion data over forested area had R 2 values for Ca, K and Mg models of 0.08, 0.31 and 0.86, respectively.The regression coefficients of each wavelength in the prediction models are shown in Figure 9.The contributions of the visible and SWIR spectral regions were significant for predicting Ca, Cu, Fe, Mn and Mg, while near infrared spectral region was more important for K compared to the other spectral regions.In general, the visible spectral region appears to be highly significant for most of the nutrient prediction models, and this finding is well-aligned with previous studies [90,91].In theory, there is no direct relationship between leaf reflectance spectra and mineral nutrient concentration, but reflectance spectra still can be used to predict nutrient concentration because of the spectrally active C-H, N-H, and O-H bonds, which are highly dependent upon the chemical composition of the leaf tissue [92].There exist inherent correlations between nutrients and leaf chemical components, such as organic matter, lignin, protein and pigments [93].Ca does not have obvious spectral signature in near infrared spectral region, and Mg is an important constituent of the chlorophyll molecule.The significant contributions of visible and SWIR spectral regions in Ca and Mg prediction could be explained by the spectral responses of organic compounds and chlorophyll content [94].Similarly, K is not spectrally active in the NIR spectral region, and the significant contribution of the NIR region could be attributed to the other chemical compounds that are sensitive to NIR [95].The inferior performance of macronutrient prediction models for Fe, Cu and Mn is mainly caused by low concentration levels [96,97].However, [98] demonstrated the potential of the NIR and SWIR spectral regions in the detection of Cu deficiency in barley leaves at early stage and [97] showed the feasibility of predicting micronutrients (Fe and Mn) in orange leaves using visible spectral region.These results were further confirmed in the current study.

Temporal Variation in Mineral Nutrient Accumulation and Physiological Indices under Ambient Ozone Concentrations
Elevated ambient ozone can negatively affect the soybean photosynthetic rate and thus, the nutrient concentration and yield.This negative effect differs between soybean genotypes [34].Overall, the nutrient concentrations of soybean leaves were significantly affected by elevated ambient ozone, and nutrient concentrations showed clear separation among different soybean genotypes from 9 July 2014 (DOY 190) towards the end of growth, which was not the case for K (Figure 10).Ca, Mg and Mn concentrations decreased until 17 July 2014 (DOY 198); Ca and Mn steadily increased towards the end.The trend was similar for all soybean genotypes, while PI88788 had consistently low Ca, Mg and Mn concentrations starting from 2 July 2017 (DOY 183).Overall, a decrease in nutrient concentration is related to a dilution effect caused by a fast increase in biomass [99].The decline in Ca and Mg at the beginning of the reproductive state could be attributed to the corresponding increase in K, which intensifies the dilution effect and inhibits Ca and Mg uptake [100].However, Ca and Mn reached their minimum values on 17 July 2014 (DOY 198) and continued to increase afterwards, regardless of the changes in K. Ca is considered to be a second messenger in plant metabolic regulation, and it increases in response to various biotic and abiotic stressors [101].Furthermore, there is a strong positive correlation between Ca and Mn [102,103]; Hirschi et al. [104] reported that tobacco leaves with increased Mn became tolerant to it.Additionally, Jiménez et al. [105] found accumulated Mn towards the end of the reproductive cycle.The fluctuation in Mg could be explained by varying influences of pods and seeds are known as sinks of soybean during reproductive stages.This unstable Mg concentration is consistent with the findings of [106].However, PI88788 showed a dramatic decrease in Mg and reached the minimum value on 2 July 2014 (DOY 183), which was followed by a small increase but was still consistently lower than other genotypes.The low Mg concentration of PI88788 can be associated with decreased chlorophyll, as most Mg is found in chlorophyll molecules [107].This result was also supported by the decreased chlorophyll index calculated from leaf spectra for the same growth stages [41].Fe and Cu fluctuated around their mean values, possibly because they are immobile and not easy to move between different parts of plants [107].Fe is necessary for biosynthesizing chlorophyll, and Fe deficiency is the main cause of chlorosis [108].Fe showed a consistently higher concentration in PI88788, and this is almost the opposite of the Mg variation trend.The increased Fe could be attributed to the need for more chlorophyll to carry out normal photosynthetic activities.Here, it is worth mentioning that ambient ozone damage on sensitive cultivar became visible after 13 August 2014 (DOY 225).About half of all available Cu in plant tissue is located in the chloroplasts, where Cu plays an important role in electron transport reactions [109].Thus, it is not surprising that the Cu concentration in PI88788 had a similar trend to Mg but was much smaller in range.Due to the activation of protective photosystem apparatus to prevent damage in high light conditions, the variations in the magnitude of fluorescence emission and photosynthesis were most significant in the afternoon when plants are exposed to high light conditions [110].On 13 August 2014, we started fluorescence and photosynthesis measurements on AK-HARROW at 10 a.m.Each sampling took about 10-15 min, with six samples per genotype, and ended at 4 p.m.One may expect higher temperature, non-photochemical quenching and therefore, lower fluorescence for PANA and WILLIAMS82 as they were measured during the peak of high temperature stress.This may be one of the factors that complicated the inter-comparison of fluorescence between the cultivars.In fact, under low light unstressed conditions (no non-photochemical protection mechanisms activated) a negative correlation exists between fluorescence and photochemistry.Most studies have observed that in the presence of plant stress and high light conditions, fluorescence declines with photosynthesis (i.e., positive correlation) as a result of protective mechanisms (e.g., deactivation of the antenna, activation of the xanthophyll cycle and non-photochemical protection) which take place in the leaf to prevent damage caused by harmful radicals formed in such stress conditions [111].

Conclusions
This paper evaluated the responses of five soybean genotypes to ambient ozone in an open agricultural experiment using hyperspectral data and plant mineral nutrient concentrations.A suite of statistical, spectral, and modeling approaches, including NDSI and PLSR were utilized.The major conclusions include the following: (1) fluorescence, photochemical quenching, electron transport, and harvest index varied significantly among the genotypes.The lowest photochemical quenching, fluorescence yield, harvest index, and terminal yield were found with the most ozone sensitive genotype PI88788, which was followed by AK-HARROW, the second most ozone sensitive genotype tested.The opposite was true for the most ozone tolerant genotypes: DWIGHT, PANA, and WILLIAMS82.(2) Ambient ozone affects how plants regulate mineral nutrient uptake and tissue element composition.When plants were exposed to significant ozone damage, the foliar concentrations of K and Mn were lower and Fe and Mg concentrations were higher for ozone-sensitive genotypes-AK-HARROW and PI88788-than the rest of the genotypes.(3) The analysis between the spectral data and leaf biophysical variables collected from the top, middle and bottom of a plant canopy demonstrated that the upper and middle leaf measurements are the best for capturing plant physiology in response to ozone concentrations.This finding is important as most optical remote sensing sensors measure upper canopy.(4) Unlike water stress, which is strongly associated with changes in stomatal conductance, fluorescence yield was found to be the most correlated to ozone stress.(5) Among the 27 indices found in the literature, the two indices developed in this study, i.e., PRI 519 [R 531 , R 519 ] and NDSI [R 419 , R 2371 ], showed the strongest correlations with fluorescence yield, electron transport, and photochemical quenching, followed by PRI 512 [R 531 , R 512 ], Anthocyanin (Ant gmn ), Carotenoid Index (Car gtln ), and Fluorescence Ratio Index1 (FRI 1 = R690/R600).
This contribution represents one of the very few studies involving the characterization of ambient ozone concentration on soybean crops in an open agricultural environment.Future research should focus on scaling up the results to satellite observations by conducting field-airborne-satellite simultaneous campaigns integrated with radiative transfer simulations.The SWIR wavelength 2371 nm was selected due to having the highest R square values among all bands.However, there are a number of high R-square hotspots in the VNIR region with R-square values very close to the highest R-square of the SWIR band.These VNIR regions, available from satellite sensors (Sentinel-2, Landsat, Worldview-2/3), provide wavelengths for satellite spectral indices.

Figure 2 .
Figure 2. Genotype-average fluorescence yields (∆F/Fm'), photochemical quenching (qP) and photosynthetic rates (µmol CO 2 m −2 s −1 ) measured on 13 August 2014, along with harvest indexes and crop yields.The lowest crop yield, harvested from PI88788, corresponds to the lowest values of ∆F/Fm' and qP.Ozone tolerant genotypes produced overall higher yields than the sensitive genotypes.

Figure 3 .
Figure 3. Variation in mineral nutrient concentration among soybean species on 17 July 2014.Whiskers indicate the 5th and 95th percentiles.Solid lines inside the box represent the median; the colored areas represent the 25th and 75th percentiles.The significant level is also shown (*** p < 0.05).The results indicate that K and Fe are less significant among soybean cultivars at the early growth stage.When plants were not yet exposed to significant ozone damage, on 17 July 2014, foliar concentrations of Ca, Fe, Mn were lower and K concentration was higher for ozone-sensitive cultivars AK-HARROW and PI88788 than the rest of the cultivars.

Figure 4 .
Figure 4. Variation in mineral nutrient concentration within soybean genotypes on 13 August 2014.Whiskers indicate the 5th and 95th percentiles.Solid lines inside the box represent the median; the colored areas represent the 25th and 75th percentiles.The significant level is also shown (*** p < 0.05).At this growth stage, Ca and Fe variations are less significantly different among soybean genotypes.When plants were exposed to significant ozone damage on 13 August 2014, the foliar concentrations of K and Mn were lower and the Fe, Mg concentrations were higher for ozone-sensitive genotypes AK-HARROW and PI88788 than for the rest of the genotypes.

Figure 5 .
Figure 5. Variation in mineral nutrient concentration within five soybean genotypes.Differences were calculated for major elements by subtracting the concentration values on 13 August 2014 from measurements on 17 July 2014.

Figure 6 .
Figure 6.Coefficient of determination (R 2 ) (a) and corresponding p-values (b) between the fluorescence yield (∆F/Fm') and NDSI (R i ,R j ).(c) shows regions of interests that have large R 2 and high significance.The NDSI was calculated using complete spectrum reflectance combinations of two wavebands at R i and R j (n = 30).

Figure 7 .
Figure7.The relationship between fluorescence yield and PRI, calculated with varying reference bands in the range, 400-1000 nm.For all the leaf positions (upper, middle, lower and average leaf spectra), the relationship is strong when the reflectance of 519 nm wavelengths are used to formulate PRI.

Figure 8 .
Figure 8.Comparison of correlations between the spectral indices and the physiological parameters.* denotes the indices proposed in this study.

Figure 10 .
Figure 10.Changes in nutrient concentrations (mean and ± standard error of the mean) by soybean genotypes over time.

Table 2 .
Spectral indices used in this study.

Table 3 .
Results of the optimal model's performance in predicting mineral nutrients.