SSSAJ Grow Your Career with SSSA
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Lebron, I.
Right arrow Articles by Lesch, S. M.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Lebron, I.
Right arrow Articles by Lesch, S. M.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Lebron, I.
Right arrow Articles by Lesch, S. M.
Related Collections
Right arrow Soil Mineralogy
Right arrow Geochemical Processes
Right arrow Time Domain Reflectometry, TDR
Published in Soil Sci. Soc. Am. J. 68:1549-1559 (2004).
© 2004 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA

DIVISION S-1—SOIL PHYSICS

The Dielectric Permittivity of Calcite and Arid Zone Soils with Carbonate Minerals

I. Lebrona,c,*, D. A. Robinsonb, S. Goldberga and S. M. Lescha

a USDA-ARS, George E. Brown, Jr. Salinity Lab., 450 W. Big Springs Rd, Riverside, CA 92507
b Dep. of Plants, Soils and Biometerology, Utah State Univ., Logan, UT 84322
c Currently, Dep. of Plants, Soils and Biometerology, Utah State Univ., Logan, UT 84322

* Corresponding author (inma_lebron{at}yahoo.com)


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Measurement of soil dielectric properties, {epsilon} are widely used to estimate water content in soils from remote sensing and from in situ soil sensors such as time domain reflectometry (TDR). The mineral permittivity also plays an important role in geochemical dissolution and precipitation. Models used to estimate water content from soils often assume a value of 5 for the mineral permittivity {epsilon}s. However, calcite (CaCO3), a major constituent of some arid and semi-arid soils, has a permittivity of 8 to 9, nearly twice the permittivity of quartz ({epsilon}s = 4.6). We studied four soils, with micaceous mineralogy, but with two soils having approximately 40% calcite. We also measured the permittivity of Iceland Spar calcite ({epsilon}s = 9.1) and a microcrystalline calcite ({epsilon}s = 8.3), and use atomistic modeling to account for differences in permittivity based on the crystal density. We found permittivities for our soils to be in the range of 5.8 to 6.6, higher calcite contents resulting in increased permittivity. The estimated permittivity of the calcite in the soils was 7.4 to 7.9, lower than the highly crystalline samples. We estimate, for a soil with a porosity of 0.5 that assuming a permittivity of 5 instead of 6.6 will result in an overestimation of water content of about 1% at saturation. This demonstrates that a large quantity of pedogenic calcite (40%) in soil is unlikely to cause substantial error in the determination of water content using standard calibration equations. However, the lower dielectric permittivity predicted for pedogenic calcite may have consequences for the interpretation and understanding of geochemical processes.

Abbreviations: TDR, time domain reflectometry


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
THE DEVELOPMENT OF electromagnetic sensor technology toward the end of the last century increased interest in the measurement of soil properties using electromagnetic techniques. The determination of soil water content from the measurement of soil dielectric permittivity has become a widely accepted method (Topp et al., 1980; Topp and Ferre, 2002; Robinson et al., 2003). Investigations have also explored the use of dielectric measurement to estimate other soil properties such as cation exchange capacity (Fernando et al., 1977), soil aggregation (Miyamoto et al., 2003), and particle-size analysis (Starr et al., 2000). Understanding the dielectric properties of soils is also required for the improved interpretation of data from methods such as ground penetrating radar (Huisman et al., 2001; Davis and Annan, 2002) and microwave remote sensing (McNairn et al., 2002). Modeling of the dielectric properties of granular materials such as soils cannot only improve calibration models and techniques to monitor properties, but also lead to new methods of material characterization. When modeling a composite material such as soil it is required that the dielectric permittivity of all phases (solid, liquid, and gas) are known, especially if the model is to be used to predict water content or porosity. The permittivity of air and water are well documented (Lide, 1999). However, a value of 5 reflecting the permittivity of quartz (4.19 to 5.00; Keller, 1989) is often used as an input parameter for the permittivity of the solid for modeling (Friedman, 1998).


    Geochemical Importance of the Mineral Permittivity
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Water is an excellent solvent, primarily because of its high dielectric constant, approximately 80. For a continuum scale analysis, the inverse square Coulomb force between two charged atoms is reduced by a factor determined by the permittivity of the medium surrounding them. The free energy, w(r) of the Coulomb interaction between two charges Q1 and Q2 is:

