# Construction of Spectral Discoloration Model for Red Lead Pigment by Aging Test and Simulating Degradation Experiment

Authors: Jinxing Liang; Xiaoxia Wan; Chan Li; Junfeng Li; Qiang Liu; Qian Cao
Journal: Journal of Spectroscopy
 (2016)

---

## Abstract

The construction of spectral discoloration model, based on aging test and simulating degradation experiment, was proposed to detect the aging degree of red lead pigment in ancient murals and to reproduce the spectral data supporting digital restoration of the ancient murals. The degradation process of red lead pigment under the aging test conditions was revealed by X-ray diffraction, scanning electron microscopy, and spectrophotometer. The simulating degradation experiment was carried out by proportionally mixing red lead and lead dioxide with referring to the results of aging test. The experimental result indicated that the pure red lead was gradually turned into black lead dioxide, and the amount of tiny particles of the aging sample increased faced with aging process. Both the chroma and lightness of red lead pigment decreased with discoloration, and its hue essentially remains unchanged. In addition, the spectral reflectance curves of the aging samples almost started rising at about 550 nm with the inflection moving slightly from about 570 nm to 550 nm. The spectral reflectance of samples in long- and in short-wavelength regions was fitted well with the logarithmic and linear function. The spectral discoloration model was established, and the real aging red lead pigment in Dunhuang murals was measured and verified the effectiveness of the model.


---

## Body

## 1. Introduction

Lead-based pigments have been widely employed and confirmed to exist in various fields of art [1–3]. Among the lead-based pigments, red lead pigment, as the bright one, was chosen by most ancient painters to enrich the tones on a large variety of artworks, including medieval manuscripts [4], wood-block print [5], canvas paintings [6], or wall paintings [7]. In most cases, however, the red lead pigment painted layers have been turned into brown or complete black, which is caused by humidity, light (mainly ultraviolet), gas environment, and other aspects [8–10]. The artistic effect of the ancient artworks was seriously affected by the degradation of red lead pigment, and people could not appreciate the original status of these precious heritages. At the same time, the discoloration of the ancient artworks also hindered the excavation of the historical information about the ancient societies. Thus, it is necessary and valuable to develop the method to digitally recover the darkened artworks to their original status, wherein the restoration of the darkened artworks should be based on the discoloration model and the discoloration model needs to satisfy the principle that the reproduced data are consistent with the real aging red lead pigment on artworks as much as possible.Modern analysis devices and technique provide the potential to systematically study the degradation process of the red lead pigment and to construct the discoloration models in different aging conditions. Over the past few decades, many researches have focused on revealing the mechanism of red lead degradation [2, 4, 5, 8, 9], but only a few studies have been involved in construction of the red lead discoloration model [11–13]. Shi and Lu firstly introduced the colorimetric and chemical models based aging simulation of Dunhuang murals for restoration and for audience reexperiencing the aging process of the murals [11]. And lately, they further formalized the empirical knowledge from artists in their models to improve the simulation results [12]. Recently, Zhao et al. studied the mechanism of red lead degradation and constructed the corresponding time and color functions for the aging process but did not validate the established functions with real aging red lead in murals [13]. The main purpose of constructing the discoloration model is to restore the fading artworks to their original colors in correct way. However, current discoloration models for red lead pigment are all based on the colorimetric theory, which were constructed in RGB or Lab color spaces and inevitably influenced by metameric problem [14].Spectral reflectance, which is robust to illumination and observer variations, is regarded at the “fingerprint” of the object color as the color parameters are calculated from spectral reflectance under any illuminations and observer conditions based on colorimetric theory. At the same time, spectral reflectance also reflects the chemical and physical properties of the object as the absorption/reflectance features of an object depend primarily on the chemical properties of its components and secondly on the physical properties. The technology of spectral imaging developed rapidly in the field of cultural heritages protection for the past decades because of its high-quality imaging abilities and its safety and effectiveness in pigment identification [15–17]. If the spectral discoloration model was constructed for red lead pigment, the virtual restoration of the ancient murals would be more scientific using the spectral images. In this study, the degradation process of red lead pigment in Dunhuang murals was studied by aging test and simulating degradation experiment. The spectral discoloration model of red lead pigment was constructed, and the validity of the model was verified by calculating and predicting the spectral reflectance of the real aging red lead pigment in Dunhuang murals.


## 2. Materials and Methods

### 2.1. Material

Generally, murals in Dunhuang grottoes consist of four layers: the coarse clay layer, the fine clay layer, the plaster layer, and the pigment layer. The coarse clay layer contains the mixture of straw and sand in an appropriate ratio to fill up the cliff. The fine clay layer contains moderate cotton or flax fibers to fill the cracks of the dry coarse clay layer as a flat plane. The plaster layer is used to cover the fine clay layer for paintings, and the pigment layer is lastly painted on the plaster layer by Chinese goat-hair brush to complete the construction of the wall paintings, in which the pigment is usually dispersed in animal glue (mostly in cowhide glue). In this study, the industrial white cardboard is applied as substrate to simulate the plaster of the murals since the chemical properties of plaster layer are stable. The pure red lead pigment is well-dispersed in the cowhide glue and painted carefully on the substrate layer to create the experimental samples. The pure red lead pigment, the animal gelatin, and the tools painting mural sample were supplied by Dunhuang Academy (Dunhuang, China). In addition, they also imparted the skills of preparing the animal gelatin solution and the techniques of painting the mural sample. The pure lead dioxide used for simulating degradation experiment was purchased from the Shenshi Chemical Co., Ltd. (Wuhan, China).


### 2.2. Apparatus

The ultraviolet (UV) light and humidity are the main factors causing the discoloration of red lead in Dunhuang murals [8, 10–13]. The UV aging experiment box ZN-P with high pressure mercury lamp and the constant-temperature-and-humidity tested box DHS-500 were chosen to carry out the aging test. The wavelength of the high pressure 40 w mercury lamp ranges from 200 to 1400 nm with the dominant emission wavelength at 365 nm. The humidity and temperature ranges of the two types aging equipment could be controlled at 30~98% RH and 0~150°C. The X-ray diffraction (XRD) analyzer Bruker Smart Apex CCDII was adopted to analyze chemical phases of the aging samples. The microstructure analysis of the aging samples was completed using Quanta 200 scanning electron microscope (SEM). The spectral reflectance of the sample ranging from 360 to 750 nm at 10 nm intervals was measured using GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in specular component excluded (SCE) model. Furthermore, the electronic analytical balance METTLER TOLEDO AB135-S/FACT, measurement accuracy of which is 0.01/0.1 mg, was used to measure the mass of the gelatin and red lead pigment when preparing the gelatin solution and painting a mural sample. The real aging red lead pigment points were measured by the PR705 spectroradiometer with the standard whiteboard produced by Photo Research Inc. (American).


### 2.3. Spectral Measurement and Color Expression

