Contrast in Terahertz Images of Archival Documents — Part I: Influence of the Optical Parameters from the Ink and Support

This study aims to objectively inform curators when terahertz time-domain (TD) imaging set in reflection mode is likely to give well-contrasted images of inscriptions in a complex archival document and is a useful non-invasive alternative to current digitisation processes. To this end, the dispersive refractive indices and absorption coefficients from various archival materials are assessed and their influence on contrast in terahertz images from historical documents is explored. Sepia ink and inks produced with bistre or verdigris mixed with a solution of Arabic gum or rabbit skin glue are unlikely to lead to well-contrasted images. However, dispersions of bone black, ivory black, iron gall ink, malachite, lapis lazuli, minium and vermilion are likely to lead to well-contrasted images. Inscriptions written with lamp black, carbon black and graphite give the best imaging results. The characteristic spectral signatures from iron gall ink, minium and vermilion pellets between 5 and 100 cm − 1 relate to a ringing effect at late collection times in TD waveforms transmitted through these pellets. The same ringing effect can be probed in waveforms reflected from iron gall, minium and vermilion ink deposits at the surface of a document. Since TD waveforms collected for each scanning pixel can be Fourier-J transformed into spectral information, terahertz TD imaging in reflection mode can serve as a hyperspectral imaging tool. However, chemical recognition and mapping of the ink is currently limited by the fact that the morphology of the document influences more the terahertz spectral response of the document than the resonant behaviour of the ink.

its path. As a result, the textual information from a specific layer would be difficult to recover on the image obtained from the amplitude of the transmitted pulse, as shadows of inscriptions from other layers within the document would be superimposed to the information of interest [9]. In reflection mode, however, such shadows or occlusions on the image obtained from the amplitude of the pulse reflected on the n th layer can be removed by using the image obtained from the amplitude of the pulse reflected on the (n − 1) th layer. The reflection configuration therefore seems more appropriate to image archival documents.
Finally, some pigments have been shown to possess a dispersive behaviour: their terahertz spectral signatures can be recovered from the Fourier-transform of TD waveforms [16,17]. Terahertz TD imaging can therefore also serve as a hyperspectral imaging tool. Combined with a highly sensitive detection system, it has been used in transmission mode as a chemical imaging tool to discriminate ink inscriptions from stains on the document [16]. Whether the reflection mode can be used for this purpose is yet to be determined. The use of terahertz TD imaging for chemical imaging is relevant to paper and book conservators as it could help distinguishing between inert and corrosive inks, such as carbon-based and iron gall inks [18], in a non-invasive way, and treat or store the document accordingly.
The aims of this study are to assess the optical properties in the terahertz region 0.15-3 THz from materials commonly found in archives and to use this information to better understand and interpret contrast in terahertz images from ink inscriptions scanned in reflection mode. This study further exploits these terahertz optical properties to refine the use of terahertz TD technology for chemical mapping.

Pigments and Binders
The selection of pigments and inks used in this study was based on the frequency of their use in archival documents [19][20][21], their predisposition to degrade the support, the availability of their components and the feasibility of their production in a laboratory. The nature of pigments acquired from art suppliers was confirmed with Fourier-transform infrared (FTIR) spectroscopy, acquired either in attenuated total reflectance (ATR) or transmission mode. Measurements were performed using the Alpha-P Bruker FTIR spectrometer with OPUS software, version 6.5. Spectroscopic data was acquired within the wavenumber range 375-4000 cm −1 , where characteristic absorption peaks from both organic and inorganic compounds can be found. Sixty-four scans were performed at 4 cm −1 spectral resolution. These spectra are supplied in Online Resource 1.
It should be noted that ivory black was acquired as Bgenuine ivory black^, while bone black was acquired from the same supplier as Bivory black^. A discussion with the supplier revealed that bone black was more likely to have been produced from calcined bones than ivory. This was confirmed by the presence of vibrations attributed to cyano stretches in cyanamide, characteristic of burnt bone, at 697 and 2012 cm −1 [22,23] in the ATR-FTIR spectrum of the Bivory black^supply, while these bands were not present in the spectrum of the Bgenuine ivory black^supply. Furthermore, neither the ATR nor transmission FTIR spectra from vermilion showed spectral feature, which is consistent with the literature [23,24]. However, Derrick et al. remind us that the red pigments cadmium red and realgar also show no spectral band in the mid-infrared region [23]. Still, the terahertz spectral signatures of these three red pigments are different [17]: vermilion shows sharp spectral bands around 1.1, 1.3 and 2.6 THz, which is in agreement with the terahertz spectrum of the vermilion supply presented in Fig. 7. It can therefore be concluded that the Bvermilion^supply is constituted of pure vermilion. The Bminium^supply showed intense spectral features at 387, 453 and 532 cm −1 , as previously reported for minium and lead(II, IV) oxide (Pb 3 O 4 ) [25,26]. It also contains minute amounts of lead carbonates, as attested by the weak bands at 685 and 1410 cm −1 [27,28]. The broad and weak absorption bands at 489 and 511 cm −1 may attest the presence of litharge and massicot [29], respectively, two polymorphs of lead(II) oxide (PbO) [20]. The supplier of lapis lazuli was not able to attest the provenance of this pigment (Chilean or Afghan). Yet, studies report that Afghan lapis lazuli can be differentiated from Chilean lapis lazuli and synthetic ultramarine thanks to the presence of an absorption band at approximately 2340 cm −1 [23,30,31]. This vibrational band was observed in none of the three ATR-FTIR spectra from the lapis lazuli powder investigated here: it can therefore be deduced that our lapis lazuli supply is not of Afghan origin. Its ATR-FTIR spectrum showed intense overlapping bands at 960-1000 cm −1 also found in the lazurite sample from the Forbes Pigment Database [32], but shows no band associated to calcium carbonate (RRUFF ID: R040070) or pyrite (RRUFF ID: R050190) from the RRUFF online database [33]. It can therefore be assumed that the supply is largely composed of pure lazurite. Finally, FTIR spectra from the supplies of bistre, malachite, verdigris, rabbit skin glue and Arabic gum confirmed that they are respectively constituted of chars from burnt lignin (a polymeric aromatic organic substance in wood) [34], copper carbonate (Cu 2 CO 3 (OH) 2 ) (RRUFF ID: R050508) [33,35], copper(II) acetate monohydrate (Cu(CH 3 COO) 2 ⋅H 2 O) [27,36], proteins [23] and polysaccharides also found in Arabic gum of industrial grade [23,37].

