key: cord-1003623-tun72v1u authors: Di Lodovico, Eliana; Marchand, Lilian; Oustrière, Nadège; Burges, Aritz; Capdeville, Gaelle; Burlett, Régis; Delzon, Sylvain; Isaure, Marie-Pierre; Marmiroli, Marta; Mench, Michel J. title: Potential ability of tobacco (Nicotiana tabacum L.) to phytomanage an urban brownfield soil date: 2021-10-18 journal: Environ Sci Pollut Res Int DOI: 10.1007/s11356-021-16411-y sha: 150ae332e0c7289f97e72b38768c311d1251110d doc_id: 1003623 cord_uid: tun72v1u The ability of tobacco (Nicotiana tabacum L. cv. Badischer Geudertheimer) for phytomanaging and remediating soil ecological functions at a contaminated site was assessed with a potted soil series made by fading an uncontaminated sandy soil with a contaminated sandy soil from the Borifer brownfield site, Bordeaux, SW France, at the 0%, 25%, 50%, 75%, and 100% addition rates. Activities of sandblasting and painting with metal-based paints occurred for decades at this urban brownfield, polluting the soil with metal(loid)s and organic contaminants, e.g., polycyclic aromatic hydrocarbons, in addition to past backfilling. Total topsoil metal(loid)s (e.g., 54,700 mg Zn and 5060 mg Cu kg(−1)) exceeded by seven- to tenfold the background values for French sandy soils, but the soil pH was 7.9, and overall, the 1M NH(4)NO(3) extractable soil fractions of metals were relatively low. Leaf area, water content of shoots, and total chlorophyll (Chl) progressively decreased with the soil contamination, but the Chl fluorescence remained constant near its optimum value. Foliar Cu and Zn concentrations varied from 17.8 ± 4.2 (0%) to 27 ± 5 mg Cu kg(−1) (100%) and from 60 ± 15 (0%) to 454 ± 53 mg Zn kg(−1) (100%), respectively. Foliar Cd concentration peaked up to 1.74 ± 0.09 mg Cd kg(−1), and its bioconcentration factor had the highest value (0.2) among those of the metal(loid)s. Few nutrient concentrations in the aboveground plant parts decreased with the soil contamination, e.g., foliar P concentration from 5972 ± 1026 (0%) to 2861 ± 334 mg kg(−1) (100%). Vulnerability to drought-induced embolism (P50) did not differ for the tobacco stems across the soil series, whereas their hydraulic efficiency (Ks) declined significantly with increasing soil contamination. Overall, this tobacco cultivar grew relatively well even in the Borifer soil (100%), keeping its photosynthetic system healthy under stress, and contaminant exposure did not increase the vulnerability of the vascular system to drought. This tobacco had a relevant potential to annually phytoextract a part of the bioavailable soil Zn and Cd, i.e., shoot removals representing here 8.8% for Zn and 43.3% for Cd of their 1M NH(4)NO(3) extractable amount in the potted Borifer soil. [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11356-021-16411-y. Despite some progress made in the management of contaminated sites by EU Member States and cooperating countries, the use of phytotechnologies to remediate contaminated brownfield soils is still not widespread (Payá- Pérez and Peláez-Sánchez 2017) . Increasing the frequency of their use remains a huge task. The US EPA listed 9 million ha of possibly contaminated land and 1327 abandoned, worst hazardous waste sites on its National Priority List (USEPA 2021). Chinese contaminated soils would reach 20 million ha (Sun et al. 2019) . Urban brownfields represent at least 150,000 ha in France. Between 2010 and 2016, 95 sites have been remediated (65 completed, 30 in progress) by the French agency of ecological transition, representing 491 ha (ADEME et al. 2018) . The lack of specific policy instruments and guidelines on the reuse of contaminated soils can impede sustainable remediation even in countries with well-developed policies and regulatory frameworks (Reinikainen et al. 2016; Moreira et al. 2021) . In parallel, prospects and threats regarding future land use, e.g., options of site redevelopment and the potential long-term liabilities, in line with economic and social uncertainties, influence the strategy for risk management and the desired remediation level, the trend being to do more than required after risk assessment (Bardos et al. 2020; Moreira et al. 2021) . Progress toward sustainable land management of contaminated soils is driven by several generic elements: e.g., increased recognition of environmental impacts and benefits of remediation, stakeholders' demand for more sustainable practices, lack of regulatory requirements, economic considerations, and institutional pressure that promotes sustainable practices (Hou et al. 2014; Rizzo et al. 2016 ; Bardos et al. 2018 Bardos et al. , 2020 . The European Parliament has just adopted a resolution in this way for the soil protection and the remediation of contaminated soils (European Parliament 2021). Two main priorities for the topic of phytoremediation and phytomanagement of polluted soils are (1) to guide the collection of case study data for identifying and assessing methods and (2) to develop tools to guide remediation of soil ecological functions underlying ecosystem services (Macci et al. 2020) . Indeed, each brownfield site has specific environmental challenges and management plan. With this in mind, the Borifer brownfield site is a part of the Parc aux Angéliques located in the Bordeaux downtown, south-west France, with contaminated soil requiring remediation. The topsoil is mostly barren at the Borifer site, due to anthropogenic soil contamination by metal(loid)s and organic xenobiotics, soil compaction, low fertility, and low water-holding capacity (Marchand and Mench 2015) . Some vegetation patches are colonizing this brownfield. Thirty-one plant species were identified in May 2018, Bromus sterilis L. and Vulpia myuros (L.) C.C.Gmel. being dominant. Most of plant colonists had a low shoot biomass, and at several patches, plants displayed some visible symptoms of phytotoxicity (e.g., leaf chlorosis, dwarfism, and sometimes purple coloration on aboveground plant parts). Afte r so il l oosen ing , a t emporary phytomanagement could cover the soil surface, prevent dust from flying, stimulate the xenobiotic biodegradation, mitigate metal(loid) migration with natural agents, and improve soil properties. This would correspond with a future land use to complete a green corridor on the right bank of the Garonne River, crossing the Bordeaux downtown. High total soil Cu and Zn and many organic contaminants were evidenced in the Borifer topsoils (Table 1 , Marchand and Mench 2015) . This can affect the microbial community decreasing its density, diversity, and vitality, especially due to the toxicity of the Cu ionic form (Cu 2+ ) (Bourceret et al. 2018; Xue et al. 2018) . Soil Zn-and Cu-contamination can affect plant growth as well (Chaney 1993; Verdejo et al. 2016; Küpper and Andresen 2016; Kolbas et al. 2018) . Relevant tools to remediate and contain soil contamination include the Gentle Remediation Options (GRO), which aim at reducing the pollutant linkages through transformation, stabilization, degradation, or extraction of contaminants using fungi (myco-), plant (phyto-), and/or bacteria-based methods, with or without amendments or chemical additives (Cundy et al., 2016) . Cultivation of fibrous, aromatic, and other non-food crops is an option, e.g., flax (Linum usitatissimum L.), jute (Corchorus capsularis L.), hemp (Cannabis sativa L.), ramie (Boehmeria nivea (L.) Gaudich., kenaf (Hibiscus cannabinus L), nettles (Urtica dioïca L.), (Herzig et al. 2014; Kidd et al. 2015; Kolbas et al. 2018 Kolbas et al. , 2020 Mench et al. 2017 Mench et al. , 2018 Ndubueze 2018) . Tobacco is a plant species frequently evaluated in Western Europe and also in China for phytomanaging contaminated soils because it can extract and store metals in its tissues: Zn and Cd in the shoots and Cu and Pb mainly in the roots (Herzig et al. 2014; Thijs et al. 2018; Yang et al. 2017; Rehman et al. 2019; Moreira et al. 2021) . It can stimulate the biodegradation of organic xenobiotic by the rhizosphere microorganisms (Rezek et al. 2004; Azaizeh et al. 2011; Prouzova et al. 2012; Bisht et al. 2015; Petrová et al., 2017) . Processes that valorize tobacco biomass were reported by Asad et al. (2017) , Yang et al. (2017) , Rehman et al. 2019, and Kolbas et al. (2020) , e.g., pyrolysis, hydrothermal oxidation, fermentation, and gasification. The main objectives to be reached by the soil phytomanagement using tobacco at the Borifer site would be to progressively remove the phytoavailable soil Cd and Zn by harvesting shoots, to immobilize in the roots the phytoavailable fraction of several metal(loid)s, e.g., Cu and Pb, to stimulate the dissipation of organic pollutants through the rhizodeposition and microbial biodegradation, and to improve soil properties, e.g., sequestration of carbon. Consequently, this study aimed at assessing the ability of tobacco to grow on a potted soil series with an increasing proportion of the contaminated soil of the Borifer site, to evaluate its behavior toward metal(loid)s in this sandy soil with mixed contamination, and to assess its vulnerability to drought-induced embolism in such soil conditions as periods of intense spring and summer droughts become more frequent with climate change. In addition to the characterization of soil properties, fresh weight (FW) and dried weight (DW) biomass of aboveground plant parts, leaf area (LA), vulnerability to drought-induced embolism (P50) and hydraulic efficiency (Ks) of tobacco stems, the foliar and stem ionomes (Salt et al. 2008) , pigment content, and chlorophyll fluorescence of the leaves were determined. The former Borifer site (44°51′ 08.6′′ N, 0°33′ 29.5′′ W, altitude 15 m) is located in the Bordeaux downtown (Gironde County, South-West France), on the right bank of the Garonne River, near the Parc aux Angéliques and Chaban-Delmas Bridge ( Fig. 1 ). At this site (1.1 ha), activities of sandblasting with abrasive glass and iron grits and painting with metal-based paints of detached metallic pieces and naval parts took place for decades (Marchand and Mench 2015) . A historical soil contamination is also due to the backfilling of the Garonne River banks. The bedrock is constituted by the Mattes clays, these bluish to grayish clays with peaty past being present on both sides of the Garonne River. Since the City of Bordeaux decided to include this site in the urban park, the soil analysis was performed showing high concentrations of metal(loid)s and xenobiotics (PAH, PCB, aliphatic hydrocarbons, etc.) (Marchand and Mench 2015 , Table S1 ). The soil texture is sandy and the metal(loid) concentrations in the topsoil largely exceed (seven-to tenfold) the background values for French sandy soil, depending on the plots. The Borifer soil (B) was collected (100 kg DW) in the 0-30-cm depth layer of the S41 plot. It is mainly sandy backfill with pebbles. The uncontaminated control soil (Ctrl) was sampled in a kitchen garden, Gradignan, France. Both sandy soils were air-dried and sieved at 4 mm. Soil texture and physico-chemical parameters were determined with standard methods and a quality scheme by INRA LAS (2018), Arras, France, i.e., inductively coupled plasma/atomic emission spectroscopy (ICP-AES) for metals after wet digestion (HF and HCIO 4 ) and hydride-generation for As after wet digestion in H 2 SO 4 /HNO 3 (2/1) with V 2 O 5 at 100°C (3 h). Soil particles were separated into a light (44% of the soil mass) and heavy (56% of the soil mass) fractions by densimetry using a sodium polytungstate liquor with a density of 2.9 (Morgun and Makarov, 2011) . Samples were ground and sieved at 100 μm. In the bulk soil, a fraction was non-crushable and was composed of shiny metallic and black particles ( Fig. 2A) . These noncrushable particles were separated into a magnetic fraction including ferromagnetic/ferrimagnetic phases and a non-magnetic fraction using a hang magnet (Fig. 2B , C). The bulk soil and fraction mineralogy was investigated by X-ray diffraction (XRD) using a D2 diffractometer (Brucker) operating with a Cu Ka source. X-ray powder diffraction patterns were recorded from 5 to 90°2 θ with a 0.02°step and a 5-s counting time per point. Soil properties are given in Tables 1 and 2. The Borifer soil was faded with the Ctrl soil, at rates varying from 0% (Ctrl) to 25% (B25), 50% (B50), 75% (B75), and 100% (B), and filled (1000 g air-dried soil) in plastic pots, with four replicates. Potted soils were irrigated with tap water, and let to react for 2 weeks (60% soil humidity, laboratory conditions). The tobacco plants cv. Badischer Geudertheimer (BaG) were firstly sowed and grown in the Ctrl soil under a greenhouse for 4 months. One plant was thereafter transplanted in each pot containing various proportions of Borifer soil on July 8. Pots were randomly placed in a laboratory with glass wall (due to the Covid 19 lockdown and to prevent heatwaves, the greenhouses being not daily accessible), their position being changed every 2-3 days. From the first-day, maximum shoot length (height), the number of leaves, and chlorophyll content (chlorophyll meter "Opti-sciences CCM-200") were determined. The measurements were taken in days 1, 23, 41, 51, and 64 for following the development of the leaf stages. On day 23, some mineral NPK fertilizers were added (Blaukorn classic, 12-8-16 (3-25)) at 40 kg N (58% NH 4 , 42% NO 3 ), 26.7 kg P 2 O 5 , 53 kg K 2 O, 10 kg MgO, 83 kg SO 3 , 0.06 kg B, 0.2 kg Fe, and 0.04 kg Zn per ha. The plants were daily watered (up to 70% of the water holding capacity) considering few aspects: the soil water capacity (10% of the soil total weight), tobacco evapotranspiration, and the increase of 5-g fresh weight (FW)/plant/week due to the gain in biomass. Between the 3rd and the 4th weeks from the transplantation day, the plants started to show the effects of the contaminant exposure with some leaves displaying interveinal discoloration (plants grown in the B soil). On the days 57 and 62, the plants were placed in a dark room overnight, and the chlorophyll fluorescence was measured (Chlorophyll fluorometer PAM-210 Walz, Germany) in the morning before the sunrise. The minimal fluorescence level in the dark-adapted state ( F 0 ) w a s m e a s u r e d u s i n g a m o d u l a t e d p u l s e (<0.05 μmol m −2 s −1 for 1.8 μs). Maximal fluorescence in this state (Fm) was measured after applying a saturating actinic light pulse of 15,000 μmol m −2 s −1 for 0.7 s. The value of Fm was recorded as the highest average of four consecutive points. Values of variable fluorescence (Fv = Fm−F 0 ) and maximum quantum efficiency of PSII photochemistry (Fv/Fm) were calculated from F 0 and Fm. This Fv/Fm ratio correlates with the number of functional PSII reaction centers, and darkadapted values of Fv/Fm can be used to quantify photoinhibition (Cambrollé et al. 2012 ). At sampling time (day 66), all leaves and stems were weighed for determining their total fresh weight (FW). After measuring the plant traits, plant parts were washed in MilliQ water and put in a ventilated oven (50°C, till constant weight) to determine the total dried weight (DW). Fresh stem samples were kept to determine their vulnerability to drought-induced embolism and hydraulic conductivity (see below). Plant samples were labeled according to the soil treatments and the addition rate of B soil (i.e., Ctrl, B25, B50, B75, and B). Four discs (0.8-cm diameter, giving 1 cm 2 in duplicate) were taken at the mid-length of the 6th leaf (starting from the top) of each plant, from both mid-rib sides. Two discs were placed in 3 mL cold DMF (N,N-dimethylformamide) (in duplicate). The extraction lasted 48h at 4°C, then the content of chlorophyll A, B, and carotenoids was computed after measuring the absorbance at 470, 664, and 647 nm with a spectrophotometer (Libra S22 of Biochrom). Formulae and extinction coefficients for computing chlorophyllous pigments were published by Lichtenhaler and Wellburn (1983) , and Blanke (1990) . All the leaves were cut from the plants and placed between filter papers. Each leaf was labeled with the treatment, the replicate number, and the leaf number. The leaves nos. 4, 8, and 9 of each plant were used to determine the leaf surface using a scanner (Epson Expression 10000 XL) and the WinFOLIA software. The plant stems were used to characterize the vulnerability to drought-induced embolism at the Caviplace laboratory (PHENOBOIS platform, INRAE-University of Bordeaux, France) with the Cavitron technique (Cochard et al. 2005) . Centrifugal force was used to establish negative pressure in the xylem and to provoke water stress-induced cavitation, using a custom-built honeycomb rotor (DGMeca, Gradignan, France) mounted on a high-speed centrifuge (Sorvall RC5 plus, MSE Scientific, London, UK). This technique enables to measure both the hydraulic conductance of a branch under negative pressure and the vulnerability of stem xylem to water stress-induced embolism caused by air seeding (Delzon et al. 2010) . Each stem was cut at the base, the foliage removed, then wrapped in a wet cloth and placed in a plastic bag to maintain it hydrated with a high moisture level. The stems were then stored in a fridge at 5°C before carrying on the vulnerability to embolism measurements. Prior to measurement, stems were re-cut under water to a standard length of 27 cm. Samples were infiltrated with a reference ionic solution of 10 mM KCl and 1 mM CaCl 2 in deionized ultrapure water. Centrifugal force was used to generate negative pressure into the xylem and induce cavitation. Maximum conductance of stem (K max in m 2 MPa −1 s −1 ) was calculated under low xylem pressures (close to zero) (P in MPa). Then, rotation speed of the centrifuge was gradually increased by 0.5 or 1 MPa, to lower xylem pressure. The percentage loss of hydraulic conductance (PLC) of the stem was determined at each pressure step to build a vulnerability curve corresponding to PLC versus xylem pressure (P) (Delzon et al. 2010 ). For each sample, a sigmoid function (Pammenter and Willigen, 1998) was fitted to the vulnerability curve using proc NLIN in SAS 9.4 according to the equation: where P 50 (MPa) is the xylem pressure inducing 50% loss of hydraulic conductance and S (% MPa −1 ) is the slope of the vulnerability curve at the inflection point (Urli et al. 2013) . The xylem pressures inducing 12% (P 12 , xylem air entry point; MPa) and 88% (P 88 ; MPa) loss of hydraulic conductivity were calculated as follows: P 12 = 50/S + P 50 and P 88 = −50/S + P 50 (Urli et al. 2013) . The specific hydraulic conductivity (Ks, m 2 MPa −1 s −1 ) of the xylem was calculated by dividing the maximum hydraulic conductivity measured at low speed by the sapwood area of the sample. These parameters were averaged for tobacco plants of each soil treatment. After the measurements, stem samples were oven-dried and weighed. Dried leaf and stem samples were ground (<1.0-mm particle size, Fritsch Pulverisette 19). Weighed aliquots (0.5 g DW) were put in TFM® (modified poly(tetrafluoroethylene) tubes with 5 mL suprapure 14 M HNO 3 and left overnight. Then, 2 mL 30% (v/ v) H 2 O 2 not stabilized by phosphates and 1 mL MilliQ water were added in each tube and samples wet-digested under microwaves (CEM Marsxpress 1200 W: from 0 to 10 min till 120°C, 11 to 15 min to reach 180°C, 15 min at 180°C, cooling during 1 h). Digested samples were filtered through ash-free papers (Dutscher 1440110B, porosity 5-8 μm) and adjusted to 100 mL with milliQ water. Certified reference material (BIPEA maize V463) and blank reagents were included in all series. Mineral composition (As, Ca, Cu, Cr, Fe, K, Mg, Mn, Na, Ni, P, and Zn) in digests was determined by ICP-MS (Thermo X series 200, INRAE USRAVE laboratory, Villenave d'Ornon, France). All elements were recovered (>95%) according to the standard values, and standard deviation for replicates was <5%. The bioconcentration factor (BCF) was computed as foliar element concentration vs. total soil element. The shoot biomass, ionome, and all other plant parameters were tested using the one-way ANOVA (analysis of variance), and when significant differences occurred between the treatments, the post-hoc Tukey test was used to make multiple comparisons of mean values. Besides, a redundancy analysis (RDA) was performed, to make a multivariate comparison among all the parameters considered. Statistical analyses were performed using the R Project for Statistical Computing package 3.0.3 (Foundation for Statistical computing, Vienna, Austria). For the specific hydraulic conductivity, the Student-Newman-Keuls test and the REG procedure were used (SAS/STAT® 13.1). Soil properties (Tables 1 and 2) The B soil displayed high metal(loid) concentrations, especially for Zn, Cu, Pb, Cr, and Ni, which largely exceeded their background concentrations in the French sandy soils and the French screening values (Table 1, Baize 1997) . Only 0.3% of total soil Cu and 0.05% of total soil Zn were present in the 1M NH 4 NO 3extractable fraction of the B soil, likely due to metal speciation and alkaline soil pH. The extractable soil Zn exceeded the indicative value (0.5 mg Zn kg −1 that guarantees full soil fertility) of the Swiss soil protection act (Herzig et al. 2014) . Extractable soil Pb and Cd were relatively low in the B soil. The Ctrl soil was not contaminated. Both soils had the same pH (7.9). The bulk Borifer soil was mainly composed of quartz (SiO 2 ), calcite (CaCO 3 ), microcline (KAlSi 3 O 8 ), fayalite (Fe 2 SiO 4 ), wüstite (FeO), and potentially pyrite (FeS 2 ) (Fig. 3a) . Densimetric and magnetic fractionation allowed to identify other crystallized minerals. In addition to fayalite and wüstite, magnetite (Fe 3 O 4 )/maghemite (Fe 2 O 3 ) and zinc and copper chromite (ZnCr 2 O 4 and CuCr 2 O 4 , respectively) were detected in the heavy fraction (Fig. 3b) . The non-crushable magnetic fraction was composed of metallic iron and magnetite/maghemite present as black and reddish particles (Figs. 2B and 3c) while the non-magnetic fraction displayed FeCu 4 likely corresponding to the shiny golden particles observed in Fig. 2C . The identified metal-containing phases, i.e., metallic iron, iron silicates, and oxides as well as Zn, Cu chromite, and FeCu 4 likely resulted from materials used at the brownfield site for abrasion and painting. They suggested a low metal phytoavailability. Organic C and total N contents were lower in the B soil than in the Ctrl soil, but similar to the values of the French Fluviosol soil series (Table 1 ; Richer de Forges 2020). The organic matter (OM) content and C/N value were high as compared to the French sandy soils (Baize 1997) . The B soil was contaminated by PAH, fluoranthene, and then pyrene (to a lesser extent) being the major compounds (Table 2 ). For 8 out of 16 PAH chemicals, their concentration in the B soil exceeded the environmental health guideline values for agricultural land use regarding either soil quality guideline or soil and food ingestion based on non-carcinogenic effects (Table 2, Environment Canada 2010). However none exceeded the environmental health guideline value for a residential/parkland use. High soil OM in the B soil may restrict the PAH bioavailability (Harvey et al. 2002) . The shoot biomass (SDW), the maximum stem length (H), and the leaf area (LA) followed the same trend. They slightly decreased along with the increased proportion of contaminated B soil (Table 3) , albeit this drop was only significant for the leaf area (−31% at 100% of B soil). The number of leaves did not significantly differ across the soil series. However some leaves displayed interveinal discoloration on their adaxial side for plants at the highest B soil rates. The water content of shoots significantly fell at the 75% and 100% addition rates of B soils (−28% and −40%, respectively). Despite the high contamination of the B soil, the impacts on the morphological parameters of tobacco plants were quite low. The plant exposure to PAHs can occur via particle-phase deposition on the waxy leaf cuticle and in the stomata, by the uptake from the gas phase through stomata, and the root uptake from the soil solution and then the liquid phase transfer in the transpiration streams (Gworek et al. 2016 ). In our pot experiment, the atmosphere was not contaminated by dust from the B soil, and continuously renewed, likely avoiding foliar exposure. The total chlorophyll density (Chl) significantly decreased (−26%) for the plants grown in the B soil, as well as the Chl A (−30%) and Chl B (−19%) (Table 4 ). This matched with the drops of leaf area and water content of the shoots (Table 3) and with some significant changes in the leaf ionome (e.g., increase of foliar Zn, Cu, Cd, and As concentrations; decrease of foliar P and Mg concentrations, Table 5 ). The excess metal(loid)s and lesser P and Mg concentrations in leaves would impact the chlorophyll biosynthesis (Guo et al. 2016; Kolbas et al. 2018) . Instead, the chlorophyll fluorescence remained constant across the soil series; Values are means ± SD (n=4); letters indicate significant differences between treatments for each parameter at p <0.05 (post hoc Tukey HSD test) The foliar and stem concentrations of several elements were affected by the increased proportion of B soil in the soil series (Tables 5 and 6 ). Foliar As, Cd, Cu, Mo, Pb, and Zn concentrations progressively increased with the addition rate of B soil (Table 5 ). Highest increases in foliar concentrations were for As (14.2 fold), Zn (7.5 fold), and Mo (4.2 fold) in the B plants as compared to the Ctrl ones. The root uptake of As and Mo oxyanions would be favored by the alkaline soil pH. In contrast, foliar Mg concentrations decreased since the 25% addition rate, and foliar Ca and P concentrations dropped to the 100% one. Foliar Al, Cr, Co, Fe, Ni, and K concentrations did not differ across the soil series. Similarly, stem As, Cd, Cu, Mo, Ni, Pb, and Zn concentrations increased along the soil series and peaked for the B plants, the largest increases being for As (tenfold) and Zn (eightfold) ( Table 6 ). Conversely, stem P concentration decreased for the B plants, stem Na concentration dropped for the 25%, 50%, and 75% addition rates, while stem Al, Ca, Cr, Co, Fe, Mg, Mn, and K concentrations did not significantly vary across the soil series. Foliar As and Zn concentrations were respectively 3 and 1.6 times higher than the stem ones. These changes in leaf and stem ionomes were likely induced by the increase in metal(loid) exposure and the shoot biomass trend to progressively decrease, this one being correlated to the addition rate of B soil (Y = 0.00008 X 2 -0.038 X + 8.84, R 2 = 0.96; the mineral elements were less diluted in the reduced biomass of aboveground plant parts). Our foliar and stem Cd concentrations were relatively low across the soil series (0.5-1.7 and 0.2-0.5 mg Cd kg −1 DW, respectively, Tables 5 and 6) despite the high total soil Cd in the B soil. This likely reflected the low concentration of 1M NH 4 NO 3 -extractable soil Cd in this B soil (Table 1 ). In Erdem et al. (2012) , the shoot DW yield of tobacco was only decreased over 4.5-6.3 mg Cd kg −1 in the shoots depending on genotypes. Cadmium concentration in the leaves and the stem of N. tabacum cv. PBD6 respectively reached 6.8 and 1.2 mg kg −1 and 164.5 and 28.6 mg kg −1 in a control (0.44 mg Cd kg −1 DW soil) and Cd-spiked (5.44 mg Cd kg −1 DW soil) acid sandy-clay soil (typical of the Bergerac area, France, and usually cultivated with tobacco) without visible phytotoxicity symptom (Mench et al. 1989) . Foliar Cu concentration of tobacco plants across the soil series just reached its upper critical threshold range (15-30 mg Cu kg −1 , Macnicol and Beckett 1985) . Stem Cu concentrations (Table 6 ) matched with the range of shoot Cu concentrations reported in Kolbas et al. (2020) for the BaG tobacco cultivar grown in field plots with a Cu/PAH-contaminated soil. According to Majsec et al. (2016) , Cd and Cu can influence each other uptake, with Cu reducing Cd translocation to tobacco shoots. Cadmium accumulation in the tobacco leaves would depend on a coordination of Cd transport, including less cell wall binding, weaker impediment by the Casparian strip, and efficient xylem loading (Huang et al. 2021) . Foliar Zn concentrations for the tobacco plants grown from 50% of B soil (Table 5 ) exceeded the physiological non-toxic levels reported for most plant species (20-150 mg kg −1 ). At the Bettwiesen site, Switzerland (total soil Zn 400-55,000 mg Zn kg −1 , extractable soil Zn 16 mg Zn kg -1 ), the BaG tobacco shoots reached up to 617 mg Zn kg −1 DW, which was higher than the foliar Zn concentration of tobacco grown on the B soil (Table 5 ), leading to 20.1 kg Zn ha −1 yr −1 for shoot Zn removal and a cleanup time of 10 years for the bioavailable Zn stripping (Herzig et al. 2014) . At the Lommel site, Belgium (sandy soil, soil pH 6.4-6.8), the shoot concentrations of BaG tobacco reached 250-475 mg Zn and 10-27 mg Cd kg −1 DW, Cd being more uptake than at the Borifer site (Thijs et al. 2018) . Some papers reported on the root-to-shoot translocation of Zn and Cd in tobacco and Zrt-/Irt-like proteins (ZIPs), which carry Zn from extracellular space and intracellular compartments to the cytosol (Kozak et al. 2019) . High Zn and Cd exposures induced changes in root part-specific expression pattern of such membrane transporters, i.e., NtZIP1-like, NtZIP5-like, NtZIP8, NtZIP11, NtIRT1, and NtIRT1-like (Palusińska et al., 2020) . The NtZIP11 plasma membrane protein, being highly expressed in tobacco leaves and upregulated by high Zn concentrations (Kozak et al. 2019) , may contribute to accumulate Zn in leaves in conditions of Zn excess as in the case of B soil. With high Zn exposure in hydroponic, Zncontaining biogenic calcite and other Zn compounds are produced through the tobacco trichomes, and this would be a mechanism involved in Zn detoxification (Sarret et al. 2006) . The BCF values based on leaf ionome and total soil content culminated for K and P, and then for Mg and Ca (Table 5) . For metal(loid)s, the BCF values peaked for Cd, then As and Mo, but in general, these values remained low (< 1) for tobacco plants cultivated in the B soil as compared for instance with values reported by Yang et al. (2017) . Nicotiana tabacum genotypes are Cd accumulator in general, for example, as compared to N. rustica L. (Mench et al. 1989; Herzig et al. 2014; Yang et al. 2017; Huang et al. 2021) , but here, the alkaline soil pH in the B soil likely limited the Cd phytoavailability as indicated by the low extractable soil Cd. The BCF value of Zn was lower than that of Cd, but this was due to the high value of total soil Zn and its chemical soil speciation. A multivariate analysis based on the data of the stem ionome and the other plant parameters (RDA, Fig. 4A ) evidenced that on the first component (70% variability explained) were most of the metal(loid) concentrations and the morpho-physiological parameters, even though both categories of vectors were in opposite directions. Along the second component (8% variance explained), stem Mg, Ca, and P concentrations and leaf area followed the same direction. The leaf area and stem P concentration were positively correlated, and stem Ca and Mg concentrations as well. This analysis confirmed that as the root exposure to metal(oid)s increased in the soil, so the shoot biomass, the shoot water content, the total chlorophyll, and the leaf area diminished. In contrast, stem P, Ca, and Mg concentrations followed a different trend in comparison to the metal(oid)s, probably because these essential nutrients were present in the soil in sufficient concentrations to satisfy the physiological needs of the plants. The decrease in the ratio of leaf Mg concentration to stem Mg concentration as the soil contamination increased corresponded with the decrease in leaf area and chlorophyll pigments, which could explain a lower demand in Mg, which is a co-factor in the protochlorophyllide, an intermediate in the biosynthesis of chlorophyll a. The same multivariate analysis for the tobacco leaves (RDA, Fig. 4B ) indicated similarities in the behavior of the vectors: opposite directions along the first axis (66% variance explained) of most of the metal(loid)s (except Cr) and the morpho-physiological parameters. Conversely, foliar Ca, P, and Mg concentrations were once again mostly along the second axis (8% variance explained) or in between both axes. Most foliar metal(loid) concentrations were negatively correlated with foliar P and Ca concentrations. As for the stems, all the plant parameters were negatively correlated with most foliar metal(loid) concentrations and were positively correlated with the foliar P and Ca concentrations. The position or the Cr vector would reflect its almost constant foliar concentration. Vulnerability to drought-induced embolism (P50), which is commonly used to characterize drought tolerance across plant species (Choat et al. 2012) , can be altered by ion-mediated changes in xylem pit membrane porosity . The objective was to quantify the extent to which increased soil contamination and metal(loid) uptake can affect drought tolerance in tobacco. The P50 values varied from −1.45 to −2.15 MPa and did not show any significant differences across the soil series (Fig. 5A, p=0.1263) . However, A Vulnerability to drought-induced embolism (hydraulic safety, P50, MPa, differences between mean values were not significant at the 5% level). B Specific hydraulic conductance (hydraulic efficiency, Ks, m 2 MPa −1 s −1 , 95% confidence limits, R 2 =0.43) of tobacco stems depending on soil treatments (% of B soil) high P50 values at the 75% and 100% addition rates matched with an increasing trend for stem Mg, Mn, and Na concentrations of the B75 and B plants as compared to those of the B50 plants (Fig. 5A , Table 6 ). The same results (no significant differences) were obtained for the other three hydraulic safety traits (P 12 , P 88 , and slope) with slightly higher slope values for increasing soil contamination. On the contrary, hydraulic efficiency (Ks) declined significantly with increasing soil contamination (Fig. 5B) . These patterns can be explained by the fact that a salt solution lead to an immediate increase in hydraulic conductance (Van Leperen et al. 2000) while sap ionic composition does not alter xylem vulnerability to embolism . Increasing soil contamination in this soil series, including excess metal(loid)s and organic xenobiotics, progressively affected the following tobacco parameters: leaf area, shoot water content, total Chl, Chl A, and Chl B densities, and the Chl A/ Chl B and total Chl/carotenoid ratios, notably for the B plants. The capacity of tobacco to keep its photosynthetic system efficient depends on the level of metal(loid) exposure and uptake, Zn stress having less effect than Cd stress on the photosynthetic function of tobacco leaves (Zhang et al. 2020) . Here, the Chl fluorescence remained constant and close to the optimum value, which means that the PSII efficiency was almost correct and that this tobacco cultivar can manage the stress caused by the contaminant exposure, maintaining the functionality of its physiological system, while some morphological traits were more disadvantaged. It is difficult to evidence the most stressful soil contaminants involved, especially since synergies/antagonisms between bioavailable contaminants taken up are possible in the cascade of biological processes. On one hand, the most disturbing changes evidenced in the leaf ionome were probably the increase in Cu and Zn concentrations, the drop of K/Na ratio, and the decrease in P and Mg concentrations (Table 5) . Interestingly, the phytoavailability of most metal(loid)s was low in the B soil regarding their total soil concentration. On the other hand, exposure to PAH excess can affect the anatomical structure and architecture of roots and water uptake (Harvey et al. 2002; Kummerova et al. 2013 ). In the B soil, fluoranthene and then pyrene, benzo(b)fluoranthene, and benzo (g,h,i)perylene to a lesser extent, had the highest concentrations in the B soil (Table 2) . Gworek et al. (2016) reported that root uptake from the soil solution occurs for the four-and five-ring PAHs (e.g., pyrene, fluoranthene), and that for three-ring PAHs (e.g., fluorene, phenanthrene, and anthracene) and dicotyledonous plants, both uptake by roots and by leaves would occur. Such PAHs can impact the photosynthetic processes, carbohydrate allocation, and ionome of plant parts and induce an oxidative stress (Kummerova et al. 2013; Dupuy et al. 2015; Tomar and Jajoo 2017) . Here, with the aging of organic contaminants in the B soil, it could be assumed that their negative influence is not evident, especially since overall tobacco growth was little affected (Table 3) . However future trials at this Borifer site will need to determine the potential PAH transfer to the aboveground biomass. This BaG tobacco cultivar can contribute to phytoextract the phytoavailable soil Zn and Cd from the B soil. Shoot Zn and Cd removals computed with the biomass (DW) of leaves and stem and their respective Zn and Cd concentrations (Tables 5 and 6) were 2.35 mg Zn and 7.6 μg Cd plant −1 . This represented 9% for Zn and 44% for Cd of their 1M NH 4 NO 3 extractable amount in the potted B soil. Based on cultural practices, i.e., plant density of 85,000-100,000 ha −1 (Sheen 1983) , the tobacco phytoextraction would be 199-234 g Zn and 0.65-0.76 g Cd ha −1 . Such values were low as compared to the BaG tobacco harvested at the Bettwiesen site (i.e., 20.1 kg ha −1 yr −1 , with 32.5 t DW ha −1 for the shoot biomass, Herzig et al. 2014 ). Here, due to the 2-month exposure and limited volume of potted B soil that can be explored by the root system, the shoot biomass was only equivalent to 6.04 t DW ha −1 . For comparison purposes, soil pH was neutral at the Bettwiesen site (7.1) but 0.1M NaNO 3 -extractable concentrations (in mg kg −1 ) were 0.15-55 for Zn and 0.003-0.01 for Cd (Herzig et al. 2014) , therefore of the same magnitude order as compared to the Borifer soil (Table 1) . Indeed, such bioavailable Zn and Cd stripping and its potential influence on pollutant linkages (e.g., herbivory, Grignet et al. 2020 ) must be in situ investigated at the Borifer site to get a better assessment. In such sandy soil, it is also possible to harvest the root systems of tobacco, which generally display high Pb and Cu concentrations (Kolbas et al. 2020) . Herbivory exposure with tobacco was not evidenced in the field trial at the St-Médard d'Eyrans site (Gironde County, France), except for snails and slugs at the 6-7 leaf stage just after transplanting on site (Kolbas et al. 2020) . The resupply of the phytoavailable soil fraction will be regulated by the metal(loid)-bearing solid phases and the soil properties, including the soil pH and Eh. Other Zn/Cd-phytoextraction options are reported in Europe: e.g., Salix viminalis L. and Arabidopsis halleri (L.) O'Kane & Al-Shehbaz alone and in inter-row co-cropping at the urban Montataire site, Oise County, France (Grignet et al. 2020) , Noccaea caerulescens (J. Presl & C. Presl) F. K. Mey (Jacobs et al., 2019) , short rotation coppice of poplars and willow, and sunflower (Thijs et al. 2018) . However, the growth of these plant species would be likely limited by the hot and dry summer weather conditions in southwest France. The black poplars (Populus nigra L.) implemented at the Chaban-Delmas site in the Parc aux Angéliques, which is adjacent to the Borifer site, suffered from this and must be irrigated in the summer (Marchand et al. 2016) . The willow and poplar leaves, which accumulate Zn and Cd, are also less easy to harvest. To implement such phytomanagement at this Borifer brownfield, removal of large stones and pieces of concrete, soil loosening, compost incorporation into the soil, and crop irrigation will be crucial points to consider. Winter crops with metal(loid)-excluders, e.g., white clover, alfalfa, and Festuca sp., between tobacco cultures would also promote the carbon sequestration, soil fertility, increase in microbial activity, and the biodegradation of organic contaminants (Kidd et al. 2015; Bourceret et al. 2018) . Metal(loid) bioavailability may however change with the biodegradation of organic contaminants. Therefore, both must be monitored in parallel to changes in microbial, mesofaunal, and animal communities. Technologies for processing shoot biomass of tobacco harvested at phytomanaged sites are listed in Asad et al. (2017) , Yang et al. (2017) , Rehman et al. (2019) , and Kolbas et al. (2020) and not detailed here. Digestates from anaerobic digestion should not exceed 1000 mg Zn and 1.5 mg Cd kg −1 (Arrêté du 8 août, 2019), so this point should be considered regarding the processing of our tobacco biomass. Currently, compost incorporation into the soil combined with hay transfer to promote both the microbial and plant communities is also an option in situ investigated at the Borifer site. As tobacco is grown from May to October in southwest France, soil cover should be provided by a winter crop to avoid soil erosion by natural agents, if possible sequestering carbon and promoting the rhizodegradation of organic contaminants (e.g., alfalfa, white clover, Marchand et al. 2016; Bourceret et al., 2018) . If the climatic conditions become drier in summer due to global warming, the cultivation of aromatic plants adapted to the Mediterranean climate should be evaluated (Raveau et al. 2020 (Raveau et al. , 2021 . Growing non-food crops is an option for the phytomanagement and remediation of contaminated soils. Here, the potential ability of tobacco (cv. Badischer Geudertheimer) for phytomanaging a contaminated soil at an urban brownfield site was assessed with a potted soil series made using the fading technique. Despite the high total soil Zn, Cu, Pb, and PAH (e.g., fluoranthene), 2 months after transplantation, the decrease in shoot biomass was only 33% for the plants grown on the undiluted contaminated soil and was not significant as compared to the uncontaminated control soil (Table 3 ). In contrast, leaf area, water content of shoots, total chlorophyll, and hydraulic efficiency of stem progressively decreased with the soil contamination (Table 3 and Fig. 5) . Overall, both foliar and stem Zn and Cd concentrations increased, while the shoot biomass production was roughly maintained, allowing to phytoextract a relevant fraction of the phytoavailable soil Zn and Cd. Few nutrient concentrations in aboveground plant parts decreased with the soil contamination, e.g., foliar and stem P concentrations. This phytomanagement option needs further long-term evaluation on plots at the site, as well as agronomic and bioaugmentation practices that can optimize it and its consequences on soil ecological functions (e.g., changes in microbial and mesofaunal communities, herbivory, organic matter, and contaminant cycles), without underestimating the potential impact of increasingly frequent droughts and heat waves and water availability being a crucial factor for this sandy soil at this brownfield site. Supplementary Information The online version contains supplementary material available at https://doi.org/10.1007/s11356-021-16411-y. La reconversion des friches polluées au service du renouvellement urbain : enseignements technico-économiques -Bilan des opérations aidées dans le cadre du dispositif ADEME d'aide aux travaux de dépollution pour la reconversion des friches polluées approuvant deux cahiers des charges pour la mise sur le marché et l'utilisation de digestats de méthanisation agricole en tant que matières fertilisantes. JORF n°0221 du 22 septembre 2019 Pretreatment of trace element-enriched biomasses grown on phytomanaged soils for bioethanol production Biodegradation of organic xenobiotic pollutants in the rhizosphere Teneurs totales en éléments traces métalliques dans les sols (France), INRA Editions The development and use of sustainability criteria in SuRF-UK's sustainable remediation framework Integrated and sustainable management of post-industrial coasts Bioremediation of polyaromatic hydrocarbons (PAHs) using rhizosphere technology Determination of chlorophyll with DMF High PAH degradation and activity of degrading bacteria during alfalfa growth where a contrasted active community developed in comparison to unplanted soil Zinc tolerance and accumulation in the salt-marsh shrub Halimione portulacoides ed) Zinc in soils and plants. Developments in plant and soil sciences Global convergence in the vulnerability of forests to drought Evaluation of a new centrifuge technique for rapid generation of xylem vulnerability curves The effects of sap ionic composition on xylem vulnerability to cavitation Brownfields to green fields: realising wider benefits from practical contaminant phytomanagement strategies Mechanism of waterstress induced cavitation in conifers: bordered pit structure and function support the hypothesis of seal capillary-seeding Morphological and physiological responses of maize (Zea mays) exposed to sand contaminated by phenanthrene Canadian Environmental Quality Guidelines. Carcinogenic and other polycyclic aromatic hydrocarbons (Environmental and Human Health Effects) Effect of cadmium stress on growth and mineral composition of two tobacco cultivars European Parliament resolution of 28 April 2021 on soil protection (2021/2548(RSP)) Urban soil phytomanagement for Zn and Cd in situ removal, greening, and Zn-rich biomass production taking care of snail exposure Magnesium deficiency in plants: an urgent problem Comparison of PAHs uptake by selected monocotyledones and dicotyledones from municipal and industrial sewage sludge Phytoremediation of polyaromatic hydrocarbons, anilines and phenols Feasibility of labile Zn phytoextraction using enhanced tobacco and sunflower: results of five-and one-year field-scale experiments in Switzerland The adoption of sustainable remediation behaviour in the US and UK: a cross country comparison and determinant analysis Differential cadmium translocation and accumulation between Nicotiana tabacum L. and Nicotiana rustica L. by transcriptome combined with chemical form analyses Sols. Méthodes Applicables Aux Sols Influence of edaphic conditions and nitrogen fertilizers on cadmium and zinc phytoextraction efficiency of Noccaea caerulescens Are nettle fibers produced on metalcontaminated lands suitable for composite applications? Agronomic practices for improving gentle remediation of trace-elementcontaminated soils Morphological and functional responses of a metal-tolerant sunflower mutant line to a copper-contaminated soil series Field evaluation of one Cu-resistant somaclonal variant and two clones of tobacco for copper phytoextraction at a wood preservation site NtZIP11, a new Zn transporter specifically upregulated in tobacco leaves by toxic Zn level Root response in Pisum sativum and Zea mays under fluoranthene stress: morphological and anatomical traits Mechanisms of metal toxicity in plants Determinations of total carotenoids and chlorophylls a and b of leaf extracts in different solvents Monitoring of a long term phytoremediation process of a soil contaminated by heavy metals and hydrocarbons in Tuscany Critical tissue concentrations of potentially toxic elements Integrative approach gives new insights into combined Cd/Cu exposure in tobacco Parc des Angéliques -Quai de Brazzasite BORIFER -33000 BORDEAUX. Diagnostic de contamination Plant responses to a phytomanaged urban technosol contaminated by trace elements and polycyclic aromatic hydrocarbons Cadmium bioavailability for Nicotiana tabacum L., Nicotiana rustica L. and Zea mays L. grown in soil added or non-added with cadmium nitrate Phytomanagement and remediation of Cucontaminated soils by high yielding crops at a former wood preservation site: sunflower biomass and ionome Phytomanagement of metal(loid)-contaminated soils: options, efficiency and value Use of sodium polytungstate in the granulo-densimetric fractionantion of soil material Potential of five plant species for phytoremediation of metal-PAH-pesticide contaminated soil. Electronic Thesis and Dissertation Repository Zn/Cd status-dependent accumulation of Zn and Cd in root parts in tobacco is accompanied by specific expression of ZIP genes A mathematical and statistical analysis of the curves illustrating vulnerability of xylem to cavitation European achievements in soil remediation and brownfield redevelopment Preliminary study of phytoremediation of brownfield soil contaminated by PAHs The effect of plants and natural compounds on bacterial population in the contaminated soil The culture of clary sage shaped the rhizospheric bacterial communities more strongly than mycorrhizal inoculation in the trace elementcontaminated soil -a two-year monitoring field trial Lounès-Hadj Sahraoui A (2021) In situ cultivation of aromatic plant species for the phytomanagement of an aged-trace element polluted soil: plant biomass improvement options and techno-economic assessment of the essential oil production channel Opportunities and challenges in the remediation of metal-contaminated soils by using tobacco (Nicotiana tabacum L.): a critical review Finnish policy approach and measures for the promotion of sustainability in contaminated land management The effect of vegetation on decrease of PAH and PCB content in longterm contaminated soil Base de données du Référentiel Régional Pédologique du département de la Gironde à 1/250 000 Comparison of international approaches to sustainable remediation ): A potential candidate for phytoremediation? Biological and economical points of view Jute: a potential candidate for phytoremediation of metals-a review Ionomics and the study of the plant ionome Trichomes of tobacco excrete zinc as zinc-substituted calcium carbonate and other zinc-containing compounds Biomass and chemical-composition of tobacco plants under high-density growth Soil contamination in China: current priorities, defining background levels and standards for heavy metals Tobacco, sunflower and high biomass SRC clones show potential for trace metal phytoextraction on a moderately contaminated field site in Belgium PSI becomes more tolerant to fluoranthene through the initiation of cyclic electron flow Xylem embolism threshold for catastrophic hydraulic failure in angiosperm trees Van Ieperen W, Van Meeteren U, Van Gelder H (2000) Fluid ionic composition influences hydraulic conductance of xylem conduits Thresholds of copper toxicity to lettuce in field-collected agricultural soils exposed to copper mining activities in Chile First evidence of large-scale PAH trends in French soils Functional activity and functional gene diversity of a Cu-contaminated soil remediated by aided phytostabilization using compost, dolomitic limestone and a mixed tree stand Phytoextraction of cadmium-contaminated soil and potential of regenerated tobacco biomass for recovery of cadmium Toxic effects of heavy metal Cd and Zn on chlorophyll, carotenoid metabolism and photosynthetic function in tobacco leaves revealed by physiological and proteomics analysis Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgements The City of Bordeaux (Bordeaux Metropole) has given access to the site and its history. Tobacco seeds were initially provided by Dr. R. Herzig, Phytotech, Bern, Switzerland. The UMR BIOGECO is member of the INRAE ecotoxicology network (ECOTOX, https://www6.inrae.fr/ecotox_eng/) and the COST Action 19116 PlantMetals (https://plantmetals.eu/plantmetals-home. html). Eliana Di Lodovico has received a 3-month grant for Traineeship (SMT) from the Erasmus Plus Mobility, The European Commission.Data availability The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Ethics approval and consent to participate Not applicable Competing interests The authors declare no competing interests.