A Semi-Theoretical Model for Water Condensation: Dew Used in Conservation of Earthen Heritage Sites

: Dew is a common but important phenomenon. Though water is previously considered to be a threat to earthen heritage sites, artiﬁcial dew is showing potential in relic preservation. A model of dew prediction on earthen sites will be essential for developing preventive protection methods, but studies of dew formation processes on relics are limited. In this study, a two parameter model is proposed. It makes approximations according to the features of earthen heritage sites, assuming that a thin and steady air layer exists close to the air–solid interface. This semi-theoretical model was based on calculations of the mass transfer process in the air layer, and was validated by simulations of laboratory experiments (R > 0.9) as well as ﬁeld experiments. Additionally, a numerical simulation, performed by the commercial software COMSOL, conﬁrmed that the difference between ﬁtting parameter δ and the thickness of assumed mass transfer ﬁeld was not signiﬁcant. This model will be helpful in developing automatic environmental control systems for stabilizing water and soluble salts, thus enhancing preventive protection of earthen heritage sites.


Introduction
Dew is a common phenomenon and is becoming increasingly important nowadays [1][2][3][4]. It forms out of water in the air once the temperature of a solid surface in contact with the air drops below the dew point.
In the conservation of cultural relics, dew has long been considered as a potential cause of damage. For example, in porous structures such as bricks, murals, sculptures and earthen sites, repeated condensation and evaporation of water can activate soluble salts. With water transportation, salts may migrate, crystallize and swell in pores, causing stress and resulting in destruction [5,6]. In addition, dew may cause corrosion of metals and worsen biological corrosion [7].
Extensive dew and water evaporation are especially serious problems for earthen heritage sites. Various kinds of materials have been developed to consolidate soil surface, such as potassium-silica solution, iron-based bioproducts, calcium hydroxide nanoparticles and organosilicones [8][9][10][11][12]. However, there is no successful case reported in humid regions after long-term preservation. The reason is clear: no matter how high the strength would be, it is just a matter of time before the periodic migration and crystallization of soluble salts destroy the structures from the inside [13]. Inappropriate materials may even cause preservation damage [14]. In the Kuahuqiao Museum, Hangzhou, China, the exhibition area is below the water table. Soluble salts have caused damage on both consolidated and unconsolidated soil surfaces (Figure 1), with diseases rates of 30.54% and 13.65%, respectively [15]. Artificial dew is a more promising method for preserving earthen heritage sites [16], though natural dew have caused problems. In recent years, preventive protection methods have received more and more recognition in the conservation of relics [10,17], and the key is minimizing potential risks by controlling and stabilizing surrounding environments. However, the example of the Kuahuqiao Museum shows that air conditioning without stabilizing inner water cannot stop salt damage. [15] As an improvement, it is promising to produce artificial dew on soil surface, using an automatic monitoring and controlling system of environment. By controlling dew amount and position, equilibrium of water potential can be built to stabilize water inside soil and prevent migration of salts. This preventive protection method will be gentle and can provide a fundamental solution for salt destruction.
In a previous study, artificial dew was produced on a soil surface. [16] Experiments in lab showed that dew formation was affected by temperature and humidity, while in situ experiments showed that dew could take salts downward. However, a rational model of dew formation rate was absent, which is essential for developing an environmental control system. Dew prediction is a common problem but has not been well studied in terms of preserving relics. Up to now, many models have been developed to predict the duration and amount of water condensation on soil [18][19][20][21][22] and have been well summarized [23]. Most of them were developed using energy and mass transport equations, or were based on assessments of latent-heat flux in the atmosphere [24]. However, none of them were developed for relic conservation uses. Alterations and simplifications can be made to produce a new model according to the features of earthen heritage sites.
In this study, we report the construction of a semi-theoretical model for preservation of earthen sites. The model is based on the mass transfer process in an assumed steady air layer, and the main driven force is the differential pressure of water vapor. The model contains three measurable variables and two fitting parameters. Simulations of dew fomation in both laboratory and field experiments from published studies are used to confirm its validity and application potential.