Pellet Preparation
For terahertz spectroscopic analysis, dry powders were analysed, and the solutions of sepia ink and iron gall inks were therefore left to dry on a ceramic plate, scraped off and ground into a powder with a pestle and mortar before analysis. Pure powders of binders, inks or pigments were pressed into three 400-mg pellets with a 10% mass concentration of sample to high density polyethylene (HDPE) powder and a diameter of 13 mm, following the optimal procedure described elsewhere [38].

Supports
Historically, cotton and linen rags have been widely used as a source of cellulosic fibrous materials for paper making. However, paper sheets containing only cellulosic fibres are water absorbent: aqueous inks would diffuse into the porous fibre mesh and the resulting ink lines would be less defined. Most paper sheets are therefore sized, i.e. impregnated with substances to prevent ink from spreading into paper fibres [39,40]. Paper sheets can be heavily sized, lightly sized or not sized at all [40]. Parchment is prepared from wet unhaired animal skin soaked in a lime solution (Ca(OH) 2 ). After the removal of hair from the limed skin, the drying stretched skin is rubbed with a pumice stone to smoothen its surface. To improve the smoothness and suitability of the parchment skin for writing purposes, its surface can further be coated with egg white, oil, glue or a solution made of soaked parchment scraps, or can be treated with hot water to form a superficial layer of size [41][42][43]. This process further ensured that the writing ink would not penetrate too deeply in the parchment fibre mesh. To improve the adhesion of ink on parchment and prevent it from running, the surface may be dusted with powdered chalk (CaCO 3 ), calcium sulphate dihydrate (CaSO 4 , 2H 2 O) and quicklime (calcium oxide, CaO) [44]. Six rag paper sheets from the nineteenth century and one parchment sheet from 1752, from the Institute for Sustainable Heritage Historic Reference Material Collection, were used as supports. To observe the influence of paper coating on terahertz images, three paper sheets were additionally sized by immersing them in a gelatine solution for less than 10 s, then placing them between two sheets of polyester and pressing them between wooden boards. The preparation of the size followed the instructions from Kolar and Strlic [45], mixing 12,500 mg of gelatine from bovine skin (Product number: G9382, Sigma Aldrich) and 1425 mg of alum (aluminium sulphate hydrate, Product number: 368458, Sigma Aldrich) in 0.5 L of deionised water. Areas of 2 × 2 cm 2 were cut from each sheet of rag paper and from three non-adjacent areas of the parchment sheet for terahertz spectroscopic analysis.

Terahertz Time-Domain Spectroscopy
Pellets and supports were analysed in transmission mode using the TPS Spectra 3000 (TeraView Ltd., Cambridge, UK). Samples were placed in a nitrogen-purged sample compartment, at the focal plane of a focused terahertz beam; 1800 scans were acquired at a scan rate of 30 Hz and co-averaged to minimise random noise. The working spectral range of the frequency-domain (FD) signals obtained through the Fourier-transform of the waveforms was 5-100 cm −1 , with a spectral resolution of 1.2 cm −1 . The signal transmitted through a 360-mg pellet of pure HDPE powder was used as a reference measurement for the analysis of pellets, while the signal transmitted through the empty sample compartment was used as a reference for the analysis of supports. Constant and frequency-dependent refractive indices, together with frequency-dependent absorption coefficients, were calculated as described elsewhere [38]. Additional pulses resulting from multiple internal reflections on the sheet surfaces were detected in the waveform transmitted through each sheet of support. Such phenomenon is often called a Fabry-Perot or Betalon^effect [46]. These subsequent pulses overlapped with the main pulse transmitted through the sheet and therefore changed the time delay between the main peak of the raw waveform transmitted through the sample sheet and the main peak of the reference waveform. To avoid such errors in readings of the time delay, also inducing errors in the calculation of the constant refractive index, deconvolution [11,47] was used to separate the main transmitted pulse from the pulses originating from the etalon effect. A bespoke deconvolution algorithm with a double Gaussian filter function was implemented in Matlab [48] for this purpose.