[1]
where {epsilon}e is the permittivity of the surrounding or environment, {epsilon}o is the permittivity of free space (8.854 pF m–1) and r is the distance between the two charges. This indicates the importance of the permittivity of the environment in determining how the charges feel each other. A step further takes us to the Born or solvation energy of an ion (µi). The Born energy, a positive quantity, gives the electrostatic free energy of an ion in a medium of permittivity {epsilon}e:

[2]
where Q is the charge of the ion and, a is the radius of the ion. This can now be used to calculate the change in free energy ({Delta}µi, joules) of moving an ion from a medium of low dielectric constant, {epsilon}e to a medium of high dielectric constant {epsilon}b.


[3]

The transfer of the ion from {epsilon}e to {epsilon}b results in negative values and thus it is energetically favorable. For a full derivation of the above equations see Israelachvili (1992). This is the core to explaining why ionic solids dissolve in polar liquids. Molecular models describing calcite, as represented by Fisler et al. (2000), represent predominantly ionic bonding and include refinements for covalency and electronic polarization. These atomistic models of carbonate minerals still use a modified form of the Born model to describe atomic interaction within the crystal (Fisler et al., 2000; Cygan et al., 2002). Thus the contrast between the permittivity of the calcite crystal and the bathing solution may have a role in determining mineral solubility.

Many arid soils contain large quantities of salts, while the more soluble ones are prone to dissolving and being transported, more insoluble salts like calcite precipitate. Over time this can amount to large quantities when the water has excess in Ca and bicarbonates. Calcite in soils does not precipitate by crystal growth but by heterogeneous nucleation (Lebron and Suarez, 1996). That is the reason why in soils, well-developed pedogenic crystals of calcite are not encountered. The presence of water-soluble organic ligands and ions such as PO4–3 act as precipitation inhibitors by blocking crystal growth sites (Paquette et al., 1995; Lebron and Suarez, 1996, 1998). The presence of other cations like Mn, common in soil environments, may cause coprecipitation or formation of solid solutions with the result of mixed crystals of variable chemical composition (Lebron and Suarez, 1999). Consequently, it is expected that the dielectric permittivity of pedogenic calcite will differ from the permittivity of perfectly crystalline calcite. Crystals formed by crystal growth in the absence of impurities and without isomorphic substitutions provide pristine crystals where the ions are translationally symmetric and ordered without imperfections or abrupt interruptions. Tabulated permittivity values for calcite generally refer to pure crystalline samples.


    Measurement of Solid Phase Permittivity
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Permittivity values for other commonly occurring soil minerals are difficult to find and are often estimated using two-phase mixing models rather than measured. Olhoeft (1979) pioneered much of the estimation of rock mineral permittivity using two-phase mixing models to provide tables of estimated mineral permittivity values (Olhoeft, 1981; Carmichael, 1982; Nelson et al., 1989; Nelson and You, 1990; Nelson, 1992). These tabulated values are very useful in indicating that the permittivity of some minerals is much greater than others. However, the method is unreliable in terms of obtaining absolute values. It relies on applying a two-phase (air and solid) dielectric mixing model to calculate the permittivity of a granular sample repacked in air, with the permittivity extrapolated to zero porosity. This method pre-assumes the correctness of the mixing model. Soils prove difficult materials for accurately estimating their solid permittivity. Predicted permittivity values for the same mineral can vary considerably depending on the choice of mixing model. Particle shapes and sizes and the presence of clay minerals that may be partially hydrated can influence the accuracy of the estimates. An additional complexity in soil systems is the variable chemical composition of minerals. Isomorphic substitutions and solid solutions are intrinsic to pedogenic conditions. Consequently, the chemical formula for clays and minerals are not unique.

In recent work, Robinson and Friedman (2003) presented a method, using a TDR probe, to measure the solid dielectric permittivity of granular materials. This method measures the effective permittivity of granular samples when immersed in different dielectric liquids. The liquids are chosen to obtain values of effective permittivity lower and higher than the expected solid permittivity. The permittivity value for the solid is obtained when the fluid permittivity matches the effective permittivity of the mixture or by interpolation to that point. The method is very effective and produces permittivity values for granular quartz in close agreement with measurements presented in the literature for single crystals. Recently the method was refined for clay minerals and fine-grained materials (Robinson, 2004a).