Theory
This model assumes that there is a steady air layer close to the air-solid interface ( Figure 2). It is thin and not disturbed by turbulent flows in the air. Water vapor would pass this layer before condensation on a solid surface. The mass flow of water vapor in this thin layer is driven mainly by the differential pressure of water vapor, while slightly Artificial dew is a more promising method for preserving earthen heritage sites [16], though natural dew have caused problems. In recent years, preventive protection methods have received more and more recognition in the conservation of relics [10,17], and the key is minimizing potential risks by controlling and stabilizing surrounding environments. However, the example of the Kuahuqiao Museum shows that air conditioning without stabilizing inner water cannot stop salt damage. [15] As an improvement, it is promising to produce artificial dew on soil surface, using an automatic monitoring and controlling system of environment. By controlling dew amount and position, equilibrium of water potential can be built to stabilize water inside soil and prevent migration of salts. This preventive protection method will be gentle and can provide a fundamental solution for salt destruction.
In a previous study, artificial dew was produced on a soil surface. [16] Experiments in lab showed that dew formation was affected by temperature and humidity, while in situ experiments showed that dew could take salts downward. However, a rational model of dew formation rate was absent, which is essential for developing an environmental control system. Dew prediction is a common problem but has not been well studied in terms of preserving relics. Up to now, many models have been developed to predict the duration and amount of water condensation on soil [18][19][20][21][22] and have been well summarized [23]. Most of them were developed using energy and mass transport equations, or were based on assessments of latent-heat flux in the atmosphere [24]. However, none of them were developed for relic conservation uses. Alterations and simplifications can be made to produce a new model according to the features of earthen heritage sites.
In this study, we report the construction of a semi-theoretical model for preservation of earthen sites. The model is based on the mass transfer process in an assumed steady air layer, and the main driven force is the differential pressure of water vapor. The model contains three measurable variables and two fitting parameters. Simulations of dew fomation in both laboratory and field experiments from published studies are used to confirm its validity and application potential.

Theory
This model assumes that there is a steady air layer close to the air-solid interface ( Figure 2). It is thin and not disturbed by turbulent flows in the air. Water vapor would pass this layer before condensation on a solid surface. The mass flow of water vapor in this thin layer is driven mainly by the differential pressure of water vapor, while slightly affected by other factors such as temperature, wind, salt content, soil suction, height and vegetation [25,26]. affected by other factors such as temperature, wind, salt content, soil suction, height and vegetation [25,26]. For earthen heritages under preservation, their environments are supposed to be stable. For example, wind should be minimized, vegetation should be removed, desalting should be conducted, heat of sunlight should be controlled and variations of soil properties at a site are often limited. As a result, many factors in other models can be neglected in this model, and the assumption of a steady air layer is reasonable.
The model solves for the mass flow in the air layer. Firstly, differential pressure of water vapor should be the main driving force of water condensation. This will cause molecular diffusion. The process can be described using Fick's laws of diffusion and their derivatives [27]: (1) Here, Ndiff (mol m −2 s −1 ) is the diffusion rate of water vapor, D (m 2 s −1 ) is the molecular diffusivity of water vapor, R = 8.314 J mol −1 K −1 is the ideal gas constant, Ts (K) is the absolute temperature of interface, pv1 (Pa) is the partial pressure of water vapor in the atmosphere, pv2 (Pa) is the partial pressure of water vapor upon the soil surface and δ (m) is the thickness of the mass transfer field. Then, alterations to Equation (1) must be made to represent the influence of temperature, wind, etc. Therefore, parameter a is added: Here, N is the total rate of dew formation. D, pv1, and pv2 can be easily calculated using the three variables: soil surface temperature Ts, air temperature Ta and air relative humidity RH. All variables are easy to measure.
The diffusivity of water vapor in atmosphere near the soil surface can be computed by [28]: For earthen heritages under preservation, their environments are supposed to be stable. For example, wind should be minimized, vegetation should be removed, desalting should be conducted, heat of sunlight should be controlled and variations of soil properties at a site are often limited. As a result, many factors in other models can be neglected in this model, and the assumption of a steady air layer is reasonable.
The model solves for the mass flow in the air layer. Firstly, differential pressure of water vapor should be the main driving force of water condensation. This will cause molecular diffusion. The process can be described using Fick's laws of diffusion and their derivatives [27]: Here, N diff (mol m −2 s −1 ) is the diffusion rate of water vapor, D (m 2 s −1 ) is the molecular diffusivity of water vapor, R = 8.314 J mol −1 K −1 is the ideal gas constant, T s (K) is the absolute temperature of interface, pv 1 (Pa) is the partial pressure of water vapor in the atmosphere, pv 2 (Pa) is the partial pressure of water vapor upon the soil surface and δ (m) is the thickness of the mass transfer field. Then, alterations to Equation (1) must be made to represent the influence of temperature, wind, etc. Therefore, parameter a is added: Here, N is the total rate of dew formation. D, pv 1 , and pv 2 can be easily calculated using the three variables: soil surface temperature T s , air temperature T a and air relative humidity RH. All variables are easy to measure. The diffusivity of water vapor in atmosphere near the soil surface can be computed by [28]: where D 0 = 2.29 × 10 −5 m 2 s −1 .
Here, p s,T is the saturated vapor pressure of water at temperature T. Over the temperature range 273~323 K, A = 58.431, B = −6750.4 and C = −4.8668 [29]. According to the definition of relative humidity, the vapor pressure in the atmosphere (pv 1 ) is given by: where RH is the relative humidity in atmosphere and T a is the air temperature. At the air-soil interface, water vapor condenses and pv 2 can be calculated considering phase equilibrium. For desalinized soil, the vapor pressure on the soil surface equals the saturated vapor pressure: Here, T s is the surface temperature of soil. If soluble salts were added to the soil, the condensed water would form a solution, and pv 2 would be reduced according to Raoult's law. In the phase equilibrium of saturated Na 2 SO 4 solution and water vapor, relative humidity of air is a function of temperature. For the temperature range 0~20 • C, the following linear equation was used [30]: Then, pv 2 can be given using the combination of Equations (4) and (7). If the soil property is unknown, Equation (6) can be used. However, this may cause deviations because of different soluble salts.

