## Servicios Personalizados

## Articulo

## Indicadores

- Citado por SciELO

## Links relacionados

- Citado por Google
- Similares en SciELO
- Similares en Google

## Bookmark

## Latin American applied research

*versión impresa* ISSN 0327-0793

### Lat. Am. appl. res. v.38 n.4 Bahía Blanca oct. 2008

**A study of ultrasonic wave propagation in bones**

**P. M. Gauzellino ^{1}, F. I. Zyserman^{1,2} and J. E. Santos^{1,2,3}**

^{1} *Facultad de Ciencias Astromómicas y Geofísicas, UNLP. Paseo del Bosque s/n, B1900FWA -La Plata, Argentina (gauze,zyserman,santos)@fcaglp.unlp.edu.ar*

^{2}

*CONICET*

^{3}

*Department of Mathematics, Purdue Unviersity, 150 N. University Sreet,*

West Lafayette, Indiana, 17907-2067, USA.

West Lafayette, Indiana, 17907-2067, USA.

*Abstract* — In the present work the propagation of ultrasonic waves in human bones is modeled by Biot equations introducing in them viscoelasticity so as to take into account attenuation mechanisms; different models for the latter are introduced and compared. Also, Biot equations are numerically solved using two dimensional finite elements. How different porosities affect elastic moduli, phase velocities, viscodynamic and intrinsic attenuation is also studied.

*Keywords* — Cortical and Trabecular Bone. Low Frequency. Biot Theory. Finite Elements.

**I. INTRODUCTION**

Osseous tissue forms the variety of bones that make up the skeletal system. The matrix is partially organic (collagen fibres in which osteocytes are embedded) and partially inorganic (mineral salts with calcium as the most important mineral). Because of its structure, bone tissue is found in two different types, Compact Bone and Cancellous or Trabecular Bone. Compact bone is very dense and strong, found on the outmost portions of all bones as a protective layer. It contains numerous osteons or Haversian Systems with a central channel through which blood vessels and nerves pass. Surrounding this channel there are multiple concentric layers of tissue known as lamellae. On the other hand, Cancellous Bone is lightweight and spongelike; the interconnected bony pieces forming this characteristic structure are called trabeculae. Trabecular bones have relatively large hollows filled by bone marrow and adipose tissues. Therefore, the relative volume fraction of solid in both structures is different; bone portions with volume fraction of solids below 70% are classified as trabecular, and over 70 % as cortical (Cowin, 1999; Smit *et al.*, 2002; Barkmann *et al.*, 2003; Bossy *et al.*, 2004; Buchanan *et al.*, 2004; Bossy *et al.*, 2005). Biot (1956 a,b,c) established a theory of propagation of elastic waves in a medium composed of a porous elastic solid saturated by a simple and viscous fluid phase. He showed the existence of a shear wave and two compressional waves, a high velocity one (wave of the first kind or Type I wave) corresponding to in-phase motion of solid and fluid and a low velocity one (wave of the second kind or Type II wave) associated to out of phase motion between solid and fluid. These compressional waves had been observed in laboratory experiments with bones, see for example Lakes *et al.* (1986); McKelvie and Palmer (1991); Kacsmarek *et al.* (2002); Lee *et al.* (2002); Cardoso *et al.* (2003).

Numerical simulation of ultrasonic wave propagation in bones using Biot theory has been reported since several years ago by many authors -see Haire and Langton (1999); Hughes *et al.* (2003); Lee *et al.* (2003); Buchanan *et al.* (2004); Fella *et al.* (2004); Wear *et al*. (2005) to name some of them-, but still there are open questions to be answered, as the reader can see below. The aim of this report is to analyze the propagation of ultrasonic waves in human bones and study how different porosities affect elastic moduli, velocities and intrinsic attenuation. Biot theory is used for two dimensional numerical simulation where bones are considered as a biphasic poroviscoelastic material. It is important to remark that the upper bound for the involved frequencies is reached when the corresponding wavelength becomes of the order of the porous size. Experiments involving frequencies beyond this limit must be treated within a different theoretical frame, as stated by the author of the theory (Biot, 1956c).

Following Mellish *et al.* (1989) and Hughes *et al.* (2003), the typical pore radii (mean values and standard deviations) are 0.285 ± 0.050 mm for young normal bone and 0.455 ± 0.130 mm for aged and osteoporotic bone. Bossy *et al*. (2005) reports an average trabecular spacing ranging from 0.5 mm to 2.0 mm; therefore Biot theory cannot be used with bones for frequencies greater than 750 kHz.