Our objectives for this work are: (i) to quantify the contribution of pedogenic calcite to the solid dielectric permittivity of soils with calcite, and draw implications for estimating water content, and (ii) to compare the permittivity of the pedogenic calcite with pristine calcite specimens of different origin, and draw implications for soil geochemical processes.


    THEORY
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Two-Phase Single Solid Dielectric in Dielectric Solution
The immersion of a granular material in a dielectric fluid gives an effective permittivity, {epsilon}eff somewhere between the permittivity of the background {epsilon}e, and the permittivity of the granular inclusions {epsilon}s. The Maxwell-Garnett (1904) mixing model based on the Lord Rayleigh (1892) formula is commonly used for describing a two-phase mixture (Sihvola, 1999) and best suited to the present discussion:

[4]
where, f is the fraction of the solid, f = (1 – {phi}), {phi} is porosity. The model does not account for the effects of particle close packing but this is minimal when the contrast ({epsilon}e/{epsilon}s) between the two phases is <3 (Robinson and Friedman, 2002). The model assumes that the solid spherical inclusions ‘see’ only the permittivity of the background. In practice in a densely packed granular material the background seen by the solid is some combination of solid and fluid and its respective permittivity. The Maxwell-Garnett model provides an upper bound where the background has a higher permittivity than the inclusion and a lower bound when this is reversed.

Figure 1 shows the relationship between the effective permittivity, {epsilon}eff, and the fluid permittivity, {epsilon}e, for a granular sample with {epsilon}s = 5, packed to three different porosities. The 1:1 line is plotted on the diagram and it can be seen that all the lines converge at the value of 5 when the solid and fluid have the same permittivity. Below this point the effective permittivity of the mixture is higher than the solution and above this point the effective permittivity is lower than the permittivity of the immersion fluid. The main point demonstrated by the modeling is that for a fixed porosity the effective permittivity and the solution permittivity have the same value at a unique point equivalent to the permittivity of the solid.



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 1. Modeled effective permittivity (Eq. [1]) of a two-phase mixture plotted vs. the permittivity of the background solution at different porosities. The unique crossing point occurs when the permittivity of the fluid matches the permittivity of the solid.

 
Two-Phase Binary Solid Dielectrics in Dielectric Solution
Another aspect to consider when dealing with soils is how to account for the effective permittivity of a material with inclusions of different volume fractions and different permittivities. Sihvola (1999) presented a simple Maxwell-Garnett based model for spherical inclusions for this problem. He demonstrated that the permittivity of the solid was not simply the arithmetic average of the solid permittivities and their respective volume fractions:

[5]
where, f1 is the fraction of solid with permittivity {epsilon}1 and f2 is the volume fraction of solid with a permittivity {epsilon}2. Robinson and Friedman (2002) presented work on the dielectric properties of mixtures of glass beads and quartz sand. They demonstrated that in the case of low contrasts ({epsilon}2/{epsilon}1 < 3) between the inclusion permittivity values, as is the case with soil minerals, negligible error was obtained by referring to a single {epsilon}s, which is the arithmetic average of the permittivity of the minerals and their volume fractions.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Soils
Four micaceous/illitic soils were chosen for this study: Ramona and Clarence, from the USA, and two soils from the Ebro Basin (Spain). Classification and dominant mineralogy are shown in Table 1, particle-size distribution, Fe and Al oxide, organic matter, and calcite content are shown in Table 2. Granulometric determinations were made in a Sedigraph 5000ET (Micromeritics, Norcross, GA1), Fe and Al oxides were extracted by the method of Jackson et al. (1986). Inorganic and organic C was analyzed by a CO2 coulometer (UIC Corp., Joliet, IL). Mineralogical identification was performed by x-ray powder diffraction. We removed the calcite from a portion of the two soils from the Ebro Basin according to the method described in Kunze and Dixon (1986). Mineralogical identification was performed by x-ray powder diffraction to ensure total calcite removal (Fig. 2) .


View this table:
[in this window]
[in a new window]
 
Table 1. Classification and dominant mineralogy of the four soils studied, Ramona and Clarence soils according to U.S. Taxonomy.

 

View this table:
[in this window]
[in a new window]
 
Table 2. Particle-size distribution and other components of the soils in this study.

 


