Next Article in Journal
Biotic versus Abiotic Controls on Bioavailable Soil Organic Carbon
Next Article in Special Issue
Soil Processes, Pedofeatures and Microscale Metal Distributions: Relevant Study of Contaminant-Dynamics Calls for Pedology-Based Soil-Depth Sampling Strategies
Previous Article in Journal
Nuclear Magnetic Resonance Analysis of Changes in Dissolved Organic Matter Composition with Successive Layering on Clay Mineral Surfaces
Previous Article in Special Issue
Altered Bacterial Communities in Long-Term No-Till Soils Associated with Stratification of Soluble Aluminum and Soil pH
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Soil Moisture Status on Soil Cadmium Phytoavailability and Accumulation in Plantain (Plantago lanceolata)

Institute of Agriculture & Environment, Massey University Manawatu, Private Bag 11 222, Palmerston North 4442, New Zealand
*
Author to whom correspondence should be addressed.
Soil Syst. 2018, 2(1), 9; https://doi.org/10.3390/soils2010009
Submission received: 29 November 2017 / Revised: 28 January 2018 / Accepted: 7 February 2018 / Published: 12 February 2018
(This article belongs to the Special Issue Soil Processes Controlling Contaminant Dynamics)

Abstract

:
The effect of fluctuating soil moisture cycles on soil cadmium (Cd) phytoavailability was investigated in a pot trial with two contrasting soils (Kereone (Allophanic), total Cd 0.79 mg kg−1; and Topehaehae (Gley), total Cd 0.61 mg kg−1) that were either sown with plantain (Plantago lanceolata) or left unseeded. Varying soil moisture contents were established using contrasting irrigation regimes: “flooded” (3 days flooded and then 11 days drained); or “non-flooded” (irrigation to 70% of potted field capacity every 7 days). Overall, there was no significant difference in mean 0.05 M CaCl2 soil extractable Cd concentrations or plant tissue Cd concentrations between flooded and non-flooded irrigation. However, there was a consistent trend for an increase in soil extractable Cd concentrations following irrigation, regardless of the irrigation regime. Mean soil extractable Cd and plant tissue Cd concentrations were significantly greater (approximately 325% and 183%, respectively) for the Topehaehae soil than the Kereone soil, despite the lower soil total Cd concentration of the Topehaehae soil. These results indicate that Cd solubility is sensitive to increases in soil moisture following periods of soil drainage, but insensitive to short-term periods of soil saturation. Plant tissue Cd concentrations in Cd-sensitive forage crops such as plantain are likely to be greater following large rainfall events over summer and autumn. This has the potential to increase animal dietary Cd exposure and rate of liver/kidney Cd accumulation.

1. Introduction

Cadmium (Cd) is a non-essential element that is readily taken up by plants. Below phytotoxic concentrations, plants have little ability to regulate Cd uptake and accumulation [1]. Because of this, plant tissue Cd concentrations tend to increase with increasing soil Cd phytoavailability [1,2,3].
Many soil factors are known to influence soil Cd phytoavailability. These include total soil Cd concentration, pH, organic matter content, mineralogy (i.e., clay types and abundance), the concentration of other anions and cations in soil solution, and soil redox potential [4,5,6,7,8,9]. Of these factors, soil pH is thought to be the most dominant factor controlling soil Cd phytoavailability, largely due to its overriding influence on the surface charge of variable charge clays and organic matter that dictates the soil’s net negative charge and therefore Cd2+ sorption potential [4,5,8,10,11,12,13].
Many authors have shown that moisture-driven changes in soil redox potential can also have a large influence on metal solubility [14,15,16,17,18]. Most research related to soil moisture-mediated changes in Cd phytoavailability has focused on prolonged periods of soil saturation (e.g., cultivation of rice), where the development of strongly-reducing soil conditions leads to a decrease in Cd solubility through the formation of insoluble Cd-sulphide complexes [19,20,21]. However, prolonged soil saturation is atypical for the majority of agricultural soil and crop management scenarios. For example, agricultural/horticultural production systems in New Zealand are based on either naturally well-drained soils or soils modified with artificial-drainage systems to control soil moisture within a desired range [22,23]. As a consequence, these soils are subject to fluctuating soil moisture and only temporary periods (i.e., several days) of soil saturation that occur when rainfall (and/or irrigation) exceeds the rate of soil drainage.
Unlike under prolonged anaerobic conditions, little is currently known regarding the behaviour of Cd under fluctuating soil moisture conditions punctuated by only short periods of saturation. Numerous studies have reported increases in the solubility of other metals following soil saturation and the associated onset of weakly- to moderately-reducing soil conditions [16,24,25]. This phenomenon is related to reduction–oxidation processes that govern the dissolution and precipitation of manganese (Mn) and iron (Fe) hydrous oxides. Under well-drained aerobic soil conditions, relatively insoluble Fe3+ and Mn4+ hydrous oxides prevail [22,24], providing sorption surfaces for many metals, including Cd [4,26,27,28]. However, under anaerobic soil conditions, Mn4+ and Fe3+ become reduced to their more soluble Mn2+ and Fe2+ oxidation states [22,24]. With the consequent dissolution of these Mn and Fe hydrous oxides, co-sorbed metal ions are released into soil solution, increasing metal phytoavailability. Notably, Ponnamperuma [29] reported that reducing conditions can be established very rapidly (within a few days), suggesting that even short periods of soil saturation can influence the solubility and phytoavailability of metals sorbed to Fe/Mn hydrous oxides. This observation is supported by the work of Bhatti et al. [30], who attributed an increase in plant tissue As concentration to the dissolution of Fe hydrous oxides for plants grown under a cyclical 3-day flooded/11-day drained irrigation regime, relative to plants grown under irrigated but non-saturated soil conditions.
Due to the affinity between Cd2+ and Fe/Mn hydrous oxides, there is a possibility that Cd phytoavailability could also be influenced by short-term periods of soil saturation, similar to the effects on As described by Bhatti et al. [30]. This hypothesis facilitates the successful examination of soil extractable Cd concentration differences over time (in addition to differences in plant Cd accumulation), since these should directly relate to the contrasting soil hydrological conditions (and redox potentials) imposed. The research described in this paper was therefore undertaken to quantify the effect of fluctuating soil moisture content and short-term periods of soil saturation on soil Cd phytoavailability and plant tissue Cd accumulation. Improving our understanding of the impact of soil moisture fluctuations on soil Cd phytoavailability may provide greater insight into appropriate strategies to better manage Cd accumulation in agricultural plant species.

