- Research
- Open access
- Published:

# Influence of cancellous bone microstructure on ultrasonic attenuation: a theoretical prediction

*BioMedical Engineering OnLine*
**volume 18**, Article number: 103 (2019)

## Abstract

### Background

Quantitative ultrasound has been used for the assessment of cancellous bone status. The attenuation mechanisms of cancellous bone, however, have not been well understood, because the microstructure of cancellous bone is significantly inhomogeneous and the interaction between ultrasound and the microstructure of cancellous bone is complex. In this study, a theoretical approach was applied to investigate the influence of the microstructure of cancellous bone on ultrasonic attenuation.

### Results

The scattering from a trabecular cylinder was significantly angle dependent. The dependencies of the ultrasonic attenuation on frequency, scatterer size, and porosity were explored from the theoretical calculation. Prediction results showed that the ultrasonic attenuation increased with the increase of frequency and decreased linearly with the increase in porosity, and the broadband ultrasound attenuation decreased with the increase in porosity. All these predicted trends were consistent with published experimental data. In addition, our model successfully explained the principle of broadband ultrasound attenuation measurement (i.e., the attenuation over the frequency range 0.3–0.65 MHz was approximately linearly proportional to frequency) by considering the contributions of scattering and absorption to attenuation.

### Conclusion

The proposed theoretical model may be a potentially valuable tool for understanding the interaction of ultrasound with cancellous bone.

## Background

As a major public health problem, osteoporosis is a metabolic skeletal disease where the risk of fragility fracture increases because the bone becomes porous and brittle [1]; it is the most common reason for a broken bone among the elderly [2]. Dual-energy X-ray absorptiometry is considered the gold standard for bone status assessments [3]; however, quantitative ultrasound may be a favorable tool for community-based screening of the high-risk population because it is a cheap, easy-to-access and non-ionizing method [4,5,6]. One of the quantitative ultrasound parameters is ultrasonic attenuation, which has been reported both in vitro and in vivo [4, 7, 8]. Bone mineral reduction is the characteristic of osteoporosis, and previous research indicated that the bone mineral content was associated with ultrasonic attenuation [4, 9]. Langton and colleagues [4] found a quasi-linear dependence of attenuation on frequency on human cancellous bone over the approximate frequency range of 0.3–0.65 MHz; based on this finding, the broadband ultrasound attenuation, which is the slope of the ultrasonic attenuation coefficient with respect to frequency, has received great attention to assess the bone status [4, 10]. Although the broadband ultrasound attenuation has been successfully applied to evaluate the bone status, the attenuation mechanisms of ultrasound propagation in bone are still not well understood, especially for cancellous bone.

The overall attenuation is the combination of absorption and scattering. Cancellous bone is a porous composite material composed of marrow and trabecular networks [11, 12]. Marrow is a fluid-like viscous medium; the ultrasound energy will be absorbed due to its viscous material properties. Compared with marrow, trabeculae have much higher acoustic impedance (\(7.9 \times 10^{6}\) vs \(1.3 \times 10^{6} {\text{kg}}/\left( {{\text{s}}\;{\text{m}}^{2} } \right)\) for common values); moreover, the complex microstructure of cancellous bone, such as the size of trabeculae, bone volume fraction, and orientations of trabecular network, significantly influences the ultrasonic propagation in cancellous bone [13]. The large mismatch of acoustic impedance and complex microstructure of cancellous bone may lead to strong scattering of an incident ultrasonic wave. Therefore, both the absorption and scattering may play important roles as ultrasonic waves propagate in cancellous bone; however, few theoretical models consider energy loss from both scattering and absorption.