Keeping in mind that Biot's theory is applicable at wavelengths much longer than the pore size (corresponding to frequencies within the 0–750 kHz range), and that the Type I and shear waves have a behaviour similar to those in an elastic solid, with high phase velocities, low attenuation and little dispersion, it is necessary to establish the critical frequency *f*_{c} that determines whether or not Poiseuille flow occurs. When the frequency range is below *f*_{c}, the relative motion of the fluid in the porous is of the Poiseuille type; the coupling between solid and fluid is preeminently viscous, leading to a slow wave of diffusive nature. The assumption of Poiseuille is not established for the higher frequency range *f*_{c} ≤*f*; in this regime inertial coupling dominates over the viscous one leading to a truly propagating slow compressional wave. Johnson and Plona (1982) showed that the latter cannot propagate if the viscous skin depth is less than the porous radius and developed a second frequency limit called the viscous frequency, *f*_{v}. The viscous skin depth depends on the frequency and the fluid properties (viscosity and density) as the next section will make clear. Results obtained by Hughes *et al.* (2003) indicate mean values of *f*_{c}=4.5 Hz and *f*_{c}=705 Hz for normal bone with water or marrow at 20^{o}C, respectively. For osteoporotic bone these same values are now *f*_{c}=1.4 Hz and *f*_{c}=231 Hz. Therefore, the viscodynamic correction must be introduced in the ultrasonic range, as has been reported in Fella *et al.* (2004). In the present paper, the theory of dynamic permeability and tortuosity presented by Johnson *et al.* (1987) is considered for calculating this frequency correction function.

Other point that is taken into account in the present work is that, when numerically studying how the wave propagation depends on the values of different model parameters, their variation must be done in a consistent fashion. For example, if the bone porosity is diminished, the bulk moduli must be accordingly modified, otherwise phase velocities will be under or overestimated. The Hashin-Shtrikman bounds (Mavko *et al.*, 1998) are used in this work to estimate the dependence of the moduli with porosity.

**II. BIOT THEORY**

Consider an isotropic, porous solid saturated by a single-phase, compressible and viscous fluid. Let and , *i* =1, · ,*d* be the averaged displacement vectors of the solid and fluid phases, respectively, where *d* indicates the Eucliden dimension, i.e., *d* =2, 3. Also, let

(1) |

be the average relative fluid displacement per unit volume of bulk material, where *Φ* means the effective porosity. Set *u* =(*u*^{s},*u*^{f} ) and recall that

(2) |

represents the change in fluid content.

Let *ε*_{ij}(*u*^{s}) be the strain tensor of the solid. Also, let *φ*_{ij}, *i,j* =1, ·· ,*d*, and *p*_{f} denote the stress tensor of the bulk material and the fluid pressure, respectively. Following (Biot, 1962), the stress-strain relations can be written in the form:

(3) |

The coefficient *µ* is equal to the shear modulus of the bulk material, considered to be equal to the shear modulus of the dry matrix. Also

(4) |

with *K*_{c} being the bulk modulus of the saturated material. Following Santos *et al.* (1992), Gassmann (1951) the coefficients in (3) can be obtained from the relations

(5) |

where *K*_{s}, *K*_{m} and *K*_{f} denote the bulk modulus of the solid grains composing the solid matrix, the dry matrix and the saturant fluid, respectively. The coefficient α is known as the effective stress coefficient of the bulk material.

**A. The modified moduli for viscoelasticity**

It is well established that frequency-dependent propagation velocity and attenuation are the most significative parameters to determine bone quality by ultrasound (Siffert and Kaufman, 2007). Since experimental measurements yield greater attenuation than the predicted by the lossless Biot theory, an attempt is here made to remedy this misbehaviour by introducing attenuation models into the latter. The chosen ones are popular within the geophysics community, and assume that the material is linear, dissipative and causal. The modified moduli *K*_{m}, *K*_{s} and *µ* are introduced by application of Biot's correspondence principle (Biot, 1956c; 1962). The poroelastic coefficients of the stress-strain relations are changed by new complex coefficients that depend on the frequency and satisfy the same relations that are satisfied by the elastic moduli. It is proposed to parametrize attenuation in bones through constant quality factor models, since there exist measurements that indicate an almost linear dependence of attenuation with frequency and the quality factor is nearly constant over a given frequency band. The following formulae allow to calculate the viscoelastic moduli for different models. The form of the moduli in the frequency-domain is obtained by Fourier transform of the stress-stain relations. For a detailed description of the differences among the models here presented, see Mavko and Nur (1979); Johnston *et al*. (1979).

**B. Liu model**

In the linear viscoelastic model proposed by Liu *et al.* (1976) the following formula allows to calculate the viscoelastic moduli:

(6) |

Here *M* represents any of the mentioned moduli and the coefficient *M*_{re} corresponds to associated elastic modulus (Bourbie *et al.*, 1987). The coefficients *R*_{M}(*ω*) and *T*_{M}(*ω*), associated with a continuous distribution of relaxation times, characterize the viscoelastic behaviour of the material and are given by (Liu *et al.*, 1976)

(7) | |

(8) |

Another measurement of the viscoelastic behaviour of the materials is

(9) |

being *Q*(*ω*) the so-called quality factor.

*T*_{1,M} and *T*_{2,M} are chosen such that the quality factor *Q*(*ω*) is approximately equal to the average value in the range *T*_{1,M} ≤*ω*≤*T*_{2,M}. In general, small values of (on the order of ten) correspond to very dissipative materials, and values of on the order of hundred correspond to quasielastic materials.

An advantage of this model over the others here considered is that the complex modulus behaves appropriately within a wide frequency range including the quasiestatic regime, i.e. *ω*→ 0.

**C. White model**

This approach introduces a perfect constant *Q* for all frequencies (White, 1965). It is a simple model where the complex moduli are calculated by

(10) |

**D. Zener or standard linear solid model**

This model, usually mentioned as SLS, is suitable to represent a distribution of single relaxation peaks centered at each frequency. The viscoelastic moduli are obtained by

(11) |

Also, it is possible consider only a single standard linear solid centered at a given frequency. In this case, the quality factor *Q* is frequency dependent.

Recall that in all cases the complex moduli are calculated by applying the different models to *K*_{s}, *K*_{m} and *µ*.

**E. The equations of motion**

The equations of motion for a porous medium saturated by a single-phase in the space-frequency domain can be written, following Santos *et al.* (1992), as

(12) | |

(13) |

where *ρ*_{s} and *ρ*_{f} are the mass densities for solid and fluid, respectively and *f*^{s} and *f*^{f} are solid and fluid sources respectively. The total mass density of the solid-fluid or bulk density is

(14) |

*g* is the mass coupling coefficient between the solid matrix and the fluid phase. This coefficient represents inertial effects. The parameter *b* is the viscosity coupling coefficient between both phases. The relationship among *b*, the fluid viscosity *η*, and the Darcy's coefficient of permeability *k* is given by

(15) |

while *g* is related with the tortuosity factor, *S* as

(16) |

It was mentioned that the assumption of Poiseuille flow breaks down when the frequencies exceed a value *ω*_{c} (Biot, 1956c; Johnson *et al.*, 1987; Carcione, 2001). This value can be accurately estimated according to the formula

(17) |

