EDX/XRD-based identification of micrometer-sized domains in scanning electron micrographs of fired clay

The thermal and mechanical properties of bricks are strongly dependent on both the chemical composition and the microstructural features of the used fired clay material. Focussing on the latter, we here identify, in terms of volume fraction, shape, and orientation characteristics, one-to-several micrometer-sized subdomains (“material phases”) within the SEM-imaged microstructure of two raw clays fired at 880 and 1100 centigrades: (1) quartz grains, (2) muscovite, (3) Fe–Mg mica, (4) feldspar grains, (5) decarbonated dolomite, (6) pores, or (7) binding matrix. This identification rests on the simultaneous use of Scanning Electron Microscopy (SEM) and Energy-dispersive X-ray spectroscopy (EDX), with correspondingly obtained data entering statistical analyses based on the Otsu algorithm, and complemented by minimum grain size and grain shape requirements, as well as by logical exclusion criteria. Crystalline and amorphous phase shares were additionally confirmed by X-ray powder diffraction measurements (PXRD). As for the investigated clays, an increased firing temperature results in dehydroxylation of muscovite, and in a reduced appearance of feldspar grains.


Introduction
Clays and clay minerals may provide a wide range of chemical and physical properties [1][2][3], which renders them suitable for many different applications. Of particular interest in civil engineering and architecture are the thermal conductivity and the mechanical properties of bricks and clay blocks, as they define their sustainability and competitiveness, in Electronic supplementary material The online version of this article (https://doi.org/10.1617/s11527-020-01531-7) contains supplementary material, which is available to authorized users.
comparison to other construction materials [4]. During production, the composition of the raw clay and the firing temperature govern the resulting mineralogy and microstructure, which in turn, dictate the thermal and mechanical properties of the fired ceramic body [5][6][7].
This has motivated the characterization of the chemical and mineralogical composition of fired clay by means of X-ray fluorescence spectrometry and X-ray diffractometry [8]. In addition, the elemental mappings gained by SEM-EDX measurements allow for a distinction between the main phases preliminary observed in the XRD, a method which has been used successfully for concrete, precisely blended-cement pastes [9][10][11][12]. A particularly comprehensive multitechnique experimental characterisation of fired clay was provided by Krakowiak et al. [13], who presented XRD spectra, typical EDX composite maps, scanning electron micrographs, mercury intrusion-based pore size distributions, and massive grid nanoindentation data for two types of solid brick. Based on these measurements, Krakowiak et al. [13] propose a hierarchical scheme for the multi-level and multi-component morphology of fired clay. In this context, a material volume of some hundred micrometers characteristic length (called ''primary brick'') comprises a ''glassy'' matrix with ''silt inclusions'' and ''micropores''. This observation scale is the very focus of the present paper, where we aim at a both more detailed and more quantitative description of the aforementioned ''primary brick''. More precisely, we introduce and identify one-to-several micrometer-sized (pseudohomogeneous) subdomains within a piece of fired clay, as seen in scanning electron microscopy. Such subdomains are called ''material phases'' in the framework of continuum micromechanics or random homogenization theory [14][15][16]. These material phases govern, through their intrinsic ''universal'' properties as well as through their shapes and volume fractions, the thermal and/or mechanical properties of the aforementioned piece of material. This material phase concept has been very successful in deciphering (micro-)structure -(mechanical) property relations in other civil and bioengineering materials, such as concrete [17,18], wood [19,20], or bone [21][22][23].
Accordingly, we are particularly interested in the volume fractions of these subdomains within a piece of material, and in their geometrical description in terms of aspect ratios and spatial orientations of spheroids approximating the ''material phases''. Therefore, we here present a novel analysis method combining SEM, EDX, and PXRD data.
The remainder of the paper is organized as follows: after introducing the basic concept in Sect. 2.1, the mineralogical composition of the two raw clays is described in Sect. 2.2. This is followed by a detailed information on clay sample preparation for the SEM-EDX analysis (Sect. 2.3). A short discussion on the analysis of the mineralogy of the fired clay samples determined by PXRD measurements is given in Sect. 2.4, providing, in combination with XRF data and the results from the SEM-EDX image analysis, a detailed insight in the chemical composition of the clay matrix. Details on the SEM-EDX analysis are presented in Sect. 2.5, followed by a description of the procedure for a statistical evaluation of the SEM-EDX data, as outlined in Sect. 2.6. The evaluation of threedimensional morphological properties from two-dimensional geometrical data is comprised in Sect. 2.7. Finally, the results are summed up in Sect. 3, and discussed thereafter (Sect. 4).

Basic concept
The aim of the present paper is the identification of the morphology and the mutual position of the material phases which govern the macroscopic thermal and mechanical behaviour of the clay, at the scale of an electron micrograph with a resolution of 195 nm, and a characteristic length of several hundred micrometers. The targeted material phases are: quartz grains, feldspar grains, Fe-Mg mica, muscovite, decarbonated dolomite, pores, and a contiguous amorphous matrix with finely dispersed quartz and feldspar particle inclusions. Quantitative access to these phases stems from careful combination of data from backscattered scanning electron microscopy (SEM), energy dispersive X-ray spectroscopy (EDX), and X-ray powder diffraction (PXRD): Namely, backscattered scanning electron micrographs allow for the identification of decarbonated dolomite, pores, and the aforementioned amorphous matrix, as well as of mineral grains made of quartz or feldspar, but not for the discrimination between the two latter minerals. Similarly, SEM images allow for the discrimination of the needle-type morphology of Fe-Mg mica and muscovite, from that of all other phases, but not for the discrimination between Fe-Mg mica and muscovite. However, quartz and feldspar can be discriminated through EDX analyses, as the presence of silicon indicates quartz; the same is true for Fe-Mg mica and muscovite, as the former is suggested by the presence of iron. The comparison of the phase volume fractions as obtained from SEM-EDX, with results from X-ray powder diffraction allows for identification of the shares of quartz and feldspar occurring in grain type and in small particle type, respectively, as well as for the determination of the chemical composition of the amorphous clay matrix. In addition to the phase volume fractions, our quantitative analysis also comprises phase shape and orientation. Therefore, all phases except the amorphous matrix are approximated as spheroids with normals which are axisymmetrically distributed around the anisotropy direction of the transversely isotropic fired clay material. These spheroids appear as ellipses in planes through the anisotropy direction, and the shape and orientation of the latter allows us to re-construct phase-specific aspect ratios orientation distributions of the spheroids.

Mineralogical analysis of two raw clay types
Two raw clays of different origin, currently used for brick production, have been investigated in course of the present paper (thereafter named Clay A and Clay B). Clay A stems from a young tectonic burial and sedimentary basin. It has been deposited in Neogene (Pannonian) in a shallow basin in brackish water environment. Due to the low salt content, fossils from shells (Congeria and Limnocardium) and snails (Melanopsis) are abundant in the sediment and sometimes make up fossil benches. The deposit consists of rather homogenous carbonate bearing clays (dolomite, calcite) with high silt content; the excavation depth of clays currently used for brick production is about 25 m. Clay B originates from an allochthonous fluviatile-swamp-pond type from Pliocene-Peistocene-Holocene (fluvio-lacrusine sedimentation). Most of the deposit consists of almost horizontal clay beds with layers of grey-blue silt with pebbles at the basis (subartesian aquifer). The major part of clays is an alteration of clays of varying color, overlaying karstified limestone bedrock. The total thickness is about 11 m.
PXRD measurements on those raw clay samples used in brick production were performed on a Bruker D8 ADVANCE Eco goniometer, using CuK a as the radiation source, with a voltage of 40 kV and a current of 25 mA. The diffraction patterns were recorded for angles 2h ranging from 5 to 85 , with a step resolution of 0.01 . A second mineralogical analysis was undertaken on the \2l-fraction in order to allow for a distinction between different clay minerals. Therefore, the raw material was blended with hydrogen peroxide to remove the organic parts. After dispersing the suspension with an ultrasonic homogenizer, the separation process was conducted by a centrifuge. Various groups of clay minerals were distinguished according to a set of methods following Schultz [24]. Tempering the sample at 550 C allowed for a separation of the main peaks of Chlorite and Kaolinite in the diffraction diagram i.e. the 7Å peak of Chlorite remains, when tempered, while the 7Å peak of Kaolinite disappears.

Fired clay sample preparation
From each of the two clay types, four bricks with dimensions of 30 Â 15 Â 125 mm 3 were extruded. Two samples of each clay type were fired at 880 C (typically used in brick production), and the other two at 1100 C (upper limit in tile production, extensive sintering process). For the sample preparation a distilled water-cooled low speed saw (Isomet, Buehler, USA) was used to cut out, from the center of each of the four fired clay bricks, cuboids measuring 6 Â 5 Â 15 mm 3 with edges being parallel or orthogonal to the extrusion direction. The latter were then embedded in resin (Agar low viscosity resin, Agar Scientific, UK), followed by oven drying at 60 C for 48 h. Thereafter, each cuboid was cut into three specimens with dimensions ranging from 0:5 Â 1 Â 4 mm 3 to 2 Â 3 Â 4 mm 3 , yielding a total of 12 samples, which were then glued, by means of a cyan acrylate adhesive, onto glassy object slides, such that the extrusion direction was oriented either perpendicular to the slides plane, or parallel to one of the long edges of the slide, see Fig. 1, where x relates to the extrusion direction, and x, y, and z form a right-handed base frame. Each sample is indicated by the type of clay, firing temperature, and the direction which is normal to the image plane. For further sample preparation, the main objective was to obtain a surface of the sample as smooth as possible, which is due to various reasons: obviously, a smooth surface allows for high quality pictures, and ensures that no material phases are quarried out during polishing and get lost during analysis. Furthermore, after scanning the samples with the SEM-EDX, a huge set of measurements will be carried out (and discussed in a subsequent series of papers) on the prepared samples to identify the elastic and thermal properties of the determined material phases by means of nanoindentation [25] and scanning thermal microscopy [26], both methods requiring a very smooth surface to ensure high quality results [27,28]. Therefore, the samples underwent an extensive polishing protocol (realised with machine PM5, Logitech, Scotland). First, the specimens were ground by means of silicone carbide papers with grain sizes of 35 microns and 9 microns, respectively. These coarse grinding steps were followed by polishing the specimens with polycrystalline diamonds (DP-Spray P, Struers, Denmark). In all these polishing steps, ethylene glycol was used as a lubricant. After each polishing step, the treated surfaces were examined with a light microscope (Zeiss Axio Imager, Carl Zeiss, Germany), in order to (1) monitor the progressing polishing success in terms of a succession of line-type traces which eventually vanish, and to (2) check the presence of abrasive grains of continuously decreasing size. In order to clean the fired clay specimens from the latter abrasives, they were first rinsed with a mild detergent dissolved in tap water, then rinsed with ethanol [29], then manually polished with a napped cloth (Microcloth, Buehler, USA), and finally rinsed with ethanol. Also the polishing machine was thoroughly cleaned after each polishing step. The fixed parts of the machine were cleaned with paper cloths saturated with distilled water and ethanol, while the removable parts underwent an ultrasonic bath in distilled water. This protocol led to the realization of the smoothest possible surface, characterized by a roughness of only 5 nm, as measured by scanning probe microscopy (TI900 Triboindenter, Hysitron, USA), allowing for high quality pictures and further experiments to derive the mechanical and thermal properties of the material phases (Table 1).

Fig. 1
For each clay and each firing temperatures three specimens were prepared taking the different spatial directions into account; the x-axis indicates the extrusion direction, the z-axis is the direction of material anisotropy

Mineralogical analysis of fired samples
The four fired samples (namely two of each clay type fired at two different temperatures) were analysed by PXRD (X'Pert, Philips, The Netherlands, CuK a radiation source, angles 2h ranging from 5 to 120 , step size of 0.02 ). The diffraction pattern entered a quantitative phase analysis (QPA) provided by the software TOPAS (Bruker, USA), being based on the Rietveld method, with the amorphous phase (see Table 4) becoming accessible through the use of Korund (A16SG 1200 C=48 h, ALCOA, USA) as an internal standard [30,31]. Such is necessary to denote the observed structures from the SEM-EDX pictures to the distinct mineral phases. Additionally, a complete identification of the material phases, including the mineralogy of the glassy matrix phase, can only be accomplished with the information of the PXRD analysis, as the SEM-EDX cannot resolve structures below a certain threshold. The fine mineral particles below this threshold are, as indicated in Sect. 2.6, accounted to the matrix phase, which consists also of an amorphous glassy fraction and mineral structures decomposed during burning. The PXRD is therefore necessary as an advanced information preliminary to the SEM-EDX elemental mapping. In order to allow for directly comparing the PXRD results with those obtained from SEM-EDX (described in Sects. 2.5 to 2.6), the XRD-derived mass fractions, MF, need to be converted to volume fractions, vf, according to the following formula, for component i: with q as the mass density of the investigated clay sample, determined by Archimedes' principle [32]; and q i as the real mass density of component i.

SEM-EDX analysis
The surfaces of the polished samples were investigated by means of Scanning Electron Microscopy (Quanta 200 FEG, FEI, USA). In this context, backscattered electrons were used, and the following settings were applied: an acceleration voltage of 20 kV, a spot size of either 3 or 3.5, and the chamber pressure was set to 0:5 mbar. For each specimen, a representative area of 271:56 Â 212:16 lm 2 was chosen, and further examined by Energy-Dispersive X-ray spectroscopy (Pegasus XM4, EDAX, USA) with a dwell time of 1000 ls and 64 frames. In the course of EDX [33], the investigated samples are bombarded with a focussed beam of electrons. The latter create electron vacancies in the inner orbital shells of atoms. This provokes atom-specific X-rays (also called chemical element-specific ''characteristic X-rays''), which are emitted when electrons from the orbits further out, ''fall'' into the aforementioned vacancies. This, in turns, results in chemical elementspecific photons which are quantified in terms of their energies, i.e. 6:391 keV for Fe (iron), 1:253 keV for Mg (magnesium), and 1:740 keV for Si (silicon), the three chemical elements which are specifically targeted in the present application. The number of such photons is then visualised, pixel by pixel, in terms of an 8bit brightness value [34], yielding elemental maps associated to the representative areas of the investigated specimens. The count of photons is influenced by bremsstrahlung (continuum background X-rays) [34], spectrometer dead time effects [35], and matrix or interelement effects [36]. The latter appears in mixtures of elements due to the differences in elastic and inelastic scattering processes and in the propagation of X-rays through the sample to reach the detector. In order to account for these various effects, the automated background and matrix correction features of the commercial software Genesis (EDAX, USA) were employed. The identification of material phases, or quasi-homogeneous subdomains in terms of continuum micromechanics to which material properties can be assigned by means of a coupled SEM-EDX analysis, is described in the following.

Phase identification based on SEM-EDX data
Based on the mineralogical information from the XRD analysis as well as the aforementioned SEM images and elemental maps, consisting of 512 Â 400 pixels, seven different material phases were identified, in the following way: -Quartz grains: The Si elemental map was converted from a brightness value image to a binary image, according to Otsu's method [37]. Accordingly, the brightness value histogram was considered as a bimodal distribution, and the threshold separating pixels of one class (belonging to the background, and representing anything but quartz grains) from those of the other class (belonging to the foreground, and representing quartz) was determined such that the inter-class variance becomes a maximum. Thereafter, agglomerations of fewer than ten foreground pixels are considered as single Si spots within the matrix, and were therefore deleted by means of the particle analyser plug-in of ImageJ [38]. The resulting binary image was finally checked against the corresponding backscattered scanning electron micrograph, in order to identify grains which were ''falsely glued'' by the Otsu method. Such grains were then separated, in the binary image, from each other. This was done by means of introducing a one pixel interface of ''separation material'', consisting of the matrix phase, which is described further below.
Comparing the volume fractions from the PXRD analysis (Table 4) with results from the SEM-EDX ( Fig. 12) indicates that fine dispersed quartz particles exist in the matrix. -Feldspar grains (anorthite): The orthoclase and the anorthite identified by the PXRD cannot be distinguished from each other within the backscattered scanning electron micrograph. Plus, the feldspar grains and the quartz grains cannot be visually separated. This allows for identifying the feldspar grains as all grain-type objects which were not qualified as quartz grains according to the method described before. Fine anorthite grains mineralizing due to higher burning temperatures are, as some quartz particles, dispersed in the matrix and therefore not captured by the SEM-EDX. -Iron-magnesium mica (Fe-Mg mica): All plateshaped features in the backscattered scanning electron micrographs, which contained iron, as indicated by the corresponding elemental maps, were identified as Fe-Mg mica. Results from the PXRD and the SEM-EDX analysis show similar results with respect to the volume fractions. -Light mica (muscovite): All plate-shaped features in the backscattered scanning electron micrographs, which did not contain iron, as indicated by the corresponding elemental maps, were identified as light mica, namely muscovite [39]. -Decarbonated dolomite: This phase was identified as comprising all non-platy features containing magnesium. As the dolomite starts to decompose between $ 700 C and 800 C, the identified grains are already at the lower burning temperature of 880 C not anymore identified by the XRD, and the grains themselves are mainly remaining structures of the dissolving process of the dolomite. Therefore, the Mg elemental map was converted from a brightness value image to a binary image, according to Yen's algorithm [40]. Accordingly, a cost function representing the discrepancy between the original map and the binary map, as well as the number of bits required for the binary image, was minimized. Afterwards, Mg pixel agglomerations of fewer than ten pixels as well as with a circularity below two, were deleted by means of the particle analyser plug-in of ImageJ [38]. Within the remaining structures of the decarbonated dolomite, gehlenite, and, with higher burning temperature, diopside as identified by the PXRD may occur. -Pores: The backscattered scanning electron micrograph was converted from a brightness value image to a binary image, according to Otsu's method [37]. Accordingly, the brightness value histogram was considered as a bimodal distribution, and the threshold separating pixels of one class (belonging to the background, representing anything but pores) from those of the other class (belonging to the foreground, representing pores) was determined such that the inter-class variance becomes a maximum. However, corresponding foreground objects, appearing as dark objects on the SEM, may also contain parts of the decarbonated dolomite phase, so that all sections of the foreground which overlap with the previously determined carbonated dolomite phase, are subtracted from what is considered as the pore space. The remaining foreground pixels form the pore phase, which thereafter undergoes a visual check, in order to split ''artificially glued'' pores. As before, those pores are separated by a ''splitting material'' belonging to the matrix phase with a thickness of one pixel. Such is necessary for a correct identification of the phase morphology. The latter concludes the phase identification description. -Matrix: Structures that could not be uniquely assigned to one of the aforementioned phases were considered to be part of the fired clay matrix. This matrix phase consists of the amorphous fractions, small amounts of minerals, such as haematite and diopside, which cannot be assigned to self-contained structures observable with the SEM-EDX, as well as of fine dispersed mineral particles (like quartz and feldspar). Plus, remaining structures of clay minerals forming new crystal or melting phases, are assigned to the matrix phase. However, all this constituent phases within the glassy matrix phase cannot further be resolved, the microheterogeneous matrix is therefore considered as a quasihomogeneous subdomain in the framework of continuum mechanics.

Phase orientation distribution
In sections through the anisotropy axis of fired clay (zaxis in Fig. 1), all phases except the matrix phase are approximated as ellipses (see Fig. 2). We interpret the latter as projections of spheroids with phase-specific aspect ratios (a n =R), a n and R being the spheroids' half-thicknesses and the radii, respectively, and with the orientation distribution of their normal vectors, defined through an angle h from the anisotropy direction. In the course of such a projection, the orientation angle h of the ellipses obtained on the cross sections, being the angle between the anisotropy direction z and the projected normal vector n of the spheroids orientation, is dependent on h as well as on / (see Fig. 3

.)
If transversal isotropy is assumed, then the normal vectors n of a group of material phases have to be equally distributed along / 2 0; 2p f g. The projected orientation angles, observed on the cross sections (with a plane x ¼ 0, resp. y ¼ 0, see Fig. 3), are given to: For the case of equally distributed angles /, the mean value of the orientation distribution by means of the projected orientation angles of a group of material phases can then be calculated to: thereby giving access to the solid orientation angle h. When calculating the weighted orientation distribution by means of the standard deviation, defined via: where N is the number of inclusions of a material phase observed in a a cross section, A i =A being the fraction of the inclusion and h Ã the weighted mean value of the orientation angles, then the standard deviation of the projected orientation angles h correlates linearly with the standard deviation of the solid orientation angles h, if h Ã ¼ 0: To derive the aspect ratio of a spheroid, projected to the plane x ¼ 0, we consider a spheroid with of the form: where x, y, z are local coordinates in the direction of local, orthogonal unit vectors v 1 ; v 2 ; n (see Fig. 4). If we rotate this coordinate system to another system x; y; z in the plane x ¼ 0, so that the projection of the spheroid's normal vector n coincides with the z direction, then the normal vector of the spheroid in this new system can be written in the form: n ¼ sinðhÞ 0 cosðhÞ 2 6 4 3 7

5: ð7Þ
The vector v 2 can be arbitrarily chosen as a unit vector perpendicular to n, for example: Resulting for the third unit vector to be: Having the local coordinate system of basis vectors we can immediately write for from which we derive, using orthonormality: x y z Using this transformation rule, we can transform Eq. (6) into the rotated system. Considering the projection into plane x ¼ 0, we finally get: The obtained aspect ratio a 1 =a 3 has to be identical to those obtained from the 2D images (see Table 7), thereby giving access to the morphology of the spheroids (listed in Table 8).  Fig. 4 Cutting a spheroid (brown) with radius R and normal n deviating from the anisotropy direction z by angle h, by means of a plane through the z-axis, yields an ellipse (yellow) whose longer half-axis preserves the length of the radius, while its shorter half-axis can be computed from Eq. (11). (Color figure online) 3 Results and discussion

Raw clay samples
According to the PXRD results on the unfired samples, clay A contained carbonates (calcite and dolomite, see Table 2), while clay B is free of such, containing more quartz. Additionally, both show some shares of phyllosilicates (chlorite, muscovite and kaolinite) as well as some feldspars (albite and microcline), see Table 2.
The PXRD on the \2l-fraction reveals that clay A appears to be dominated by smectite and illite, containing also chlorite and kaolinite. In contrast, Clay B (poor in Ca), is dominated by vermiculite, with some constituents of illite, chlorite and kaolinite (see Table 3).

Fired clay samples
Results from the PXRD analysis for the fired clay are given in Table 4.
Clay A, representing a clay rich in carbonates, shows, in agreement with literature [42], a wide compositional variability with respect to the carbonate-poor clay B (Table 4, Fig. 12). This can be explained by the breakdown of chlorite [43] at $ 500 C, as well as calcite and dolomite at $ 850 C, releasing Ca and Mg, which produce a compositional array of new forming minerals with the surrounding clay minerals [44]. With higher burning temperature, clay A shows a significant decrease in the quartz content, something which can be typically expected for carbonate-containing clays [45,46]. This occurs due to a sudden onset of a melting phase at temperatures of $ 1080 C, which can be observed specifically at the boundaries of the quartz grains (see Fig. 5a). Riccardi et al. [42] reported that such transformation at a wider reaction layer around quartz, beginning at temperatures of $ 1050 C, can be observed in Ca-rich clays only, and can be explained with the nucleation of newly forming anorthite round the quartz grains due to the presence of diffused calcite in the clay matrix. Besides, the high-temperature modification cristobalite is forming at a burning temperature of 1100 C. The highest amount of newly forming crystals within clay A is represented by anorthite, which has also been previously reported by other authors [45], followed by small amounts of diopside and the aforementioned cristobalite. As the raw material contains dolomite and calcite, which decarbonate at temperatures of $ 700 C to 800 C, these minerals are not detected by the PXRD within our samples. Also, the PXRD results for the fired clay at 880 C do not clearly indicate the presence of Mg and Ca containing minerals, oxides or hydroxides, a problem which has already been adressed by Cultrone et al. [47]. However, the original microtextural site can still be observed at a burning temperature of 880 C in the SEM, significantly enriched in Mg and Ca (see Fig. 12a). With higher burning temperatures, and the former dolomite and calcite grains fully reacted, only ring-like textures remain still visible (see Fig. 12), a  phenomenon which has also been reported by Riccardi et al. [42]. The calcite forms, in addition with siliciumdioxide of clay minerals, anorthite and gehlenite, while gehlenite as an intermediate forming mineral converts with SiO 2 from clay minerals or quartz grains to anorthite with its higher silicid acid content during the sinter reaction [45,48]. Muscovite, paragonite and kaolinite, along with the microcline, are partly dissipated during the firing process, but still observable as muscovite and orthoclase at a burning temperature of 880 C. The muscovite and orthoclase then dissolve, as expected from literature [47,49], at high burning temperatures. Regarding the constant amount of hematite, the following can be stated: all calciumaluminum silicates accept varying amounts of iron within their mineral structures [45,50], which is why only little amounts of pure hematite are found, and no new crystallization can be observed at higher temperatures. Finally, the absence of albite at burning temperatures exceeding 850 C is caused by its change of the chemical composition by reaction with the surrounding clay matrix [44], when a rim of  5 SEM micrographs of (a) quartz grains melting due to higher burning temperatures at 1100 C within clay A and (b) newly formed mullite crystals in the ''glassy'' vitrified matrix phase of clay B due to higher burning temperatures at 1100 C K-feldspars is constantly growing, reducing the Na content in the former albite crystal. Clay B, which is poor in carbonates, shows a substantial different behaviour. As iron is not trapped within the lattice of calcium-aluminum silicates in clay B, an increase in hematite content can be found with higher burning temperature [43]. As with clay A, quartz and muscovite dissolve (partly) at burning temperatures of 1100 C, while the only newly formed mineral is mullite at $ 1050 C [46], observed also in SEM pictures (see Fig. 5b). Cultrone et al. [49] indicate that large formations of mullite out of the melt phase can be observed. Such would explain the decrease of Al 2 O 3 within the matrix at a burning temperature of 1100 C, see Table 6. Compared to clay A, a more distinctive melt formation can already be observed at the lower burning temperature of 880 C. This is caused by the decomposition of carbonates in clay A, promoting the fast crystallization of calcium-aluminum silicates, at the expense of the formation of a melt phase [43].
As expected, the volume fraction of SiO 2 within the amorphous part of the matrix is increasing with a higher burning temperature due to the formation of a melting phase (see Fig. 5a) and the decomposition of clay minerals. While the content of Al 2 O 3 in the matrix of clay A stays constant, the newly forming of mullite crystals in the matrix (see Fig. 5b) reduces the fraction of Al 2 O 3 in clay B. The small reduction of MgO 2 in clay A can be explained by the newly forming diopside, while the reduction of Fe 2 O 3 in clay B can be assigned to the increasing mineral content of hematite.
To derive the chemical composition of the matrix, containing the amorphous fraction, mineral structures decomposed during burning as well as fine dispersed mineral grains, the PXRD data can be used to estimate the chemical composition. Besides the physical properties, the RRUFF TM Project database [41] contains the chemical composition of minerals, which we can use in combination with the results from the PXRD to compute the mass fractions of molecules recorded by the XRD. Comparison of the results with the XRF data (see Table 5) reveals the chemical composition of the amorphous fraction, given in Table 6.

SEM-EDX image analysis
Otsu's method applied to the Si elemental maps, combined with the requirement of 10 pixel minimum grain size, and after ''splitting artificially glued grains'' according to an SEM-based visual check, allows for satisfactory identification of the quartz grains in all investigated samples, compare images labelled with (a) and (c), respectively, in Fig. 6 (see also supplementary material). Grain-type objects, not identified as quartz grains by Otsu's method applied on the Si elemental maps, thereafter were assigned as feldspar grains (see Fig. 7), following the protocol in Sect. 2.6. Overlaying plate-shaped features in the backscattered scanning electron micrographs (see Fig. 8a) with the corresponding Fe elemental maps  Fig. 8b) allows for satisfactory identification of Fe-Mg mica and of muscovite, respectively, see cyan and blue feature in Fig. 8c. Yen's algorithm applied to the Mg elemental maps, combined with the requirement of 10 pixels minimum grain size and of 0.2 minimum circularity, allowed for the identification of the decarbonated dolomite phase in all investigated samples, compare images labelled with (a) and (c), respectively, in Fig. 9 (see also supplementary material). Otsu's method applied to the backscattered scanning electron micrographs, and subsequent subtraction of all foreground features which are already assigned to the decarbonated dolomite phase, followed by splitting ''artificially glued pores'' according to an SEM-based visual check, allows for satisfactory identification of the pore phase, see Fig. 10 (see also  supplementary material). As the aforementioned binary maps are representative for the entire (spatial) clay samples, their mean grey values give direct access to the phase volume fraction, see Fig. 11. At 880 C firing temperature, the quartz grain volume fraction obtained from SEM-EDX are consistently lower than the quartz volume fraction from XRD, compare  Table 4. This indicates that at this lower firing temperature, quartz appears not only in grain form, but also as fine dispersed particles within the matrix phase. At 1100 C firing temperature, the quartz grain volume fractions obtained from SEM-EDX agree very well with the XRD-derived quartz volume fraction, compare Fig. 11c, d and Table 4. This indicates that at this higher firing temperature, quartz appears virtually exclusively in the form of several micrometer-sized grains. At 880 C firing temperature, the volume fractions of muscovite obtained from SEM-EDX agree fairly well with those derived from XRD, compare Fig. 11a, b and Table 4. At 1100 C firing temperature, the SEM-EDX analysis confirms the disappearance of muscovite, as already known from the PXRD results. The higher porosity of clay A (compare Fig. 11a, c to 11b, d) indicates the higher carbonate content in its raw state [49], cf. 2.3. Carbonate-bearing clays are known to be rather heterogeneous [45], which is consistent with clay A exhibiting indeed more material phases than clay B (compare Fig. 11a, c to 11b, d). 1100 C ignition-induced disappearance of pores smaller than 2l m (see Fig. 14) is consistent with the vitrification processes described by Cultrone et al. [49] and Zouaoui et al. [51]. However, micropores of size below the resolution of the SEM-images with 0:192 lm pixel size may still prevail. This is suggested by an intravoxel microCT-image analysis which was experimentally validated by weighing tests and mercury intrusion porosimetry [52]. Correspondingly obtained porosity values are about 8% higher than those obtained from present SEM-EDX results of Fig. 11a, c. These additional pores of clearly submicrometer size are part of the matrix phase. Simultaneous vanishing of the muscovite phase is probably due to dehydroxylation and weakening the long range structural organization of the latter, as reported by Gridi-Bennadji et al. [53] and Barlow and Manning [54].