2. Materials and Methods

The research described in this paper was conducted as a pot-trial, under controlled conditions, over the period April–September 2015. Two irrigation/soil moisture regimes were overlain on two contrasting soils (Kereone fine sandy loam and Topehaehae sandy clay loam; Table 1), which were either sown in plantain (Plantago lanceolata, cultivar “Ceres Tonic”) or left as an unplanted control. Five replicates of each treatment were used (i.e., 40 pots in total). Plantain was chosen as the Cd phytoavailability indicator species in this trial based on previous research [31], which recorded high tissue Cd concentrations in this cultivar. Although chicory was shown to accumulate even greater tissue Cd concentrations than plantain [31] “Ceres Tonic” plantain is known to have greater cool season growth rate [32]—an important factor given that this trial was carried out over winter.
The imposed irrigation strategies (based on the methodology of Bhatti et al. [30]) were (i) “non-flooded”, where soil moisture was maintained at 70% potted field capacity (FcP) utilising a 7-day irrigation cycle; and (ii) “flooded”, where soil moisture was maintained at 110% FcP for 3 days, subsequent to which drainage was permitted, allowing aeration of the soil over the ensuing 11 days, before the 14-day irrigation cycle was repeated.

2.1. Pot Design

To facilitate the proposed irrigation/soil moisture strategies, square white plastic pots (16.5 × 16.5 × 19 cm) were prepared based on the methodology of Bhatti et al. [30]. Briefly, for flood-irrigated pots, a 7 mm diameter hole was drilled into the side of the pot approximately 2 cm above the base. Into this, a 30 cm length of rubber tube (external diameter 8 mm and internal diameter 4 mm) was inserted such that 15 cm of tube extended across the internal width of the pot and 15 cm extended externally (Figure 1). To facilitate flooded soil conditions for 3 days, the external length of rubber tube was bent over itself and clamped to prevent drainage. During the subsequent 11-day drained period, the tube was unclamped to facilitate soil drainage and aeration. Pots used for non-flooded conditions had two 3 mm diameter holes drilled into the base of each side (i.e., eight holes per pot) to permit drainage and soil aeration.
A 20 × 20 cm section of non-woven mulch mat was placed in the bottom of each non-flooded treatment pot. This water-permeable fabric extended approximately 2 cm up the internal wall of each pot, preventing soil loss through or blockage of the drainage holes located near the base of each side. For flooded pots, the section of rubber tube inside the pot had a narrow slit (12 cm × 0.3 cm) cut out along the length of the tube. This internal rubber tube was then wrapped in mulch mat to prevent potential sediment blockage or loss. The bottom of all pots was then filled with 400 g of gravel (6–13 mm) to facilitate aeration and water movement around the drainage holes. All pots were then weighed to determine mass without soil added.

2.2. Soils

Soils were collected from a commercial dairy farm in the Waikato region of New Zealand. Classification and chemical characteristics of the two soils is provided in Table 1. At the collection sites, the top 3 cm of soil was scraped off to remove the dense root mat of the ryegrass/white clover pasture sward. Soil from the A-horizon was then retrieved from 3 to 15 cm depth and sieved (field-moist) to 4 mm. Each pot was packed with 3500 g of field moist soil, with a subsample of each soil taken for moisture content determination (24 h at 105 °C). The dry-soil gravimetric moisture content at potting was calculated at 0.42 and 0.43 for the Topehaehae and Kereone soils, respectively. On this basis, the oven-dry corrected mass of soil in each pot was 2464 and 2451 g for the Topehaehae and Kereone soils, respectively.
All pots were then watered until the soil was saturated. Once drainage had ceased, pots were re-weighed to determine the mass of water in each pot at FcP. This mass of water at FcP was then multiplied by either 70% or 110%, and added to the combined mass of the oven-dry soil and pot apparatus, so as to determine the target total pot mass for each soil type by irrigation treatment.

2.3. Experimental Layout and Management

All pots were laid out in 4 rows of 10 pots in completely randomised design (CRD). The trial was carried out under a clear Perspex-roofed shelter at ambient air temperature (ranging from 0.1 °C to 21.7 °C) in Tauranga, Bay of Plenty. Daily mean, minimum, and maximum air temperature data for Tauranga over the experimental period is shown in Figure 2. Throughout the duration of the trial, pots were rotated on a weekly basis to ensure there was no environmental bias for any replicate.
Plantain seeds were sown in the appropriate pots on 20 April 2015 (day −20). Divergence in irrigation strategy occurred on 10 May 2015 (day 0), once seedlings had germinated, become well established, and had been thinned to one healthy plant per pot. The trial was completed on 13 September 2015 (day 126).

2.4. Soil Sampling and Plant Harvest

Immediately prior to the initiation of the different irrigation strategies (day 0), a single soil core (100 mm depth, 18 mm diameter) was extracted from each pot. Further soil samples were taken 3 days later, timed in alignment with the end of the 3-day “saturated” phase of the flooded treatments. Flooded treatments were allowed to drain over the subsequent 11 days, before this “pre-irrigation” and “+3 day” sampling strategy was repeated, coinciding with each 14-day flood irrigation cycle. In total, this fortnightly soil sampling cycle was repeated nine times (i.e., 18 paired samples). Soil removed at each sampling was backfilled with surplus soil from each soil type that had been sieved and kept in storage. Soil sampling sites were rotated clockwise around each pot, with the location of backfilled soil sampling sites marked with a toothpick.
Plantain was harvested from all pots when “grazing-maturity” was reached within the treatments exhibiting the greatest growth rates (non-flooded Kereone and Topehaehae treatments). This occurred on day 126. Plant tissue samples were harvested 20 mm above the soil surface, weighed, washed, and oven-dried (60 °C for 3 days) for dry matter content (% DM), dry matter yield (g DM pot−1) and tissue Cd (mg kg−1 DM) assessment.

2.5. Chemical Analysis

