New Approach for Calculating the Effective Dielectric Constant of the Moist Soil for Microwaves

Microwave remote sensing techniques are used, among others, for temporally and spatially highly-resolved observations of land-surface properties, e.g., for the management of agricultural productivity and water resource, as well as to improve the performances of numerical weather prediction and climate simulations with soil moisture data. In this context, the effective dielectric constant of the soil is a key variable to quantify the land surface properties. We propose a new approach for the effective dielectric constant of the multiphase soil that is based on an arithmetic average of the dielectric constants of the land-surface components with damping. The results show, on average, better agreement with experimental data than previous approaches. Furthermore, the proposed new model overcomes the theoretical limitation of previous models in the incorporation of non-physical parameters to simulate measured data experimentally with satisfactory accuracy. For microwave remote sensing such as SMAP (Soil Moisture Active Passive), SMOS (Soil Moisture and Ocean Salinity) and AMSR-E (Advanced Microwave Scanning Radiometer for EOS), the physical-based model in our study showed a 23–35% RMSE (root-mean-square error) reduction compared to the most prevalent refractive mixing model in the prediction of the dielectric constant for the real and imaginary part, respectively. Furthermore, in radiowave bands used in portable soil sensors such as TDR (time-domain reflectometer) and GPR (ground-penetrating radar) the new dielectric mixing model reduced RMSE by up to 53% in the prediction of the dielectric constant. We found that the permittivity over the saturation point (porosity of dry soil) has a very different and varying pattern compared to that measured in the unsaturated condition. However, in our study, this pattern was mathematically derived from the same mixing rule applied for the unsaturated condition. It is expected that the new dielectric mixing model might help to improve the accuracy of flood monitoring by satellite.


Introduction
In contrast to radiation in the visible and IR spectral regions, microwaves generally penetrate non-precipitating clouds and the gaseous atmosphere without strong interaction (with the exception only of some microwave regions at which absorption by some gases takes place).Therefore, microwaves allow for reliable information about the soil properties and temperatures of the terrestrial surface, as well as of vegetation independent of cloud cover.Furthermore, microwave remote sensing data enable us to reduce the uncertainties of remote sensing observations made at other wavelengths regarding the land surface for improved estimations of atmospheric trace gases, aerosols and clouds.
In summary, we can thus state that the understanding of the interaction of microwaves is important, not only for land surface remote sensing, but also for observations of the Earth system as a whole.
In modern space-borne microwave remote sensing [1], several bands are used: L-band (1-2 GHz), C-band (4-8 GHz), X-band (8-12 GHz) and K-band .The accurate computation of the effective dielectric constant is essential for both passive and active microwave sensors.Furthermore, it is also critical for the analysis of materials in materials sciences [2].The first dielectric mixing formulas were proposed for cavities, which are (hypothetical) spheres [3,4], monodispersed spheres [5], polydispersed spheres [6], non-spherical [7], as well as for non-spherical nanoporous media and nanoparticles [8,9].However, these mixing models contain an innate limitation for complex multiphase materials such as moist soil: the practical design of the dielectric mixing model based on a microgeometry approach relies on empirical adjustments [2].When a material is exposed to an electric field, its dielectric constant describes the interaction.Thus, the remote sensing of land-surface properties, such as soil moisture, requires an effective operator for the computation of the effective dielectric constant (Figure 1).aerosols and clouds.In summary, we can thus state that the understanding of the interaction of microwaves is important, not only for land surface remote sensing, but also for observations of the Earth system as a whole.
In modern space-borne microwave remote sensing [1], several bands are used: L-band (1-2 GHz), C-band (4-8 GHz), X-band (8-12 GHz) and K-band .The accurate computation of the effective dielectric constant is essential for both passive and active microwave sensors.Furthermore, it is also critical for the analysis of materials in materials sciences [2].The first dielectric mixing formulas were proposed for cavities, which are (hypothetical) spheres [3,4], monodispersed spheres [5], polydispersed spheres [6], non-spherical [7], as well as for non-spherical nanoporous media and nanoparticles [8,9].However, these mixing models contain an innate limitation for complex multiphase materials such as moist soil: the practical design of the dielectric mixing model based on a microgeometry approach relies on empirical adjustments [2].When a material is exposed to an electric field, its dielectric constant describes the interaction.Thus, the remote sensing of land-surface properties, such as soil moisture, requires an effective operator for the computation of the effective dielectric constant (Figure 1).

