Delayed Inhibition of Photosynthetic Performance—Three Linden Species in an Urban Environment

: The presented study concerns the delayed response of photosynthetic performance to summer drought. Increasing drought events are thought to be one of the most detrimental and expensive natural disasters. For that reason, we studied the condition of urban trees, which are especially susceptible to water deﬁcit. Three species of linden trees ( T. tomentosa, T. × europaea and T. cordata ) growing in the city area of Krakow, Poland, were investigated for two consecutive years to compare their response to a precipitation deﬁcit. For the ﬁrst time, the physiological traits of linden urban trees were compared in years with and without a natural summer drought. For evaluation of the photosynthetic performance, we used measurements of the chlorophyll a ﬂuorescence concerning the OJIP transients and derived parameters. The contents of photosynthetic pigments, such as chlorophylls and carotenoids, were analysed, and stress-sensible parameter ratios, such as chl a/b and chl/carot, were calculated. The most common method of research conduction assumes examination close to the presence of a stressor. We chose an innovative method of prolonged investigation after relief from the stress conditions. While there was no evident reaction of PSII to meteorological drought during the stress occurrence, all of the investigated species of linden revealed delayed reaction of PSII to precipitation deﬁcit in summer. Changes in PSII were connected with the decrease of active reaction centres (RCs), which was evident in the increase of the ABS/RC, DI 0 /RC and TR 0 /RC, and the decrease of ET 0 /RC. T. × europaea demonstrated a more intense decrease of certain parameters during late senescence compared with other species. Despite the close relation of T. cordata and T. × europaea , species showed different physiological traits during senescence after a summer drought. Our results underline the need for monitoring tree conditions not only during stress but also in the wider perspective.


Introduction
Drought has a broad range of definitions, which includes its drivers, sequence, and consequences of such events. 'Meteorological drought' is connected with the lack of precipitation, 'hydrological drought' with poor soil moisture, 'agroecological drought' concerns lower crop productivity and 'socioeconomic drought' is an effect of combined hydrological and agroecological drought [1,2]. Due to its complexity and lack of a unique definition, it is difficult to establish the exact beginning of a drought event and, thus, its duration [3]. This matter is highly important in many fields of study as increasing drought events are considered to be the most detrimental and expensive of natural disasters [4].
Due to climate change and global warming, droughts are predicted to be a more frequent and severe phenomenon connected with both widespread warming and precipitation decline [5]. The severity and frequency of drought are, therefore, linked with progressive climate change and also with thee specific features of a particular environment [6]. One of the specific environments where drought occurs with high intensity is urban areas. The and length) around day 160 of the year (June), and they are fully mature from that moment. In our work, the leaves collected in June were treated as fully developed.
From mid-September (around the 260th day of the year), changes related to aging were noticed. While the main leaf fall was observed from mid to late October, single leaves fell even at the beginning of November [21]. At this time, the disappearance of chlorophyll pigments and carotenoids were observed, with the former decomposing faster. For this reason, the car/chl ratio increased towards the end of the season. With reference to this information, in the presented research, the leaves collected in September were described as early senescence and those collected in October as late senescence.
In addition, changes in the chlorophyll a fluorescence were found to be noticeable during the growing season [22]. Five oak species showed an increase in the maximum quantum yield of the primary photochemistry of PSII and the performance index on an absorption basis at the beginning of the growing season, followed by stabilisation of these parameters and their reduction in October. The shape of the fluorescence transient curve also changed both at the beginning of leaf development and during senescence.
Somewhat differently, in the case of Fagus sylvatica L., the authors did not record changes in the maximum quantum yield of PSII (Fv/Fm) during the growing season [23]. The influence of artificial drought on this parameter was also not demonstrated. Interestingly, withholding irrigation decreased the PI total parameter, but after a few weeks of the duration of drought and more than 3 weeks later, the reduction in the intensity of photosynthesis started; the drought started in June and a decrease in the PI parameter was recorded in mid-August.
The information presented shows that the photosynthetic activity and chlorophyll a fluorescence of the leaves change during the growing season. Therefore, in the presented work, we demonstrate the dynamics of changes in the year with average precipitation and in the year with precipitation deficit.
The presented study concerns urban trees, which, as quoted, are particularly exposed to drought stress. It is known that linden trees show an anisohydric reaction to drought. Such behaviour maintains relatively high transpiration rates and carbon assimilation under water deficit in contrast to isohydric reactions. Such adaptation results in delayed growth reduction, evident in the year following drought events [24]. Such a reaction would not limit the carbon reactions themselves, while the CO 2 would be assimilated through stomata.
However, due to observed precipitation limitation (see results) and high temperatures, we assumed the disturbance of photosynthesis light reactions. As presented by Gillner et al. [25] anisohydric species, such as Corylus or Tilia, can maintain high efficiency of photosystem II during drought. While there was no evident reaction of PSII to meteorological drought in summer, we decided to evaluate if there were any further consequences of a summer drought event for PSII functioning. The study aimed to investigate how a precipitation deficit in summer affects the photosynthetic performance at two stages of senescence. To meet this objective, we chose an innovative method of prolonged investigation after the relief of the stress conditions. We formulated two main hypotheses: The PSII reaction can be delayed and revealed few weeks after the stress factor appearance.