Biot’s theory [14, 15], which predicts that absorption is due to viscous losses at internal interfaces, has been applied to cancellous bone [16]. Numerous parameters are required to be known in this theory, and generally these parameters are difficult to be acquired, which limits its application in cancellous bone. In addition, this theory underestimates the experimentally observed attenuation in cancellous bone [17]. Meanwhile, some other scholars suggested that scattering might be a crucial contribution of attenuation for cancellous bone besides the absorption. Chaffai et al. [18] used the weak scattering model, which was based on a random description of the scattering tissue. Strelitzki et al. [19] adopted the binary mixture model, which was based on the assumption that the scattering was proportional to the mean fluctuations of velocity and density. The density fluctuation was ignored in their binary mixture model for cancellous bone. Although the weak scattering model and the binary mixture model have been successfully applied to model scattering from soft tissue, they require isotropy and small contrasts in acoustic properties between the solid trabeculae and the fluid marrow, which is not satisfied because of the large mismatch of acoustic impedance between trabeculae and marrow. Wear [20] proposed the Faran cylinder model, which assumed that trabecular cylinders with diameters much smaller than incident wavelength were the dominant source of scattering in bone. He further suggested that absorption was likely to be significant component of attenuation than scattering because of the discrepancy of attenuation at diagnostic frequencies between the theoretical prediction and experiment data. Thus, to well understand the attenuation mechanisms of cancellous bone, it is necessary to consider the energy loss from both the absorption and scattering attenuation mechanisms. The objective of this article was to investigate the dependence of ultrasonic attenuation on the cancellous bone microstructure by considering both absorption and scattering from a theoretical prediction.

## Results

Figure 1a, b shows typical scattering patterns of the scattered wave from a trabecular cylinder at wavelength \(\lambda = 8\pi a\) and \(\lambda = 4\pi a\). The scattering was significantly angle dependent; the wave was scattered almost uniformly in the backward directions, and the intensity amplitude of scattered wave was much smaller in the forward direction than that in the backward direction. Figure 1c demonstrates the normalized quantity \(\gamma^{\text{sca}}\), which represents the total scattered intensity, as a function of \(ka = 2\pi a/\lambda\). The total intensity scattered varied nonlinearly with \(ka\).

Figure 2 shows the attenuation coefficients as a function of frequency for different trabecular sizes at the porosity of 80%. Six experiment data from literature were also plotted in the figure. It was noted that the attenuation increased with the increase in frequency, and this trend was consistent with previous experimental reports [7, 10, 21,22,23]. Theoretical calculation also revealed that the attenuation was higher for larger trabecular size.

Figure 3 demonstrates the attenuations due to absorption and scattering over the frequency range of 0.3–0.65 MHz at a trabecular size of 100 µm and a porosity of 80%. Both attenuation due to absorption and scattering increased as the increasing of frequency. The attenuation was greater due to absorption than due to scattering below the frequency of 0.6 MHz; however, the attenuation due to scattering, which grew in the form of a power law function, increased faster than that due to absorption, which grew approximately linearly. It was noted that the power exponent was found to be 2.83 (close to 3.0) by fitting the attenuation coefficient due to scattering with a power law function. Although the attenuation from scattering was strongly non-linear, the overall attenuation approximately linearly increased with the frequency, which might be mainly caused by the compensation of the attenuation from absorption.

The theoretical attenuation coefficients versus porosity at different trabecular sizes are shown in Fig. 4. The ultrasound attenuation coefficient decreased linearly with the increase in porosity, and a similar trend was observed from Lee’s experiment data [10]. The broadband ultrasound attenuation coefficient at different trabecular sizes is illustrated in Fig. 5. It was observed that the broadband ultrasound attenuation coefficient decreased with the increase in porosity from the theoretical calculation, and a similar trend was found from Hdgskinson’s and Lee’s experimental studies [8, 10].

## Discussion

In this study, a theoretical model was applied to investigate the influence of microstructure of cancellous bone on ultrasonic attenuation by considering attenuation from both scattering and absorption. This model did not fully account for the trabecular microstructure; however, it successfully predicted the trends between ultrasonic attenuation coefficient and trabecular size, bone volume fraction, and frequency. The dependence of broadband ultrasound attenuation on porosity was also successfully revealed.