Figure 1.
Measuring principles of bare soil properties (black boxes) by satellite and airborne remote sensing in the microwave region (A) and TDR (time-domain reflectometer) and GPR (ground-penetrating radar) (B).The dielectric constant (hexagon) is the link between the measured parameters of the sensors (brightness temperature (TB) and refractive index, respectively) and the targeted soil properties.TDR (time-domain reflectometer) and GPR (ground-penetrating radar) measure the refractive index, compute the effective dielectric constant and quantify the soil water content by using the separately-obtained temperature and soil texture information.In contrast to this, airborne and space-borne remote sensing instruments measure the brightness temperature TB (Box A in Figure 1).To retrieve the soil moisture content from the measured TB, the traditional retrieval method requires ancillary information to account for soil temperature and soil texture, as well as for the effects of vegetation.Studies have attempted to obtain this ancillary data for vegetation from the vegetation parameter b [10], multi-frequency microwave sensor measurements [11][12][13], the NDVI (normalized differenced vegetation index) [14] or the MPDI (microwave polarization difference index) [15][16][17].However, the measurement of land-surface properties is already very uncertain for bare soil without vegetation.Hence, for this simpler situation, the ancillary information about soil temperature and soil texture is already quite critical for retrieving the soil moisture content.As a consequence, we TDR (time-domain reflectometer) and GPR (ground-penetrating radar) measure the refractive index, compute the effective dielectric constant and quantify the soil water content by using the separately-obtained temperature and soil texture information.In contrast to this, airborne and space-borne remote sensing instruments measure the brightness temperature TB (Box A in Figure 1).To retrieve the soil moisture content from the measured TB, the traditional retrieval method requires ancillary information to account for soil temperature and soil texture, as well as for the effects of vegetation.Studies have attempted to obtain this ancillary data for vegetation from the vegetation parameter b [10], multi-frequency microwave sensor measurements [11][12][13], the NDVI (normalized differenced vegetation index) [14] or the MPDI (microwave polarization difference index) [15][16][17].However, the measurement of land-surface properties is already very uncertain for bare soil without vegetation.Hence, for this simpler situation, the ancillary information about soil temperature and soil texture is already quite critical for retrieving the soil moisture content.As a consequence, we focus in this study on the bare soil.The simulation results are validated with TDR measurements.Future studies may extend the new dielectric mixing model presented here to ground vegetation and canopy layers.
The dielectric mixing model introduced here can be integrated into land surface models.All of the information that is necessary for calculating TB is available in the land surface model (such as soil temperature, soil moisture [18] and soil texture [19]).Thus, the new mixing model can be used in the future within a forward operator for determining TB.Recently, new "physically-based" radiative transfer models were proposed [14,16,20] for the retrieval of near surface soil moisture from passive microwave measurements.These can be better combined with data assimilation [21][22][23], especially for the estimation of the root zone soil moisture using SMOS [24,25] and SMAP [26,27].By combining our dielectric mixing model with the radiative transfer models, the TB can be obtained more accurately in the microwave spectral region, which is of course beneficial not only for the data assimilation scheme but also the field campaign for the calibration and validation such as SMOSREX [28] or SMAPVEX [29].
This paper is organized as follows: In Section 2, we propose a new multiphase model for the effective dielectric constant; the damping factor and the dielectric constant of the bound water are discussed.In Section 3, we apply the new model and compare the results with experimental data.Finally, a discussion and the conclusion of this study are given in Sections 4 and 5.