2.
The species differ in the course of aging during normal and drought years.

Study Sites and Objects
The study was performed in Krakow (19 • 57 E, 50 • 03 N), the second-largest city in Poland, where intensive urbanisation has occurred in the last 50 years. Due to rapid city development, an increase in the impermeable surface coverage, which contributes to surface runoff and the worsening of water management, has been observed [26]. The experiment was conducted for three years (2013-2015) in ten different locations (two consecutive years for every location). The investigated trees were located within the urban area of Krakow. During the years of the experiment, the examined trees were not irrigated, and the Municipal Management of Urban Green Areas in Krakow did not conduct any irrigation program.
The soil pH KCl measured for every location ranged from 5.5 to 7.3 in the spring and from 4.9 to 7.2 in the autumn (data not published). The electrical conductivity of the soil in the investigated location ranged from 160 to 483 µS·cm −1 in the spring and from 82 to 402 µS·cm −1 in the autumn of 2014 (data not published). Taking into account the results of the soil analysis, we assumed that there was no excessive salinity or excessively high pH, which could be separate stress factors. However, climatic and growth environment features in urban areas are highly variable.
For that reason, we chose strict criteria for the selection of the studied trees. To ensure the credibility of the obtained results, we assumed three selection criteria. All of the examined trees were from (i) one genus-to avoid inter-generic differences; (ii) a sufficient number of trees in one location  individuals, planted at the same time)-to ensure a representative research sample; (iii) with a similar age of individuals-to avoid differences between young and mature trees (all selected plants were up to 15 years, planted on study sites no less than 2 years before the study).
The investigated trees were located along roads or in the near vicinity, planted as a strip with trees of the same age and species. Our research included 93 trees from the Tilia genus: 16 individuals of T. cordata (Mill.), 36 of T. tomentosa (Moench.) and 41 of T. × europaea (for more details see Table S1 and Figure S1 from Supplementary Materials). All of the investigated species have similar distributions. T. cordata and T. tomentosa have European origins; however T. tomentosa originates from drier regions. T. × europaea is a hybrid species of T. cordata and T. platyphyllos, which are both native in Poland.
The morphology of the species leaves is different, concerning mostly the size, thickness and hair coverage. Among the investigated species, T. cordata has the smallest leaf blades, and they are also less hairy. T. tomentosa has the largest leaf blades with dense hairs below [27]. Due to those morphological features, and its origin in relatively dry regions, T. tomentosa is thought to be better suited for urban areas compared with other species. However, while T. cordata is widely planted in European cities its suitability, as well as its hybrids is more often perceived [28,29].

