Evaluation of Bio-Based Extraction Methods by Spectroscopic Methods

: New technologies are in development regarding the preservation of waterlogged archaeological wood items contaminated with Fe / S species. To this purpose, a bio-based treatment to extract these harmful species before further damages occur is presented. Thiobacillus denitriﬁcans and desferoxamine were employed based on their speciﬁc properties to solubilize iron sulﬁdes and uptake iron. The biological treatment was compared with oxidizing and complexing agents (sodium persulfate and ethylene diamine tetraacetate) traditionally used in conservation-restoration. Mock-ups of fresh balsa as well as fresh and archeological oak and pinewood were prepared to simulate degraded waterlogged wood by immersion in corrosive Fe / S solutions. The e ﬃ ciency of both biological and chemical extraction methods was evaluated through ATR-FTIR and Raman spectroscopies and validated by statistical approach. Results showed that treatments did not a ﬀ ect the wood composition, meaning that no wood degradation was induced. However, the chemical method tended to bleach the samples and after treatment, reduced sulfur species were still identiﬁed by Raman analyses. Finally, statistical approaches allowed validating ATR-FTIR results.


Introduction
Anoxic environments promote the preservation of waterlogged archaeological wooden objects. When artefacts are, for example, buried in sediments, the surrounding environment allows the wood structure to perdure [1]. The cellulosic content is however subjected to a slow deterioration by erosion bacteria (EB) as primary wood degraders under these conditions [2,3].
First, the lignocellulose structure of the cell walls is broken by EB [4,5]. The wood surface is degraded allowing the bacterial invasion to reach the cell lumen and continuing through the cellulose-rich layer. In this layer, the crystallized cellulose is converted into an amorphous substance [1,6,7] that is then used by secondary wood degraders such as sulfate-reducing bacteria (SRB) as carbon and energy source. In parallel, in sulfate-rich environments, hydrogen sulfide H 2 S is produced by the reaction of SRB with organic carbohydrates of wood, according to Equation (1) [8]: 2(CH 2 O) + SO 4 2− + 2H + → H 2 S + 2H 2 O + 2CO 2 (1) H 2 S diffuses within the wood and reacts with iron ions, if present (i.e., as nails or in the burial environment). Iron sulfides are then formed and accumulated within the wood timbers according to Equation (2) [4,9]: Transfer of the artefacts from an anoxic to an oxic environment reveals severe conservation issues. Indeed, iron sulfides are oxygen-sensitive [10]. Acidification and salts efflorescence occur leading to irreversible physical and chemical damages. The pH drop enhances the degradation of cellulose through an acidic hydrolysis resulting in a loss of strength. In addition, formed iron oxides, iron sulfates and/or elemental sulfur induce cracks of the cell structure as these salts are more voluminous than their precursor iron sulfides [11]. Moreover, iron itself catalyzes the oxidation of sulfur compounds and cellulose degradation, accelerating the decay process [12,13].
Iron and sulfur accumulation within waterlogged archaeological wood is of great concern for conservation and preservation. It occurs frequently in marine environments but could also occur on sweet water finds. Several important wooden cultural heritage artefacts such as the Vasa warship in Sweden or the Mary Rose in United Kingdom present acidification and salt precipitation issues and require a long-term monitoring [14].
Usually, polyethylene glycol (PEG) is the main consolidation agent employed, either by spraying or immersion [5,12]. However, this method does not prevent from acidification and salts efflorescence. Moreover, iron keeps corroding in contact with moist PEG inducing even more degradation of the wooden structure [15]. Complexing agents such as derivatives of ethylene diamine tetra acetic (EDTA), employed to extract iron, present high extraction yields but do not avoid wood swelling and are time-consuming [16]. In addition, neutralizing methods are also performed to remediate the acidification. Alkaline baths and ammonia gas are the most efficient methods employed so far but can induce holocellulose degradation [4].
Alternatively, less aggressive methods are currently being developed. For instance, bio-based technologies have proven to be sustainable solutions for the preservation of cultural artefacts [17]. Microbial processes can be used for consolidation, cleaning or inhibition of reactive corrosion products while reducing the toxicity related to traditionally used treatments [18]. For metal conservation, Beauvaria bassiana fungal strain produces oxalic acid in presence of toxic metals such as copper. The copper is then transformed into copper oxalates. Applied on outdoor bronze sculptures, this biotechnology allows forming a protective and stable layer preventing the object from further degradation and leaching [19]. Microorganisms can also be employed to stabilize corroded layers of archaeological iron artefacts and akaganeite β-FeOOH, a chlorine-containing iron oxyhydroxide, can be converted into chemically stable vivianite Fe 3 (PO 4 ) 2 ·8H 2 O, siderite FeCO 3 or magnetite Fe 3 O 4 . As well, few studies related bio-based conservation methods on waterlogged archaeological wood and the removal of mackinawite FeS and pyrite FeS 2 , the main iron sulfides detected in such artefacts [20].
The goal of the present research is to develop an innovative, green and sustainable extraction method using the specific capacities of some microorganisms, such as Thiobacillus denitrificans, to extract harmful Fe/S species. In fact, T. denitrificans, a facultative anaerobe bacterium, can employ nitrate NO 3 − as electron donor to oxidize sulfur compounds at neutral pH. Several case studies showed its faculty to dissolve iron sulfide compounds [14,21]: In addition, siderophores were investigated for the extraction of iron on organic materials [22]. In fact, siderophores are produced by many bacteria under iron-deficient conditions and they present the highest iron affinity, allowing their use in pharmacology, agriculture, medicine and other applications [23]. For example, the application of siderophores was investigated on paper and wooden artefacts to remove iron staining [24]. The results showed that siderophores were the most effective regarding the extraction of iron in comparison to commonly used complexing agents. They represent a promising green alternative to current chemical iron extraction methods [25].
The efficiency of the proposed bio-based treatment regarding Fe/S species extraction is evaluated by comparing with the main chemical methods usually employed on different wood mock-ups. Further consolidation with PEG and freeze-drying is also foreseen to check the compatibility of bioremediation with the following consolidation and drying steps. Attenuated total reflectance infrared (ATR-FTIR) and Raman spectroscopies analyses were performed to evaluate the wood degradation and Fe/S extraction rates. First results showed promising results with the biological treatment as the appearance of biologically treated samples was less impacted than the chemically treated samples. Moreover, ATR-FTIR analyses did not show any wood degradation and Raman revealed that reduced sulfur compounds were still present only after chemical extraction. Statistical analyses were also performed on the different datasets to validate the previous observations achieved.

