Leaf morpho-physiological traits of Populus sibirica and Ulmus pumila in different irrigation regimes and fertilizer types

Site description

The present study is one of the research outputs of a big project in Mongolia. Hence, the data were collected as previously described in our previous study (Byambadorj et al., 2021b). The study took place in Lun soum, Tuv Province, Mongolia, specifically located at latitude 47°52′15.43″N and longitude 105°10′46.4″E. Lun soum is situated approximately 135 kilometers west of Ulaanbaatar, and its elevation is 1,130 m above sea level. The research site is defined by the presence of Middle Khalkha dry steppe, an area that has undergone significant degradation due to extensive grazing activities by livestock (Ulziykhutag, 1989). The dominant vegetation in the area consists of xerophytic and mesoxerophytic grass species (e.g., Stipa sareptana subsp. krylovii Roshev. and Cleistogenes squarrosa Trin. (Lavrenko, Karamysheva & Nikulina, 1991)). The study area is characterized by Kastanozems soil, specifically of the Loamic type, which is notable for its absence of distinct horizon differentiation (IUSS Working Group WRB, 2015). The surface soil and underlying soil hardness measure 4.5 kg cm⁻² and 1.7 kg cm⁻², respectively (Batkhishig, 2016).

Typically, the growing season in Mongolia commences in either May or June and concludes with the arrival of the initial frost in September (NAMEM, 2021). Throughout the duration of the study, the average temperature recorded was 0.5 °C, while the total rainfall ranged between 93.1 mm and 298.7 mm (2010–2021). Environmental data was gathered at the study site using the HOBO climate data logger. The data revealed that the highest average temperature (21.3 °C) and precipitation (71.8 mm) occurred during the months of July and August. Conversely, the lowest average temperature (10.5 °C) and precipitation (7 mm) were observed in May and September, respectively (Fig. 1).

Experimental materials and design

This study was conducted in a 12-year-old of Populus sibirica and Ulmus pumila during two growing seasons (i.e., July 2021 and 2022). We started the experiment in mid-May 2011 and started to measure leaf morpho-physiological traits from June to September of 2021 and 2022. Because the environment in Mongolia is quite varied, we were able to better identify seasonal trends or changes in the characteristics we analyzed by involving two growth seasons. In 2011, a total of 512 two-year-old seedlings were utilized for tree planting. These seedlings had initial root collar diameter (RCD) and height measurements of 0.51 ± 0.02 mm and 68 ± 2.94 cm, and 0.33 ± 0.01 mm and 51 ± 1.14 cm, respectively. The selected species possessed initial root lengths ranging from 21.50 to 25.9 cm. The planting process involved placing the seedlings into circular holes, which were approximately 60–70 cm deep and 50–60 cm wide. Furthermore, the seedlings received equal amounts of irrigation across all treatments during the acclimatization period, allowing them time to adapt to their new environment. Following the initial month, the study implemented four distinct irrigation regimes: a control group with no irrigation (0 L h⁻¹), a group with irrigation at a rate of 2 L h⁻¹ (0.25 mm m⁻²), a group with irrigation at a rate of 4 L h⁻¹ (0.5 mm m⁻²), and a group with irrigation at a rate of 8 L h⁻¹ (1.0 mm m⁻²). Irrigation was applied every four days for a duration of five hours each time. An irrigation hose system connected to a 25-ton water tank facilitated the watering process.

Additionally, two types of fertilizers were used, namely NPK and compost (COMP), in different quantities. Prior to transplanting, these fertilizers were mixed with the natural soil to fill the holes. The NPK fertilizer was composed of solid granules containing a nitrogen, phosphorus, and potassium mixture with a ratio of 12:16:4. The quantities of fertilizers used were 120 g tree-1 for NPK and 500 g tree-1 for compost.