Experimental studies have shown that the attenuation coefficient is a quasi-linear function of frequency over the typical range of 0.3–0.65 MHz [4, 24, 25], and this quasi-linear assumption is the basis of broadband ultrasound attenuation measurement. From the perspective of scattering, however, the quasi-linear assumption is difficult to be explained. As in the low-frequency limit, the scattered intensity from a cylinder is proportional to the frequency cubed [26]. Strong non-linear scattering would occur. Wear [20] investigated the frequency dependence of ultrasonic backscatter from human cancellous bone. He observed the discrepancy of attenuation at diagnostic frequencies between the theoretical prediction and experimental data, and then proposed that the absorption was likely to contribute more to attenuation than scattering in the low-frequency range. In our calculation, we found that attenuation due to scattering was well fitted by a power law function with a power index very close to 3. Although the scattering portion of the intensity was almost frequency-cubed dependent, we found that the attenuation due to absorption portion was comparable to the attenuation due to scattering and was more linear in trend. The non-linear trend of attenuation due to scattering was compensated by the introduction of absorption, which caused the quasi-linear behavior of the overall attenuation. Good agreement between the theoretical prediction and linear fitting was observed in the frequency range of 0.3 to 0.65 MHz. Thus, the approach presented here theoretically validates the quasi-linear assumption of attenuation coefficient for cancellous bone over a typical frequency range of 0.3–0.65 MHz.

There are several limitations in our approach. First, the trabeculae were assumed to be cylinders with uniform radius in the model, and the geometrical heterogeneity of trabecular structure was not fully accounted for. This simple assumption may influence the accuracy of our theoretical predictions and further hinder the clinical application. The real shape, dimension, and orientation of trabeculae are highly complex. Future studies are required to obtain a more exact prediction of attenuation from the cancellous bone using more realistic bone microstructure. Cancellous bone may be considered as a composite material composed of matrix and inclusions. In the field of micromechanics, the ellipsoidal is used to define the inclusion shape and the orientation distribution function is used to define the orientation of inclusion in the matrix [27, 28], which may be applied to characterize the microstructure of cancellous bone. Second, multiple scattering was ignored in this study and only single scattering was considered. However, our assumption may be appropriate, since Wear [20] experimentally investigated the frequency dependence of ultrasonic backscatter from human trabecular bone and suggested that multiple scattering effects on the average were insignificant. Despite these limitations, the model successfully predicted the trends of ultrasonic attenuation in terms of frequency, trabecular size, and porosity, and provided a theoretical explanation for the quasi-linear assumption of attenuation coefficient for cancellous bone over a typical frequency range of 0.3–0.65 MHz. It is expected that this model may be a potentially valuable tool for understanding the interaction of ultrasound with bone.

## Conclusions

In conclusion, both absorption and scattering were considered for the prediction of ultrasonic attenuation of cancellous bone. The model successfully explained the quasi-linear behavior of the attenuation coefficient at low frequencies. The dependency of the ultrasonic attenuation on frequency, scatter size, and porosity was also predicted, and corresponding trends were consistent with previous reports. It is expected that our results may provide a valuable insight into ultrasonic propagation in human cancellous bone and aid in the bone assessment for the diagnosis of osteoporosis in clinics.

## Methods

A theoretical method was adopted from the literature [29] in which ultrasonic attenuation coefficient was assumed to be the summation of absorption and scattering from both the matrix and inclusions. This theory models the independent scattering and wave attenuation in two-phase composites and takes into consideration the attenuation property of the matrix and the scattering and absorption cross sections of a single inclusion. The exact formulation for the overall attenuation coefficient is given by

where \(\alpha\) is the total attenuation coefficient, \(\emptyset\) is the inclusion volume fraction, \(\alpha_{1}\) is the attenuation coefficient of matrix, \(n_{\text{s}}\) is the number density of inclusions, \(\gamma^{\text{sca}}\) is the scattering cross section of a single inclusion in the infinite matrix, and \(\gamma^{\text{abs}}\) is the absorption cross section of a single inclusion in the infinite matrix.

In this study, the cancellous bone was modeled as a two-phase composite material with inclusions of identical trabecular cylinders with radius \(a\) embedded in the homogeneous matrix of marrow. The typical trabecular thickness is approximately in the range of 80 to 440 µm [30,31,32,33]. Moreover, the trabecular cylinders were assumed to be long relative to the ultrasound beam cross section and perpendicular to the incident ultrasonic wave. The microstructure of cancellous bone in this model was consistent with one type of microstructures of cancellous bone [20, 34], and only rod shaped trabeculae were used. The absorption part from the trabecular cylinders was ignored due to small bone volume fraction and large difference of impedance in solid trabecular and fluid marrow. Under these assumptions, the total ultrasonic attenuation coefficient of cancellous bone can be written as

