Profiles of relative permittivity and electrical conductivity from unsaturated soil water content models ANNALS OF GEOPHYSICS, 59, 3, 2016, G0320; doi:10.4401/ag-6990 G0320 Profiles of relative permittivity and electrical conductivity from unsaturated soil water content models Roberta Porretta*, Fabio Bianchi Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy ABSTRACT A mathematical model of water diffusion in the vadose zone has been im- plemented for different types of soil textures in order to determine the soil water content (SWC) profiles in dependence of depth and time. From these profiles, obtained for different soils, we derived the characteristic electrical parameters, such as relative permittivity (fr ) and electrical con- ductivity (v), and their variation in time, employing empirical relations available in the scientific literature. The simulation through mathemati- cal models has been performed taking into account different types of soils characterized by the percentage composition of sand, clay and silt in the textural triangle, which provides some physical and chemical properties that affect the water retention in the soil. The resulting simulated profiles of SWC and consequently permittivity and conductivity profiles, span over a certain range of values suggesting the best techniques and the lim- its in geophysical investigation. Moreover this a-prior knowledge helps in the elaboration and interpretation of permittivity and conductivity data obtained by the measurements. Permittivity and conductivity profiles are particularly useful in some environmental applications when the soil tex- tures are (or supposed to be) known as in the typical case of landfill leachate dispersion. Since the soil textures beneath or nearby a disposal waste are characterized by a SWC, the values of fr and vat various depth can be directly inferred. 1. Introduction Electrical conductivity (v) and relative dielectric per- mittivity (fr ) are two independent physical properties that characterize the behavior of soil when excited by electric currents or electromagnetic fields [Settimi 2011]. The electrical conductivity and dielectric permittivity of the constituent minerals are the main controlling factors of the electrical properties of soils [Kriaa et al. 2014]. The derivation of these quantities related to the water con- tent provides useful information regarding the practical use of geophysical investigation (GPR, TDR, etc.), ap- plied to unsaturated subsurface vadose zone of medium. This zone has relevant importance in different ambits and especially in geophysical prospection [Huisman et al. 2003]. The unsaturated portion of soil concerns the first meters beneath the ground level down to the water table, where the saturated zone begins. In this section hydro- geological properties (physical and chemical) are impor- tant since they dominate the electric properties of soil, described by the two quantities v and fr . These last quan- tities can be directly measured through several techniques or theoretically deduced (constitutive relations) in a ho- mogeneous medium. In this paper these parameters are estimated through soil hydrologic models [Huisman et al. 2003]. In the latter case the SWC in vadose zone plays an important role to determine both relative permittivity and electrical conductivity since they are strongly de- pendent on it. The different soil textures can retain dif- ferent percentage of water, due to its porosity, under the action of two agents: the gravity potential and the ma- trix potential. The diffusive equation is driven by these quantities and cannot be solved without models. These models need in input physical parameters that quantify the textures of soil. In the past decades authors pro- posed different numerical parameters [Mualem 1976, van Genuchten 1980] to solve the so-called Richards equation that well describes the water diffusion in the vadose zone. Once the input parameters for the particular type of soil are established it is possible to numerically calculate the SWC, which is function of time and depth. SWC is the most important hydro-geological parameters that jointly to salinity and porosity affects the electrical properties in the unsaturated soils. This paper deals with simulations of soil model in order to evaluate the SWC and consequently through Topp and Archie empirical relationships to determine the electrical parameters of soil. In this specific case, we obtained 3D profiles that show temporal and spatial variation of the electrical properties. The study focuses on soil investigation theoretical approach, that could be a useful tool in some environmental and geophysi- Article history Received February 19, 2016; accepted May 17, 2016. Subject classification: Soil texture, Diffusion, Water content, Conductivity/permittivity profiles. cal applications. Furthermore it could represent a pre- liminary phase which precedes the measurements and facilitates the interpretation of the results (GPR, TDR, ERT and other similar geophysical investigation tech- niques e.g. after raining survey or in landfill control). 2. Water content models Water flows in saturated porous medium and water content retention in different kind of soils has been sub- ject of many studies. To describe this phenomenon in the vadose zone Richards proposed an equation which is a combination between the flux equation and mass con- servation, and it characterizes the SWC trend over time [Richards 1931]. The software developed to derivate SWC presented in this paper it’s based on well-known litera- ture models such as van Genuchten, Mualem, Brooks and Corey. The introduction of these models was neces- sary, due to the impossibility to analytically solve the Richards’ equation. The water diffusion equation into unsaturated soil in one-dimensional terms can assume the following form: where i is the SWC [m3/m3], H represents the hy- draulic head: H = } + z , being } the matrix potential and z the elevation above a reference plan (level). K is hydraulic conductivity in (m/s). The above is the so- called Richard’s equation that has no closed form solu- tion because K is depend dependent on i. The substitution of the hydraulic head H with }+ z , leads to the so-called Richard’s mixed form: (2) Here please note that K is the hydraulic conductivity in [cm/day] and i is the soil water content [cm3/cm3]. It describes the movement of the fluid (fluid motion) in a non-saturated and porous medium. It is a non lin- ear partial differential equation and, as such, it has no analytical closed form solution and must be integrated numerically. Furthermore, K(i) is not a constant in the strict sense, as it depends on the same i, hence the problem need an appropriate model to be treated. Here it is pre- sented in one-dimensional term along the vertical z. Equations (2) and (3) have a two dependent variables i and H (or }), so the resolution of the two forms re- quires the definition of a constitutive equation that binds i and } as in Figure 1 [Mualem 1976]. From USDA (United States Departement of Agriculture) triangle textures shown in Figure 2 typical van Genuchten pa- rameters are reported in Table 1 [Leij et al. 1996]. As mentioned earlier K is not a constant, because of its dependence on i and for this reasons, a model for each different soil textures is required. This demands a constitutive equation that takes into account the water retention curve. Such a curve depends on the soil pososity, or, even better depends on the effective saturation Se [Leij et al. 1997]. The latter can be expressed as: (4) where: is is the water content in saturated conditions and ir is the residual water content. It means that if Se = 0 the fluid can be found only in gaseous phase, t Kz2 2 2 2 2 2i = S Hz2 2 X t z z z z z K zK 12 2 2 2 2 2 2 2 2 2 2 2i W W = ++ =S SX X" "% % t z K2 2 2 2 i i = Q" z 12 2 i W +SV X% Se ri i= -Q V rsi i-Q V PORRETTA AND BIANCHI 2 (1) (3) Figure 1. Water retention curve that links the two variables i or } depending on the type of soil [Tuller and Or 2005]. Figure 2. USDA (United States Departement of Agriculture) tex- tures triangle. 3 while if Se = 1 we are in saturated conditions. Because the hydraulic conductivity K depends on Se, therefore the introduction of models [Brooks and Corey 1964, Mualem 1976] that take into account such a depend- ence K→K (Se) is necessary. In this work the following model proposed by van Genuchten has been adopted: (5) Such a model [van Genuchten 1980] is widely em- ployed in this field. In this work was also introduced the dependence on the matrix potential like: (6) with tabulated parameters a, n = 1 − 1/m (7) Equation (3) after some few arrangements be- comes Equation (8) is numerically solved in Matlab pro- gram through the van Genuchten’s model, given in re- lations (5), (6) and (7). This allows to obtain both i(z,t) and }(z,t) variables in function of time and depth. Since the implemented model it’s based on equations having validity in vadose zone the simulated depth can only reach the water table. Limits on the depth param- eter are intrinsic in the matrix potential that cannot act on more than 10 m above the water table in ideal con- ditions, which reflects in a 2-5 m action in real condi- tions depending on soil texture. Also the heterogeneity of the soil has been strongly simplified and the missing data interpolated with a statistical algorithm. 3. Empirical relations to derive relative permittivity of the soil Many physical and empirical models have been suggested for the evaluation of the i−fr relationship in the literature. The relative permittivity is the main elec- tric quantity employed to define SWC in the soil since it can be easily measured through various techniques. The relation that ties these two quantities was empirically established [Topp et al. 1980]. This equation is valid for a wide range of mineral soils and independent from soil bulk density, ambient temperature, and salt content. This led many authors to the use of the term “universal” for this equation with appropriate caveat that “in or- ganic soils or heavy clay soils problems arise which may require site-specific calibration” [Cosenza et al. 2003, and the reference therein]. In this specific study the Topp’s formula [Topp et al. 1980] has been used . From well-known water content profiles, obtained by a mod- eling software [Bianchi et al. 2015], Topp’s relation (9) allows to derivate the relative permittivity values of soil from the SWC. This empirical model was generated using time-domain reflectometry (TDR) at a frequency between 1 MHz and 1 GHz to measure fr for several mineral soils. The estimated error in this model is 0.013. General Topp’s model can be expressed: (9) where fr is the relative permittivity and i is the soil vol- umetric water content. They also provided another in- verse relation as follows: z K zt 12 2 2 2 2 2 2 2 i i W W W = +Q SV X" % nm h 1 n mV s r nV 1 1 2 2i a} i i a a W + = - - + -Q Q QV mV nV Se 1 1= +Q n mVa} mV K S K Se e=Q V S1 /e m mW1 2%- R mW" 2% . . .3 03 9 3 146 76 7r 2 3 f i i i= + + - ELECTRICAL PROFILES FROM SWC MODELS (8) Textural class N ir [cm3/cm3] is [cm3/cm3] a [1/cm] n Sand 126 0.058 0.37 0.035 3.19 Loamy sand 51 0.074 0.39 0.035 2.39 Sandy loam 78 0.067 0.37 0.021 1.61 Loam 61 0.083 0.46 0.025 1.31 Silt 3 0.123 0.48 0.006 1.53 Silt loam 101 0.061 0.43 0.012 1.39 Sandy clay loam 37 0.086 0.40 0.033 1.49 Clay loam 23 0.129 0.47 0.030 1.37 Silty clay loam 20 0.098 0.55 0.027 1.41 Silty clay 12 0.163 0.47 0.023 1.39 Clay 25 0.102 0.51 0.021 1.20 Table 1. Typical van Genuchten model with parameters (a, n) including residual (ir) and saturated (is). (10) A specific calibration is needed for soils with higher water content or organic matter: (11) for organic soil; (12) for 450 µm glass beds [Mukhlisin and Saputra 2013]. The model relationships works better for frequen- cies around 100 MHz. At higher frequencies and mois- ture contents near to saturation (i~0.4) the Topp-model over-predicts relative permittivity by up to 20%. At very low water contents the Topp-model doesn’t perform well, particularly for soils with a large clay content [van Dam et al. 2005]. However, later studies have shown the dependency of the i−fr relationship on clay content [Persson et al. 2000, Bouksila et al. 2008] and mineral- ogy [Cosenza and Tabbagh 2004], organic matter and porosity or soil density [Malicki et al. 1996 and Persson et al. 2002], and soluble salt content [Dalton 1992, Nadler et al. 1999, Persson et al. 2000]. Ledieu et al. [1986] have shown another form of empirical equations to describe the i−fr relationship, that can be used to expand the Topp model for higher water content [van Dam et al. 2005] of the form: (13) where b0 and b1 are two empirical parameters depend- ing on soil type. In the same paper Ledieu et al. [1986] found the following relationship: (14) This relationship appears to work better for most mineral non-magnetic soils over a range of frequencies between 1 MHz and 10 GHz [Hamed et al. 2006]. Roth et al. [1990] for mineral soil proposed this empirical relationship: (15) While for organic soil and material is: (16) the error estimations of these equations for mineral soil and organic soil are 0.015 and 0.035 cm3/cm3, respec- tively [Mukhlisin and Saputra 2013]. Although a comparison of all the above relation- ships showed a similar trend (Figure 3), Topp’s formula (blue curve) resulted the relationship that best fitted the input data, since is calibrated for soils that can hold a water content up to 50%. Indeed in unsaturated zone the pores in the soil can be filled not only by water but also by air. 4. Empirical relations to derive conductivity of the soil Another major controlling factor, of the soil elec- trical proprieties estimated in this study, is electrical . . . . w 5 3 10 2 92 10 5 5 10 4 3 10 r r r r 2 2 4 6 # # # # i f f f f =- + - + - - - - Q V . .1 74 0 34 135 55r 2 3 f i i i= - + - .. .13 57 31 7 1 68 2r 2 3 f i i i= - + - bbr 0 1f = + i .0 1138 .r 0 5 i f= .0 17585 - . . . .0 0728 0 044 0 00195 0 0000361r r r 2 3 i f f f=- + - + . . . .0 0233 0 0285 0 000431 0 00000304r r r 2 3 i f f f=- + - + PORRETTA AND BIANCHI 4 Figure 3. Empirical relations comparison. 5 conductivity. Although usually geoelectric investiga- tions are focused on the measurement of soil resistivity, in this paper conductivity is preferred to resistivity (i.e. the inverse quantity), since prevalently we refer to rain water data given in terms of conductivity entering in the following equations. The empirical relation to esti- mate the electrical conductivity both for saturated and unsaturated soil is the Archie’s law [Shah and Singh 2005] also spread in the literature. This last one connect electrical conductivity of sat- urated rocks (vo ) to the conductivity of the electrolyte pore of soils (vw ) , it can be expressed as: where vqo is the conductivity contribution which char- acterizes high clay percentage soils. It can only be neg- lected for high saline water content filling high porosity soils. F is the formation factor which is equal to: (18) The parameter z is the porosity, and m is the ce- mentation exponent which increases with compaction, cementation and consolidation; it varies between 1.3 and 2.5. Unconsolidated sands have values in the range between 1.3 and 1.5 [Dannowski and Yaramanci 1999, and the reference therein]. Generally, the resistivity value is greatly influenced by basic soil characteristics variation such as fraction of solid, air and water. Ac- cording to [Griffiths and King 1981, Telford et al. 1990], resistivity value was highly influenced by pore fluid and grain matrix of geomaterials. In compact condition, it was found that at high compaction values corresponds lower resistivity value (or higher electrical conductivity). The volumes of pore in compact condition were reduced and cause the current to easily propagate especially dur- ing the existing of water [Abidin et al. 2013]. Another form has been proposed, known as Archie’s second law: which has been developed as an extension for unsatu- rated rocks and soils including the I factor (20) The saturation index I depends on the degree of saturation Sw and the saturation exponent n. This latter was observed to be about 2 for consolidated rocks and to range from 1.3 to 2 for unconsolidated sands [Lesmes 2005, and the reference therein]. Archie’s second law (19) has been chosen and implemented in the modeling software to obtain profiles of electrical conductibility trend over time which characterizes the vadose zone. This physical quantity depends on a number of pa- rameters as soil texture, water holding capacity, organic matter, salinity and ions exchange capacity. The con- ductivities of water filling pores is an important factor in the process of electric current flow through the soil; especially when a quantity of salts are dissolved in it. 5. Simulations and results The software estimation of the effective water content, in a porous soil is determined by Richard’s equation (Equation 3). The latter is a partial differential equation that can only be solved with the implementa- tion of the van Genuchten soil models. Hence water content is strongly affected by physical property of soil material (porosity, capillarity, etc.) and fluid properties (viscosity). Indeed at equal initial conditions, retention water content only depends on matrix potential in op- position to the effect of gravity, preventing leaching through the ground till the water table. In the per- formed software simulation, specific initial condition (infiltration velocity), boundary conditions (water table depth, ground level and saturated zone hydraulic head) and simulation parameters as time and depth have been taken into account. Infiltration is the volume of water passing into the soil per unit of area per unit of time. Simplifying it has the dimensions of velocity [m3 m-2 s-1]. Depth is in- tended not only as the depth of the water table, or line of saturation, but also the points in the modeled spaceF I 1 1 o wv v= w qov+ I Sw n= - F 1 o wv v= + oqv F mz= - ELECTRICAL PROFILES FROM SWC MODELS (17) (19) Soil type Porosity [cm³/cm³] Infiltration velocity [cm/day] vw [mS/m] Cementation exponent Saturation exponent Water table depth [cm] Time [day] Sand 0.375 4.010 5 1.4 2 200 4 Silt 0.489 1.688 5 1.4 2 200 4 Clay 0.459 1.603 5 2.5 2 200 4 Table 2. Simulation parameters. where the solution is calculated. Hydraulic head is an equivalent measure of the pressure expressed in height of water column, and it is the sum of pressure head (fluid internal pressure) and elevation head (pressure due to gravity). Following a summarizing Table 2 of the main simulation coefficients, parameters and physical quan- tities that were used to determinate the electrical con- ductivity, permittivity and water content relative to the soil textures object of this paper. The simulation through mathematical models of the WC profiles has been reported, respectively for the most common texture soils types (Figures 4, 5, 6). The picture obtained by the modeling software, shows range of water content up to 50%, and the relative per- mittivity and electrical conductivity values linked to it, plotted in a 3D space. A time of four days has been cho- sen because sand soils after that amount of time lose all the water retained in it, and resulted completely dry. Due to this quite fast change in SWC, the electrical properties are also subject to a rather remarkable vari- ation. So the results of the simulation show how taking a measure on a certain soil texture (after-rain survey) could be misleading if the soil retention is not taken in consideration. From Figure 4 it is possible to see how, in a typical clay soil, water content in the topsoil layer changes very little due to its high degree of retention. In the other hand, sand soil (Figure 5) with medium- size texture have a low moisture holding capacity, which results in a seep through the ground downwards. Silt take up an intermediate behavior between sand and clay soils (Figure 6). Topp’s derivated relative permit- tivity as shown below for clay, sand and silt (Figures 7, 8, 9) and Table 3 shown an average of relative permit- tivity values (fr ). PORRETTA AND BIANCHI 6 Figure 4. Water content section in clay soil textures. Figure 5. Water content section in sand soil textures. 7 ELECTRICAL PROFILES FROM SWC MODELS Figure 6. Water content section in silt soil textures. Figure 7. Relative permittivity section in clay soil textures. Figure 8. Relative permittivity section in sand soil textures. Since relative permittivity is an electrical physical characteristic which exclusively depends on inner water content, for high degree of moisture content, relative permittivity in simulation test increase. Indeed, clay soils texture are characterized by higher relative permittivity values than sand soils, which gener- ally shows much lower relative permittivity. Silt soils, show a range of values between 0 and 24. Profiles for the electrical conductivity have been present in according to Archie’s second law (Equation 19). Also in this case, the profiles were obtained simulating a time of four days as shown in Figures 10, 11 and 12. Since electrical conductivity is linked to total dis- solved solids (TDS) and increase considerably in func- tion of this quantity. To simplify the model, a fixed total solute concentration has been considered, which results PORRETTA AND BIANCHI 8 Figure 9. Relative permittivity section in silt soil textures. i [cm3/cm3] fr 0.10 5.343 0.15 7.451 0.20 10.116 0.25 13.281 0.30 16.889 0.35 20.881 0.40 25.201 0.45 29.791 0.50 34.592 Figure 10. Electrical conductivity section in clay soil textures. Table 3. Relative Permittivity and corresponding SWC values ac- cording to Topp relation. 9 in a typical rain water electrical conductivity value of about 5 mS/m. Furthermore, in the specific case of clay model, the contribution of the surface grains or particle conductivity (of about 90.94 mS/m) has been taken into consideration. Clay surface’s mineral parti- cles influence considerably the electrical conductivity of this specific texture which results in an increase of the total current flow through soil as shown in Figure 10 and Table 4. 6. Conclusions This paper dealt with the estimation of relative permittivity and electrical conductivity profiles inferred from SWC models for different soil textures. As well known the different percentages of the soil components, namely sand, clay and silt, described trough a numeri- cal model allows to solve the Richards diffusion equa- tion and determine the SWC in dependence of depth and time. Such quantity is mainly responsible for the val- ues assumed by the relative permittivity and conductivity along depth and time. Because of seldom geophysicists have a prior knowledge of the subsoil electrical parame- ters even in the case of homogeneous soil’s composi- tion these simulation yields of SWC along the vadose zone once recognized some characteristic properties of soil textures. Conversely, more commonly, this infor- mation can be used in the opposite sense inferring the SWC from the relative permittivity and conductivity. The prior knowledge of the permittivity and conduc- tivity modelled profiles could help geophysicists and op- erators to better interpret the field measurement’s results as in application of landfill monitoring where the soil tex- tures are known. The modelled SWC and corresponding electric parameters are useful to choice the suitable tech- ELECTRICAL PROFILES FROM SWC MODELS Figure 11. Electrical conductivity section in sand soil textures. Figure 12. Electrical conductivity section in silt soil textures. niques to be employed and the optimal mathematical inversion algorithm to retrieve the electric profiles. The results of these simulations as reported in Tables 3 and 4, and Figures 7, 8 and 9 show the strong dependence of the permittivity from the SWC. Other soil parameters have negligible influence on this quantity. Different is the case in which we consider the conductivity. In fact SWC dependence is important even if salinity, temperature, etc. are less relevant but still important for soil conductivity. The estimation of this two profiles have relevance especially when the soil tex- tures are well characterised. In such a case this allows to simulate SWC through the above discussed model and consequently estimate these two electric quantities in function of depth and time in unsaturated soil. A com- parison with others studies carried out previously [Zhou et al. 2001, Michot et al. 2003, Schwartz et al. 2008], showed the validity of the inverse method to convert soil moisture content to electrical resistivity using ERT with a modified form of the Archie’s law. Starting from a given water content, electrical resistivity trend in time can be traced. Application in landfill monitoring is one the fruitful application since leachate dispersion can be profitably inferred through the measurements of the soil electric parameters fr and v. References Abidin, M.H.Z., F. Ahmad., D.C. Wijeyesekera., R. Saad and M.F.T. Baharuddin (2013). Soil Resistivity Measurements to Predict Moisture Content and Density in Loose and Dense Soil, Applied Mechan- ics and Materials, 353/356, 911-917. Bianchi, F., M. Chiappini and R. Giordano (2015). For- mazione, diffusione e distribuzione del percolato: modelli matematici integrati del progetto SIGLOD (Formation and distribution of leachate: integrated mathematical models of SIGLOD project), Quader- ni di Geofisica, 128, 25 p., ISSN 1590-2595, (in Italian). Bouksila, F., M. Persson, M.R. Berndtsson and A. Bahri (2008). Soil water content and salinity determina- tion using different dielectric methods in saline gyp- siferous soil, Hydrological Sciences Journal, 53 (1), 253-265. Brooks, R.H., and A.T. Corey (1964). Hydraulic prop- erties of porous media, Hydrology Papers, 3, Col- orado State University, Fort Collins, Colorado. Cosenza, P., C. Camerlynck and A. Tabbagh (2003). Dif- ferential effective medium schemes for investigating the relationship between high-frequency relative di- electric permittivity and water content of soils, Water Resources Research, 39 (9). Cosenza, P., and A. Tabbagh (2004). Electromagnetic determination of clay water content: role of the mi- croporosity, Applied Clay Science, 26, 21-36. Dalton, F.N. (1992). Development of time-domain re- flectrometry for measuring soil water content and bulk soil electrical conductivity, In: G.C. Topp and W.D. Reynolds (eds.), Advances in Measurement of Soil Physical Properties: Bringing Theory into Prac- tice, Soil Science Society of America Spec. Publ. no. 30, Madison, Wisconsin, USA, 143-167. Dannowski, G., and U. Yaramanci (1999). Estimation of water content and porosity using combined radar and geoelectrical measurements, Technical Univer- sity of Berlin, Dept. of Applied Geophysics, Acker- str. 71-76, D-13355 Berlin, Germany. Griffiths, D.H., and R.F. King (1981). Applied Geophysics for Geologist and Engineers, The Element of Geo- physical Prospecting, Pergamon Press, Oxford. Hamed, Y., S. Ghada and M. Persson (2006). Evaluation of the WET sensor compared to time domain Re- flectometry, Hydrological Sciences Journal, 51 (4), 671-681. Huisman J.A., S.S. Hubbard., J.D. Redman and A.P. PORRETTA AND BIANCHI 10 i [cm3/cm3] v Sand [mS/m] v Silt [mS/m] v Clay [mS/m] 0.1 9.010E-05 7.680E-05 9.094E+01 0.15 2.026E-04 1.728E-04 9.094E+01 0.2 3.603E-04 3.072E-04 9.094E+01 0.25 5.629E-04 4.800E-04 9.094E+01 0.3 8.106E-04 6.912E-04 9.094E+01 0.35 1.103E-03 9.409E-04 9.094E+01 0.4 1.441E-03 1.229E-03 9.094E+01 0.45 1.824E-03 1.555E-03 9.094E+01 0.5 2.252E-03 1.920E-03 9.094E+01 Table 4. Electrical conductivity and corresponding SWC values according to Archie relation. 11 Annan (2003). A Measuring Soil Water Content with Ground Penetrating Radar: A Review, Vadose Zone Journal, 2, 476-491. Kriaa, A., M. Hajjib, F. Jamoussib and A.H. Hamzaouib (2014). Electrical Conductivity of 1:1 and 2:1 Clay Minerals, Surface Engineering and Applied Electro- chemistry, 50 (1), 84-94. Ledieu, J., P. De Ridder, P. De Clerk and S. Dautrebande (1986). A Method of measuring soil moisture by Time Domain Reflectometry, Journal of Hydrology, 88, 319-328. Leij, F.J., W.J. Alves., M.Th. van Genuchten and J.R. Williams (1996). The UNSODA unsaturated hy- draulic database. EPA/600/R-96/095, U.S. Environ- mental Protection Agency, Cincinnati, OH. Leij, F.J., W.B. Russell and S.M. Lesch (1997). Closed- Form expression For Water Retention and conduc- tivity data, Ground Water, 35 (5), 848-858. Lesmes, D.P., and P. Shmulik Friedman (2005). Rela- tionships between the electrical and hydrogeologi- cal properties of rocks and soils, In: Y. Rubin and S.S. Hubbard (eds.), Hydrogeophysics, Springer, Dordrecht, The Netherlands, Ch. 4, 87-128. Malicki, M.A., R. Plagge and C.H Roth (1996). Im- proving the calibration of dielectric TDR soil mois- ture determination taking into account the solid soil, European Journal of Soil Science, 47, 357-366. Michot, D., Y. Benderitter., A. Dorigny, B. Nicoullaud, D. King and A. Tabbagh (2003). Spatial and tempo- ral monitoring of soil water content with an irrigated corn crop cover using electrical resistivity tomogra- phy, Water Resources Research, 39, 1138. Mualem, Y. (1976). A new model for predicting the hy- draulic conductivity of unsaturated porous media, Water Resources Research, 12 (3), 513-522. Mukhlisin, M., and A. Saputra (2013). Performance Evaluation of Volumetric Water Content and Rela- tive Permittivity Models, The ScientificWorld Jour- nal, Vol. 2013, Article ID 421762, 7 p. Nadler, A., A. Gamliel and I. Peretz (1999). Practical as- pects of salinity effect on TDR-measured water con- tent: a field study, Soil Science Society of America Journal, 63, 1070-1076. Persson, M., R. Berndtsson, S. Nasri, J. Albergel, P. Zante and Y. Yumegaki (2000). Solute transport and water content measurements in clay soils using time domain reflectometry, Hydrological Sciences Jour- nal, 45, 833-847. Persson, M., B. Sivakumar, R. Berndtsson, O.H. Jacob- sen and P. Schjønning (2002). Predicting the dielec- tric constant-water content relationship using artificial neural networks, Soil Science Society of America Journal, 66, 1424-1429. Richards, L.A. (1931). Capillary conduction of liquids through porous mediums, Physics, 1 (5), 318-333. Roth, K., R. Schulin, H. Fluhler and W. Attinger (1990). Calibration of time domain reflectometry for water content measurement using a composite dielectric approach, Water Resources Research, 26, 2267-2273. Schwartz, B.F., M.E. Schreiber, P.S. Pooler and J.D. Rim- stidt (2008). Calibrating access-tube time domain reflectometry soil water measurements in deep het- erogeneous soils, Soil Science Society of America Journal, 72 (4), 917-930. Settimi, A. (2011). Performance of electrical spectroscopy using a RESPER probe to measure salinity and water content of concrete and terrestrial soil, Annals of Geophysics, 54 (4), 400-413. Shah, P.H., and D.N. Singh (2005). Generalized Archie’s Law for Estimation of Soil Electrical Conductivity, Journal of ASTM International, 2 (5), 1-20. Telford, W.M., L.P. Geldart and R.E. Sheriff (1990). Ap- plied Geophysics, 2nd ed., Cambridge University Press, Cambridge, 792 p. Topp, G.C., J.L. Davis and A.P. Annan (1980). Electro- magnetic determination of soil water content: Measurements in coaxial transmission lines, Water Resources Research, 16, 574-582. Tuller, M., and D. Or (2005). Retention of water in soil and the soil water characteristic curve, In: D. Hillel (ed.), Encyclopedia if Soils in the Environment, El- sevier. van Dam, R.L., B. Borchers and J.M.H. Hendrickx (2005). Methods for prediction of soil dielectric properties: a review, In: Proceedings of SPIE - The International Society for Optical Engineering 5794, 188-197. van Genuchten, M.Th. (1980). A closed-form equation for predicting the hydraulic, conductivity of unsat- urated soils, Soil Science Society of America Journal, 44, 892-898. Zhou, Q.Y., J. Shimada and A. Sato (2001). Three-dimen- sional spatial and temporal monitoring of soil water content using electrical resistivity tomography, Water Resources Research, 37 (2), 273-285. * Corresponding author: Roberta Porretta, Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy; email: roberta.pr@live.com. © 2016 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. ELECTRICAL PROFILES FROM SWC MODELS << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles false /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Warning /CompatibilityLevel 1.3 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.1000 /ColorConversionStrategy /LeaveColorUnchanged /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams true /MaxSubsetPct 100 /Optimize false /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments false /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile (None) /AlwaysEmbed [ true /AndaleMono /Apple-Chancery /Arial-Black /Arial-BoldItalicMT /Arial-BoldMT /Arial-ItalicMT /ArialMT /CapitalsRegular /Charcoal /Chicago /ComicSansMS /ComicSansMS-Bold /Courier /Courier-Bold /CourierNewPS-BoldItalicMT /CourierNewPS-BoldMT /CourierNewPS-ItalicMT /CourierNewPSMT /GadgetRegular /Geneva /Georgia /Georgia-Bold /Georgia-BoldItalic /Georgia-Italic /Helvetica /Helvetica-Bold /HelveticaInserat-Roman /HoeflerText-Black /HoeflerText-BlackItalic /HoeflerText-Italic /HoeflerText-Ornaments /HoeflerText-Regular /Impact /Monaco /NewYork /Palatino-Bold /Palatino-BoldItalic /Palatino-Italic /Palatino-Roman /SandRegular /Skia-Regular /Symbol /TechnoRegular /TextileRegular /Times-Bold /Times-BoldItalic /Times-Italic /Times-Roman /TimesNewRomanPS-BoldItalicMT /TimesNewRomanPS-BoldMT /TimesNewRomanPS-ItalicMT /TimesNewRomanPSMT /Trebuchet-BoldItalic /TrebuchetMS /TrebuchetMS-Bold /TrebuchetMS-Italic /Verdana /Verdana-Bold /Verdana-BoldItalic /Verdana-Italic /Webdings ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 150 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.10000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 150 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.10000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.08250 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile (None) /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName (http://www.color.org) /PDFXTrapped /Unknown /CreateJDFFile false /SyntheticBoldness 1.000000 /Description << /ENU (Use these settings to create PDF documents with higher image resolution for high quality pre-press printing. The PDF documents can be opened with Acrobat and Reader 5.0 and later. These settings require font embedding.) /JPN /FRA /DEU /PTB /DAN /NLD /ESP /SUO /NOR /SVE /KOR /CHS /CHT /ITA >> >> setdistillerparams << /HWResolution [2400 2400] /PageSize [595.000 842.000] >> setpagedevice