relationship of Moisture content and its dielectric constant

0 Introduction

http://www.soilmanagementindia.com/soil-water/how-to-measure-soil-moisture-soil-water-relationship-soil-science/4514

Drying the moist soil to a constant weight in a drying oven is controlled at 110° ± 5°C. Samples should be dried for at least 24 h. Weight of soil remaining after oven drying is used as the weight of soil solids. Moisture content expressed as a per cent is equal to the weight of water divided by weight of soil solids all times 100.

Volumetric water content (%) can be calculated if soil bulk density (BD) is known.

Soil moisture content can be expressed in depth of water (cm) per unit depth of soil.

Depth of water per unit soil depth = Volumetric water content x Soil depth.

1 Topp's model for soils

Ref:http://www.decagon.com.br/micro/dr-topp/

The equation specifically is: 

Θ (m3/m3) = 4.3 X 10-6 * ε3 - 5.5 X 10-4 * ε2 + 2.92 X 10-2 * ε -5.3 X 10-2

Topp and his team showed that the volumetric water content of soil can be determined from the apparent dielectric constant of the soils, independent of soil type, soil temperature, and soluble salt content. 

https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2003WR002589

ksi=1;i=0;a=[];b=[];moiscontent=[];
for ksi=10:5:40
    i=i+1;
%     disp('dieletric constant')
%     ksi
a(i)=ksi;
    pho=1.8;% g/cm3 dry weight
%moiscontent(i)=1*(4.3*10^(-6)*ksi^3-5.5^10^(-4)*ksi^2+2.92*10^(-2)*ksi-0.053); %by volume
    moiscontent(i)=(sqrt(ksi)-0.819-0.6187*pho-0.159*pho^2)/(7.17+1.18*pho);
   % disp('moisture content%')
   b(i)= 100* moiscontent*1/((1-moiscontent)*pho+moiscontent);% water content 
   
end
figure
yyaxis left
plot (a,moiscontent,'b-')
ylabel('volumetric moisture content')
hold on
yyaxis right
plot(a,b,'r')
ylabel('soil moisture')
xlabel('dietricity constant')

Ref:http://quebec.hwr.arizona.edu/classes/hwr432/2014/roth92-soil-dielectric-constant-vs-water-content.pdf

2 For pavement

ksi=1
for ksi=10:1:16
    ksi
   % moiscontent=4.3*10^(-6)*ksi^3-5.5^10^(-4)*ksi^2+2.92*10^(-2)*ksi-0.053;%by volume
    moiscontent=-6.216+2.383*ksi-0.0598*ksi^2+0.0006*ksi^3
   
   %
    disp('moisture content%')
   100* moiscontent*1/((100-moiscontent-3)*2.6+moiscontent)
end

density of pavement=2.6g/cm3

Result

Compared with Dr. Muller's result. They are almost the same.