View larger version (35K):
[in this window]
[in a new window]
 
Fig. 2. X-ray diffraction patterns for the Ebro Basin soils EB1 and EB2 before and after the removal of calcite. Observe the disappearance of the calcite peak after treatment.

 
Calcite Minerals
The properties of the calcite minerals used in this study are shown in Table 3. Particle densities were determined using the excluded volume method described in Flint and Flint (2002). The calcite powder was obtained from Pfizer (Multiflex MM from Pfizer, Los Angeles, CA), the Iceland Spar was from Ward's (Ward's Natural Science Establishment Inc., New York). Characterization of the morphology, size, and crystallinity was done using microscopic and x-ray techniques. Figure 3a shows a micrograph of the Iceland Spar using a scanning electron microscope (AMRAY, 3200, AMRAY Inc., Bedford, MA). Crystallinity was analyzed by selected area electron diffraction for Iceland Spar (Fig. 3b) (SAED, Philips Analytical, Eindhoven, the Netherlands). Figure 3c shows the micrograph of the microcrystalline calcite using a transmission electron microscope, and the x-ray diffraction pattern using energy dispersive X-ray (EDX) (Fig. 3d) (Philips CM 300, Philips Analyitical, Eindhoven, the Netherlands).


View this table:
[in this window]
[in a new window]
 
Table 3. Physical and chemical properties of Iceland Spar and microcrystalline calcite (CaCO3).

 


View larger version (178K):
[in this window]
[in a new window]
 
Fig. 3. (a) Scanning electron micrograph of Iceland Spar. (b) Diffraction pattern for Iceland Spar. (c) Transmission micrograph for microcrystalline calcite. (d) Diffraction pattern for microcrystalline calcite.

 
Dielectric Immersion Fluids
We used the method described in Robinson (2004a) to measure the dielectric permittivity of fine-grained granular materials. The critical aspect in applying this method is finding fluids whose permittivities lie in a range between 1 and 15 to provide the dielectric background. Furthermore, the characteristics of the background fluids are important, they must not exhibit strong relaxation at the frequency of measurement (0.1–1 GHz). Robinson and Friedman (2003) used air, {epsilon} = 1; penetrating oil (WD-40, WD-40 Company, San Diego, CA), {epsilon} = 2.3; methylene chloride, {epsilon} = 8.8; and acetone, {epsilon} = 20.8. They also found that the alcohols were unsuitable because of relaxation in the frequency bandwidth of the TDR. In this work mixtures of fluids were used to obtain the low permittivity background dielectric solutions. The main fluids used were: air, {epsilon} = 1; corn oil, {epsilon} = 2.4' dicholoromethane, {epsilon} = 9.1; and acetone, {epsilon} = 20.8, at 25°C. The acetone and corn oil are miscible, which allowed us to prepare dielectric background solutions with a permittivity between 2.4 and 20.8.

Measurement of the Effective Permittivity and TDR Probe Construction
A Tektronix (1502C) TDR cable tester and a custom made coaxial TDR probe were used throughout the experiments to measure the effective permittivity of soil and clay suspensions. The TDR was connected to a personal computer, which was used to collect and analyze waveforms using software developed by Heimovaara and de Water (1993). Cable tester and probe were connected via a 0.75 m, (50-ohm RG 58) coaxial cable. The probe had an internal stainless steel electrode (A in Fig. 4) , 0.003175 m in diameter and 0.16 m in total length with 0.152 m projecting above the chemical resistant delrin probe head (Dupont, Wilmington, DE). The outer part of the TDR probe was a brass cylinder (B in Fig. 4) 0.156 m long with an internal diameter of 0.085 m. This cylinder had a brass collar soldered onto it at one end so that a screw fitting (C in Fig. 4) could tighten the coaxial cell onto the microwave electrical connector (E in Fig. 4). The central electrode had a male shaped end that fitted into the female microwave connector (E in Fig. 4). The microwave electrical connector was high quality and rated to pass frequencies of up to 11 GHz. The coaxial TDR probe was designed with a screw and push fit to allow separation and cleaning. A small probe was developed to minimize the amount of clay mineral required for the measurement. The volume of the cell was 8.0134 cm3. The TDR probe was calibrated for effective length (0.154 m) using air and deionized water (Heimovaara, 1993; Robinson et al., 2003).