Single-Phase Dielectric Mixing Model
For land surface remote sensing, two mixing theories have been proposed, namely the 'dielectric average' originally proposed by Brown [30]: and the 'refractive average' originally proposed by Birchak [31]: Both approaches relate the effective dielectric constant ε eff of a material to the dielectric constants ε j of its m different components weighted according to their volumetric fractions v j in the mixture.The refractive index n eff measured by TDR or GPR is related to the time duration of electromagnetic wave propagation in a medium, in our case in the soil-water mixture [32], according to: where c is the speed of light, t is the traveling time along the probe rod and L is the length of the probe rod.A further refinement proposed by Ansoult et al. [33] and by Dobson et al. [34] is the empirical modification of Equation ( 1) with an exponent α, the so-called 'shape factor', according to: The linear relationship of the soil moisture with the refractive index, similar to Equation (2), has been widely used in calibration models [35][36][37][38].This approach has further evolved into various power-law-based models: because the soil moisture estimation following Equation (2) could not meet the required accuracy for different soil textures and frequencies of electromagnetic waves, various values have been proposed for α [33,34,[39][40][41][42][43][44][45][46].Furthermore, a nonlinear relationship between water content and the dielectric constant has been proposed for empirical calibration models [47][48][49].Figure 2 gives an overview of these models together with the different values used for α.In Appendix A, we discuss that the effective refractive index n eff should preferably be calculated from the effective dielectric constant ε eff based on the arithmetic mean (parallel mixing) according to Equation (1), and not (2) or (4).

Multi-Phase Dielectric Mixing Model
The volumetric fractions v j in Equation (1) can be described by the soil porosity p and volumetric fraction of water, w, as proposed by Wang and Schmugge [52], ε eff = (1 − p)ε soil + wε water + (p − w)ε air (5) Because the dielectric constant follows the superposition rule of polarizability (see Appendix A), this constant, and not the refractive index, is linearly proportional to the volume of the physical material.Therefore, the arithmetic average of the dielectric constant by volumes of soil water, of soil particles and of air as shown in Equation ( 5) is the preferred approach from a theoretical point of view.
The volumetric ratios of the main soil minerals v sand , v silt and v clay enable us to derive a sub-phase model for ε soil according to: The effective dielectric constant of water in the land surface is the arithmetic average of the dielectric constant of free and of bound water according to: In contrast to the mineral component in the soil, the computation of the dielectric constant of soil water ε water is more challenging because v bound and v free are unknown and cannot be measured; these parameters can only be determined by comparisons of the measured refractive indices with simulated data.
Both the real and the imaginary parts of the complex dielectric constant now change in such a way with the soil water fraction w.For low values of w, when the dry soil particles become moist and w remains below the wilting point, the volumetric soil water content consists only of bound water, and thus, v bound = 1 in Equation ( 9) and ε water = ε bound in Equation (8).Then, Equation (5) becomes: When w is larger than the wilting point w wp , but remains smaller than the porosity p, the dielectric constant of water in the soil is a composite of the dielectric constants of bound water and free water, both of which show very different values for the dielectric constant [55][56][57].The dominant force of bound water is the charge of soil particles, which attracts water molecules and causes them to adhere to the surface of soil particles forming thin layers around the mineral particles (see also Figure 3).The presence of increasing amounts of water in the soil serves to thicken the thin layer of water molecules, which causes the van der Waals force between these molecules [58] to dominate due to the surface charges of the soil particles; as a result, the water is able to move freely.This transition point, the wilting point w wp , depends on the size and the characteristic electric charge of the soil particle.Based on this phase transition, we can model v bound assuming a linear decrease with the total water increase for w wp < w < p.This assumption of a linear relation can be expressed as: When the total soil water w increases from w wp to p, v bound decreases from 1 and 0. We thus get from the assumption of a linear relation: b = − 1 p − w wp (12) and: a = p p − w wp (13) By substituting Equations ( 12) and ( 13) into Equation ( 11), we can write v bound in relation to w wp as: (14) and then get according to the conservation law of Equation ( 9): Thus, we can calculate the effective dielectric constant ε eff in dependence of w by replacing the unknown v bound and v free in the equation for ε water with the known constants w wp , p, ε free .In summary, we obtain as relation between ε eff and w: Next, we investigate the case w > p.After heavy precipitation, not all water infiltrates the soil medium quickly.In this case, we must consider standing water on top of the soil within the observed volume.With increasing w, the fraction of standing water over the soil layer increases, and the fraction of saturated soil medium (free water and soil particles) decreases.We can express the measurable effective dielectric constant with: Similar to the derivation of Equation ( 14), the standing water fraction v standing can be computed using the linear approximation proposed in Equation ( 11) between v standing and w.Then, we obtain: and: We see from Equation ( 18) that if w = p, then v standing = 0.If w = 1, then v standing = 1.Equation (19) shows that if w = p, then v saturated = 1; if w = 1, then v saturated = 0. Equation ( 17) is thus equivalent to: It is interesting to note that the effective dielectric constant of wet soil with w > p can simply be calculated from ε soil , ε free and the measured w regardless of the wilting point and porosity.
The distinct dielectric properties of bound water and free water cause the effective dielectric constant of wet soil (Figure 3b) to become a nonlinear function of w with three domains according to Equations ( 10), ( 16) and (20), the boundaries of which are given by w wp and p, as shown in Figure 4.   10), ( 16) and (20), where w is the wilting point and p is the porosity, both of which are specific parameters for a given soil texture.The gray solid and dotted curves represent the free and bound water contributions, w v ε and w v ε to ε , respectively.10), ( 16) and (20), where w is the wilting point and p is the porosity, both of which are specific parameters for a given soil texture.The gray solid and dotted curves represent the free and bound water contributions, w v ε and w v ε to ε , respectively.10), ( 16) and (20), where w wp is the wilting point and p is the porosity, both of which are specific parameters for a given soil texture.
The gray solid and dotted curves represent the free and bound water contributions, w v free ε free and w v bound ε bound to ε eff , respectively.
For the computation of effective conductivity, it is popular to apply the Bruggeman model [6,59] and Maxwell Garnett model [5,60] similarly to the calculation of the effective dielectric constant.However, in several studies, the computation of the effective conductivity is based on the arithmetic average of the volumetric mixing ratio [2,61,62].Similar to Equations ( 10), ( 16) and (20), we can develop a new multiphase mixing model also for the effective conductivity and get: for w ≤ w wp ( 21) and: The further derivation of σ eff , which consists of soil σ soil , bound water σ bound , free water σ free , salinity σ salt and air σ air , is provided in Appendix B.
The complex dielectric constant ε * consists of the effective polarization term ε eff as the real part, the effective dielectric loss ε eff and effective conductivity σ eff according to: For w wp < w < p and p < w, the effective dielectric constants for both the real part ε eff and imaginary part ε eff are computed by Equations ( 10), ( 16) and (20).