Prior to the formal experiment, the spectral measurement accuracy and repeatability as well as the consistency between two types of spectral measurement devices were evaluated to ensure the validity of the spectral measurement results. Six saturated standard color samples of GretagMacbeth ColorCkher Chart (red, green, blue, yellow, magenta, and cyan) were used to carry out the evaluation procedure. The root-mean-square error (RMSE) was calculated between the measured spectral and the corresponding certified values. Each measurement was repeated fifty times at the interval of twenty seconds. Statistical results of the accuracy, repeatability, and consistency are summarized in Table1. The mean RMSE for all the color samples measured by 7000A and PR705 is less than 0.49% and the maximum RMSE is less than 0.81%. These results indicate that there is good accuracy and repeatability for both spectral measurement devices. The mean and maximum RMSE between 7000A and PR705 are 0.53% and 0.85%, respectively, which is slightly higher than the accuracy and repeatability of each device alone but still maintains a fairly good consistency.Table 1
Statistical results of spectral measurement accuracy, repeatability, and consistency of Color-Eye 7000A spectrophotometer and PR705 spectroradiometer.
RMS (%)
Red
Green
Blue
Yellow
Magenta
Cyan
7000A
Mean
0.34
0.35
0.44
0.35
0.46
0.37
Max
0.49
0.54
0.69
0.72
0.71
0.61
Min
0.22
0.22
0.22
0.31
0.34
0.21
PR705
Mean
0.37
0.38
0.49
0.42
0.49
0.42
Max
0.56
0.66
0.68
0.69
0.81
0.76
Min
0.26
0.21
0.24
0.33
0.37
0.23
PR705 versus 7000A
Mean
0.44
0.42
0.47
0.42
0.53
0.41
Max
0.79
0.81
0.78
0.82
0.85
0.78
Min
0.36
0.34
0.37
0.32
0.39
0.36After the spectral reflectance was acquired, the color parameter of samples was calculated under the CIE standard illuminant D50 and the CIE 1931 standard observer. The spectral reflectance was firstly converted to the CIEX
Y
Z tristimulus according to (1), where x
(
λ
), y
(
λ
), and z
(
λ
) represent color-matching functions, E
(
λ
) represents the spectral power distribution of the CIE standard illuminant, S
(
λ
) is the measured spectral reflectance of the samples, and k is the adjustment factor. And then the chromaticity values of the samples like Y
x
y, L
∗
a
∗
b
∗, or sRGB were calculated from the CIE X
Y
Z tristimulus, and the transformation from X
Y
Z to Y
x
y, L
∗
a
∗
b
∗, and sRGB values can be referenced in the literature [18]:(1)
X
=
k
∫
λ
x
λ
E
λ
S
λ
d
λ
,
Y
=
k
∫
λ
y
λ
E
λ
S
λ
d
λ
,
Z
=
k
∫
λ
z
λ
E
λ
S
λ
d
λ
,
k
=
100
∑
λ
y
λ
E
λ
d
λ
.Besides, the color of the experiment samples can be presented by Matlab toolbox OptProp as the sRGB values of the samples can be rendered as uniformity color patches by the function of OptImage [19]. The color difference between the samples was calculated by the CIEDE2000 (DE00) formula [20].


### 2.4. Parameters Setting for Aging Test