Climatic Data
Detailed climatic data (temperature) were collected from the climatic computer located outside of the University of Agriculture's glasshouse. Precipitation measurements for Krakow and data comparing each month to precipitation norm based on measurements from the 1971-2000 period were gathered from Bulletins of the Hydrological and Meteorological Service [30][31][32]. The temperature and precipitation are presented in a Gaussen-Walter diagram.
This kind of climatic data presentation allows the estimation of a precipitation deficit or its excess, according to the specific presentation of the temperature and precipitation record. In the original presentation of this diagram, 10 • C corresponds to 20 mm of precipitation. However, as reported by Treder et al. [33] for the climatic conditions of Poland, 40 mm (or even 45 mm) would be more appropriate to counterbalance the temperature of 10 • C. For that reason, we applied the modifications proposed by Łukasiewicz [34] with a ratio of 1 • C to 4 mm between the mean temperature of a particular month and the amount of precipitation.

Leaves Collection and Analysis
Plant material was collected in the morning hours, from the outside part of the crown, and the eastern side of the tree. Such procedures eliminated the differences connected with non-uniform illumination of different parts of the crown [35]. Leaf samples were gathered from five branches parts from each tree (about 25-50 leaves from one individual). Leaves were carefully protected from drying and mechanical damage and immediately transported to the laboratory for further analysis.
All parameters were compared in three terms of collection. Leaves collected in June were described as "full development" (according to Lichtentaler and Babani [36]); leaves collected in September were "early senescence"; and, in October, the leaves were "late senescence". Only leaves visually recognised as healthy (without chlorosis, necrosis or foraging pests), were taken for further investigation. In early senescence, only green leaves were chosen for analysis, where possible. In the late senescence stage, leaves have visible symptoms of natural senescence, such as discoloration and marginal necrosis.
The parameters of a given species in a selected period of a dry or normal year were considered as a treatment (three species × three developmental phases × 2 years with different precipitation).