Wood Samples
Three wood species were selected: balsa, oak and pine. Balsa wood (Ochroma pyramidale) is a light and porous wood, with a structure similar to the one of waterlogged archaeological wood. It has already been used in some research projects in mock-ups to model waterlogged wood [26]. Oak (Quercus sp.) and pine (Pinus sp.) are the main wood species found on waterlogged wood artefacts and the most described cases in literature [16,27]. Fresh oak wood has been also reported as a mock-up to mimic waterlogged wood [13,28].
Fresh and archaeological wood samples were used. Fresh balsa was purchased from Coop Brico + Loisirs (Neuchâtel, Switzerland). Fresh oak and pine were provided from a carpenter, J.F. Liabeuf (Cernier, Switzerland). Concerning archaeologic samples, poles and boards of Neolithic oak and pine (respectively) were provided by the Archaeological Service of Bern Canton (ADB) and two Neolithic oak poles by the Swiss National Museum (SNM).
Cubes samples of 2 × 2 × 2 cm 3 were cut. Concerning archaeological wood, the samples were cut on the external layer of the poles (earlywood) and following the wood ring to have homogeneous samples sets. All the samples present the transversal (Tv), Tangential (Tg) and Radial (Rd) cutting sections.
For each wood species and type, 10 sets of 19 samples were prepared. For each set, one sample was kept as reference, on which no artificial contamination or extraction was performed. The set distribution is described in Table 1. All samples were stored in degassed deionized water as stock solution. Sets A1, A2, B, C, E and F were artificially contaminated with iron and sulfur species according to the protocol described in Albelda-Berenguer et al. [14]. The purpose of this impregnation is to form iron sulfides compounds such as mackinawite (FeS) and/or pyrite (FeS 2 ). A water pre-immersion treatment under vacuum (600 mbars) was also applied for sets A1, A2, B and E before impregnation to enhance the desired wood degradation and to simulate the state of degradation observed on real waterlogged wood artefacts. Regarding D1, D2, G1 and G2 sets, spot tests to detect Fe 2+ , Fe 3+ , sulfides and sulfates were performed and revealed the presence of iron and sulfates and in some parts sulfides. The tests were performed with potassium thiocyanate 33% and barium chloride to detect iron and sulfates, respectively. Machery-Nagel ® test paper n • 90761 (Fisher Scientific AG, Reinach, Switzerland) was employed for sulfides detection [29].