Ink Preparation and Application
For terahertz imaging experiments, ink solutions were prepared by dispersing the pigment, with no further grinding, in a solution of Arabic gum. Although a broad range of ratios of binder over water can be found in the literature, the same proportion recommended by Kolar and Strlic [45] was consistently used (i.e. 625 mg of Arabic gum with 25 mL of deionised water heated at 70°C and left to cool down before the addition of the pigment), as it gives suitable viscosity and increases control over the experiments. Inks were prepared with a mass ratio of pigment to gum of 1:2. Mass ratios were preferred over volume ratios, as the volume of a pigment would be difficult to measure accurately due to its inter-and intra-particle pores, and as the penetration properties of the binder into the pigment pores are not known. To demonstrate the influence of ink on image contrast in a document, any indentations left by the writing tool should be avoided, as they can influence the image contrast [49]. Application with a brush or using screen printing was also avoided as the volume of ink and pigment concentration applied could not be uniformly distributed, which would have increased the complexity of the experiment and interpretation. Application with an ink jet was not considered, due to the colloidal nature of the pigment suspensions, which would have risked clogging the ink jet heads. For these reasons, 2 mL of each fresh ink was applied on historical rag paper and parchment sheets using pipettes. The resulting ink deposits had a diameter of 2 ± 0.5 mm, depending on the affinity of the ink to the support. Although the diameter of ink inscriptions written with most contemporary pens is usually smaller than 2 mm, some inscriptions written with quills, nib pens or brushes can be wider than 2 mm: the diameter of these ink droplets is therefore not unrealistic with regards to the width of archival ink inscriptions. Additionally, the ink spot was significantly larger than the step size and associated spatial resolution of the imaging set-up (see section below), which ensured that observation of the dot in the terahertz image was possible.
Despite the care taken in controlling the amount of pigment and ink applied on the support, it should be noted that the colloidal nature of these ink suspensions often led to a non-uniform deposition of particles during evaporation of water from the ink drop, creating ring-shaped solid residues (also known as Bcoffee ring effect^ [50]). Additionally, application of ink with a pipette sometimes introduced microscopic bubbles in the droplet, creating local depletion of particle concentration. Finally, ink may penetrate non-uniformly into fibres constituting the support. All these phenomena led to the non-uniform distribution of ink thickness and pigment concentration within some ink dots: these variations of pigment concentration and topography have been recorded by acquiring micro-photographs of the ink dots (displayed alongside terahertz images of the same ink dots in Figs. 2, 3, 4 and 5) and accounted for in the interpretation of terahertz images.

Set-up
Terahertz images were acquired in reflection mode by fitting a reflection imaging module in the nitrogen-purged sample compartment of TPS Spectra 3000. This module comprises of a software-controlled motorised scanning table with a set of ellipsoidal mirrors underneath, focusing the generated terahertz pulses onto the holed scanning stage at an angle of incidence of 30°. Samples were placed on the scanning table and maintained steady and flat by placing them between two ring-shaped magnets taped to the stage.

Image Resolution
We define spatial resolution as the minimum distance between two separate features on a captured object leading to discernible features in the image of this object. Theoretical criteria, such as the Rayleigh or the Buxton criteria, are commonly applied to estimate the spatial resolution of an optical system [51,52]. These criteria use the shape of the 2D lateral distribution of the irradiance or amplitude of the electromagnetic field emitted by a point source, also called the point spread function (PSF). For example, the Rayleigh criterion assumes that the PSF of the irradiance is a first-order Bessel function in order to evaluate the limit of resolution r = (1.22 * f * λ) / D, radius of the first dark Airy ring of the PSF. The optics of the terahertz imaging system used in this study are such that the ratio of the focal length f and aperture diameter D of the imaging system is equal to 1. It would follow that the limit of resolution is approximately 180 μm at 2 THz. However, to confirm the distribution governing the PSF of the amplitude of the terahertz pulsed electromagnetic field and better select the relevant criterion for spatial resolution, the emission profile from the terahertz source would need to be recorded at the focal plane of the imaging system. While placing a CCD camera at the focal plane would only enable to obtain an image of the irradiance profile (and not the amplitude of the pulse at a specific time or frequency) [53], an alternative is to place a terahertz detector on a translation stage at the focal plane to record the terahertz pulse point by point along a line [54]. The subsequent Fourier-transform of each of these pulses recorded at different positions would enable to recover the discrete PSF of the terahertz electromagnetic field amplitude at specific frequencies. However, since a commercial terahertz set-up with fixed optics and antennas was used in our study, such experiment could not be performed. Another alternative consists in recording the minimum distance between two separate features on a 1951 USAF resolution test target (Product ID: R3L3S1P and R3L3S1N, Thorlabs Inc) leading to discernible features in the image of this test target. The amplitude at 2 THz of the Fourier-transformed reflected time domain waveform recorded at each location in a regularly spaced 2D grid was used to produce a 2D image of the target. At the limit of resolution defined by Rayleigh, two features on the grid are not resolved if their composite intensity distribution, sum of the PSF from the two features, exhibits a central depression corresponding to a decrease in intensity of 19% [51]. For the image of the grid at 2 THz acquired with a lateral step size of 100 μm, such central depression is observed for features spaced by 180 μm. It is interesting to note that this figure coincides with the diffraction limit calculated above using the Rayleigh criterion. The following assumptions may therefore hold true in our case: (i) the PSF from the pulsed terahertz beam follows a first-order Bessel distribution, and (ii) the resolution of the system is only limited by the diffraction phenomenon when using a small scanning step size. Using the maximum amplitude of the reflected time domain waveform to produce an image of the grid, the limit of resolution as defined by a central depression in amplitude of 19% was observed for features spaced by 280 μm for a lateral step size of 100 μm. It should finally be reminded that this limit of 19% was considered satisfactory by Lord Rayleigh because it matched what the eye could detect through telescopes or microscopes [55]. Considering the resolution and size of today's computer monitors, this limit might be obsolete and might be better defined as the spacing at which a central depression similar to the standard deviation of the amplitudes collected from a blank area of the grid is reached. On the image of the grid produced from the maximum amplitude of the pulse, this new limit is reached for features spaced by 180-200 μm, for a lateral step size of 100 μm,