Six red lead samples were painted for the aging test. The temperature and humidity were settled at the level of 60°C and 80% RH for the aging tests. The test boxes were maintained in a stable working state. The red lead pigment samples were placed alternately in two test boxes with the interval of 12 hours starting from the constant-temperature-and-humidity test box. From beginning of aging test, one sample was picked out every five days, until the aging test was completed. Samples of 0 days, 5 days, 10 days, 15 days, 20 days, and 25 days were denoted as sample 1 (1#), sample 2 (2#), sample 3 (3#), sample 4 (4#), sample 5 (5#), and sample 6 (6#). Each sample was measured and analyzed using the professional equipment mentioned above. The spectral reflectance of the samples was measured in a constant-temperature-and-humidity optical laboratory with temperature at 25°C and humidity at 40% RH.


## 2.1. Material

Generally, murals in Dunhuang grottoes consist of four layers: the coarse clay layer, the fine clay layer, the plaster layer, and the pigment layer. The coarse clay layer contains the mixture of straw and sand in an appropriate ratio to fill up the cliff. The fine clay layer contains moderate cotton or flax fibers to fill the cracks of the dry coarse clay layer as a flat plane. The plaster layer is used to cover the fine clay layer for paintings, and the pigment layer is lastly painted on the plaster layer by Chinese goat-hair brush to complete the construction of the wall paintings, in which the pigment is usually dispersed in animal glue (mostly in cowhide glue). In this study, the industrial white cardboard is applied as substrate to simulate the plaster of the murals since the chemical properties of plaster layer are stable. The pure red lead pigment is well-dispersed in the cowhide glue and painted carefully on the substrate layer to create the experimental samples. The pure red lead pigment, the animal gelatin, and the tools painting mural sample were supplied by Dunhuang Academy (Dunhuang, China). In addition, they also imparted the skills of preparing the animal gelatin solution and the techniques of painting the mural sample. The pure lead dioxide used for simulating degradation experiment was purchased from the Shenshi Chemical Co., Ltd. (Wuhan, China).


## 2.2. Apparatus

The ultraviolet (UV) light and humidity are the main factors causing the discoloration of red lead in Dunhuang murals [8, 10–13]. The UV aging experiment box ZN-P with high pressure mercury lamp and the constant-temperature-and-humidity tested box DHS-500 were chosen to carry out the aging test. The wavelength of the high pressure 40 w mercury lamp ranges from 200 to 1400 nm with the dominant emission wavelength at 365 nm. The humidity and temperature ranges of the two types aging equipment could be controlled at 30~98% RH and 0~150°C. The X-ray diffraction (XRD) analyzer Bruker Smart Apex CCDII was adopted to analyze chemical phases of the aging samples. The microstructure analysis of the aging samples was completed using Quanta 200 scanning electron microscope (SEM). The spectral reflectance of the sample ranging from 360 to 750 nm at 10 nm intervals was measured using GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in specular component excluded (SCE) model. Furthermore, the electronic analytical balance METTLER TOLEDO AB135-S/FACT, measurement accuracy of which is 0.01/0.1 mg, was used to measure the mass of the gelatin and red lead pigment when preparing the gelatin solution and painting a mural sample. The real aging red lead pigment points were measured by the PR705 spectroradiometer with the standard whiteboard produced by Photo Research Inc. (American).


## 2.3. Spectral Measurement and Color Expression

Prior to the formal experiment, the spectral measurement accuracy and repeatability as well as the consistency between two types of spectral measurement devices were evaluated to ensure the validity of the spectral measurement results. Six saturated standard color samples of GretagMacbeth ColorCkher Chart (red, green, blue, yellow, magenta, and cyan) were used to carry out the evaluation procedure. The root-mean-square error (RMSE) was calculated between the measured spectral and the corresponding certified values. Each measurement was repeated fifty times at the interval of twenty seconds. Statistical results of the accuracy, repeatability, and consistency are summarized in Table1. The mean RMSE for all the color samples measured by 7000A and PR705 is less than 0.49% and the maximum RMSE is less than 0.81%. These results indicate that there is good accuracy and repeatability for both spectral measurement devices. The mean and maximum RMSE between 7000A and PR705 are 0.53% and 0.85%, respectively, which is slightly higher than the accuracy and repeatability of each device alone but still maintains a fairly good consistency.Table 1
Statistical results of spectral measurement accuracy, repeatability, and consistency of Color-Eye 7000A spectrophotometer and PR705 spectroradiometer.
RMS (%)
Red
Green
Blue
Yellow
Magenta
Cyan
7000A
Mean
0.34
0.35
0.44
0.35
0.46
0.37
Max
0.49
0.54
0.69
0.72
0.71
0.61
Min
0.22
0.22
0.22
0.31
0.34
0.21
PR705
Mean
0.37
0.38
0.49
0.42
0.49
0.42
Max
0.56
0.66
0.68
0.69
0.81
0.76
Min
0.26
0.21
0.24
0.33
0.37
0.23
PR705 versus 7000A
Mean
0.44
0.42
0.47
0.42
0.53
0.41
Max
0.79
0.81
0.78
0.82
0.85
0.78
Min
0.36
0.34
0.37
0.32
0.39
0.36After the spectral reflectance was acquired, the color parameter of samples was calculated under the CIE standard illuminant D50 and the CIE 1931 standard observer. The spectral reflectance was firstly converted to the CIEX
Y
Z tristimulus according to (1), where x
(
λ
), y
(
λ
), and z
(
λ
) represent color-matching functions, E
(
λ
) represents the spectral power distribution of the CIE standard illuminant, S
(
λ
) is the measured spectral reflectance of the samples, and k is the adjustment factor. And then the chromaticity values of the samples like Y
x
y, L
∗
a
∗
b
∗, or sRGB were calculated from the CIE X
Y
Z tristimulus, and the transformation from X
Y
Z to Y
x
y, L
∗
a
∗
b
∗, and sRGB values can be referenced in the literature [18]:(1)
X
=
k
∫
λ
x
λ
E
λ
S
λ
d
λ
,
Y
=
k
∫
λ
y
λ
E
λ
S
λ
d
λ
,
Z
=
k
∫
λ
z
λ
E
λ
S
λ
d
λ
,
k
=
100
∑
λ
y
λ
E
λ
d
λ
.Besides, the color of the experiment samples can be presented by Matlab toolbox OptProp as the sRGB values of the samples can be rendered as uniformity color patches by the function of OptImage [19]. The color difference between the samples was calculated by the CIEDE2000 (DE00) formula [20].


## 2.4. Parameters Setting for Aging Test

Six red lead samples were painted for the aging test. The temperature and humidity were settled at the level of 60°C and 80% RH for the aging tests. The test boxes were maintained in a stable working state. The red lead pigment samples were placed alternately in two test boxes with the interval of 12 hours starting from the constant-temperature-and-humidity test box. From beginning of aging test, one sample was picked out every five days, until the aging test was completed. Samples of 0 days, 5 days, 10 days, 15 days, 20 days, and 25 days were denoted as sample 1 (1#), sample 2 (2#), sample 3 (3#), sample 4 (4#), sample 5 (5#), and sample 6 (6#). Each sample was measured and analyzed using the professional equipment mentioned above. The spectral reflectance of the samples was measured in a constant-temperature-and-humidity optical laboratory with temperature at 25°C and humidity at 40% RH.


## 3. Result and Discussion

### 3.1. Aging Test

The color appearance of red lead pigment changed from bright oranges to dim blacks as shown in Figure1. The difference in color appearance between the two adjacent samples from 1# to 6# is gradually getting smaller from visual perception. The detailed information for characterization and analysis of the chemical composition, microstructure, spectral reflectance, and chromaticity values of the aging samples were described in the following parts.Figure 1
Samples in different aging degrees: 0 days (1#), 5 days (2#), 10 days (3#), 15 days (4#), 20 days (5#), and 25 days (6#).


#### 3.1.1. Chemical Composition

As illustrated in Figure2, XRD results indicated that the aging product of red lead under the aging test environment was lead dioxide, explained by the increased diffraction peaks of PbO2 and the decreased diffraction peaks of Pb3O4. This is consistent with the previous researches [2, 4, 5]. Furthermore, the diffraction intensity ratio of PbO2 to Pb3O4 in Table 2 and the discrimination of color between the adjacent samples in Figure 1 indicated the strong color covering power of lead dioxide. A little amount of lead dioxide may darken the sample remarkably, but when the amount of lead dioxide reaches a certain degree, its color covering ability will reach the maximum.Table 2
The diffraction intensity ratio of PbO2 to Pb3O4 of the aging samples.
Sample
1#
2#
3#
4#
5#
6#
Intensity ratio (I
P
b
O
2/I
P
b
3
O
4)
0
0.048
0.072
0.138
0.258
1.503Figure 2
XRD results for the aging samples:(
) represents the diffraction peak of PbO2. ∗ represents the diffraction peak of Pb3O4.


#### 3.1.2. Microstructure

The microstructure image of the aging samples captured with SEM is shown in Figure3. It demonstrated that the amount of the large particle is decreasing and the amount of tiny particle is increasing faced with the aging process, but there are still some large particles wrapped by tiny particles. The grain size distribution of the aging samples measured through analysis software was summarized in Table 3. The statistics shows that the maximum, minimum, and mean particle size are reduced as the aging time passes, and the minimum particle size is ultimately stable at the level of 0.21 μm. Based on the particle size changing trend described above, the minimum particles in the aging samples may be the lead dioxide. The phenomenon of the large particles wrapped by tiny particles may be the reason why the previous researchers observed the incomplete oxidized red lead pigment grains in the cross section of the pigment layer painted by red lead pigment [7, 8].Table 3
Particles distribution of the aging samples.
Sample
1#
2#
3#
4#
5#
6#
Max/(μm)
8.21
5.41
3.62
2.52
2.41
2.12
Min/(μm)
0.63
0.5
0.4
0.38
0.21
0.21
Mean/(μm)
2.12
1.78
1.36
1.03
0.79
0.66Figure 3
Microstructure of the six aging samples.Moreover, the optic path of the incident lights becomes increasingly more complex owing to the reduced large particles and the increased tiny particles of the aging process. The complexity of the sample’s microstructure enhanced the optic path of the incident lights by increasing its frequency of reflection according to light propagation theory. Coupled with the strong light absorption of PbO2, the amount of emergent light was inevitably reduced and the samples dramatically grew darker over aging time as shown in Figure 1.


#### 3.1.3. Spectral Reflectance and Chromaticity

The spectral reflectance of the aging samples was measured with GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in SCE model. As illustrated in Figure4(a), the spectral reflectance of the aging samples of long-wavelength from 550 to 750 nm apparently decreased with the aging time, while the decreased value of the spectral reflectance in the same aging interval becomes smaller and smaller. It is suggested that the darkening trends of samples become slower, but this does not reveal that the chemical composition transition also obeys the same law. The intensity ratio of PbO2 to Pb3O4 in Table 2 shows that the chemical composition transition rates exhibit a thoroughly inverse trend. The inconsistency between sample darkening trends and chemical composition transition rate from red lead to lead dioxide may be mainly caused by the strong color covering power of lead dioxide and its proportion in the aging samples as described above. According to the spectral reflectance change trends, we preliminary infer that the logarithmic function could be well applied to construct the red lead pigment discoloration model under the aging conditions.Figure 4
Distribution of the spectral reflectance of the aging samples (a) and the corresponding first derivate (b).
(a)
(b)In addition, it was easy to find that the spectral reflectance of the aging samples almost has the same rising point in about 550 nm, and the faint decrease in short-wavelength regions is from 360 to 550 nm. As the spectral reflectance of the aging sample is primarily determined by the chemical properties of its components, the spectral features of the aging samples could be explained from the chemical components aspect as the aging samples share the same chemical components except sample 1#. Red lead absorbs the light from 360 to 550 nm and reflects the rest in visible light region, while lead dioxide absorbs almost all visible light. The different ratios of red lead and lead dioxide contained in the aging samples lead to the almost same rising point and the low reflection in short-wavelength region. Furthermore, as shown in Figure4(b), the different ratios of red lead and lead dioxide contained in the aging samples also influence the inflection of the spectral curve. With the amount of lead dioxide increased in the aging samples, the spectral inflection of samples is slightly moved from about 570 to 550 nm.The chromaticity and color difference of the aging samples were summarized in Table4. The brightness value of the aging samples Y and L
∗ as well as the chroma value C
∗ decreased rapidly at the beginning of the experiment, but the speed slows down later in the experiment. Similarly, the color difference DE00 increased rapidly at the beginning of the experiment and the growth slackened later in the experiment. The changing trends of the chromaticity values consisted with the trend of the spectral reflectance in the long-wavelength regions. The reason for such phenomenon was that the chromaticity is linearly calculated from spectral reflectance multiplied by spectral power distributions of the illumination and by color-matching functions [19].Table 4
Chromaticity and color difference of the aging samples.
Sample
Color coordinate
Chroma (C
∗)
D
E
00
Δ
L
∗
Δ
a
∗
Δ
b
∗
Δ
D
E
00
Y
x
y
L
∗
a
∗
b
∗
1#
30.935
0.577
0.371
62.452
52.263
50.530
72.71
0
—
—
—
—
2#
16.951
0.535
0.374
48.198
38.825
36.977
53.62
14.32
14.254
13.438
13.553
14.32
3#
11.300
0.503
0.374
40.081
28.550
25.495
38.28
24.37
8.117
10.275
11.482
10.05
4#
8.241
0.474
0.367
34.479
18.997
16.983
23.59
31.44
5.602
9.533
8.512
7.07
5#
6.454
0.451
0.363
30.530
13.478
12.879
18.64
35.62
3.949
5.499
4.104
4.18
6#
5.920
0.439
0.361
29.209
12.961
12.133
17.75
36.67
1.321
0.517
0.746
1.05The hue of the aging samples ina
∗
-
b
∗ planes and in CIE 1931 chromaticity diagram was plotted in Figure 5, where a
∗ represents the red-green response and b
∗ represents the yellow-blue response from human visual system in CIE L
∗
a
∗
b
∗ uniform color space, respectively. Figure 5(a) indicated that the hue changed trends can be fitted very well with the linear model, and Figure 5(b) shows the same varying tendencies. The values of Δ
a
∗, Δ
b
∗, and Δ
L
∗ were also summarized in Table 4, which shows the same changing trends as ΔDE00. Based on the data changing trends of Δ
a
∗, Δ
b
∗, Δ
L
∗, and C
∗ in Table 4 and the hue distributions of aging samples plotted in Figure 5, it could be concluded that, during the red lead pigment aging process, the main changing chromaticity feature was the brightness and chroma value, and the hue of the aging samples almost remains unchanged.Figure 5
Hues of aging samples ina
∗
-
b
∗ planes (a) and CIE 1931 chromaticity diagram (b).
(a)
(b)In the next part, the simulating degradation experiment was carried out with referring to the result of aging test. And the spectral reflectance and chromaticity of the aging samples were further analyzed for construction of the spectral discoloration model.


### 3.2. Simulating Degradation Experiment

In this experiment, eleven experimental samples were painted using the painting method simulating Dunhuang murals. The samples were mixed by red lead and lead dioxide in different mixing ratios with total mass controlled at 0.6 g as the quantity is enough for painting a sample in this study, and the mass change step of the pure red lead and lead dioxide was set at 0.06 g. The name and the corresponding mixing ratio of the simulating degradation samples are listed in Table5.Table 5
Names and mixing ratios of the simulating degradation samples.
Sample
S1
S2
S3
S4
S5
S6
S7
S8
S9
S10
S11
Pb3O4/g
0.60
0.54
0.48
0.42
0.36
0.30
0.24
0.18
0.12
0.06
0.00
PbO2/g
0.00
0.06
0.12
0.18
0.24
0.30
0.36
0.42
0.48
0.54
0.60The spectral reflectance of the simulating degradation samples was also measured by GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in SCE model, and the chromaticity of each simulating degradation sample was calculated using Matlab toolbox OptProp. The color of each simulating degradation sample is rendered in Figure6. The spectral reflectance and the corresponding first derivate of the simulating degradation samples are plotted in Figure 7, and the chromaticity value of the simulating degradation samples is summarized in Table 6.Table 6
Chromaticity features of the simulating degradation samples.
Sample
Color coordinate
Chroma (C
∗)
D
E
00
Δ
L
∗
Δ
a
∗
Δ
b
∗
Δ
D
E
00
Y
x
y
L
∗
a
∗
b
∗
S1
30.935
0.577
0.371
62.452
52.263
50.530
72.71
0
—
—
—
—
S2
16.521
0.458
0.366
47.651
26.696
25.834
37.14
17.35
14.801
25.567
24.696
17.35
S3
13.228
0.426
0.361
43.104
19.287
18.466
26.69
23.68
4.547
7.409
7.368
6.33
S4
10.788
0.399
0.354
39.227
14.140
12.822
19.08
29.05
3.877
5.147
5.644
5.37
S5
10.139
0.385
0.351
38.094
11.601
10.349
15.54
31.00
1.133
2.539
2.473
1.95
S6
8.425
0.377
0.349
34.858
9.590
8.676
12.92
34.22
3.236
2.011
1.673
3.22
S7
7.402
0.361
0.345
32.712
6.923
6.237
9.31
36.78
2.146
2.667
2.439
2.56
S8
6.625
0.350
0.343
30.944
4.977
4.770
6.89
38.66
1.768
1.946
1.467
1.88
S9
6.272
0.336
0.338
30.099
3.108
2.701
4.12
40.14
0.845
1.869
2.069
1.48
S10
5.650
0.327
0.336
28.512
1.782
1.633
2.42
41.65
1.587
1.326
1.068
1.51
S11
5.225
0.316
0.332
27.375
0.276
0.345
0.44
43.03
1.137
1.506
1.288
1.38Figure 6
Simulating degradation samples.Figure 7
Distribution of the spectral reflectance of the simulating degradation samples (a) and the corresponding first derivate (b).
(a)
(b)It indicated that there is a strong similarity between the aging test and simulating degradation experiment about the changing trend of the spectral reflectance, the first derivate of spectral reflectance, and the chromaticity values. However, there is slight difference between the spectral reflectance curves shapes of the two experiments according to the results of gelatin degradation in Dunhuang murals carried out by Li [21]. The humidity and light are two important factors causing the aging of the gelatin which lead to the gelatin solution from almost transparent status gradually becoming light yellow, so it can be inferred that the slight difference between the spectral reflectance curves shapes may be caused by the aging of the gelatin in the samples. However, the similarity of spectral and chromaticity features of the samples between the aging test and simulating degradation experiment inferred that the aging of the gelatin just has a very small influence on the aging trends of red lead pigment. The hue distribution of the simulating degradation samples in a
∗
-
b
∗ planes and CIE 1931 chromaticity diagram was shown in Figure 8. The hue distribution also maintains high similarity with the aging test result as it can be fitted very well using the linear model. The similarity of the two experimental results may also prove that the mechanism of the red lead pigment discoloration in Dunhuang murals is caused by the process of red lead altered in lead dioxide.Figure 8
Hue distribution of the simulating degradation samples ina
∗
-
b
∗ planes (a) and CIE 1931 chromaticity diagram (b).
(a)
(b)


### 3.3. Spectral Discoloration Model Construction

The decreased trend of spectral reflectance in long-wavelength regions of the simulating degradation samples was fitted with some simple linear or nonlinear functions. Here the mass ratio of lead dioxide in total mass was used to monitor the change of spectral reflectance. Figure9 shows the decreasing trend of the spectral reflectance at 700 nm. It can be seen from Figure 9 that the decreased trend could be fitted very well using the logarithmic model and that R
2 value of the fitting equation is high up to 0.99515. Moreover, the average R
2 value from 560 to 750 nm is higher than 0.97, and the minimum R
2 value is 0.95329.Figure 9
Decreasing trend of the spectral reflectance of the simulating degradation samples at 700 nm.The other simple functions were also adopted to fit the changing trend of the simulating degradation samples.R
2 values of these fitting functions and the logarithmic fitting function are summarized in Table 7. R
2 values in Table 7 show that the logarithmic function has the highest accuracy which is consistent with the previous hypothesis. Furthermore, the decreasing trend in short-wavelength regions from 360 to 550 nm acquires the best fitting accuracy using the linear function. The fitting result shows that the average R
2 value from 360 to 550 nm is higher than 0.97 and the minimum R
2 value is 0.95993. Thus the spectral discoloration model could be constructed with the abovementioned results. Besides, the spectral reflectance is regarded as the “fingerprint” of the object color and all the chromaticity values are calculated from the spectral reflectance under different illuminating and viewing conditions, so there is no need to construct the discoloration model in color space.Table 7
R
2 values of different fitting functions at 700 nm.
Functions
R
2
Exponential
0.92245
Linear
0.63268
Logarithmic
0.99515
Polynomial(
2
)
∗
0.84378
Power
0.82483
∗Note: “polynomial (
2
)” represents the two-order polynomial.Based on the above analysis, the spectral discoloration model of red lead pigment was constructed as follows:(2)
Y
i
=
A
x
+
B
,
if
1
≤
i
<
22
,
A
ln
⁡
x
+
B
,
if
22
≤
i
≤
40
,where i represents the number of the sampling points from 360 to 750 nm at 10 nm intervals; x is mass ratios of lead dioxide in the simulating degradation samples; Y
i denotes the spectral reflectance in different sampling point for the specific x; and A and B contain the model parameters in each sampling wavelength from 360 to 750 nm at 10 nm interval; namely, (3)
A
=
A
1
⋮
A
i
⋮
A
40
,
B
=
B
1
⋮
B
i
⋮
B
40
.After the spectral discoloration model was constructed, its prediction accuracy was firstly examined by the simulating degradation samples with the known aging degree ofx. The results show that the constructed spectral discoloration model is effective as the average spectral root-mean-square error, the average spectral goodness-of-fit coefficient [22], and the average color difference DE00 of the eleven simulation aging samples are all kept in a fairly good level at 0.0591, 0.9979, and 1.14, respectively.


### 3.4. Real Aging Red Lead Pigment in Dunhuang Murals

The red lead pigment is widely used in Dunhuang murals, especially for the flesh of the Buddha in murals. Most of the murals painted with red lead pigment are black today. But there is still some incompletely discolored red lead pigment identified by chemical analysis of previous researchers [8, 23]. Several red lead pigment points in different aging degree in two typical grottoes from different dynasties were measured with the in situ nondestructive measurement method. The validity of the model was verified by calculating and predicting the spectral reflectance of the real aging red lead pigment.The PR705 spectroradiometer with the standard whiteboard produced by American Photo Research was used to measure the spectral reflectance of the real aging red lead pigment in Dunhuang murals, and the actual measured work scene is shown in Figure10. The insert in Figure 10 represents the geometrical conditions of illumination and measurement. The points measured by PR705 spectroradiometer were marked in Figure 11. The measured spectral reflectance of the five marked points was plotted in Figure 12.Figure 10
Spectral reflectance measured using PR705 spectroradiometer in cave 465 (Yuan Dynasty).Figure 11
Five real aging red lead pigment points in Dunhuang murals: (a) Yuan Dynasty cave 465 and (b) Early Tang Dynasty cave 220.
(a)
(b)Figure 12
Distribution of spectral reflectance of the marked five points in the figure.Because of the larger measurement noise in the short-wavelength below 400 nm and the long-wavelength above 700 nm, the effective wavelength range of the sampled points was restricted to 400~700 nm as plotted in Figure12. It can be seen that there is a very strong similarity between the spectral reflectance shapes of the real aging red lead pigment and the above two experiment samples. The hue distribution of the real aging red lead points also can be fitted well with the linear model with the slope, the intercept, and R
2 value being equal to 0.9692, 0.5621, and 0.9597. Based on the constructed spectral discoloration model, the simulated aging degree of the five real aging red lead pigment points was calculated as 0.8028%, 31.4346%, 5.2151%, 9.7909%, and 73.4961% using the inversed spectral discoloration model.After the detection of aging degree, the spectral reflectance of the five points was predicted using the proposed spectral discoloration model as plotted in Figure13 with the dash-dot line. The prediction accuracy was also evaluated by the metrics of spectral root-mean-square error, the spectral goodness-of-fit coefficient, and the color difference DE00. The predicted results of the five points are summarized in Table 8. The average spectral root-mean-square error and color difference DE00 are 1.21% and 2.64, respectively. The average spectral goodness-of-fit coefficient is high up to 99.75%. The predicted results indicated that the constructed spectral discoloration model could well predict the spectral reflectance of the real aging red lead pigment in Dunhuang murals based on the simulated aging degree. It can be concluded that the constructed spectral discoloration model for red lead pigment could work well in detecting the aging degree of red lead pigment in ancient murals and reproducing the spectral data for specific aging degree.Table 8
The predicted accuracy of the five real aging red lead pigment points.
Point number
RMS (%)
DE00
GFC (%)
1
1.50
1.99
99.93
2
1.30
3.24
99.59
3
1.28
2.61
99.88
4
1.14
2.87
99.86
5
0.81
2.48
99.47
Mean
1.21
2.64
99.75Figure 13
Measured (solid line) and predicted (dash-dot line) spectral reflectance of the five real aging red lead pigment points using the spectral discoloration model.


## 3.1. Aging Test

The color appearance of red lead pigment changed from bright oranges to dim blacks as shown in Figure1. The difference in color appearance between the two adjacent samples from 1# to 6# is gradually getting smaller from visual perception. The detailed information for characterization and analysis of the chemical composition, microstructure, spectral reflectance, and chromaticity values of the aging samples were described in the following parts.Figure 1
Samples in different aging degrees: 0 days (1#), 5 days (2#), 10 days (3#), 15 days (4#), 20 days (5#), and 25 days (6#).


### 3.1.1. Chemical Composition

As illustrated in Figure2, XRD results indicated that the aging product of red lead under the aging test environment was lead dioxide, explained by the increased diffraction peaks of PbO2 and the decreased diffraction peaks of Pb3O4. This is consistent with the previous researches [2, 4, 5]. Furthermore, the diffraction intensity ratio of PbO2 to Pb3O4 in Table 2 and the discrimination of color between the adjacent samples in Figure 1 indicated the strong color covering power of lead dioxide. A little amount of lead dioxide may darken the sample remarkably, but when the amount of lead dioxide reaches a certain degree, its color covering ability will reach the maximum.Table 2
The diffraction intensity ratio of PbO2 to Pb3O4 of the aging samples.
Sample
1#
2#
3#
4#
5#
6#
Intensity ratio (I
P
b
O
2/I
P
b
3
O
4)
0
0.048
0.072
0.138
0.258
1.503Figure 2
XRD results for the aging samples:(
) represents the diffraction peak of PbO2. ∗ represents the diffraction peak of Pb3O4.


### 3.1.2. Microstructure

The microstructure image of the aging samples captured with SEM is shown in Figure3. It demonstrated that the amount of the large particle is decreasing and the amount of tiny particle is increasing faced with the aging process, but there are still some large particles wrapped by tiny particles. The grain size distribution of the aging samples measured through analysis software was summarized in Table 3. The statistics shows that the maximum, minimum, and mean particle size are reduced as the aging time passes, and the minimum particle size is ultimately stable at the level of 0.21 μm. Based on the particle size changing trend described above, the minimum particles in the aging samples may be the lead dioxide. The phenomenon of the large particles wrapped by tiny particles may be the reason why the previous researchers observed the incomplete oxidized red lead pigment grains in the cross section of the pigment layer painted by red lead pigment [7, 8].Table 3
Particles distribution of the aging samples.
Sample
1#
2#
3#
4#
5#
6#
Max/(μm)
8.21
5.41
3.62
2.52
2.41
2.12
Min/(μm)
0.63
0.5
0.4
0.38
0.21
0.21
Mean/(μm)
2.12
1.78
1.36
1.03
0.79
0.66Figure 3
Microstructure of the six aging samples.Moreover, the optic path of the incident lights becomes increasingly more complex owing to the reduced large particles and the increased tiny particles of the aging process. The complexity of the sample’s microstructure enhanced the optic path of the incident lights by increasing its frequency of reflection according to light propagation theory. Coupled with the strong light absorption of PbO2, the amount of emergent light was inevitably reduced and the samples dramatically grew darker over aging time as shown in Figure 1.


### 3.1.3. Spectral Reflectance and Chromaticity

The spectral reflectance of the aging samples was measured with GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in SCE model. As illustrated in Figure4(a), the spectral reflectance of the aging samples of long-wavelength from 550 to 750 nm apparently decreased with the aging time, while the decreased value of the spectral reflectance in the same aging interval becomes smaller and smaller. It is suggested that the darkening trends of samples become slower, but this does not reveal that the chemical composition transition also obeys the same law. The intensity ratio of PbO2 to Pb3O4 in Table 2 shows that the chemical composition transition rates exhibit a thoroughly inverse trend. The inconsistency between sample darkening trends and chemical composition transition rate from red lead to lead dioxide may be mainly caused by the strong color covering power of lead dioxide and its proportion in the aging samples as described above. According to the spectral reflectance change trends, we preliminary infer that the logarithmic function could be well applied to construct the red lead pigment discoloration model under the aging conditions.Figure 4
Distribution of the spectral reflectance of the aging samples (a) and the corresponding first derivate (b).
(a)
(b)In addition, it was easy to find that the spectral reflectance of the aging samples almost has the same rising point in about 550 nm, and the faint decrease in short-wavelength regions is from 360 to 550 nm. As the spectral reflectance of the aging sample is primarily determined by the chemical properties of its components, the spectral features of the aging samples could be explained from the chemical components aspect as the aging samples share the same chemical components except sample 1#. Red lead absorbs the light from 360 to 550 nm and reflects the rest in visible light region, while lead dioxide absorbs almost all visible light. The different ratios of red lead and lead dioxide contained in the aging samples lead to the almost same rising point and the low reflection in short-wavelength region. Furthermore, as shown in Figure4(b), the different ratios of red lead and lead dioxide contained in the aging samples also influence the inflection of the spectral curve. With the amount of lead dioxide increased in the aging samples, the spectral inflection of samples is slightly moved from about 570 to 550 nm.The chromaticity and color difference of the aging samples were summarized in Table4. The brightness value of the aging samples Y and L
∗ as well as the chroma value C
∗ decreased rapidly at the beginning of the experiment, but the speed slows down later in the experiment. Similarly, the color difference DE00 increased rapidly at the beginning of the experiment and the growth slackened later in the experiment. The changing trends of the chromaticity values consisted with the trend of the spectral reflectance in the long-wavelength regions. The reason for such phenomenon was that the chromaticity is linearly calculated from spectral reflectance multiplied by spectral power distributions of the illumination and by color-matching functions [19].Table 4
Chromaticity and color difference of the aging samples.
Sample
Color coordinate
Chroma (C
∗)
D
E
00
Δ
L
∗
Δ
a
∗
Δ
b
∗
Δ
D
E
00
Y
x
y
L
∗
a
∗
b
∗
1#
30.935
0.577
0.371
62.452
52.263
50.530
72.71
0
—
—
—
—
2#
16.951
0.535
0.374
48.198
38.825
36.977
53.62
14.32
14.254
13.438
13.553
14.32
3#
11.300
0.503
0.374
40.081
28.550
25.495
38.28
24.37
8.117
10.275
11.482
10.05
4#
8.241
0.474
0.367
34.479
18.997
16.983
23.59
31.44
5.602
9.533
8.512
7.07
5#
6.454
0.451
0.363
30.530
13.478
12.879
18.64
35.62
3.949
5.499
4.104
4.18
6#
5.920
0.439
0.361
29.209
12.961
12.133
17.75
36.67
1.321
0.517
0.746
1.05The hue of the aging samples ina
∗
-
b
∗ planes and in CIE 1931 chromaticity diagram was plotted in Figure 5, where a
∗ represents the red-green response and b
∗ represents the yellow-blue response from human visual system in CIE L
∗
a
∗
b
∗ uniform color space, respectively. Figure 5(a) indicated that the hue changed trends can be fitted very well with the linear model, and Figure 5(b) shows the same varying tendencies. The values of Δ
a
∗, Δ
b
∗, and Δ
L
∗ were also summarized in Table 4, which shows the same changing trends as ΔDE00. Based on the data changing trends of Δ
a
∗, Δ
b
∗, Δ
L
∗, and C
∗ in Table 4 and the hue distributions of aging samples plotted in Figure 5, it could be concluded that, during the red lead pigment aging process, the main changing chromaticity feature was the brightness and chroma value, and the hue of the aging samples almost remains unchanged.Figure 5
Hues of aging samples ina
∗
-
b
∗ planes (a) and CIE 1931 chromaticity diagram (b).
(a)
(b)In the next part, the simulating degradation experiment was carried out with referring to the result of aging test. And the spectral reflectance and chromaticity of the aging samples were further analyzed for construction of the spectral discoloration model.


## 3.1.1. Chemical Composition

As illustrated in Figure2, XRD results indicated that the aging product of red lead under the aging test environment was lead dioxide, explained by the increased diffraction peaks of PbO2 and the decreased diffraction peaks of Pb3O4. This is consistent with the previous researches [2, 4, 5]. Furthermore, the diffraction intensity ratio of PbO2 to Pb3O4 in Table 2 and the discrimination of color between the adjacent samples in Figure 1 indicated the strong color covering power of lead dioxide. A little amount of lead dioxide may darken the sample remarkably, but when the amount of lead dioxide reaches a certain degree, its color covering ability will reach the maximum.Table 2
The diffraction intensity ratio of PbO2 to Pb3O4 of the aging samples.
Sample
1#
2#
3#
4#
5#
6#
Intensity ratio (I
P
b
O
2/I
P
b
3
O
4)
0
0.048
0.072
0.138
0.258
1.503Figure 2
XRD results for the aging samples:(
) represents the diffraction peak of PbO2. ∗ represents the diffraction peak of Pb3O4.


## 3.1.2. Microstructure

The microstructure image of the aging samples captured with SEM is shown in Figure3. It demonstrated that the amount of the large particle is decreasing and the amount of tiny particle is increasing faced with the aging process, but there are still some large particles wrapped by tiny particles. The grain size distribution of the aging samples measured through analysis software was summarized in Table 3. The statistics shows that the maximum, minimum, and mean particle size are reduced as the aging time passes, and the minimum particle size is ultimately stable at the level of 0.21 μm. Based on the particle size changing trend described above, the minimum particles in the aging samples may be the lead dioxide. The phenomenon of the large particles wrapped by tiny particles may be the reason why the previous researchers observed the incomplete oxidized red lead pigment grains in the cross section of the pigment layer painted by red lead pigment [7, 8].Table 3
Particles distribution of the aging samples.
Sample
1#
2#
3#
4#
5#
6#
Max/(μm)
8.21
5.41
3.62
2.52
2.41
2.12
Min/(μm)
0.63
0.5
0.4
0.38
0.21
0.21
Mean/(μm)
2.12
1.78
1.36
1.03
0.79
0.66Figure 3
Microstructure of the six aging samples.Moreover, the optic path of the incident lights becomes increasingly more complex owing to the reduced large particles and the increased tiny particles of the aging process. The complexity of the sample’s microstructure enhanced the optic path of the incident lights by increasing its frequency of reflection according to light propagation theory. Coupled with the strong light absorption of PbO2, the amount of emergent light was inevitably reduced and the samples dramatically grew darker over aging time as shown in Figure 1.


## 3.1.3. Spectral Reflectance and Chromaticity

The spectral reflectance of the aging samples was measured with GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in SCE model. As illustrated in Figure4(a), the spectral reflectance of the aging samples of long-wavelength from 550 to 750 nm apparently decreased with the aging time, while the decreased value of the spectral reflectance in the same aging interval becomes smaller and smaller. It is suggested that the darkening trends of samples become slower, but this does not reveal that the chemical composition transition also obeys the same law. The intensity ratio of PbO2 to Pb3O4 in Table 2 shows that the chemical composition transition rates exhibit a thoroughly inverse trend. The inconsistency between sample darkening trends and chemical composition transition rate from red lead to lead dioxide may be mainly caused by the strong color covering power of lead dioxide and its proportion in the aging samples as described above. According to the spectral reflectance change trends, we preliminary infer that the logarithmic function could be well applied to construct the red lead pigment discoloration model under the aging conditions.Figure 4
Distribution of the spectral reflectance of the aging samples (a) and the corresponding first derivate (b).
(a)
(b)In addition, it was easy to find that the spectral reflectance of the aging samples almost has the same rising point in about 550 nm, and the faint decrease in short-wavelength regions is from 360 to 550 nm. As the spectral reflectance of the aging sample is primarily determined by the chemical properties of its components, the spectral features of the aging samples could be explained from the chemical components aspect as the aging samples share the same chemical components except sample 1#. Red lead absorbs the light from 360 to 550 nm and reflects the rest in visible light region, while lead dioxide absorbs almost all visible light. The different ratios of red lead and lead dioxide contained in the aging samples lead to the almost same rising point and the low reflection in short-wavelength region. Furthermore, as shown in Figure4(b), the different ratios of red lead and lead dioxide contained in the aging samples also influence the inflection of the spectral curve. With the amount of lead dioxide increased in the aging samples, the spectral inflection of samples is slightly moved from about 570 to 550 nm.The chromaticity and color difference of the aging samples were summarized in Table4. The brightness value of the aging samples Y and L
∗ as well as the chroma value C
∗ decreased rapidly at the beginning of the experiment, but the speed slows down later in the experiment. Similarly, the color difference DE00 increased rapidly at the beginning of the experiment and the growth slackened later in the experiment. The changing trends of the chromaticity values consisted with the trend of the spectral reflectance in the long-wavelength regions. The reason for such phenomenon was that the chromaticity is linearly calculated from spectral reflectance multiplied by spectral power distributions of the illumination and by color-matching functions [19].Table 4
Chromaticity and color difference of the aging samples.
Sample
Color coordinate
Chroma (C
∗)
D
E
00
Δ
L
∗
Δ
a
∗
Δ
b
∗
Δ
D
E
00
Y
x
y
L
∗
a
∗
b
∗
1#
30.935
0.577
0.371
62.452
52.263
50.530
72.71
0
—
—
—
—
2#
16.951
0.535
0.374
48.198
38.825
36.977
53.62
14.32
14.254
13.438
13.553
14.32
3#
11.300
0.503
0.374
40.081
28.550
25.495
38.28
24.37
8.117
10.275
11.482
10.05
4#
8.241
0.474
0.367
34.479
18.997
16.983
23.59
31.44
5.602
9.533
8.512
7.07
5#
6.454
0.451
0.363
30.530
13.478
12.879
18.64
35.62
3.949
5.499
4.104
4.18
6#
5.920
0.439
0.361
29.209
12.961
12.133
17.75
36.67
1.321
0.517
0.746
1.05The hue of the aging samples ina
∗
-
b
∗ planes and in CIE 1931 chromaticity diagram was plotted in Figure 5, where a
∗ represents the red-green response and b
∗ represents the yellow-blue response from human visual system in CIE L
∗
a
∗
b
∗ uniform color space, respectively. Figure 5(a) indicated that the hue changed trends can be fitted very well with the linear model, and Figure 5(b) shows the same varying tendencies. The values of Δ
a
∗, Δ
b
∗, and Δ
L
∗ were also summarized in Table 4, which shows the same changing trends as ΔDE00. Based on the data changing trends of Δ
a
∗, Δ
b
∗, Δ
L
∗, and C
∗ in Table 4 and the hue distributions of aging samples plotted in Figure 5, it could be concluded that, during the red lead pigment aging process, the main changing chromaticity feature was the brightness and chroma value, and the hue of the aging samples almost remains unchanged.Figure 5
Hues of aging samples ina
∗
-
b
∗ planes (a) and CIE 1931 chromaticity diagram (b).
(a)
(b)In the next part, the simulating degradation experiment was carried out with referring to the result of aging test. And the spectral reflectance and chromaticity of the aging samples were further analyzed for construction of the spectral discoloration model.


## 3.2. Simulating Degradation Experiment

In this experiment, eleven experimental samples were painted using the painting method simulating Dunhuang murals. The samples were mixed by red lead and lead dioxide in different mixing ratios with total mass controlled at 0.6 g as the quantity is enough for painting a sample in this study, and the mass change step of the pure red lead and lead dioxide was set at 0.06 g. The name and the corresponding mixing ratio of the simulating degradation samples are listed in Table5.Table 5
Names and mixing ratios of the simulating degradation samples.
Sample
S1
S2
S3
S4
S5
S6
S7
S8
S9
S10
S11
Pb3O4/g
0.60
0.54
0.48
0.42
0.36
0.30
0.24
0.18
0.12
0.06
0.00
PbO2/g
0.00
0.06
0.12
0.18
0.24
0.30
0.36
0.42
0.48
0.54
0.60The spectral reflectance of the simulating degradation samples was also measured by GretagMacbeth X-Rite Color-Eye 7000A spectrophotometer in SCE model, and the chromaticity of each simulating degradation sample was calculated using Matlab toolbox OptProp. The color of each simulating degradation sample is rendered in Figure6. The spectral reflectance and the corresponding first derivate of the simulating degradation samples are plotted in Figure 7, and the chromaticity value of the simulating degradation samples is summarized in Table 6.Table 6
Chromaticity features of the simulating degradation samples.
Sample
Color coordinate
Chroma (C
∗)
D
E
00
Δ
L
∗
Δ
a
∗
Δ
b
∗
Δ
D
E
00
Y
x
y
L
∗
a
∗
b
∗
S1
30.935
0.577
0.371
62.452
52.263
50.530
72.71
0
—
—
—
—
S2
16.521
0.458
0.366
47.651
26.696
25.834
37.14
17.35
14.801
25.567
24.696
17.35
S3
13.228
0.426
0.361
43.104
19.287
18.466
26.69
23.68
4.547
7.409
7.368
6.33
S4
10.788
0.399
0.354
39.227
14.140
12.822
19.08
29.05
3.877
5.147
5.644
5.37
S5
10.139
0.385
0.351
38.094
11.601
10.349
15.54
31.00
1.133
2.539
2.473
1.95
S6
8.425
0.377
0.349
34.858
9.590
8.676
12.92
34.22
3.236
2.011
1.673
3.22
S7
7.402
0.361
0.345
32.712
6.923
6.237
9.31
36.78
2.146
2.667
2.439
2.56
S8
6.625
0.350
0.343
30.944
4.977
4.770
6.89
38.66
1.768
1.946
1.467
1.88
S9
6.272
0.336
0.338
30.099
3.108
2.701
4.12
40.14
0.845
1.869
2.069
1.48
S10
5.650
0.327
0.336
28.512
1.782
1.633
2.42
41.65
1.587
1.326
1.068
1.51
S11
5.225
0.316
0.332
27.375
0.276
0.345
0.44
43.03
1.137
1.506
1.288
1.38Figure 6
Simulating degradation samples.Figure 7
Distribution of the spectral reflectance of the simulating degradation samples (a) and the corresponding first derivate (b).
(a)
(b)It indicated that there is a strong similarity between the aging test and simulating degradation experiment about the changing trend of the spectral reflectance, the first derivate of spectral reflectance, and the chromaticity values. However, there is slight difference between the spectral reflectance curves shapes of the two experiments according to the results of gelatin degradation in Dunhuang murals carried out by Li [21]. The humidity and light are two important factors causing the aging of the gelatin which lead to the gelatin solution from almost transparent status gradually becoming light yellow, so it can be inferred that the slight difference between the spectral reflectance curves shapes may be caused by the aging of the gelatin in the samples. However, the similarity of spectral and chromaticity features of the samples between the aging test and simulating degradation experiment inferred that the aging of the gelatin just has a very small influence on the aging trends of red lead pigment. The hue distribution of the simulating degradation samples in a
∗
-
b
∗ planes and CIE 1931 chromaticity diagram was shown in Figure 8. The hue distribution also maintains high similarity with the aging test result as it can be fitted very well using the linear model. The similarity of the two experimental results may also prove that the mechanism of the red lead pigment discoloration in Dunhuang murals is caused by the process of red lead altered in lead dioxide.Figure 8
Hue distribution of the simulating degradation samples ina
∗
-
b
∗ planes (a) and CIE 1931 chromaticity diagram (b).
(a)
(b)


## 3.3. Spectral Discoloration Model Construction

The decreased trend of spectral reflectance in long-wavelength regions of the simulating degradation samples was fitted with some simple linear or nonlinear functions. Here the mass ratio of lead dioxide in total mass was used to monitor the change of spectral reflectance. Figure9 shows the decreasing trend of the spectral reflectance at 700 nm. It can be seen from Figure 9 that the decreased trend could be fitted very well using the logarithmic model and that R
2 value of the fitting equation is high up to 0.99515. Moreover, the average R
2 value from 560 to 750 nm is higher than 0.97, and the minimum R
2 value is 0.95329.Figure 9
Decreasing trend of the spectral reflectance of the simulating degradation samples at 700 nm.The other simple functions were also adopted to fit the changing trend of the simulating degradation samples.R
2 values of these fitting functions and the logarithmic fitting function are summarized in Table 7. R
2 values in Table 7 show that the logarithmic function has the highest accuracy which is consistent with the previous hypothesis. Furthermore, the decreasing trend in short-wavelength regions from 360 to 550 nm acquires the best fitting accuracy using the linear function. The fitting result shows that the average R
2 value from 360 to 550 nm is higher than 0.97 and the minimum R
2 value is 0.95993. Thus the spectral discoloration model could be constructed with the abovementioned results. Besides, the spectral reflectance is regarded as the “fingerprint” of the object color and all the chromaticity values are calculated from the spectral reflectance under different illuminating and viewing conditions, so there is no need to construct the discoloration model in color space.Table 7
R
2 values of different fitting functions at 700 nm.
Functions
R
2
Exponential
0.92245
Linear
0.63268
Logarithmic
0.99515
Polynomial(
2
)
∗
0.84378
Power
0.82483
∗Note: “polynomial (
2
)” represents the two-order polynomial.Based on the above analysis, the spectral discoloration model of red lead pigment was constructed as follows:(2)
Y
i
=
A
x
+
B
,
if
1
≤
i
<
22
,
A
ln
⁡
x
+
B
,
if
22
≤
i
≤
40
,where i represents the number of the sampling points from 360 to 750 nm at 10 nm intervals; x is mass ratios of lead dioxide in the simulating degradation samples; Y
i denotes the spectral reflectance in different sampling point for the specific x; and A and B contain the model parameters in each sampling wavelength from 360 to 750 nm at 10 nm interval; namely, (3)
A
=
A
1
⋮
A
i
⋮
A
40
,
B
=
B
1
⋮
B
i
⋮
B
40
.After the spectral discoloration model was constructed, its prediction accuracy was firstly examined by the simulating degradation samples with the known aging degree ofx. The results show that the constructed spectral discoloration model is effective as the average spectral root-mean-square error, the average spectral goodness-of-fit coefficient [22], and the average color difference DE00 of the eleven simulation aging samples are all kept in a fairly good level at 0.0591, 0.9979, and 1.14, respectively.


## 3.4. Real Aging Red Lead Pigment in Dunhuang Murals

The red lead pigment is widely used in Dunhuang murals, especially for the flesh of the Buddha in murals. Most of the murals painted with red lead pigment are black today. But there is still some incompletely discolored red lead pigment identified by chemical analysis of previous researchers [8, 23]. Several red lead pigment points in different aging degree in two typical grottoes from different dynasties were measured with the in situ nondestructive measurement method. The validity of the model was verified by calculating and predicting the spectral reflectance of the real aging red lead pigment.The PR705 spectroradiometer with the standard whiteboard produced by American Photo Research was used to measure the spectral reflectance of the real aging red lead pigment in Dunhuang murals, and the actual measured work scene is shown in Figure10. The insert in Figure 10 represents the geometrical conditions of illumination and measurement. The points measured by PR705 spectroradiometer were marked in Figure 11. The measured spectral reflectance of the five marked points was plotted in Figure 12.Figure 10
Spectral reflectance measured using PR705 spectroradiometer in cave 465 (Yuan Dynasty).Figure 11
Five real aging red lead pigment points in Dunhuang murals: (a) Yuan Dynasty cave 465 and (b) Early Tang Dynasty cave 220.
(a)
(b)Figure 12
Distribution of spectral reflectance of the marked five points in the figure.Because of the larger measurement noise in the short-wavelength below 400 nm and the long-wavelength above 700 nm, the effective wavelength range of the sampled points was restricted to 400~700 nm as plotted in Figure12. It can be seen that there is a very strong similarity between the spectral reflectance shapes of the real aging red lead pigment and the above two experiment samples. The hue distribution of the real aging red lead points also can be fitted well with the linear model with the slope, the intercept, and R
2 value being equal to 0.9692, 0.5621, and 0.9597. Based on the constructed spectral discoloration model, the simulated aging degree of the five real aging red lead pigment points was calculated as 0.8028%, 31.4346%, 5.2151%, 9.7909%, and 73.4961% using the inversed spectral discoloration model.After the detection of aging degree, the spectral reflectance of the five points was predicted using the proposed spectral discoloration model as plotted in Figure13 with the dash-dot line. The prediction accuracy was also evaluated by the metrics of spectral root-mean-square error, the spectral goodness-of-fit coefficient, and the color difference DE00. The predicted results of the five points are summarized in Table 8. The average spectral root-mean-square error and color difference DE00 are 1.21% and 2.64, respectively. The average spectral goodness-of-fit coefficient is high up to 99.75%. The predicted results indicated that the constructed spectral discoloration model could well predict the spectral reflectance of the real aging red lead pigment in Dunhuang murals based on the simulated aging degree. It can be concluded that the constructed spectral discoloration model for red lead pigment could work well in detecting the aging degree of red lead pigment in ancient murals and reproducing the spectral data for specific aging degree.Table 8
The predicted accuracy of the five real aging red lead pigment points.
Point number
RMS (%)
DE00
GFC (%)
1
1.50
1.99
99.93
2
1.30
3.24
99.59
3
1.28
2.61
99.88
4
1.14
2.87
99.86
5
0.81
2.48
99.47
Mean
1.21
2.64
99.75Figure 13
Measured (solid line) and predicted (dash-dot line) spectral reflectance of the five real aging red lead pigment points using the spectral discoloration model.


## 4. Conclusion

A spectral discoloration model for the aging red lead pigment in Dunhuang murals was tentatively constructed to detect the aging degree of red lead pigment and to reproduce the spectral reflectance based on aging test and simulating degradation experiment. The chemistry and physical attributes of red lead pigment during its aging process under the aging conditions were characterized, where the red lead was identified changing into lead dioxide and the amount of tiny particles increasing in the painted layer faced with the aging process. The spectral reflectance and chromaticity of the aging red lead pigment were analyzed and the spectral discoloration model was constructed to describe the red lead pigment aging process. The test results indicate that the constructed spectral discoloration model could be well applied to calculate the aging degree and to predict the spectral reflectance of red lead pigment in Dunhuang murals. Based on the constructed spectral discoloration model, the virtual restoration of the color of the darkened murals painted by red lead pigments will be more accurate and scientific.



---

*Source: 3517496-2016-09-01.xml*