Water availability and nutrient-status of soils play crucial roles in seedling establishment and plant survival in coal-spoiled areas worldwide. Restoration of spoils pertains to the application of proper doses of nutrients and water, and selection of particular plant species for efficient revegetation. This study aimed at examining the potential effects of different combinations of soil-water and fertilizers (nitrogen, N and phosphorus, P) on morpho-physiological and biochemical attributes of Amorpha fruticosa grown in coal-mined spoils. Three factors five-level central-composite-design with optimization technique response surface methodology (rsm) was used to optimize water irrigation and fertilizer application strategies.
Our results revealed a strong correlation between experimental data and predicted values developed from the rsm model. The best responses of A. fruticosa in terms of plant height, stem diameter, root length, and dry biomass were observed under a high-water regime. Low-water regime caused a notable reduction in growth-associated parameters, and fertilization with either N or P did not show positive effects on those parameters, indicating that soil-water was the most influential factor for growth performance. Leaf water potential, gas-exchange parameters, and chlorophyll content significantly increased under high levels of soil-water, N and P, suggesting a synergistic effect of these factors for the improvement of photosynthesis-related parameters. At low soil-water contents and N-P fertilizer application levels, enhanced accumulation of malondialdehyde and proline indicated that A. fruticosa suffered from oxidative and osmotic stresses. Amorpha fruticosa also responded to oxidative stress by accelerating the activities of superoxide dismutase, catalase, and peroxidase. The effects of both fertilizers relied on soil-water, and fertilization was most effective under well-watered conditions. The maximum growth of A. fruticosa was observed under the combination of soil-water, N-dose and P-dose at 76% field capacity, 52.0 mg kg− 1 and 49.0 mg kg− 1, respectively.
Our results demonstrate that rsm effectively designed appropriate doses of water and N-P fertilizer to restore coal-spoiled soils. Furthermore, A. fruticosa responded to low-water and fertilizer-shortage by upregulating defensive mechanism to avoid damage induced by such deficiencies. Finally, our findings provide effective strategies for revegetation of coal-contaminated spoils with A. fruticosa using appropriate doses of water and N-P fertilizers.
Water-shortage is considered as one of the major ecological limiting factors, affecting the restoration of vegetation and plant productivity. Simultaneously, the soil-water deficiency is usually connected with the unobtainability of nutrients in arid and semi-arid areas of the world . Drought along with nutrient deficiency limits vegetation establishment, and the relationships between vegetation restoration and soil-moisture availability have recently gathered significant attention because of their importance in re-establishment of an ecosystem in the drought-prone areas [2,3,4]. The northwest-part of China considerably covers a big region of the world’s dry-lands . Seedling establishment and plant growth in this region largely depend on natural precipitation . Moreover, this region comprises of many coal-mine areas, contributing around 70% of China’s coal production. Coal-mining leads to serious ecological and environmental problems, such as removal of the earth-surface cover, damage of plant roots, destruction of vegetation and turning of coal-mining areas into deserts [6, 7]. Thus, the coal-mined spoils of northwest-part of China could be an appropriate resource to explore the possibility of ecosystem restoration using particular plant species and appropriate doses of water and fertilizers.
Artificial ecological engineering (revegetation) is the appropriate method for restoration of ecology in an abandoned coal gob piles in order to improve the environmental conditions. However, restoration of such kind of degraded areas is very challenging and considered as a daunting task for the environmentalists . An adequate amount of soil-moisture and nutrients is the key to support plant growth and development in dry-land areas . Moreover, the development of a self-sustaining revegetation ecosystem also depends on the selection of suitable plant species, which, together with effective treatment strategies using soil-water and fertilizers can help increase restoration effectiveness . The plants which are originated in the desert-environment have better adaptable photosynthetic systems compared to plants that originated in other normal environments. To continue life cycles in the desert areas, plants develop several survival mechanisms associated with morphological, physiological and biochemical adaptations [11, 12]. Amorpha fruticosa is a deciduous shrub (1–6 m tall) and able to grow in diverse environmental conditions . It is a commonly used shrub for revegetation in Loess Plateau of China because of its drought resistance and adaptation capabilities under infertile conditions. A. fruticosa is frequently cultivated in areas after landslides and escarpments, and for stabilization of the soils on railway embankments, as it’s roots play important roles in the reduction of soil erosion. It can also effectively dwells in degraded-environments in the form of brushes with its seed and suckers .
Ecological problems due to coal mining and the vegetation reestablishment in coal-mining areas are very important environmental issues, urging in-depth research for the sustainable environment [8, 15]. Numerous findings reported the impacts of soil-water and nutrients on the growth and development of various crops in coal-mined areas [16,17,18]. The results suggested that the coupling effects of soil-water and fertilizers played vital roles on plant growth for the restoration of vegetation. To our knowledge, the combinatorial effects of soil-water and fertilizers on desert shrubs in coal-mined spoils have rarely been performed. Hence, it is worth investigating the growth response of A. fruticosa planted in coal-mined spoils concerning different strategies of water and fertilizer applications. Thus, the objective of this study was to evaluate the growth and development of A. fruticosa in response to different combination of water (W) and fertilizers (nitrogen, N and phosphorus, P) using central composite design (CCD) with optimization technique response surface methodology (rsm), and to find out optimum water and fertilizer doses for vegetation restoration in a coal-contaminated desert area in China. For this purpose, we have examined the effect of W, N and P on the morphological, physiological and biochemical responses of A. fruticosa by considering plant height, stem diameter, root length, dry biomass, root/shoot (R/S) ratio of biomass, leaf water potential (LWP), photosynthesis rate (Pn), transpiration rate (Tr), stomatal conductance (Gs), water use efficiency (WUE), the levels of chlorophyll (Chl), malondialdehyde (MDA), and proline (Pro), and the activities of antioxidant enzymes superoxide dismutase (SOD), catalase (CAT) and peroxidase (POD).
Effects on morphological growth responses of A. fruticosa
The maximum increase of plant height (101.23 cm), stem diameter (7.26 mm), root length (73.12 cm) and dry biomass (42.37 g) were recorded under the treatments of W (80% Field capacity (FC)), N (60 mg kg− 1) and P (90 mg kg− 1) (Table 1). The minimum root-shoot (R/S) biomass ratio of 1.51 was observed when the plants were treated with the highest dose of W (80% FC), and medium doses of N (60 mg kg− 1) and P (90 mg kg− 1) fertilizers (Table 1). The lowest increment of plant height (52.04 cm) was recorded under the combination of W, N and P at 48.1% FC, 95.68 mg kg− 1 and 36.49 mg kg− 1, respectively (Table 1). The increase of stem diameter was 4.83 mm when the plants exposed to medium W (60% FC) and N (60 mg kg− 1), and highest P doses (180 mg kg− 1) (Table 1). On the other hand, the lowest root length (44.96 cm) was recorded when the plants were treated with medium W (60% FC), N (60 mg kg− 1) and P (90 mg kg− 1) doses (Table 1). The minimum dry biomass (24.0 g) and highest R/S biomass ratio (2.09) were observed at W (40% FC), N (60 mg kg− 1) and P (90 mg kg− 1) doses (Table 1). Positive values of the parameter coefficients W (a1) and N (a2) exhibited positive effect of W, and N on plant height, stem diameter, and root length, and the relative magnitude was: W > N (Table 3). In contrast, P (a3) had no positive effect on the above-mentioned parameters (Table 3). For plant dry biomass parameter coefficients, W, N, and P had a positive effect with W having the greatest effect, followed by N and P (Table 3). The coefficient P was positive for R/S biomass ratio in contrast W and N had a negative effect (Table 3).
Effects on leaf water potential, gas-exchange parameters, and chlorophyll contents of A. fruticosa
The maximum increase and decrease of leaf water potential (LWP) (− 1.08 MPa and − 4.81 MPa, respectively) were obtained at medium N (60 mg kg− 1) and P (90 mg kg− 1), but at highest and lowest W (80 and 40% FC) level, respectively (Table 1). The maximum values obtained for photosynthesis rate (Pn) and transpiration rate (Tr) were 13.52 (μmol m− 2 s− 1) and 5.90 (mmol m− 2 s− 1), respectively, under the treatments of W at 71.9% FC, N at 95.68 mg kg− 1 and P at 143.51 mg kg− 1 (Table 1). In contrast, the lowest value of Pn (10.78 μmol m− 2 s− 1) and (Gs) (0.17 mol m− 2 s− 1) were recorded under the combinations of W, N, and P (48.1% FC, 24.32 mg kg− 1 and 143.51 mg kg− 1, respectively) (Table 1). The treatment with W (48.1% FC), N (95.68 mg kg− 1) and P (143.51 mg kg− 1) resulted in minimum Tr (3.64 mmol m− 2 s− 1) and maximum water use efficiency (WUE) (3.48 μmol mmol− 1) (Table 1). All coefficients for LWP, Pn, and Tr were positive (Table 3). Although W and N positively influenced Gs, WUE was effectively increased by N only. For Gs and WUE, P did not show any positive effect (Table 3).
The minimum levels of chlorophyll (Chl) a (2.16 mg g− 1 fresh weight (FW)), Chl b (1.06 mg g− 1 FW) and total Chls (3.22 mg g− 1 FW) were observed in response to the treatments with W at 60% FC, N at 0 mg kg− 1 and P at 90 mg kg− 1 (Table 2). The highest levels of Chl a (2.69 mg g− 1 FW) and total Chl (4.69 mg g− 1 FW) were obtained in response to the treatments with W at 80% FC, N at 60 mg kg− 1 and P at 90 mg kg− 1 (Table 2). Chl a, Chl b and total Chl contents in A. fruticosa leaves significantly increased by the addition of W, N, and P, where N had the maximum positive impact followed by W and P (Table 3).
Effects on malondialdehyde and proline accumulation, and the activities of antioxidant enzymes of A. fruticosa
The malondialdehyde (MDA) content of A. fruticosa seedlings was suppressed by all 3 factors; however, proline (Pro) content was considerably decreased by W and P (Table 3). Among the factors, W was most influential for decreasing MDA and Pro contents in the leaves of A. fruticosa (Table 3). The minimum MDA content of 23.34 μmol g− 1 FW was recorded with W at 71.9% FC, N at 95.68 mg kg− 1 and P at 36.49 mg kg− 1 (Table 2). In contrast, the maximum MDA content was observed under the combinations of W, N, and P (48.1% FC, 95.68 mg kg− 1 and 143.51 mg kg− 1, respectively) (Table 2). The highest W (80% FC) and medium N (60 mg kg− 1) and P (90 mg kg− 1) doses represented minimum accumulation of Pro (1.54 μmol g− 1 FW) in the leaves of experimental seedlings; while the maximum Pro content of 5.15 μmol g− 1 FW was obtained with W at 48.1% FC, N at 95.68 mg kg− 1 and P at 36.49 mg kg− 1 (Table 2).
The superoxide dismutase (SOD) (154.68 U g− 1 FW) and catalase (CAT) (196.15 U min− 1 g− 1 FW) activities considerably increased under the treatment with W at 40% FC, N at 60 mg kg− 1 and P at 90 mg kg− 1. The minimum SOD activity (89.22 U g− 1 FW) was observed at medium W (60% FC), N (60 mg kg− 1), and P (90 mg kg− 1) levels, whereas the minimum CAT activity (152.97 U min− 1 g− 1 FW) was recoreded with W at 71.9% FC, N at 24.32 mg kg− 1 and P at 143.51 mg kg− 1. The lowest peroxidase (POD) activity (132.75 U min− 1 g− 1 FW) in the leaves of A. fruticosa was observed with the combinations of W, N, and P (71.9% FC, 95.68 mg kg− 1 and 36.49 mg kg− 1, respectively). However, the POD activity significantly increased (254.66 U min− 1 g− 1 FW) in response to the treatments with W at 48.1% FC, N at 95.68 mg kg− 1 and P at 143.51 mg kg− 1 (Table 2). The activities of SOD and POD were decreased by W and N; while CAT activity was considerably decreased by W and P (Table 3).
Interactive effect of process variables on morphological parameters
The mutual effect of W, N, and P over plant height and stem diameter displayed that plant height and stem diameter increased linearly with W rising from 40 to 80% FC. At low W dose (40% FC), increment of N and P doses showed a little positive effect on the elevation of plant height and stem diameter, whereas at high W dose (80% FC), plant height and stem diameter significantly increased with the addition of N and P fertilizers (Fig. 1a, b, Additional file 1: Figure S1). The root length initially decreased with the decrease of W dose (from 80 to 60% FC), but after reaching a certain point (at 60% FC), it started increasing with the further reduction of W dose (Fig. 1c, Additional file 1: Figure S1). Under low W dose (40% FC), N fertilizer did not express any significant effects, while at high W dose (80% FC) addition of N fertilizer increased root length slowly (Fig. 1c, Additional file 1: Figure S1). Dry biomass of A. fruticosa seedlings significantly influenced in response to W dose regardless of the doses of N and P fertilizers. The positive influence of N on plant dry biomass was detected at a high W dose (80% FC) (Fig. 1d). In contrast, under low (40% FC) and high (80% FC) W doses, the addition of P-fertilizer slightly increased dry biomass, but continuous P-addition resulted in negative impacts (Fig. 1e). The R/S biomass ratio dramatically enhanced with the decline of W level, whereas P had no obvious effects on R/S biomass ratio (Fig. 1f, Additional file 1: Figure S1).
Interactive effect of process variables on LWP, gas-exchange parameters, and chlorophyll contents
It was observed that LWP significantly declined at low W condition (around 40 to 48.1% FC), and the addition of N and P fertilizers did not show any positive effects on the recovery of LWP (Fig. 2a, Additional file 2: Figure S2). A mutual effect of W with N-P fertilizers over Pn and Tr were observed, and the data showed that both Pn and Tr increased with the addition of W, N and P, and the maximum Pn and Tr were obtained at W (71.9% FC), N (95.68 mg kg− 1) and P (143.51 mg kg− 1) doses (Fig. 2b-d, Additional file 2: Figure S2). The Gs initially increased with W and N doses up to a certain level (W at 60% FC and N at 60 mg kg− 1), whereas it started to decline with the increase of W and N above this level (Fig. 2e). Under low W condition (48.1% FC), WUE significantly increased with the increase of N (95.68 mg kg− 1) and P doses (143.51 mg kg− 1) (Fig. 1f, Additional file 2: Figure S2). The combined effect of W and N, and W and P on Chl a, Chl b, and total Chls indicated that under both low and high W levels, total Chl contents dramatically increased with the addition of N, whereas the interactive effect of W and P had no significant influence on the levels of total Chl (Fig. 2g-h, Additional file 2: Figure S2).
Interactive effect of process variables on MDA and pro contents, and antioxidant enzyme activities
Under low W condition (48.1% FC), the addition of N (95.68 mg kg− 1) and P (143.51 mg kg− 1) fertilizer dramatically increased MDA content in our experimental seedlings (Fig. 3a-b, Additional file 3: Figure S3). The highest W level (80% FC) and medium N-P fertilizer (60 mg kg− 1 and 90 mg kg− 1, respectively) content was beneficial for the decrease of Pro content (Fig. 3c, Additional file 3: Figure S3).
The medium dose of W (60% FC) and N-P fertilizer (60 mg kg− 1 and 90 mg kg− 1, respectively) was suitable for the decrease of SOD activity (Fig. 3d, Additional file 3: Figure S3). When W level was 71.90% FC, CAT activity significantly decreased with the increase of P fertilizer (Fig. 3e, Additional file 3: Figure S3). Under low W and N doses (48.1% FC and 24.32 mg kg− 1, respectively), POD activity considerably increased. On the other hand, POD activity exhibited a decreasing trend in response to the increase of either W or N-P fertilizers (Fig. 3f, Additional file 3: Figure S3).
Optimization of soil-water, nitrogen and phosphorus fertilizer
According to CCD result, the optimum water and fertilizer rate to obtain maximum growth performance of A. fruticosa was determined by Derringer’s desired function approach as follows: W of 75.68% FC, N of 51.64 mg kg− 1, P of 49.49 mg kg− 1 (Fig. 4). Under this condition, maximum plant growth was obtained (Fig. 4). Accurateness of the optimized condition was identified by comparing the average experimental values and the predicted values developed by the model. The desirability ramp obtained from optimal points is represented in Fig. 4.
Our study demonstrated that maximum plant height, stem diameter, root length, and dry biomass were obtained at high soil-water content; and in this situation, the addition of N-P fertilizers significantly improved these growth associated attributes (Fig. 1a-e, Additional file 1: Figure S1). Low water levels caused a significant decline in plant dry biomass, and the effect of N-P fertilizers also depended on the availability of soil-water. At low water-level, increases in the doses of N-P fertilizers exhibited a little effect on the improvement of plant biomass, and the highest biomass was obtained at high water level (Fig. 1d-e; Additional file 1: Figure S1; Table 1). These findings indicated that due to scarcity of soil-water, nutrients present in the soil did not mobilize well, and plant roots could not uptake it adequately, which ultimately prevented water and nutrients flow through the xylem to the surrounding cells, as observed in other plant species under water-shortage conditions [19,20,21]. Thus, the deficiency of nutrients and water eventually led to the inhibition of cell growth and elongation , as was also manifested in this study in terms of plant height and stem diameter (Table 1). However, the mutual influence of soil-water and N-P fertilizers considerably improved the functions of each other, which ultimately positively increased the plant height, stem diameter, root length and dry biomass of A. fruticosa, which were in parallel with previous studies [22, 23]. The previous study reported that soil-water was the most important factor for plant morphological responses, and the effect of fertilizers like N and P on plant growth depends on soil-water content . In this experiment, we observed that R/S biomass ratio significantly increased under low water regime (Fig. 1f; Additional file 1: Figure S1; Table 1), indicating that under water and nutrient shortage conditions, carbon translocated from leaves to roots, thereby increasing R/S biomass ratio . The increased R/S biomass ratio during stress conditions could also be explained as one of the adaptive mechanisms of plants to make greater root surfaces for absorbing much water and nutrients from soils . However, other study reported that a greater R/S biomass ratio under stress conditions was resulted due to decreased shoot growth .
LWP acts as a powerful indicator of plants’ water status, and a reduced level of water has been known to be involved in declining LWP in plant species . Similarly, we observed that low level of water caused a severe reduction of LWP in A. fruticosa (Fig. 2a; Additional file 2: Figure S2; Table 1). Gas-exchange parameters are considered as important physiological parameters, having significant roles in plant growth and development. The status of gas-exchange parameters also correlated with the adaptive responses of the plants to cope under environmental assaults . We also observed that interaction of high level of water with high-level fertilizers significantly increased Pn, Tr, and Gs, whereas water and fertilizer limitations led to a considerable decline of those gas-exchange parameters (Fig. 2b-f; Additional file 2: Figure S2; Table 1). Under adverse environmental conditions, stomata of leaves are closed and Pnis inhibited, and thus, declined levels of Gs and Pn have been considered as indications of stomatal limitation . Furthermore, a noticeable decrease in Gs during stress conditions indicate an efficient adaptive mechanism to control Tr . The severe water-shortage condition is extremely unfavourable for absorbing and utilizing nutritional components from the soils by the plant root-systems . Under well-watered conditions, a reasonable coupling model of water and fertilizers could be adopted to utilize N and P . Our results corroborated with the findings of several reports, and indicated that high soil-water content and fertilizers positively increased Pn, Gs and Tr; however, plants are used to reduce these parameters in order to survive under water-shortage conditions [28, 33]. Our experimental results exhibited that during water-shortage conditions, WUE increased with the addition of N-P fertilizer, and the minimum WUE was detected under high soil-water and N-P fertilizer doses, which were in parallel with a previous report .
The status of Chl content is considered as the tolerance index of the plants under stressful conditions . Total Chl (Chl a + Chl b) content considerably decreased with the increasing level of water-shortage (Fig. 2g-h; Additional file 2: Figure S2; Table 2), indicating that A. fruticosa mounted an adaptive mechanism to protect themselves from the damage of low-water induced photo-inhibition . Gallé and Feller  also reported that continuous drought led to decrease Chl content, reduce stomatal conductance, and finally reduce Pn in the leaves of Fagus sylvatica. However, together with N, soil-water significantly increased Chl content, whereas P did not show any positive effect on Chl content. Since N is an integral part of Chl molecules, N-shortage conditions were associated with the reduced level of Chl and chlorosis-related symptoms in plants [38, 39].
We observed a considerable increase of lipid peroxidation product MDA and osmoprotectant Pro under water-shortage condition compared to the well-watered conditions (Fig. 3a-c; Additional file 3: Figure S3; Table 2). These results suggest that water deficiency led to the production of ROS, which ultimately caused lipid peroxidation of the membrane, leading to greater production of MDA. However, MDA content decreased gradually with the addition of soil-water and N-P fertilizer, and the results showed that high W and N-P fertilizer doses contributed to drought tolerance by inhibiting the damage of cell membranes and enhancing osmoregulation in A. fruticosa seedlings. Saneoka et al.  also reported that MDA content in bentgrass (Agrostis palustris) was higher under drought conditions than those found in control condition, and MDA content decreased with the addition of N to the growth medium. Similarly, a water-deficient state resulted in osmotic stress, and A. fruticosa accumulated a high level of Pro to overcome the adverse effects of osmotic stress. Our results corroborated the findings of others in mulberry (Morus alba) and moth beans (Vigna aconitifolia) [41, 42], and suggested that water-stressed plants increased osmoprotective adaptation by accumulating high levels of osmoprotectants like Pro.
To protect from the deleterious effects of ROS produced under stressful conditions, plants normally increase their protective mechanisms against ROS by enhancing the activities of several enzymes, including SOD, CAT and POD . A wealth of reports suggest that exposure of plants to drought led to the enhancement of the activities of antioxidant enzymes, such as SOD, CAT, and POD, and increased enzymatic responses represent a greater antioxidant capacity for providing drought resistance in plants [43,44,45]. A remarkable increase in SOD activity in the leaves of low W and high N-treated seedlings (Fig. 3d; Additional file 3: Figure S3; Table 2) indicated that A. fruticosa plants improved their antioxidative defence status to protect themselves from ROS-mediated damage. A significant rise in CAT and POD activities in the leaves of A. fruticosa plants under low W, N, and P treatments also implied that A. fruticosa plants responded well to water and nutrient shortage by activating antioxidant capacity (Fig. 3e-f; Additional file 3: Figure S3; Table 2).
The maximum growth performance of A. fruticosa was observed in response to the treatments with W at 75.68% FC, N at 51.64 mg kg− 1 and P at 49.49 mg kg− 1 (Fig. 4). In accordance with our results, Wang et al.  also demonstrated that microbial inoculum with inorganic fertilizer (N at 61.73 mg kg− 1 and P at 40.0 mg kg− 1 doses) significantly improved morphological, physiological, and antioxidant enzyme responses in Medicago sativa when grown in coal-mined spoils.
The levels of doses developed for W, N and P using rsm correlated with the growth performance of A. fruticosa, suggesting that rsm could be used as an effective tool to design the perfect combination of water and fertilizers for revegetation of coal-spoiled soils. Our results revealed that application of various doses of W, N and P significantly altered growth-associated attributes, including shoot height, stem diameter, root length and dry biomass of A. fruticosa. Although all the three factors (W, N and P) exhibited significant impacts on morphological, physiological and biochemical traits of A. fruticosa, W was more prominent in influencing those traits when compared with N and P. The maximum effects on growth-contributing parameters were observed when high level of W (76% FC) was applied with 52.0 mg kg− 1 of N and 49.0 mg kg− 1 of P. This combination of W, N and P was also highly effective in enhancing the levels of LWP, gas exchange parameters and photosynthetic pigment contents. However, under low doses of W, which might have induced drought stress, A. fruticosa enhanced the accumulation of Pro and the activities of antioxidant enzymes SOD, POD and CAT. These results indicated that A. fruticosa employed osmoprotective effects of Pro to overcome osmotic stress and antioxidative roles of SOD, CAT and POD to reduce oxidative stress in order to survive under water-shortage conditions. Collectively, the findings of this study provided valuable and useful information that could be used as effective measures for revegetation and ecological restoration of coal-contaminated subsidence areas in China, and perhaps in other arid and semi-arid areas of the world. However, field studies with a range of plant species using suggested doses of W, N and P (76% FC, 52.0 mg kg− 1, 49.0 mg kg− 1, respectively) would be beneficial to ensure improved management and sustainable restoration of coal-spoiled areas.
Plant material and growth conditions
The coal-mined spoils with topsoil were collected from the Yangchangwan coal mining area (106° 35′ ~ 106° 38′E, 37° 59′ ~ 38° 03′N) of Lingwu city, Ningdong, Ningxia province, China. The sample was naturally dried, mixed, and finally sieved to 5 mm. Fourteen (14) kg of soil was taken to fill up each pot (upper diameter of 32 cm, bottom diameter of 27 cm, and height of 30 cm). One-year-old A. fruticosa seedlings with similar height (33.0 ± 3.0 cm) and stem diameter (2.8 ± 0.3 mm) were collected from Hengwang Seedlings Greening Company, Alashan County, Inner Mongolia, China. The plant species were identified based on its characteristics and was approved by the Jilantai Forestry Bureau, Alashan county, Inner Mongolia, China before collection. The collection of coal-mined spoils with topsoil and plant specimens were done in a public right of way, which means no specific permission was required for collecting soil and plant samples, and we have not deposited any voucher specimen. After transplanting one A. fruticosa seedling to each pot in March 2018, all pots were placed in an artificial shed at Northwest A&F University, (N 34° 16′, E 108° 4′), Yangling. The mean annual precipitation and temperature of the experimental site (outside of the artificial shed) were 650 mm and 13.7 °C, respectively. In the first month, sufficient water was provided to all pots to induce seedling recovery from transplantation stress. One month later, when the water content was close to the lower level of the water threshold, experimental treatments were applied, up to October 2018 according to the experimental design obtained from rsm. The amount of moisture lost from each pot through transpiration and evaporation was determined through a weighing method, and irrigation was done every day during the entire trial period. Fertilizer was applied in holes near the root zone of the seedlings. Nitrogen (N) fertilizer (urea-46% N) was applied in four times (¼ N, May 6th; ¼ N, June 6th; ¼ N, July 6th; ¼ N, August 6th), whereas phosphorus (P) fertilizer (Triple super phosphate-46% P2O5) was applied in two equal half (1/2 May 6th and 1/2 July 6th). Before starting the pot experiment, various soil properties were calculated, and recorded as FC (11.177%), soil bulk density (1.48 g cm− 3), soil pH (7.83), total N (0.157 g kg− 1), available N (6.14 mg kg− 1), available P (1.98 mg kg− 1) and available K (66.88 mg kg− 1).
Measurement of morphological growth
Changes in plant height and stem diameter were measured using measuring-tape and slide calliper, respectively, at every 15-days-interval throughout the growing seasons. Increase of growth in plant height and stem diameter were calculated by subtracting initial reading from final reading. For the measurement of plant biomass, plant samples were carefully washed with water for removing debris, and different plant parts (leaves, stems, and roots) were separated accordingly. Plant samples were initially dried in an oven for 20 min at 105 °C, followed by further drying at 80 °C, and finally dry biomass of different plant parts was determined. The maximum root length also measured using a measuring-tape. The root-shoot (R/S) biomass ratio was measured by below-ground biomass divided by above-ground biomass.
Determination of leaf gas-exchange parameters and leaf water potential
Leaf gas-exchange parameters, including net photosynthesis rate (Pn), transpiration rate (Tr), stomatal conductance (Gs) and water use efficiency (WUE) were measured on sunny days between 8:30 am and 11:30 am using a Portable Photosynthesis System CIRAS-3 (PP Systems, Amesbury, MA, U.S.A). Three seedlings from each treatment were used for determining gas-exchange parameters. At the time of measurement, photosynthetic active radiations were maintained at 1000 μmol m− 2 s− 1, CO2 at 300 μmol mol− 1, leaf temperature at 25 °c and relative humidity at 70% (provided by a built-in red LED light source). Leaf water potential (LWP) of the selected seedlings was determined by PMS-Model 1000 plant water potential meter (PMS Instrument Company, Albany, OR, USA) before 06:00 am.
Estimation of chlorophyll contents
Chlorophyll (Chl) was extracted from completely expanded fresh leaf samples (0.1 g) using a 10 mL mixture of ethanol, acetone and distilled water (4.5:4.5:1) following the method of Chen and Chen . The absorbance of the extracts was taken at 645 and 663 nm, and the contents of Chl a, Chl b and total Chl were calculated as mg g− 1 fresh weight (FW) by the formula described by Arnon .
Measurement of malondialdehyde and proline contents
The method of Wang et al.  was adopted to extract, and determine the contents of malondialdehyde (MDA) in leaf samples of A. fruticosa. Extract (2 mL) was mixed with an equal volume of 0.5% 2-thiobarbituric acid (TBA, dissolved in 15% trichloroacetic acid), and the mixture was heated for 30 min at 100 °C followed by cooling in an ice bath. The mixture was then centrifuged at 10,000 rpm for 10 min, and the absorbance of the supernatants was read at 450, 532 and 600 nm using a spectrophotometer. The MDA content was calculated and expressed as μmol g− 1 FW.
Free proline (Pro) from fresh leaf samples was extracted using sulphosalicylic acid solution, and measured with ninhydrin according to Bates et al. . The fresh leaf samples (0.1 g) were extracted with 3% aqueous sulfosalicylic acid (10 mL) and placed in water-bath for boiling for 10 min followed by filtering with Whatman No. 2 filter paper. Plant extract (2 mL) was mixed with 2 mL ninhydrin solution and 2 mL glacial acetic acid, and the mixture was heated for 30 min at 100 °C, followed by cooling in an ice bath. Toluene (4 mL) was used to extract proline from the mixture; the content of Pro was measured after taking the absorbance at 520 nm using a spectrophotometer.
Determination of antioxidant enzyme activities
Fresh leaf samples (0.3 g) were homogenized in a mortar and pestle (ice-cold conditions) using 8 mL of 50 mM sodium phosphate buffer (pH 7.8), and the homogenates were centrifuged at 10,000 rpm for 20 min at 4 °C. The supernatant was separated and used for the estimation of enzyme activities.
The activity of superoxide dismutase (SOD, EC 188.8.131.52) was determined following the method of Zhang et al. . The reaction mixture contained 3 mL of phosphate buffer (pH 7.8), 0.6 mL of 130 mM methionine buffer, 0.6 mL of 750 μM nitroblue tetrazolium buffer, 0.6 mL of 100 μM EDTA-Na buffer, 0.6 mL of 20 μM riboflavin and 0.2 mL of enzyme extract. The photoreduction of nitroblue tetrazolium was determined at 560 nm using a spectrophotometer. One unit (U) of SOD activity was equal to the amount of enzyme needed to produce 50% inhibition of the color reaction.
Catalase (CAT, EC 184.108.40.206) activity was determined following the method reported in Beers and Sizer . The reaction mixture (2.6 mL) contained 100 mM phosphate buffer (pH 7.0), 20 mM H2O2 and 100 μL enzyme extract. The reduction in H2O2 was monitored at 240 nm. Enzyme activity was expressed as CAT U min− 1 g− 1 FW.
Peroxidase (POD, EC 220.127.116.11) activity was carried out using the guaiacol oxidation method described by Ekmekci and Terzioglu . Enzyme extract (20 μL) was mixed with 3 mL of 50 mM phosphate buffer (pH 6.0) having 20.1 mM guaiacol, 12.3 mM H2O2, and guaiacol oxidation was measured based on absorbance increase at 470 nm for 3 min. Enzyme activity was expressed as units of guaiacol oxidized min− 1 g− 1 FW.
Experimental design for statistical analysis and optimization
In our study, response surface methodology (RSM) based on central composite design (CCD) method was used with statistical software design expert (Trial version 11, Stat-Ease Inc., USA). Ranges of independent variables of our study were chosen based on some previous literature [54,55,56] and presented in Table 4. Different levels (+α, + 1, 0, − 1, − α) of each factor were taken for operating this study. The distances from the centre to axial points are represented by α and calculated by the following equation.
Where k represents the number of factors, and therefore, α value in our study was 1.682. The CCD comprises 20 experimental runs with eight factorial points, six axial points and six replicates at the center points. A total of 80 seedlings subjected to twenty treatments with four repetitions were arranged randomly. Independent variables effects on responses were analyzed by using a second-order polynomial equation:
Where, Y is the response variables; a0 is regression constant-coefficient; a1, a2 and a3 are linear terms; a4, a5, and a6 are interaction terms; a7, a8, and a9 determine quadratic terms; x1, x2, and x3 = represents the coded value of W, N, and P.
RSM was carried out to optimize the process variables on the growth and development of A. fruticosa, and the experimental results are presented in Tables 1 and 2. The relationship between the response variables and process variables was obtained by applying multiple regression analysis. The final models, in terms of coded factors were represented in Table 3. A Positive and negative sign in front of the terms indicated synergistic and antagonistic effects, respectively. To determine the significance and fitness of the model for response variables, analysis of variances (ANOVA) was carried out. The model F and P-value were found for different parameters indicated that the model was highly significant. To determine the model precision, the results of the experiment were examined with coefficient of variation (CV), adequate precision (AP) and R2 values. In our present experiment, the adjusted R-squared (R2a) value was detected to be very near and a little bit lesser then R2, which designated a strong correlation between the experimental and predicted values. For reasonable agreement, the difference between predicted R-squared (R2p) value and R2a value should be less than 0.2 , and we found that all difference values were less than 0.2 which confirmed that the form of the model chosen to explain the relationship between the factors and the response was well-correlated. 3-D response surface plots obtained from statistical software were designed to find out the relationship between process variables and to assess the response variables over process variables. At the time of creating 3-D response surface plot, we retained our third factor as a central level (as our model comprises three factors).
Availability of data and materials
The dataset generated and analyzed during the study are available from the corresponding author on reasonable request.
Root to shoot biomass ratio
Leaf water potential
Water use efficiency
Central composite design
Response surface methodology
Valladares F, Dobarro I, Sánchez-Gómez D, Pearcy RW. Photoinhibition and drought in Mediterranean woody saplings: scaling effects and interactions in sun and shade phenotypes. J Exp Bot. 2005;56:483–94.
Liu Y, Gao M, Wu W, Tanveer SK, Wen X, Liao Y. The effects of conservation tillage practices on the soil water-holding capacity of a non-irrigated apple orchard in the loess plateau, China. Soil Tillage Res. 2013;130:7–12.
Xiao W, Fu Y, Wang T, Lv X. Effects of land use transitions due to underground coal mining on ecosystem services in high groundwater table areas: a case study in the Yanzhou coalfield. Land Use Policy. 2018;71:213–21.
Brigić A, Vujčić-Karlo S, Kepčija RM, Stančić Z, Alegro A, Ternjej I. Taxon specific response of carabids (Coleoptera, Carabidae) and other soil invertebrate taxa on invasive plant Amorpha fruticosa in wetlands. Biol Invasions. 2014;16:1497–514.
Kong J, Pei Z, Du M, Sun G, Zhang X. Effects of arbuscular mycorrhizal fungi on the drought resistance of the mining area repair plant Sainfoin. Int J Min Sci Technol. 2014;24:485–9. https://doi.org/10.1016/j.ijmst.2014.05.011.
Devnarain N, Crampton BG, Chikwamba R, Becker JVW, O’Kennedy MM. Physiological responses of selected African sorghum landraces to progressive water stress and re-watering. South African J Bot. 2016;103:61–9.
Andivia E, Fernández M, Vázquez-Piqué J. Autumn fertilization of Quercus ilex ssp. ballota (Desf.) Samp. Nursery seedlings: effects on morpho-physiology and field performance. Ann For Sci. 2011;68:543–53.
Chiappero J, del Rosario Cappellari L, Sosa Alderete LG, Palermo TB, Banchio E. Plant growth promoting rhizobacteria improve the antioxidant status in Mentha piperita grown under drought stress leading to an enhancement of plant growth and total phenolic content. Ind Crops Prod. 2019;139:111553. https://doi.org/10.1016/j.indcrop.2019.111553.
Deng B, Du W, Liu C, Sun W, Tian S, Dong H. Antioxidant response to drought, cold and nutrient stress in two ploidy levels of tobacco plants: low resource requirement confers polytolerance in polyploids? Plant Growth Regul. 2012;66:37–47.
Hessini K, Ghandour M, Albouchi A, Soltani A, Werner KH, Abdelly C. Biomass production, photosynthesis, and leaf water relations of Spartina alterniflora under moderate water stress. J Plant Res. 2008;121:311–8.
Keshavarz Afshar R, Chaichi MR, Moghadam H, Ehteshami SMR. Irrigation, phosphorus fertilizer and phosphorus solubilizing microorganism effects on yield and forage quality of turnip (Brassica rapa L.) in an arid region of Iran. Agric Res. 2012;1:370–8.
Pan J, Liu Y, Zhong X, Lampayan RM, Singleton GR, Huang N, et al. Grain yield, water productivity and nitrogen use efficiency of rice under different water management and fertilizer-N inputs in South China. Agric Water Manag. 2017;184:191–200.
Munshaw GC, Ervin EH, Shang C, Askew SD, Zhang X, Lemus RW. Influence of late-season iron, nitrogen, and seaweed extract on fall color retention and cold tolerance of four bermudagrass cultivars. Crop Sci. 2006;46:273–83.
Saneoka H, Moghaieb REA, Premachandra GS, Fujita K. Nitrogen nutrition and water stress effects on cell membrane stability and leaf water relations in Agrostis palustris Huds. Environ Exp Bot. 2004;52:131–8.
Molinari HBC, Marur CJ, Daros E, De Campos MKF, De Carvalho JFRP, Filho JCB, et al. Evaluation of the stress-inducible production of proline in transgenic sugarcane (Saccharum spp.): osmotic adjustment, chlorophyll fluorescence and oxidative stress. Physiol Plant. 2007;130:218–29.
Shao H, Liang Z, Shao M. Changes of anti-oxidative enzymes and MDA content under soil water deficits among 10 wheat (Triticum aestivum L.) genotypes at maturation stage. Colloids Surf B Biointerfaces. 2005;45:7–13.
Bhatt D, Negi M, Sharma P, Saxena SC, Dobriyal AK, Arora S. Responses to drought induced oxidative stress in five finger millet varieties differing in their geographical distribution. Physiol Mol Biol Plants. 2011;17:347–53.
Wang J, You Y, Chen W, Xu Q, Wang J, Liu Y, et al. Optimal hypobaric treatment delays ripening of honey peach fruit via increasing endogenous energy status and enhancing antioxidant defence systems during storage. Postharvest Biol Technol. 2015;101:1–9. https://doi.org/10.1016/j.postharvbio.2014.11.004.
Zhao R, Guo W, Bi N, Guo J, Wang L, Zhao J, et al. Arbuscular mycorrhizal fungi affect the growth, nutrient uptake and water status of maize (Zea mays L.) grown in two types of coal mine spoils under drought stress. Appl Soil Ecol. 2015;88:41–9.
Huang J, Yu H, Lin H, Zhang Y, Searle EB, Yuan Z. Phosphorus amendment mitigates nitrogen addition-induced phosphorus limitation in two plant species in a desert steppe, China. Plant Soil. 2016;399:221–32.
Mourabet M, El Rhilassi A, El Boujaady H, Bennani-Ziatni M, El Hamri R, Taitai A. Removal of fluoride from aqueous solution by adsorption on Apatitic tricalcium phosphate using box-Behnken design and desirability function. Appl Surf Sci. 2012;258:4402–10.
The authors express their special gratitude to all the funding sources for the financial assistance and also want to thank Feng shulin, Boyuan Li, Xu Ma and Guoli Lv for extending help during sample analysis. We are very grateful to the editor and anonymous reviewers, for very valuable comments that have greatly helped us to further improve the quality of our manuscript.
This research was financed by National Key Research and Development Program of China “Eco-security technology for coal mining bases in the Northwestern arid desert regions in China” (2017YFC0504400); “Studies on the vegetation rehabilitation and conservation in abandoned coal mining land research” (2017YFC0504402). National Natural Science Foundation Project of China (31670713) and Yulin City Forestry Science and Technology Plan Projects (K403021528). The supporters did not play any role in the design of the study, data collection, and analysis of the data or writing the manuscript.
Authors and Affiliations
College of Natural Resources and Environment, Northwest A&F University, Yangling, Shaanxi, 712100, People’s Republic of China
Rana Roy, Jinxin Wang & Ashim Sikdar
Department of Agroforestry & Environmental Science, Sylhet Agricultural University, Sylhet, 3100, Bangladesh
Rana Roy & Ashim Sikdar
Department of Biochemistry and Molecular Biology, Bangabandhu Sheikh Mujibur Rahman Agricultural University, Gazipur, 1706, Bangladesh
Mohammad Golam Mostofa
Key Laboratory of Plant Nutrition and the Agri-Environment in Northwest China, Ministry of Agriculture, Yangling, 712100, People’s Republic of China
School of Economics and Finance, Xi’an Jiaotong University, Xi’an, 710049, People’s Republic of China
RR and JXW conceived and designed the experiment; RR performed the experiments; RR and MGM analyzed the data; AS and TS assisted with the experiments; RR, JXW and MGM wrote and revised the manuscript. All authors have read and approved the final manuscript.
Response surface plots showing the effect of soil-water (W), nitrogen (N) and phosphorus (P) on the plant height (a-b), stem diameter (c-d), root length (e-f), dry biomass (g), and root-shoot (R/S) biomass ratio (h-i).
Response surface plots showing the effect of soil-water (W), nitrogen (N) and phosphorus (P) on the leaf water potential (LWP) (a-b), photosynthesis rate (Pn) (c), transpiration rate (Tr) (d-e), stomatal conductance (Gs) (f-g), water use efficiency (WUE) (h-i), chlorophyll a (Chl a) (j-l), chlorophyll b (Chl b) (m-o) and total chlorophyll (Total Chl) content (p).
Response surface plots showing the effect of soil-water (W), nitrogen (N) and phosphorus (P) on the malondialdehyde (MDA) (a) and proline (Pro) (b-c) contents, and activities of superoxide dismutase (SOD) (d-e), catalase (CAT) (f-g), and peroxidase (POD) (h-i).
Rights and permissions
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://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Roy, R., Mostofa, M.G., Wang, J. et al. Improvement of growth performance of Amorpha fruticosa under contrasting regime of water and fertilizer in coal-contaminated spoils using response surface methodology.
BMC Plant Biol20, 181 (2020). https://doi.org/10.1186/s12870-020-02397-1