The compost used in the study was derived from thoroughly decomposed sheep manure, and it possessed the following chemical characteristics: pH of 7.4, organic matter content ranging from 18.0% to 25.0%, nitrogen content between 5.0 g kg⁻¹ and 7.0 g kg⁻¹, and total calcium (Ca), magnesium (Mg), potassium (K), and sodium (Na) contents of 9.29 g kg⁻¹, 7.02 g kg⁻¹, 9.18 g kg⁻¹, and 0.05 g kg⁻¹, respectively.To carry out the experiment, a total of twelve plots, each measuring 20 m × 10 m, were established for each tree species at the study site. These plots were divided as follows: four plots for the irrigation-only treatments, four plots for the irrigation combined with NPK fertilizer treatments, and four plots for the irrigation combined with compost treatments (Fig. 2). The plots were arranged with a spacing of 2.5 m between trees and 2.5 m between rows. For the irrigation-only treatments, a total of 32 seedlings were used, while for the irrigation combined with fertilizer or compost treatments, 16 trees were utilized.

Measurements of leaf morphological traits

A total of 1,200 leaves, representing three replicates for each of the ten treatments and two tree species (20 leaves × three replicates × 10 treatments × two species), were collected for the determination of leaf area (LA) and specific leaf area (SLA). Leaves were not collected from the controlled plots as, by the end of the experiment, all seedlings in those plots were either dead or in a deteriorated condition. The collected leaf samples were initially placed in zip lock bags and stored in a cold container for transportation to the laboratory. Leaves were scanned using an HP LaserJet scanner (M1132 MFP, Palo Alto, CA, USA) at a resolution of 600 dpi. The image analysis software, ImageJ, was employed to analyze the leaf area (LA) based on the procedures outlined in Schneider, Rasband & Eliceiri (2012), with modifications as described in Hernandez et al. (2020). For the determination of specific leaf area (SLA), the leaves were dried in an oven at 65 °C for three days. Subsequently, the dried leaves were weighed using a high precision electronic scale (d = 0.001 g, Discovery Semi-Micro and Analytical Balance, Ohaus Corp., Switzerland). The SLA was calculated using the equation SLA = LA (cm²) divided by the leaf dry mass (mg), following the procedure outlined by Li, Shoup & Chen (2005).

Measurements of leaf physiological traits

Chlorophyll concentration

A random selection of 66 trees (n = 9) was made, and leaves that were healthy, fully expanded, and exposed to sunlight were harvested for the purpose of chlorophyll content analysis. The MC-100 leaf-clip chlorophyll meter from Apogee Instruments, Inc. (Logan, UT, USA) was used to perform the measurements. The chlorophyll content index (CCI), as described by Padilla et al. (2018), was utilized for evaluating the chlorophyll levels.

Stomatal conductance

The measurement of stomatal conductance involved using the SC-1 leaf porometer on sun-exposed leaves from the mid crown of 66 trees. The measurements were conducted between 08:00 and 11:00 h, with a sample size of nine trees. The SC-1 leaf porometer from Decagon Devices, Inc. (Pullman, WA, USA) was chosen for its capability to assess leaf conductance rates while considering the influence of surrounding environmental factors. The calibration of the SC-1 leaf porometer ensured accurate and reliable measurements. To ensure accurate measurements, the calibration of the SC-1 leaf porometer was conducted before conducting field observations. This calibration was necessary as the leaf clip needed to reach temperature equilibrium with the surrounding environment. The readings from the leaf porometer were taken on the same leaf as the infrared thermometry (IRT) temperature readings to maintain consistency in the measurements. The SC-1 leaf porometer was used to measure leaf stomatal conductance automatically, providing readings in mmol m⁻² s⁻¹. The measurements were taken over a period of 30 s, during which the porometer also recorded air temperature and humidity. Stomatal conductance values near 0 mmol m⁻² s⁻¹ indicate significant stress, while values near 500 mmol m⁻² s⁻¹ suggest no stress on the leaves.

Chlorophyll fluorescence