Following the procedure of Li et al. [33], 6 g of field moist soil was extracted in 30 mL of 0.05 M CaCl2 with 24 h end-over-end shaking. The samples were then centrifuged at 11,953× g for 10 min before filtering (Whatman No. 42). Plant tissue samples were crushed and homogenised prior to analysis for Cd, following the methodology described by Jeyakumar et al. [34] for plant metal analysis. Briefly, 0.1 g of each sample was digested using 10 mL of 69% HNO3. Cadmium in the soil extracts and plant tissue digests was assessed using a graphite furnace atomic absorption spectrophotometer (GFAAS; Perkin-Elmer AAnalyst 600, Waltham, MA, USA).
The pH of 0.05 M CaCl2 soil extracts was assessed for all samples at sampling days 0, 28, 84, and 112. This was undertaken using a Eutech Instruments CyberScan pH 310 fitted with an Ag/AgCl oxidation reduction potential reference electrode.
For each treatment, the mean soil extractable Cd concentration was determined using replicate data across all sampling occasions. Treatment-mean soil extractable Cd concentrations were also independently assessed for samples collected at pre-irrigation and +3 day sampling timings. This also allowed the mean change in soil extractable Cd concentrations between pre-irrigation and +3 day samplings to be determined on a treatments basis, to test the hypothesis that this variance should be greater in flooded than non-flooded treatments.

2.6. Quality Control

All chemicals used in the experiments were of analytical grade. The GFAAS limit of detection for Cd in digest solution was 0.002 mg L−1. Analytical accuracy was assessed through parallel analysis of a certified standard reference material (NIST 1573a, tomato leaves). The analysed concentration differed by less than 5% of the certified value. GFAAS calibration rescale was performed every 30 samples.

2.7. Statistical Analysis

Statistical analysis was carried out using IBM SPSS Statistics for Windows, version 22 (IBM Corp., Armonk, NY, USA). Significant outliers were identified using Tukey’s outlier filtering method [35]. Analysis of variance (ANOVA) was undertaken to assess differences in treatment means, and Tukey’s honestly significant difference (HSD) was run as a post-hoc analysis to identify differences between individual treatment means.

3. Results

3.1. Soil Extractable Cd

The mean soil extractable Cd concentration across the duration of the trial was significantly (p < 0.001) different between the two soils (Table 2), with soil extractable Cd concentrations being approximately four-times greater in the Topehaehae soil (85 µg kg−1) than the Kereone soil (20 µg kg−1). However, there was no significant effect (p > 0.05) of irrigation regime or of plant on extractable Cd concentration within the same soil.
When samples were separated into those taken at pre-irrigation and at the +3 day timings (Table 2), there was a consistent trend whereby mean soil extractable Cd concentrations were greater for the +3 day samples relative to the pre-irrigation samples; however, these differences only reached significance (p < 0.05) for treatments within the Kereone soil. Within the Kereone soil, the difference in soil extractable Cd concentrations between pre-irrigation and +3 day sampling timings was of greater significance for flooded treatments (p < 0.01) than the non-flooded treatments.

Soil Extractable Cd over Time

Figure 3 presents a graphical comparison of the mean soil extractable Cd concentration at each sampling date for the flooded and non-flooded treatments of the Topehaehae planted, Topehaehae unplanted, Kereone planted, and Kereone unplanted treatments (Figure 3a–d, respectively). The trend for greater soil extractable Cd concentrations at +3 day samplings relative to the pre-irrigation samplings in both flooded and non-flooded treatments is illustrated in these graphs.
Within all treatments, there was a general increase in soil extractable Cd concentrations over the first 4–6 weeks of the trial. This phase was then followed by a decrease in soil extractable Cd concentration, with the Topehaehae soil in particular showing an abrupt 50–60% decrease over days 42–56 (Figure 3a,b). Although not as abrupt, a similar trend was seen for the Kereone soil (Figure 3c,d), where an approximate 60% decrease in soil extractable Cd concentration occurred over the latter half of the trial. This decline was not an effect of plant uptake, since soil extractable Cd concentrations in non-planted treatments showed the exact same patterns as the planted treatments throughout the duration of the trial.

3.2. Soil Solution pH

Periodic assessment of soil solution pH throughout the duration of the trial indicated that with the exception of the Kereone planted flooded treatment at day 28, there was no consistent significant (p > 0.05) effect of irrigation regime on solution pH (Table 3). Overall, there was a significant (p < 0.001) increase in soil solution pH throughout the course of the trial.

3.3. Plantain Tissue Cd Concentration, Yield, and Cd Uptake

Significant (p < 0.001) differences in overall treatment mean plantain tissue Cd concentrations were apparent (Table 2), with mean tissue Cd concentrations decreasing in the order Topehaehae non-flooded > Topehaehae flooded > Kereone flooded > Kereone non-flooded. However, the main effect was that of soil type (p < 0.001), with plant tissue Cd concentrations for the Topehaehae soil being approximately double those for the Kereone soil. In contrast, the main effect of irrigation regime on tissue Cd concentration was not significant (p > 0.05).
Dry matter yield of plantain was significantly (p < 0.001; Table 2) influenced by irrigation regime, however there was no relationship (p > 0.05) to soil type. Within the Kereone soil, there was an inverse relationship between plant yield and tissue Cd concentration. As a consequence, total Cd uptake was not significantly (p > 0.05) different for the flooded and non-flooded treatments. In contrast, in the Topehaehae soil, plant tissue Cd concentrations were not influenced (p > 0.05) by differences in plant yield between the two irrigation regimes. As a consequence, total Cd uptake was significantly (p < 0.05) greater in the higher-yielding non-flooded treatment.

4. Discussion

4.1. Effect of Fluctuating Soil Moisture on Soil Cd Phytoavailability

The results indicate that the solubility of soil Cd is sensitive to initial increases in soil moisture content following periods of soil drainage, but insensitive to short-term periods of soil saturation. This is shown by the consistent increases in soil extractable Cd concentrations measured at the +3 day sampling timing relative to pre-irrigation, regardless of irrigation strategy (Table 2 and Figure 3). In addition, there was a lack of difference in plant tissue Cd concentrations between flooded and non-flooded irrigation regimes in our trial (Table 2), providing further evidence that Cd phytoavailability is insensitive to short-term periods of soil saturation. This is an important finding, given that previous research [16,21,24,25,26] has indicated that changes in soil redox potential can influence Cd solubility.
Although we did not directly assesses soil redox potential, it is known that soil redox potential can change within a few days of saturation [29]. As change in redox potential is expected to influence soil solution pH [29,36,37,38,39], the lack of difference in soil solution pH between irrigation regimes in our study (Table 3) suggests that either (i) both soils buffered potential pH changes brought about by differences in redox processes driven by the different irrigation regimes; or (ii) that the 3-day flooded period was insufficient for generating reducing soil conditions. Notably, Bhatti et al. [30] employed a very similar irrigation methodology to that used in our study, and reported significant increases in plant As concentration under flooded compared to non-flooded irrigation treatments, which they attributed to the dissolution of Fe hydrous oxides under the short-term flooded conditions. Although Bhatti et al. [30] did not directly measure soil redox potential, their results indicate that the solubility of As is much more sensitive than Cd to short periods of soil saturation.
Potentially, the lack of difference in Cd solubility between the two irrigation regimes in this study could also be related to soil pH. It is known that the dissolution of Mn hydrous oxides can occur even under oxidising conditions when soil pH is mildly acidic [36]. Given that both the Topehaehae and Kereone soils had mildly acidic soil pH at the time of soil collection (pH 5.6 and 6.3, respectively; Table 1), it is possible that dissolution of Mn hydrous oxides (and concurrent release of Cd) could have occurred with increased soil moisture, regardless of irrigation regime.