Data Collection
Data including environmental variables and dew amounts were collected from published articles, and simulations were conducted to validate the proposed model. One piece of literature is our previous study [16] and the other is an important study conducted by Jacobs et al. [18,19,31]. Experimental protocols of our previous study are given in the following sections, because they were originally in Chinese.

Laboratory Experiments
Water 2021, 13, 52 5 of 12 overfitting and large deviations beyond tested conditions. In this study, experimental data were reprocessed from the beginning using the proposed semi-theoretical model.

Field Experiments in the Kuahuqiao Museum
In the previous study, field experiments were conducted to preliminarily prove that artificial dew can take surface salts downward. [16] Monitored environmental variables and cumulative dew amounts were simulated using a new model in this study.
At the Kuahuqiao Site Museum, an 8000-year old dugout canoe is on display. As preservation methods of the boat and the surrounding earthen site, the museum firstly built a ceiling to prevent rain wash, and then installed an air-conditioning system to prevent microorganism growth and excessive water condensation. Consolidation of soil surface using epoxy resin was also carried out [32]. However, these methods prevented water transportation and salt migration uniaxially. After the construction of ceiling and environmental control system, soluble salts moved upward along with water evaporation, crystallized on soil surfaces, and caused damage, [15] as shown in Figure 1. This suggests that factors such as water transportation must be considered in controlling surrounding environments, as a part of preventive protection. The experimental setup used is shown in Figure 4. Five 0.5 m × 0.5 m square boards pasted with electrothermal films were spliced and arranged to form a cubic space. A humidifier composed of a water bath, a plastic cover and a fan was placed in the cubic space. The walls were then heated and moisture was produced for 6 h. The temperature of the electrothermal films and the water bath were kept at 45 and 50 • C, respectively. The electrical power was then switched off and the top board removed for 1 h so that the soil could cool off naturally. The protocol was carried out four times. Water content inside the soil (0~5 cm deep) was measured using a soil temperature and humidity measuring instrument.
prevent microorganism growth and excessive water condensation. Consolidation of soil surface using epoxy resin was also carried out [32]. However, these methods prevented water transportation and salt migration uniaxially. After the construction of ceiling and environmental control system, soluble salts moved upward along with water evaporation, crystallized on soil surfaces, and caused damage, [15] as shown in Figure 1. This suggests that factors such as water transportation must be considered in controlling surrounding environments, as a part of preventive protection.
The experimental setup used is shown in Figure 4. Five 0.5 m × 0.5 m square boards pasted with electrothermal films were spliced and arranged to form a cubic space. A humidifier composed of a water bath, a plastic cover and a fan was placed in the cubic space. The walls were then heated and moisture was produced for 6 h. The temperature of the electrothermal films and the water bath were kept at 45 and 50 °C, respectively. The electrical power was then switched off and the top board removed for 1 h so that the soil could cool off naturally. The protocol was carried out four times. Water content inside the soil (0~5 cm deep) was measured using a soil temperature and humidity measuring instrument.  [16]. Water will condense on the soil surface where the temperature was below dew point.