Chlorophyll a fluorescence measurement (n = 9) were performed on leaves collected from a total of 66 trees. The measurements were carried out on a sunny morning, between 08:00 and 11:00 h, prior to the midday depression of photosynthesis. The chlorophyll a fluorescence was quantified using a Handy PEA+ fluorimeter (Plant Efficiency Analyzer, Hansatech Instruments Ltd., Pentney, UK). This instrument allows for the measurement of changes in chlorophyll a fluorescence over a duration ranging from ten microseconds to one second. To prepare the leaves for measurement, a dark adaptation process was employed using specialized leaf clips. Following a dark-adapted period of 30 min, the leaves were exposed to saturating red light pulses with a peak wavelength of 650 nm and intensity of 3,000 µmol (photon) m⁻² s⁻¹. The signal gain for the season was set at × 1.0, as described in the works of Strasser, Srivastava & Tsimilli-Michael (2000) and Yusuf et al. (2010). Table 1 provides a detailed description of the measurements related to the chlorophyll fluorescence OJIP transients.

Leaf water potential

Using the methods outlined in Scholander et al. (1964), the predawn (ψ_p) and midday (ψ_m) leaf water potentials were measured in healthy, fully expanded, and sun-exposed apical leaves. A pressure chamber (Model 1505D EXP, PMS Instrument Company, Albany, OR, USA) was employed for this purpose. The measurements were taken at a height of approximately 1.5 m above the ground.

Statistical analysis

To evaluate the impact of irrigation and fertilization treatments on various morpho-physiological traits, such as chlorophyll concentration, stomatal conductance, chlorophyll fluorescence, and leaf water potential, a one-way to three-way analysis of variance (ANOVA) was conducted for each species. To further analyze and compare the treatments, Duncan’s multiple range test (DMRT) was employed. All statistical analyses were performed using the Statistical Analysis Software (SAS) package developed by SAS Institute Inc. in 2014.

Effects different irrigation regimes and fertilizer types on leaf morphological traits of Populus sibirica and Ulmus pumila

Results of the present study revealed significant variation in leaf area (LA) between control and NPK/COMP-treated seedlings of P. sibirica and U. pumila (Table 2; Table S1). As expected, the control group generally yielded smaller LA compared with NPK/COMP-treated seedlings, suggesting that control seedlings’ leaf growth were negatively affected. The effects of the treatments also differed by growing season. In P. sibirica, the 4–8 L h⁻¹ without fertilizer had larger LA than control (CONT), NPK, and COMP treatments, which all had similar effects on LA, during the first growing season. A different pattern in LA was observed during the second growing season, i.e., 4 L h⁻¹ + COMP resulted in the highest LA compared with the other treatments. In the case of U. pumila, the 2 L h⁻¹ without fertilizer and 2 L h⁻¹ + NPK had larger LA compared with CON and other NPK/COMP treatments for 2021 and 2022 growing seasons, respectively. Moreover, the results of our study have revealed an intriguing interaction effect between year, irrigation, and fertilization on LA for both P. sibirica and U. pumila (Table 2; Table S2). Over time, such a combination resulted in a substantially worse effect in the LA in the control groups lacking irrigation and fertilizer.

The SLA also varied significantly between control and NPK/COMP-treated seedlings of P. sibirica and U. pumila, and the pattern differed by growing season (Table 2; Table S2). We also found a significant egffect of year × irrigation × fertilization interaction on SLA of U. pumila (Table S2). Among treatments, P. sibirica seedlings at 4 L h⁻¹ without fertilizer had the highest SLA for 2021 and at 2 L h⁻¹ + NPK for 2022. A different pattern was observed in U. pumila, i.e., CONT and 4 L h⁻¹ had the highest SLA than the other treatments for 2021 and 2022 growing seasons, respectively.

Effects different irrigation regimes and fertilizer types on leaf physiological traits of Populus sibirica and Ulmus pumila