Signal to Noise Ratio
The signal to noise ratio is described here as the maximum value of the dynamic range. The dynamic range was calculated using the logarithmic ratio of the Fourier-transform of the TD signal reflected from a gold-coated mirror and the frequency-domain noise floor [56,57]. The signal to noise ratio was measured to be ≈ 54 dB at 0.6 THz (r Rayleigh = 609 μm). At 2 THz, the dynamic range value lowers to 36 dB, until reaching its minimal value of 0 dB at 3.2 THz (r Rayleigh = 115 μm).
The maximum amplitude of the pulse in the time domain represents the average reflected power across the dynamic range [58]. It should therefore be noted that terahertz set-ups with a narrower dynamic range (i.e. poorer dynamic range values in the 1-3-THz region) than the system used in this study would lead to a poorer resolution in images produced from the maximum pulse amplitude than the values reported in the above section.

Signal Processing
Terahertz images from features at a specific interface within a stratified object can be obtained by identifying the pulse corresponding to the layer of interest in the reflected waveform, and producing a map of its amplitude throughout the scanned area [59][60][61][62]. For a layer of ink applied on a support 80 to 300 μm thick, two to three interfaces may be considered, depending on the thickness of the ink layer and penetration of the ink into the fibres of the support: the air-ink interface, the inkfront surface of the support interface and the back surface of the support-air interface. Depending on the thickness of the ink and support layers with respect to the bandwidth of the incident terahertz pulse, pulses reflected from each of these interfaces may overlap. This overlap would be greatly influenced by the optical properties of the ink and the support and would change the amplitude, time position and shape of the peaks observed in the collected waveform. Using the amplitude of a specific peak in the collected waveform without further signal processing to objectively compare the contrast in terahertz images from different ink dots on different supports would therefore be unsound. The comparison should rather be based on imaging data obtained from waveforms processed to minimise the influence of the overlap between pulses reflected from different interfaces.
Deconvolution attempts to retrieve the theoretical impulse response function and thereby improves the separation between the two (or three) pulses [11,47]. Deconvolution of the collected waveform therefore enables the retrieval of a signal in which the peak amplitudes are little or no longer affected by the pulses overlap. Terahertz images from ink dots presented in this study were therefore produced using the maximum amplitude of the deconvolved signals, rather than the raw signals, calculated for each scanning position (i.e. each pixel) using the built-in deconvolution tool in the TVL Imaging Suite software (TeraView Ltd). This deconvolution tool uses a double Gaussian filter [47]. The same high-frequency and lowfrequency cut-offs were consistently used for all scanned ink dots.

Normalisation of Images and Calculation of Contrast
To objectively compare the image contrast between ink dots, the maximum amplitude value M(i,j) of the deconvolved waveform obtained for each pixel (corresponding to the amplitude of the pulse reflected at the air-front surface of the sheet interface) was normalised using the following formula: with M reduced as the cropped image to obtain an image size of 3.6 × 3.6 mm 2 (area within the user-defined red rectangle in Fig. 1) and sSUB avg the average of the pixel values corresponding to the support (area outside of the user-defined yellow circle in Fig. 1). All grey scales were set from 1 (black) to 1.5 (white): while the support always appears dark in the normalised terahertz images, the more reflective the ink dot is in the terahertz region, the brighter it appears in the image. Fig. 1 Terahertz image of a vermilion ink dot obtained, after deconvolution, from the maximum amplitude of the reflected waveform. The image is overlaid by a schematic diagram describing the boundaries used for the image normalisation and calculation of the contrast value C. Red lines correspond to the cropped area to obtain a 3.6 × 3.6-mm 2 image (corresponding here to 48 pixels, with one pixel corresponding to a 75-μm scanning step). All pixels within the blue circle, the diameter of which corresponds to the minimum of the distances between blue horizontal or vertical lines, correspond to the region of interest within the ink dot. All pixels outside the yellow circle, the diameter of which corresponds to the maximum of the distances between yellow horizontal or vertical lines, correspond to the region of interest of the support To enable the simultaneous comparison of contrast between terahertz images of 43 different ink dots on supports, a contrast value C was calculated from all terahertz images and summarised in Table 2: with sROI avg as the average of the pixel values corresponding to the inked area (area inside the user-defined blue circle in Fig. 1). Regions of interest excluded the vicinity of the ink dot edges, as upon drying, some ink droplets created ring-shaped solid residues inducing unpredictable scattering of the terahertz beam. Large differences of reflectivity between the ink and the support lead to large C values.