4.2. Change in Soil Extractable Cd as a Function of Time

This study is the first to evaluate changes in Cd solubility over time under contrasting wetting/drying soil moisture cycles. The study recorded an increase in soil extractable Cd concentration over the first 4–6 weeks of the trial, followed by a large (50–60%) decrease thereafter, which was consistent for both irrigation regimes and soil types (Figure 3). Data presented in Figure 3 indicate that in the initial 4–6 weeks of the trial, the increase in soil extractable Cd concentrations that occurred between each pre-irrigation and +3 day sampling was not reversed (or completely reversed) over the subsequent drained phase, resulting in successively greater soil extractable Cd concentrations at each subsequent pre-irrigation sampling. With relatively dry soil conditions at the time of soil collection (gravimetric moisture content (θm) for both soils of ~0.42 g g−1, compared to θm at FcP of 0.64 and 0.73 g g−1 for the Topehaehae and Kereone soils, respectively), it is possible that the trend for increasing soil extractable Cd concentrations over the first 4–6 weeks of the trial was a consequence of the abrupt increase in soil moisture. This is consistent with previous research [40] that showed a large increase in soil solution Cd2+ and total Cd concentrations (approximately 100% and 50%, respectively) in the first 48 h following rewetting of air-dried soil samples, which the authors attributed to the establishment of anaerobic soil microbial respiration processes. Furthermore, as soil organic matter is an important sink for Cd in New Zealand agricultural soils [41] it is also possible that mineralisation of soil organic matter following sieving/potting may have reduced the soils’ Cd sorption potential, thereby increasing labile-Cd during the first 4–6 weeks of this study.
The consistent decline in soil extractable Cd concentrations that occurred beyond days 42–56 coincided with a small but significant (p < 0.001) increase in the pH of 0.05 M CaCl2 soil extracts that occurred relative to that measured at day 0 (Table 3). With an increase in pH, specific sorption of Cd would be expected to increase with a corresponding decrease in Cd phytoavailability [4,8,11,13]. The increase in soil solution pH is therefore consistent with the decrease in soil extractable Cd concentrations measured mid-way through the trial. This does not contradict the increase in soil extractable Cd concentrations observed during the initial 4–6 weeks (Figure 3), since soil solution pH during this period did not change significantly (p > 0.05) (Table 3).
Alternatively, the decline in soil extractable Cd concentrations mid-way through the trial may have been temperature-related. The decline aligned with the coldest period of winter, when overnight temperatures regularly approached 0 °C (Figure 2). With a decrease in temperature, a decrease in the decomposition rate of soil organic carbon is expected [42]. As a consequence of this decrease in energy availability, lower respiratory oxygen demand is likely to result in a decrease in microbially-induced reduction [43]. Hindersmann and Mansfeldt [24] reported an inverse relationship between temperature and the time taken to transition from oxidising to weakly reducing soil conditions, noting an increase from 10 to 160 h at temperatures of 15 °C and 7 °C, respectively. In our trial, with the maximum soil saturation period being only 3 days, it may be that oxidising conditions remained prevalent over the cooler winter period in all treatments. Consequent precipitation of Mn and Fe oxides under oxidising soil conditions would result in greater specific sorption of Cd [24,26], possibly explaining the overall decline in soil extractable Cd concentrations beyond day 42. Direct assessment of microbial respiration rate and/or redox potential was outside the scope of this project, being designed to follow irrigation/drainage methodology similar to that employed in previous research [30] to explore changes in metal phytoavailability with fluctuating soil moisture.

4.3. Effect of Soil Type on Cd Phytoavailability and Uptake

Despite soil total Cd concentration being 30% greater in the Kereone soil than the Topehaehae soil (0.79 vs. 0.61 mg kg−1, respectively; Table 1), soil extractable Cd concentrations were around four times greater in the Topehaehae soil (Table 2). Correspondingly, tissue Cd concentrations were also greater (approximately double) for plantain grown in the Topehaehae soil than the Kereone soil (Table 2).
Although an inverse relationship between tissue Cd concentration and plant yield is frequently reported in the literature [44,45,46], this effect was only evident in the Kereone soil, resulting in equivalent plant Cd uptake between low-yielding flooded (1.020 µg Cd pot−1) and high-yielding non-flooded (1.001 µg Cd pot−1) irrigation regimes. In contrast, the greater extractable Cd concentrations of the Topehaehae soil sustained plant tissue Cd concentration irrespective of the differing plant growth rate under the two irrigation regimes, resulting in much greater plant Cd uptake for the high-yielding non-flooded treatment (2.602 µg Cd pot−1) than the low-yielding flooded treatment (0.880 µg Cd pot−1) (Table 2). Our observation that an inverse relationship between plant tissue Cd concentration and plant yield only occurred at low Cd phytoavailability is consistent with the finding of Mortvedt et al. [45].
The differences in Cd phytoavailability between the two soil types can be explained by differences in soil characteristics. Although cation exchange capacity, Mn oxide, and Fe and Al amorphous oxide content have each been shown to be correlated with soil Cd desorption, soil total Cd, pH, and organic matter content are considered to be the dominant factors regulating Cd solubility and desorption [28]. Of these three factors, total Cd is positively correlated while soil pH and organic matter are negatively correlated [9,28]. The overriding effect of soil pH on metal phytoavailability is regularly cited in the literature—particularly for soils with high variable charge capacity, such as those in this study [4,12,13]. Notably, the pH of the Topehaehae soil was 0.7 pH units lower than that of the Kereone soil at the outset of the trial. The greater pH of the Kereone soil, coupled with its greater organic matter content and predominance of variable charge clays (especially allophane and Mn oxides) is likely to explain the much lower extractable Cd concentrations of this soil type.
The importance of soil pH and organic matter content in regulating soil Cd solubility was demonstrated by Gray et al. [28], who developed an algorithm to predict soluble Cd concentrations based on these three factors, which accounted for 81% of the variation in soil solution Cd concentrations of 29 New Zealand soils. This algorithm was used to predict Cd solubility for the soils in our study, using the soil pH, carbon, and total Cd values in Table 1 (carbon being converted to soil organic matter using the Van Bremmelen factor of 1.724 [47]). Soil solution Cd concentrations were predicted to be almost three times greater for the Topehaehae soil than the Kereone soil. In relative terms, this is very similar to the approximate three-fold difference in mean soil extractable Cd concentrations measured for the Topehaehae and Kereone soils (i.e., ~87 µg kg−1 vs. ~20 µg kg−1, respectively; Table 2), demonstrating the importance of soil pH and organic matter content in regulating the phytoavailability of soil total Cd.