where \(a\) is the radius of trabecular cylinder and \(\emptyset\) is the bone volume fraction. The value of \(\left( {1 - \emptyset } \right)\) is defined as the porosity. The attenuation coefficient of marrow was determined using a power law equation as follows:

where \(f\) is the frequency in MHz, \(a_{0}\) and \(n\) were determined to be 3.15 and 1.15 by fitting Gammel’s experimental data [35], respectively.

The scattering cross section of an ultrasound wave scattered from a solid cylinder, \(\gamma^{\text{sca}}\), is given by [36]

where \(\varPi_{\text{s}}\) is the total power scattered by the cylinder per unit length, \(I\) is the incident intensity, \(k\) is the wavenumber defined by \(2\pi /\lambda\) (\(\lambda\) is the wavelength in the marrow surrounding the trabecular), \({\varepsilon_m}\) is the Newman constants (\({\varepsilon_0} = 1\) and \({\varepsilon_m} = 2\) for \(m > 0)\), and \(\eta_{m}\) is the phase-shift angle of the *n*th scattered wave. \(\eta_{m}\) can be calculated from [37]

where

and

where \(J_{m}\) and \(N_{m}\) correspond to Bessel and Newman functions, respectively; \(\rho\) and \(\rho_{1}\) are the density of marrow and trabecular, respectively, and \(\upsilon\) is the Poisson’s ratio of trabecular bone. The angle \(\xi_{m} \left( {x_{1} ,\nu } \right)\) is given by

where \(x_{1} = k_{1} a = \frac{\omega }{{c_{1} }}a\), \(x_{2} = k_{2} a = \frac{\omega }{{c_{2} }}a\), and \(c_{1}\) and \(c_{2}\) correspond to the longitudinal and shear wave velocity in the trabeculae, respectively.

Substitution of Eqs. (5)–(11) into Eq. (4) yields the scattering cross section of an ultrasound wave scattered from a solid cylinder, \(\gamma^{\text{sca}} .\) The overall ultrasonic attenuation coefficient of the cancellous bone, \(\alpha ,\) can be calculated by further substitution of Eqs. (3) and (4) into Eq. (2). It can be observed from the above formulations that the attenuation coefficients of cancellous bone depend on the incident frequency, the material properties of the matrix and inclusions, and volume fraction and geometry of the inclusions. For very long wavelength (i.e., value of \(ka\) goes to zero), the scattering cross section can be simply written as

In the following section, we will apply these formulations to predict the influence of frequency, scatter size, and bone volume fraction on the ultrasonic attenuation of cancellous bone. In the calculation, the values of density, longitudinal wave speed, and the Poison’s ratio of solid trabeculae were adopted from literature [38, 39] and set to 2000 kg/m^{3}, 3929 m/s, and 0.24, respectively. The density and wave speed in the marrow were assumed to be 900 kg/m^{3} and 1400 m/s, respectively [40].

## Availability of data and materials

Not applicable.

## References

Chan C, Mohamed N, Ima-Nirwana S, Chin K-Y. A review of knowledge, belief and practice regarding osteoporosis among adolescents and young adults. Int J Environ Res Public Health. 2018;15(8):1727.

Svedbom A, Ivergård M, Hernlund E, Rizzoli R, Kanis JA. Epidemiology and economic burden of osteoporosis in Switzerland. Arch Osteoporos. 2014;9(1):187.

Srichan W, Thasanasuwan W, Kijboonchoo K, Rojroongwasinkul N, Wimonpeerapattana W, Khouw I, Deurenberg P. Bone status measured by quantitative ultrasound: a comparison with DXA in Thai children. Eur J Clin Nutr. 2015;2008(9):1497–507.

Langton C, Palmer S, Porter R. The measurement of broadband ultrasonic attenuation in cancellous bone. Eng Med. 1984;13(2):89–91.