Step Size
When using a scanning step size (i.e. a pixel size) half the width of the ink dot, there is always a pixel which is fully inscribed within the dot and faithfully renders the signal from the ink, showing maximum contrast, regardless of the position of the dot with regards to the pixel grid. Using a scanning step size below half the width of the ink line does not, in theory, improve the contrast of the dot. Yet a smaller step size may enable to render the width and shape of the ink dot more accurately and to appreciate small variations in peak amplitude due to changes in topography and/or concentrations of the ink within the dot. To illustrate this concept, terahertz images of the same vermilion and minium ink dots on rag paper were acquired with two different step sizes: 75 and 250 μm (images provided in Online Resource 2). The difference in calculated contrast values C between images with a different step size was 0.05 for both the vermilion and the minium ink dot (C vermilion, 75 μm = 1.48, C vermilion, 250 μm = 1.43, C minium, 75 μm = 1.37, C minium, 250 μm = 1.32). In comparison, the difference in image contrast recorded for two dots of vermilion ink applied on the same support and scanned with the same scanning step size can be as large as 0.22. The influence of the step size on the image contrast is therefore not considered to be significant. Still, the rendering of the width, shape and ink distribution of the ink dot is poor in images acquired with a 250-μm step size. To recover it as accurately as possible, images can be processed by increasing their sampling rate through the use of an interpolation filter. The Matlab function imresize has been used for this purpose [48].
The resulting images are shown in Online Resource 2 and have the same number of pixels as the images acquired with a 75-μm step size. The width, shape and ink distribution of the ink dot in these processed images closely resemble those observed in images acquired with a 75-μm step size, although a blur can still be observed. Additionally, the change in contrast value compared to unprocessed images does not exceed 0.02 (C vermilion, 250 μm, image resizing = 1.43 and C minium, 250 μm, image resizing = 1.30): image resizing has therefore no significant influence on the image contrast value. Reducing the step size to 75 μm therefore seems to add little to the information already enclosed in images acquired with a 250-μm step size. From the discussion in the above sections on image resolution and signal to noise, we learned that the diameter of the 2D PSF of the terahertz pulsed beam can be estimated to be above 230 μm (i.e. r Rayleigh, 3.2 THz * 2). It follows that with a step size of 250 μm, the amount of information on the interaction between the incident pulse and the sample that can be collected is maximised. With a step size of 75 μm, this information is therefore oversampled. This is why image resizing yields similar results as reducing the step size but with a slightly poorer rendering of the width, shape and ink distribution of the dot. In this study, scans were acquired with a step size between 75 and 150 μm, with no further image resizing. However, for future terahertz studies of documents, to avoid acquiring large size data and to reduce data acquisition time while rendering the width, shape and ink distribution of the ink inscriptions as accurately as possible, a large scanning step size (provided this scanning step size is less than half the width of the ink line) combined with image resizing can therefore be used.

Categorisation Based on the Constant Refractive Index
Constant refractive indices from all pellets and supports are displayed in Table 1. Interestingly, there seems to be a correlation between the refractive index values from pigments and binders with their electrical conductivity properties. According to Coutaz et al. [63], dry dielectrics have a low refractive index (below ≈ 2), intrinsic semiconductors have a larger refractive index (between ≈ 2 and 5), while metals have a very large refractive index. Coutaz et al. further remind that metals in the terahertz range can be considered as perfect electrical conductors, i.e. as having electrons unbound to atoms and free to move under the influence of the applied terahertz electromagnetic field. The Drude model can therefore be applied to describe the electronic and optical properties of metals in the terahertz range. This model helps to evidence that metals absorb much more strongly terahertz radiation than they reflect it. Paper and parchment are dielectrics. Arabic gum constituted of carbohydrates (as in paper) and rabbit skin glue constituted of collagen (as in parchment) can also be considered as dry dielectrics. It is therefore not surprising that these materials exhibit the lowest refractive indices in Table 1.
Graphite is a semi-metal (i.e. a conductor along the plane of each carbon sheet) [64]. It is thus expected to have a high refractive index and an even higher absorption coefficient: this would explain why no signal is transmitted through a pellet of graphite with a 10% mass concentration. Carbon black comes from the incomplete combustion of gaseous hydrocarbons and, as pointed out by Hauptman et al. [65], consists of graphite-like crystalline and amorphous carbon domains. The authors also indicate that the major factor that determines the electrical conductivity of carbon black particles is the graphitic character giving delocalized π-electrons [65]. Lamp black, ivory black and bone black also result from the combustion of carbonaceous materials [20] and may also contain some graphitic and some amorphous carbon. Yet, ivory black and bone black also contain inorganic calcium phosphate and calcium carbonate: they would therefore have a lower mass fraction of carbon than graphite, carbon black or lamp black within the 10% mass concentration pellet. As a result, they are likely to contain a lower mass fraction of graphitic carbon too, and would be less conductive, so would have a lower refractive index and absorption coefficient. This would explain why a transmitted signal is measured through pellets of ivory black and bone black. Bistre contains aromatic carbons connected to functional groups containing hydrogen and oxygen atoms, resulting from the charred lignin [34]: there is therefore an even smaller mass fraction of pure carbon to potentially form graphitic structures in bistre pellets than in ivory black or bone black pellets. It could therefore be assumed that bistre is less conductive than ivory black and bone black, which would be why it shows a lower refractive index in Table 1. Minium (Pb 3 O 4 ) and vermilion (or cinnabar, α-HgS) can be considered as semiconductors [66,67]. Minium and vermilion therefore have a larger refractive index (close to 2 as reported in Table 1) than dielectrics such as Arabic gum or rabbit skin glue. The image contrast of ink inscriptions on support depends on the difference between the refractive index of the ink and the refractive index of the support. Inks with a refractive index similar to the refractive index of the surrounding support are unlikely to give high image contrast and are presented as roman in Table 1. Inks with a refractive index much larger than the refractive index of the surrounding support are very likely to give a good image contrast and are presented as bold in Table 1. All other inks with refractive indices larger than the refractive indices of supports, but still allowing the signal transmitted through the pellet to be detected, are likely to give a well-contrasted terahertz image. This categorization of pigments and inks based solely on the semi-quantitative comparison of their refractive indices with those from supports may appear subjective and needs to be confirmed by observing the contrast in terahertz images from ink inscriptions made with those same pigments and binders and applied on those same supports.