where *a*_{p} = 2(5*k*/*Φ*)^{1/2} is the pore size parameter as function of the permeability *k* and the porosity *Φ* (Johnson, 1982; Bear, 1972; Scheidegger, 1974; Hovem and Ingram, 1979). For frecuencies higher than *ω*_{c}, the coupling coefficients depend on the frequency as

(18) | |

(19) |

The function *F*(*ω*)=F_{r}(*ω*)+ *i F*_{i}(*ω*) is a correction factor proposed by Biot (1956c) as an universal function representing effects of frequency for different size of the pores and their geometry. Johnson *et al.* (1987) provided the following modified version that is used in this article,

(20) |

and ζ can be calculated from

(21) |

The parameter ζ has the dimensions of length and it is a geometrical parameter that represent a characteristic of the porous medium.

A second limit called viscous frequency, *ω*_{v}, was suggested by Johnson and Plona (1982) as

(22) |

being *ζ* a scaling constant with an empirically derived value of the order of 0.01.

In order to completely define the problem to be studied, boundary conditions must be attached to equations (12-13). Considering from now on a two dimensional medium, let Γ denote the frontier of the domain Ω, with *v* being the unit outer normal at the boundary and its unit tangent such that {*v*,} is an orthonormal system on Γ. The absorbing boundary condition chosen is

(23) |

The matrix D is positive definite and it is given by , where , and

(24) |

**F. Phase velocity and attenuation**

Plane waves can be used for obtaining the characteristics of waves propagating in a medium. Because the medium is isotropic, without loss of generality the plane waves propagating in the *x*-direction are

(25) |

where is the complex wave vector. Substituting these expressions into equations of motion (12-13) and zeroing the right-hand yields

(26) |

(27) |

Taking the determinant of this system of equations equal to zero and considering the complex velocity *v*_{c}=*ω*/*k*_{c}, the following dispersion relation is obtained:

(28) |

where *H*_{c}= λ_{c} +2*µ*. The two roots of this second-order equation in *v*_{c}^{2} correspond to the fast and slow compressional waves. The phase velocity and the attenuation factor can be calculated by

(29) |

respectively.

**G. The effective moduli and compressional velocities**

In order to theoretically predict the effective elastic moduli, *K*_{m} and *µ*=*µ*_{m} of a mixture of a porous solid material and a fluid, generally, it is necessary to specify the volume fractions and the moduli of the constituents. The Hashin-Shtrikman expressions allow to predict the upper and lower bounds of the elastic bulk and shear moduli of the dry matrix without being forced to consider details of the geometries of the phases (Hashin and Shtrikman, 1963; Mavko *et al.*, 1998). A similar approach to calculate upper and lower bounds for other Biot parameters can be found in Norris (1989). For the upper bounds the Hashin-Shtrikman formulae are

(30) | |

(31) |

The lower bounds are calculated by interchanging moduli and volume fractions; the separation between both bounds depends on how different the phases are. Because Biot's effective medium model consider a homogeneous modulus, a possibility is to take the average of the upper and lower bounds, and . Notice that there exist other ways of calculating these bounds, e.g. Gassmann's relations, but they are not used here because they are appropriated for quasi-static regimes. **Figure 1** shows the separation between the upper and average bounds vs. porosity assuming *K*_{s} = 16.74 *GPa*, *µ*_{s} =6.25 *GPa* for bone (Smit *et al.*, 2002) and *K*_{air} =0.01 *GPa*, *µ*_{air} =0 *GPa* for the air. Therefore, the lower bound for *K*_{m} is in the range 0.01 *GPa* to 0.25 *GPa* and for *µ*_{m} is zero. The fast and slow compressional velocities can be estimated from the above results. Several considerations concerning the pore size, the porosity, the permeability, the tortuosity and some geometric details in addition to previous elastic moduli have been taken into account to plot velocities vs. porosity in Fig. 2. This figure illustrates the behavior of elastic and anelastic velocities for different viscoelastic models with the quality factor = 30 if the bone porous matrix is saturated by water. Note that both type I and type II wave velocities tend to the speed of sound in water (1563 ± 25 m/s), Nicholson *et al.* (2000) as the porosity increases (lower solid fraction). All Biot viscoelastic velocities are greater than elastic ones and the major difference is about 10 % for the Liu model. The behaviour of the attenuation coefficient as a function of the porosity is shown in Fig. 3. While the attenuation of the fast wave increases slowly when the solid fraction decreases, the attenuation of the slow wave presents high values for cortical bone porosities. **Figure 4** represents the attenuation coefficient as a function of the frequency for trabecular and cortical structures with *Φ* =0.8 and *Φ* =0.05, respectively. The curves correspond to = 30, *T*_{1}_{,M} =1 *s* and *T*_{2}_{,M} =1.10^{-12} *s*. It must be noticed that in these models the type of dispersion is anomalous, i.e., the phase velocity increases with frequency. In addition, it is possible to analyze how the attenuation of fast and slow waves depends on . Considering a cancellous bone, Fig. 5 displays curves of attenuation vs. frequency for = 30 and = 10(Liu model). The attenuation is notably stronger for the fast wave; and for this model these results are close to the experimental results of Serpe and Rho (1996).

**Figure 1**: Effective moduli (measured in Pa) of the porous matrix of dry bone. Bottom: Bulk modulus. Top: Shear modulus.

**Figure 2**: Compressional Biot velocities as a function of the porosity for different attenuation models. The single SLS is centered at f0 = 300kHz. Bottom: Type I. Top: Type II.

**Figure 3**: Attenuation coefficients, measured in dB/cm, as a function of the porosity and different attenuation models for Type I wave (bottom) and Type II wave (top). Again, the single SLS is centered at f0 = 300 kHz.

**Figure 4**: Attenuation coefficients, measured in dB/cm, for Type I and Type II waves in trabecular and cortical bones as a function of the frequency, for different attenuation models. As in previous figures, the single SLS is centered at f0 = 300 kHz.

**Figure 5**: Attenuation coefficients (dB/cm) as a function of the frequency for Type I wave (bottom) and Type II wave (top). Liu model with = 30 and = 10.

**III. NUMERICAL APPROXIMATION**

Naturally, the space–frequency domain is chosen to describe phenomena of dissipation and anelastic attenuation suffered by the different type of waves here considered. Moreover, this domain allows to compute the solution at the current time without the knowledge of the time history of the system, that is, the solution is computed for each temporal frequency and a finite number of them are required to obtain the solution in the space-time domain via inverse Fourier transform (Douglas *et al.*, 1994; 1993; Gauzellino *et al.*, 2001).

As already mentioned, the numerical resolution of Eq. (12-13-23) is performed using piecewise linear finite elements.

**A. Wave propagation modeling**

A numerical experiment to analyze the influence of porosity variations in the observed traces within a piece of bone is designed. The following simple source is chosen:

(32) |

being *f*_{0} = 300kHz the principal frequency of the source, with *ξ* =8 *f*_{0} and t_{0} =1.25/*f*_{0}. A finite number of frequencies between 0 kHz and 600 kHz is used in order to approximate *u*^{s}(*x*, *t*) and *u*^{f}(*x*, *t*). Table 1 shows estimated parameters at different porosities to be used in the simulation of wave propagation. The fluid phase is water and its parameters *K* =2.3*GPa*, *ρ* =1*gcm*^{-3} and *η* =0.01*Poise* are taken from (McKelvie and Palmer, 1991; Buchanan *et al*., 2004). The bone specimen is represented by a square sample with sides of 30 mm. The signal transmitted through the bone is recorded by a receiver located at 30 mm from the source. **Figure 6** illustrates the signals acquired, *u*^{s}(1−*Φ*)+ *u*^{f} *Φ* for different porosities. It can be seen that the slow wave amplitude is relatively important in cancellous structures contrarywise to what happens in a compact structure. In trabecular bone, the velocities of the fast and slow waves are similar due to the relatively high percentage of fluid. Therefore, when the porosity is very high, what is recorded is a superposition of the two waves. It has been already established that osteoporotic bone can be also characterized by low density. For the fast wave, if the density decreases to 20 %, the velocity increases about 8 % and the amplitude increases too; on the contrary, for the slow wave the changes are not significant.

Table 1: Estimated parameters for cortical and trabecular bone. CB: Cortical bone, NTB: Normal trabecular bone, OTB: Osteoporotic trabecular bone.

**Figure 6**: Signals acquired, *u*^{s}(1−*Φ*)+ *u*^{f} *Φ*, for samples with different porosities; the bottom one corresponds to a cortical bone, the top one to different trabecular bones.

**IV. CONCLUSIONS**

It is long established that the Biot theory for wave propagation in a fluid saturated porous medium may be a useful tool to model ultrasound probing of both trabecular and cortical bone. In this paper an upper bound for the frequency spectrum of the acoustic source has been introduced; this upper bound is far lower than frequencies used in the literature assuming the validity of Biot theory.

An open question is which underlying loss process within the bone leads to wave attenuation. In the present effort this mechanism is accounted for assuming viscoelastic behaviour, modeled by several different mechanisms. The propagation of fast and slow waves travelling through both cancellous and cortical structures is fairly well followed by all implemented possibilites.

Is is also observed that the measured signal, when the target is trabecular bone is an overlap of both slow and fast waves, rendering meaningless the idea of extracting information by analyzing just one of the two involved propagation modes.

**ACKNOWLEDGEMENTS**

This work was funded in part by grants PIP 5126/04 (CONICET) and PICT 03-13376 (ANPCyT)

**REFERENCES**

1. Barkmann, R., W. Timm, R. Scheffzcyk, E. Lochmüller and C. Glüer, "Fast wave propagation in the epiphysis of human finger phalanges. Associations between ultrasound variables and structure of trabecular bone", *Proc. WCU*, París, 1173-1176 (2003). [ Links ]

2. Bear, J., *Dynamics of fluids in porous media*. Dover Publications, New York (1972). [ Links ]

3. Biot, M.A., "Theory of deformation of a porous viscoelastic anisotropic solid", *J. Appl. Phys.*, **27**, 459- 467 (1956a). [ Links ]

4. Biot, M.A., "Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low frequency range", *Journal of the Acoustical Society of America*, **28**, 168-178 (1956b). [ Links ]

5. Biot, M.A., "Theory of propagation of elastic waves in a fluid-saturated porous solid. II. High frequency range", *Journal of the Acoustical Society of America*, **28**, 179-191 (1956c). [ Links ]

6. Biot, M.A., "Mechanics of deformation and acoustic propagation in porous media", *J. Appl. Phys*., **33**, 1482-1498 (1962). [ Links ]

7. Bossy, E., M. Talmant and P. Laugier, "Three-dimensional simulations of ultrasonic axial transmission velocity measurement on cortical bone models", *Journal of the Acoustical Society of America*, **115**, 2314-2324 (2004). [ Links ]

8. Bossy, E., F. Padilla, F. Peyrin and P. Laugier, "Three-dimensional simulation of ultrasound propagation through trabecular bone structures measured by synchrotron microtomography", *Physics in Medicine and Biology*, **50**, 5545-5556 (2005). [ Links ]

9. Bourbie, T., O. Coussy and B. Zinszner, *Acoustics of Porous Media,* Editions Technip, Paris (1987). [ Links ]

10. Buchanan, J., R. Gilbert and K. Khashanah, "Determination of the parameters of cancellous bone using low frequency acoustic measurements", *Journal of Computational Acoustics*, **2**, 99-126 (2004). [ Links ]

11. Carcione, J.M., "Wave fields in real media: Wave propagation in anisotropic, anelastic and porous media", *Handbook of Geophysical Exploration*, **31** (2001). [ Links ]

12. Cardoso, L., A. Meunier and C. Oddou, "Theoretical and experimental ultrasonic characterization of anisotropic properties in cancellous bone", *Proc. WCU*, París, 1169-1172 (2003). [ Links ]

13. Cowin, S., "Bone poroelasticity", *Journal of Biomechanics*, **32**, 217-238 (1999). [ Links ]

14. Douglas, Jr. J., J.E. Santos, D. Sheen and L. Bennethum, "Frequency domain treatment of one- dimensional scalar waves", *Math. Models Meth. Appl. Sci.*, **3**, 171-194 (1993). [ Links ]

15. Douglas, Jr. J., J.E. Santos and D. Sheen, "Approximation of scalar waves in the space-frequency domain", *Math. Models Meth. Appl. Sci*., **4**, 509-531, (1994). [ Links ]

16. Fella, Z.E.A., J. Chapelon, S. Berger, W. Lauriks and C. Depollier, "Ultrasonic wave propagation in human cancellous bone: Application of biot theory", *Journal of the Acoustical Society of America*, **116,** 61-73 (2004). [ Links ]

17. Gassmann, F., "Über die elastizität poröser medien", *Vierteljahrsschrift der Naturforschenden Gessellschaft,* Zurich, **96**, 1-23 (1951). [ Links ]

18. Gauzellino, P.M., J.E. Santos and D. Sheen, "Frequency domain wave propagation modelling in exploration seismology", *Journal of Computational Acoustics*, **9**, 941-955 (2001). [ Links ]

19. Haire, T.J. and C.M. Langton, "Biot theory: A review of its application to ultrasound propagation through cancellous bone", *Bone*, **24**, 291-295 (1999). [ Links ]

20. Hashin, Z. and S. Shtrikman, "A variational approach to the elastic behavior of multiphase materials", *Journal Mech. Phys. Solids*, **11**, 127-140 (1963). [ Links ]

21. Hovem, J.M. and G.D. Ingram, "Viscous attenuation of sound in saturated sand", *J. Acoust. Soc. Amer.*, **66**, 1807-1812 (1979). [ Links ]

22. Hughes, E., T. Leighton, G. Petley, P. White and R. Chivers, "Estimation of critical and viscous frequencies for biot theory in cancellous bone", *Ultrasonics*, **41**, 365-368 (2003). [ Links ]

23. Johnson, D.L., "Elastodynamics of gels", *Journal of Chemical Physics*, **77**, 1531-1539 (1982). [ Links ]

24. Johnson, D.L. and T.J. Plona, "Acoustic slow waves and the consolidation transition", *J. Acoust. Soc. Am*., **72**, 556-565 (1982). [ Links ]

25. Johnson, D.L., J. Koplik and R. Dashen, "Theory of dynamic permeability and tortuosity in fluid-saturated porous media", *Journal of Fluid Mechanics*, **176**, 379- 402 (1987). [ Links ]

26. Johnston, D.H., M.N. Toksösz and A. Timur, "Attenuation of seismic waves: Part II", *Geophysics*, **44,** 691- 711 (1979). [ Links ]

27. Kacsmarek, M., J. Kubik and M. Pakula, "Short ultrasonic waves in cancellous bone", *Ultrasonics*, **40**, 95-100 (2002). [ Links ]

28. Lakes, R., H. Yoon and L. Katz, "Ultrasonic wave propagation and attenuation in wet bone", *J. Biomed. Eng*., **8**, 143-148 (1986). [ Links ]

29. Lee, T., R. Lakes and A. Lal, "Investigation of bovine bone by resonant ultrasound spectroscopy and transmission ultrasound", *Biomech. and Modeling in Mechanobiology*, **1**, 165-175 (2002). [ Links ]

30. Lee, K.I., H. Roh and S.W. Yoon, "Acoustic wave propagation in bovine cancellous bone: application of the modified Biot-Attemborough model", *Journal of the Acoustical Society of America*, **114**, 2284-2293 (2003). [ Links ]

31. Liu, H.P., D.L. Anderson and H. Kanamori, "Velocity dispersion due to anelasticity; implications for seismology and mantle composition", *Geophys. J. R. Astr. Soc.*, **l47**, 41-58 (1976). [ Links ]

32. Mavko, G. and A. Nur, "Wave attenuation in partially saturated rocks", *Geophysics*, **44**, 161-178, (1979). [ Links ]

33. Mavko, G., T. Mikerji and J. Dvorkin, *The Rock Physics Handbook, tools for seismic analysis in porous media*, Cambdridge University Press (1998). [ Links ]

34. McKelvie, M. and S. Palmer, "The interaction of ultrasound with cancellous bone", *Physics in Medicine and Biology*, **36,** 1331-1340 (1991). [ Links ]

35. Mellish, R., N. Garrahan and J. Compston, "Age-related changes in trabecular width and spacing in human iliac crest biopsies", *J. Bone Miner.*, **6**, 331-338 (1989). [ Links ]

36. Nicholson, P., R. Strelitzki, R. Cleveland and M. Bouxsein, "Scattering of ultrasound in cancellous bone: predictions from a theorical model", *Journal of Biomechanics*, **33**, 503-506 (2000). [ Links ]

37. Norris, A.N., "Stoneley wave attenuation and dispersion in permeable formations", *Geophysics*, **54**, 330-341 (1989). [ Links ]

38. Santos, J.E., J.M. Corberó, C.L. Ravazzoli and J.L. Hensley, "Reflection and transmission coefficients in fluid-saturated porous media", *Journal of the Acoustical Society of America*, **91**, 1911-1923, (1992). [ Links ]

39. Scheidegger, A.E., "The physics of flow through porous media", University of Toronto (1974). [ Links ]

40. Serpe, L. and J. Rho, "The nonlinear transition period of broadband ultrasound attenuation as bone density varies", *Journal of Biomechanics*, **29**, 963-966 (1996). [ Links ]

41. Siffert, R.S. and J.J. Kaufman, "Ultrasonic bone assessment: The time has come", *Bone*, **40**, 5-8 (2007). [ Links ]

42. Smit, T., J. Huyghe and S. Cowin, "Estimation of the poroelastic parameters of cortical bone", *Journal of Biomechanics*, **35**, 829-835 (2002). [ Links ]

43. Wear, K., A. Laib, A. Stuber and J. Reynolds, "Comparison of measurements of phase velocity in human calcaneus to Biot theory", *Journal of the Acoustical Society of America*, **117**, 3319-3324 (2005). [ Links ]

44. White, J.E., *Seismic Waves, Radiation, Transmission and Attenuation*, McGraw-Hill (1965). [ Links ]

**Received: August 27, 2007. Accepted: September 20, 2008. Recommended by Subject Editor: Alberto Cutiño**