View larger version (132K):
[in this window]
[in a new window]
 
Fig. 4. Photograph of the equipment specifically made for soil dielectric measurements. Internal stainless steel electrode 0.003175 m in diameter and 0.16 m in total length with 0.152 m projecting above the chemical resistant delrin probe head at the base (A); outer brass tube, 0.156 m long with an internal diameter of 0.085 m (B); retaining nut fixing B to E (C); male connector (D); microwave electrical connector (E); coaxial cable (F); syringe assembly for mixing suspension(G); mixing rod (H); syringe used to draw fluid out of the coaxial cell, a similar syringe was used with a piece of tygon tube to connect to G, to apply tension to de-air the suspension (I).

 
Measurement Procedure
The measurement procedure relies on making a suspension of the mineral in the background dielectric fluid. In the first step, soils and mineral samples were oven dried at 150°C overnight to remove any hygroscopic water. The samples were removed from the oven and allowed to cool in a desiccator. A 10-g subsample from the dessicator was removed and accurately weighed in a plastic weighing boat on a four-point balance. The sample was then gently poured into the coaxial TDR probe until it was flush with the top, thus forming a loose packing, measurements were then made. The weighing boat was then reweighed with the remaining material to determine the mass of sample poured into the probe. The solid was then removed from the coaxial probe and the probe was dismantled and cleaned.

The second set of measurements was then made, initially by measuring the permittivity of the liquid. The solid–liquid suspension was then made by removing the fluid from the probe using a large syringe with a 0.16-m tube attached. The suspension was made using this fluid, maintained at the same temperature in a 10-mL syringe (Fig. 4, G). The nozzle of the syringe was blocked with a stopper and 4 mL of liquid was placed in the 10-mL syringe. The weighed mineral sample was poured into the syringe and mixed using a stainless steel stirring rod (Fig. 4, H). The suspension sample was now made up to a mark on the syringe representing 8.0 cm3. The stopper was removed from the nozzle and another syringe attached to the syringe with the suspension via a tygon tube. A small tension was applied using this syringe to draw out any remaining air. As bubbles formed, the syringe with the suspension was gently tapped to remove these bubbles. This suspension was then injected into the coaxial TDR probe; the suspension had the same solid fraction as the sample measured with air as the background. If the sample was slightly short of reaching the top of the coaxial probe a little more background fluid was added until flush with the top of the probe. Measurements were again made with the TDR and after this the coaxial probe was dismantled and cleaned thoroughly. This procedure was then performed several more times with different background dielectric fluids. Statistical analysis was performed on the grouped data set to determine the 95% confidence interval for each estimate of permittivity. The methodology for this is described in the appendix.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Dielectric Permittivity of Calcite Minerals
The measured permittivity values for Iceland Spar and microcrystalline calcite in the different fluids are presented in Fig. 5 . We used an empirical model of the form {alpha}{epsilon}eß (where {alpha} and ß are constants) to interpolate between the measured values to determine the permittivity for the Iceland Spar calcite ({epsilon}s = 9.1), and for the microcrystalline calcite powder ({epsilon}s = 8.3). Figure 5 shows also the uncertainty around the fitted model with the 95% confidence interval. The appendix presents detailed information about the quantification of the uncertainty. The error bars on the microcrystalline calcite show a large uncertainty, however, one of the measurements is located directly underneath the estimate (not seen in Fig. 5) and a fifth measurement was made at a slightly higher permittivity (seen in Fig. 5 underneath the estimate for the Iceland Spar) giving a good correspondence between the measurements and the model. We will also demonstrate that we would expect the solid permittivity value for the micro-crystalline calcite to be lower than for the Iceland Spar as the crystals are less dense. The obtained values are compared with a range of values from the literature in Table 4. The value for Iceland Spar was in good agreement with values obtained from atomistic modeling (Fisler et al., 2000) but was a little higher than values measured using the immersion method for a single calcite crystal (Schmidt, 1902) or values reported in Keller (1989). The values for both the permittivity and particle density quoted from Olhoeft (1981) are in poor agreement with a permittivity value of 8 to 9 and the more generally accepted value of 2.71 g cm–3 for the density of crystalline calcite. This may reflect the method of repacking to estimate the permittivity and also the quality of the samples, which were probably natural and may have contained many impurities. The values for the particle density would certainly indicate this. We are confident of our measured Iceland Spar value since the calcite grains were 1 mm in diameter ensuring no air entrapment, which could have biased results. The value for the microcrystalline calcite falls within previously documented permittivity ranges for calcite.