The chlorophyll concentration, stomatal conductance (g_s), and chlorophyll fluorescence varied significantly between treatments in both P. sibirica and U. pumila (Figs. 3–6; Table 2; Table S3). In this study, we found no significant effects of year × irrigation × fertilization interaction on leaf physiological traits of both species, except in stomatal conductance of U. pumila (Tables S4; S5). However, the year × fertilization interaction significantly affected F_m, F_v, F_v/F_m, and PI_ABS in both species (Table S5). The effects of irrigation × fertilization interaction on the chlorophyll fluorescence parameters, except in F_m, were also significant in the case of U. pumila. In P. sibirica, the chlorophyll concentration in the first year was generally higher in CONT, watering-alone treatments, and 2 L h⁻¹+NPK compared with the other treatments (Fig. 3). This contradicts our hypothesis that the leaf morpho-physiological characteristics of control seedlings would be significantly worse than those of treated seedlings. A different pattern was observed in U. pumila, i.e., the 2 L h⁻¹ without any fertilizer yielded a significantly higher chlorophyll concentration compared with other treatments during the second year of the experiment. Notably, the NPK-alone treatment in U. pumila had the lowest chlorophyll concentration among the treatments.

The stomatal conductance (g_s) varied significantly between treatments, and the pattern was different by year in both P. sibirica and U. pumila (Fig. 4). In P. sibirica, CONT had a significantly greater g_s than the other treatments during the first year. A contrasting pattern was observed during the second year of the experiment, i.e., the highest g_s was found in 2 L h⁻¹+NPK, while the lowest values were generally recorded in CONT, all watering-alone treatments, and in 4 L h ⁻¹+ COMP. In the case of U. pumila, the CONT with NPK or COMP generally had a higher g _scompared with the other treatments in both experimental periods. The lowest values were generally recorded in 8 L h⁻¹+ NPK in both experimental periods, although all water + fertilizer treatment combinations decreased relative to CONT with and without NPK or COMP.

The whole OJIP curves of P. sibirica and U. pumila under different watering and fertilization treatments are shown in Figs. 5 and 6. Here, we found different O-J-I-P kinetic curves patterns for only water treatments, water + NPK treatments, and water + COMP treatments for both species and experimental periods. The general pattern is that the fluorescence transients at CONT were characterized by lower values of curves than those at either NPK/COMP-treated ones. Similar O (origin) phase curve (minimal fluorescence) was observed across treatments and species in both experimental periods. In the P (peak) phase of P. sibirica, the 2 L h⁻¹ generally had the highest fluorescence, while the lowest was observed in CONT for only water treatments, water + NPK treatments, and water + COMP treatments, particularly in 2021 experimental period. Different patterns were observed in U. pumila depending on treatments. For only water treatments, 2 L h⁻¹ generally had the highest fluorescence in both experimental periods. The CONT and/or 2 L h⁻¹ + NPK yielded the highest fluorescence in 2021, but a different pattern was observed in 2022 (i.e., 2 L h⁻¹ + NPK = 4 L h⁻¹ + NPK = 8 L h⁻¹ + NPK > CONT) for water + NPK treatments. A contrasting pattern was observed in the case of water + COMP treatments, i.e., 4 L h⁻¹+ COMP and/or 8 L h⁻¹+ COMP generally had the highest fluorescence among the treatments for U. pumila.

We observed changes in chlorophyll a fluorescence transient parameters across different treatments, and the patterns differed by year depending on the type of fertilizer and amount of irrigation (Figs. S1; S2). In both species, no significant changes in other parameters (e.g., F_v/F_m, F_o, TR_o/RC) were detected during the study periods (Table 3). Here, PI_ABS and F_m were generally lower in CONT with/without NPK or COMP for both P. sibirica and U. pumila (Table 3). In P. sibirica, PI_ABS and F_mwere generally greater in 2 L h⁻¹ + NPK or COMP throughout the first experimental period. Different patterns in the treatment with the highest value were detected throughout the second year for both PI_ABS (i.e., only 8 L h⁻¹, 8 L h⁻¹ + NPK, 2 L h⁻¹ + COMP) and F_m (i.e., only 2 L h⁻¹, 2 L h⁻¹+ NPK, 2/4/8 L h⁻¹+ COMP).