Extraction Methods
Each set was photographically documented before and after extraction ( Figure 1). Before extraction, the samples were put in 400 mL of deionized water and sonicated for 2 and then 4 min. The water was changed between the two-sonication baths. The sonication allows removing the chemicals employed for impregnation. Two baths were necessary to have a clear solution. Six samples per set were used for each extraction method tested and labelled as biologically treated (BT) and chemically treated (CT) samples respectively (see Sections 2.2.1 and 2.2.2). In parallel, six samples remained in deionized water as control and labelled as non-treated (NT) samples. The volume of extraction solution was set to 200 mL. At the end of treatment, BT and CT samples were put in degassed deionized water and several cycles of 2 × 20 min vacuum (1 mbars) followed by 1 h sonication were applied on the samples until constant pH was reached and storage solution conductivity was in the range 10-30 µS/cm −1 (as the rinsing water).

Biological Extraction
T. denitrificans (DSMZ 12475) was obtained from the German Collection of Microorganisms and Cell Cultures (DMSZ, Braunschweig, Germany). The strain was cultivated using standard anaerobic techniques at 30 • C in the dark. Desferoxamine (DFO) was purchased as Desferal ® (Novartis Pharma Schweiz AG, Basel, Switzerland). The BT samples were immersed in a solution of DFO 84 mM for 10 days. Samples were kept at room temperature and with orbital agitation (100 rpm). After rinsing with deionized water, the samples were immersed with T. denitrificans for 20 days. Samples inoculated with the bacteria were incubated at 30 • C in darkness without agitation.

Chemical Extraction
The applied protocol was based on a previous study [26]. The samples were immersed in sodium persulfate 0.1 M for 1 day to oxidize sulfur species. After rinsing with deionized water, the samples were then immersed in EDTA 0.125 M for 7 days to complex free iron ions. During the whole process, the samples were kept at room temperature without agitation.