Images from Peak Amplitude
The first thing to observe when reading Table 2 and comparing Figs. 2, 3 and 4 is that lamp black ink dots systematically lead to a very good image contrast. Additionally, for a given support, sepia, bistre and verdigris ink dots systematically belong to the category of lower image contrast values. The other ink dots show image contrasts between these two extremes. Categories listed in Table 1 are therefore confirmed when looking at image contrast, and conclusions about the likelihood of different pigments leading to a good image contrast are reinforced.
Additionally, comparing Figs. 2 and 4, it can be observed that image contrasts are systematically poorer when imaging ink dots on parchment compared to rag paper sheets. This can be explained by the fact that the refractive index of the parchment sheet used as a support is larger (1.71 ± 0.02) than the refractive indices of rag paper sheets, with or without additional sizing (1.29 ± 0.02 and 1.33 ± 0.06). The differences in refractive index values for inks and the surrounding support are therefore reduced in the case of parchment, and it follows that the differences in the amplitude of the pulse reflected at the air-front surface of the sheet interface with or without ink present are lower for parchment. As for rag paper sheets with and without additional sizing, no significant difference in the constant refractive  . 2 Photographs of ink dots on original rag paper and the associated terahertz images acquired in reflection mode. Terahertz images were produced using the maximum amplitude of the reflected waveforms after deconvolution. All images were normalised and are displayed on the same grey scale Fig. 3 Photographs of ink dots on additionally sized rag paper and the associated terahertz images acquired in reflection mode. Terahertz images were produced using the maximum amplitude of the reflected waveforms after deconvolution. All images were normalised and are displayed on the same grey scale Fig. 4 Photographs of ink dots on parchment and the associated terahertz images acquired in reflection mode. Terahertz images were produced using the maximum amplitude of the reflected waveforms after deconvolution. All images were normalised and are displayed on the same grey scale index is observed. However, ink dots applied on the additionally sized rag paper sheet tend to have a broader diameter than ink dots on the original rag paper sheets, while the applied volume of the ink remained the same (2 mL). This implies that sizing modifies the surface properties of paper in such a way that the angle of contact between the ink and the sheet is increased. With the exception of bone black (which seems to have undergone no or little spread in the presence of sizing), contrast values tend to be poorer for the rag paper sheet with additional sizing than for the original rag paper sheet. Pigment concentration therefore plays a role in the image contrast: a larger concentration of pigment on the paper surface leads to an increased amplitude of the reflected waveform. This is further confirmed by the variations of contrast within an ink dot (mostly) coinciding with variations of pigment concentration, as clearly illustrated by the minium ink dot on additionally sized paper in Fig. 3.
When not caused by variations in concentration, variations in contrast within an ink dot can arise from scattering. The configuration of the optics within the imaging set-up is such that the illuminating source is at 30°from the normal to the surface, and a detector is also placed at 30°from the normal to the surface. An edge can be considered as a diffuse reflector: the beam would be reflected in all directions, but the intensity of the radiation (and therefore also the amplitude of the electric field) reflected from this edge would be larger in the direction of the normal to the surface, according to Lambert's cosine law. A schematic diagram describing the distribution of the intensity of the electromagnetic field reflected from the surface of an ink dot is available in Online Resource 3. As a result, the amplitude of the electric field collected from an edge at an oblique angle would be reduced compared to a specular reflection. This would appear dark in the terahertz image produced from the amplitude of the pulse reflected from the front surface of the sheet. This phenomenon is particularly observable for bone black on any support, or for iron gall ink and bistre on the original rag paper: due to scattering, the coffee-stain-like ring left upon drying appears darker than the centre of the ink dot, even though it may represent an area with larger pigment concentration.

Repeatability
To acknowledge the difference between variations of C values inherent to minor variations in shape, concentration or topography of the dots, and variations of C values caused by the different refractive properties of different inks, standard deviation values have been calculated after applying the same ink suspensions of vermilion (colour pigment likely leading to a good image contrast) and bistre (black pigment unlikely to lead to a good image contrast) in triplicate droplets on the same support and imaging them with terahertz TD imaging. In Fig. 5 and Table 2, it can be observed that all bistre ink dots systematically lead to a poor image contrast, while the contrast of vermilion ink dots can vary greatly between triplicates on the same support, likely due to scattering and distribution of ink concentration. Similarly, the local accumulation of minium pigment in the ink dot displayed in Online Resource 2 seems to be at the origin of the larger contrast value reported in Table 2, compared to the contrast value reported for the minium ink dot with an apparent uniform pigment distribution, displayed in Fig. 2. This is the reason why the influence of the refractive index from pigments and supports on the image contrast can only be qualitatively, and not quantitatively, assessed using trends and categories rather than comparing the contrast from ink dots pairwise.