The effects of applied irrigation and fertilizer on diurnal leaf water potential are shown in Figs. 7–8. Among the treatments, predawn and midday leaf water potentials of P. sibirica and U. pumila were generally the lowest in CONT, followed by 2 L h⁻¹+NPK/COMP during the first growing season. For P. sibirica during the first year of the experiment, leaf water potential began to drop from 5:00 to 13:00 (midday), especially for water + COMP treatments, and began to increase as early as 9:00 to 23:00 for only water and water + NPK treatments. During the second growing season for P. sibirica, an obviously different pattern was observed, i.e., the 8 L h⁻¹+COMP had the lowest midday leaf water potential among the treatments. For U. pumila during the second growing season, leaf water potential at only NPK treatment had declined much lower (more negative) than only water or only COMP treatments, although it showed a tendency to increase starting from 19:00.

The pattern of treatment effects on predawn (ψ_p) and midday (ψ_m) leaf water potentials of P. sibirica and U. pumila differed between 2021 and 2022 growing seasons (Table 4). In P. sibirica, CONT generally had a significantly lower (more negative) ψ_p and ψ_m compared with NPK + 2L h⁻¹/4L h⁻¹ irrigation, particularly in the second growing season. The water + NPK treatments showed a tendency to have a higher ψ_p and ψ_m than water + COMP treatments, regardless of the amount of irrigation for P. sibirica. The ψ_p was generally the lowest (more negative) at NPK without water among NPK and COMP treatments for U. pumila. Lastly, the ψ_m of U. pumila was the highest in either NPK or COMP alone treatments than NPK/COMP + water treatments.

Effects of different fertilization and watering regimes on the leaf area and specific leaf area of Populus sibirica and Ulmus pumila

As expected, the leaves of P. sibirica displayed notable responses to varying fertilization and watering conditions, aligning with the findings of previous studies (Larwanou, Adamou & Abasse, 2014; Guo et al., 2021). In the initial growing season, significantly higher leaf area (LA) and specific leaf area (SLA) were observed in the 4 L h⁻¹ and/or 8 L h⁻¹ without fertilizer treatments compared to the control (CONT) and other treatment groups. This suggests that the presence of water alone may have been sufficient for stimulating leaf growth in P. sibirica, without the need for additional NPK or COMP fertilizers. Given the semi-arid nature of the study site, it is possible that P. sibirica seedlings initially prioritized energy investment in expanding leaf area by utilizing water resources more efficiently during the early stages of development. This preference for water over fertilizer may be attributed to the higher water and energy requirements associated with the absorption and uptake of dissolved nutrients via mass flow, resulting in lower leaf area as observed in seedlings treated with water in combination with NPK or COMP, particularly in the first growing season. Furthermore, during the initial growth and development stages, cells typically possess meristematic properties and have higher water demands compared to fertilizer requirements. The process of plant cell expansion involves the concurrent uptake of water into the cells (Cosgrove, 2014). This phenomenon can help explain the significant impact of the interaction between year, irrigation, and fertilization on the leaf area (LA) of P. sibirica seedlings. During the initial stages of growth, there was a greater demand for irrigation compared to fertilizer, which may have intensified as the seedlings continued to develop. The need for fertilizer could have arisen from a gradual increase in growth over time. This finding further elucidates the lower LA observed in the control group (no irrigation), as water likely played a crucial role in stimulating leaf area growth.