4.4. Implication of These Findings for Cd Management in Agricultural Soils

The trend shown in our trial for increases in soil extractable Cd concentration with increases in soil moisture is a new and potentially important finding with regard to Cd management in agricultural soils. This observation suggests that after a period of drainage or soil drying, soil Cd phytoavailability is extremely sensitive to subsequent rainfall or irrigation inputs and indicates that Cd phytoavailability increases with increasing soil moisture prior to the point of saturation and the associated onset of reducing soil conditions.
This new knowledge may provide greater insight into appropriate soil Cd management strategies for minimising livestock dietary Cd exposure. For example, under typical New Zealand non-irrigated pastoral agriculture conditions, chicory and plantain crops are planted to provide quality feed when ryegrass-white clover-based pastures are compromised by hot, dry summer to early-autumn conditions [48]. As chicory and plantain have been shown to be ‘Cd-accumulator’ species [31] it is likely that their tissue Cd concentrations will be enriched following rainfall events that substantially increase soil moisture during the predominantly dry summer and autumn soil conditions. This effect has been observed in recent field trials investigating Cd accumulation in chicory (currently unpublished). Consistent across two contrasting soils, mean tissue Cd concentrations were significantly (p < 0.05) greater at the second harvest (3.56–4.67 mg kg−1 DM) relative to the first (1.45–1.53 mg kg−1 DM) and third (1.81–2.03 mg kg−1 DM) harvests. Notably, a large rainfall event (23 mm rainfall within 24 h) occurred during the second growth cycle that rewet the topsoil. In contrast, rainfall was infrequent and of low intensity during the first growth cycle, while only 1 mm of rainfall arrived during the entirety of the third growth cycle.
Practically, there is little ability to avoid grazing these crops following periods of increased soil moisture. Hence, where these Cd-sensitive crops are to be sown in Cd-enriched soil, management interventions should be explored to minimise potential increases in Cd phytoavailability that will be aligned with periods of increased soil moisture.
As soil pH is considered to be one of the single most important factors influencing Cd sorption–desorption behaviour and solubility [4,5,10,13,28], managing soil pH at the upper end of the optimum range for these crops would be a logical starting point. However, the effect of increasing soil pH on Cd accumulation in plantain and chicory has yet to be quantified in field trials. Furthermore, although the sorption of Cd is known to increase with increasing soil pH with a concurrent shift towards less labile, specifically-sorbed Cd [4,11,49], the assumption that this will reduce the size of the soil Cd fraction mobilised during periods of increased soil moisture needs to be proven. Further research in this area is warranted, as we do not yet fully understand the mechanisms behind the increases in Cd phytoavailability that were seen in our study that paralleled increases in soil moisture. Further to this, direct quantification of redox potential, soluble Mn/Fe concentrations, soil Cd fractionation, microbial respiration rate, and the influence of soil temperature is recommended for future research. This may help to explain the variation in soil extractable Cd concentrations observed in this research.

5. Conclusions

The results from this pot trial indicate that temporary fluctuations in soil moisture content can influence soil Cd phytoavailability, with increases in soil extractable Cd concentration paralleling increases in soil moisture content. However, soil saturation for short periods (3 days) did not appear to exert any influence on soil Cd phytoavailability, relative to irrigated but non-saturated soil conditions. The data from our trial suggest that plant tissue Cd concentrations are likely to be enriched in Cd-accumulator crops (such as chicory and plantain) during periods where soil moisture is abruptly increased following a period of drying. Our results also reinforce that soil pH and organic matter appear to be key factors regulating soil Cd phytoavailability, and need to be incorporated into Cd management frameworks where the objective is to manage risk of plant Cd accumulation.

Acknowledgments

We wish to acknowledge the Fertiliser Association of New Zealand (FANZ) for funding this research.

Author Contributions

