1. Introduction
Among Middle Eastern countries, Iran faces unprecedented climate change. For example, during the past 40 years, the precipitation decrease was 2.56 mm per year in northwest Iran [
1]. An increase in temperature of at least 2 °C is expected in the next decades. This increase in temperature is predicted to be accompanied by a severe decline (~35 %) in precipitation and an anomalous increase in upward longwave radiation, which may be further enhanced by CO
2 emissions that are the largest among all Middle Eastern countries [
2].
Although Iran is mostly an arid or semi-arid country, part of the country with a temperate climate near the Caspian Sea is becoming semi-arid because of the climate change [
3]. Hyperspectral images and vegetation indexes indicate that climate change and human population dynamics have reduced the forest area in Iran by about 43% from 1972 to 2009 [
4,
5]. Reforestation is required to increase landscape water-holding capacity and also to increase biomass accumulation and, therefore, carbon storage [
6].
Many research articles have evaluated the positive effects of hydrogel on selected forest tree species [
7,
8,
9]. Hydrogel is a synthetic polyacrylamide with a large capacity for water retention and storage. Application of hydrogel can counter conditions of low precipitation, high evapotranspiration demand, and low soil water-holding capacity.
Trees and other plants grown in hydrogel-treated soil showed improved physiological and morphological traits and increased survival, water-use efficiency, and dry matter production [
10,
11]. By forming a flexible envelope, hydrogel mimics the effects of mucilage naturally exuded by roots in order to maintain ion and water exchange processes between the root system and the rhizosphere. Root exudates containing amino acids and saccharides positively influence the soil microorganisms and therefore increase the availability of organic and inorganic nutrients to roots [
12]. The leaves of hydrogel-treated beech and eucalyptus trees contained higher contents of nutrients (N, K, and Mg) than the leaves of untreated trees [
13,
14]. Hydrogel-treated plants generally require 20% less fertilizer than untreated plants [
13]). By increasing the uptake of Ca
+, Na
+, and Cl
−, hydrogel may help plants survive the combination of drought and salt stress [
15].
Positive effects of hydrogel have been reported by many studies, but several studies have reported non-significant effects of hydrogel on plant growth [
13]. Germination seems to be the only trait negatively influenced by hydrogel application [
8,
12]. Hydrogel might suppress germination because the polymer is able to withhold large amounts of water when excessive irrigation is applied [
10]. When used alone, hydrogel may also harm seedlings by reducing root aeration [
12,
16]. Therefore, the application of an organic substrate (vermiculite, bark, etc.) along with hydrogel is highly recommended to increase soil porosity and thus facilitate root and rhizosphere gas exchange. For the same reason, hydrogel may also be less effective in clay or loamy soils than in sandy soils [
17].
Although many studies have documented positive effects of hydrogel on plant growth, a comprehensive study with different mixture treatments with hydrogel content and with several tree species in a semi-arid ecosystem is lacking. The goal of the current research was to conduct such a study with a focus on the morphological and physiological traits of the studied species. In a replicated field experiment, we compared the effects of hydrogel treatments with a non-treated control on the performance of 20 common tree species that differ in their drought resistance.
2. Materials and Methods
The research was conducted at 15 km east of Tehran in 2008 (35.819 N, 51.511 E; 1857 m a.s.l.). The average slope was 8%, and slopes were oriented to the southeast. The average annual precipitation was 208 mm, and most precipitation occurred in the winter. Mean annual rainfall ranged from 51 to 1835 mm [
18]. The soil was a sandy loam with a pH of 7.5 and an extremely low N content (0.02% on a mass base).
The experiment used 2-year-old saplings that had been grown in a plantation in Karaj, which is 50 km from Tehran. The experiment occupied 2500 m2 and was protected against livestock and rodents by a 1.5-m-high metal mesh fence. During the spring of 2008, the saplings were planted 2.0 m from each other at the study site. The following seven treatments were applied immediately after planting:
- (1)
CO (control)
- (2)
OM (hydrogel and sawdust)
- (3)
HM (hydrogel and organic fertilizer)
- (4)
HC (hydrogel and compost)
- (5)
HK (hydrogel and wheat straw)
- (6)
HS (hydrogel and subsoil)
- (7)
HP (hydrogel and subsoil with a cobbles cover)
The seven treatments were applied in a factorial design to 20 tree species (
Table 1) with 10 replicate saplings for each combination of treatment and tree species, giving a total of 1400 plants (7 treatments × 20 species × 10 replicates).
A 1250-g quantity of the hygroscopic polymer (trade name STOCKOSORB® 300; Evonik Nutrition & Care GmbH, Essen, Germany) was added to 100 L of water. After it congealed, 2 L of the mixture (plus added materials) was applied to the soil in the root zone of each sapling. The amount of added material (wheat straw, compost, subsoil, sawdust) was approx. 1/3 of the total volume that was added to the roots. The fertilizer in the form of capsules (10 g, with N, P, K nutrients) was administered to the roots together with the hydrogel. The control plots received no fertilizer.
The following parameters were evaluated 6 months after the saplings were planted and treated: height, health (on a scale from 1 to 4, where 1 was a dead plant and 4 was a vigorously growing plant), and leaf osmotic potential at 4 pm (measured during 1 week in 2008). In total, six leaves per treatment and tree species were measured. The vapor-pressure osmometer Wescor Model 5500 (Wescor, UT, USA) was used for measuring osmotic potential. To estimate drought vulnerability, previous reports of P
50 values for each species were used. P
50 is the leaf water potential at which 50% of hydraulic conductivity is lost (
Table 1). Survival was assessed 1 year after planting.
Biochemical analyses of leaf material were performed 6 months after planting. Biochemical parameters were analyzed in control and HS (hydrogel with subsoil) treatments only. The methodology is briefly described in the following paragraphs.
2.1. Photosynthetic Pigments
Photosynthetic pigments were extracted from an assimilation apparatus in acetone. The extract was passed through filter paper and centrifuged. The pigment content was calculated using the following equations [
24]:
Concentration of chlorophyll a [c(chl a)] and b [c(chl b)] and carotenoids c(car) is calculated from the pigment absorbance at different wavelengths (A663, 646, 470) measured spectrophotometrically (DR 3900, Hach Company, Loveland, CO, USA).
2.2. Rubisco
The activity of ribulose-l,5-bisphosphate carboxylase oxygenase (Rubisco) was measured with a spectrophotometer according to Barta et al. 2011 [
25].
2.3. Nitrogen
The N content in dried (48 h, 80 °C) and weighed leaves or needles was measured with the automatic elemental analyzer (CNS-2000, LECO Corp., St. Joseph, MI, USA). Commercial standards (sulfamethazine and alfalfa) were used for calibration.
2.4. Electrical Conductivity
The electrical conductivity of stem samples was measured with a Milwaukee Mi180 multimeter after the samples were cooled to 5 °C, warmed to 20 °C for 1 h, and then autoclaved. Electrical conductivity was measured as the ratio of values before and after autoclaving.
2.5. Prolin
Prolin in an assimilation apparatus was estimated with the spectrophotometer (DR 3900, HACH Company, Loveland, MI, USA) based on the reaction between amino acid ninhydrin (2,2-dihydroxyindane-1,3-dione) [
26].
2.6. Abscisic Acid
ABA was measured by GS-MS according to Okamoto et al. [
27].
2.7. Total Protein
Protein concentration in assimilation was determined according to Bradford using bovine serum albumin as the standard [
28].
2.8. Zinc
Zn content in assimilation was measured with an atomic absorption spectrophotometer (AA-670, Shimadzu, Columbia, MD, USA) according to the methods of Diatloff and Rengel (2001) [
29].
2.9. Relative Water Content
For determination of RWC, the weight of leaves or needles was measured in their fresh state (FW), after immersion for 3 h in de-ionized water (TW), and after drying in an oven (DW). RWC was calculated according to the following equation:
In order to make the statistical results clearer, the findings were expressed for two treatments only: control one and all treatments with hydrogel addition were embodied into one hydrogel treatment (although no pure hydrogel treatment was used). All treatments with hydrogel addition (e.g., OM, HM, HC, HK, HS, HP) were joined together into one group. Finally, two groups were compared—control and hydrogel treatments. Nevertheless, the detailed results of statistical differences among all tree species and all seven treatments are part of
Appendix A.
The selected trees included 15 broadleaf and 5 coniferous species that are common in Western Asia (
Table 1). For better visibility of key findings, we divided tree species into two groups. Drought-tolerant (TOL) and drought sensitive (SEN) tree species were selected on the basis of P
50 and differences in biochemical parameters. TOL species comprise
Olea,
Eleagnus, and both
Pinus species, the rest of the investigated tree species was part of the SEN group.
Statistical analysis was conducted using R software environment, version 3.6.1 (R Core Team, 2019). Plots were made using packages ggplot2 [
30] and fmsb [
31].
Prior to statistical analyses, data were checked for normality and log-transformed when appropriate.
To assess differences among treatments and species combinations, a generalized least squares model was applied (gls function from nmle package) [
32], which allowed to model specific variance structures. Possible variance structures were fitted, and candidate models were compared by Akaike’s Information Criterion (AIC).
The model with the lowest AIC value was used. The normality of the distribution of the residuals was assessed using QQ-plots.
The model can be written in a form:
where log(height)
ijkl is log-transformed height increase of ith observation, for jth species, kth treatment and lth site; species is categorical variable with two levels (sensitive and tolerant) and treatment is also categorical variable with two levels (control and hydrogel); ε
ij as an error term of ith observation at lth site. Variance of errors was specified as varIdent(~1|site), allowing estimation of different variances at each of 10 sites and varExp(form = ~height) combined together with varComp function. Residuals are normally distributed with mean 0 and var(ε
il) =
×
:
where δ is the estimated parameter.
Post-hoc pairwise comparisons were performed by the Tukey HSD method implemented in emmeans function from emmeans package. As we worked with log-transformed data, reported values represent back transformed means.
For survival data, prop.test function from base R was used to perform Pearson’s chi-squared test statistic to evaluate the significance of differences among hydrogel treatments. Subsequently, post-hoc pairwise comparison was performed using
p-values adjusted by Holm’s method for correction for multiple testing as implemented in pairwise.prop.test function from base R. Detailed information is part of
Table A1,
Table A2 and
Table A3.
3. Results
Relative to the control, each of the treatments with hydrogel significantly increased the survival rate and positively influenced the height increases and health status of all tree species (
Figure 1 and
Figure 2).
Both drought-tolerant (TOL) and drought-sensitive (SEN) were in very poor condition (no excellent health status) without hydrogel. Survival of all species was 37% in control and 81% in all hydrogel treatments together. The application of hydrogel together with subsoil and cobbles on the surface was especially beneficial for
Eleagnus,
Robinia, and
Juniperus species. The increase in height for the treatment with cobbles was 35%–45% greater than for the other hydrogel treatments (
Table A4). Health status was highest for
Pinus nigra,
Olea, and
Erythrostemon species when they were treated with hydrogel and organic fertilizer or compost (
Pinus and
Olea), or with hydrogel and sawdust or subsoil (
Olea and
Erythrostemon). For
Erythrostemon, a health status of four (i.e., excellent) was recorded for the hydrogel treatment with sawdust, and 70% of all trees were considered to be in perfect shape. For the untreated control, no
Erythrostemon tree had a health status of four, and 70% had a health status less than three. No matter what kind of hydrogel treatment was used
Olea,
Cercis,
Celtis, and
Erythrostemon species showed the highest rate of survival compared to other tree species with any hydrogel treatment. On the other hand, survival of
Platanus,
Acer, and
Juglans was similar among untreated control trees and those treated with hydrogel (
Table A4).
The key biochemical traits were influenced in the HS treatment (hydrogel with subsoil) compared to control one (
Figure 3 and
Table A5). Nevertheless, no statistically significant differences in biochemical traits between control and hydrogel treatments were found.
4. Discussion
Hydrogel application has been reported to enhance the aboveground growth of the plant, support survival, and significantly improve the health status [
10,
11,
14,
16]. Moreover, our findings indicate that the magnitude of the response to hydrogel strongly depends on the drought-sensitivity or drought- tolerance of the tree species. In response to hydrogel treatment, the increases in leaf relative water content (RWC) and leaf dry mass (LDM) were expected. Both of these traits were previously recognized as common effects of hydrogel treatment [
11]. Similar to the results of M’barki [
33] hydrogel increased growth, dry weight, RWC because of stronger uptake of the micro- and macronutrients, especially nitrogen and potassium [
34]. With a higher accumulation of nutrients, the increased amounts of total proteins, photosynthetic pigments and Rubisco could be explained [
35]. Similarly, a higher quantum yield (Fv/Fm) or higher photosynthetic rate [
8] is usually reported for hydrogel-treated plants. One of the main benefits of hydrogel treatment is an increase in prolin. In both the SEN and TOL group, higher prolin concentration was recorded. From this point of view, no correlation with osmotic potential or RWC is surprising because the osmolytic properties of prolin usually explain higher water potential in hydrogel-treated plants [
36]. In addition to the osmoprotective function, it acts as a metal chelator or antioxidant [
37]. Moreover, the positive effect of hydrogel treatment was proved by lower concentration of malonaldehyde (MDA) or higher activity of superoxide dismutase (SOD) [
38]. Higher concentration of MDA stands for membrane destruction, and the presence of SOD suggests the presence of reactive oxygen species—both signs of drought stress that hydrogel treatment can alleviate.
Our findings imply that hydrogel efficacy depends on tree species susceptibility to environmental stress. Based on our results, the hydrogel application is more beneficial for drought-sensitive tree species compared to drought-tolerant ones. Drought-sensitive trees typically inhabit deep soils that contain sufficient organic matter and with sufficient supplies of water and nutrients. Their P
50 values do not exceed −2.0 MPa (see
Table 1).
Cupressus, in particular, is considered as one of the most drought-tolerant tree species [
19]. Moreover, hydrogel should be applied with additional substances as straw or subsoil that can increase soil porosity and, therefore, oxygen supply to the roots. In this respect, an expanded form of vermiculite can be effective [
16]. Added substrate porosity seems to be crucial for any hydrogel amendment treatments as poor physiological tree growth may be related to the reduced root aeration [
12]. The different mixtures of various substances with hydrogel probably explains the differences among hydrogel treatments in their effects on morphological traits and why hydrogel with subsoil was often the best treatment for any particular morphological trait. The treatment with hydrogel, subsoil, and cobbles was also favorable for height increase. This confirms that the traditional practice of cobbles placement in the Anasazi Rio Grande landscape in the USA improves local water management by trapping seasonal runoff and reducing evaporation [
39].