However, a different trend was observed in U. pumila seedlings, where a lower water amount (2 L h⁻¹) resulted in higher leaf area (LA) compared to the other treatments during the first growing season. This suggests that U. pumila has a lower water requirement for enhancing leaf area growth compared to P. sibirica, despite both species being known for their drought resistance. The potentially higher drought resistance of U. pumila can be inferred from the significant effect of the interaction between year, irrigation, and fertilization on LA and specific leaf area (SLA) of P. sibirica. Additionally, U. pumila is classified as an invasive noxious tree, capable of reproducing through wind-dispersed seeds and sprouts, and able to thrive in regions where other trees struggle, even in dry mesic areas (USDA, 2017). It is possible that P. sibirica seedlings had already adapted to the prevailing environment after the first growing season, leading to a contrasting pattern in LA (i.e., 4 L h⁻¹ + COMP > CONT) and SLA (i.e., 2 L h⁻¹ + NPK > CONT) during the second growing season. The findings indicate that the combination of water with COMP and/or NPK is necessary to support the growth of the seedlings, especially in the second year after the initial growth stage, in order to achieve their maximum growth potential. This combination of water and fertilizers likely stimulated the activity and movement of biological organisms present in the compost, leading to improved soil aeration, moisture retention, and nutrient availability. Specifically, the application of organic fertilizer enhanced the nitrogen content in the soil, a crucial element that influences the vegetative development of plants by promoting increased leaf fresh weight, improved quality, and larger leaf area (Mojeremane et al., 2015). Furthermore, during the second growing season, the treatment of 2 L h⁻¹ + NPK resulted in a higher leaf area (LA) compared to other treatments in P. sibirica. This outcome can be attributed to the higher drought tolerance of U. pumila, where the addition of COMP did not provide significant benefits in terms of improving soil moisture and nutrient availability. Species with high invasiveness, like U. pumila, possess the ability to not only thrive in low-resource environments but also enhance macronutrient accumulation, soil properties, and soil microbial communities (Souza-alonso & Novoa, 2014; Ye et al., 2019). This explains the observed larger specific leaf area (SLA) of U. pumila in the control treatment and 4 L h⁻¹ without fertilizer during the first and second growing seasons, respectively.

Effects of different fertilization and watering regimes on physiological traits of Populus sibirica and Ulmus pumila

In this study, the leaf physiological traits of both species were significantly influenced by different water and fertilizer treatments, which aligns with the findings of previous experiments (Sumaira, Tao & Zhou, 2017; Byambadorj et al., 2021b; Byambadorj et al., 2021a; Guo et al., 2021). Interestingly, higher chlorophyll concentration was observed in the control treatment, all treatments with water only, and the 2 L h⁻¹ + NPK treatment of P. sibirica seedlings compared to the other treatments, particularly during the first growing season. Similarly, the 2 L h⁻¹ treatment without any fertilizer resulted in a higher chlorophyll concentration compared to the other treatments in U. pumila. Despite the semi-arid conditions at the study site, our results indicate that both species are capable of synthesizing chlorophyll and maintaining high levels in their leaves even without the addition of fertilizer or irrigation. This is noteworthy considering the significantly lower fluorescence parameters (PI_ABS and F_m) observed in the control treatment. While many studies have reported a significant decline or alteration in chlorophyll content in plants subjected to long-term drought stress, most of these studies focused on herbaceous plants (Li et al., 2006; Khayatnezhad & Gholamin, 2012; Chen et al., 2016). Therefore, the observed stability of chlorophyll content in P. sibirica and U. pumila could serve as a potential indicator of their adaptation to drought conditions.