Frequency-Dependent Dielectric Relaxation Model
The frequency-dependent dielectric constant for free and bound water can be derived by Debye's relaxation formula, substituting: τ free = (1.1109×10 −10 + 3.824 × 10 −12 T + 6.938 × 10 −14 T 2 − 5.096 into Equations ( 10), ( 16) and (20).ε max free is the static dielectric constant of free water computed according to Klein and Swift [63], and ε min free is obtained according to Lane and Saxton [64].τ free is the dielectric relaxation time for free water given in Equation (30) based on [65].T and S used in these equations are soil temperature ( • C) and salinity ( ).The bound water relaxation is frequency-dependent (omega in Equations ( 25) and ( 26)).It also changes according to soil texture [66,67], the variation of which is reflected in the bound water relaxation model ( 25) and ( 26) by applying a linear equation varying with v clay to the estimation of the parameters ε max bound as shown in (32): with the relaxation time of the bound water (33): The model for the effective dielectric constant for the bound water ( 25) and (26) does depend on wavelength and soil texture, and the optimal fitting parameters, a and b, will be provided in Section 3.

Bulk Dielectric Mixing Model
For the dielectric constant of a three-dimensional medium, its geometrical structure and thickness are relevant.For example, the effective value of the polarization density P should be reduced by 1/3, a number termed the depolarization factor for the spherical shape [3].The non-spherical shape effect of soil particles in saturated and unsaturated porous media was reviewed by [68].Recent studies showed that by increasing the thickness of a medium, the effective value of the dielectric constant decreases [44,[69][70][71].In short, a reduction of the bulk dielectric constant compared to the flat 2D dielectric constant is attributable to both the 3D microgeometry and the thickness of the medium.In the following, we idealize the complex morphology of the land surface properties found in nature and consider the soil as isotropic.Thus, we focus on the damping effect within the sampling depth, such as the penetration depth of TDR or GPR or an emission depth of SMOS [72] or SMAP [73].
The sampling depth describes that part of the medium that is responsible for the surface emission; thus, the dielectric constant of this layer is relevant [74] because only the energy emitted from it is detected by remote sensing instruments [75].In other words, the signal detected over a surface is not only related to the arithmetic average of the dielectric constant of the soil properties, but also determined by the sampling depth.With increasing depth, the contribution of the signal from that depth to the average signal intensity becomes smaller and smaller.Without considering attenuation, the quantification of soil moisture content from the surface emission alone is underestimated.
The computation of the surface emission from the sampling depth is related to a nonlinear weighting function, such as a decaying function with depth [76].Wilheit [77] proposed a layered model for the electric field in the distinct interface between air (free space) and the soil surface.For each layer, the surface polarization density is related to the surface effective susceptibility χ eff and the electric field E(z) in that layer according to: If the medium is homogenous, Beer-Lambert's law yields: (35) with the effective penetration depth δ eff , which depends on the imaginary part of the effective refractive index.By integrating the surface polarization density P eff (z) up to the penetration depth z, we obtain: The same can also be written with the effective electric susceptibility of a bulk substance χ 3D eff as: Therefore, we get: where: By inserting the attenuation factor H z into Equation ( 24), we can obtain the real part of the 3D complex effective dielectric constant according to: and for the imaginary part: Furthermore, we get for the effective conductivity: If the sampling of the soil layer is obtained with penetration depth δ eff , one needs to insert δ eff for z in Equation (39), and H z becomes the constant H δ (independent of wavelength and extinction coefficient) As discussed later, we found that this simple approach of a constant damping factor for the three-dimensional nature of the soil provides excellent results compared to experimental data.

Comparison with Other Approaches
In the following, we compare the general expression of the dielectric mixing and refractive mixing model: ε eff = aε soil + bε bound + cε free + dε air (44) n eff = an soil + bn bound + cn free + dn air (46) where a, b, c and d are the volumetric mixing ratios of each phase, which is a function of soil properties and empirical parameters, as listed in Table 1.To ensure correspondence between the simulation of the effective dielectric constant and the actual measurements, Wang and Schmugge [53] proposed the transition moisture w t , which is larger than the wilting point.If w t is utilized to compute the bound water fraction b and free water fraction c according to Equations ( 14) and ( 15), b increases and c decreases, thereby causing an underestimation of the effective dielectric constant.Therefore, the use of w t requires the empirical parameter Υ to obtain a fit between the predicted and measured effective dielectric constants.Dobson et al. [34] suggested for α in Equation (45) an empirical value of 0.65 and w β with β, which is a function of soil texture for c.This choice leads to a total volumetric ratio larger than 1 in some cases, which is of course unrealistic.In the approach of Mironov et al. [53,78], a negative volumetric ratio appears in d, the volumetric mixing ration of air, in some cases.Moreover, the free water dielectric constant reaches quite high values of up to 100 and becomes furthermore an empirical function of the clay-mixing ratio, which seems unrealistic.In addition, this approach does not use the porosity information in the refractive mixing model (see Equation ( 46)) and, thus, does not show a connection between the soil texture and the effective dielectric constant, which seems unrealistic, as well.Therefore, while previous models simulate measured data experimentally with satisfactory accuracy, these models have to assume unphysical parameters in order to achieve this, which is unsatisfactory from a theoretical point of view.We also investigated a simple calibration model based on the linear relationship between the refractive index √ ε eff and the soil water content w : which is widely utilized for GPR applications, where a = 0.1168, b = −0.19 in [36], and for TDR applications, where a = 0.1138, b = −0.1758 in [35], a = 0.1181, b = −0.1841 in [37], and a = 0.14, b = −0.2 in [38].This approach has been suggested as the refractive-index mixing approach, which only considers the soil water content.We can express the effective dielectric constant ε eff as a polynomial function of the soil water content w according to: Equation ( 48) needs various empirical fitting parameters for different soil textures.This quadratic polynomial calibration model was also proposed for microwave remote sensing by Hallikainen [50] with: Compared to Equation (48), Equation ( 49) is not only a function of the soil water content w, but also of the soil texture.Therefore, the model is able to simulate the effective dielectric constant of various soil textures with good agreement to the measurements; however, it requires quite a large number of nine empirical fitting parameters.
One can see from in Equation ( 50), which is actually the same formula of Equation ( 16), that the fitting parameters of the empirical polynomial models of Equations ( 48) and ( 49) are a combination of the dielectric constants ε bound and ε free and the physical soil properties w wp and p.As a result, one can estimate the soil water content from the dielectric measurements without the need for empirical parameters or different calibrations for different soil conditions; the only unknown parameter, which cannot be resolved, is ε bound .However, this single parameter can be determined by comparing the mixing model with observations.In this study, we compare our model with above models including a modified Mironov model.The model of Mironov et al. (2009) is the originally proposed refractive mixing model [53,78], which uses a static dielectric constant of free water (ε max free ) of 100 in Equations ( 25) and ( 26).In contrast to this, " Mironov et al. (2009), modified" is the model of Mironov et al. (2009) modified with ε max free = 80, which is the same value as used in all other approaches.This value can either be measured [79] or modeled [63] at 20 • C with the surface pressure level of 0.1 MPa (see also Equations ( 25) and ( 26)).The performance scores of these models are evaluated with: and:

Summary of New Model
In summary, we can calculate the wave-dependent complex effective dielectric constant for various soil textures using the following formulas: for w ≤ w wp : for w wp < w ≤ p: and for p < w:

Experimental Reference Data
In the prediction of the dielectric constant, the model proposed in this study requires the direct gravimetric measurement of soil moisture, temperature, volumetric mixing ratio of clay, silt, sand and salinity measurements.The fact that the organic material volume information is not taken into account is likely to be an error factor in our model.Unlike soil samples (A-I in Table 2), soil water contents measured in radio waves (J-X in Table 2) do not contain salinity information.Although it is obviously an error factor in the imaginary part of the dielectric constant, the verification was performed with the dielectric measurements in the range of the radio wave band because it is not a significant error factor in the computation of the real parts of the dielectric constant.
The verification of the dielectric constants predicted by the model with the observed soil parameters as an input is carried out based on the dielectric constant measured by microwave probes or TDR.We used the Roth et al. data [48] as shown in Table 2 on the assumption that instrumental error was minimized.The total number of samples was verified with a total of 131 point measurements (69 from L-, 41 from C-, 21 from K-bands) at 9 sample sites in the microwave and 148 point measurements at 15 sample sites in 30-MHz measurements.Table 2 shows the values of the volumetric mixing ratio of soil texture for v sand , v silt and v clay , salinity and temperature.In order to find proper values for the wilting point w wp and porosity p, we first identify the soil texture by matching the measured sand, silt and clay mixing ratios of the references to the USDA (United States Department of Agriculture) soil texture classification.According to the identified soil texture, we obtained w wp and p from the input table for the soil parameters in NOAH-MP (Noah land surface model with multi-parameterization options) (Niu et al. 2011b), which are based on the STATSGO (State Soil Geographic) Data Base (Soil Survey Staff).Table 3 lists the values for w wp and p following this approach for different soil textures.ε bound , ε free ; ε bound and ε free are calculated with Equations ( 25) and (26).σ bound and σ free are calculated with Equations (A13) and (A14), respectively, using Table 4; σ 25 • C and ϕ are calculated with Equations (A10) and (A11), respectively, using the salinity and temperature given in Table 2; σ soil follows Equation (A5) with the data of Table 4 and σ air = 0.  [80] are obtained by weight.Because the specific gravity of soil minerals is almost the same [80], their weight percentages are equivalently utilized as the volumetric mixing ratio and the values of salinity notated with * were assumed as there was no information in the references about which value fits best to the experimental data.Table 4. Electrical conductivity of different soil textures [82].
In order to find the fitting parameters a and b of Equation ( 32), the optimum ε max bound reaching the least RMSE (root-mean-square error) compared to the measured effective dielectric constant are determined from the wet soil A-I according to various authors from Table 2.The optimum value is the best combination of values for unknown bound water parameters ε max bound , assuming it is probably closest to the true soil condition.The ε max bound is determined when the RMSE becomes the minimum in the given clay contents as follows.
Regardless of the availability of bound water information, we are able to complete the prediction of the effective soil dielectric constant based on the measurements of volumetric clay contents owing to the proposed bound water model, Equations ( 25) and (26).Finally, our model can be tested based on soil texture, temperature and salinity information with the relevant formulas, as summarized in Table 5.

Results
The performance of the new approach was assessed by comparison with experimental data, as well as with previous dielectric mixing models.The Wang and Schmugge model has a smaller dielectric constant for bound water compared to other approaches (see Table 6).The model of Dobson et al. does not deal separately with bound and free water, but integrates them into a dielectric constant of water.For the comparison with Mironov et al. [53], we simulated the effective dielectric constant with Mironov et al. (2009) and Mironov et al. (2009), modified, which are defined at the end of Section 2.5.Because the purpose is the operational use of microwave remote sensing for soil moisture estimation on a global scale, the consistent performance of the dielectric prediction with high mean accuracy is decisive.Therefore, we validated the new approach with experiments in the L-band (Figure 5), C-band (Figure 6), K-band (Figure 7) and their RMSE analysis in Figure 8. Furthermore, we performed the prediction of the effective dielectric constant in a radiowave band at 30 MHz (Figure 9) over various soil textures and provided the RMSE results in Figure 10.The average values are also listed in Table 7. Dobson's model [34] is based on the evaluation and calibration of five soil texture values from sandy loam to silt clay for which the sand mixing ratio ranged from 5-51%.The simulation of pure sand with this model shows a too high dielectric constant in the real part (see Figure 5a), e.g., 27 instead of 18 for w = 0.275.Other approaches, such as that of Wang and Schmugge [52] and that of Hallikainen et al. [49], generally underestimate the dielectric constant in such a situation and, therefore, result in too large values of soil moisture from the dielectric measurements.The values obtained with our model are in good agreement (e.g., only a 3% underestimation for w = 0.275) with the dielectric constant measurements of wet sand, also capturing the measurements above the porosity point.The models of Mironov et al. (2009) [53] and Hallikainen et al. [49] overestimate the imaginary part of the dielectric constant (e.g., 67% and 144% overestimated in w = 0.275, respectively (see Figure 5a), whereas other models including ours are closer to the measurements (e.g., an 18% overestimation for w = 0.275).When we compare the values obtained from the various models along the gradient of increasing clay content (see Figure 5a-i), the dielectric predictions for the real part are similar; in general, all of the predictions are mostly in reasonable agreement with the measured values of the dielectric constants.However, the models show significantly different performances in the prediction of the imaginary part.Considering the modeled curve for the imaginary part shown in Figure 5b,g, in the higher clay contents (Figure 5g), the predictions of various models diverge.The new model is in best agreement for both silt and silty clay loam with absolute RMSE values of only 0.34 and 0.48, respectively.For the clay case (Figure 5i), the new model shows an overestimation.One of the possible reasons might be a swelling process of heavy clay in the oversaturated soil condition.A significantly larger porosity of the sampled clay soil than the value that we applied in our model (0.500, see Table 4) might explain this disagreement.1 and 2) from our approach, Wang and Schmugge, Dobson et al., Hallikainen et      (a-c) Same as Figure 5, but for 18 GHz (K-band) with measured data from [34] and [66].
Our new approach is more accurate even than the previous semi-empirical and empirical approaches (see the averaged analysis in Table 7 and Figure 8).(a-c) Same as Figure 5, but for 18 GHz (K-band) with measured data from [34] and [66].

Discussion
The range of w > p is not covered in the existing models.Using the newly-proposed model, it is possible to improve the accuracy of the flood monitoring using the satellite and of the soil moisture estimation by the soil moisture probes even in the challenging circumstance such as the swelling or ponding condition.For example, the average RMSE of the new model in the samples A, B, D, G, J, K in Figures 5 and 9, considered that the case w > p is included in those samples, is 36%, 77% and 42% lower than Mironov et al. (2009) in the prediction of the real and imaginary parts in the microwave band and of the real part in the radio wave band, respectively.However, the present version has been ideally designed according to twelve soil texture classification in which the wilting point and porosity of soil medium are also restricted to only twelve classes.In a real soil medium, because of its heterogeneous nature, not all soil particles reach the wilting point and porosity identically in the bulk of soil medium.For example, if the total water is measured as much as dry porosity, some of the portion can still be occupied with bound water defined in the formula for w < p, and others already can exceed the porosity point expanding water portion w > p.However, as is stated, the total amount of the soil water will be p.If this heterogeneous nature of the soil medium is successfully reflected by the error function, a smoothing effect will appear where the discontinuity has been shown in the current model.However, the approach based on this strong assumption requires a verification based on the purposefully-designed experimental observations.Through this kind of effort, the model proposed in this study can be improved as a continuous function mathematically differentiable at all points, which is an important requirement in the minimization process in the soil moisture retrieval algorithm of SMAP, SMOS and AMSR-E.

Conclusions
The physically-correct averaging method for the effective dielectric constant is the arithmetic mean because it follows the superposition rule of polarizability.We propose a new model for the effective dielectric constant of bare soil, which considers the soil water phase as consisting of free and bound water and the bulk behavior by including a damping factor.The resulting effective dielectric constant of bare soil shows a complex nonlinear behavior as a function of the volumetric fraction of water.In addition, our model allows for the calculation of the effective dielectric constant in the oversaturated soil condition.This should facilitate investigating surface runoff and the infiltration capacity of the surface and, thus, the tracking and forecast of flooding events.

Figure 1 .
Figure 1.Measuring principles of bare soil properties (black boxes) by satellite and airborne remote sensing in the microwave region (A) and TDR (time-domain reflectometer) and GPR (ground-penetrating radar) (B).The dielectric constant (hexagon) is the link between the measured parameters of the sensors (brightness temperature (TB) and refractive index, respectively) and the targeted soil properties.

Figure 3 .
Figure3.Illustration of the rescaling process for the volumetric fractions in moist soil.(a) Total soil water content w is smaller than the wilting point wwp.Here, only bound water is present; (b) w greater than wwp, but smaller than the soil porosity p. Part of the water is bound and part is free; (c) w greater than p.There is free water inside the soil and standing water above the soil, but no bound water in our mixing model.The dielectric constant of free and standing water is the same.

Figure 4 .
Figure 4. Multiphase behavior of the effective dielectric constant ε for different values of the soil water fraction w according to the new approach of Equations (10), (16) and(20), where w is the

Figure 3 . 30 Figure 3 .
Figure3.Illustration of the rescaling process for the volumetric fractions in moist soil.(a) Total soil water content w is smaller than the wilting point w wp .Here, only bound water is present; (b) w greater than w wp , but smaller than the soil porosity p. Part of the water is bound and part is free; (c) w greater than p.There is free water inside the soil and standing water above the soil, but no bound water in our mixing model.The dielectric constant of free and standing water is the same.

Figure 4 .
Figure 4. Multiphase behavior of the effective dielectric constant ε for different values of the soil water fraction w according to the new approach of Equations (10), (16) and(20), where w is the

Figure 4 .
Figure 4. Multiphase behavior of the effective dielectric constant ε eff for different values of the soil water fraction w according to the new approach of Equations (10), (16) and(20), where w wp is the wilting point and p is the porosity, both of which are specific parameters for a given soil texture.The gray solid and dotted curves represent the free and bound water contributions, w v free ε free and w v bound ε bound to ε eff , respectively.

Figure 5 .
Figure 5. (a-i) Effective dielectric constant at 1.4 GHz (L-band) for nine soil samples from sand to clay (see Tables1 and 2) from our approach, Wang and Schmugge, Dobson et al.,Hallikainen et  al. and Mironov et al. in comparison to measured data from [34,52,80,81]: vertical lines indicate wwp and p.

Figure 5 .
Figure 5. (a-i) Effective dielectric constant at 1.4 GHz (L-band) for nine soil samples from sand to clay (see Tables1 and 2) from our approach, Wang and Schmugge, Dobson et al.,Hallikainen et  al. and Mironov et al. in comparison to measured data from [34,52,80,81]: vertical lines indicate w wp and p.

Figure 5 .
Figure 5. (a-i) Effective dielectric constant at 1.4 GHz (L-band) for nine soil samples from sand to clay (see Tables1 and 2) from our approach, Wang and Schmugge, Dobson et al.,Hallikainen et  al. and Mironov et al. in comparison to measured data from [34,52,80,81]: vertical lines indicate wwp and p.

Figure 8 .Figure 8 .
Figure 8. Absolute RMSE for the real part (a) and imaginary part (b) and relative RMSE (in %) for the real part (c) and imaginary part (d) of the effective dielectric constant for different soil textures in the L-, C-and K-band for the models H [49], D [34], M1 [53], M2 (Mironov et al. (2009), modified), W [52] and P (this work).The last values show the overall RMSE.

Table 1 .
(44) of volumetric mixing ratios in the general expression of the dielectric and refractive mixing formulas.The following equations were used:(44)for Wang and Schmugge and this paper (Park et al.), (45) for Dobson et al., and (46) for Mironov et al.

Table 2 .
Volumetric mixing ratio of soil texture, salinity and temperature from the references.

Table 3 .
Wilting point and porosity according to soil texture from the NOAH land surface model with multi-parameterization options (NOAH-MP).
[53]ing point and porosity from the State Soil Geographic (STATSGO) data (cm 3 /cm 3 )[53]applied to the NOAH-MP land surface model; * increased witling point and porosity to match the measurement points.

Table 5 .
Required formulas and data to compute the complex effective dielectric constant Equations (53)-(58).

Table 6 .
Complex dielectric constant for land surface properties.The range of bound water dielectric constant computed from sand to clay; * same value as Mironov et al. (2009); value for porous Soil A and for solid Soil B.
al. and Mironov et al. in comparison to measured data from [34,52,80,81]: vertical lines indicate w wp and p.

Table 8 .
Absolute RMSE over different soil textures in the radiowave band (30 MHz) for the models D (Dobson et al. 1985), M1 (Mironov et al. 2009), W (Wang and Schmugge 1980) and P (the proposed model in this study).Gray boxes mark the best performance.