Dielectric Properties of Glass Beads with Talc as a Reference Material for Calibration and Verification of Dielectric Methods and Devices for Measuring Soil Moisture.

This paper presents dielectric measurements of talc, glass beads, and their mixtures under different moisture and salinity levels. The measurements were conducted using a prototype seven-rod probe (15 mm long central rod) connected to a single port of vector network analyzer. The samples were moistened with distilled water and KCl solutions in order to obtain six different moisture content levels. The complex dielectric permittivity was determined from vector network analyzer reflection-coefficient measurements based on the open-water-liquid calibration procedure. Next, the fitting of volumetric water content–real part of dielectric permittivity calibration curves was performed for each material at selected frequencies, and the obtained relations were compared with well-known calibration equations. Additionally, a salinity index for the tested materials was calculated. It was concluded that pure talc is not an optimal material for the calibration and verification of dielectric methods. The calibration curves obtained for glass beads and the mixtures of glass beads with talc gave results close to well-known reference calibration functions. Additionally, the addition of talc caused the data points to be less scattered. Moreover, the values of the salinity index for the tested materials were in a good agreement with literature data for sand. The obtained results indicated that glass beads with the addition of talc can be used as a reference material for the calibration and verification of dielectric methods and devices for soil moisture measurement.


Introduction
The monitoring of soil water content provides essential information that is necessary to mitigate the impact of natural disasters like floods or droughts, and it is also a valuable tool for a precise irrigation, which, in turn, can bring reductions in water consumption [1,2]. There are several methods used to measure soil moisture that can be divided into in-situ methods and remote methods. Among the where the real (ε ) and imaginary (ε ) parts of dielectric permittivity depend on frequency ( f ) of the applied electromagnetic wave. It is recognised that ε in soil mainly depends on volumetric water content (θ), and, to a lesser degree, on other soil parameters such as mineralogical properties and particle-size distribution. The imaginary part ε depends on the dielectric loss (ε d ) and on the bulk electrical conductivity (σ b ) according to the following formula: where ε 0 ≈ 8.85 ·10 −12 F m −1 is vacuum dielectric permittivity. In the literature, there are a number of theoretical and empirical models describing the relation among the θ, σ b , and electrical conductivity of soil solution. One of them is the salinity index model developed in [6], where SI is defined as: where ε a is the apparent dielectric permittivity of soil. The application of dielectric techniques in order to estimate soil moisture content has resulted in the development of an increasing number of sensors working at various frequencies from kHz up to several GHz [4]. TDR sensors, operating in a frequency band of 1-2 GHz, measure dielectric permittivity with a high accuracy. However, they are still too expensive to be commonly used. The FDR and capacitance methods represent an alternative to the TDR technique because of their lower cost. Nevertheless, these sensors usually operate at a single frequency in the range of 20-300 MHz. At these frequencies, soil complex permittivity can exhibit several polarization phenomena, including bound water, double-layer, and Maxwell-Wagner effects [7,8]. These phenomena depend not only on water content and salinity but also on other properties of soils such as mineralogy and particle-size distribution.
From a user point of view, a dielectric sensor should be able to operate in different types of soil and should be capable of adapting to changing environmental conditions. The manufacturers usually offer commercial sensors with the calibration equation that relates θ to either the dielectric permittivity or the output magnitude unit of the sensor. However, in many cases, a soil-specific calibration is recommended, even though it is challenging for many users to perform. In general, dielectric sensors can be calibrated in two ways. One way, called one-step calibration, is to directly calibrate the sensor against soil water content like in the gravimetric method of soil moisture determination [9]. However, applying a direct calibration is often too labour-intensive and time-consuming. The use of a two-step calibration seems to be a more promising method [10]. The first step of this calibration relies on measurements of reference (calibration) materials with known values of ε * . The second step is the application of the relations between θ and ε * in the form of an empirical [11] or semi-theoretical model. Materials 2020, 13, 1968 3 of 12 This approach allows for the determination of θ. However, the results obtained from different dielectric sensors may vary.
The electromagnetic sensors are usually evaluated using the determination of permittivity in the soil over a range of water contents [12]. These studies are useful for the demonstration of the general water content measurement capability of specific soils. The obtained results are often misleading, and the effects arising from bound water or salinity can be disguised by used soil-specific calibration. The response of the sensors is varied among soil type, temperature, and salinity [13]. The reliable comparison of the performance of dielectric sensors requires the use of repeatable and stable reference materials with well-known dielectric properties, reflecting the dielectric spectrum of soil. The "ideal" reference materials should be non-toxic, inexpensive, and easily accessible. Due to the fact that many soils (e.g., especially fine-textured soils) exhibit multiple relaxation phenomena including Maxwell-Wagner effects at low frequencies, the reference materials should also reproduce dielectric relaxation and ionic conductivity behaviour [10]. Several authors have developed sensors calibration techniques by using calibration media such as air and liquids like water, 2-isopropoxyethanol, dioxane, ethylene glycol, methanol, and ethanol [10,[14][15][16]. The use of liquids as calibration media has several advantages: (i) no air gaps and density variations, (ii) the ability to calibrate multiple sensors in a wide range of permittivity (from 2 to 35), and (iii) the possibility to separate sensor-and soil-specific effects (air or water gaps, metal electrode polarization, and dielectric dispersion at a specific frequency) [17]. Nevertheless, most of the calibration liquids, such as dioxane and 2-isopropoxyethanol, are harmful to health. Additionally, liquids do not represent all of the loss mechanisms and relaxations that occur in the soil in the range of the minimum to the maximum moisture level. An alternative is the use of grainy solids or powdery materials that reflect the dielectric permittivity spectrum of a natural soil, e.g., minerals (kaolinite, bentonite), glass beads, talc, sand, and calcium carbonate under different moisture conditions [14,18]. However, the knowledge of the dielectric complex permittivity spectrum of these materials under different moisture conditions, including different salinity levels, is required. The application of glass beads and talc as reference materials has several advantages: (i) they regular spherical particles with a known and repeatable diameter, (ii) they have easy access, (iii) they allow for their easy recovery from a mixture, and (iv) a mixture of these materials with water has low hazard potential and there are no effects on the sensor [19]. Additionally, it is worth noting that talc is characterized by low dielectric losses [20].
The aim of the study was to measure the complex dielectric permittivity of talc, glass beads, and their mixtures under different moisture and salinity levels. This paper presents the following:

•
Electromagnetic simulations of the seven-rod probe with different dimensions using the Ansys HFSS software (Version 19.2.0, ANSYS, Inc. Southpointe, Canonsburg, PA, USA).

•
Measurements of the ε and ε values of talc, glass beads, and glass beads with 5% and 10% of talc under different moisture and salinity conditions (electrical conductivities of 0.5, 1.0, and 1.5 S·m −1 ) using the short seven-rod probe connected to a PNA-L N5230C VNA (Agilent Technologies, Santa Clara, CA, USA) in the operational frequency range of the probe.

•
Comparison between ε and θ with well-known calibration equations at a several frequencies.

•
Estimation of repeatability of the materials under test.

•
Evaluation of the salinity index model.

Seven-Rod Probe
In this study, complex dielectric permittivity spectra were measured using a short seven-rod probe with 15 and 20 mm long inner and outer rods, respectively, with a measurement volume of about 2800 mm 3 . The probe was manufactured at the Laboratory of Dielectric Spectroscopy of the Institute of Agrophysics PAS, Lublin, Poland (Figure 1a,b). The probe used in this work (a short seven-rod probe) was a modified version of the probe described in [21]. The construction of the seven-rod probe Materials 2020, 13,1968 4 of 12 enabled us to measure soil moisture in a well-defined volume. Because the operating frequency range of the probe shown in [21] was 200 MHz, in order to extend its operational frequency range, another construction of probes with shorter central and outer rods was proposed. Probes with 25 and 30 mm long central and outer rods and 15 mm and 20 mm long central and outer rods, respectively, were tested by the numerical simulations using the finite element method (FEM) in the Ansys HFSS software [22]. The simulated probes were surrounded by tested material in the form of a cylinder of 70 mm in height and 15 mm in diameter filled with calibration media (air, water, and ethanol) and isopropanol as a verification liquid. More details about the simulation space and the applied model were presented in [21]. The extraction of complex permittivity spectrum from the simulated reflection coefficient (S 11 ) was performed using the following bilinear equation: where c 1 , c 2 , and c 3 are calibration constants determined in the open-water-liquid (OWL) calibration [23,24].
Materials 2020, 13, x FOR PEER REVIEW 4 of 12 range, another construction of probes with shorter central and outer rods was proposed. Probes with 25 and 30 mm long central and outer rods and 15 mm and 20 mm long central and outer rods, respectively, were tested by the numerical simulations using the finite element method (FEM) in the Ansys HFSS software [22]. The simulated probes were surrounded by tested material in the form of a cylinder of 70 mm in height and 15 mm in diameter filled with calibration media (air, water, and ethanol) and isopropanol as a verification liquid. More details about the simulation space and the applied model were presented in [21]. The extraction of complex permittivity spectrum from the simulated reflection coefficient ( ) was performed using the following bilinear equation: where , , and are calibration constants determined in the open-water-liquid (OWL) calibration [23,24]. The value of the of isopropanol (verification liquid) obtained from Equation (4) was based on the simulation data. Next, the obtained of isopropanol obtained from the two probes with different dimensions were compared with the data shown in [21] for a probe with 40 mm long inner rod and with the data of isopropanol found in the literature [25].
The measurements of dielectric properties of the mixture of glass beads with 5% of talc with distilled water were performed in three independent repetitions. Additionally, six samples with different moisture contents were obtained for talc. The measurements were conducted using the short sevenrod probe (Figure 1) connected to a PNA-L N5230C VNA (Agilent Technologies) over the frequency range from 1 to 500 MHz ( Figure 2) at 20 °C. At those frequencies, an OWL calibration and the complex dielectric permittivity extraction procedure described above were performed using methanol as the calibration liquid. The weight of both air-dry glass beads and talc did not change after drying at 105 °C for 24 h. This suggested that these materials do not absorb water. A mass of 60 g material (talc, glass beads or glass beads with talc) was mixed with a designated amount (in gram) of liquids in order to obtain 6 different values of moisture contents from 2% to 3% of volumetric water content to near saturation. Next, the sample was transferred to a vial (50 mL) where the volume of each measured sample was about 30800 mm 3 . After this experiment, all samples were dried in the laboratory oven in 105 °C in order to obtain the dry mass of the samples [26]. The volumetric water content was calculated for all samples. In the next step, the -calibration curves in the form of linear functions of dielectric permittivity were fitted to the obtained data. Additionally, the salinity index (SI) was calculated. The value of the ε of isopropanol (verification liquid) obtained from Equation (4) was based on the simulation data. Next, the obtained ε of isopropanol obtained from the two probes with different dimensions were compared with the data shown in [21] for a probe with 40 mm long inner rod and with the data of isopropanol found in the literature [25].

Materials Characteristics and Experiment
The glass beads, 90-106 µm (Fuji Manufacturing Industries, Japan), and mixtures of glass beads with 5% and 10% of talc (Sigma Aldrich) were moistened with distilled water (2 × 10 −4 S·m −1 ) and with three potassium chloride (KCl) solutions of electrical conductivities of (σ KCl ) 0.5, 1.0, and 1.5 S·m −1 . The measurements of dielectric properties of the mixture of glass beads with 5% of talc with distilled water were performed in three independent repetitions. Additionally, six samples with different moisture contents were obtained for talc. The measurements were conducted using the short seven-rod probe (Figure 1) connected to a PNA-L N5230C VNA (Agilent Technologies) over the frequency range from 1 to 500 MHz ( Figure 2) at 20 • C. At those frequencies, an OWL calibration and the complex dielectric permittivity extraction procedure described above were performed using methanol as the calibration liquid. The weight of both air-dry glass beads and talc did not change after drying at 105 • C for 24 h. This suggested that these materials do not absorb water. A mass of 60 g material (talc, glass beads or glass beads with talc) was mixed with a designated amount (in gram) of liquids in order to obtain 6 different values of moisture contents from 2% to 3% of volumetric water content to near saturation. Next, the sample was transferred to a vial (50 mL) where the volume of each measured sample was about 30,800 mm 3 . After this experiment, all samples were dried in the laboratory oven in 105 • C in order to obtain the dry mass of the samples [26]. The volumetric water content was calculated for all samples. In the next step, the θ-ε calibration curves in the form of linear functions of dielectric permittivity were fitted to the obtained data. Additionally, the salinity index (SI) was calculated.

Results and Discussion
The shortening the rods of the probe with respect to the solution from [21] enabled us to extend the measurement frequency range. The simulation of the values of isopropanol obtained from the short seven-rod probe with the 15 mm central rod agreed with data obtained from the literature [25] up to 500 MHz (Figure 2), which was also confirmed by the low value of root-mean-square error (0.024) between the simulated value and the known reference permittivity of isopropanol. The similar root mean square error (RMSE) (0.025) was obtained up to 200 MHz for the probe with a 40 mm central rod, as presented in [21]. The sample spectra of selected materials obtained using the short seven-rod probe over the frequency range from 1 to 500 MHz are presented in Figure 3.

Results and Discussion
The shortening the rods of the probe with respect to the solution from [21] enabled us to extend the measurement frequency range. The simulation of the ε values of isopropanol obtained from the short seven-rod probe with the 15 mm central rod agreed with data obtained from the literature [25] up to 500 MHz (Figure 2), which was also confirmed by the low value of root-mean-square error (0.024) between the simulated ε value and the known reference permittivity of isopropanol. The similar root mean square error (RMSE) (0.025) was obtained up to 200 MHz for the probe with a 40 mm central rod, as presented in [21].
The sample ε spectra of selected materials obtained using the short seven-rod probe over the frequency range from 1 to 500 MHz are presented in Figure 3.

Results and Discussion
The shortening the rods of the probe with respect to the solution from [21] enabled us to extend the measurement frequency range. The simulation of the values of isopropanol obtained from the short seven-rod probe with the 15 mm central rod agreed with data obtained from the literature [25] up to 500 MHz (Figure 2), which was also confirmed by the low value of root-mean-square error (0.024) between the simulated value and the known reference permittivity of isopropanol. The similar root mean square error (RMSE) (0.025) was obtained up to 200 MHz for the probe with a 40 mm central rod, as presented in [21]. The sample spectra of selected materials obtained using the short seven-rod probe over the frequency range from 1 to 500 MHz are presented in Figure 3.  Despite the fact that the simulation gave the results up to 500 MHz that were consistent with the literature, in the case of the measurements, the upper working frequency of the probe was limited to 400 MHz because a slight increase of was observed above 400 MHz for several samples, which is not a physical behaviour. Additionally, some scattering of the values was observed at low frequencies. This scattering was more visible in the case of samples with higher water and salt contents. As expected, for each sample, the value of increased with an increasing water content [21,27,28]. Among all the samples, the highest value of for the saturation point was obtained for talc. According to the literature data, talc and glass beads exhibit hydrophobic character. This means that water does not adsorb inside the particles of these materials. However, water can be absorbed on the particles' surface due to the forces cohesion and adhesion [29]. The highest value of for the saturation point of talc can be explained by its smaller particle size as opposed to those of the glass beads.
In the next step, the -calibration curves were fitted in the form of a linear function of the square root of dielectric permittivity, as given by the following equation: The fitting equations were obtained for each of the materials, except for talc, at selected frequencies for data including distilled water and all sodium chloride solutions. For talc, it was only obtained for distilled water. Figure 4 shows the relations between values for the materials under test and an appropriate reference equation available in literature: Figure 3. The ε values of (a) talc, (b) glass beads, glass beads with (c) 5% and (d) 10% of talc moisturized with distilled water, and glass beads with 5% of talc moisturized with KCl solutions of (e) 0.5 and (f) 1.5 S·m −1 conductivities over the frequency range from 1 to 500 MHz (selected data). The legends present the θ (volumetric water content) value of samples.
Despite the fact that the simulation gave the results up to 500 MHz that were consistent with the literature, in the case of the measurements, the upper working frequency of the probe was limited to 400 MHz because a slight increase of ε was observed above 400 MHz for several samples, which is not a physical behaviour. Additionally, some scattering of the ε values was observed at low frequencies.
This scattering was more visible in the case of samples with higher water and salt contents. As expected, for each sample, the value of ε increased with an increasing water content [21,27,28]. Among all the samples, the highest value of ε for the saturation point was obtained for talc. According to the literature data, talc and glass beads exhibit hydrophobic character. This means that water does not adsorb inside the particles of these materials. However, water can be absorbed on the particles' surface due to the forces cohesion and adhesion [29]. The highest value of ε for the saturation point of talc can be explained by its smaller particle size as opposed to those of the glass beads.
In the next step, the θ-calibration curves were fitted in the form of a linear function of the square root of dielectric permittivity, as given by the following equation: The fitting equations were obtained for each of the materials, except for talc, at selected frequencies for data including distilled water and all sodium chloride solutions. For talc, it was only obtained for distilled water. Figure 4 shows the relations between ε values for the materials under test and an appropriate reference equation available in literature: Materials 2020, 13, x FOR PEER REVIEW 6 of 12 Figure 3. The values of (a) talc, (b) glass beads, glass beads with (c) 5% and (d) 10% of talc moisturized with distilled water, and glass beads with 5% of talc moisturized with KCl solutions of (e) 0.5 and (f) 1.5 S·m −1 conductivities over the frequency range from 1 to 500 MHz (selected data). The legends present the (volumetric water content) value of samples.
Despite the fact that the simulation gave the results up to 500 MHz that were consistent with the literature, in the case of the measurements, the upper working frequency of the probe was limited to 400 MHz because a slight increase of was observed above 400 MHz for several samples, which is not a physical behaviour. Additionally, some scattering of the values was observed at low frequencies. This scattering was more visible in the case of samples with higher water and salt contents. As expected, for each sample, the value of increased with an increasing water content [21,27,28]. Among all the samples, the highest value of for the saturation point was obtained for talc. According to the literature data, talc and glass beads exhibit hydrophobic character. This means that water does not adsorb inside the particles of these materials. However, water can be absorbed on the particles' surface due to the forces cohesion and adhesion [29]. The highest value of for the saturation point of talc can be explained by its smaller particle size as opposed to those of the glass beads.
In the next step, the -calibration curves were fitted in the form of a linear function of the square root of dielectric permittivity, as given by the following equation: The fitting equations were obtained for each of the materials, except for talc, at selected frequencies for data including distilled water and all sodium chloride solutions. For talc, it was only obtained for distilled water. Figure 4 shows the relations between values for the materials under test and an appropriate reference equation available in literature:  [4,11,28,30,31], dashed lines represent the fitted functions according to Equation (5), where the √ ε linear functions were presented in ε coordinates.
Wet2 sensor calibrations for mineral and clay soils at 20 MHz [4]: HydraProbe (HP) calibration for mineral and loam soils at 50 MHz [4]: ML3 ThetaProbe (TP) calibration for mineral soils at 100 MHz [31]: Calibration function for an FDR sensor obtained in the frequency range of 390-480 MHz [30]: Topp's equation obtained for TDR [11]: Calibration equation for 10 soils at 20, 50, and 100 MHz [28]: The application of pure talc as a reference material was under consideration. However, as can be seen in Figure 4a, talc stood out from the rest of considered materials and calibration curves. Moreover, it was difficult to remove the short seven rod probe from the talc sample, especially in the case of higher water content. The abandonment of further measurements of talc mixed with KCl solutions was caused by those two reasons. For the same value of θ, the value of ε for tested materials increased according to ε talc < ε glass beads < ε mixture glass beads and talc . The addition of talc to the glass beads caused a diminished air space volume between the glass beads, which, in turn, decreased the water content at saturation for the mixture of glass beads with talc. Additionally, both glass beads and the mixture of glass beads with talc gave results closer to well-known reference calibration functions. In the case of glass beads, the data points were more scattered in comparison to glass beads with the talc addition. Additionally, there were no significant differences between application 5% and 10% of talc. The fitted parameters of Equation (5), including the R 2 and RMSE values obtained for each material under different moisture and conductivity conditions for selected frequencies, are presented in Table 1. Generally, the materials under test gave high values of R 2 , which increased with an increase of frequency. The highest value of R 2 was observed in the glass beads with 10% of talc. It also occurred that the maximum value of the RMSE was obtained at 20 MHz. For each material, the RMSE decreased with an increase in frequency. The smallest differences between the maximum and minimum values of the RMSE among the tested materials was observed for glass beads with 10% of talc.
The ε spectra of selected materials are presented in Figure 5. Generally, the materials under test gave high values of R 2 , which increased with an increase of frequency. The highest value of R 2 was observed in the glass beads with 10% of talc. It also occurred that the maximum value of the RMSE was obtained at 20 MHz. For each material, the RMSE decreased with an increase in frequency. The smallest differences between the maximum and minimum values of the RMSE among the tested materials was observed for glass beads with 10% of talc.
The spectra of selected materials are presented in Figure 5.
Materials 2020, 13, x FOR PEER REVIEW 9 of 12 Figure 5. The (imaginary electric permittivity) values of (a) talc, (b) glass beads, glass beads with (c) 5% and (d) 10% of talc moisturized with distilled water, and glass beads with 5% of talc moisturized with KCl solutions of (e) 0.5 and (f) 1.5 S·m −1 conductivities over the frequency range from 1 to 500 MHz (selected data). The legends present the value of samples.
As expected, the addition of salt increased the value of the samples. Talc exhibited the lowest dielectric loss, which is characteristic of this material. In the next step, was obtained using the Levenberg-Marquardt algorithm as implemented in the MATLAB (Version 7.12.0.635, MathWorks Inc., Natick, MA USA) lsqnonlin function for nonlinear least square fitting. Then SI values for measurement materials were calculated using Equation (3). Following [32], was substituted with ′. For the calculation, the samples with < 3.5 were excluded from analysis, because the salinity As expected, the addition of salt increased the value of the ε samples. Talc exhibited the lowest dielectric loss, which is characteristic of this material. In the next step, σ b was obtained using the Levenberg-Marquardt algorithm as implemented in the MATLAB (Version 7.12.0.635, MathWorks Inc., Natick, MA, USA) lsqnonlin function for nonlinear least square fitting. Then SI values for measurement materials were calculated using Equation (3). Following [32], ε a was substituted with ε . For the calculation, the samples with ε < 3.5 were excluded from analysis, because the salinity index model does not work in the case of dry samples. The relations between σ b and ε at 400 MHz, presented in Figure 6 for the tested materials, were linear for all applied solutions with high value of R 2 > 0.95. The highest R 2 was obtained for glass beads with the addition of talc. Figure 5. The (imaginary electric permittivity) values of (a) talc, (b) glass beads, glass beads with (c) 5% and (d) 10% of talc moisturized with distilled water, and glass beads with 5% of talc moisturized with KCl solutions of (e) 0.5 and (f) 1.5 S·m −1 conductivities over the frequency range from 1 to 500 MHz (selected data). The legends present the value of samples.
As expected, the addition of salt increased the value of the samples. Talc exhibited the lowest dielectric loss, which is characteristic of this material. In the next step, was obtained using the Levenberg-Marquardt algorithm as implemented in the MATLAB (Version 7.12.0.635, MathWorks Inc., Natick, MA USA) lsqnonlin function for nonlinear least square fitting. Then SI values for measurement materials were calculated using Equation (3). Following [32], was substituted with ′. For the calculation, the samples with < 3.5 were excluded from analysis, because the salinity index model does not work in the case of dry samples. The relations between and at 400 MHz, presented in Figure 6 for the tested materials, were linear for all applied solutions with high value of R 2 > 0.95. The highest R 2 was obtained for glass beads with the addition of talc. Figure 6. Relation between the bulk electrical conductivity ( ) and ′ of the measured (a) glass beads, (b) glass beads with 5% of talc, and (c) glass beads with 10% of talc moistened with distilled water and KCl solutions. The value of (S·m −1 ) is given in the legend. The linear regression equation and R 2 for each material is also presented. Figure 7 shows the relation between SI and for materials under test. The relations were highly linear for each considered material, which was represented by high R 2 (Table 2), where l is the slope of a relation between SI and the . For the glass beads with 5% of talc, l = 0.0124 ± 0.000194, which matched the l = 0.0128 value that was calculated using TDR for sand [6]. Similar results were also obtained for glass beads with 0.26 mm [33].  Figure 7 shows the relation between SI and σ KCl for materials under test. The relations were highly linear for each considered material, which was represented by high R 2 (Table 2), where l is the slope of a relation between SI and the σ KCl . For the glass beads with 5% of talc, l = 0.0124 ± 0.000194, which matched the l = 0.0128 value that was calculated using TDR for sand [6]. Similar results were also obtained for glass beads with 0.26 mm [33].
Materials 2020, 13, x FOR PEER REVIEW 10 of 12 Figure 7. Relation between SI (salinity index) and of the measured glass beads, glass beads with 5% of talc, and glass beads with 10% of talc moistened with distilled water and KCl solutions. The value of is given in the legend. Glass beads 0.0109 1.48·10 -4 0.9998 Figure 7. Relation between SI (salinity index) and σ KCl of the measured glass beads, glass beads with 5% of talc, and glass beads with 10% of talc moistened with distilled water and KCl solutions. The value of σ KCl is given in the legend. Next, the repeatability of the dielectric properties of the most promising reference material was examined. In this case, the mixture of glass beads with 5% of talc moistened with distilled water was measured in three independent repetitions. It was calculated that the differences in the value of ε' between repetitions did not exceed 1%. This showed that the glass beads with talc met repeatability requirements, which is the one of important parameter for a reference dielectric material.

Conclusions
This paper presented the dielectric properties of glass beads and mixtures of glass beads with talc. It was concluded that the addition of talc to the glass beads provided more repeatability and stability of the material, which was observed through the diminished scattering of the data points in comparison to the results from the glass beads alone. The fitting equations obtained for the mixture of glass beads with talc, especially for the 5% addition of talc, gave results close to well-known reference calibration functions. It was also observed that pure talc stood out from the rest of the tested materials and the calibration curves, which makes pure talc a less optimal calibration material than the aforementioned mixtures. Additionally, the addition of 5% of talc also provided the highest linear correlations between σ b and ε , which is an important assumption of the salinity index model. It was observed that the slope l of SI versus the σ KCl relation was similar to the slope obtained for sand. In summary, the obtained data indicated that the glass beads with the addition of 5% of talc can be a promising material for the calibration and verification of dielectric methods and devices for measuring soil moisture. The addition of clay material to a mixture of glass bead with 5% of talc is planned in future studies in order to provide better modelling of low-frequency dielectric relaxation phenomena.