View larger version (27K):
[in this window]
[in a new window]
 
Fig. 5. Effective permittivity of coarse-grained (500 µm) natural calcite sample (Iceland Spar) obtained from a single calcite crystal that was crushed. The calcite powder was a synthetically made sample that was very fine-grained (0.1 µm). The solid symbols represent the estimated permittivity with 95% confidence estimates, Iceland Spar {epsilon}s = 9.1; Micro-crystalline calcite {epsilon}s = 8.3. The forth point for the micro-crystalline calcite lies directly under the estimate, giving perhaps greater confidence in the value of 8.3 than the confidence estimate might suggest.

 

View this table:
[in this window]
[in a new window]
 
Table 4. Measured and modeled values for the permittivity of calcite obtained from the literature.

 
Dielectric Permittivity of Soils
We measured the permittivities of the soil suspensions and used an empirical power function to interpolate between the data points. The obtained curve was then used to calculate the interception of the fitting function with the line y = x, this value provides the dielectric permittivity of the soil mineral phase (Robinson and Friedman, 2003). We obtained {epsilon}s = 6.0 for the Ramona soil and {epsilon}s = 5.8 for the Clarence soil, values in close agreement with their mineralogy (Fig. 6) and values previously measured for illite (Robinson, 2004a).



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 6. Effective permittivity of Clarence (Cl) and Ramona (Ra) soils immersed in different dielectric fluids. The solid symbols represent the estimated permittivity with 95% confidence estimates, Cl {epsilon}s = 5.8; Ra {epsilon}s = 6.0.

 
Figures 7a and 7b present the results for the two Ebro Basin soils. The graphs show the permittivity measured for the unaltered soil mineralogy ({epsilon}s = 6.6 for EB1, {epsilon}s = 6.3 for EB2) and for the soil after treatment to remove the calcite. The permittivity of the soil mineral phase after the removal of the calcite was {epsilon}s = 5.7 for EB1 and 5.8 for EB2, in good agreement with their micaceous mineralogy. However, the permittivity of the soil mineral phase with calcite was not as high as we would expect considering that they have 35 and 41% calcite.



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 7. Effective permittivity of the Ebro Basin soils before (EB1C and EB2C) and after (EB1NC and EB2NC) treatment with acid to remove the calcite. The solid symbols represent the estimated permittivity with 95% confidence estimates, EB1C {epsilon}s = 6.6; EB2C {epsilon}s = 6.3; EB1NC {epsilon}s = 5.7; EB2NC {epsilon}s = 5.8.

 
In Table 5 we present results where we predicted the permittivity that we would expect for the whole soil using a simple arithmetic averaging, {epsilon}s = {epsilon}1f1 + {epsilon}2f2 where {epsilon}s is the permittivity of the whole soil, f1 is the fraction of solid with permittivity {epsilon}1 (micaceous clay), and f2 is the volume fraction of solid with a permittivity {epsilon}2 (calcite). This assumed that the calcite in the soil had (i) the permittivity of the Iceland Spar, and (ii) the permittivity of the microcrystalline calcite. Finally, for (iii) we calculated the permittivity that the soil calcite must have had to give the permittivity value measured for the whole soil. These calculated values for the pedogenic calcite were both less than 8 and much lower than expected.


View this table:
[in this window]
[in a new window]
 
Table 5. Dielectric permittivity of the Ebro Basin soils. Calculations of the soil dielectric permittivity ({epsilon}s) were made based on: {epsilon}s = {epsilon}1f1 + {epsilon}2f2 where {epsilon}s is the permittivity of the whole soil, f1 is the fraction of solid with permittivity {epsilon}1 (soil without the calcite), and f2 is the volume fraction of solid with a permittivity {epsilon}2 (calcite).{dagger}

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Differences Between Calcite Permittivity Values
In the introduction (Eq. [1]–[3]) we proposed that the permittivity plays a role in determining the solubility of a crystal in a background fluid. As the permittivity contrast between the crystal and in our case water becomes greater so the change in free energy becomes more negative, and so the crystal is likely to become more soluble. Particle size has been shown to play an important role in the solubility of calcite. Chave and Schmalz (1966) demonstrated that as size decreased so the apparent surface energy increased. Their work indicated that calcite particles under 0.1 µm would be unstable in aqueous systems, tending to be more soluble. They suggested that as particles became smaller the edges and corners play a greater role in the surface energy. They also considered that grinding of calcite could cause a phase transition. It has been shown that calcite can be inverted to aragonite by grinding with a pestle and mortar (Jamieson and Goldsmith, 1960).