A.D. Stafford and C.W.N. Anderson conceived and designed the experiments; A.D. Stafford performed the experiments; P.J. provided laboratory analytical support; A.D. Stafford analysed the data; A.D. Stafford wrote the paper with support from C.W.N. Anderson, P.J. and M.J. Hedley.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adriano, D.C. Trace Elements in the Terrestrial Environment; Springer: New York, NY, USA, 1986; 533p. [Google Scholar]
  2. Smolders, E.; Mertens, J. Cadmium. In Heavy Metals in Soils: Trace Metals and Metalloids in Soils and Their Bioavailability; Alloway, B.J., Ed.; Springer: Dordrecht, The Netherlands, 2012; pp. 283–311. [Google Scholar]
  3. Williams, C.H.; David, D.J. Some effects of the distribution of cadmium and phosphate in the root zone on the cadmium content of plants. Aust. J. Soil Res. 1977, 15, 59–68. [Google Scholar] [CrossRef]
  4. Loganathan, P.; Vigneswaran, S.; Kandasamy, J.; Naidu, R. Cadmium Sorption and Desorption in Soils: A Review. Crit. Rev. Environ. Sci. Technol. 2012, 42, 489–533. [Google Scholar] [CrossRef]
  5. Young, S.D. Chemistry of heavy metals and metalloids in soils. In Heavy Metals in Soils: Trace Metals and Metalloids in Soils and Their Bioavailability; Alloway, B.J., Ed.; Springer: Dordrecht, The Netherlands, 2012; pp. 51–95. [Google Scholar]
  6. McLaughlin, M.J.; Tiller, K.G.; Naidu, R.; Stevens, D.P. Review: The behaviour and environmental impact of contaminants in fertilisers. Aust. J. Soil Res. 1996, 34, 1–54. [Google Scholar] [CrossRef]
  7. Tiller, K.G.; Gerth, J.; Brummer, G. The relative affinities of Cd, Ni and Zn for different soil clay fractions and goethite. Geoderma 1984, 34, 17–35. [Google Scholar] [CrossRef]
  8. McBride, M.B. Reactions controlling heavy metal solubility in soils. Adv. Soil Sci. 1989, 10, 1–56. [Google Scholar]
  9. McBride, M.; Sauve, S.; Hendershot, W. Solubility control of Cu, Zn, Cd and Pb in contaminated soils. Eur. J. Soil Sci. 1997, 48, 337–346. [Google Scholar] [CrossRef]
  10. Christensen, T.H.; Haung, P.M. Solid phase cadmium and the reactions of aqueous cadmium with soil surfaces. In Cadmium in Soils and Plants; McLaughlin, M.J., Singh, B.R., Eds.; Kluwer Academic Publishers: Dordrecht, The Netherlands, 1999; pp. 65–96. [Google Scholar]
  11. Naidu, R.; Bolan, N.S.; Kookana, R.S.; Tiller, K.G. Ionic-strength and pH effects on the sorption of cadmium and the surface charge of soils. Eur. J. Soil Sci. 1994, 45, 419–429. [Google Scholar] [CrossRef]
  12. Naidu, R.; Kookana, R.S.; Sumner, M.E.; Harter, R.D.; Tiller, K.G. Cadmium sorption and transport in variable charge soils: A review. J. Environ. Qual. 1997, 26, 602–617. [Google Scholar] [CrossRef]
  13. Bolan, N.S.; Naidu, R.; Syers, J.K.; Tillman, R.W. Surface charge and solute interactions in soils. Adv. Agron. 1999, 67, 88–141. [Google Scholar]
  14. Iu, K.L.; Pulford, I.D.; Duncan, H.J. Influence of waterlogging and lime or organic matter additions on the distribution of trace metals in an acid soil—I. Manganese and iron. Plant Soil 1981, 59, 317–326. [Google Scholar] [CrossRef]
  15. Hernandez-Soriano, M.C.; Jimenez-Lopez, J.C. Effects of soil water content and organic matter addition on the speciation and bioavailability of heavy metals. Sci. Total Environ. 2012, 423, 55–61. [Google Scholar] [CrossRef] [PubMed]
  16. Charlatchka, R.; Cambier, P. Influence of reducing conditions on solubility of trace metals in contaminated soils. Water Air Soil Pollut. 2000, 118, 143–167. [Google Scholar] [CrossRef]
  17. Atta, S.K.; Mohammed, S.A.; Van Cleemput, O.; Zayed, A. Transformations of iron and manganese under controlled E(h), E(h)-pH conditions and addition of organic matter. Soil Technol. 1996, 9, 223–237. [Google Scholar] [CrossRef]
  18. Kashem, M.A.; Singh, B.R. Metal availability in contaminated soils: I. Effects of flooding and organic matter on changes in Eh, pH and solubility of Cd, Ni and Zn. Nutr. Cycl. Agroecosyst. 2001, 61, 247–255. [Google Scholar]
  19. Zheng, S.; Zheng, X.; Chen, C. Transformation of metal speciation in purple soil as affected by waterlogging. Int. J. Environ. Sci. Technol. 2013, 10, 351–358. [Google Scholar] [CrossRef]
  20. Zhang, C.; Ge, Y.; Yao, H.; Chen, X.; Hu, M. Iron oxidation-reduction and its impacts on cadmium bioavailability in paddy soils: A review. Front. Environ. Sci. Eng. 2012, 6, 509–517. [Google Scholar] [CrossRef]
  21. De Livera, J.; McLaughlin, M.J.; Hettiarachchi, G.M.; Kirby, J.K.; Beak, D.G. Cadmium solubility in paddy soils: Effects of soil oxidation, metal sulfides and competitive ions. Sci. Total Environ. 2011, 409, 1489–1497. [Google Scholar] [CrossRef] [PubMed]
  22. McLaren, R.G.; Cameron, K.C. Soil Science: Sustainable Production and Environmental Protection, 2nd ed.; Oxford University Press: Auckland, New Zealand, 1996; p. 304. [Google Scholar]
  23. Rumball, P.J. Pasture Plant, Water-Table Interactions on Manawatu Sandy Country. In Proceedings of the New Zealand Grassland Association 39, Palmerston North, New Zealand, 1977; Available online: https://www.grassland.org.nz/publications/nzgrassland_publication_1394.pdf (accessed on 9 February 2018).
  24. Hindersmann, I.; Mansfeldt, T. Trace element solubility in a multimetal-contaminated soil as affected by redox conditions. Water Air Soil Pollut. 2014, 225, 1–20. [Google Scholar] [CrossRef]
  25. Cornu, S.; Cattle, J.A.; Samouelian, A.; Laveuf, C.; Guilherme, L.R.G.; Alberic, P. Impact of redox cycles on manganese, iron, cobalt, and lead in nodules. Soil Sci. Soc. Am. J. 2009, 73, 1231–1241. [Google Scholar] [CrossRef] [Green Version]
  26. Backes, C.A.; McLaren, R.G.; Rate, A.W.; Swift, R.S. Kinetics of cadmium and cobalt desorption from iron and manganese oxides. Soil Sci. Soc. Am. J. 1995, 59, 778–785. [Google Scholar] [CrossRef]
  27. Kim, N.D.; Fergusson, J.E. Adsorption of cadmium by an aquent New Zealand soil and its components. Aust. J. Soil Res. 1992, 30, 159–167. [Google Scholar] [CrossRef]
  28. Gray, C.W.; McLaren, R.G.; Roberts, A.H.C.; Condron, L.M. Solubility, sorption and desorption of native and added cadmium in relation to properties of soils in New Zealand. Eur. J. Soil Sci. 1999, 50, 127–137. [Google Scholar] [CrossRef]
  29. Ponnamperuma, F.N. The chemistry of submerged soils. In Advances in Agronomy; Academic Press: Cambridge, MA, USA, 1972; Volume 24, pp. 29–96. [Google Scholar]
  30. Bhatti, S.M.; Anderson, C.W.N.; Stewart, R.B.; Robinson, B.H. Risk assessment of vegetables irrigated with arsenic-contaminated water. Environ. Sci. Process. Impacts 2013, 15, 1866–1875. [Google Scholar] [CrossRef] [PubMed]
  31. Stafford, A.D.; Anderson, C.W.N.; Hedley, M.J.; McDowell, R.W. Cadmium accumulation by forage species used in New Zealand livestock grazing systems. Geoderma Reg. 2016, 7, 11–18. [Google Scholar] [CrossRef]
  32. Powell, A.M.; Kemp, P.D.; Jaya, I.K.D.; Osborne, M.A. Establishment, Growth and Development of Plantain and Chicory under Grazing. In Proceedings of the New Zealand Grassland Association 69, Wairakei, New Zealand, 2007; Available online: https://www.grassland.org.nz/publications/nzgrassland_publication_144.pdf (accessed on 9 December 2017).
  33. Li, Z.; McLaren, R.G.; Metherell, A.K. Cobalt and manganese relationships in New Zealand soils. N. Z. J. Agric. Res. 2001, 44, 191–200. [Google Scholar] [CrossRef]
  34. Jeyakumar, P.; Loganathan, P.; Anderson, C.N.; Sivakumaran, S.; McLaren, R. Comparative tolerance of Pinus radiata and microbial activity to copper and zinc in a soil treated with metal-amended biosolids. Environ. Sci. Pollut. Res. 2014, 21, 3254–3263. [Google Scholar] [CrossRef] [PubMed]
  35. Tukey, J.W. Exploratory Data Analysis; Addison-Wesley: Reading, PA, USA, 1977; 688p. [Google Scholar]
  36. Gotoh, S.; Patrick, W.H. Transformations of manganese in a waterlogged soil as affected by redox potential and pH. Soil Sci. Soc. Am. J. 1972, 36, 738–742. [Google Scholar] [CrossRef]
  37. DeLaune, R.D.; Reddy, K.R. Redox Potential. In Encyclopedia of Soils in the Environment; Hillel, D., Ed.; Academic Press: Amsterdam, The Netherlands, 2005; pp. 366–371. [Google Scholar]
  38. Narteh, L.T.; Sahrawat, K.L. Influence of flooding on electrochemical and chemical properties of West African soils. Geoderma 1999, 87, 179–207. [Google Scholar] [CrossRef]
  39. Gotoh, S.; Patrick, W.H. Transformations of iron in a waterlogged soil as influenced by redox potential and pH. Soil Sci. Soc. Am. J. 1974, 38, 66–71. [Google Scholar] [CrossRef]
  40. Cornu, J.Y.; Denaix, L.; Schneider, A.; Pellerin, S. Temporal evolution of redox processes and free Cd dynamics in a metal-contaminated soil after rewetting. Chemosphere 2007, 70, 306–314. [Google Scholar] [CrossRef] [PubMed]
  41. Gray, C.W.; McLaren, R.G.; Roberts, A.H.C.; Condron, L.M. Fractionation of soil cadmium from some New Zealand soils. Commun. Soil Sci. Plant Anal. 2000, 31, 1261–1273. [Google Scholar] [CrossRef]
  42. Yamane, I.; Sato, K. Effect of temperature on the decomposition of organic substances in flooded soil. Soil Sci. Plant Nutr. 1967, 13, 94–100. [Google Scholar] [CrossRef]
  43. Lovley, D.R. Microbial reduction of iron, manganese, and other metals. Adv. Agron. 1995, 54, 175–231. [Google Scholar]
  44. Loganathan, P.; Hedley, M.J.; Gregg, P.E.H.; Currie, L.D. Effect of phosphate fertiliser type on the accumulation and plant availability of cadmium in grassland soils. Nutr. Cycl. Agroecosyst. 1997, 47, 169–178. [Google Scholar] [CrossRef]
  45. Mortvedt, J.J.; Mays, D.A.; Osborn, G. Uptake by wheat of cadmium and other heavy metal contaminants in phosphate fertilizers. J. Environ. Qual. 1981, 10, 193–197. [Google Scholar] [CrossRef]
  46. Roberts, A.H.C.; Longhurst, R.D. Cadmium cycling in sheep-grazed hill-country pastures. N. Z. J. Agric. Res. 2002, 45, 103–112. [Google Scholar] [CrossRef]
  47. Tan, K.H. Soil Sampling, Preparation, and Analysis, 2nd ed.; Taylor & Francis: Boca Raton, FL, USA, 2005; 672p. [Google Scholar]
  48. Somasiri, S.C.; Kenyon, P.R.; Kemp, P.D.; Morel, P.C.H.; Morris, S.T. Growth performance and carcass characteristics of lambs grazing forage mixes inclusive of plantain (Plantago lanceolata L.) and chicory (Cichorium intybus L.). Small Rumin. Res. 2015, 127, 20–27. [Google Scholar] [CrossRef]
  49. Bolan, N.S.; Adriano, D.C.; Mani, P.A.; Duraisamy, A. Immobilization and phytoavailability of cadmium in variable charge soils. II. Effect of lime addition. Plant Soil 2003, 251, 187–198. [Google Scholar]