The current study also revealed different patterns in g _s depending on the growing season and species, which can be attributed to the inherent adaptation of the species to the prevailing environmental conditions at the study site. In the case of P. sibirica, g_s was significantly higher in the control treatment (CONT) during the early growth stage. This higher g_s in CONT can be attributed to the smaller leaf area of seedlings under this treatment, leading to increased exposure of leaves to higher solar radiation and air temperature. Previous research has shown that despite a decrease in leaf water potential, stomatal conductance of Populus and Pinus species increases with rising temperatures (Urban et al., 2017). However, prolonged increases in g_s can result in substantial water loss through transpiration. This explains the generally lower (more negative) predawn and midday leaf water potential observed in P. sibirica under CONT. It is likely that the species avoids excessive water loss by reducing leaf area during the early growth stage. Smaller leaves with fewer stomata can minimize water loss through transpiration (Wang et al., 2019; Hernandez & Park, 2022). This also explains the higher g_s observed in the 2 L h⁻¹ + NPK treatment during the second growing season. As compared to the early growth stage in the first season, the demand for water and nutrients may have increased as the seedlings grew rapidly. Higher g_s may have allowed growing leaves to maximize photosynthetic exchange and water conductance before senescence, which can be accelerated by abiotic stresses such as drought, nutrient limitation, extreme temperatures, and oxidative stress commonly experienced in dry land environments.

In the case of U. pumila, the treatment with NPK/COMP fertilizers without irrigation resulted in the highest g_s among the treatments, despite lower ψ_p and ψ_m. This finding supports the results of Cavaleri & Sack (2010), who reported that invasive species exhibit higher stomatal conductance and sap flow rates compared to native species. The application of fertilizer may have played a role in regulating stomatal movement in U. pumila, leading to increased g_s even in the dry environment of the study site. For instance, a study demonstrated that both mesophyll and stomatal conductance were reduced in seedlings with potassium deficiency (Jin et al., 2011). Potassium deficiency can cause dryness throughout the plant and hinder shoot development and water-use efficiency (Arquero, Barranco & Benlloch, 2006). Therefore, the observed high g_s in U. pumila under the NPK/COMP treatment without irrigation may be attributed to a dehydration coping mechanism that helps maintain cell turgidity despite decreasing leaf water potential, thereby avoiding excessive water loss through transpiration.

Implications for effective selection of species potentially suitable in dry lands

When considering the selection of tree species for establishing plantations in semi-arid and arid regions, where climate conditions are projected to become warmer and drier, it is generally preferable to choose species that can conserve water. In this study, both P. sibirica and U. pumila exhibited distinct responses to the irrigation and fertilization treatments, indicating their adaptability to thrive in arid regions with limited water availability. These species have demonstrated characteristics of water conservation, suggesting that they may also be well-suited to nutrient-poor soils and can tolerate low nutrient levels. Specifically, both species displayed the ability to efficiently allocate water and nutrient resources by adjusting leaf area growth, chlorophyll concentration, and stomatal conductance in response to different irrigation and fertilizer levels. This adaptive behavior allows them to optimize photosynthetic exchange and water conductance, even under conditions of limited water and nutrient availability.

The observed higher PI_ABS and F_m in all treated seedlings of both species, regardless of fertilizer type and amount of irrigation, can support such maximization of photosynthetic exchange and water conductance, although those with lower amount of irrigation are much preferable given the arid conditions in the area. Higher PI _ABSsuggests a stable, if not enhanced, overall photosynthetic performance, in treated seedlings of both species. Higher PI_ABS can also be associated with the observed larger LA in treated seedlings (2 L h⁻¹ without fertilizer and 2 L h⁻¹ + NPK) because it also promotes development of photosynthetically active plant organs (Lepeduš et al., 2011). Monitoring PI_ABS of P. sibirica and U. pumila throughout their growing period can provide valuable insights into a plant’s response to drought stress and its adaptive mechanisms.

The response of U. pumila to the COMP fertilizer appeared to be more favorable compared to NPK, regardless of the presence or absence of irrigation. On the other hand, P. sibirica showed a higher degree of responsiveness to NPK when coupled with irrigation. Therefore, the suitability of these species in a specific semi-arid or arid environment may depend on the initial levels of soil water and nutrients. It is important to assess the initial conditions of soil moisture and nutrient content in the target planting area, or similar environments, prior to introducing seedlings of P. sibirica and U. pumila. This step is crucial to ensure long-term sustainability and achieve a high rate of plant survival.