Field Experiments in a Desert System
Jacobs et al. [18,19,31] conducted an important study on a desert system. They measured daily dew amounts as well as micrometeorological conditions and constructed a simulation model. Environmental variables and cumulative dew amounts were collected from published figures in their study.  [16]. Water will condense on the soil surface where the temperature was below dew point.

Field Experiments in a Desert System
Jacobs et al. [18,19,31] conducted an important study on a desert system. They measured daily dew amounts as well as micrometeorological conditions and constructed a simulation model. Environmental variables and cumulative dew amounts were collected from published figures in their study.

Simulations Using the Semi-Theoretical Model
Data including environmental variables and cumulative dew amounts were reprocessed using Equations (2)-(7).
First, cumulative dew curves of laboratory experiments were manually divided into two phases, because transition points were observed at similar position of curves. The relationship between cumulative dew and time was estimated using linear least square fitting with Excel software, and the slopes of the fitting lines are the water condensation rates.
Then, simulations of the dew formation rates were conducted using the proposed model. The least squares curve fit routine (lsqcurvefit) with optimization was used to estimate the values of the two parameters. Curve fitting was performed using Matlab R2016b software (MathWorks, Inc., Natick, MA, America) and the codes are shown in the supplementary materials file "phase1.m" and "phase2.m".
Last, simulations of field experiments were conducted. Dew formation rates at specific points in time were estimated. Corresponding environmental variables were used to estimate the fitting parameters of the field experiment. Equation (6) was used in fitting because salts in soil were not controlled. Combining environmental variable curves with fitting parameters, the dew formation rate curve and cumulative dew amount curve were constructed.

Numerical Simulations of Laboratory Experiments
Numerical simulations were performed to validate the assumption of a thin and steady air layer in the proposed model. In laboratory experiments, soil samples were placed in culture dishes. As a result, an equivalent air layer was constructed in numerical simulations. The simulation is axisymmetric two-dimensional including a fluid domain composed of moist air ( Figure 5). Equivalent layer thickness δ e was calculated and was compared with the fitting parameter δ to check consistency.
Last, simulations of field experiments were conducted. Dew formation rates at specific points in time were estimated. Corresponding environmental variables were used to estimate the fitting parameters of the field experiment. Equation (6) was used in fitting because salts in soil were not controlled. Combining environmental variable curves with fitting parameters, the dew formation rate curve and cumulative dew amount curve were constructed.

Numerical Simulations of Laboratory Experiments
Numerical simulations were performed to validate the assumption of a thin and steady air layer in the proposed model. In laboratory experiments, soil samples were placed in culture dishes. As a result, an equivalent air layer was constructed in numerical simulations. The simulation is axisymmetric two-dimensional including a fluid domain composed of moist air ( Figure 5). Equivalent layer thickness δe was calculated and was compared with the fitting parameter δ to check consistency.  Simulations were performed with Comsol Multiphysics version 5.4 (COMSOL, Inc., Stockholm, Sweden), which was installed on a Dell OptiPlex 3060 with Intel Core i5-8500 (3 GHz). Three interfaces including laminar flow, heat transfer in moist air and moisture transport in air were used. Natural convection and moist transportation were solved together as a stationary problem, and optimization studies were conducted to determine the thickness of the mass transfer field.
Input boundary conditions including three environmental variables are shown in Table 1. Relative humidity on the upper surface (RH), temperature on the upper surface (T a ) and temperature on the lower surface (T s ) are listed. Additionally, the relative humidity on the lower surface (RH interface ) was calculated using Equation (7).
In the optimization study, the height of the computational domain (thickness of mass transfer field δ e ) was the controlled variable and optimization objective was minimizing the difference between the water vapor flux and the measured dew formation rate (N). Mesh and other unspecified settings remained as default. (Figure 5c)

Laboratory Experiments and Numerical Simulations
Cumulative dew curves of laboratory experiments were composed of two phases, both of which were linearly fitted ( Figure 6 and supplementary materials file "Figs lab cum dew.pdf", all R-square values are above 0.99). The transition points between the two phases were about 0.8-1 kg/m 2 . Around these points, water could be observed on the soil surface and the soil samples became water-saturated. Then, simulation of the dew formation rate was conducted using Equations (2)- (7) and values of the fitting parameters were determined.
Mesh and other unspecified settings remained as default. (Figure 5c)