Figure 1. Design of pots (left, completed; right, partially complete showing the slit in the internal drainage tube) used for generating cyclical 3-day flooded, 11-day drained soil conditions.
Figure 1. Design of pots (left, completed; right, partially complete showing the slit in the internal drainage tube) used for generating cyclical 3-day flooded, 11-day drained soil conditions.
Soilsystems 02 00009 g001
Figure 2. Daily mean, maximum, and minimum temperatures recorded at Tauranga airport (NIWA “CliFlo” climate station 1615) over the duration of the trial period.
Figure 2. Daily mean, maximum, and minimum temperatures recorded at Tauranga airport (NIWA “CliFlo” climate station 1615) over the duration of the trial period.
Soilsystems 02 00009 g002
Figure 3. Mean soil extractable Cd concentrations (µg kg−1) at each sampling date under flooded and non-flooded conditions for (a) Topehaehae Planted; (b) Topehaehae Unplanted; (c) Kereone Planted; and (d) Kereone Unplanted treatments. ANOVA p-values are provided for comparison of means at each sampling date. Error bars represent the 95% confidence interval.
Figure 3. Mean soil extractable Cd concentrations (µg kg−1) at each sampling date under flooded and non-flooded conditions for (a) Topehaehae Planted; (b) Topehaehae Unplanted; (c) Kereone Planted; and (d) Kereone Unplanted treatments. ANOVA p-values are provided for comparison of means at each sampling date. Error bars represent the 95% confidence interval.
Soilsystems 02 00009 g003
Table 1. Soil classification and chemical characteristics of the two soils used in this study.
Table 1. Soil classification and chemical characteristics of the two soils used in this study.
Soil TypeKereone Fine Sandy LoamTopehaehae Sandy Clay Loam
New Zealand Soil ClassificationTypic Orthic AllophanicTypic Recent Gley
USDA TaxononyUdandAquept
Total Cd (mg kg−1 DW)0.790.61
Total P (mg kg−1 DW)15771132
Total Mn (mg kg−1 DW)1840710
Total Fe (%)2.12.4
Total C (%)5.73.2
pH6.35.6
Cation Exchange Capacity (meq 100 g−1)2617
P retention (%)6934
Table 2. Mean soil extractable Cd concentration (µg kg−1) across the 18 paired sampling dates; (i) Overall mean; (ii) Mean of pre-irrigation samplings; and (iii) Mean of +3 day samplings; in addition to (iv) Mean plant tissue Cd concentration (mg kg−1 DM); (v) Plant yield (g DM pot−1); and (vi) Plant Cd uptake (µg Cd pot−1). Within a row, treatment means with the same letter indicate that differences are not significant from one another (p < 0.05).
Table 2. Mean soil extractable Cd concentration (µg kg−1) across the 18 paired sampling dates; (i) Overall mean; (ii) Mean of pre-irrigation samplings; and (iii) Mean of +3 day samplings; in addition to (iv) Mean plant tissue Cd concentration (mg kg−1 DM); (v) Plant yield (g DM pot−1); and (vi) Plant Cd uptake (µg Cd pot−1). Within a row, treatment means with the same letter indicate that differences are not significant from one another (p < 0.05).
MeasureTreatmentp-Value
KPNKUNKPFKUFTPNTUNTPFTUF
Soil extractable Cd (µg kg−1)(i) Overall mean20.9 (a)21.1 (a)20.0 (a)18.5 (a)84.4 (b)89.3 (b)88.9 (b)86.1 (b)***
Std. Error0.70.70.80.83.43.43.53.4
n8887888486888784
(ii) Pre-irrigation mean19.3 (a)19.4 (a)17.9 (a)16.1 (a)83.0 (b)87.6 (b)85.4 (b)81.4 (b)***
Std. Error1.01.11.01.15.45.24.95.0
n4443454144444440
(iii) +3 day mean22.6 (a)22.7 (a)22.1 (a)20.7 (a)86.0 (b)91.0 (b)92.5 (b)90.4 (b)***
Std. Error1.00.91.21.14.14.34.94.6
n4444434342444344
p-value (pre-irrigation vs. +3 day)******n.s.n.s.n.s.n.s.
(iv) Plant tissue Cd (mg kg−1 DM)Mean0.208 (a)0.384 (ab)0.650 (c)0.604 (bc)***
Std. Error0.0590.0770.0540.068
n5555
Main effect—Soil type***
Main effect—Irrigation regimen.s.
(v) Plant yield (g DM pot−1)Mean5.40 (c)2.70 (ab)4.02 (bc)1.51 (a)***
Std. Error0.570.290.300.16
Main effect—Soil typen.s.
Main effect—Irrigation regime***
(vi) Plant Cd uptake (µg Cd pot−1)Mean1.001 (a)1.020 (a)2.602 (b)0.880 (a)***
Std. Error0.2240.1820.2800.056
Main effect—Soil type*
Main effect—Irrigation regime*
* p ≤ 0.05, ** p ≤ 0.01, *** p ≤ 0.001, n.s. = not significant. Treatment codes: K = Kereone, T = Topehaehae; P = Planted, U = Unplanted; N = Non-flooded, F = Flooded.
Table 3. pH of 0.05 M CaCl2 extracts assessed periodically throughout the trial. Means with the same letter indicate that differences from one another are not significant (p < 0.05). Lower-case letters represent comparisons between irrigation treatment means (n = 5) and upper-case letters represent comparisons between the overall mean pH at different time intervals (n = 20).
Table 3. pH of 0.05 M CaCl2 extracts assessed periodically throughout the trial. Means with the same letter indicate that differences from one another are not significant (p < 0.05). Lower-case letters represent comparisons between irrigation treatment means (n = 5) and upper-case letters represent comparisons between the overall mean pH at different time intervals (n = 20).
Irrigation TreatmentDay of Samplingp-Value
02884112
KPN4.044.09 (ab)4.294.44
KPF3.944.28 (b)4.384.27
TPN3.974.00 (a)4.174.23
TPF3.924.10 (ab)4.134.12
p-value0.1500.0420.2940.144
Sampling mean3.97 (A)4.11 (AB)4.24 (B)4.27 (B)<0.001