ATR-FTIR Spectroscopy
Wood degradation was evaluated by Fourier Transform InfraRed spectroscopy. Wood mock-ups were characterized in the 4000-650 cm −1 range using a Nicolet iS5 instrument (Thermo Fisher Scientific, Waltham, MA, USA) equipped with an Attenuated Total Reflectance accessory, with 16 scans, resolution of 4 cm −1 and an analyzed area of 62.5 µm. Two samples per treatment (BT, CT and NT) per set were analyzed and three spectra per five faces were recorded, for a total of 90 spectra per set. The height of vibrational bands at 1034 cm −1 assigned to holocellulose H and lignin L, 1158, 1374 cm −1 assigned to holocellulose H, cellulose C and 1506 cm −1 assigned to lignin L were used to calculated three ratios: R1 = I(1158)/I(1506), R2 = I(1374)/I(1506) and R3 = I(1034)/I(1506) according to literature [30,31]. Bands height of spectra were measured using Rstudio software (Rstudio, Boston, MA, USA) and ChemoSpec package (Bryan A. Hanson, DePauw Univeristy, https: //bryanhanson.github.io/ChemoSpec/index.html). Heights of bands were calculated from the baseline, corrected with a polynomial function, and normalization of the spectra. These ratios allow determining which wood components were affected by the treatments and thus the wood state of degradation.

Raman Spectroscopy
Compounds present on the wood samples were identified by Raman spectroscopy. The measurements were performed on wet samples, immersed in their stock solution (degassed water). This methodology allows avoiding the overwarming of the samples that would induce a possible modification of the compounds under the laser excitation.
Analyses were carried out using a LabRam Aramis Horiba instrument (HORIBA Jobin Yvon GmbH, Bensheim, Germany) and LabSpec 5 software with a 632.8 nm laser, in the range 100-1500 cm −1 , with 100× objective (numerical aperture of 0.9) and laser power of 0.99 mW at the surface of the sample (D1 filter, 10%), 1800 lines/mm grating, 1000 µm confocal pinhole and 100 µm spectrometer entrance slit. Two samples per treatment (BT, CT and NT) per set were analyzed. One spectrum at the center of three faces was recorded. Two spectra as routine ones (10 acquisitions of 5 s) and as another with a long acquisition (25 × 10 secs). A total of 18 spectra per set were recorded.
In parallel, reference compounds (mackinawite FeS, goethite α-FeOOH, hematite Fe 2 O 3 ) were synthetized according to Bourdoiseau et al. and Schwertmann [32,33]. These precipitates as well as sulfur (Sigma, α-S 8 ) and pyrite (http://www.pierrequiroule.ch/, FeS 2 ) were used to create a reference spectral dataset. Spectra were acquired with the same parameters as above. Moreover, water pre-immersed wood and archaeologic samples of different wood types were used to obtain reference spectra of the wood substrates.

Statistic Approach
Chemometric analyses were performed to validate which extraction method was the most efficient in terms of wood degradation and Fe/S extraction rates. Principal component analyses (PCA) were carried out on all spectroscopic measurements. Deeper investigations were achieved on spectroscopic data with Least Absolute Shrinkage and Selection Operator (LASSO) analyses. Rstudio software was used with the packages Chemospec, ggplot2, and glmnet. The spectra dataset was corrected (baseline, bin, normalization, Savitzky-Golay filter) before any analyses.

ATR-FTIR Spectrosocpy
Spectra recorded on untreated samples were considered as a starting point (t0) for the characterization campaign. Before extraction, differences could be observed between the spectra from fresh and archaeological wood sets (Figure 2A). Vibrational bands at 1034 and 1059 cm −1 were very intense for fresh wood and decreased in intensity for archaeological wood sets. Only set F (Figure 2A, light green line) which is lake pinewood presented a mean spectrum similar to fresh samples. These bands were assigned to C-O stretching of holocellulose and lignin [34]. Moreover, the intensity of vibrational band at 1738 cm −1 assigned to C=O stretching in xylan, a component of hemicellulose [35], decreased for archaeological wood sets. All these observations suggested a degradation of holocellulose content within the archaeological wood sets. On the contrary, archaeological samples displayed intense band at 1650 cm −1 attributed to H 2 O deformation due to the water absorbed by waterlogged wood during burial time (Figure 2A). The vibrational bands at 1237 and 1267 cm −1 were assigned to C-O stretching (lignin and xylan) and to C-H and O-H wagging (cellulose and xylan) [13,34,35]. Generally, a higher intensity of these bands was observed for fresh wood presenting differences between the wood species.
After extraction, BT and CT extraction methods presented similar spectra while NT samples showed different bands intensities. The bands at 1034 and 1059 cm −1 were still observed for fresh samples ( Figure 2B, full lines). However, the intensities of bands in the range 1300-1200 cm −1 decreased. All sets presented adsorbed water with a characteristic band at 1650 cm −1 . Archeological samples still had bands less intense in the range 1400-1200 cm −1 suggesting that the degradation of fresh wood was not comparable to in situ degradation ( Figure 2B, dashed lines). Finally, none of the sets displayed the band at 1738 cm −1 , meaning that xylan content was probably degraded for all wood samples ( Figure 2B).  In addition, ATR-FTIR ratios were calculated (Table 2). For fresh wood, we observed a general decrease of the ratios suggesting a degradation of holocellulose content after artificial degradation and contamination (i.e., sets A1, A2, B and E). In addition, pinewood samples from sets E and F showed higher ratios values.  We noticed that after extraction, whatever the treatment (i.e., BT or CT) applied, the ratios calculated were quite similar with the one before extraction (NT). This observation suggests that the extraction methods did not have any impact on the wood matrix. This is even more evident for not artificially contaminated sets (i.e., sets D1, D2, G1 and G2).

Raman Spectroscopy
First, the reference compounds were analyzed. Raman spectra obtained are displayed in Figure 3. Spectra obtained from goethite, hematite, sulfur and pyrite were similar to that described in the literature [10,36]. For the synthetized mackinawite, the bands obtained at 119(s), 235(w), 282(s), 320(w), 360(s) and 452(m) cm −1 suggest the formation of an intermediary phase, partially oxidized mackinawite Fe 1−x S [10]. Analyses of the wood samples before extraction showed that reduced sulfur compounds were present on the artificially contaminated samples. A1, A2, C, E and F sets displayed the same spectrum as reference Fe 1−x S ( Figure 4B). Elemental sulfur was identified only on set B samples (fresh oak) with characteristics bands at 152(m), 220(s) and 473(m) cm −1 ( Figure 4A). Finally, Raman spectra for sets D1, D2, G1 and G2 did not reveal the presence of any inorganic compounds. This observation is in contradiction with the previous spot tests performed on wood samples from sets D1, D2, G1 and G2. After extraction, no reduced sulfur compounds were detected on BT samples suggesting that the extraction seemed successfully achieved (Figure 4). While elemental sulfur was still identified on some CT samples.

Chemometrics
PCA analyses were carried out on the ATR-FTIR dataset before and after extraction. Before extraction, different clusters can be observed ( Figure 6A). All three species of fresh wood gathered with an overlapping of fresh balsa and oak samples. On the other side, all the archaeological samples clustered and overlapped. Only samples from set F (lake pine) differed slightly from the other sets, suggesting that this set is less degraded.
After extraction, all sets gathered, independently of the extraction method applied ( Figure 6B). Some clusters can be observed. Sets E and F (artificially contaminated pinewood) overlapped, A1 and A2 sets (balsa) also overlapped as for sets C, D1, D2 and G2 (archaeological oak wood). Only exception is set G1, overlapping with archaeological oak sets though G1 is archaeologic pinewood. It seemed that when a certain state of degradation is reached, the wood species does not influence the classification. These clusters are related to wood composition, which is identical within each cluster. LASSO regression was carried out on the ATR-FTIR dataset. All the analyses showed that the vibrational band at 1158 cm −1 , assigned to C-O-C vibration in holocellulose, is the main important variable regarding the classification of our data.
PCA was also performed for the Raman dataset. No distinction could be made between artificially contaminated and waterlogged wood ( Figure 6C) or among the extraction methods ( Figure 6D). Deeper analyses were performed on the entire dataset through LASSO regression.
At first, spectra of the reference compounds were used as model. Then, the model was applied on the Raman spectra dataset to predict its classification. However, the results obtained showed that our model was not optimal as hematite bands were the main variables employed for the classification. The interpretation of the experimental data needs further investigation.

ATR-FTIR Spectroscopy
Vibrational bands assigned to carbohydrates at 1034 and 1059 cm −1 were observed for fresh samples but their intensities were lower for archaeological samples (except set F). This is in correlation with the R1 (H/L) and R2 (C/L) ratios calculated. Indeed, the fresh wood samples showed ratios in the range 1.00-1.50 while ratios of archaeological wood (oak and pine) displayed lower ratio values (0.80-1.00). Previous studies of modelling waterlogged wood have demonstrated that these ratios decreased after artificial contamination of the fresh samples [14]. This is also observed here, with ratios of NT samples being lower than the ratios of fresh balsa, oak and pine wood species. Moreover, the ratios after extraction remained in the same range. It can then be ascertained that the decrease of the ratios was not due to the extraction but to the artificial contamination protocol applied on some wood mock-ups.
Comparing the treatments (BT, CT) with NT samples, we can observe that the ratios were in the same range for each set. The treatment applied on the samples did not affect the carbohydrates content, whether it was a chemical or a biological treatment. Regarding the degradation criterion, none of the extraction methods enhanced wood degradation suggesting that they all guarantee the preservation of waterlogged wood objects. Previous study has already demonstrated that T. denitrificans was harmless for wood substrate [14]. Iron extraction with siderophores as first extraction step did not interfere with the substrate neither.
However, R3 ratio displayed higher values for sets E (fresh pine) and F (lake pine). Both these sets are pinewood artificially contaminated. R3 ratio is based on the height of the vibrational band at 1034 cm −1 attributed to C-O stretching of holocellulose and lignin [34]. The artificial contamination did not degrade wood content as the values for NT samples in sets E and F are in the same range as fresh pine. The ratios increased after extraction treatment suggesting that lignin content was degraded. This is an extra correlation for these sets. Indeed, all mock-ups in set F showed very intense vibrational bands at 1304 and 1059 cm −1 even though this pinewood was dated from the nineteenth century. Comparing with literature, waterlogged pine wood generally present higher ratio values than oak species [38]. It seems that the impregnation protocol has a different effect depending on wood specie leading to different states of carbohydrates degradation. Intern structure of these pinewood species (porosity, presence of resin canals) may lead to different degradation processes. For example, Neolithic pine was reported to be a yellow and hard wood while Neolithic oak was described as yellow-red wood with the heartwood well preserved but the sapwood often destroyed [39]. To reach similar state of degradation, a longer water pre-immersion time or increase of vacuum, forcing the structure to collapse should be applied on fresh wood samples in further experiments.
Results obtained from the ratios imply that no cellulosic degradation occurred during extraction. However, visual observation showed an important aesthetic change, especially regarding chemically treated samples (Figure 1). If the ratios are not impacted, minor wood components such as extractives may react during extraction. Another reason could be that inorganic salts react with the persulfate, leading to their dissolution and then discoloration of the wood surface.
PCA analysis was carried out on the same dataset. The results showed that before artificial contamination fresh and archaeological wood samples can be differentiated. The clustering is based on the frequencies in the range 1010-1060 cm −1 and 3250-3270 cm −1 . These ranges affected the first loadings for the PCA analysis. The range 3250-3270 cm −1 is assigned to O-H stretching vibrations [34]. This proves that the lignin parts were not the main factor to cluster this dataset. PCA analysis is then based on the carbohydrates content and allows differentiating the sets. However, after extraction, as observed before for spectra and ratios, it was more complicated to cluster data. PCA score plots obtained showed that all the samples gathered. Sets E and F clustered, validating the observations above. This cluster proved that pine, whether fresh or slightly decayed, presents similar behavior. Oak samples (sets C, D1, D2 and G2) also cluster together. All these sets are from archaeological wood. Their confident ellipses overlapped the one of set G1 (archaeological pine) meaning that spectra obtained for all the sets present similarities. All archaeological oak wood sets present similar state of degradation as well as set G1. However, set B (fresh oak) did not cluster with the other oak groups but was closer to sets E and F (fresh and lake pine, respectively) meaning that the fresh oak degradation was not as efficient as expected. Both fresh oak and pine wood species were not as degraded as expected during the artificial contamination protocol. In addition, their spectra were similar to the one of set F samples. Application of artificial contamination of fresh oak and pine led to a degradation state similar to lake pine. The ranges from 810-850 cm −1 , 1020-1065 cm −1 and 1508-1558 cm −1 were the main frequencies to cluster the samples. After extraction, we noticed that the lignin range is as important as carbohydrates range meaning that this wood component interfered in the clustering. Lignin vibrational bands are more intense compare to before extraction. This could suggest that the carbohydrates vibrational bands are as intense as the lignin ones.
LASSO analysis validated these observations. All three factors Wood Species, Wood Type and Artificially Contamination mainly depended of the variable 1158 cm −1 , which is the vibrational band attributed to C-O-C vibration within holocellulose. The statistical approach performed showed that carbohydrates content was the main variable to differentiate, classify and interpret the dataset.

Raman Spectroscopy
In parallel, Raman spectroscopy was performed on the samples to identify the compounds present. Based on a comparison with the reference compounds (Figure 3), it can be proven that the impregnation protocol applied on some sets to perform an artificial contamination led to the formation of partially oxidized mackinawite Fe 1−x S. Only set B samples displayed different spectra with elemental sulfur S 8 detected. The precipitation of sulfur may be due to the wood matrix. Indeed, fresh oak wood presents a high content of tannin [40]. This was confirmed by analyzing the precipitates obtained from the water pre-immersion solution (data not showed). If tannins were released during the water pre-immersion phase, it is likely that they were also released in the impregnation solution. As the first step of the impregnation with ferrous chlorine solution lasted 4 h, it is possible that the tannins released reacted with iron ions present in solution to form iron-tannin complexes. Iron will then not be available to react with the sulfide ions during the second step of the impregnation protocol. On the contrary, elemental sulfur was not formed on archaeological oak sets since the tannins were absent as previously degraded during burial time. Moreover, the water solution where oak samples were pre-immersed turned orange while this was not the case for sets A1 and A2 (fresh balsa) and set E (fresh pine) whose pre-immersion water solutions remained colorless. The color of the solution for fresh balsa and fresh pine suggested that no release of tannins occurred during the water pre-immersion step and explain why iron sulfides were formed afterwards during the impregnation step. In literature, it was demonstrated that iron sulfides can be oxidized by microorganisms through two different mechanisms: thiosulfate and polysulfide pathways. These mechanisms depend of the acid solubility of the iron sulfide species [41]. Indeed, pyrite is an acid-insoluble sulfide while mackinawite or partially oxidized mackinawite are acid-soluble compounds [42]. Only acid-soluble iron sulfides species were here investigated, and the efficiency of BT assessed only for this family of compounds.
However, no iron sulfide compounds were identified on sets D1, D2, G1 and G2. The spot tests previously performed validated the presence of iron and sulfate species. Their concentration within archaeological wood may be low or the contamination heterogeneous, explaining why nothing was detected by Raman.
After extraction, iron sulfides compounds were absent from all BT samples. It seemed then that at least on surface, the extraction of Fe/S species with the proposed microbial treatment was efficient. On the contrary, CT samples still presented some reduced sulfur compounds, mainly elemental sulfur. Especially, as can be seen on Table 2, sulfur was identified on sets A2, E and F. As observed with previous ATR-FTIR results, different chemical extraction rates may be observed for pinewood. In fact, both sets of fresh and lake pine presented sulfur, but all iron seemed extracted. The chemical treatment allowed uptaking the iron, at least on surface, but the sulfur was not complexed and solubilized.
These first results showed that BT, even though it lasted 30 days against 8 days for CT, was more efficient regarding both iron and sulfur extraction. However, the analyses were only performed on surface of the samples and did not allow evaluating the complete extraction of Fe/S species within the wood structure. Further analyses should be performed in the samples core, before and after extraction, to evaluate how deep extraction of Fe/S species occurred. Moreover, contamination with iron disulfides species should be performed to evaluate if BT can be applied for any sulfide-metal extraction.
Furthermore, Raman analyses displayed bands not attributed to reduced sulfur compounds or iron oxides. The strongest bands were in the range 1200-1500 cm −1 . Investigation of these bands showed that they may be due to the wood substrate. If Raman was not considered to evaluate wood degradation, these results demonstrated that it could be employed as complementary method with ATR-FTIR spectroscopy to determine the state of degradation of wood.
PCA analyses carried out on the dataset before and after extraction did not show clusters suggesting a poor distinction among the sets. However, previous observations showed that no reduced sulfur species were identified for sets D1, D2, G1 and G2 before extraction. A plausible explanation for the clustering of all sets together may be due to the presence of the vibrational bands attributed to the wood substrate. Indeed, for all, a similar pattern was observed and assigned to the wood based on previous studies [37]. This pattern may interfere and be the main variable for the PCA analyses. However, similar analyses performed only in the range 100-600 cm −1 displayed same results with all samples clustering. It seemed that PCA was not the most convenient method to interpret this dataset. Additional analyses through Hierarchical Clustering Analysis (HCA) and spectra distance (data not shown) did not provide further information.
Deeper investigation was then performed, similar to the one applied for the ATR-FTIR dataset. The regression performed depending of the five reference compounds (i.e., partially oxidized mackinawite, sulfur, pyrite, goethite and hematite) showed that the main variable were hematite bands for Raman spectra. This result already showed that the prediction model was not well adapted for the dataset. Moreover, the application of the model to the dataset resulted in false error regarding reduced sulfur compounds. The prediction model should be optimized. As described in Section 3.2., no iron oxide compounds were observed on the surface of the samples. Their characteristic bands should not be the main variable to understand and classify the Raman dataset.
This study demonstrated that PCA-LASSO analyses allow validating our ATR-FTIR dataset. However, concerning Raman spectroscopy, these two methods did not validate the observations. Further statistical analyses such as one-way ANOVA or K nearest neighbor (KNN) analyses could be applied on the Raman dataset. Results then obtained could help to understand the ones observed so far and to provide a prediction model for further analyses.

Application to Real Artefacts
This study proves the feasibility of a biological extraction method for Fe and S species. However, as discussed in Section 4.2., the analyses were performed only on the surface of the samples and not in the core. Further investigations should be carried out to validate the efficiency of the microorganisms to extract the species within the objects. In addition, it has been reported that for the larger objects, such as the Vasa warship, iron and sulfur species were detected in the two first centimeters [43]. The proposed biological treatment would be able to reach such depth for large objects to have an optimal extraction rate if small-shaped bacteria are employed. T. denitrificans are a short rod-shaped bacteria of dimensions 0.5 × 1-3 µm [44] and the typical range of wood pores is 150-200 µm for hardwood, such as oak, [45] and 30-40 µm for softwood, such as pine [46]. They can penetrate within the wood pores and then have the capacity to oxidize iron sulfides compounds in situ.
An application on such larger objects should also be considered. In this study, the samples were cubes of 2 × 2 × 2 cm 3 . If some waterlogged archaeological wooden objects may be of similar shape and/or volume, most of them are more massive. Indeed, the main waterlogged artefacts are everyday life objects, from different periods and civilizations [47][48][49]. They range from buckets, tools to boats. For instance, over 500 boats have been excavated from Swiss lakes since 1850 [49]. In this case, a treatment per immersion requires an infrastructure adapted to such artefacts, with large immersion installation baths, as it is the case in ARC-Nucléart, Grenoble (FR) [50]. Another solution would be application by spray, as developed for the consolidation of the Vasa warship at the Mary Rose Trust, Portsmouth (UK). Another point should be the duration of the treatment. An adaptation of the time of the treatment should be considered to ascertain that the extraction of harmful species was efficient till the core of the objects.

Conclusions
With the perspective of developing an innovative green biological extraction method, ATR-FTIR and Raman spectroscopies methods were used to compare the efficiency of biological treatment over commonly used chemical extraction method. First results were promising. After extraction, reduced sulfur compounds identified by Raman spectroscopy were not identified anymore after biological treatment (BT) while sulfur was still present for fresh and lake pine after chemical treatment (CT). The wood species may have more impact on the efficiency of CT but not BT as extraction on surface seemed efficient of all BT mock-ups. However, it is worth mentioning that both treatments seemed to be quite smooth methods as no further wood degradation was observed with ATR-FTIR measurements.
The statistical approach performed on the spectroscopic dataset validated our observations regarding ATR-FTIR measurements. Different regression and prediction models should be investigated for Raman spectra, to have a better understanding from the data available. In the following months, the consolidation phase in PEG baths foreseen within the MICMAC project will be applied to assess the long-term behavior of the samples after different extraction methods. Prediction models will be then optimized to classify spectroscopic data before and after extraction, as well as after consolidation. As well, Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) is foreseen to evaluate the extraction rate of Fe and S species for each extraction method studied. Raman analyses of the samples core will also be performed to evaluate how deep samples are previously impregnated with iron and sulfur species as well as how deep the extraction is effective. Furthermore, imaging analyses (i.e., Raman and Scanning Electron Microscope imaging) of woods samples are planned before and after treatment to evaluate possible changes in the wood structure (cell collapse or swelling for example).