Laboratory Experiments and Numerical Simulations
Cumulative dew curves of laboratory experiments were composed of two phases, both of which were linearly fitted ( Figure 6 and supplementary materials file "Figs lab cum dew.pdf", all R-square values are above 0.99). The transition points between the two phases were about 0.8-1 kg/m 2 . Around these points, water could be observed on the soil surface and the soil samples became water-saturated. Then, simulation of the dew formation rate was conducted using Equations (2)- (7) and values of the fitting parameters were determined. Results are given in Table 1. Simulations of both phases produced satisfying results. Statistical tests revealed that the equations for both phases I and II explained a high proportion of the variation in the dew formation rate-i.e., both R 2 values between measured and simulated dew formation rates were >0.81. The unexplained variance can be attributed to systematic error caused by accuracy of the model as well as random error in observation.
In the proposed model, a steady air layer was assumed and the inner water vapor flux equals the measured dew formation rate N. As a result, an equivalent air layer was constructed in a numerical simulation. Once N and the other boundary conditions (RH, Ta, Ts) were set, the equivalent thickness of the assumed air layer (denoted by δe) was determined. Results are also shown in Table 1. Results are given in Table 1. Simulations of both phases produced satisfying results. Statistical tests revealed that the equations for both phases I and II explained a high proportion of the variation in the dew formation rate-i.e., both R 2 values between measured and simulated dew formation rates were >0.81. The unexplained variance can be attributed to systematic error caused by accuracy of the model as well as random error in observation.
In the proposed model, a steady air layer was assumed and the inner water vapor flux equals the measured dew formation rate N. As a result, an equivalent air layer was constructed in a numerical simulation. Once N and the other boundary conditions (RH, T a , T s ) were set, the equivalent thickness of the assumed air layer (denoted by δ e ) was determined. Results are also shown in Table 1.

Field Experiments at the Kuahuqiao Site Museum
During the field experiment, humidity, air temperatures and soil temperature were monitored (Figure 7a). The measured cumulative amount of dew (Figure 7b) was estimated using known dry soil density (1.802 g/cm 3 ) and the measured moisture content. These original data were cited, and a cumulative dew curve was constructed using Equations (2)- (6) to simulate this field experiment. Simulation was satisfied when δ = 4.58 × 10 −3 and a = 2.91 × 10 3 .
Three small concaves can be observed in the simulation curve (Figure 7b) in accordance with the fluctuations of environmental variables, but disappeared in the measured curve. This is explained by the delay in measuring dew amounts. Measured dew amounts were calculated using the water content beneath the soil surface while simulated dew amounts were on air-soil surface. When the surface dew formation stopped, free water transferred downward continuously. Thus, fluctuations were smoothed by the mass transfer process in soil.

Field Experiments at the Kuahuqiao Site Museum
During the field experiment, humidity, air temperatures and soil temperature were monitored (Figure 7a). The measured cumulative amount of dew (Figure 7b) was estimated using known dry soil density (1.802 g/cm 3 ) and the measured moisture content. These original data were cited, and a cumulative dew curve was constructed using Equations (2)-(6) to simulate this field experiment. Simulation was satisfied when δ = 4.58 × 10 −3 and a = 2.91 × 10 3 . Three small concaves can be observed in the simulation curve (Figure 7b) in accordance with the fluctuations of environmental variables, but disappeared in the measured curve. This is explained by the delay in measuring dew amounts. Measured dew amounts were calculated using the water content beneath the soil surface while simulated dew amounts were on air-soil surface. When the surface dew formation stopped, free water transferred downward continuously. Thus, fluctuations were smoothed by the mass transfer process in soil.

Field Experiments in a Desert System
Data from another field experiment [18,19,31] were also simulated. Using the proposed model, a simulation of the cumulative dew amount of dew was performed and the results were shown in Figure 8. pv2 was calculated using Equation (6) because the species of soluble salts were unknown. Fitting parameters used were δ = 3.82 × 10 −2 and a = 1.25 × 10 3 . During the time range 10.0-11.6 h, calculated dew amounts were below zero. They were replaced by zero manually, because the value of water content in soil cannot be negative. This simulation of an outdoor experiment was satisfactory, suggesting that the proposed model is potentially useful beyond heritage preservation.