Share and Cite

MDPI and ACS Style

Stafford, A.; Jeyakumar, P.; Hedley, M.; Anderson, C. Influence of Soil Moisture Status on Soil Cadmium Phytoavailability and Accumulation in Plantain (Plantago lanceolata). Soil Syst. 2018, 2, 9. https://doi.org/10.3390/soils2010009

AMA Style

Stafford A, Jeyakumar P, Hedley M, Anderson C. Influence of Soil Moisture Status on Soil Cadmium Phytoavailability and Accumulation in Plantain (Plantago lanceolata). Soil Systems. 2018; 2(1):9. https://doi.org/10.3390/soils2010009

Chicago/Turabian Style

Stafford, Aaron, Paramsothy Jeyakumar, Michael Hedley, and Christopher Anderson. 2018. "Influence of Soil Moisture Status on Soil Cadmium Phytoavailability and Accumulation in Plantain (Plantago lanceolata)" Soil Systems 2, no. 1: 9. https://doi.org/10.3390/soils2010009

APA Style

Stafford, A., Jeyakumar, P., Hedley, M., & Anderson, C. (2018). Influence of Soil Moisture Status on Soil Cadmium Phytoavailability and Accumulation in Plantain (Plantago lanceolata). Soil Systems, 2(1), 9. https://doi.org/10.3390/soils2010009

Article Metrics

Back to TopTop