# Advanced mapping of inorganic treatments on porous carbonate stones by combined synchrotron radiation high lateral µXRPD and µXRF

G. Massinelli, N. Marinoni, C. Colombo, G. D. Gatta, M. Realini, M.

Burghammer, E. Possenti

► To cite this version:

G. Massinelli, N. Marinoni, C. Colombo, G. D. Gatta, M. Realini, et al.. Advanced mapping of inorganic treatments on porous carbonate stones by combined synchrotron radiation high lateral µXRPD and µXRF. Scientific Reports, 2024, 14 (1), pp.9108. 10.1038/s41598-024-58718-z. hal-05099486

HAL Id: hal-05099486

https://hal.science/hal-05099486v1

Submitted on 5 Jun 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Advanced mapping of inorganic treatments on porous carbonate stones by combined synchrotron radiation high lateral µXRPD and µXRF

G. Massinelli<sup>1, ∗</sup>, N. Marinoni<sup>1, †</sup>, C. Colombo<sup>2</sup>, G. D. Gatta<sup>1</sup>, M. Realini<sup>2</sup>, M. Burghammer<sup>3</sup> & E. Possenti<sup>2</sup>

Understanding the effects of consolidating inorganic mineral treatments on carbonate stones of cultural heritage, and on the nature and distribution of newly formed products within the matrix, poses a significant challenge in Heritage Science and Conservation Science. Existing analytical methods often fail to deliver spatial and compositional insights into the newly formed crystalline phases with the appropriate high lateral resolution. In this study, we explore the capabilities and limitations of synchrotron radiation (SR) micro-X-ray powder diffraction (µXRPD) mapping combined with micro-X-ray fluorescence (µXRF) to give insight into compounds formed following the application of ammonium oxalate (AmOx) and diammonium phosphate-based (DAP) solutions on porous carbonate stone. Ultimately, the integration of µXRPD mapping and µXRF analysis proved itself a powerful asset in providing precise qualitative and quantitative data on the newly formed phases, in the case of both calcium oxalates (CaOxs) and calcium phosphates (CaPs), and their complex stratigraphic distribution, thus opening a new route for applications to a more comprehensive study of inorganic treatments applied to carbonate substrates.

The consolidation of carbonate stones poses a significant challenge in preserving Cultural Heritage (CH) buildings and artifacts due to material decohesion caused by decay¹. Since the mid-1990s, various inorganic treatments have been developed to restore microstructural cohesion within damaged stones matrixes²⁻⁴. In such context, ammonium oxalate (AmOx) and di-ammonium phosphate (DAP) solutions have gained prominence for carbonate stones in the past two decades⁵⁻⁷. These water-soluble salts react with the carbonate matrix of the stone substrate to produce micrometric reaction products (ranging from 0.1 to 4 µm in crystal size)⁸. Specifically, AmOx primarily produces Ca-oxalate (CaOxs) while DAP's reaction should ideally yield calcium phosphate phases (CaPs). Oxalate phases exhibit limited penetration into the stone (~1 mm) whereas DAP-based treatments offer greater consolidating capabilities in terms of treatment penetration (2–8 mm, depending on the stone substrate)⁹,¹⁰.

However, inorganic mineral treatments exhibit significant variability in terms of the newly formed phases within the treated stone materials and their effect on the stone microstructure. Factors such as their crystal-chemical composition, crystallinity, penetration, and distribution are influenced by a series of interrelated variables playing during the solution reaction with the calcite matrix (i.e., solution molarity, treatment methodology, etc.)¹⁰,¹². Consequently, studying the outcome of the reaction process, specifically the new “reaction products-calcite matrix” system, presents a significant challenge. Despite this complexity, understanding the modifications induced by inorganic treatments is crucial, given that the physical, chemical, and microstructural properties of the new system significantly impact the macroscopic interaction of the building stone with the agents of decay. The complexity of studying these transformations becomes particularly pronounced in the case of reaction products that share a similar elemental composition to those of the carbonate stone substrate, like

¹Dipartimento di Scienze della Terra "Ardito Desio", Università degli Studi di Milano, Via Botticelli 23, 20133 Milan, Italy. ²Istituto di Scienze del Patrimonio Culturale (ISPC), Consiglio Nazionale delle Ricerche (CNR), Via R. Cozzi 53, 20125 Milano, Italy. ³European Synchron Radiation Facility, 71 Avenue des Martyrs, 38000 Grenoble, France. ⁴email: giulia.massinelli@unimi.it; nicoletta.marinoni@unimi.it
in the case of CaOxs formed by AmOx solutions, or in the case of complex reaction kinetics in the presence of carbonate ions that result in not stoichiometric reaction. Examples of this last scenario are CaPs phases formed with DAP solutions, which ideally should yield crystalline hydroxyapatite (HAP, Ca₁₀(PO₄)₆(OH)₂), however, the result is complex mixtures of crystalline and poorly crystalline phosphate phases¹¹.