Chlorophyll a Fluorescence
For every location of tree growth, fluorescence measurements were made with the light intensity of 3000 µmol·m −2 ·s −1 , for 10 replications, in each of the studied months. Prior to measurement, the leaves were dark-adapted for 30 min. All fluorescence measurements were obtained using Handy PEA (Plant Efficiency Analyser), Hansatech (UK). We analysed the OJIP fluorescence transient, called the JIP-test, where O is the fluorescence origin (first measured minimal level), J and I are intermediates and P is the maximum level of the fluorescence curve. We also calculated some of the JIP-test parameters, which are connected with a particular phase of electron transport.
All calculations were made according to formulas from Stirbet and Govindjee [37] and Stirbet et al. [38]: F 0 (the minimum chlorophyll a fluorescence), Fm (the maximum chlorophyll a fluorescence after dark adaptation), Fv/F 0 (the ratio of the photochemical and non-photochemical processes in photosystem II (PSII) (the maximum efficiency of the photochemical processes of PSII), Fv/Fm (the maximum quantum yield of PSII photochemistry), Tfm (the time to reach the maximum chlorophyll fluorescence), Area (the area above the OJIP transient and Fm line), PI total (the performance index for energy conservation from absorption to the reduction of PSI acceptors), RC/ABS (the PSII reaction centre per absorbed photon flux) ABS/RC (the absorbed photon flux per PSII reaction centre (RC) or apparent antenna size of an active PSII), TR 0 /RC (the maximum trapped exciton flux per active PSII), ET 0 /RC (the flux of electrons transferred from the primary electron acceptor (QA) per active PSII reaction centre) and DI 0 /RC (the flux of energy dissipated in processes other than trapping per active PSII reaction centre). We evaluated the parameters concerning specific energy fluxes expressed per active PSII reactive centre. All of the measurements concerning the OJIP-test were taken for 10 replications for the mixed leaf sample from one location (for more details see Table S1 from Supplementary Materials).

Chlorophylls and Carotenoids Content
Assimilation pigment contents were analysed using chemical procedures. Pigment extraction proceeded with 80% acetone, according to procedures described by Wellburn [39]. Analyses were made in four replications for every location of tree growth in each of the studied months (for more details see Table S2 from Supplementary Materials). The measurements were taken on a Hitachi U-2900 (Tokyo, Japan) spectrometer with 663, 646, and 470 nm wavelengths. In addition to the chlorophyll sum, we also calculated the parameters of the chlorophyll a to b ratio (chl a/b) and chlorophyll sum to carotenoids (chl/carot), which are considered drought-stress-sensitive indicators [40,41].

Statistical Analysis
All statistical analyses were performed using Statistica v13.1, StatSoft software (TIBCO Software Inc., 2017, Palo Alto, CA, USA) from Statistica (data analysis software system, version 13.). The statistical differences were performed with the Kruskal-Wallis method and Dunn test. The test was performed separately to distinguish differences between (I) early senescence and full development and (II) late senescence and full development. Significant differences were established with p < 0.05.
The method of cluster analysis was also used to group the examined treatments based on the measured features. Prior to analysis, all data were standardised. The Ward's method and Euclidean distance were used to group treatments. The dendrogram was created, and a line marked the cut-off point that allowed for the identification of the clusters. The treatments with similar characteristics were classified into one cluster.

Climatic Conditions
While the temperature course in studied years is similar, the precipitation distribution in the same month differs greatly between years ( Figure 1). In 2013, precipitation deficits occurred in April, July and August and later in October. In 2014, from June to July, there was only a minor precipitation deficit according to the Gaussen-Walter diagram, which we recognised as negligible. In the third year of study, the months from June to August were characterised by insufficient precipitation. A similar classification of consecutive months in the years of study was observed in relation to the precipitation norm for Poland based on measurements from the 1971-2000 period ( Figure 2).
The years 2013 and 2015 can be characterised as years with summer precipitation deficits, and 2014 can be classified as normal according to Figures 1 and 2. For that reason, the measured leaf parameters are presented in groups for each species of linden trees, comparing the results from years with a precipitation deficit (D) and with precipitation in the summer described as normal (N).

OJIP-Transient
All of the investigated species showed a lower fluorescence transient during late senescence in the years with precipitation deficit in comparison to other terms of analysis ( Figure 3). Meanwhile, for T. cordata, while the O-K phases were similar to other terms of analysis, the J-I and I-P phases were lower in all of the investigated species. In T. cordata, the early senescence (D) fluorescence transient (especially the J-I and I-P phases) were lower than in other terms. For T. × europaea, the late senescence (N) curve was also lower than other terms in a normal year, which was not observed for other species. In addition, the chlorophyll a fluorescence transient for full development (D) in T. × europaea was lower than in a normal year. In general, the F P (Fm) was lower for T. × europaea, reaching a maximum of about 2000, while for T. cordata and T. tomentosa, it was about 2300 and 2500, respectively. What is interesting for T. tomentosa is that, apart from late senescence (D), all OJIP curves had a similar course.

OJIP Test and Photosynthetic Pigment Parameters during the Normal Year
As presented for all investigated species, during the year without a precipitation deficit, the parameters in early senescence did not differ drastically from those in full development. F 0 and the chl a/b ratio were significantly lower in all species during early senescence, and DI 0 /RC was lower for both T. tomentosa and T. cordata. The Fv/F 0 , Fv/Fm and chl/carot ratios were significantly higher during early senescence than in the full development stage for T. cordata and T. tomentosa. In the late senescence Area, PI total and chl a/b were lower than in the full development stage for all species, and TR 0 /RC and ABS/RC were higher for T. cordata and T. × europaea. An increase in DI 0 /RC at the late senescence stage was observed only in T. × europaea. The chlorophyll content did not differ significantly in the investigated terms for all of the investigated species.

OJIP Test and Photosynthetic Pigment Parameters during the Year with a Precipitation Deficit
The differences between parameters in the year with a precipitation deficit were much more evident for both stages of senescence. PI total was the only parameter that decreased in early senescence for all of the investigated species. The Fv/F 0 , Fv/Fm, RC/ABS, ET 0 /RC and chl/carot ratios decreased in early senescence in comparison to full development for T. × europaea and T. tomentosa. For all of the investigated species, Tfm was higher during early senescence than in full development, and F 0 , ABS/RC, TR 0 /RC and DI 0 /RC were higher in T. × europaea and T. tomentosa.
In late senescence, the Fm, Area, Fv/F 0 , Fv/Fm, RC/ABS, ET 0 /RC, PI total , chl sum and chl/carot ratios decreased in all taxons. The DI 0 /RC increased greatly in late senescence in comparison to full development for all of the investigated taxons. The Tfm parameter was higher in late senescence for T. cordata and T. tomentosa but lower in T. × europaea in comparison to full development (Figure 4).

Cluster Analysis
A cluster analysis made on the basis of all tested parameters separated as 18 treatments (three species × three developmental phases with two categories of precipitation deficit N and D) into four different clusters with cut-off point of a Euclidean distance 7.5 ( Figure 5). Cluster 1 included all species analysed during full development and early senescence from category D (precipitation deficit). The second cluster contained all combinations with category N except T. × europaea_LS_N. This treatment was classified into cluster 3 along with T. cordata_LS_D and T. tomentosa_LS_D. One treatment was assigned to the last fourth cluster (T. × europaea_LS_D). Figure 6 was prepared based on the Euclidean distance of a given development phase from full development for each species and precipitation deficit category, separately. It can be noticed that T. × europaea analysed during late senescence was classified separately than other species late senescence leaves from the specific year (D or N, respectively). Late senescence leaves during the year without precipitation deficit of T. tomentosa and T. cordata were included in cluster 2 and the case of T. × europaea in cluster 3. In addition, during the year with precipitation deficit, the late senescence leaves of T. tomentosa and T. cordata were separated in a different cluster than T. × europaea (clusters 3 and 4, respectively).
In addition, in Figure 6, the distance from the full development phase to early senescence and late senescence during N year is small, except for T. × europaea _LS_N. The distance between the full development phase and ES of plants from the year with a rainfall deficit is slightly higher. Much larger is the distance from full development to phases/treatments included in cluster 3, which emphasizes their difference compared to the full development (changes during senescence). However, the largest changes from FD to LS are observed in the case of T. × europaea during year with a precipitation deficit, where the greatest distance was recorded. Fully developed leaves were set as 1, and early and late senescence presented in relation to fully developed leaves. Statistical analysis was performed separately to distinguish statistical differences between early senescence and full development (first mark), and late senescence and full development (second mark; − means there were no statistical differences, * means there were statistical differences according to the Kruskal-Wallis method and Dunn test with p < 0.05).

Discussion
Drought is a complex phenomenon that can be measured and considered in many ways. For that cause, it is difficult to indicate the exact beginning and duration of such an event [3]. In addition to in-depth techniques testing the soil moisture and evapotranspiration rates, there are also simplified methods that enable the determination of drought events under limited climatic data [33].
To estimate the drought event appearances during years of study, we used the Gaussen-Walter climatic diagram, which facilitates estimation of the precipitation deficit

Discussion
Drought is a complex phenomenon that can be measured and considered in many ways. For that cause, it is difficult to indicate the exact beginning and duration of such an event [3]. In addition to in-depth techniques testing the soil moisture and evapotranspiration rates, there are also simplified methods that enable the determination of drought events under limited climatic data [33].
To estimate the drought event appearances during years of study, we used the Gaussen-Walter climatic diagram, which facilitates estimation of the precipitation deficit with specific chart plotting. According to this method, each group of investigated trees were examined during one year with a summer (months VII-VIII) precipitation deficit, and one with the summer precipitation classified as normal or moist (according to the presented climatic diagram and analysis of precipitation norm- Figures 1 and 2).
Stressful conditions, such as drought, severely affect photosynthesis in all its phases and at many levels. Drought can cause the degradation of chlorophyll and thylakoid membranes, hence, disturbing photosystem functioning and mainly the light-absorbing efficiency [42]. Chlorophyll a fluorescence is one of the effective measurements that provides information about PSII functioning. It has been proven suitable in measurements of tree reactions to specific stresses, such as drought [43,44].
In most of the available research, the examination of stress effects was conducted close to the presence of the stressor. We decided to apply the innovative way of prolonged investigation, after the relief from stress conditions. We did not observe an evident reaction of investigated trees PSII during stress occurrence. Therefore, we decided to evaluate if drought events can influence the autumnal senescence course.
Autumnal senescence is a programmed process that enables nutrient relocation before the dormancy period. Aging patterns are strongly connected to the plant species, and environmental conditions can alter the natural course of those reactions [45]. It is known that, during the senescence of deciduous tree leaves, the content of photosynthetic pigments, such as chlorophylls and carotenoids, decreases; however, the degradation processes take place at different rates [46]. The chlorophyll degradation traits were different among investigated species.
It has been proven that deciduous tree species (maple, oak and beech) vary in the chlorophyll degradation trend during senescence [47]. The mentioned authors also assumed that species differed in the onset of senescence processes. As presented for four different tree species (S. aucuparia, A. platanoides, B. pendula and P. padus), individual leaves maintained relatively high chlorophyll content for the whole autumnal senescence period, estimated at 2 months [48]. We underline that the fast decrease of chlorophyll was observed only about a week before leaf abscission. A similar reaction was observed in the investigated linden trees (Figure 4).
The chlorophyll sum decreased significantly in the year with a summer precipitation deficit during late senescence for all of the investigated species in comparison to full development and in early senescence for T. × europaea. However, it did not decrease in early or late senescence for any species after the summer with 'normal' precipitation. Such a result indicates disturbances in the natural chlorophyll degradation in the year with a summer precipitation deficit. This observation indicates a delayed reaction of physiological traits to summer natural drought.
It was also presented, for Tilia cordata growing in Poland, that the content of Chl a was lower in autumn in comparison to spring and summer, while the Chl b and carotenoid contents did not change significantly [49]. This is connected with the faster senescencerelated degradation of chlorophyll a over chl b and carotenoids, and this is consistent with our results during the year with normal precipitation (Figure 4). The chlorophyll to carotenoid ratio, which was decreased during senescence (D), was unaffected or even higher during senescence (N) in comparison to full development. This indicates a more drastic reaction during senescence in the year with precipitation deficit in comparison to the 'normal' year.
The mentioned changes in the rate of chlorophyll pigment degradation were not so evident in the presented chl a/b ratio. This is caused by the decrease in the chl b content in early and late senescence (D). For that reason, the difference between chl a/b during senescence (D) and full development (D) was insignificant in most cases. The general degradation of both chlorophylls was observed during year with a precipitation deficit. The reduction of total chlorophyll in LS compared to FD in the year with a precipitation deficit was about 58% in T. cordata, 73% in T. tomentosa and 95% in T. × europaea, which emphasizes the differences in reactions to drought between these species.
Under natural leaf senescence, the OJIP transient became flatter, and the values of the fluorescence in the J, I, and P decreased [50]. The responses of the chlorophyll a fluorescence transient (OJIP curve) were changed after drought in summer (D), which was especially observed during late senescence (Figure 3). The reduction in the I-P phase was observed under the drought stress of Fraxinus ornus. It correlates with active plastocyanin (PC) and P700, and reflects the ability for reducing final acceptors beyond PSI.
In addition, the reduction in the J-phase is connected with the down-regulation of PSII photochemistry as a photoprotective mechanism observed in Fraxinus ornus under drought [43]. In the presented results for all investigated Tilia species, both photosystems were disturbed in the late senescence (D) while, in the 'normal' year, this occurred only in T. × europaea (decrease in the J-phase and I-P phases). These findings are contrary to drought-stressed T. cordata growing in controlled conditions in pots [51]. In such conditions, the authors observed the increase of K and J bands of OJIP test after 28 days without plant watering.
As presented for two young (15-20 years old) deciduous tree species Acer saccharum and Quercus bicolor, a high photochemical efficiency was maintained during autumnal senescence [52]. Selig and Bohne [53] observed no differences between the control and drought-treated Tilia cordata Fv/Fm parameters. In some cases, this parameter is not sensitive enough, as was suggested previously [44,54,55]. Similar results were presented in our study for investigated linden species during the year with 'normal' precipitation in the summer (Figure 4). There was no significant decrease in the maximum photochemical efficiency of PSII (measured with Fv/Fm) during early senescence for any species. However, such a decrease was observed in the late senescence for T. × europaea, which again had a distinct reaction.
In the year with a summer precipitation deficit, there was a significant decrease of the Fv/Fm parameter during late senescence in all of the investigated species. However, the reduction of this parameter was not uniform in all the species studied, in the case of T. cordata, it was about 4%, in T. tomentosa 30% and in T. × europaea, 62%. The decrease of the Fv/Fm parameter during drought was also observed in Tilia cordata by Kalaji et al. [51]. Such a reaction was also observed in Acer saccharum leaves growing in the city park area, during autumnal leaf coloration [56]. The change from green to yellow and orange leaves resulted in the decrease of the Fv/Fm parameter from 0.75 to 0.4 and 0.25, respectively. Such a decrease of the maximum quantum yield of PSII is connected with a loss of photosynthetic capacity.
Changes in the specific fluxes were observed during natural leaf senescence. As presented for Acer platanoides leaves, during yellowing, ABS/RC, TR 0 /RC and DI 0 /RC increase and ET 0 /RC decrease was observed [50]. However, the changes in investigated Tilia leaves were more evident during senescence after natural summer drought than after a 'normal' summer. Similarly to our results, after 28 days of drought, Kalaji et al. [51] observed an increase of F 0 , ABS/RC, TR 0 /RC and DI 0 /RC as well as a decrease of Fm, Fv/F 0 , Fv/Fm and PI total in the investigated Tilia cordata leaves. An increase of DI 0 /RC is an indicator of the dissipation of absorbed photosystem energy when the photochemistry is disrupted [57,58].
This indicates that the drought stress observed in summer inhibited electron transport and stimulated the dissipation of energy in late senescence in Tilia leaves. Although Kalaji et al. [41] did not observe changes in ET 0 /RC, in our experiment, a decrease of electron transport was significant in the LS_D in all Tilia species. Absorbed radiation energy was redirected not to electron transport and photochemical reactions but to the processes of its dissipation. Furthermore, the redirection of absorbed energy was also confirmed by the decrease of Fv/F 0 parameter. This observation could be an indicator of a lower photochemical process efficiency in PSII [59]. In addition, the decrease of Fm or Fv/Fm could be the information about the reduced potential of primary acceptor QA of PSII for receiving electron due to the dissociation of the light-harvesting complex (LHC) from the PSII [51]. During the dehydration of plant tissue, the TR 0 /RC increased as a result of the increase of inactive RCs and decrease of active RCs until 30% of RWC (relative water content in tissue) [60]. The decrease of active RCs was also possible in our experiment in late senescence Tilia leaves after a summer drought.
During the year with normal precipitation, T. cordata showed no changes in the ABS/RC, TR 0 /RC and DI 0 /RC parameters during senescence, while, in the case of T. × europaea, significant changes were observed. In addition, the T. cordata PI total parameter during aging (comparison of FD with LS) was reduced by 65%; in T. × europaea, this was reduced by 74%; and, in T. tomentosa, a decrease of 31% was observed. During the year with a deficit of precipitation, changes in the ABS/RC, TR 0 /RC, ET 0 /RC and DI 0 /RC as well as in the PI total parameters were observed in all the species studied during senescence.
However, the levels of change were not similar. The PI total decreased by 73% and 75% in the case of T. cordata and T. tomentosa, while in T. × europaea decreased by 99%. The increase in energy dissipation expressed by the DI 0 /RC parameter was 2.3 fold in the case of T. cordata, 9 fold in T. tomentosa and as much as 25 fold in T. × europaea. This clearly shows the differences in the response to summer drought stress between the studied species.
Cluster analysis is a useful tool for treatment classification. The usefulness of cluster analysis in the classification of similarity of Tilia species based on morphological features was proven [61]. However, in the presented study, the most important comparison was the reaction of different species to stress conditions (summer with precipitation deficit). Therefore, the result of cluster analysis shows treatments with similar and distant responses. Four distinct groups were separated, where treatments exposed to a precipitation deficit in summer were clearly separated from those collected during years with better rainfall distribution.
Moreover, some growth phases of plants exposed to summer precipitation deficits were separated: full development and early senescence from late senescence. This indicates significant changes in the tested parameters and deterioration of the functioning state of the photosynthetic apparatus. There was no phase separation during the better distribution of rainfall (N), except for the late senescence of T. × europaea. It is known that plants differ in their 'drought resistance' adaptations, and reactions are species-dependent [62].
The results of the cluster analysis suggest that the response to senescence and summer drought stress is different for T. cordata and T. × europaea, despite the species relationship. Swoczyna et al. [63] classified T. cordata and T. × europaea into one group based on fluorescence parameters and found that, as forest species, they are sensitive to drought stress. However, those classification also included other, not-related urban tree species.
Three species of linden were examined in the presented study, while, in the studies of Swoczyna et al. [63], various tree species were compared (Tilia, Acer, Quercus, Gleditsia, Pyrus, Platanus and Ginkgo), which likely contributed to the classification of T. cordata and T. × europaea into one group. In our research, in the late senescence phase, both species are classified into different groups, which indicates that the pattern of senescence and response to drought stress in both species is different. The greatest change from full development to late senescence was seen in T. × europaea during the year with precipitation deficit.

1.
Investigated species of linden trees, i.e., Tilia tomentosa, T. × europaea and T. cordata revealed delayed reactions to precipitation deficits. The photosynthetic performance of the investigated species was disturbed in the late senescence phase after summer with precipitation deficits. Those changes were less evident during senescence after summer without precipitation deficits.

2.
Changes in the OJIP test included mostly parameters related to the number of active reaction centres. After the summer precipitation deficits, increases of the ABS/RC, DI 0 /RC and TR 0 /RC and a decrease of ET 0 /RC were related to the decrease of active RCs.

3.
Different severities of the senescence-related reactions were observed in T. × europaea. This heterotic species demonstrated a more intense decrease of certain parameters during late senescence (N) compared with other species (i.e., Fm, Fv/Fm and Fv/F 0 ). In the years with precipitation deficits, those differences were more evident. For example, compared to full development, the decrease in the maximum fluorescence (Fm) during late senescence (D) was about 30% for T. cordata and T. tomentosa and about 75% for T. × europaea. Despite the close relation of T. cordata and T. × europaea, the species showed different physiological traits during senescence after precipitation deficits in summer.

4.
The presented results underline the need for monitoring tree conditions not only during stress but also in the wider perspective, concerning a prolonged time of measurement for at least few weeks after the stress appearance. Practitioners and urban planners need to take into account that the impact of a stress factor on tree physiology can occur many weeks beyond its occurrence.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/f12060761/s1, Table S1: Locations of investigated trees with number of individuals. Figure S1: Locations of investigated trees presented on Krakow's map, red is for T. × europaea, blue is for T. tomentosa and green is for T. cordata; Table S2: Table with the means of parameters derived from OJIP-test and chlorophyll pigment content, for investigated Tilia species. In the column 'term' FD is for full development, ES is for early senescence and LS is for late senescence. In the column 'year' N is for normal year and D is for dry year (with precipitation deficit in summer). Funding: The presented study was supported by The Ministry of Education and Science, Subsidy SUB2020-050014-D011 and Subsidy SUB/2020-050012-D011.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.