Morphological analysis
From a morphological viewpoint, the overall emerging picture is the following: at 880 C firing temperature, a contiguous matrix contains slit-type pores, quartz and feldspar grains, as well as platy muscovite inclusions; in clay A, also grains of decarbonated dolomite appear, and with a burning temperature of 1100 C degenerate to pore coatings, see Figs. 12a, b and 13a, b. At 1100 C firing temperature, the pores become more ellipsoidal in shape. Regarding these SEM images, the following can be stated: the induced secondary porosity in clay A due to the decarbonation of dolomite and calcite (see [42,49]), can be observed, while porosity in clay B at lower burning temperatures seems to be mainly driven by shrinkage of the clay matrix (Fig. 13). With higher burning temperatures, the formation of a melt phase leads to more ellipsoidal and spherical pores, compared to the sharp edged pores at a burning temperature of 880 C, which is in accordance with observations by Cultrone et al. [47].
Quantitative morphological information for all distinguished material phases is given in Tables 7  and 8, and reveals the preferential orientation to be with the extrusion axis, something which has already been observed by previous authors [7,13]. In general, it can be stated that the more elongated the inclusions are, the more their preferential orientation is with the extrusion axis. While the less elongated inclusions (quartz, feldspar and decarbonized dolomite) show a huge weighted standard deviation of the preferential orientation, the elongated inclusions (phyllosilica) are aligned almost parallel to the extrusion axis.

Conclusions
Conclusively, our new phase identification protocol, stemming from the combination of SEM-EDX images with a preliminary conducted PXRD, is fully consistent with many previous investigations, while marking, at the same time, a new level of completeness, as the entire space of SEM-resolved microstructure is clearly assigned to one of up to seven phases with different shapes. Additionally, the chemical composition of the matrix phase has been calculated using obtained data by PXRD and XRF measurements, confirming the assumption of finely dispersed quartz and feldspar particles below the resolution threshold of the SEM. Plus, these results are in line with the mineralogical analysis of the unfired raw material and the fired clay at different burning temperatures. Finally, the morphological evaluation of all material phases allows for a thorough description of the morphometrics of the latter. This opens the way towards a more quantitative understanding of structure-property relations in fired clay, in particular for mechanical and thermal properties predicted with tools of continuum micromechanics or homogenization theory [15]. The latter proved as very suitable for this purpose, in the context of various construction materials. With the presented identification of the material phases and their morphology, the first step towards such a physically based multiscale material model able to estimate macroscopic elastic and thermal properties of fired clay [26] is taken. 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/.