Our samples were carefully chosen to provide two extremes of calcite particle size (approximately 1000 and 0.1 µm) without grinding. However, our measurements indicated that although both samples were calcite, the 0.1-µm particles had a lower density. We use this to provide an explanation as to the difference in measured permittivity for each calcite sample. The confidence interval on the microcrystalline calcite measurements overlapped with the Iceland Spar measurements. However, we believe the measurements accurately reflect the calcite properties, and to explain the differences between the measured values of calcite permittivity we use atomistic modeling. This allows us to understand how changes in the crystal density impact the permittivity. Highly sophisticated atomistic models are being continually updated to predict the physical properties of minerals including the permittivity (Gale and Rohl, 2003). However, a good understanding of the important mechanisms can be gained using a more simplistic, semi-physical modeling approach. Polarization mechanisms in solids can be divided into three major categories, electronic, ionic, and orientational; it is the first two of these that are of interest in calcite. The link between the mineral dielectric polarizability, {alpha}D and the measured real part of the dielectric permittivity of the solid, {epsilon}s is given by the Clausius–Mosotti equation:

[6]
where Vm is the molar volume in cubic angstroms {in many reference texts molar volume is given in cm3 and can be converted to a molar volume in cubic angstroms by: [mole mass/x-ray density (g cm–3)]/Avogadro's number] x1024}, b is defined as 4{pi}/3, and {epsilon}s is the real part of the complex dielectric permittivity of the solid (Shannon, 1993; Robinson, 2004b). This simplistic model makes several assumptions; an isotropic distribution of spherical ions in space and all the calculated values of permittivity are average values, taking no account of crystal orientation. The molar volume Vm (cubic Angstroms) is a function of the crystal density indicating the importance of this in determining the polarizability of the crystal. Molar volumes for many minerals can be found in Weast (1984) for example. Shannon (1993) presented a list of ion polarizability values derived from measurements using single crystals of oxide minerals. No measurements were made for the carbonate ion, however, O2 dominates the polarizability so we can assume that the calcite crystal polarizability is dominated by the Ca ion and three oxygen ions summed using the oxide additivity rule (Shannon, 1993). The values of ion polarizabilities used were, 3.16 cubic angstroms for Ca and 2.01 cubic angstroms for O2 giving a total polarizability of 9.19 cubic angstroms. These values were used in the Clausius–Mosotti equation to predict mineral permittivity values according to:

[7]

The model is sensitive to the crystal density, through Vm. By adjusting the density value from 2.71 (Iceland Spar) to 2.66 g cm–3 (microcrystalline calcite) the estimated permittivity using Eq. [7] drops from 8.5 to 8.0. The solid permittivity estimate of 8.5 is exceptionally good considering the simplifications in the modeling approach. The half unit, permittivity difference between the two values, caused by altering the molar volume can account fully for the difference in permittivity observed between the Iceland Spar and microcrystalline calcite sample. Thus we suggest that differences in crystal density can account for the observed differences in the crystal permittivity.

By using two confirmed calcite samples we can rule out phase transition as being the mechanism leading to increased permittivity. We would also suggest that based on the arguments presented in Eq. [1] through [3] that our observations are in keeping with measurements of calcite solubility in natural waters in arid environments. In Eq. [3] we demonstrated the importance of the contrast between the permittivity of the crystal {epsilon}e and the background solution {epsilon}b; the greater the contrast the more energetically favorable the transfer is. Thus calcite with a lower permittivity may be part of the solubility mechanism that makes it more soluble in the natural environment. Suarez (1977) and Suarez and Rhoades (1982) found supersaturation of soil solutions to be common in semiarid environments, explaining this by the simultaneous dissolution of silicate minerals and subsequent precipitation of heterogeneously nucleated calcite. Suarez (1983) also found this to be the case for water from areas of the lower Colorado River. Further work might explore in more detail the role of particle size and mineral permittivity on calcite solubility.

We know that soils are not an ideal environment for crystal growth and although micrographs of pedogenic calcite show the classic rhombohedral structure it is likely that the crystal may contain many internal dislocations and defects and these could also reduce the permittivity of the crystal. Accurate dielectric measurement of minerals might be a way of determining their crystallinity, this may help in our understanding of geochemical processes since precipitation and dissolution processes are affected by the permittivity of both the mineral and liquid phase (Israelachvili, 1992).

Implications for Water Content Determination
Equation [4] can be used to estimate the effective permittivity of a two-phase mixture of water and mineral. We found that calculating the effective permittivity of a saturated soil with a porosity of 0.5, using a value of 5.0 for the mineral permittivity then a value of 6.6, representing the soil with 41% calcite, altered the determined effective permittivity from 35.0 to 36.1. Assuming that this would be the maximum difference between permittivity, and applying a refractive index style mixing formula a difference in determined water content of about 0.01 m3 m–3 results. This indicates that even large quantities of pedogenic calcite (40%) are unlikely to substantially affect the determination of water content using a standard calibration.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Measurements of the dielectric permittivity of calcite samples and arid zone soils with pedogenic calcite are presented. Iceland Spar calcite had a measured permittivity of 9.1 and a microcrystalline calcite had a measured value of 8.3. We explain this difference using an atomistic model based on the difference between the measured density of the two samples. The presence of pedogenic calcite in soils was found to increase the permittivity from 5.7 and 5.8 without calcite to 6.3 and 6.6 for two soils with 35 and 41% calcite, respectively. Modeling calculations suggested that the permittivity of the pedogenic calcite was between 7.4 and 7.9, lower than either of the pure samples. Imperfections in the pedogenic calcite crystal, distortions due to isomorphic substitutions, lower crystal density and/or solid solutions within the crystal might explain these lower permittivity values. We speculate that the lower permittivity observed might be a mechanism that makes pedogenic calcite more soluble and perhaps in part accounts for observed supersaturation of natural waters with respect to calcite.

We also examine the implications of determining water content from dielectric mixing models assuming a permittivity of 5 for the permittivity of the solid. Our calculations suggest that assuming a solid permittivity of 5 instead of 6.6 for a soil with a porosity of 0.5 would only lead to an overestimate of water content of about 0.01 m3 m–3 at saturation.


    APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
Approximate 95% Confidence Interval Estimates
For this study, we assumed that the following multi-parameter nonlinear model could adequately describe the laboratory response data:

[A1]
where {epsilon}eij represents the jth solution permittivity measurement in the ith experiment, {epsilon}effij represents the corresponding effective permittivity measurement, {alpha}i and ßi represent the model parameters, and {eta}ij represents the random noise component (i.e., the error), which was assumed to be independent and proportional to the measured response. We then log transformed both sides of [A1] to yield:

[A2]
which is identical to a standard analysis of covariance (ANOCOVA) model. (A residual assessment performed on this model indicated that the transformed errors satisfied the usual normality assumptions). We then solved [A2] for {epsilon}e = {epsilon}eff to produce the crossover point for each experiment; that is,

[A3]

An approximate 95% confidence interval for each crossover point was then constructed using a first-order Delta Method calculation (Casella and Berger, 1990). Specifically, the approximate variance of each crossover point estimate was calculated as

[A4]
and ±2 times the square root of this variance estimate was used to create the log confidence intervals. Note that the ß parameter variance and covariance estimates used in [A4] were derived from [A2]; that is, the fitted ANOCOVA model. Back-transformed point and interval estimates were then calculated from the variance estimates produced by [A4].


    ACKNOWLEDGMENTS
 
The authors acknowledge funding provided in part by a USDA NRI grant 2002-35107-12507.


    NOTES
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 
1 Trade names are provided for the benefit of the reader and do not imply any endorsement by the USDA. Back

Received for publication October 3, 2003.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 Geochemical Importance of the...
 Measurement of Solid Phase...
 THEORY
 MATERIALS AND METHODS
 NOTES
 RESULTS
 DISCUSSION
 CONCLUSIONS
 APPENDIX
 REFERENCES
 





This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal