Microwave complex permittivity and anisotropy of conifer wood chips vs moisture content: experiments and modeling

The complex microwave permittivity—including anisotropy- of wood chips of softwood has been measured for different moisture contents in the band 0.75 to 2.5 GHz using an ultra-wide band radio transmission technique. The real and imaginary parts increase monotonically with moisture content. The wood chips are oriented by gravity, which gives anisotropic permittivity. The anisotropy ratio of the real part increases from 1.1 to 1.6 with moisture content from 0 to 120%. The anisotropy ratio of the imaginary part is around 2.5 at all moisture contents. Effective medium models were used to model the permittivity. The Bruggeman, and two versions of the Maxwell Garnett model gave good results at low moisture content (below the fiber saturation point). Above the fiber saturation point only the Bruggeman model gave results in agreement with experiments. The difference in model performance suggests that the free water does not follow the wood chips geometry.


Introduction
Wood chips are used in pulp and paper industry and as fuel in district heating. The moisture content of wood chips influences the combustion and emission of carbon dioxide [1,2] and has been found to be an important quality parameter in process control in pulp mills [3]. It may also change during storage [4]. To measure the moisture content has therefore attracted interest and indirect methods like near-infrared spectroscopy, X-ray absorption, and radio frequency and microwave methods are used [5,6].
The microwave permittivity of wood chips changes significantly with moisture content [7,8]. The measured complex permittivity of sawdust in the range 0.5 to 15 GHz was reported in [7] and varied with moisture.
Ottosson et al. [8] determined the permittivity (real part) of wood chips in volumes of several m 3 using ultra wideband (UWB) radar; the permittivity was found to be anisotropic, which was caused by the orientation of the woodchips by gravity, and the anisotropy increased with increasing moisture content.
Effective medium models or mixing models are used to calculate an effective dielectric permittivity of inhomogeneous composite materials [9,10]. The permittivity of the constituent materials as well as their volume fractions are used in the models. There are numerous effective medium models and if a specific model is suitable for a specific composite depends on various parameters, such as the geometry, volume fractions and permittivity of the constituents.
Paz et al. [11] measured the permittivity of woody biomass at 0.3-0.8 GHz and found that the Maxwell Garnett effective medium model was suitable for modeling the real part of permittivity vs moisture content. The Rönnow et al. Journal of Wood Science (2022) 68:22 inclusions were randomly oriented ellipsoidal inclusions, which result in an isotropic permittivity.
In this paper, we present measurements of the complex permittivity of wood chips of different moisture content at microwave frequencies (0.75-2.5 GHz). In particular, we present data for the anisotropy. We use an UWB radar equipment in transmission mode, a method that is suitable also for in line industrial applications and large volumes. We also use effective medium models for the complex anisotropy. To our knowledge, experimental data and effective medium modeling of the anisotropy of the complex permittivity of wood chips vs moisture content have not been reported before.

Wood chips
We analyze wood chips made of stem wood of predominantly Norway spruce (Picea abies) with possible small amounts of Baltic pine (Pinus sylvestris) from middle Sweden, i.e., conifer or softwood. The woodchips were produced for industrial purpose and in particular as fuel in district heating. It was, therefore, not possible to evaluate the permittivity and anisotropy of the wood before the it was cut into wood chips and used in the study. Figure 1 shows wood chips of stem wood of the type investigated here, where the left image is a top view and the right a side view. The side view is blurred since it is taken through transparent plastic used in the lab, in order to have the wood chips oriented as in the radio measurements (describe below). Notice that the individual wood chips have the shape of long rods in most cases. The length is in the range 1-10 cm. The short side of the individual wood chips is 10-20% of the length. The wood chips are not randomly oriented in all directions. Instead, they are oriented in such a way that the long side is predominantly horizontal (perpendicular to gravity) [8].

Experimental equipment
During the measurements the wood chips were placed in a plastic box of the type euro container (an industrial stacking container conforming to the VDA 4500 standard [12]). Four antenna housings were placed on each side of the box for radio transmission measurement (see Fig. 2). The box had the interior size of 26.8 × 36.8 × 31.5 cm 3 . The antenna housings were made of steel and had the size of 14.0 × 14.0 × 14.0 cm 3 . Each housing included one balanced Vivaldi-antenna, which could be rotated to enable the measurement of the permittivity parallel and perpendicular to gravity (cf. [13]). The transmit and receive antennas were both in the parallel or perpendicular orientation, in the respective measurement cases. The system was used in transmission mode, i.e., the antenna on one side is transmitting a signal that propagates through the wood chips and is received by the opposite antenna.
The antennas were connected to a radar system called DiRP (digital radar processor) that is based on M-sequence technology (cf. [14]), to generate the UWB signal. The sampling frequency was 25 GHz and the signal's center frequency was 2.5 GHz in a 3.5 GHz wide band. In this frequency range the signal level after propagating through wood chips is high enough to be detectable for wood chips of both high and low moisture content. Higher frequencies would be more attenuated resulting in too low signal levels. The antenna size is practical in laboratory and industry applications. Lower frequency would require larger antennas.
A conventional cement mixer (Biltema 17-686), Fig. 3 (left), was used to mix wood chips and water to get wood chips of moisture content. A scale (Dini Argeo IP67/IP68) was used to measure the weight of water and wood chips, Fig. 3 (right).

Experimental procedure
To characterize the wood chips, 100 pieces of wood chips were randomly selected and the dimensions were measured using a caliper. The wood chips were also put in a box of transparent plastic without any lid and photographed (see Fig. 1

above).
The moisturization process started by placing 8975 kg of relatively dry wood chips (starting chips) in a larger conditioning box. This amount was used in the successive moistening. We use the dry based moisture content, where m w and m d is the weight of the wet and dry wood chips, respectively.
Initially, the moisture content of the starting chips was measured by taking samples (à 500 g) that were weighted before, denoted m w1s , and after drying, denoted m d1s . The drying was performed during 24 h in + 105 °C, which is the approved industrial drying procedure in Sweden [15]. The moisture content of the starting chips, M d1 , was determined using Eq. 1. The moisture content of the starting chips was M d1 = 2.23%. The dry weight, m d1 , of the entire amount of starting chips was determined from Eq. 1, using M d1 and the measured weight of the starting chips, m w1 .
The moisture content of the wood chips was increased in a number of steps. In each step, n, the moisture content, M dn , was determined from the measured weight, m wn , and the initial weight, m w1 using Eq. 1. The wood chips were the same in all steps.
The moistening process was performed by sprinkling 1 L of water on the wood chips during mixing by the cement mixer (Fig. 3a) for each moistening cycle. Altogether, 10 moistening cycles were performed. Each succeeding moistening cycle (except for the first one) was based on preceding moistened wood chips. Each moistening cycle took about 30 min and the whole moistening experiment was carried out for 8 h (one working day). In each cycle, following acts and measures were made: 1. Wood chips and water were mixed in the cement mixer ( Fig. 3a). Water was sprinkled over the wood chips, 2. Wood chips were moved from the cement mixer to the conditioning box and mixed again, 3. The weight, m wn , of the wood chips in the conditioning box was measured (Fig. 3b),  4. Wood chips were moved from the conditioning box to the sample box and mixed again, 5. The sample box was packed by shaking the box towards the floor three times and adding more wood chips until the box was filled. Approximately half of the amount of wood chips from the conditioning box was needed to fill the sample box, 6. Radio transmission measurements with antennas in vertical and horizontal direction were performed on the sample box (Fig. 2a), and 7. The weight, m ′ wn , of the wood chips in the sample box was measured.
The mixing in both step 1 and 2 was made to achieve a homogeneous moisture distribution in the sample box.
A conventional error analysis gives the error in the measured moisture content where w is the error in the weight measuring.

Determining permittivity
We determined the complex permittivity of wood chips from transmitted radar pulses. We describe the main steps of the determination method. Details can be found elsewhere [8,13]. We model the wave propagation using the complex refractive index, ñ . The real and imaginary parts of the permittivity, ε = ε ′ + iε ′′ , are related to the complex refractive index ñ = n + iκ by or ñ = √ε . We determine ñ from measured radar pulses that have propagated through the wood chips. The time delay gives the real part, n, [8]. The imaginary part, κ, is obtained from the relative damping vs. frequency [13].
In a reference measurement with air in the test box, a reference pulse, y 0 (t) , is measured vs time, t. In a second measurement a pulse, y 1 (t) , is measured after it has been transmitted through the wood chips.
The refractive index, n, is determined from the difference in time delay between y 1 (t) and y 0 (t) . In order to make the time estimate more robust, we use a technique in which the Fourier transform of a transfer function, h(t) was calculated as [16] (2) where Y 1 (ν) and Y 0 (ν) are the Fourier transforms of y 1 (t) and y 0 (t) , respectively, ν is the frequency, and * denotes the complex conjugate; h(t) is obtained by taking the inverse Fourier transform of H(ν) . The parameter γ is small and real valued and works as noise filtering. The time of the peak value of h(t) is the difference in time between y 1 (t) and y 0 (t) , i.e., T 1 -T 0 . The refractive index is determined as [8] We estimate the error in the refractive index as: where z is the thickness of the box, and ∆z its error; ∆ (T 1 -T 0 ) is the error in the determined time difference.
The imaginary part, κ, is determined from the damping of y 1 (t) relative to y 0 (t) as described in [13]. The ratio of the signals' Fourier transforms is analyzed, where β 0 and β 1 are functions that that take into account the effects such as finite sample size and finite distance between the antennas [13]. These coefficients have small frequency dependence. We neglect all frequency dependence in H except that from the attenuation in the wood chips, which is given by the slope α = −2πκz/c 0 . Taking the natural logarithm of Eq. 7, we get that ln Y 1 (ν) Y 0 (ν) depends linearly on ν with the slope α = −2πκz/c 0 . Thus, we determine the slope, α by the least square method and then calculate κ = −αc 0 (2πz) . We estimate the error in κ as where �α is the error in the fitted coefficient of a straight line by linear regression (see e.g., [17]) and z is the error in the thickness of the wood chips in the test box.
The error in the permittivity is obtained from putting the errors in Eq. 8 into Eq. 3, To describe the anisotropy of ε , we the calculate the anisotropy ratios where ε ′ hor (ε ′′ hor ) is the real (imaginary) part of ε for horizontal E-field, and ε ′ ver (ε ′′ ver ) for vertical E-field, respectively.

Effective medium theory
Effective medium models (or dielectric mixing models) are used for calculating the complex dielectric permittivity, ε = ε ′ + iε ′′ , of inhomogeneous media of different materials, in which the length scale of the inhomogeneities is small compared to the wavelength of the electromagnetic radiation [9,10]. Effective medium models like Bruggeman and Maxwell Garnett are derived from Maxwell's equations for electromagnetic wave propagation and are quasistatic approximations. The models contain physical parameters and a comparison between modeled and experimentally determined permittivity can, hence, give insight into the physics of the inhomogeneous mixing. We use three different effective medium models that include inhomogeneities that give rise to anisotropic permittivity. All models have as input parameters the permittivity and volume fractions of the different materials (in our case air, wood and water) and the depolarization factors that model the effect of anisotropic inhomogeneities (in our case the elongation of the wood chips and the air and water in between them). We use literature values for the permittivity of wood at different moisture content in the models, together with volume fractions that we estimate from the measured masses.
The relaxations mechanism of wood with moisture is complex (see e.g., [18]). To describe the effect on the permittivity we use the concept of bounded and free water [19]. At moisture contents below the fiber saturation points the water molecules are adsorbed to the cell walls and do not interact as free water molecules with the E-field. Above the fiber saturation point, some water molecules are in the cell wall cavities or on the wood chips' surfaces and interact with the E-field as in free water. The fiber saturation point is at moisture contents of around 30% [19].
Wood chips have anisotropic permittivity. There are two physical properties that may give anisotropic effective permittivity. The first is that wood has different permittivity in different directions, or more precisely, for the E-field 1) radial (to the annual rings); 2) tangential (to the annual rings), and 3) parallel to the fiber direction or tangential to the stem [20]. The second is that (10) the wood chips that are in the shape of plates or needles are oriented by gravity such that the effective permittivity becomes different for the E-field parallel or perpendicular to gravity [8]. We use an isotropic equivalent permittivity in the models. In [21], different methods to describe an isotropic equivalent permittivity of wood were compared. It was found that the arithmetic mean of the permittivity of the three main axes differs only by a few percent from other physically motivated methods. We therefore use the arithmetic mean in the cases when we use an isotropic permittivity of wood, i.e., ε w = 2ε �,w + ε ⊥,w where ε ,w and ε ⊥,w are the permittivity of the wood with the E-field parallel and perpendicular to the fibers, respectively.
The first model that we use is a multiphase Bruggeman (BG-MPA) effective medium model for spheroidal inclusions [10,22]. This model has the advantage of being symmetric in the constituents, i.e., the analytical expressions do not distinguish between a host material, like air, and a particle material, like wood. This feature is advantageous in our case, since the volume fractions of the constituents may vary with water content. The spheroidal inclusions-in contrast to spherical ones-give anisotropic effective permittivity. The BG-MPA gives an effective complex permittivity for three directions, ε e,x , ε e,y , and ε e,z , which are obtained by solving for ε e,j for j = x, y, z . In Eq. 11, N is the number of different phases or materials, f n is the volume fraction of material n, ε n is the complex permittivity of material n, and L n,j is the depolarization factor for inclusions of material n in direction j. Notice, that f 1 + f 2 + . . . f N = 1 and that L n,x + L n,y + L n,z = 1 for each n. For isotropic inclusions, L n,x + L n,y + L n,z = 1 3 . For infinitely long needles in the z-direction, L n,x + L n,y = 0 and L n,z = 1 . When solving Eq. 11 an equation of order N in ε e,j has to be solved and consequently multiple roots are obtained. We solve the equations analytically for N = 2 and numerically for N = 3. Only one root is physically correct and we use the method in [23] to select the correct one. In the modeling below we use ε 1 for the permittivity air, ε 2 for wood, and ε 3 for water, respectively.
The second model that we use is the Maxwell Garnett model for layered spheroidal inclusions (MG-LSI) [24]. It is derived for a geometry of a spheroidal with a coating in a host material. It resembles the wood particles in air with the assumption of free water on the particles' surface. It gives the effective complex permittivity [24]: (11) N n=1 f nε n −ε e,j ε e,j + L n,j ε n −ε e,j = 0 where ε 1 is the complex permittivity of the host material (air), ε 3 of the outer layer material (water) and ε 2 of the core material (wood). f 1 , f 2 , and f 3 are the volume fractions of materials 1, 2, and 3, respectively, and, as for the Bruggeman model, f 1 + f 2 + f 3 = 1 . L j,n are as before the depolarization factors of material n in the direction j (j = x, y, or z). The third model is an anisotropic Maxwell Garnett model for spheroidal inclusions of multiple types and multiple depolarization factors (MG-MPA) [25]. The inclusion can also have different orientation and distributions of orientation. Paz et al. [11] used the Maxwell Garnett model for sawdust with randomly oriented anisotropic particles, causing an isotropic effective permittivity. We use the model assuming that all inclusions are aligned to the coordinate system of the measurements system (which is given by the orientation of the antennas, the box and the wood chips' orientation by gravity). The effective permittivity in the x,y, and z-directions is given by the diagonal elements of the matrix ε e , where ε 1 is the permittivity of the host material (air), I is the identity matrix and the matrix η n is diagonal in our case with the diagonal elements and j = 1,2,3 corresponds to the x, y, and z-direction, respectively. As before, L n,j is the depolarization factor of material n in direction j and f n is the volume fraction of material n. (12) ε e,j =ε 1 +ε 1 σ j 1 − L 3,j σ j , η n,jj = 16f n (ε n −ε 1 ) 3(ε 1 + L n,j (ε n −ε 1 )) ,

Modeling wood chips moisture content
The effective medium models above require that the volume fractions of the different materials are known. From the experiments we get the mass of the wood chips in the sample box, m ′ w , for different moisture contents, M d (we omit the index n which refers to moistening cycle). The dry weight of the wood chips in the sample box, m ′ d , can be calculated from M d and m ′ w using Eq. 1. The weight of the water is then m ′ The volume of the sample box containing the wood chips is V sb .
Below the fiber saturation point (which we take to be M d = 30% [19]) there is no free water and the mixture consists of air and wood. The volume fraction of wood becomes where ρ w is the density of the wood, which depends on M d [26] where G w is the specific gravity of wood and ρ H 2 O = 1 g cm −3 is the density of water. The density ρ w is an average of all the wood chips, similarly to the depolarization factors L n,j . The variation in ρ w will be relatively small since all wood chips are of conifer wood. The volume fraction of air is obtained as f 1 = 1 − f 2 below the fiber saturation point (i.e., we set f 3 = 0, in the effective medium model above, since it is the volume fraction of free water).
Above the fiber saturation point, we assume that the wood's volume fraction, f 2 , is the same as at the fiber saturation point (M d = 30%); the mass of the free water, denoted m ′

Wood chips dimensions
The probability density function (pdf ) of the measured dimensions of the 100 measured wood chips was calculated. The longest length is denoted x' , the second length y' , and the shortest z' , as shown in the inset of Fig. 4(a). In Fig. 4(a) the pdf of the wood chips as measured along the longest axis, x' , is shown. The length varies between 10 and 100 mm with a clear peak at 20 mm. For each wood chip the relative dimensions were calculated, and the corresponding pdf:s calculated. In Fig. 4(b) the pdf:s of the ratios y'/x' , z'/x' , and z'/y' are shown. The ratio y'/x' has a clear peak at 0.1, which indicates that the typical wood chip is 10 times longer than it is wide. The ratios z'/x' and z'/y' do not have that clear peaks, but it is clear that the wood chips are elongated.
The two-dimensional autocorrelation of the images (converted to gray scale) in Fig. 1 were calculated. In Fig. 4c the autocorrelation is shown for the top and side views in Fig. 1. In Fig. 4c, x, y and z refer to the coordinate system in Fig. 2. For the top view the x-and z-directions have mainly the same autocorrelation function, which indicates random orientation of the wood chips. For the side view the autocorrelation in the x-direction is clearly larger than in the y-direction, which indicates that the wood chips' are oriented in the x-direction. The autocorrelation length, ξ, is the length at which the autocorrelation function is ξ = 1/e = 0.37. For the top view (17) mm (x) and ξ = 23 mm (z). For the side view ξ = 25 mm (y) and ξ = 103 mm (x). Thus, the side view shows that the wood chips are oriented in the horizontal plane.

Permittivity of wood chips vs moisture
In Fig. 5a the radar signal of a reference measurementwhere the signal has passed through air-is shown with a signal that has passed through wood chips of moisture  To determine κ, the frequency range 0.7 to 2.5 GHz was used where the signal-to-noise ratio was high enough to enable the fitting of a straight line. In Fig. 6a the magnitude vs frequency of the two signals in Fig. 5a are shown. The signal that has passed through wood chips is clearly lower, particularly at frequencies above 2.5 GHz. In Fig. 6b the natural logarithm of the ratio of the signals in Fig. 6(a) is shown vs frequency. Also shown is the fitted straight line used to determine κ. Notice that the signals' levels are low below 0.7 GHz, at 2.1 GHz, and at 2.7 GHz where the antennas gain is low. Data points at 2.1 GHz and above 2.7 GHz were therefore not used in the determination of κ. Figure 7 shows ε ′ vs M d . It increases monotonically with M d for both the directions of the E-field. The values for the horizontal E-field are consistently higher than those for vertical E-field. Fig. 6 a Amplitude vs frequency for the signal in Fig. 5 and b the amplitude ratio of the signals vs frequency used to determine κ Fig. 7 The real part of the permittivity, ε ′ , vs moisture content, M d , for wood chips. Data for the E-field horizontal (x in Fig. 2) and vertical (y in Fig. 2) are shown Fig. 8 The imaginary part of the permittivity, ε ′′ , vs moisture content, M d , for wood chips. Data for the E-field horizontal (x in Fig. 2) and vertical (y in Fig. 2) are shown Fig. 9 Anisotropy ratios k' and k'' of the permittivity for horizontal E-field and for vertical E-field for ε ′ , and ε ′′ , in Figs. 7 and 8, respectively, vs moisture content, M d Figure 8 shows ε ′′ vs M d . It also increases monotonically (except one data point) with M d . As for ε ′ , the values for horizontal E-field are consistently higher than for vertical E-field.
In Fig. 9 the anisotropy ratios of the permittivity, k ′ and k ′′ , are shown for the data in Figs. 7 and 8, respectively. The real part, k ′ , increases monotonically with M d from approximately 1.1 to 1.6. For the imaginary part, k ′′ , the errors in the data are large. Any increase or decrease with M d . cannot be distinguished. However, k ′′ is clearly larger than k ′ and is in the range 2 to 3 at high M d (above 40%).

Effective medium modeling
To model the complex effective permittivity in the x-(horizontal) and y-directions (vertical), we used the models described above. We used literature values for the permittivity. The volume fractions were determined from the measured masses of the wood chips. As free parameters, used to fit the modeled permittivity to the experimental ones, we used the specific gravity, G w in Eq. 16, and depolarization factors, L n,j in the respective model for M d < 30%. For M d > 30% only the depolarization factors of free water was used as a free parameter. (The specific gravity and depolarization factors for wood and air obtained from fitting of data for M d < 30% was used also for M d > 30%.) The modeling is described in more detail in this section.
We use literature values for the permittivity of soft wood from ch. 5 in [27] to calculate ε w , as described above. Figure 10 shows some literature values for the permittivity of wood vs moisture content up to 30% [27,28]. Also shown are the values for ε w that we use in the modeling. We see in the figure that the scattering in the data is relatively large, in particular for ε ′′ w . For ε ′ w the values based on [27] are in the middle of given literature data. For ε ′′ w the values based on [27] are relatively low. We therefore also use higher values for ε ′′ w , also shown in the figure, to investigate the effect of the modeling results from different values of ε ′′ w , since the literature data show such large variations.
For the free water we use ε = 79 + i5 [29] in all the models.
Below the fiber saturation point (for M d < 30%) we model ε as a dielectric mixture of wood and air. We set f 3 = 0 in the BG-MPA, MG-LSI, and MG-MPA models above. The modeling is made in these steps: a. The volume fractions of air, f 1 , and wood, f 2 , are calculated from the experimentally determined masses, as described above. b. For the permittivity of wood, ε w , we used the values in Fig. 10. Thus, we use two values of ε ′′ w , which results in different modeled ε. c. We use the depolarization factors, L n,j , in the respective models to fit the modeled anisotropy ratio of ε ′ , k ′ , to the experimental one in Fig. 9. The anisotropy ratio of ε ′′ , k ′′ , has too large experimental errors to enable a fit (see Fig. 9). We set L n,x = L n,z � = L n,y in all the models, since the wood chips are randomly oriented in the x-z-plane. In the BG-MPA model we set L 1,j = L 2,j for j = x,y, and z. We thus reduce the number of parameters and get the more commonly used version of the model in [22]. d. We use G w to fit the modeled ε ′ to the experimental one in Fig. 7. ε ′ is used since the experimental errors are smaller than for ε ′′ and the scattering in the literature data (Fig. 10) is also smaller for ε ′ than for ε ′′ .
Notice, that for M d < 30% we use only two parameters in the models to fit the experimental data: The specific gravity, G w and depolarization factor L n,x ( L n,y and L n,z are not independent of L n,x ). The experimental values for ε ′′ , are not used in the fit.
Above the fiber saturation point (i.e., for M d > 30%) we model ε as a mixture of air, wood and water: i. The volume fractions of air, f 1 , wood, f 2 , and water, f 3 , are calculated from the experimentally determined masses, as described above. ii. For ε w we use the values in Fig. 10 at M d = 30%. For the free water we use ε = 79 + i5 [29]. iii. The value for G w obtained in the modeling of data for M d < 30% is used. iv. For wood and air we use the depolarization factors as obtained from the modeling for M d < 30%. v. The depolarization factors of water are used to fit the modeled anisotropy ratio, k ′ , to the experimental in Fig. 9. We set L 3,x = L 3,z � = L 3,y as for M d < 30%.
Thus, we use only one free parameter, the depolarization factor of water L 3,x , as free parameters for M d > 30% in addition to G w and L 1,x that were obtained from modeling for M d < 30%.
Technically, the depolarization factors of air and water could be changing with moisture content, since with increasing amount of free water the geometry of the water and air inclusions may change. To vary these

Fig. 11
Volume fraction of air, wood, and water vs moisture content for G w = 0.34

Fig. 12
Permittivity and anisotropy ratios vs moisture content, modeled using the BG-MPA model (Eq. 11) with parameters as in Fig. 10 Fig. 11 the volume fractions of air, f 1 , wood, f 2 , and water, f 3 , vs M d , for G w = 0.34 , used in the BG-MPA model are shown. Up to M d = 30% there is no free water and the volume fraction of wood is close to constant; f 2 changes marginally with M d . Above M d = 30% the amount of water increases and the amount of air decreases; air is replaced by water as the moisture content increases. The curve shapes of f 1 , f 2 , and f 3 , that at not straight lines, are given by the experimental values for M d and m ′ w . The modeled permittivity and anisotropy ratio for the BG-MPA model are shown in Fig. 12. The experimental data in Figs. 7, 8, and 9 are shown for comparison. The ratio k ′ , increases from ca 1 to 1.3 for M d going from 0 to 30%. This behavior is in good agreement with the experimental data in Fig. 9. The ratio, k ′′ , increases from 1.2 to 2.2, but the experimental errors in Fig. 9 make a comparison impossible. For ε ′ the modeled data in Fig. 12 are in good agreement with experimental data in Fig. 7. It increases from ca 1.2 to 1.8 and 2.3 for vertical and horizontal E-fields, respectively. Notice that the modeled ε ′ is not significantly affected by the ε ′′ w used in the model. The high and low ε ′′ w give the same modeled ε ′ . For ε ′′ , the modeled data in Fig. 12 are significantly different for the high and low ε ′′ w in Fig. 10. For both cases the modeled ε ′′ increases monotonically with M d . The ε ′′ values in Fig. 12 for high ε ′′ w are too high by up to 30%; The ε ′′ values for low ε ′′ w are too low by up to 30%, compared with experimental data in Fig. 8. Considering the errors in Fig. 8 the modeled data for ε ′′ are in fair to good agreement when using the high ε ′′ w in Fig. 10. For M d > 30% the modeled anisotropy ratio k ′ , in Fig. 12 increases to 1.5 at M d = 120% in agreement with Fig. 9. The anisotropy ratio, k ′′ , is practically constant at 2.2 for all M d , also in agreement with experimental data in Fig. 9. ε ′ increases up to 4.5 (horizontal) and 3.0 (vertical) at M d = 120%, in qualitative agreement with Fig. 9. For ε ′′ the low ε ′′ w in Fig. 10 gives clearly too low values for M d > 30%. For the high ε ′′ w in Fig. 10 the modeled ε ′′ in Fig. 12 increases up to 1.4 (horizontally) and 0.7 (vertically) at M d = 120% in reasonable agreement with Fig. 8. The different ε ′′ w values used affect in practice only the modeled ε ′′ values. The modeled ε ′ and anisotropy ratios, k ′ and k ′′ , are not significantly affected.
The used value for G w in Table 1 is well in line with data in the literature; for spruce the specific gravity is in the range 0.33 to 0.42 and for pine 0.34 to 0.51 in [30]. The depolarization factors are naturally different from those of isotropic inhomogeneities ( L x = L y = L z = 1 3 ).
Notice that the depolarization factors for the free water ( L 3,x , L 3,y , L 3,z ) in Table 1 are slightly closer to a sphere than those for the wood. The explanation may be that the geometry of the water does not perfectly replicate that of the wood, but instead creates shapes that are less elongated.
Thus, using the BG-MPA model, Eq. 9, we can reproduce the main features of the experimental data in Fig. 7-9. It lends credence to the BG-MPA model, since only three free parameters were used and that the agreements is good also for ε ′′ , which was not used in the fitting. Figure 13 shows the model results for the MG-LSI model. For M d < 30% the anisotropy ratios k ′ and k ′′ are similar to the BG-MPA results and are in good agreement with experimental data in Fig. 9 and 7, respectively. The modeled ε ′ in Fig. 13 is also in good agreement with the experimental data in Fig. 7. For ε ′′ the modeled data using the low ε ′′ w in Fig. 10 give clearly too small values, but the using the high ε ′′ w the agreement with experimental data in Fig. 8 is good.   Fig. 13 Permittivity and anisotropy ratios vs moisture content, modeled using the MG_LSI model (Eq. 12) with parameters as in Fig. 10  Above M d = 30% the modeled data in Fig. 13 for the anisotropy ratios k ′ and k ′′ are too high. Notice in Table 1 that the depolarization factors for wood and free water are the same. The model was very sensitive to any difference between these depolarization factors and gave results that were in poor agreement with experiential data. ε ′ in Fig. 13 shows a behavior similar to the experimental data in Fig. 7. For ε ′′ , the MG-LSI model gives values that decrease with moisture content. The results are not in agreement with the experimental data in Fig. 8. We conclude that the MG-LSI model does not well describe the complex permittivity of wood ships above the fiber saturation point. The explanation is probably that the free water does not cover the wood chips as layers of equal thickness.
The model results for the MG-MPA model are shown in Fig. 14. Up to M d = 30% the results resemble those of the MG-LSI and BG-MPA models, both for the anisotropy ratios k ′ and k ′′ , and for ε ′ . For ε ′′ , however, the results are lower than for the other models and clearly lower than the experimental values, also for the high ε ′′ w values in Fig. 10.
For M d > 30% the results for the MG-MPA model are in better agreement with experimental results for the anisotropy ratios in Fig. 9 than the MG-LSI model and similar to the BG-MPA model (Fig. 12). For ε ′ and ε ′′ the model gives too low values compared to the experimental ones in Figs. 7 and 8. ε ′′ does not change with M d , which is slightly better in agreement with experimental results than the MG-LSI model.
Thus, the MG-LSI and MG-MPA models give reasonable results for data below M d = 30%, but become increasingly poorer with increasing moisture content. The explanation is that these models assume that the material consists of inclusions in a host material (here air), and that this assumption becomes increasingly poorer as the volume fraction of air decreases (cf. Figure 10).

Discussion
We notice that the experimental ε ′ and ε ′′ increase monotonically with moisture content for horizontal and vertical E-field. The experimental data do not show any change in behavior at M d = 30%, in contrast to all the models. The explanation is likely that the limit M d > 30% for free water to appear is not exactly the same for all wood particles (cf. [26]) and that the moisture may not be perfectly homogeneously distributed within the wood chips, though care was taken to get homogeneous moisture distribution.
Experimental data for the permittivity of saw dust in the range 0.5 to 15 GHz increased monotonically with moisture content at all frequencies [7]; the data at 1 GHz are ε = 1.5 + i0.08 at M d = 15%, ε = 2.8 + i0.28 at M d = 41% and ε = 5 + i0.5 at M d = 82%. They report isotropic permittivity. The data are approximately the same as our vertical ε ′ and our horizontal ε ′′ (cf. Figure 7 and 8). Notice that the size of the wood particles in [7] were around 1 mm, i.e., significantly smaller than in our case.
MG-LSI is poor above the fiber saturation point in particular compared to BG-MPA. We suggest that the explanation is that the free water does not follow the shape of layers on the wood chip surfaces. Instead it is less uniformly distributed, which explains why the MG-MPA and in particular BG-MPA are good in agreements with the experimental data.
We notice that in Ref [11] the Maxwell Garnett model was found to model the experimental ε of woody biomass well (The biomass was mixtures of air, wood and bark). The moisture content was not as high as in our case so their results are in agreement with ours, i.e., that Maxwell Garnett can be used to model ε at low moisture content. In [11] they did not include the Bruggeman model, but models that assume that there is a host material (air) in which the inhomogeneities are located.
Our experiments combined moistening of wood chips with radio measurements; the experimental procedure had  Fig. 10 and Table 1. For comparison the experimental data in Figs. 7, 8, and 9 are also shown to be feasible for both types of experiments. The experimental time, 8 h, is shorter than, e.g., the 24 h reported in [31] for conditioning of wood chips. The individual wood chips may not have reached a homogeneous moisture distribution, which may explain some of the discrepancy between modeled and experimental permittivity. We studied industrial wood chips and in industrial applications the moistening process is not well controlled, which may affect industrial radio applications such as [8]. To our knowledge the change in permittivity with time during moistening of wood chips has not been investigated and may be the subject of future research.