Extraction of Spectral Information from Ink Inscriptions
While the constant refractive index is a good indicator of the reflective properties of a material with respect to the terahertz pulse, which encompasses all frequencies, materials may exhibit a  Photographs of three different bistre and vermilion ink dots applied on rag paper, additionally sized rag paper and parchment, and the associated terahertz images acquired in reflection mode. Terahertz images were produced using the maximum amplitude of the reflected waveforms after deconvolution. All images were normalised and are displayed on the same grey scale dispersive behaviour in the terahertz range. The frequency-dependent refractive index from pigments and dry inks is displayed in Fig. 6. Additionally, their frequency-dependent absorption coefficients are displayed in Fig. 7. The first thing to observe is that the samples of vermilion, minium, verdigris and dry iron gall ink exhibit a reproducible spectral signature between 0.15 and 3 THz. Figure 7 shows that vermilion and minium exhibit intense and narrow absorption bands at the resonance frequencies 38, 42 and 89 cm −1 (1.15, 1.27 and 2.66 THz) and 55, 62, 71 and 84 cm −1 (1.64, 1.86, 2.14 and 2.51 THz), respectively. Iron gall ink shows three broad absorption bands at 38, 51 and 64 cm −1 (1.14, 1.53 and 1.92 THz), and verdigris shows a sharp absorption band at 34 cm −1 (1.02 THz) and a broad absorption band at 76 cm −1 (2.27 THz). The specimens of minium, vermilion, verdigris and iron gall ink studied here are therefore crystalline materials with a resonant behaviour in the terahertz range 0.1-3 THz. It should also be noted that the absorption features from vermilion (cinnabar) and minium are in agreement with the spectral signatures of the supplies from Zecchi, Florence, analysed with a different terahertz TD set-up in transmission and available on the online database thzdb.org [17]. The supply of verdigris from Zecchi, analysed using Fouriertransform infrared spectrometer in the range 0.6-13 THz and available in the same online database also displays absorption bands at ≈ 1 and 2.25 THz, but the signal to noise ratio is poor below 5 THz. The resonant behaviour translates into an Banomalous^dispersion in refractive index values: the refractive index is larger than the resonant frequency at lower frequencies, and it is lower at larger frequencies. If pigments and inks exhibit larger absorption or reflection properties in specific spectral regions, an improved image contrast of the inscriptions could potentially be obtained using the amplitude of the FD signal in these specific spectral regions only. Besides, these features are characteristic of a specific pigment or ink, so the frequencydependent behaviour of the contrast in the FD image could reveal the nature of the ink or pigment. In other words, using the mathematical tool of Fourier-transform to convert TD waveforms into FD spectra, and benefiting from the strong resonant behaviour of some pigments and inks in the terahertz region, terahertz TD imaging could potentially be used as a chemical imaging tool. Figure 8 shows the TD B-scans and raw waveforms from ink dots on rag paper. An interesting phenomenon is observed in the B-scans from iron gall ink and minium ink dots on rag paper and also in the TD B-scan from vermilion shown in Fig. 9: a Bringing^effect is observed in the tail of the TD waveform reflected from an inked area, which translates into Bripples^observed after the back surface of the paper sheet on the TD B-scan, as if additional layers were present under the paper back surface. This phenomenon is systematically observed in all iron gall, minium and vermilion ink dots on rag paper and parchment sheets, although far less pronounced for iron gall ink dots, and is never observed for ink dots made of other pigments. These ink dots are the only pigments with characteristic spectral signatures in the terahertz range (with the exception of verdigris with a weak spectral feature at a low frequency and a very broad a) b) Fig. 8 a Terahertz TD B-scans along a line across an ink dot on rag paper (sepia, malachite, lamp black, iron gall ink or minium mixed with Arabic gum) and b associated raw TD waveforms extracted from pixel positions indicated with numbers. The inset in b shows the tail of the waveforms collected in positions 6 and 7 on a magnified scale. The star indicates the presence of an additional TD feature due to the presence of an ink dot at the surface of the sheet feature at a high frequency). More specifically, the Bringing^and Brippling^phenomena are very clearly observed for the two pigments with very strong and narrow absorption bands: minium and vermilion. This ringing effect is also clearly observed in the tail of the waveforms transmitted through spectroscopic pellets of minium or vermilion mixed with HDPE (data not shown here). The frequency of this Bringing^effect observed in the signal transmitted through spectroscopic pellets was 1.7 THz for pellets of minium and 1.2 THz for pellets of vermilion. These frequencies correspond to spectral regions of high absorption in the absorption coefficients of minium and vermilion shown in Fig. 7. The frequency of the Bringing^effect observed in waveforms reflected from minium and vermilion ink dots on rag paper and parchment sheets showed similar values (between 1.7 and 1.8 THz for minium and 1-1.1 THz for vermilion). It can therefore be assumed that there is a frequency-dependent absorption of the terahertz pulse transmitted through the thin layer of ink applied on paper and reflected at the inkpaper interface and that this frequency dependence is characteristic of the vibrational properties in the region 0.1-3 THz of the pigment. Unfortunately, this Bringing^effect has a weak amplitude and is significantly attenuated or lost when imaging the same ink dots in Bverso^configuration, i.e. with the incident pulse impinging on the air-support interface, then support-ink, and finally ink-air. This is due to the fact that the reflected pulse that has been partially absorbed by the ink layer on the back surface travels through a layer of support before reaching the detector in reflection mode, and is therefore attenuated. This Bringing^effect may also be attenuated or lost during the deconvolution process, depending on the filter function used [47]. Besides, the presence of this Bringing^effect in TD waveforms does not seem to significantly influence the shape of the FD spectra from ink dots on support obtained from the Fourier-transform of the reflected TD waveforms: it is difficult to observe a clear decrease in amplitude at either 1.7 or 1.2 THz for the pixels corresponding to minium and vermilion ink dots, respectively, in the FD B-scans shown in Fig. 9b, d. However, topographic features and morphology of the layered structure have a clear influence on the FD signals [49]. The presence of an ink layer introduces an increase in the optical path, which changes the time delay between the signals reflected from the front and back surfaces of the sheet. This phenomenon and the presence of a new interface create the appearance of an additional layer between the front and back surface of the sheet, i.e. the grey patches observed between 23 and 24 ps in the TD B-scans for pixels within the inked area (Fig. 9a, c). Since pigment distribution within the ink dot, thickness of the ink layer and penetration of the ink into paper fibres are not perfectly uniform across an ink dot, two grey patches are observed for the minium ink dot and three for the vermilion ink dot, instead of a continuous grey layer. The exact same grey patches are observed in the FD B-scans within a frequency range where the rest of the pixels display no (or low) signal (Fig. 9b, d). The FD response therefore seems to be largely influenced by the number of layers in the document and their thickness. To summarise, when imaging a document with terahertz TD technology in reflection mode, the terahertz pulse can undergo dispersive absorption when transmitted through a thin layer of ink on support before it is reflected on the ink-support interface, but it seems that morphological modifications due to the presence of this ink layer influence more the FD response of the document than the resonant behaviour of the materials present within the document.

Conclusions
The influence of the difference in refractive index between the ink and the support underneath and around it on the image contrast produced from the amplitude of the pulse reflected at the air-front surface of the sheet interface has been demonstrated. Sepia ink and inks produced with bistre or verdigris mixed with a solution of Arabic gum or rabbit skin glue are unlikely to lead to well-contrasted images. However, bone black, ivory black, iron gall ink, malachite, lapis lazuli, minium and vermilion are likely to lead to well-contrasted images, while lamp black, carbon black and graphite give the best imaging results. This categorisation enables to provide recommendations to paper and archives conservators on occasions when terahertz TD imaging in reflection mode is a valuable depth-imaging technology to uncover text within a document. For inks containing sepia, bistre or verdigris, the use of terahertz TD imaging in reflection mode is therefore not recommended, unless the set-up has a very large signal to noise ratio and can probe very low differences in reflective properties. In addition to simply adapt the grey scale or brightness scale of the image, image processing algorithms can further be used to improve the image contrast of inscriptions and for character recognition, such as k-means clustering algorithms [68].
While the dispersive behaviour of iron gall ink, vermilion and minium can be probed in reflection mode, the sample morphology (e.g. changes in the number and thickness of layers) appears to influence more the FD response of the document than the resonant behaviour of the materials present within the document. As a result, the use of terahertz TD imaging in reflection mode as a chemical imaging tool for the improvement of image contrast using specific spectral windows or for material identification is currently limited. If there is prior knowledge about the document, such as an estimate of the number and thickness range of the investigated layers, simulating the signal that has propagated through the document could help in recovering the spectral signature from the different materials present within it [69]. However, the number of layers in a sealed archival document is rarely known. Besides, such approach remains very complex and its complexity increases with the number of layers in the document. Terahertz TD imaging in transmission mode may be preferred for chemical imaging, as the FD response would be less influenced by the sample morphology. Yet the set-up needs to show an improved signal to noise ratio in the frequency range of interest in order to detect small changes in the dispersion of the signal caused by the ink [16,70]. Alternatively, a continuous wave imaging set-up tuned at the frequency where the material of interest shows strong absorption or reflection may be used to improve the contrast of any feature containing this material (e.g. terahertz sources tuned at 1.15, 1.27 or 2.66 THz to improve the contrast of vermilion ink lines and tuned at 1.64, 1.86, 2.14 or 2.51 THz for minium ink lines) [71][72][73]. Compact terahertz quantum cascade laser sources or Gunn diodes can be used for this purpose and can further be combined with tomography to recover some depth information.
To further inform conservators of the usefulness of terahertz TD imaging to uncover texts in documents, the terahertz optical properties of synthetic dyes and pigments found in documents from the last two centuries may also be studied, such as aniline, nigrosine, methylene blue and Prussian blue. Finally, this study combined theory, spectroscopy and imaging to provide recommendations on occasions when terahertz TD imaging is a valuable depth-imaging technology to uncover text within a document. To strengthen the reliability of these recommendations, the first step would be to objectively assess if the refractive index and absorption coefficient presented in this study faithfully render the optical properties of the associated pigment. To do so, it would be recommended to perform further analysis using a large number of spectroscopic pellets, but also using different supplies of the same pigment, as well as using different set-ups (round robin test). The accuracy of these optical coefficients would also need to be improved by taking pellet porosity and particle size into account, through the combined use of helium gas pycnometry and effective medium theory [74]. Once a more complete and accurate database of optical coefficients from ink and support materials is designed, each terahertz researcher would need to refine the contrast categories described above to the signal to noise imaging performances of his or her terahertz set-up, in order to provide customed system-related recommendations.