In accordance, a range of analytical techniques has been employed tailored to characterize carbonate stone treated with inorganic mineral treatments and gain insights into the composition and penetration depth of the newly formed phases to assert their performance. Conventional methods such as SEM-EDX, XRD, and FT-IR have been widely used but lack the specificity and lateral resolution required for characterizing phases with different degrees of crystallinity, in addition to being invasive (as many of their experimental setups require a sample) or, in some cases destructive (i.e., with a complex sample preparation)¹³⁻¹⁵. High lateral resolution micro-spectroscopy techniques, like µ-ATR-FTIR and µ-Raman mapping or spectroscopic imaging, have proven valuable for probing the spatial distribution of consolidants, providing complementary structural and compositional insights on the newly formed phases¹¹,¹⁶⁻¹⁸. However, these techniques may face challenges in distinguishing the newly formed phases within the microstructure due to overlapping vibrational bands. More recently, synchrotron radiation (SR) facilities have experienced a remarkable rise in the study of Cultural Heritage materials¹⁹⁻²². Advanced radiation imaging techniques, such as SR-based X-ray micro-computed tomography (µ-CT), have been employed to non-invasively determine the effects of the new product crystallization on the 3D stone's microstructure at a maximum voxel size of 1 µm³. However, no chemical information is provided²³⁻²⁷. Moreover, micro-X-ray diffraction computed tomography (XRDCT) measurements allow structural characterization and 3D localization of the reaction products but may lack point-specific information at good spatial resolution²⁸,²⁹. Lastly, more recently, SR-based spectroscopy techniques such as Ca K-edge 2D X-ray absorption near-edge structure (XANES), with µ-XRF mapping, have been successfully employed to spatially distinguish among amorphous calcium carbonate with a micrometric lateral resolution³⁰.

Despite the advancements in analytical techniques, it is still challenging to achieve comprehensive spatially resolved insights into the “reaction products-calcite matrix” system. Synchrotron radiation (SR) micro-X-ray powder diffraction (µXRPD) mapping stands out as a powerful asset, as it can provide spatially resolved information about the crystallographic properties of the reaction products at a lateral resolution down to the submicrometric scales. This technique, when complemented with the elemental information of micro-X-ray fluorescence (µXRF) mapping, becomes particularly well-suited for the micro-analysis of highly heterogeneous systems, especially when the new phases are present in minor fractions if compared to the substrate's minerals or are poorly crystalline. These techniques, although well-known and already used in many fields over the years³¹⁻³⁵, in CH materials are mainly used in the examination of oil paintings²²⁻³⁹ and wall paintings⁴⁰,⁴¹, ceramics⁴²,⁴⁴, paper⁴⁵, and metals⁴⁶,⁴⁷, but, to the best of our knowledge, have not been employed to study the reaction products of inorganic consolidation solutions on porous carbonate stones in CH.

In this paper, by building on what others have done in this field, we're also bringing something new to the table by addressing the challenges of dealing with more complicated systems. We explore the capabilities and limitations of this methodology approach combining SR-µXRPD and µXRF mapping, in collecting comprehensive qualitative and quantitative information about the crystal-chemical compositions and spatial distribution of the newly formed phases (CaPs and CaOxs), resulting from the reaction of AmOx and DAP solutions on porous carbonate stones, to test its applicability and contribution in the ongoing knowledge related to stone conservation in Cultural Heritage.

# Materials and methods

## Sample selection and preparation

Noto Yellowish Limestone, our primary substrate, is a building stone used for Baroque monuments in the Val di Noto (south-eastern Sicily) with a predominant composition of calcite, with minor clay minerals, quartz, and iron hydroxides. This stone shows a pore size range from 1 to 10 µm and from 1 to 0.1 µm and an open porosity ranging between 25 and 35%²⁴,²⁹.

A tailored preparation method was devised to preserve crystalline phases and microstructural information. A thin section, of approximately 30 µm thickness, was prepared traversal (from the treated surface to the bulk) from each 5 × 5 × 2 cm³ treated limestone specimen. Embedding the sample in resin and adding a piece of polycarbonate to the section were essential steps to stabilize the thin section, during the follow-up cutting procedure. The thin sections were resized to a final size of 2 × 4 mm, to optimize the mounting space on the sample holders. Such preparation is compatible with the ID13 set-up allowing for optimal transmission of X-rays and a controlled and uniform probed voxel across the two-dimensional surface.

To obtain a more objective assessment of the capabilities and limitations of the techniques employed, we subjected the limestone samples to various treatment protocols. This systematic approach allows us to assess the performance of our methodologies under different conditions. Two protocols encompassed the application of just one solution (either AmOx or DAP). Another one adopted the application of a sequential treatment (DAP followed by AmOx solution). The diverse treatments aim to replicate real-world situations encountered in the conservation of CH materials. Concentrations for DAP and AmOx as well as application methodology and condition were determined based on precedent experiments and established practices from conservation worksites.

## Combined micro X-ray fluorescence (µXRF) and micro X-ray powder diffraction (µXRPD)

The combined Micro X-ray fluorescence (µXRF) and Micro X-ray Powder diffraction (µXRPD) investigations were conducted at the ID13 beamline at the European Synchrotron Radiation Facility (ESRF) in Grenoble (Fig. 1a). The samples were mounted together on a 4 mm diameter 12-hole holder, as shown in Fig. 1b, then
<figure><img src="image_3.png" /><figcaption>**Figure 1.** (a) SR-μXRPD scanning and μXRF set-up at ID13 beamline, ESRF, (b) some of the analyzed thin section (2×4 mm) mounted on the 4 mm diameter 12-hole holder, and (c) Schematic illustration of μXRPD mapping and μXRF experiment on Noto limestone thin section. (d) μXRD patterns and μXRF spectra were acquired at each pixel of an area of 400 µm (horizontal) by 400 µm (vertical) defined over the sample surface.</figcaption></figure>

mounted vertically, perpendicular to the X-ray beam. The μXRPD branch performed crystalline phase mapping using a 2.5 × 2.5 µm² beam with an energy of 13 keV. Throughout the thin section, three regions of interest (ROIs) were systematically collected, commencing from the sample surface (as illustrated in Fig. 1c). Each ROI corresponds to a map size of 400 µm (horizontal) by 400 µm (vertical) covering an area of the sample of approximately 1.2 mm. This specific ROI size was chosen to strike a balance between achieving sufficient spatial detail and maintaining manageable data acquisition times. Considering the characteristics of the limestone samples, including potential variations in crystalline phases or microstructural features, the 400 µm by 400 µm ROI was deemed apt for capturing relevant information without the pitfalls of oversampling or undersampling. Two-dimensional (2D) diffraction patterns were acquired in transmission mode at every pixel of 2D maps using the graphical user interface (GUI) Daiquiri. These X-ray diffraction patterns were subsequently transformed into one-dimensional (1D) diffractograms through azimuthal integration, facilitated by Jupyter Notebooks based on the PyFAI software package⁴⁶. Simultaneously to each diffraction map, μXRF spectra were obtained for the same map (Fig. 1d), using a Vortex EM detector mounted orthogonally to the beam.

## Data processing

Data processing was an essential step in our study as the “reaction products–stone matrix” system exhibited Ca-baring elemental composition but with diverse crystalline phases in different locations and degrees of crystallinity. To investigate this complexity, we relied on the PyMca 5.9.2 ROI imaging software (link https://sourceforge.net/projects/pymca/)⁴⁸, which enables simultaneous visualization of μXRPD and μXRF imaging data. This parallel processing workflow maximized the synergy of the two techniques, allowing us to not only visualize elemental composition but also directly identify crystalline phases and their location.

The initial data processing step involved I/I₀ normalization of the intensity map pixels to rectify variations induced by the incoming X-ray beam. Figure 2a shows an example of the false color appearance of the μXRF map before normalization, emphasizing the importance of this step to produce an accurate representation of elemental intensity (in Fig. 2b the elemental map after the normalization is shown).

Afterward, the normalized μXRF map underwent a fitting process for specific elements composing our “reaction products-calcite matrix” system. The fitting process generated distribution maps for these elements, enabling effective comparisons of their presence and distribution. An illustrative example Fig. 2c displays the elemental distribution map for calcium in one of the samples.

Regarding the μXRPD imaging, the software displayed a correlated μXRPD pattern representing the averaged profile of all acquired patterns across the 2D map. Raw data preprocessing involved 1D background correction, addressing the contributions of the amorphous polycarbonate layer in the sample preparation, which can make it more challenging to identify and analyze the diffraction peaks associated with the limestone system. This improvement in peak-to-background ratio facilitates more reliable peak identification, especially for weak or overlapping peaks, like in the case of the new crystalline phases. The exclusion of low-angle regions was done to not result in the loss of important information associated with low angles, like some amorphous or poorly
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Data processing on the PyMca ROI imaging software of one of the samples. (a) Original µXRF map of an ROI; (b) I/I₀ normalization of the µXRF map; (c) RGB correlation of the elemental distribution maps of calcium, iron, and magnesium; (d) Original µXRPD map of the same ROI and (e) µXRPD map after the 1D background correction.</figcaption></figure>

crystalline phases of CaPs and CaOxs. In Fig. 2d,e, the effect on µXRPD maps of the background correction process is shown.

In the post-preprocessing, the identification of crystalline phases within µXRPD patterns, exported from PyMca, employed the X'Pert Highscore Plus software 5.2 (link https://www.malvernpanalytical.com/en/suppo rt/product-support/software/highscore-software-update). The analysis complemented the data processing using PyMca. Leveraging the strengths of PyMca, we conducted tailored data post-processing to address the diverse analytical inquiries arising in the study of the effects of inorganic treatments on carbonate stones. The subsequent post-processing methodologies will be described in the upcoming sections.

PyMca’s “RGB Correlator” tool was mainly used as it visually highlights distribution patterns of elements alongside crystalline phases, creating µXRPD and µXRF correlation maps. This capability proved indispensable in unraveling relationships between element distribution and identifying correlations or anticorrelations in calcium-based crystalline phases within a predominantly calcite matrix.

Finally, the quantification of XRD data involved image analysis using Fiji, Open-Source software ImageJ (OS version, link https://imagej.net/ij/), applied to diffraction maps corresponding to each identified phase. Utilizing the threshold tool, the software extracted the areal percentage of each phase within the specified region of interest. This process provided a quantitative distribution of each phase. To ensure the robustness of our findings, average percentages and their standard deviations were calculated for each treatment. Also, to mitigate subjectivity associated with threshold tools, small squares were inserted inside the areas designated for segmentation, and all the data were plotted and described by a function. The data is well-fitted with a Gaussian function. The limits of this Gaussian function were employed during the segmentation. Furthermore, all images processed with this method were acquired under the same experimental conditions and saved with consistent dimensions and resolutions to maintain constant image quality.

**Results and discussion**

**Qualitative phase characterization of the "reaction products-calcite matrix" system**

Initially, we explored the combined µXRPD and µXRF mapping capabilities to characterize the different reaction products (CaOxs and CaPs) within the newly formed system for all treatment methodologies (single treatment and sequential treatments).

To begin, we analyzed the average µXRPD patterns for all Noto Limestone samples, both untreated and treated. The average diffraction patterns of the untreated limestone samples show predominantly the peaks associated with calcite, with occasional quartz. This trend remains consistent across treated samples' diffraction patterns, as illustrated in Fig. 3. In these patterns, the main peaks of the reaction products were weak. This result aligns with prior literature and reflects how the reaction products of both solutions are in a low wt% fraction if compared to the abundant calcite matrix¹⁷.

Focusing specifically on the reaction products from the AmOx solution (Fig. 4a and b), the CaOxs, low-intensity peaks appear in some surface-proximate diffraction patterns (as shown in Fig. 4c), where the highest crystallization activity for these phases usually occurs, suggesting the presence of substances, but not a
<figure><img src="image_5.png" /><figcaption>**Figure 3.** (a) Average µXRPD patterns extracted from representative surface maps of an untreated Noto limestone sample (NT) and three samples treated by AmOx, DAP, and DAP + AmOx; (b) A zoom (position ~ 8–20 2Theta) to highlight the marker peaks of calcite (CAL) and quartz (QTZ).</figcaption></figure>

conclusively identification. The challenge in their identification arises due to the nano-crystalline nature of the main CaOx phases—whewellite (CaOx monohydrate) and weddellite (CaOx dihydrate)—together with their low fractions if compared to the calcite matrix⁴⁴. Even with our high-resolution XRPD patterns, discerning these CaOxs in all patterns remains difficult.

The sole examination of µXRF imaging does not facilitate their identification either, as CaOxs share a similar elemental composition to the calcite matrix (Ca, C, O). However, in the case of the AmOx-treated samples, the high-resolution µXRF Ca maps have shown, as seen in Fig. 4d, low-intensity (~ 70–100 in the scale bar) and high-intensity areas (corresponding to ~ 170–250 in the scale bar). An unusually big gap in the intensity of Ca through the elemental map was here registered. These intensity variations throughout the Ca maps could, therefore, be correlated with different calcium contents in crystalline phases within our “reaction products-calcite matrix” system. Given that CaOxs contain less calcium (around 38% CaO) than calcite (approximately 55% CaO)⁵², these intensity fluctuations in the elemental map could suggest a different distribution of oxalate and calcite phases within the stone system, corresponding respectively to the low and high-intensity area.

Combining the µXRF maps with µXRPD data, we could generate *area-specific* average diffraction patterns by selecting pixels from areas exhibiting these calcium intensity variations. These patterns reveal peaks attributed to CaOxs with good resolution, as shown in Fig. 4e. The *area-specific* diffraction pattern exhibited a high signal-to-noise ratio and a consistently flat baseline, which helps us confirm not only the crystallization but also facilitates the unambiguous identification of whewellite and weddellite.

The situation becomes even more complex in the case of DAP-treatment samples (Fig. 4f and g). Like CaOxs, CaPs are nanometric-sized products and low in concentration⁴⁹. Additionally, the complex reaction of DAP solutions with calcite does not yield only hydroxyapatite (HAP) but results in a mixture of by-products with varying Ca/P ratios in concentration distributed within the matrix based on reaction conditions⁵⁰. These by-product phases share similarities in elemental composition and crystalline structure⁵¹, making their phase identification even more complex. Consequently, map average diffraction patterns of DAP-treated samples lack identifiable peaks even in our diffraction patterns, as shown in Fig. 4h, which hinders the confirmation of the presence of CaPs.

Moving on to µXRF imaging, a situation like the AmOx-treated sample was observed in the elemental maps of calcium distribution in samples treated with DAP. This is coupled with localized high phosphorus intensities (Fig. 4i). Given the known high Ca/P ratio in HAP (around 1.76), correlating calcium intensity variations with areas of high phosphorus concentration implies the potential CaP presence. Confirming this involves associating the µXRF imaging with the µXRPD data, generating *area-specific* diffraction patterns for pixels showing Ca/P elemental correlation. By generating specific diffraction patterns for these correlated areas, the identification of CaPs was achievable. The resulting diffraction patterns showed well-resolved HAP peaks that supported its visualization and identification. These patterns exhibit also good resolution at the low angular range where by-products of the reaction reside. The result is that this combined technique aids in identifying traditionally difficult-to-characterize²⁰ by-products such as octacalcium phosphate (OCP) and other CaP phases associated with HAP crystallization, such as Ammonium Dihydrogen Phosphate ADP, and carbonate-substituted HAP (C-HAP), as shown in Fig. 4l, providing crucial information about the composition of the DAP “reaction products-calcite matrix” system.

The combination of techniques proves effective for qualitative phase analysis of both CaOxs (whewellite and weddellite) and CaP (HAP and by-products), even in the more complex scenario of sequential treatment samples (DAP + AmOx). Despite the additional challenges posed here because of the simultaneous presence of CaOxs and CaPs, this approach facilitated the visualization and characterization of the newly crystallized phases (oxalate and phosphate) within the complex systems generated by these treatments, as shown in Fig. 5.

In the end, focusing solely on the average µXRPD pattern, the identification of reaction products was limited mainly because of their low concentration compared to the dominant calcite matrix. In contrast, examining µXRF data alone resulted challenging in samples lacking elemental markers for identifying CaOxs and could only hint
<figure><img src="image_6.png" /><figcaption>DAP-treated Noto Limestone sample</figcaption></figure>

<figure><img src="image_7.png" /></figure>

**Figure 4.** Comparison between single treated samples. For the AmOx treated: (a) Optical image of the cross-section. Selection of the investigated ROI (green) and the investigated pixels (red); (b) Representative average surface µXRF map of the ROI; (c) Map average µXRPD pattern extracted from the ROI; (d) Elemental map of calcium showing variation in intensity; (e) Area-specific average µXRPD pattern for the red area exhibiting low Ca intensity. A zoom of the diffraction patterns (position ~8–20 2Theta) is shown to highlight the marker peaks of whewellite (WHE) and weddellite (WED). For the DAP treated: (f) Optical image of the cross-section. Selection of the investigated ROI (green) and the investigated pixels (red); (g) Representative surface µXRF map of the ROI; (h) Map average µXRPD pattern extracted from the ROI (i) RGB correlation map of the elemental distribution of calcium (blue) and phosphorous (red); (l) Area-specific average µXRPD pattern for the red area exhibiting correlation Ca/P. The patterns are presented in a zoom (position ~16–25 2Theta) to highlight the marker peaks hydroxyapatite (HAP), Ammonium Hydrogen Phosphate (ADP), and carbonate-substituted HAP (C-HAP).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** DAP + AmOx treated Noto Limestone sample (a) Optical image of the cross-section. Selection of the investigated ROI (blue) and the investigated pixels (red); (b) Area-specific average µXRPD pattern for the red area and a zoom (position ~8–21 2Theta) to highlight the marker peaks of hydroxyapatite (HAP), whewellite (WHE), and weddellite (WED).</figcaption></figure>

at the crystallization of CaPs when the markers were present. Integrating both techniques enabled the production of *area-specific* diffraction patterns associated with those zones exhibiting calcium intensity variations or high Ca/P correlation, allowing for the identification of specific crystalline phases characteristic of the newly formed “reaction products-calcite matrix” system.

## Phases' penetration depth and spatial distribution of the reaction products

The combination of techniques used in this study yields an invaluable advantage in generating comprehensive µXRPD maps. *Phase-specific* distribution maps were generated by isolating distinct peaks associated with the identified crystalline phases in the diffraction patterns, providing detailed insight into the spatial distribution of each phase. SR-µXRPD mapping allows micrometric individual mapping of each reaction product distribution within the “reaction products-calcite matrix” systems, offering good spatial resolution. Such data not only visualizes the penetration depth of the crystallization but also localizes where the chemical reactions take place.

Our analysis consistently reveals *phase-specific* color maps accurately depicting the areal distribution of CaOxs and CaPs crystalline phases across all examined samples, without artifacts.

In the context of AmOx-treated samples, the areal distribution data through RGB correlation maps reveals a distinct distribution pattern of CaOx phases: a gradient from the surface, with weddellite (WED) overlaying whewellite (WHE). This observation underscores the spatial relationship between these phases shown in Fig. 6a. The data underscored a crystalline network following the solution diffusion path within the stone's microstructure, originating from the surface interface. These findings highlight the CaOx areal arrangement given by the treatment. Such data offer, for the first time, an opportunity to establish correlations between the phases of CaOx formed and their spatial distribution, allowing interpretations of the consolidation effect. Currently, literature only suggests that the distribution and penetration of distinct CaOxs significantly affect resultant system properties⁴⁹.

The application of µXRPD mapping of DAP reaction products has revealed unprecedented details regarding the areal distribution of CaP phases within the stone, particularly highlighting hydroxyapatite (HAP) and its associated by-products (ADP, OCP, and C-HAP). As shown in Fig. 6b, HAP distribution follows a gradient from the surface, less homogeneous than the CaOxs network. The by-products, less identified in prior studies, exhibit a distribution pattern linked to HAP, extending to similar depths within the stone and localized in spots within its porous network. It is important to acknowledge that the representation of these CaP phases in the color maps might be subject to potential inaccuracies (either underestimated or overestimated) because of the mixture of phases within our systems. CaP peaks are collocated in the diffraction pattern close to the peaks of other phases complicating their isolation and mapping. Despite this, the ability to pinpoint the spatial distribution of both HAP and its by-products remains crucial. The literature already suggests correlations between the presence of these by-products, such as ADP indicating acidification, and changes in the stone's condition¹⁰. Therefore, having data on CaP distribution allows for a deeper exploration of these correlations, facilitating a better understanding of the localized effects of CaP phases on the stone's chemical changes and conditions.

In the case of sequential treatment, the mapping of CaOxs and CaP phases, and the examination of phase coexistence through RGB correlation maps, highlight distinct variations in both penetration depth and areal distribution compared to single treatments. The data (Fig. 6c) showed that, in sequential treatments, HAP and by-products are dispersed throughout the stone bulk rather than follow a gradient from the stone surface. CaOxs seemly exhibit different penetration depths compared to single treatments, although they are still primarily located on the surface. Additionally, the maps illustrate sparse areas of coexistence between CaOxs
<figure><img src="image_9.png" /><figcaption>**Figure 6.** (a) RGB correlation of the µXRPD areal distribution maps of whewellite (WHE) and weddellite (WED) in an AmOx-treated sample; (b) RGB correlation of the µXRPD areal distribution maps of the marker peaks of identified HAP and by-products (ADP and OCP) in a DAP-treated sample; (c) RGB correlation of the µXRPD areal distribution maps the main CaOxs (whewellite and weddellite) and CaPs (hydroxyapatite) in a DAP + AmOx treated sample. All three figures are an ROI of 400 × 400 µm² taken from the surface of the samples.</figcaption></figure>

and CaPs. These observations strongly indicate interactions between the solutions' reactions. The significance of obtaining these data lies in its capacity to shed light on these interactions, thereby enhancing our understanding of combining treatment effects on stone microstructures absent in the literature.

## Quantitative phase analysis and orientation of the reaction products

The high resolution of µXRPD maps also enables us to obtain information on quantitative phase concentration and orientation at the scale of the studied area (ROI = 400 × 400 µm²). Acquiring such data at this scale, with a clear view of the phases' distribution, provides us with essential insights for assessing localized effects and evaluating treatment performance.

For quantitative assessment of newly crystallized reaction products, image analysis of phase-specific µXRPD maps was used to quantify the areal fraction of reaction product observed in the 2D map (as explained in the Methods paragraph). The % ± Standard Deviation results are shown in Table 1. In the case of AmOx-treated samples, this quantification confirmed existing literature by highlighting predominant whewellite crystallization, while uncovering unexpectedly high content of weddellite in some surface areas. Similarly, in samples treated with DAP solution, through quantitative evaluation, we categorized the main reaction products (HAP and OCP, 22.2% and 7.6%, respectively) as well as reaction by-products into minor (between 1% and 0.1%) and trace (less than 0.1%) products, based on their concentrations. The localized quantitative data derived from µXRPD map image analysis, of the sequentially treated samples, revealed that the % of crystallization of the reaction compounds is, generally, in line with the one in the single treatments.

While generating µXRPD distribution maps by selecting Bragg peaks corresponding to the crystalline phase, we observed variations among maps for different diffraction peaks. Each peak in the diffraction pattern represents specific crystallographic planes, mapping the areal distribution of distinct crystal planes. RGB correlation of maps associated with different peaks of a single phase was employed to reveal the crystal orientation of the reaction products. The lack of color merging in specific areas in RGB correlation maps indicated the crystal orientation of selected crystalline planes. For instance, in the case of CaOxs (both whewellite and weddellite), they exhibit randomly distributed crystallites within the stone matrix, showing up as a single cohesive yellow color on the maps (Fig. 7a). On the other hand, CaPs, primarily HAP, in DAP-treated samples show preferential orientations in some areas, as distinct green (HAP $d_{310}$) and red (HAP $d_{002}$) color distributions show up on the correlation maps (Fig. 7b). This observation becomes particularly intriguing when studying sequential treatments where no significant orientation of CaPs can be detected (Fig. 7c). Here the data show how CaPs crystallized differently.

<table><caption>**Table 1.** Quantitative assessment based on image analysis of the reaction products areal fraction (Mean % ± Standard deviation) on surface maps for various treatment methodologies.</caption><thead><tr><th rowspan="2"></th><th colspan="6">Reaction products</th></tr><tr><th>WHE</th><th>WED</th><th>HAP</th><th>C-HAP</th><th>OCP</th><th>ADP</th></tr></thead><tbody><tr><th>AmOx-treated</th><td>25 ± 2.7</td><td>1.5 ± 0.5</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><th>DAP-treated</th><td>-</td><td>-</td><td>22 ± 3.9</td><td>1.5 ± 0.6</td><td>7.6 ± 1.4</td><td>0.3 ± 0.1</td></tr><tr><th>DAP + AmOx treated</th><td>21 ± 3.5</td><td>2.8 ± 0.6</td><td>16 ± 2.2</td><td>1.2 ± 0.4</td><td>5.4 ± 1.7</td><td>0.5 ± 0.2</td></tr></tbody></table>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** (a) RGB correlation of the µXRPD areal distribution maps of two marker peaks of whewellite in an AmOx-treated sample; (b) RGB correlation of the µXRPD areal distribution maps of two marker peaks of hydroxyapatite in a DAP-treated sample; (c) RGB correlation of the µXRPD areal distribution maps of two marker peaks of hydroxyapatite in a DAP + AmOx treated sample.</figcaption></figure>

This proves an interaction between the solutions (the DAP-stone system is disturbed by the crystallization of CaOxs,) and the reactivity of the CaPs phases to be disturbed by the interaction.

Overall, our experimental setup prioritized high spatial resolution and simultaneous µXRF and µXRPD acquisition, but this compromised the features of the collected X-ray diffraction patterns hindering, for example, a proper Rietveld analysis. The preclusion of Rietveld's full-profile fit hinders any structural analysis of the crystalline phases, along with a phase quantification (in a multi-phase system) based on the refined scale factor of each crystalline component. However, this limitation is here efficiently overcome, as 1) the crystal structure of the co-present crystalline phases, in the system under investigation, is already well known and 2) the relative (areal) fraction of the crystalline components can be obtained through the µXRPD map analysis, along with their average crystallites orientation. In this light, the lack of the Rietveld full-profile fit does not diminish the description of the complex polycrystalline system under investigation and opens a new route for other scenarios in which the application of full-profile fit analysis is precluded.

## Conclusion

Our research effectively demonstrates how the synergy of µXRF and µXRPD at high resolution leads to an efficient characterization of the newly formed conservation products (CaOxs and CaPs), overcoming the analytical challenge of their low wt% concentration compared to the dominant calcite matrix. While µXRPD alone struggled, due to their low abundance, µXRF provided elemental data, leading to *area-specific* diffraction patterns, and ultimately enabling phase identification in complex multiphase mixtures.

Moreover, these techniques provide invaluable insights into the penetration depth and spatial distribution of reaction products, enhancing our understanding of the treated stone system. The data derived from µXRPD mapping proved, for the first time, to be crucial in assessing reaction product concentrations. Additionally, orientation analysis of crystalline phases on the µXRPD maps offered a deeper understanding of phase formation and preferential orientations within the treated samples, and of the phase interaction in the case of sequential treatment.

While acknowledging some limitations, mainly due to the need for micro-sampling, the combination of the techniques provided high-quality data on CH stone consolidation. Moreover, the obtained positive outcomes not only contribute to our fundamental understanding of these complex systems and inorganic mineral treatments but also underline their potential applications in further studies.

In conclusion, this study underscores the significance of combining advanced synchrotron-based X-rays 2D mapping techniques for comprehensive and spatially resolved analyses of complex CH multiphase systems, which broadens the applicability of these advanced synchrotron-based approaches in the study of inorganic mineral treatments in Conservation Science and Heritage Science.

## Data availability

The data used and analyzed during the current study are available from the corresponding author upon reasonable request.

Received: 2 February 2024; Accepted: 2 April 2024

Published online: 20 April 2024

## References

1. Sena da Fonseca, B. Current trends in stone consolidation research: An overview and discussion. *Buildings* **13**(2), 403 (2023).

2. Giorgi, R., Baglioni, M., Berti, D. & Baglioni, P. New methodologies for the conservation of cultural heritage: Micellar solutions, microemulsions, and hydroxide nanoparticles. *Acc. Chem. Res.* **43**, 695-704 (2010).
3. Hansen, E. et al. A review of selected inorganic consolidants and protective treatments for porous calcareous materials. *Stud. Conserv.* **48**, 13–25 (2003).

4. Matteini, M. Inorganic treatments for the consolidation and protection of stone artifacts and mural paintings. *Conserv. Sci. Cult. Herit.* **8**, 13–27 (2008).

5. Celik, S. E., Gulen, J. & Viles, H. A. Evaluating the effectiveness of DAP as a consolidant on Turkish building stones. *Constr. Build. Mater.* **262**, 120765 (2020).

6. Menningen, J., Sassoni, E., Sobott, R. & Siegesmund, S. Constraints of the durability of inorganic and organic consolidants for marble. *Environ. Earth Sci.* **80**, 370 (2021).

7. Sassoni, E., Delhomme, C., Forst, S., Graziani, G. & He'nin, J., Masi, G., Palazzo, A., Rolland, O., and Verge's-Belmin, V., Phosphate treatments for stone conservation: a 3-year field study in the Royal Palace of Versailles (France). *Mater. Struct.* **54**, 140 (2021).

8. Possenti, E. et al. Diammonium hydrogen phosphate for the consolidation of building materials. Investigation of newly-formed calcium phosphates. *Constr. Build. Mater.* **195**, 557–563 (2019).

9. Matteini, M., Rescic, S., Fratini, F. & Botticelli, G. Ammonium phosphates as consolidating agents for carbonatic stone materials used in architecture and cultural heritage: Preliminary research. *Int. J. Architect. Herit.* **5**, 717–736 (2011).

10. Possenti, E. et al. New insight on the interaction of diammonium hydrogen phosphate conservation treatment with carbonatic substrates: A multi-analytical approach. *Microchem. J.* **127**, 79–86 (2016).

11. Possenti, E., Colombo, C., Realini, M., Song, C. L. & Kazarian, S. G. Time-resolved ATR-FTIR spectroscopy and macro-ATR-FTIR spectroscopic imaging of inorganic treatments for stone conservation. *Anal. Chem.* **93**(44), 14635–14642 (2021).

12. Conti, C. et al. Phase transformation of calcium oxalate dihydrate-monohydrate: Effects of relative humidity and new spectroscopic data. *Spectrochim. Acta Part A Mol. Biomol. Spectrosc.* **128**, 413–419 (2014).

13. Sassoni, E., Graziani, G., Franzoni, E. & Scherer, G. W. Calcium phosphate coatings for marble conservation: Influence of ethanol and isopropanol addition to the precipitation medium on the coating microstructure and performance. *Corros. Sci.* **136**, 255–267 (2018).

14. Sassoni, E., Ugolotti, G. & Pagani, M. Nanolime, nanosilica, or ammonium phosphate? Laboratory and field study on consolidation of a byzantine marble sarcophagus. *Constr. Build. Mater.* **262**, 120784 (2020).

15. Possenti, E. et al. Grazing incidence synchrotron X-ray diffraction of marbles consolidated with diammonium hydrogen phosphate treatments: Non-destructive probing of buried minerals. *Appl. Phys. A* **124**, 383 (2018).

16. Osticoli, I. et al. Micro-Raman analysis on the combined use of ammonium oxalate and ammonium phosphate for the consolidation and protection of carbonate stone artifacts. *J. Raman Spectrosc.* **48**(7), 966–971 (2017).

17. Conti, C. et al. Investigation of calcium and magnesium phosphate crystals in stones treated with diammonium hydrogen phosphate conservation product: potential of micro-raman spectroscopy. *Crystals* **13**(8), 1212 (2023).

18. Calore, N. et al. High resolution ATR µ-FTIR to map the diffusion of conservation treatments applied to painted plasters. *Vib. Spectrosc.* **98**, 105–110 (2018).

19. Cotte, M., Genty-Vincent, A., Janssens, K. & Susini, J. Applications of synchrotron X-ray nano-probes in the field of cultural heritage. *Cr. Phys.* **19**, 575–588 (2018).

20. Bertrand, L. et al. Development and trends in synchrotron studies of ancient and historical materials. *Phys. Rep.* **519**(2), 51–96 (2012).

21. Janssens, K., Van der Snickt, G., Vanmeert, F., Legrand, S., Nuyts, G., Alfeld, M. & De Wael, K. Non-invasive and non-destructive examination of artistic pigments, paints, and paintings by means of X-ray methods. *Anal. Chem. Cult. Herit.* 77–128 (2017).

22. Cotte, M. et al. New opportunities offered by the ESRF to the cultural and natural heritage communities. *Synchrotron Radiat. News* **35**(5), 3–9 (2022).

23. Slavíková, M. et al. X-ray radiography and tomography for monitoring the penetration depth of consolidants in Opuka—the building stone of Prague monuments. *J. Cult. Herit.* **13**(4), 357–364 (2012).

24. Massinelli, G. et al. 4D imaging and quantification of the effect of a phosphate-based treatment on building materials using synchrotron X-ray Microtomography: New insights into the dynamic of the consolidation process. *Constr. Build. Mater.* **397**, 132348 (2023).

25. Conti, C. et al. Investigation of ammonium oxalate diffusion in carbonatic substrates by neutron tomography. *J. Cult. Herit.* **19**, 463–466 (2016).

26. Realini, M. et al. Development of neutron imaging quantitative data treatment to assess conservation products in cultural heritage. *Anal. Bioanal. Chem.* **409**, 6133–6139 (2017).

27. Marinoni, N., Voltolini, M., Mancini, L. & Cella, F. Influence of aggregate mineralogy on alkali-silica reaction studied by X-ray powder diffraction and imaging techniques. *J. Mater. Sci.* **47**, 2845–2855 (2012).

28. Possenti, E., Conti, C., Gatta, G. D., Marinoni, N., Merlini, M., Realini, M. & Colombo, C. Synchrotron X-ray diffraction computed tomography to non-destructively study inorganic treatments for stone conservation. *Iscience*, **25**(10) (2022).

29. Possenti, E., Marinoni, N., Conti, C., Realini, M., Vaughan, G. B., & Colombo, C. (2024). Synchrotron radiation X-ray diffraction computed tomography (XRDCT): A new tool in cultural heritage and stone conservation for 3D non-destructive probing and phase analysis of inorganic re-treatments. *Analyst*, **139**(1), 143–152 (2014).

30. Monico, L. et al. Synchrotron radiation Ca K-edge 2D-XANES spectroscopy for studying the stratigraphic distribution of calcium-based consolidants applied in limestones. *Sci. Rep.* **10**(1), 14337 (2020).

31. Seidel, R. et al. Mapping fibre orientation in complex-shaped biological systems with micrometer resolution by scanning X-ray microdiffraction. *Micron* **39**(2), 198–205 (2008).

32. Lemke, H., Burghammer, M., Flot, D., Rössle, M. & Riekel, C. Structural processes during starch granule hydration by synchrotron radiation microdiffraction. *Biomacromolecules* **5**(4), 1316–1324 (2004).

33. Riekel, C., Burghammer, M. & Schertler, G. Protein crystallography microdiffraction. *Curr. Opin. Struct. Biol.* **15**(5), 556–562 (2005).

34. Turyanskaya, A. et al. Detection and imaging of gadolinium accumulation in human bone tissue by micro-and submicron-XRF. *Sci. Rep.* **10**(1), 6301 (2020).

35. De Samber, B. et al. Three-dimensional elemental imaging by means of synchrotron radiation micro-XRF: Developments and applications in environmental chemistry. *Anal. Bioanal. Chem.* **390**, 267–271 (2008).

36. Ghirardello, M. et al. Application of synchrotron radiation-based micro-analysis on cadmium yellows in Pablo Picasso's femme. *Microsc. Microanal.* **28**(5), 1504–1513 (2022).

37. de Mecquenem, C. et al. A multimodal study of smalt preservation and degradation on the painting “Woman doing a Libation or Artemisia” from an anonymous painter of the Fontainebleau School. *Eur. Phys. J. Plus* **138**(2), 1–8 (2023).

38. Gonzalez, V. et al. Lead (II) formate in rembrandt's night watch: Detection and distribution from the macro, Äto the micro, Äëscale. *Angew. Chem.* **135**(16), e02216478 (2023).

39. Avranovich Clerici, E. et al. Multi-scale X-ray imaging of the pigment discoloration processes triggered by chlorine compounds in the upper Basilica of Saint Francis of Assisi. *Molecules* **28**(16), 6106 (2023).

40. Oriols, N. et al. Carbonation of fresco mural paintings with a dolomitic mortar. *Cem. Concr. Res.* **157**, 106828 (2022).

41. Sciau, P., Goudeau, P., Tamura, N. & Dooryhee, E. Micro scanning X-ray diffraction study of Gallo-Roman Terra Sigillata ceramics. *Appl. Phys. A* **83**, 219–224 (2006).
42. Leon, Y. et al. The nature of marbled Terra Sigillata slips a combined µXRF and µXRD investigation. *Appl. Phys. A* **99**, 419–425 (2010).

43. Wang, T. et al. Synchrotron radiation-based multi-analytical approach for studying underglaze color: The microstructure of Chinese Qinghua blue decors (Ming dynasty). *Anal. Chim. Acta* **928**, 20–31 (2016).

44. Guarneri, N. et al. Imaging and micro-invasive analyses of black stains on the passepartout of Codex Atlanticus Folio 843 by Leonardo da Vinci. *Sci. Rep.* **13**(1), 4902 (2023).

45. Kergourlay, F. et al. Stabilization treatment of cultural heritage artifacts: In situ monitoring of marine iron objects dechlorinated in alkali solution. *Corros. Sci.* **132**, 21–34 (2018).

46. Grousset, S. et al. In situ monitoring of corrosion processes by coupled micro-XRF/micro-XRD mapping to understand the degradation mechanisms of reinforcing bars in hydraulic binders from historic monuments. *J. Anal. At. Spectrom.* **30**(3), 721–729 (2015).

47. Ševčík, R. et al. Synthetic calcium carbonate improves the effectiveness of treatments with nanolime to contrast decay in highly porous limestone. *Sci. Rep.* **9**(1), 15278 (2019).

48. Kieffer, J., Valls, V., Blanc, N. & Hennig, C. New tools for calibrating diffraction setups. *J. Synchrotron Radiat.* **27**, 558–566 (2020).

49. Possenti, E., Conti, C., Gatta, G. D., Realini, M. & Colombo, C. Diammonium hydrogenphosphate treatment on dolostone: The role of Mg in the crystallization process. *Coatings* **9**(3), 169 (2019).

50. Suzuki, O., Shiwaku, Y. & Hamai, R. Octacalcium phosphate bone substitute materials: Comparison between properties of biomaterials and other calcium phosphate materials. *Dent. Mater. J.* **39**(2), 187–199 (2020).

51. Dorozhkin, S. V. Calcium orthophosphates. *J. Mater. Sci.* **42**(4), 1061–1095 (2007).

52. Conti, C. et al. Synthesis of calcium oxalate trihydrate: New data by vibrational spectroscopy and synchrotron X-ray diffraction. *Spectrochim. Acta Part A Mol. Biomol. Spectrosc.* **150**, 721–730 (2015).

## Acknowledgements

For the beamtime received, thanks are expressed to the European Synchrotron Radiation Facility ID13 beamline (https://doi.org/https://doi.org/10.15151/ESRF-ES-981402187). The authors also acknowledge the TS Lab & Geoservice snc (Pisa, Italy) for developing an ad hoc effective protocol to prepare polished thin sections suitable for the investigations. The authors acknowledge support from the University of Milan through the APC initiative.

## Author contributions

All authors reviewed the manuscript. G.M., E.P., and N.M. wrote the main manuscript text. G.M. and M.B. performed the analysis and G.M., E.P., N.M., and C.C. interpreted the results. Preparation of Figures: N.M. and G.M. All authors (G.M., E.P., C.C., D.G., M.R., M.B., N.M.) assisted in the revision of the submitted manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to G.M. or N.M.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024