Truscott JG. Quantitative ultrasound—assessment of osteoporosis and bone status. In: Njeh CF, Hans D, Fuerst T, Glüer C-C, Genant HK, editors. Ultrasound in medicine and biology. London: Martin Dunitz; 1999. ISBN 1-85317-679-6, 420 pp. 2000. Vol. 26, No. 9. p. 1539.

Zaini WM, Md AA. Bone mineral density assessment in pre- and postmenopausal women: comparison between T-scores by heel QUS and DXA in HRPZII. Med J Malays. 2012;67(5):487–90.

Wear KA. Ultrasonic attenuation in human calcaneus from 0.2 to 1.7 MHz. IEEE Trans Ultrason Ferroelectr Freq Control. 2001;48(2):602–8.

Hodgskinson R, Njeh CF, Whitehead MA, Langton CM. The non-linear relationship between BUA and porosity in cancellous bone. Phys Med Biol. 1996;41(11):2411.

Chaffaı S, Peyrin F, Nuzzo S, Porcher R, Berger G, Laugier P. Ultrasonic characterization of human cancellous bone using transmission and backscatter measurements: relationships to density and microstructure. Bone. 2002;30(1):229–37.

Lee KI, Heui-Seol R, Suk Wang Y. Acoustic wave propagation in bovine cancellous bone: application of the Modified Biot–Attenborough model. J Acoust Soc Am. 2003;114(1):2284–93.

Currey JD. The mechanical adaptations of bones. Princeton: Princeton University Press; 2014.

Hutmacher DW, Schantz JT, Lam CXF, Tan KC, Lim TC. State of the art and future directions of scaffold-based bone engineering from a biomaterials perspective. J Tissue Eng Regen Med. 2007;1(4):245–60.

Nicholson P, Strelitzki R, Cleveland R, Bouxsein M. Scattering of ultrasound in cancellous bone: predictions from a theoretical model. J Biomech. 2000;33(4):503–6.

Biot MA. Theory of propagation of elastic waves in a fluid‐saturated porous solid. I. Low‐frequency range. J Acoust Soc Am. 1956;28(2):168–78.

Biot MA. Theory of propagation of elastic waves in a fluid‐saturated porous solid. II. Higher frequency range. J Acoust Soc Am. 1956;28(2):179–91.

Fellah ZEA, Chapelon JY, Berger S, Lauriks W, Depollier C. Ultrasonic wave propagation in human cancellous bone: application of Biot theory. J Acoust Soc Am. 2004;116(1):61–73.

Haire TJ, Langton CM. Biot theory: a review of its application to ultrasound propagation through cancellous bone. Bone. 1999;24(4):291–5.

Chaffaı S, Roberjot V, Peyrin F, Berger G, Laugier P. Frequency dependence of ultrasonic backscattering in cancellous bone: autocorrelation model and experimental results. J Acoust Soc Am. 2000;108(5):2403–11.

Strelitzki R, Nicholson P, Paech V. A model for ultrasonic scattering in cancellous bone based on velocity fluctuations in a binary mixture. Physiol Meas. 1998;19(2):189.

Wear KA. Frequency dependence of ultrasonic backscatter from human trabecular bone: theory and experiment. J Acoust Soc Am. 1999;106(6):3659–64.

Chaffaï S, Padilla F, Berger G, Laugier P. In vitro measurement of the frequency-dependent attenuation in cancellous bone between 0.2 and 2 MHz. J Acoust Soc Am. 2000;108(1):1281–9.

Pakula M, Padilla F, Laugier P. Influence of the filling fluid on frequency-dependent velocity and attenuation in cancellous bones between 0.35 and 2.5 MHz. J Acoust Soc Am. 2009;126(6):3301–10.

Padilla F, Laugier P. Prediction of ultrasound attenuation in cancellous bones using poroelasticity and scattering theories. In: 2001 IEEE ultrasonics symposium proceedings an international symposium (Cat No01CH37263); 7–10 Oct. 2001. vol. 1202. 2001. p. 1201–4.

