Ethylenediurea (EDU) effects on hybrid larch saplings exposed to ambient or elevated ozone over three growing seasons

Ground-level ozone (O3) pollution is a persistent environmental issue that can lead to adverse effects on trees and wood production, thus indicating a need for forestry interventions to mediate O3 effects. We treated hybrid larch (Larix gmelinii var. japonica × L. kaempferi) saplings grown in nutrient-poor soils with 0 or 400 mg L−1 water solutions of the antiozonant ethylenediurea (EDU0, EDU400) and exposed them to ambient O3 (AOZ; 08:00 − 18:00 ≈ 30 nmol mol−1) or elevated O3 (EOZ; 08:00 − 18:00 ≈ 60 nmol mol−1) over three growing seasons. We found that EDU400 protected saplings against most effects of EOZ, which included extensive visible foliar injury, premature senescence, decreased photosynthetic pigment contents and altered balance between pigments, suppressed gas exchange and biomass production, and impaired leaf litter decay. While EOZ had limited effects on plant growth (suppressed stem diameter), it decreased the total number of buds per plant, an effect that was not observed in the first growing season. These results indicate that responses to EOZ might have implications to plant competitiveness, in the long term, as a result of decreased potential for vegetative growth. However, when buds were standardized per unit of branches biomass, EOZ significantly increased the number of buds per unit of biomass, suggesting a potentially increased investment to bud development, in an effort to enhance growth potential and competitiveness in the next growing season. EDU400 minimized most of these effects of EOZ, significantly enhancing plant health under O3-induced stress. The effect of EDU was attributed mainly to a biochemical mode of action. Therefore, hybrid larch, which is superior to its parents, can be significantly improved by EDU under long-term elevated O3 exposure, providing a perspective for enhancing afforestation practices.