Field Experiments in a Desert System
Data from another field experiment [18,19,31] were also simulated. Using the proposed model, a simulation of the cumulative dew amount of dew was performed and the results were shown in Figure 8. pv 2 was calculated using Equation (6) because the species of soluble salts were unknown. Fitting parameters used were δ = 3.82 × 10 −2 and a = 1.25 × 10 3 . During the time range 10.0-11.6 h, calculated dew amounts were below zero. They were replaced by zero manually, because the value of water content in soil cannot be negative. This simulation of an outdoor experiment was satisfactory, suggesting that the proposed model is potentially useful beyond heritage preservation.
Water 2021, 13, x FOR PEER REVIEW 10 of 13 Figure 8. Measured cumulative amount of dew collected in a desert system [18] and simulated curve.

Model Validation Using Numerical Simulation
The validity of the proposed model in laboratory experiments has been preliminarily shown in 4.1, by R values over 0.9. Here, the accordance between δ and δe provides more evidence for the assumptions in the proposed model. The difference between fitting parameter δ and thickness of the equivalent air layer δe was analyzed using a two-tailed t-test (shown in Table 2, α = 0.05). For both phases, δ was in the confidence intervals of δe, indicating that there was no significant difference between the two groups. The two parameters were not statistically different.

Model Validation Using Numerical Simulation
The validity of the proposed model in laboratory experiments has been preliminarily shown in 4.1, by R values over 0.9. Here, the accordance between δ and δ e provides more evidence for the assumptions in the proposed model. The difference between fitting parameter δ and thickness of the equivalent air layer δ e was analyzed using a two-tailed t-test (shown in Table 2, α = 0.05). For both phases, δ was in the confidence intervals of δ e , indicating that there was no significant difference between the two groups. The two parameters were not statistically different. Additionally, significant differences can be observed between δ e values in two phases (p = 0.01), indicating that the thickness of the air layer was changed after saturation of soil. This rationalized the division of the two phases of dew formation.

Application Potential and Unexplained Parameters
Setting of the three experiments varied largely. The space of dew formation was centimeter-level in the laboratory experiment, decimeter-level in the Kuahuqiao Museum and open space in desert. The wind causing turbulent flow was zero in the laboratory experiment, strong in the Kuahuqiao Museum (approx. 10 m/s according to product description) and weak in the desert (approx. 2.7 m/s mean annual wind speed [19]). The property of soil was controlled in lab and unknown in situ. Only the water condensation process was tested in our previous study [16], while evaporation was included in desert. As a result, the proposed model has shown application potential in various situations.
The parameters used for model fitting are listed in Table 3. When this model is used to predict dew in a new environment, a preliminary experiment is recommended to determine the values of parameters δ and a. It should be mentioned that the assumption of the proposed model was not fully met in field experiments, where parameter a were extremely high and the maximum differences of vapor pressure were 5176 and 268 Pa. In these cases, neglected factors such as wind and soil suction, may be the driving force of water condensation. As a result, the physical meaning of δ must be reconsidered. It is interesting that simulations using the proposed model still worked well, and the variations of parameter δ seem to be explicable. On the one hand, δ was much smaller at the Kuahuqiao Site. This is reasonable because the electric fan caused intensely turbulent flows. On the other hand, δ was lower in laboratory experiments compared with in a desert; the reason should be the limited space in a desiccator.
Additionally, parameter a, simulating Phase II of laboratory experiments, was close to zero. This can be attributed to the water saturation of soil-the influence of soil suction was prevented after water emerged on the soil surface.
In general, the application potential of the proposed model has been shown, though the physical meanings of parameters in field experiments are obscured. Further study will be needed to give more particular explanations.

Conclusions
In this study, a semi-theoretical model is constructed to calculate dew formation rates. This model assumes a steady air layer at air-solid interface and mass transfer in it is calculated. The model contains two fitting parameters: δ represents thickness of the transfer field and a represents overall influence of other factors. Three variables of the model are easy to measure: relative humidity in air (RH), air temperature (T a ) and the temperature of the soil surface (T s ).
The model yielded satisfactory results in fitting laboratory experiments (R > 0.9), as well as in fitting published field experiments. Moreover, the physical interpretation of parameter δ in laboratory experiments was studied using a numerical simulation.
Models such as the one presented here will be useful for preservation of earthen heritage sites. Further research should be conducted to make the model more practical. For example, more accurate interpretations of parameters can be studied. As an expectation, an automatic system using a practical model can be developed to monitor and control surrounding environments of cultural relics. Artificial dew could be produced to prevent soil desiccation in arid regions or salt accumulation in humid regions, while superfluous water condensation can be avoided before it causes damage.