Langton CM, Njeh CF, Hodgskinson R, Currey JD. Prediction of mechanical properties of the human calcaneus by broadband ultrasound attenuation. Bone. 1996;18(6):495–503.

Tavakoli MB, Evans JA. Dependence of the velocity and attenuation of ultrasound in bone on the mineral content. Phys Med Biol. 1991;36(11):1529.

Morse PM, Ingard KU. Theoretical acoustics. Princeton: Princeton University Press; 1968.

Mura T. Micromechanics of defects in solids. New York: Springer; 1987.

Feng C, Wang Y, Yang J. Effects of reorientation of graphene platelets (GPLs) on Young’s modulus of polymer composites under bi-axial stretching. Nanomaterials. 2018;8(1):532.

Biwa S. Independent scattering and wave attenuation in viscoelastic composites. Mech Mater. 2001;33(11):635–47.

Majumdar S, Genant H, Grampp S, Newitt D, Truong VH, Lin J, Mathur A. Correlation of trabecular bone structure with age, bone mineral density, and osteoporotic status: in vivo studies in the distal radius using high resolution magnetic resonance imaging. J Bone Miner Res. 1997;12(1):111–8.

Woo DG, Won YY, Han SK, Lim D. A biomechanical study of osteoporotic vertebral trabecular bone: the use of micro-CT and high-resolution finite element analysis. J Mech Sci Technol. 2007;21(4):593–601.

Nicholson PH, Müller R, Cheng XG, Rüegsegger P, Van DPG, Dequeker J, Boonen S. Quantitative ultrasound and trabecular architecture in the human calcaneus. J Bone Miner Res. 2010;16(10):1886–92.

Singh I. The architecture of cancellous bone. J Anat. 1978;127(Pt 2):305.

Singh I. The architecture of cancellous bone. J Anat. 1978;127(2):305–10.

Gammell PM, Croissette DH, Le Heyser RC. Temperature and frequency dependence of ultrasonic attenuation in selected tissues. Ultrasound Med Biol. 1979;5(3):269–77.

Morse PM, Ingard KU. Theoretical acoustics. Princeton: Princeton University Press; 1986.

Faran JJ Jr. Sound scattering by solid cylinders and spheres. J Acoust Soc Am. 1951;23(4):405–18.

Antich PP, Anderson JA, Ashman RB, Dowdey JE, Gonzales J, Murry RC, Zerwekh JE, Pak CY. Measurement of mechanical properties of bone material in vitro by ultrasound reflection: methodology and comparison with ultrasound transmission. J Bone Miner Res. 2010;6(4):417–26.

Moss WC, King MJ, Blackman EG. Skull flexure from blast waves: a mechanism for brain injury with implications for helmet design. Phys Rev Lett. 2009;103(10):108702.

Kubo T, Fujimori K, Cazier N, Saeki T, Matsukawa M. Properties of ultrasonic waves in bovine bone marrow. Ultrasound Med Biol. 2011;37(11):1923–9.

## Acknowledgements

Not applicable.

## Funding

This work was supported by Natural Science Foundation of Zhejiang Province (LQ15H180002), Scientific Research Staring Foundation for the Returned Overseas Chinese Scholars of Ministry of Education of China, Health Foundation for Creative Talents in Zhejiang Province, and Project Foundation for the College Young and Middle-aged Academic Leader of Zhejiang Province.

## Author information

### Authors and Affiliations

### Contributions

JL and YY designed the research. JL, LL, and JZ developed code. JL drafted the manuscript. JZ and LL collected data from literature. LL and JZ modified figures. JL, LL, JZ, and YY revised the article. All authors reviewed the final manuscript. All authors read and approved the final manuscript.

### Corresponding author

## Ethics declarations

### Ethics approval and consent to participate

Not applicable.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

## Additional information

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

**Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## About this article

### Cite this article

Liu, J., Lan, L., Zhou, J. *et al.* Influence of cancellous bone microstructure on ultrasonic attenuation: a theoretical prediction.
*BioMed Eng OnLine* **18**, 103 (2019). https://doi.org/10.1186/s12938-019-0724-4

Received:

Accepted:

Published:

DOI: https://doi.org/10.1186/s12938-019-0724-4