Introduction
Larches (Pinaceae), unlike the vast majority of conifers, are deciduous trees (Gower and Richards 1990;Givnish 2002;Osawa et al. 2010). They cover permafrost ecosystems in eastern Siberia (Osawa et al. 2010), and dominate nutrientpoor areas of the boreal zone and various high-latitude Abstract Ground-level ozone (O 3 ) pollution is a persistent environmental issue that can lead to adverse effects on trees and wood production, thus indicating a need for forestry interventions to mediate O 3 effects. We treated hybrid larch (Larix gmelinii var. japonica × L. kaempferi) saplings grown in nutrient-poor soils with 0 or 400 mg L −1 water solutions of the antiozonant ethylenediurea (EDU0, EDU400) and exposed them to ambient O 3 (AOZ; 08:00 − 18:00 ≈ 30 nmol mol −1 ) or elevated O 3 (EOZ; 08:00 − 18:00 ≈ 60 nmol mol −1 ) over three growing seasons. We found that EDU400 protected saplings against 1 3 peatlands of the northern hemisphere (Givnish 2002). Their dominance in the latter habitat is often described as a paradox because evergreen vegetation would be anticipated to dominate (Givnish 2002). The advantage of larches over most other deciduous and evergreen trees can be partly explained by the higher nutrient retranslocation efficiency of larch (Matyssek 1986;Givnish 2002;Kitaoka et al. 2020) and advantageous water relations (Berg and Chapin 1994).
Owing to their unique ecophysiological characteristics, larches have a very wide distribution throughout the northern hemisphere (Givnish 2002;Usoltsev et al. 2002;Ryu et al. 2009). Larches are major trees for wood and biomass production in large regions of the northern hemisphere, thus being of considerable economic and environmental (CO 2 sink) value (Zhang et al. 2000;Usoltsev et al. 2002;Osawa et al. 2010;Kurinobu 2015;Ohta et al. 2019). Hence, important breeding programs targeted the improvement of larches by creating hybrids in Europe and Asia (Kurinobu 2015).
Japanese larch (Larix kaempferi (Lamb.) Carr.) has been a major afforestation species in northeast Asia, including Japan, since a long time ago (Ryu et al. 2009;Kurinobu 2015). However, it displays a number of disadvantages that can result in reduced growth and timber production, such as susceptibility to root rot and shoot blight diseases, deer and redback voles grazing, wind, and snow (Ryu et al. 2009;Kurinobu 2015). Therefore, long-lasting breeding programs have led to a new hybrid (F 1 ) with improved characteristics (Kurinobu 2015). This hybrid between Dahurian larch (L. gmelinii var. japonica) and Japanese larch (L. kaempferi) was found to be promising and planted widely in northern Japan (Ryu et al. 2009;Kurinobu 2015). Seedlings of hybrid larch are planted into natural soils in the field for afforestation, without any fertilization. However, while hybrid larch displays heterosis under various environmental conditions, it appears to be susceptible to ground-level ozone (O 3 ) (Sugai et al. 2018;Wang et al. 2018).
Ozone concentrations are elevated throughout the globe, relative to the pre-industrial concentrations, and O 3 is one of the most puzzling-to-control air pollutants nowadays Gao et al. 2020;Diaz et al. 2020). Current ambient O 3 concentrations are potentially toxic for agricultural crops, trees/shrubs, and other types of vegetation in many areas of the northern hemisphere (Singh et al. 2018;Käffer et al. 2019;Sicard et al. 2020;Proietti et al. 2020;Mukherjee et al. 2021), especially in the Asian region (Koike et al. 2013;Izuta 2017). Ozone is taken up by plants during gas exchange with the atmosphere via stomata, thus generating the production of harmful chemical species and activating stress defense mechanisms (Cieslik et al. 2009;Kitao et al. 2015;Moura et al. 2018;Bellini and De Tullio 2019;Kinose et al. 2020;Xu et al. 2021). These responses lead to inhibitory effects above taxon-specific stress thresholds, including inhibition of photosynthesis and photosynthates translocation, growth, yields, and productivity (Fiscus et al. 2005;Matyssek et al. 2010). Ozone pollution can affect a multitude of ecological processes and functions (Ueno et al. 2016(Ueno et al. , 2020Fuhrer et al. 2016;Yue et al. 2017;Unger et al. 2020;Blande 2021), as well as reduce the yields of agricultural crops (Mills et al. 2018;Feng et al. 2019a andHu et al. 2020). Moreover, O 3 pollution can significantly reduce annual wood production (Sacchelli et al. 2021), creating a need to protect trees that are important for wood production.
Ethylenediurea (EDU) is a synthetic chemical that was reported to protect bean plants from O 3 injury in the 1970s, and has been extensively studied for its phytoprotective mode against O 3 injury since then (Manning et al. 2011;Singh et al. 2015;Agathokleous 2017;Tiwari 2017). Although EDU improves various physiological, biochemical, growth, and productivity traits of plants under O 3 stress in the framework of hormesis, there is high intraspecific variation and no universal mechanism explaining its phytoprotective mode of action (Feng et al. 2010;Oksanen et al. 2013;Singh et al. 2015;Agathokleous 2017;Tiwari 2017;Ashrafuzzaman et al. 2018;Pandey et al. 2019;Surabhi et al. 2020). Besides, EDU studies concerning tree species are few Feng et al. 2010;Agathokleous 2017;Giovannelli et al. 2019) relative to EDU studies concerning non-tree species. Importantly, there have been only two research programs studying the efficacy of EDU in protecting tree species over several years; one with poplars and one with Japanese larch (Giovannelli et al. 2019;Agathokleous et al. 2021). Therefore, the effects of EDU on plants under O 3 stress over several years remain largely unexplored, yet it is essential to confirm whether EDU can protect various trees against O 3 for multiple years, indicating that long-term responses should be evaluated.
In the framework of a previous study, two independent experiments were performed to study EDU-O 3 effects on Japanese larch (hybrid larch's pollen father), showing that EDU significantly protected Japanese larch plants against O 3 -induced phytotoxicty ). In the first experiment, the response of seedlings (size and biomass) to increasing doses of EDU was studied in a free-air O 3 -concentration enrichment (FACE) system, by exposing seedlings to ambient or elevated O 3 for one growing season. In the second experiment, saplings were treated with EDU while exposed to only ambient O 3 in two growing seasons and only elevated O 3 in the following two growing seasons, and their size and biomass were studied ). In the present study, we evaluated the effects of EDU (0 or 400 mg L −1 ) on hybrid larch saplings exposed to either ambient or elevated O 3 for three growing seasons using a contemporary FACE facility. An advantage of this experiment, compared to the vast majority of studies on O 3 effects on plants, is that a non-fertilized soil was used to create near real-world conditions. Several traits of gas exchange, photosynthetic pigments, growth, and biomass were studied, as well as litter decomposition. We hypothesized that elevated O 3 would negatively affect gas exchange and photosynthetic pigments, thus also inhibiting growth and biomass production, effects that EDU can alleviate (Feng et al. 2010;Oksanen et al. 2013;Singh et al. 2015;Agathokleous 2017;Tiwari 2017). Because of the decreased energy reservoirs under severe O 3 stress, we hypothesized that plants may show altered allocation of resources toward buds, which define next year (future) re-growth, and became interested in whether EDU can protect against such effects. As elevated O 3 decreases foliage quality, especially causing visible injuries on the leaves (e.g., necrosis, chlorosis, or discoloration), we also hypothesized that such low-quality leaf litter may lead to changes in litter decay, and that EDU may protect against such O 3 -induced effects by protecting leaves against O 3 toxicity. Elevated O 3 can also impair leaf N resorption , which may affect litter decay rate and plant growth in long-term. For example, if the decay rate of potentially N-enriched litter is faster in elevated O 3 than in ambient O 3 , potentially inhibited plant growth due to impaired leaf N resorption might be alleviated to some extent. By showing that EDU can significantly improve hybrid larch seedlings under O 3 stress, this study offers a perspective to enhance the success of afforestation with hybrid larch by preconditioning plants to O 3 via the exogenous application of EDU. This experiment differs from the previous experiments with larches exposed to EDU-O 3  in that this experiment includes (1) hybrid larch, (2) both ambient and elevated O 3 treatments in three growing seasons, and (3) assessment of physiological mechanisms as well as growth over time.

Experimental site
The experiment was conducted at the Sapporo Experimental Forest (43°0' N, 141°2' E, 15 m a.s.l.) of the Field Science Center of Northern Biosphere (FSC) of Hokkaido University, Sapporo, Japan, in three growing seasons (2014 − 2016). In this area the snow-free season commonly spans from late April to mid-November. Data of meteorological parameters were logged by a nearby monitoring station (43°03.6′ N, 141°19.7′ E) operated by the Japan Meteorological Agency (http:// www. jma. go. jp/ jma/ indexe. html). The average values (± indicates standard deviation for environmental data and standard error for experimental/statistical results hereafter) of the main parameters were calculated for the months April-October of 2014, 2015, and 2016. The monthly average, daily maximum, and daily minimum air temperature was 16.2 ± 5.2, 20.8 ± 5.1, and 12.5 ± 5.6 °C, respectively, while the wind speed was 3.69 ± 0.53 m s −1 . The relative humidity, monthly total sunshine duration, and monthly precipitation were 67% ± 7%, 1877 ± 37 h, and 106 ± 63 mm, respectively.

Plant material
Hybrid larch two-year-old saplings were stored in dark under low temperature, at the nursery of Forestry Research Institute, Hokkaido Research Organization, Bibai (near Sapporo), Japan. Forty-eight uniform saplings were selected and transported to the campus of Hokkaido University in Sapporo, on May 15, 2014, and kept in an incubator at a constant temperature of 3 ± 0.5 °C under light. The saplings were planted in 15-L pots on June 5, 2014. The pots contained an 1: 1 mixture of two commercial soils (pH = 5.92 ± 0.02), viz. well-weathered volcanic ash Akadama soil and Kanuma pumice soil (DCM Homac CO., LTD., Sapporo, JP), which contain no organic matter and are poor in nutrients (Agathokleous et al. 2016a). Soils are originated from Kanuma town of Tochigi prefecture, central Japan. Both soils are originated from volcanic ash (a kind of Vitric Andosols), and are deficient in phosphorus and poor in nitrogen, and commonly occur in Hokkaido (Schmincke 2004 The pots were placed under ambient conditions at a sunlight-exposed ambient site of the experimental nursery of the FSC. A fully randomized design was employed, with a pot-to-pot distance of 30 cm. At planting, the average values of height and fresh weight were 13.8 ± 0.4 cm and 1.7 ± 0.1 g. The average diameter at the base, middle, and uppermost part of the stem was 2.0 ± 0.06, 1.5 ± 0.04, and 1.2 ± 0.03 mm, respectively. The saplings had on average 16.8 ± 0.5 buds and 20.2 ± 1.5 needles. The diameter was assessed with a digital caliper with an accuracy of 0.01 mm, while the height was evaluated using a measuring tape with a 1-mm gradient. The fresh weight was measured after roots were immersed in deionized water so to have uniform amount of humidity across saplings.
On 6 June, 2014, the pots were transferred into the FACE plots. Eight randomly selected saplings were placed at the south side of each FACE plot (the schematic diagram of the plots was presented elsewhere , and all plots remained under ambient conditions until the O 3 treatment began. Henceforth, and in each growing season, the pots position in each FACE plot was re-adjusted biweekly to minimize edge effects. Furthermore, the saplings of each O 3 treatment were rotated across the three FACE plots with same O 3 treatment bimonthly. Apart from the EDU and O 3 treatments, plants were not treated with any other means. Therefore, plants were also rain-fed. However, on 20 May, 2016, 950 mL of tap water were given to each plant due to unusual increases of temperature during those days (> 25 °C). Furthermore, 400 mL of water were given to each plant the night before gas exchange measurements.

EDU treatments
Saplings were treated with soil drench containing either 0 or 400 mg EDU L −1 (200 mL per pot per application). The concentration of 400 mg L −1 was selected as the upper limit of the range of EDU concentrations offering sufficient phytoprotection from O 3 toxicity based on the broad literature (Feng et al. 2010;Agathokleous 2017). This concentration also effectively protected Japanese larch (hybrid larch parent) seedlings and saplings against O 3 phytotoxicity (Agathokleous et al. 2021). Before applying EDU in 2014, half of the saplings in each FACE plot were assigned the one EDU treatment and the other half were assigned the other EDU treatment; the assignment was random. Henceforth, each sapling was receiving the same EDU treatment every 9 days during the growing seasons of 2014, 2015, and 2016. EDU was first applied on July 29, April 14, and April 24 of 2014, 2015, and 2016, respectively, and repeated 12, 22, and 17 times in each of the growing seasons. Hence, each sapling was treated with EDU 51 times in total, equalling to 10.2 L of water solution. The 9-day application interval was selected considering the known effective persistence of EDU in the leaf apoplast, and based on the widely successful results that EDU produced in the literature Agathokleous 2017). There were 4 saplings per EDU treatment per FACE plot, i.e., 12 saplings per EDU treatment per O 3 treatment.

Ozone treatments
For O 3 treatments, the FACE system with tree communities of the FSC was used, which consists of 6 plots scattered around the experimental forest. Details of the technical characteristics of this system can be found in supplemental materials published along with a recent article  was inactive, plants in the plots were exposed to ambient O 3 concentrations. Ambient O 3 concentrations were logged every 1 min (TUV-1100; Tokyo Industries Inc. Tokyo, Japan). The mean daily 10-h (08:00 − 18:00 Japan Standard Time, JST) ambient O 3 mixing ratio was 22.3 ± 3.3, 34.3 ± 5.5, and 32.3 ± 5.7 nmol mol −1 in the growing seasons of 2014, 2015, and 2016, respectively. The mean 10-h O 3 concentrations in the EOZ plots were 60.1, 72.1, and 55.8 nmol mol −1 in the growing seasons of 2014, 2015, and 2016, respectively. Further details about O 3 exposure were presented in earlier publications .

Growth
Growth of plants was monitored through time. Plant height, crown span (from two farthest points of the crown), stem basal diameter (average of two crosswise measurements per plant), and total number of buds per plant were measured. Plant height and crown span were measured (cm) with a measuring tape with an accuracy of 1 mm. Stem diameter was measured (mm; 2-decimals accuracy) with a digital caliber. Plant height, crown span, stem diameter, and number of buds of all plants were measured before the beginning of treatments (7/

Photosynthetic pigments
Photosynthetic pigments were measured in sun-lit needles collected on October 31, 2014, and July 25, 2015, and September 2, 2015. This taxon has dimorphic shoots/branches, with single needles and several buds on long shoots and dense clusters of needles (20 − 50) and a single bud on short shoots, and the normal needle senescence of the experimental plants began in November.
Some needles, usually single needles on long shoots, are not shed and overwinter, especially during the seedling stage (Kurahashi 1988), thus helping plant regrowth in the following spring-summer. For the first assessment (2014), single needles on the main shoot (stem) were analyzed to examine potential influence to be carried over in the next growing season. For the second and third assessments (2015), mature needles from dense clusters of needles on short shoots were analyzed to examine the effect in the current growing season. The sampling of July (2015) was conducted during the peak of the growing season, when severe O 3 visible injury existed in the plants treated with EDU0 and EOZ (it appeared from early June); no such visible injury existed in plants treated with EDU400 and EOZ (Fig. 1). During the September (2015) sampling, however, the severely injured needles of plants treated with EDU0 and EOZ had already been shed prematurely. Therefore, the sampling conducted in September (2015) concerns less-affected needles, in terms of O 3 -induced visible injury. Plant senescence was recorded after careful visual observation of each plant and its surrounding soil surface. All the plants were observed, and binary data were collected (yes or no) to calculate how many plants per treatment were (prematurely) senescing. There was no defoliation induced by any other factor over the course of the experiment.
For each plant, needles were put in 4 mL of DMSO and placed in an incubator at 65 °C for 4 h, while shaking from time to time (Shinano et al. 1996). Before inserting them into DMSO, the needles were scanned (Canon LIDE 40, Tokyo), and the area of the tissue (single-side projection) was measured using the LAI32 v.0.377 software (© Kazukiyo Yamamoto). Extracts were evaluated (GeneSpec III; Hitachi Genetic Systems; MiraiBio, Alameda, CA) immediately after incubation, and the absorbance was recorded at the optical densities of 470, 543, 648, and 665 nm. The values at 543 and 470 nm were used to assess carotenoids, according to the method of Lichtenthaler 1987. The values at 648 and 665 nm were used to assess chlorophyll a, chlorophyll b, and total chlorophylls (chlorophyll a + b), according to a previously described method (Barnes et al. 1992;Shinano et al. 1996). The ratios of chlorophyll a to chlorophyll b and carotenoids to total chlorophylls were also calculated.

Gas exchange
Gas exchange was measured in sun-lit needles on clear days between June 23 and 30, 2015 (all plants were measured). Mature needles of long-shoot were selected randomly; O 3 -symptomatic needles of random O 3 -injury degree were collected from EOZ. A soil drench of 400 mL water was given to each seedling the evening before the measurements.
The needles were adapted to the chamber for 20 min before starting the A/Ci curve program. After photosynthesis measurements, the measured needles were sampled for area estimation. Gas exchange was recorded with different LiCor instruments (LI-6400, Li-Cor Inc., Lincoln, NE, USA) measuring simultaneously. Plants of different experimental conditions were rotated between LiCor instruments and over time. Net photosynthetic rate (A 380 ), stomatal conductance (gs 380 ), and transpiration rate (E 380 ) were measured at ambient CO 2 (380 μmol mol −1 ), 55% − 65% relative air humidity, 25 °C leaf temperature, and a light saturated photosynthetic active photon flux (PDF) of 1500 μmol m −2 s −1 (Koike et al. 2012). The A/Ci curve (the curve between net CO 2 assimilation rate and intercellular CO 2 concentration) was developed using fluctuating CO 2 levels (12 levels) in the leaf chamber fluorometer, with constant PDF, relative humidity, and temperature. The LI-6400 AutoProgram was used. Based on the A/ Ci curve and using a C3 plant biochemical model (Farquhar et al. 1980;Long and Bernacchi 2003), net photosynthetic rate at saturated CO 2 of 1500 μmol CO 2 mol −1 (A max ), maximum rate of carboxylation (Vc max ), and maximum rate of electron transport rate (J max ) were estimated. Gas exchange measurements were conducted during morning hours (completed by noon).

Litter decay
On 10 October, 2015, mesh was installed on each pot to collect the leaf litter; senescence had not started yet (or restarted for EOZ plants that had premature senescence earlier; see Results). Leaf litter was collected from each plant on October 19, 2015, put in air-dry oven for 48 h (65 °C), and stored in zip-lock bags at room temperature and dark conditions until further use in the next snow-free season. In September 2016, collected needles of all plants of each treatment were pooled and mixed to give a pool of fallen needles for each treatment (AOZ × EDU0, AOZ × EDU400, EOZ × EDU0, and EOZ × EDU400). From this pool of fallen needles, random samples were taken in order to prepare 12 samples per treatment. A customized plastic mesh bag of standardized size, which was same for all bags, was prepared for each sample. A standardized amount of needles was added in each litter bag, and the bag was closed using stainless steel staples. While authors remember that the initial amount of needles was standardized to 1.5 g, no record of this value could be found in the experimental notes. However, only the relative difference among treatments was of interest in this study. On September 7, 2016, the prepared litter bags were implanted in soil (8 cm depth). Four experimental plots were established, scattered across the experimental forest, under ambient conditions. Three litter bags of each treatment were implanted in each plot, with a distance of 20 cm × 20 cm between litter bags. The experimental design was fully randomized in each plot. Tea bags were retrieved from all the plots 74 days after the implanting (20 November), i.e., at the end of the snow-free period. The needles of each tea bag were dried in an oven (60 °C) to constant mass, and the mass was measured.

Biomass
At the end of the experiment, all plants were harvested (October 2, 2016). Each root system was retrieved as an intact unit (no root losses), followed by removal of soil particles and alien roots by hand and gently washing with tap water. Each organ (root, needles, and stem + branches) was dried in an oven to constant mass at 60 °C, and the shoot (needles + stem + branches) and total plant (shoot + root) dry mass was calculated. The ratio of root mass to shoot mass (R/S) was also calculated.

Data analyses
The threshold of statistical significance for all the tests was set a priori to an alpha level of 0.05. Because some of the data distributions did not follow a Gaussian pattern, the data of each response variable were transformed with a Box-Cox power transformation, following a method described previously (Agathokleous et al. 2016b). While photosynthetic pigments were evaluated three times in total, a dependentsample design with time being within effects factor was not employed because measurements were conducted with different kinds of needles from different years. Therefore, data of photosynthetic pigments, gas exchange, biomass (and ratio of total number of buds per plant at harvest to branches dry mass at harvest), and litter mass following retrieval from soil were analyzed with a General Linear Model (GLM) adjusted with Method I Sum of Squares of Overall and Spiegel (Overall and Spiegel 1969;Howell and Mcconaughy 1982). Ozone and EDU were fixed factors and plot was random factor. This "stepdown analysis" is ideal for balanced analyses of variance (equal n, orthogonal), and the specification of the effects follows a proper hierarchical order in the model statement. Growth traits were analyzed with a three-factor GLM, with time within-effects factor and O 3 and EDU between-effects factors. Time factor had 8, 8, 6, and 3 levels for plant height, crown span, stem diameter, and number of buds, respectively. Posthoc multiple comparisons, where needed, were examined with Bonferroni test. We also calculated the Cohen's delta (δ) for AOZ-EDU0 versus EOZ-EDU0 for both the number of buds and the biomass of branches per plant, as described by Agathokleous et al. (2016b). Processing of data and statistical analyses were done with EXCEL 2010 (Microsoft, Redmond, CA, USA) and STATISTICA (v.10, StatSoft, Tulsa, OK, USA).

Visible foliar injury
Visible foliar injury was monitored in the second growing season (2015), in which almost all physiological measurements were also conducted. Premature senescence occurred in 33.3% of the EOZ × EDU0 plants in early June; none of the plants of other treatments exhibited senescence at this point. At this point of time, all the EOZ × EDU0 plants exhibited O 3 -induced visible foliar injury (e.g., Fig. 1), affecting on average 28.3% ± 6.5% of the foliage, and the injury was low to moderate. Conversely, 50% of the EOZ × EDU400 plants exhibited O 3 -induced visible foliar injury, and the injury was negligible, affecting only on average 1.0% ± 0.4% of the foliage. Premature senescence of the few needles of these EOZ × EDU400 plants began in early July. At this point of time, on average 2.8% ± 0.9% and 48.4% ± 10.0% of the foliage of EOZ × EDU400 and EOZ × EDU0 plants was injured (Fig. 1). By August 28, the premature senescence had stopped in O 3 -injured plants because injured needles had already been shed; the normal senescence in other plants had not started before mid. October. No O 3 -induced visible injury was observed in AOZ plants throughout the experiment.

Leaf photosynthetic pigments
In October, 2014, O 3 was not a significant main factor for any of the photosynthetic pigment traits (Figs. 2a-d, 3a-b). EDU was a marginally significant main factor for chlorophyll b (Fig. 2b), chlorophyll a + b (Fig. 2c), and chlorophyll a to b ratio (Fig. 3a). In particular, saplings treated with EDU400 had 43.9% and 32.0% higher chlorophyll b (Fig. 2b) and chlorophyll a + b (Fig. 2c) content, and 39.0% lower chlorophyll a-b ratio than plants treated with EDU0. Similarly, the O 3 × EDU interaction was significant for chlorophyll b (Fig. 2b), chlorophyll a + b (Fig. 2c), and chlorophyll a-b ratio (Fig. 3a). However, the only significant differences in multiple comparisons among means (O 3 × EDU interaction) were between EDU0 and EDU400 in EOZ for chlorophyll b (Fig. 2b) and chlorophyll a + b (Fig. 2c). Specifically, plants treated with EDU400 had 50.2% higher chlorophyll b content (Fig. 2b) and 27.4% higher chlorophyll a + b content (Fig. 2c) than plants treated with EDU0 in EOZ.
In September, 2015, O 3 had a significant main effect on all the pigment traits (Figs. 2i-j, 3e-f) except chlorophyll a + b (Fig. 2k) and carotenoids (Fig. 2l). EOZ, relative to AOZ, decreased chlorophyll a (Fig. 2i) by 33.5% and increased chlorophyll b (Fig. 2j) by 42.8% (Fig. 2c). Furthermore, plants had 2.3 and 1.4 times larger chlorophyll a to b (Fig. 3e) and chlorophyll to carotenoid (Fig. 3f) ratios in Fig. 2 Means (± se) of photosynthetic pigments in hybrid larch plants grown in ambient ozone (AOZ) or elevated ozone (EOZ) and treated with 0 mg L −1 (EDU0) or 400 mg L −1 (EDU400) ethylenediurea (EDU) over three growing seasons. Photosynthetic pigments were evaluated in the first two growing seasons. The data were analyzed with a General Linear Model (GLM). Different letters above error bars indicate statistically significant differences according to Bonferroni posthoc test for significant interaction main effects. For nonsignificant interaction main effect, posthoc test was not performed, thus the means are marked with the same letters above error bars. The level of statistical significance was α = 0.05 AOZ than in EOZ (EDU treatments pooled). EDU, as a main factor, had a significant effect on carotenoids only: EDU400 increased carotenoid content by 20.3% compared with EDU0 (Fig. 2l). The O 3 × EDU interaction was non-significant for all the pigment traits.

Leaf gas exchange
Ozone was a significant main factor for A 380 , gs 380 , E 380 , and A max in late June, 2015 (Fig. 4). EOZ decreased A 380 , gs 380 , E 380 , and A max by 63.4, 41.4, 47.0, and 69.3%, respectively, compared with AOZ. EDU was a significant main factor for only gs 380 , which was 23.6% higher in plants treated with EDU400 than in plants treated with EDU0. V Cmax and J max were not significantly affected by any factor, and none of the gas exchange traits exhibited a significant O 3 × EDU interaction.

Plant growth
Ozone was not a significant main factor for height (Fig. 5a), crown span (Fig. 5b), and number of buds (Fig. 5d); however, it was a significant (marginally) main factor for stem diameter (Fig. 5c). EOZ caused a 9.8% inhibition of stem diameter compared with AOZ. EDU was a non-significant main factor for all the growth traits, all of which varied significantly through time. The O 3 × EDU interaction was not significant for any of the traits. For height (Fig. 5a) and crown span (Fig. 5b), none of the interactions was significant. For stem diameter, O 3 × Time interaction was marginally significant; however, there was no difference of biological relevance as to the treatments effect in the multiple comparisons (Fig. 5c). For the number of buds, all interactions were significant except O 3 × EDU (Fig. 5d). Regarding O 3 × Time interaction (significant at P < 0.05), however, there was no relevant significant difference among means of AOZ and EOZ within each time point when tested with the "strict" Bonferroni test. Regarding EDU × Time interaction, there was no significant difference between EDU0 and EDU400 either before treatments (first growing season; time point 0) or in the spring of the second growing season (time point 1) (Fig. 5d). However, at the end of the third growing season (time point 2), plants treated with EDU400 had 16.6% more buds than plants treated with EDU0. Further examination of the O 3 × EDU × Time interaction revealed that this difference at the end of the third growing season appeared because EDU400 (average in EOZ = 223.3 buds) significantly protected (+ 67.9%) buds of plants against EOZ, compared with EDU0 (average in EOZ = 133.0 buds). Moreover, EDU400 led to 31.4% more buds in plants exposed to EOZ compared to plants treated with EDU400 and exposed to AOZ. No significant difference was observed between AOZ-EDU0 and AOZ-EDU400 or between AOZ-EDU0 and EOZ-EDU400 treatments at the end of the third growing season (Fig. 5d). Neither was a significant difference observed among any treatments before the beginning of the treatments in the first growing season and one year later (spring of the second growing season). In order to understand if the number of buds was due to smaller plants overall or because of fewer buds per unit biomass, the ratio of the total number of buds per plant at harvest to the dry mass of branches at harvest was analyzed. EDU was not a significant main factor, but O 3 and O 3 × EDU interaction were significant main factors (Fig. 6). Conversely to the preceding analysis with the total number of buds per plants only, this further analysis per unit of branches biomass revealed that plants in AOZ (mean = 45.5 ± 1.1 buds, n = 6) had 38.2% fewer buds per unit of branches biomass than plants in EOZ (mean = 62.7 ± 5.1 buds, n = 6). The only significant differences among means were between EOZ-EDU0 vs. AOZ-EDU0 (− 59.9%) and AOZ-EDU400 (− 55.0%) (Fig. 6).

Plant biomass
Ozone main effect was not significant for stem + branches (Fig. 7a), needles (Fig. 7c), and shoot (Fig. 7d) biomass. However, EOZ significantly decreased root biomass by 52.1% (Fig. 7b), total plant biomass by 33.0% (Fig. 7e), and root/shoot ratio by 32.5% (Fig. 7f). EDU was a significant factor only for needles biomass (Fig. 7c) and shoot biomass (Fig. 7d), which were 54.9% and 31.9% higher in plants treated with EDU400 than in plants treated with EDU0. The O 3 × EDU interaction was significant for all the biomass traits (Fig. 7a-e) except root/shoot ratio (Fig. 7f). The same pattern was observed in the biomass traits with significant O 3 × EDU interaction: EDU400 protected plants against EOZ-induced inhibition. Specifically, EDU400 led Means (± se) of gas exchange traits of hybrid larch plants grown in ambient ozone (AOZ) or elevated ozone (EOZ) and treated with 0 mg L −1 (EDU0) or 400 mg L −1 (EDU400) ethylenediurea (EDU) over three growing seasons. Net photosynthetic rate (A 380 ), stomatal conductance (gs 380 ), transpiration rate (E 380 ), net photosynthetic rate at 1500 μmol mol −1 (A max ) in the form of CO 2 , maximum rate of carboxylation (V Cmax ), and maximum rate of electron transport rate (J max ) were evaluated in late June, 2015 (second growing season). The data were analyzed with a General Linear Model (GLM). Different letters above error bars indicate statistically significant differences according to Bonferroni posthoc test for significant interaction main effects. For non-significant interaction main effect, posthoc test was not performed, thus the means are marked with the same letters above error bars. The level of statistical significance was α = 0.05 to 2.1, 1.8, 3.3, 2.4, and 2.2 times higher stem + branches (Fig. 7a), root (Fig. 7b), needles (Fig. 7c), shoot (Fig. 7d), and total plant (Fig. 7e) biomass, respectively, compared with EDU0, in EOZ. The mean of EOZ-EDU400 was not significantly different from the means of AOZ-EDU0 and AOZ-EDU400, except for needles biomass where the mean of EOZ-EDU400 was 1.6 times higher than the mean of AOZ-EDU400 (Fig. 7c).

Leaf litter decay
Ozone, EDU, and O 3 × EDU interaction were marginally significant for leaf litter decay (Fig. 8). Leaf litter from plants exposed to EOZ for two growing seasons had 35.5% less biomass than leaf litter from plants exposed to AOZ, following soil incubation for 74 days. Furthermore, leaf litter from plants treated with EDU400 in three growing seasons had 50.8% more biomass than leaf litter from plants treated with EDU0 (Fig. 8). The O 3 × EDU interaction had a pattern similar with that of most biomass traits, where EDU400 protected plants against EOZ impact. In particular,  Mean (± se) of ratio of the total number of buds per plant at harvest to the dry mass of branches at harvest in hybrid larch plants grown in ambient ozone (AOZ) or elevated ozone (EOZ) and treated with 0 mg L −1 (EDU0) or 400 mg L −1 (EDU400) ethylenediurea (EDU) over three growing seasons. The data were analyzed with a General Linear Model (GLM). Different letters above bars indicate statistically significant differences according to Bonferroni posthoc test for significant interaction main effects. The level of statistical significance was α = 0.05 the mean of EOZ-EDU0 was 2.2, 2.4, and 2.4 times smaller than the means of AOZ-EDU0, AOZ-EDU400, and EOZ-EDU400, respectively (Fig. 8). There was no significant difference among means of AOZ-EDU0, AOZ-EDU400, and EOZ-EDU400.

The overall picture of EDU protection against EOZ-induced phytotoxicity
This research revealed significant negative effects of EOZ on hybrid larch plants, which appeared as extensive visible foliar injury, premature senescence, suppressed gas exchange, inhibited stem diameter growth, suppressed biomass production, and impaired leaf litter decay. These findings suggest that hybrid larch is susceptible to elevated O 3 , in agreement with previous studies showing that hybrid and Japanese larches are susceptible to elevated O 3 (Koike et al. 2012;Agathokleous et al. 2017Agathokleous et al. , 2021. However, this research also revealed for the first time that EDU can protect against many of the negative effects of O 3 on hybrid larch, with the potential to fully protect productivity and enhance the capacity for the next-year growth (based on number of buds). Positive effects of EDU on plants under O 3 stress have been found in various trees and other types of perennial plants Agathokleous 2017;Tiwari Fig. 7 Means (± se) of biomass traits of hybrid larch plants grown in ambient ozone (AOZ) or elevated ozone (EOZ) and treated with 0 mg L −1 (EDU0) or 400 mg L −1 (EDU400) ethylenediurea (EDU) over three growing seasons. Biomass was evaluated at the end of the experiment. The data were analyzed with a General Linear Model (GLM). Different letters above bars indicate statistically significant differences according to Bonferroni posthoc test for significant interaction main effects. For non-significant interaction main effect, posthoc test was not performed, thus the means are marked with the same letters above error bars. The level of statistical significance was α = 0.05. Note: needles biomass (C) does not include all the needles shed during the growing season, as severely-injured needles of EOZexposed plants were prematurely shed earlier than the normal senescence (

O 3 and EDU effects on photosynthetic pigments
Photosynthetic pigments are key response endpoints for assessing oxidative stress in plants (Agathokleous 2021), widely shown to be affected by O 3 (Fiscus et al. 2005;Karnosky et al. 2007;Xin et al. 2016;Tiwari and Agrawal 2018;Feng et al. 2019b;Grulke and Heath 2020). For example, a meta-analysis of O 3 effect on poplar revealed that elevated O 3 (70-106 nmol mol −1 ) and current ambient O 3 (40-50 nmol mol −1 ) decreased chlorophyll content by 34% and 13%, compared with charcoal-filtered air (7-27 nmol mol −1 ) (Feng et al. 2019b). Another meta-analysis of O 3 effects on woody plants also demonstrated that elevated O 3 (88-98 nmol mol −1 ) decreased chlorophyll a, chlorophyll b, chlorophyll a + b, and carotenoids by 17%, 20%, 17%, and 15%, respectively; chlorophyll a/b ratio was not significantly changed . A large suppression of photosynthetic pigment content and a small change in chlorophyll a/b ratio were also found in Japanese and hybrid larch tall saplings grown in different soil types and exposed to ambient air or elevated O 3 (≈66 nmol mol −1 ) over two growing seasons . Chlorophyll degradation follows palisade cells injury, resulting in more visible chlorotic flecks that are indicative of O 3 -stress symptomology in trees (Karnosky et al. 2007). Despite photosynthetic pigments exhibit a high temporal variability, it remains poorly understood how O 3 affects photosynthetic pigments throughout time. The effect of EDU on woody plants grown under either ambient air or elevated O 3 throughout time remains obscure.
The results of this study show that O 3 and EDU effects on photosynthetic pigments varied with time. When measured in long-shoot needles in the fall of the first growing season, chlorophyll b was more responsive to O 3 and EDU than chlorophyll a. EDU400 increased chlorophyll b content (50%) and chlorophyll a + b content (27%), compared with EDU0, in EOZ. Moreover, EDU appeared to affect chlorophyll b independently of O 3 , as EDU400 led to a lower chlorophyll a to b ratio (39%, relative to EDU0), which indicates that it increased more chlorophyll b content relative to chlorophyll a content. Increased chlorophyll b content (13%, relative to EDU0) due to EDU (500 mg L −1 ; applied once) was also found in Pinto bean plants (Phaseolus vulgaris L.); O 3 had no significant effect on chlorophyll b content (Agathokleous et al. 2016b). These results together indicate that enhanced chlorophyll b, an accessory pigment occurring only in peripheral light-harvesting complexes and aiding in absorbing higher photons/blue light, may be an important mechanism of adaptive response by which EDU maximizes or optimizes light harvesting capacity and protect plants. By enhancing chlorophyll b, plants expand their light-harvesting antenna to attain a higher light-absorption efficiency per leaf unit, thus optimizing carbon balance and enhancing production capacity (Eggink et al. 2001) in Japanese larch seedlings treated with different light levels (Ji et al. 2015). That chlorophyll b was sustained at relatively high levels under EOZ-induced stress when plants were treated with EDU400, allowing hybrid larch to better use light, suggests a consequently higher competitive capacity under light competition as found in Japanese larch grown with different light levels (Ji et al. 2015). Since several of these needles overwinter under snow to help plants regrow in the next growing season, this may be a mechanism to increase light absorption efficiency per leaf unit to aid plants increase their regrowth potential in the following spring. Physiologically, chlorophyll b might be sustained at relatively high levels under EOZ-induced stress when plants were treated with EDU400 in the framework of an overcompensation response to protect chlorophyll a and efficient energy transfer between chlorophyll a molecules, control the antenna size, ensure the absorbance of high energy quantities, and/or dissipate excess energy (Sakuraba et al. 2010).  (2015). The data were analyzed with a General Linear Model (GLM). Different letters above bars indicate statistically significant differences according to Bonferroni posthoc test for significant interaction main effects. The level of statistical significance was α = 0.05 Chlorophyll degradation occurs in the framework of autumn senescence in deciduous trees (Mattila et al. 2018). An increase in chlorophyll a to b ratio can also be an indicator of leaf senescence (Lu et al. 2001). Ozone did not significantly affect (either alone or in interaction with EDU) chlorophyll a to b ratio in October, 2014, and July, 2015, and reduced chlorophyll a to b ratio in September, 2015 Considering also that most of the EOZ-damaged needles were already shed by the time the needles were sampled for pigment analyses, there is a possibility that younger needles of EOZ-exposed plants had enhanced stress-coping skills with decelerated senescence to protect plants. Furthermore, as EDU had a significant overall effect on chlorophyll a to b ratio (O 3 treatments pooled) in October, 2014, EDU might have some effect on chlorophyll degradation. One speculation is that EDU might delay leaf senescence via a mode of action similar to N fertilizer (Van Heerwaarden et al. 2003;Fu et al. 2020). However, although EDU contains N, abundant literature suggests that it does not act as a N fertilizer in plants (Agathokleous 2017), indicating a need for reconsideration of its basic mode of action in protecting plants against O 3 -induced toxicity. Here, a question is also raised whether N resorption in leaves of EOZ-exposed plants (Dong-Gyu et al. 2015) was protected by EDU, maybe due to keeping leaf longevity; however, N content in leaves and resorption was not evaluated in this research. As decomposition in chlorophyll-protein complexes has a relevant role in leaf N resorption, leaves of EOZ-exposed plants might also exhibit a delay in leaf N resorption. It may be postulated that not only early-shedding, but also a delay in the decomposition of chlorophyll-protein complexes might impair leaf N resorption under elevated O 3 , which should be further studied and confirmed experimentally.
In contrast to the long-shoot needles that showed relatively smaller effects of EOZ in the fall of the first growing season, short-shoot needles in mid. summer of the second growing season exhibited stronger effects of EOZ, as indicated by 70, 78, 71, 36, and 43% lower chlorophyll a, chlorophyll b, chlorophyll a + b, carotenoids, and chlorophyll to carotenoid ratio, respectively, in EOZ than in AOZ. However, chlorophyll a, chlorophyll b, and chlorophyll a + b contents were fully protected by EDU400 against EOZ. These findings are similar with those in short-shoot needles in early fall of the second growing season, although EOZ had a relatively smaller effect than earlier in the summer, as indicated by 34 and 43% decreased chlorophyll a and chlorophyll b and a non-significantly affected chlorophyll a + b content. The lower effect compared to summer might be also explained by the fact that the most O 3 -damaged needles had already been shed prematurely (see visible foliar injury in Fig. 1). Lower chlorophyll a to b and chlorophyll to carotenoid ratios due to EOZ suggest that EOZ affected plants by decreasing more chlorophyll a than chlorophyll b (EOZ tended to decrease chlorophyll a but increase chlorophyll b) and more chlorophylls than carotenoids. While plants had 2.3 and 1.2 times lower chlorophyll a to b and chlorophyll to carotenoid ratios in EOZ than in AOZ, EDU did not protect against such changes. Such alterations in photosynthetic pigments ratios have been found in several studies with O 3 and other abiotic stresses (Knudson et al. 1977;Saitanis et al. 2001;Sayyad-Amin et al. 2016). Besides, EDU had a significant effect (EDU400: + 20% relative to EDU0) on carotenoids only. Overall, across all time points, EOZ widely affected photosynthetic pigments and pigments ratios balance, commonly suppressing the light harvesting potential, and EDU400 alleviated many of the EOZ negative effects, but not all as its effect varied with trait and time point. A set of reasons may be involved in the variation in the effects from measurement time to measurement time, including (1) temporal variation of gas exchange physiology, (2) different age/type of leaf tissue analyzed, (3) time-dependent sensitivity of plants to O 3 due to e.g. ontogenic stage and senescence processes, and (4) different O 3 doses taken up by tissues at different time points, with the potential of accumulated oxidative stress over time.

O 3 and EDU effects on gas exchange
In agreement with the EOZ effect on photosynthetic pigments in the summer of the second growing season, EOZ suppressed A 380 , gs 380 , E 380, and A max by 63, 41, 47, and 69%, respectively, compared with AOZ, but not Vc max and J max . However, EDU400 only increased (24%) gs 380 , compared with EDU0 (O 3 treatments pooled), and had no significant effect on other gas exchange traits. These findings may indicate a conservative strategy of coordination between V Cmax and J max to counteract photoinhibition in the case of limiting carboxylation at the cost of maximizing photosynthesis at light-limiting conditions (Walker et al. 2014). These results also suggest that EDU mode of action against O 3 damage was rather related to biochemical processes. However, it should be noted that the mode of action of EDU appears to be not only species-specific but also cultivar-or genotypespecific (Singh et al. 2015;Jiang et al. 2018;Pandey et al. 2019;Fatima et al. 2019).

O 3 and EDU effects on growth and productivity
In this experiment, EOZ had only few effects on hybrid larch growth, which mainly appeared as inhibited stem diameter growth. Such limited effects of elevated O 3 are in agreement with a recent study with Japanese larch saplings grown in the same soils and exposed to O 3 in four growing seasons . These findings highlight that experiments with woody plants grown in fertile/fertilized soils, i.e., the vast majority of studies in the literature of O 3 effects on plants (see also Körner 1995Körner , 2006, may overestimate effects of elevated O 3 , especially for larches that are naturally grown in nutrient-limited conditions (Schulze et al. 2019). However, EOZ significantly inhibited the total number of buds per plant in the third growing season (but not earlier), an effect that was minimized by EDU400. This novel finding that EOZ reduces the total number of buds in woody plants after multiple growing seasons, suggests that EOZ can adversely affect hybrid larch plant growth and competitive capacity in the next growing season. However, this effect could not be observed in the first growing season, suggesting that experiments lasting only one growing season may not be ideal to depict such effects. In support to this, other experiments with willow (Salix sachalinensis Fr. Schmidt) and Japanese larch exposed to O 3 conditions similar with the ones in this study -but for only one growing season-did not find significantly changed number of buds due to O 3 (Agathokleous et al. 2016a. To understand if the decreased number of plants was due to only smaller plants in EOZ (with no affected number of buds per unit of size), we standardized the total number of buds per plant (at harvest; the measurement conducted in the third growing season) according with the biomass of branches per plant (at harvest). Interestingly, this analysis revealed a significantly increased number of buds per unit of branches biomass due to EOZ, an effect that was mediated by EDU400. This response of EOZ-exposed plants treated with EDU0 may indicate an ecological strategy for investment to more buds per unit biomass to enhance the growth potential and competitive capacity in the next growing seasons. We did not find a paper reporting this kind of response to O 3 , regarding investment to buds, thus, this is the first study documenting this mechanism under O 3 stress. Here, it should also be mentioned that the number of buds exhibited a larger effect (δ = 10.11) relative to branches biomass (δ = 2.82). This result cancels out the possibility that branches biomass exhibited a stronger effect of EOZ compared to the number of buds. In larch, the differentiation and development of buds as well as their responses to environmental changes have long been studied mainly for the artificial management of flowering promotion and the ecological managements under climate change (e.g., Heitmuller and Melchior 1960;Owens andMolder 1979, 2011;Sieber and Luscher 1995;Powell 1995;Ziauka and Kuusienë 2006;Owens 2008), suggesting that the plastic behavior of bud formation in response to environmental changes would be a unique characteristic of this larch. Our study further confirms this, using O 3 as a different type of environmental stress (air pollutant), and provides a description of processing responses of individual leaves that may provide a perspective for further research to better understand the physiological mechanisms of such plastic responses of bud formation.
As a result of suppressed light-harvesting potential and photosynthetic capacity, EOZ inhibited root and total plant biomass, but not stem + branches, needles, and shoot biomass. These phenomena might be partly explained by the non-fertilized soil, simulating natural conditions, which led to limited effects of EOZ on plant aboveground growth. As a result, EOZ significantly decreased R/S ratio, which indicates that root was more affected than shoot in these non-fertilized soils, and such an effect may suggest negative effects for woody plant cultivation in the long term as a disproportional biomass of roots may not support well aboveground growth. The sensitivity of root biomass to elevated O 3 has been shown by many studies Pan et al. 2019Pan et al. , 2020Mrak et al. 2019Mrak et al. , 2020Tisdale et al. 2021), commonly with fertile soils, with a reduced allocation to the roots (compared to the shoots) occurring in many cases (Grantz et al. 2006). EDU400 did not mediate the EOZinduced effect on the R/S ratio; however, it minimized EOZ effect on stem + branches, root, needles, shoot, and total plant biomass. These results suggest that EDU400 could minimize EOZ-induced biomass losses, but it offered a higher degree of protection to shoot (compared to root), even if it was applied to plant root (soil drench). Furthermore, biomass allocation responded to EOZ even under EDU400, but EOZ-induced accelerated leaf senescence was effectively circumvented by EDU400. The preferable biomass allocation into shoot under O 3 stress, irrespective of EDU, may suggest that plants perceived O 3 as a signal of preferential accumulation to shoot, even in the presence of EDU400. Enhancement of shoot growth by O 3 was also found in studies with different tree species . Based on the results, EDU might have not prevented O 3 uptake in the plant tissue interior, but might have enhanced the stability of some function regarding leaf senescence. EDU might contribute to maintaining membrane integrity, which increases tolerance to O 3 (i.e., preventing leaf senescence), as it has been reported for volatile organic compounds function or drought-enhanced tolerance to heat (Zhang et al. 2019). However, further studies are needed to test these hypotheses.

O 3 and EDU effects on leaf litter decay
Making a further step, this research evaluated for the first time the potential effect of EDU on the decay of leaf litter of EOZ-exposed plants. Although the incubation of litter bags was relatively short, the results showed a significantly lower mass of leaf litter following a 74-day decay when litter was exposed to EOZ (compared to AOZ), suggesting that leaf litter from EOZ-exposed plants exhibited a higher/faster decay than leaf litter from AOZ-exposed plants. This study did not aim to investigate the underpinning mechanisms of the difference in decay but used the amount of decayed biomass as an index of health that may translate to effects on the community and ecosystem. Therefore, the reason(s) for the higher decay of leaf litter from EOZ compared to leaf litter from AOZ remain(s) unknown. For example, it might be even due to difference in the type of leaf litter as EOZexposed plants (in particular those of EDU0 group) lost most of the severely O 3 -injured needles before litter collection. Another explanation is that EOZ-EDU0 leaves of hybrid larch might have higher litter N (Dong-Gyu et al. 2015), although not analyzed in this study. A few studies showed before that O 3 can alter the decomposition process using various plant species (Boerner and Rebbeck 1995;Baldantoni et al. 2013), albeit there is a discrepancy in the findings (Parsons et al. 2004(Parsons et al. , 2008Lindroth 2010;Vitale et al. 2019). Ozone slowed down (Parsons et al. 2008;Baldantoni et al. 2013;Fu et al. 2018), enhanced (Fenn and Dunn 1989;Boerner and Rebbeck 1995;Parsons et al. 2008;Fu et al. 2018), or did not significantly affect (Boerner and Rebbeck 1995;Fu et al. 2018) biomass loss of leaf litter. These differences can be attributed to various factors such as species, incubation duration, and O 3 exposure characteristics. Similar to our findings, more mass remained in leaf litter of yellow poplar (Liriodendron tulipifera L.) exposed to ambient or elevated O 3 (relative to charcoal-filtered air) in the first 10 months of incubation, but these differences no longer existed by the end of the year (Boerner and Rebbeck 1995). Likewise, elevated O 3 decelerated the leaf litter (Ginkgo biloba L.) decomposition early in the incubation period and slightly accelerated it in later stages of the incubation (Fu et al. 2018). While the difference in decay rate may be temporal, it may have implications to various ecological processes, including nutrient cycling processes. That EDU400 fully protected from the EOZ-induced alteration of leaf litter decay suggests for the first time that positive effects of EDU on plants may have further protection of soil community/ ecosystem.

Conclusion
To conclude, EOZ led to various adverse effects on the hybrid larch plants, spanning from the leaf level (pigments, gas exchange) to the level of individual (biomass, number of buds) as well as leaf litter decay, and EDU400 minimized most of the EOZ-induced effects. A novel ecological strategy of plants under O 3 stress was also revealed, in the framework of which hybrid larch saplings increased the investment toward enhanced development of buds per unit of branches biomass in EOZ, an effect that was mediated by EDU400. This strategy may be attributed to an effort of plants to increase their growth potential and competitiveness in the next growing season. The protection of plants by EDU400 was mainly due to a biochemical mode of action. All in all, EDU400 can be considered as a tool to improve this heterosis-displaying hybrid larch with superior characteristics under elevated O 3 conditions in the future.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.