BioMed Research International
Volume 2018 (2018), Article ID 4786329, 7 pages
https://doi.org/10.1155/2018/4786329

Research Article

Variability in Ultrasound Backscatter Induced by Trabecular Microstructure Deterioration in Cancellous Bone

Xingxing Chou,¹ Feng Xu,¹ Ying Li,¹ Chengcheng Liu,² Dean Ta,1,3,4 and Lawrence H. Le⁵

¹Department of Electronic Engineering, Fudan University, Shanghai 200433, China
²Institute of Acoustics, Tongji University, No. 1239 Siping Road, Shanghai 200092, China
³State Key Laboratory of ASIC and System, Fudan University, Shanghai 200433, China
⁴Key Laboratory of Medical Imaging Computing and Computer Assisted Intervention (MICCAI) of Shanghai, Shanghai 200032, China
⁵Department of Radiology and Diagnostic Imaging, University of Alberta, Edmonton, AB, Canada

Correspondence should be addressed to Chengcheng Liu; chengchengliu@tongji.edu.cn and Dean Ta; tda@fudan.edu.cn

Received 26 October 2017; Accepted 1 January 2018; Published 29 January 2018

Academic Editor: Weibao Qiu

Copyright © 2018 Xingxing Chou et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

To determine the relationship between the ultrasonic backscatter parameters and trabecular microstructural variations in cancellous bone, three erosion procedures were performed to simulate various changes in the cancellous bone microstructure. The finite difference time domain (FDTD) method was used to simulate the backscatter signal in cancellous bone. Ultrasonic backscatter properties were derived as functions of the porosity when the ultrasound incident directions were perpendicular and parallel to the major trabeculae direction (MTD), respectively. The variability in the apparent backscatter coefficient (ABC) and apparent integrated backscatter (AIB) due to the trabecular microstructure was revealed. Significant negative correlations between the backscatter parameters (ABC and AIB) and the porosity of the cancellous bone were observed. The simulations showed that the ABC and AIB were influenced by the direction of the trabecular microstructural variations. The linear regressions between the ultrasonic backscatter parameters (ABC and AIB) and the porosity showed significantly different slopes for three erosion procedures when they are ultrasonically perpendicular (for ABC, −1.22 dB, −0.98 dB, and −0.46 dB; for AIB, −0.74 dB, −0.69 dB, and −0.25 dB) and parallel (for ABC, −1.87 dB, −0.69 dB, and −0.51 dB; for AIB, −0.9 dB, −0.5 dB, and −0.34 dB) to the MTD. This paper investigated the relationship between ultrasonic backscatter and cancellous bone microstructure deterioration and indicated that the ultrasonic backscatter could be affected by cancellous bone microstructure deterioration direction.

1. Introduction

Osteoporosis is a multifactorial skeletal disease characterized by decreased bone mass and deteriorated microarchitecture that leads to an increased risk of fracture [1]. Early detection and treatment of osteoporosis are essential for decreasing the risk of fracture. Lots of situations can lead to osteoporosis, such as age-related and microgravity-related situations [2–6]. Both aggravating trend of aging population and the aerospace development indicate the importance of early detection and treatment of osteoporosis.

Ultrasonic backscatter has shown great advantages and potential as a noninvasive tool for cancellous bone assessment [7–15]. Compared with ultrasonic through-transmission measurement, the backscatter measurement can be performed in pulse-echo mode with a single transducer and has easier access to skeletal sites such as the hip and spine. Hosokawa has studied the changes of the ultrasonic through-transmission signal [16]. In theory, the backscatter signal could provide more microstructural information; the backscatter signal is closely related to cancellous bone properties, including the bone mineral density (BMD), bone volume fraction (BV/TV), trabecular separation (Tb.Sp), ultimate strength, and Young's modulus [7–14, 17–25].

Bone is a tissue undergoing continuous construction and degradation; the location of cancellous bone in people’s body and different bone loss and growth processes determine the various cancellous bone microstructure. The trabecular orientation of bone tissue changes in response to mechanical stimuli; the process of bone loss destruction and reconstruction is anisotropic [26, 27]. In age-related osteoporosis, the trabecular elements perpendicular to the major trabecular direction (MTD) are more strongly lost than those parallel to the MTD. Because the weakly oriented trabecular elements to which large loads are not usually applied are the first to disappear, and the porosity of the bone increases (and the state of the osteoporosis progresses) [28, 29]. For spaceflight-induced bone loss, both the weak and strong oriented trabecular elements to which loads are not applied disappear rapidly [1, 4, 5, 28].

The ultrasonic backscattering and propagation are substantially affected by the cancellous bone microstructure [19]; thus the reliability needs to be further improved in ultrasonic backscatter apparatuses, especially assessment of the bone mass changes during bone loss and growth, because the relationship of the ultrasonic backscatter and cancellous bone microstructure parameters is not yet clearly understood. However, a detailed investigation on the relationship between the ultrasonic backscatter and cancellous bone microstructure is difficult because of the various cancellous bone microstructure in the bone loss process. Some image erosion algorithms have been used to simulate the degradation of cancellous bones [16, 29, 30]. Hosokawa realized various cancellous bone microstructure using image erosion methods [16, 29]. Three erosion procedures correspond to age-related, microgravity-related, and other reasons related to bone loss [16, 29]. This paper cited their algorithms to simulate the degradation of cancellous bone in normal bone loss, weightlessness, or microgravity environment bone loss. The ultrasonic backscatter parameters, such as apparent backscatter coefficient (ABC) and apparent integrated backscatter (AIB), are generally measured from fixed region of interest in the ultrasonic backscatter signal [21]. A detailed investigation on the relationship between the ultrasonic backscatter and cancellous bone microstructure is needed.

The objective of this study is to investigate the variability in ultrasonic backscatter induced by different deteriorations of trabecular microstructure. Image erosion methods were used to simulate the deteriorations of trabecular microstructure, and three erosion procedures were performed to realize deteriorations in the cancellous bone microstructure. The FDTD method was used to simulate the backscatter signal in cancellous bone. The variability in the ABC and AIB due to the deteriorations of trabecular microstructure was revealed.

2. Methods

The reconstruction of the 3D microcomputed tomographic (μ-CT) images is useful for the numerical analysis of cancellous bone [32, 33]. The finite difference time domain (FDTD) method is useful for simulating the ultrasound propagation in cancellous bone [29, 34]. The cancellous bone model for the FDTD simulation was realized by the reconstruction of the 3D microcomputed tomographic (-CT) images from the cancellous bone [35].

2.1. Cancellous Bone Erosion or Deteriorations

A cancellous bone specimen (approximately  mm) was sawed from a bovine distal tibia, and the trabecular image was provided by a μ-CT system (skyscan1076, Bruker micro-CT, Belgium) with a spatial resolution of 36.4 μ.

The binary image was obtained from the gray image by the automatic threshold function in the MATLAB to clearly distinguish between the trabeculae and bone marrow. The trabecular structure is with a MTD in most normal cancellous bone [36]. The 2D trabecular structures of the cancellous bone model are defined in - plane. As shown in Figure 1, a MTD along the -direction can be observed.

Figure 1: Deteriorations in trabecular microstructure simulated by three erosion procedures.

An image erosion technique was used to erode the edges of the trabeculae in the cancellous bone model and from which to simulate the bone loss process [30, 35]. The erosion procedure was to transform the solid bone into bone marrow [16]. The porosity increased by an increment approximately 2% at the same time. The trabeculae was eroded by three erosion procedures named A, B, and C; each of the erosion procedures was applied in different direction of the trabecular edges, from which to realize distinct changing processes of the trabecular microstructure in different bone loss processes. In erosion procedure A, the erosion was randomly distributed in every direction [16]. In the other two procedures B and C, the erosions were distributed in the - and -direction, respectively. Procedure C realized the age-related bone loss, procedure B realized the spaceflight-induced bone loss, and procedure A is for any other reason.

The direction of the erosion distributed was set in the erosion function. To the three erosion procedures, as an example, Figure 1 shows the different changes in trabecular structure induced by them. In Figure 1, the image in the left with a 61% porosity is original image before erosion, and the porosity is increasing with the erosion of original image. As an example, when the porosity is 68%, the images of three procedures are shown in the middle and their porosities are the same; compared with the original image, in the image of procedure A (top), the erosion is randomly distributed in every direction. In the other two procedures B (middle) and C (below), the erosions were distributed in the - and -direction, respectively. The solid bones (white) decrease in the corresponding direction. When the porosity is 81% (right), the difference of three erosion procedures is more obvious. It appears that the trabeculae changes of the three erosion procedures are different in directions. Thus, various trabecular microstructures of different causes could be realized by the three erosion procedures, and the trabecular orientation (or the pore orientation) in the -direction becomes stronger in the cancellous bone model. The porosities of all the cancellous bone models eroded in the three procedures are shown in Figure 2. Each erosion procedure is performed 16 times; porosity of cancellous bone before erosion is 61%, increased by an increment approximately 2% at the same time and increased from 60% to 90% based on the general range.

Figure 2: Porosity of cancellous bone with respect to erosion times.

2.2. Ultrasonic Backscatter Simulations

Figure 3 shows the FDTD simulation model, with a total region of  mm, for the ultrasonic backscatter measurement, and a cancellous bone model ( mm) was placed in the center [21]. The transmitting surfaces were with a diameter of 3.64 mm. The physical parameters of the simulation model are listed in Table 1.

Table 1: Physical parameter values of cancellous bone [31].

Trabeculae
Bone Marrow

First Lamé coefficient (GPa)
14.8
2.2

Second Lamé coefficient (GPa)
8.3
0

Density (kg/m³)
1960
1000

Normal resistance coefficient (s⁻¹)

75

Shear resistance coefficient (s⁻¹)

0

Figure 3: The geometry of the simulation model.

As shown in Figure 3, the ultrasonic propagation was along the -direction. With the cancellous bone model rotated by 90 degrees, the ultrasonic propagation is perpendicular to the MTD of the cancellous bone when the MTD is in -direction.

A Gaussian-modulated sinusoidal pulse was used as the ultrasound pressure source [21]:

where is the bandwidth and is central frequency. In the simulation, the parameters are defined as follows: MHz; the central frequency was set to 1 MHz; the space step was set to 36.4 μm, corresponding to the voxel size of the cancellous bone image; and the time step was 5  [37].

2.3. Backscatter Signal Analysis

Figure 4 shows a typical simulated backscatter signal at 1 . The backscattered signal of interest (SOI) was selected by a rectangular window of in length of the backscatter signal, where  μs [20].

Figure 4: Simulated backscatter signals and signal of interest selection.

ABC and AIB were defined as follows [7, 11, 38, 39]:

where is the amplitude spectrum of the backscatter SOI, is the reference spectrum of the backscatter signal reflected by a standard steel plate, and and correspond to the −6 dB effective frequency band. The central frequency is used in the calculation of the ABC.

3. Results

3.1. Backscatter Properties for Ultrasonic Propagation Perpendicular to MTD

Figure 5 shows the ABC versus the porosity induced by three different erosion procedures for ultrasonic propagation perpendicular to MTD. The ABC shows significant negative correlations with the porosity of the cancellous bone in all procedures (procedure A: ; procedure B: ; procedure C: ). In the three different erosion procedures, the cancellous bone microstructure undergoes different changes, and the ultrasonic backscatter signals are different. The ABC of procedure B are the smaller than those of the other two procedures, and in procedure B the trabecular microstructural variations are parallel to MTD. The linear fittings for ABC versus the porosity of the three erosion procedures are listed in Table 2. Significant differences are observed between the slopes for the three erosion procedures (procedure A: −1.22 dB; procedure B: −0.98 dB; procedure C: −0.46 dB). The absolute values of the slope for procedure C are the smallest.

Table 2: Linear fitting for the ultrasonic backscatter parameters (ABC and AIB) versus the porosity of the three erosion procedures perpendicular and parallel to .

Perpendicular to MTD
Parallel to MTD

Procedure A
Procedure B
Procedure C
Procedure A
Procedure B
Procedure C

ABC

    Intercept (dB)
51.59
31.16
0.02
96.27
13.10
285

    Slope (dB)
−1.22
−0.98
−0.46
−1.87
−0.69
−0.51

AIB

    Intercept (dB)
34.23
28.85
0.78
44.82
16.21
5.94

    Slope (dB)
−0.74
−0.69
−0.25
−0.90
−0.50
−0.34

of the values are below 0.01.

Figure 5: Relationships between ABC and cancellous bone porosity induced by different erosions for ultrasonic propagating perpendicular to the MTD.

The AIB is an important parameter of ultrasonic backscattering. The AIB versus the porosity of the cancellous bone induced by the different erosion procedures perpendicular to MTD is shown in Figure 6. Significant negative correlations with the porosity of the cancellous bone are observed in all procedures (procedure A: ; procedure B: ; procedure C: ). The liner regressions between the AIB and porosity showed significant differences in the slopes for the three erosion procedures (procedure A: −0.74 dB; procedure B: −0.69 dB; procedure C: −0.25 dB).

Figure 6: Relationships between AIB and cancellous bone porosity induced by different erosions for ultrasonic propagating perpendicular to the MTD.

3.2. Backscatter Properties for Ultrasonic Propagation along MTD

The ABC and AIB results versus the porosity of the cancellous bone induced by the different erosion procedures along the MTD are shown in Figures 7 and 8, respectively. Both the ABC and AIB show significant negative correlations with the porosity of the cancellous bone in all procedures (ABC: , −0.84, −0.71; AIB: , −0.90, −0.92). For the three different erosion procedures, with the increase in the cancellous bone porosity, the cancellous bone microstructure undergoes various changes. The changes in ABC and AIB are different for the same porosity. Based on the regression results listed in Table 2, the slopes for the three erosion procedures are different. Compared with the results perpendicular to the MTD, the values of ABC and AIB are smaller, and the corresponding slopes are different.

Figure 7: Relationships between ABC and cancellous bone porosity induced by different erosions for ultrasonic propagation along the MTD.

Figure 8: Relationships between AIB and cancellous bone porosity induced by different erosions for ultrasonic propagation along the MTD.

4. Discussion

ABC and AIB reflect the frequency-related intensity of the backscatter signal. The signal strength is mainly affected by the scattering cross section and attenuation. The reflected signal energy increases with the scattering cross section. When the ultrasonic incident direction is along the MTD, the scattering cross section is smaller than that for ultrasonic perpendicular to MTD and the values of ABC are smaller. Besides, the attenuation increases with trabecular bone length in the transmission direction. Thus, the ABC and AIB not only are influenced by bone mass or porosity but also will be influenced by trabecular microstructure of cancellous bone.

The simulations show that ABC and AIB vary differently for the three procedures with the increasing of porosity. With the notations , is used to represent ABC, AIB under procedure , the simulation results show that when the ultrasonic propagation is perpendicular to MTD, both ABC and AIB show significant negative correlations with the porosity of the cancellous bone in all procedures. The slopes of the liner fitting results for the three erosion procedures (ABC: −1.22 dB, −0.98 dB, and −0.46 dB; AIB: −0.74 dB, −0.69 dB, and −0.25 dB) are different. To the same porosity, the and are smaller than those of the other two procedures. In procedure B, the trabecular microstructure deterioration along the MTD, the length of the trabeculae in the propagation direction, is larger than those of the other procedures, which simulate the spaceflight-induced bone loss. The results indicate that the ABC and AIB may be smallest in the case of the trabecular microstructure deterioration along the MTD when propagation is perpendicular to MTD and with the increasing porosity. The slopes of procedure C are the smallest, which simulates the age-related bone loss. Besides, the difference values of the three procedures also increase for both ABC and AIB.

The and for the three erosion procedures are also different when the propagation is along MTD, but compared to the results of propagation perpendicular to MTD, the values of ABC and AIB are smaller, the corresponding slopes of procedures A, B, and C (perpendicular MTD: ABC: −1.22 dB, −0.98 dB, −0.46 dB; AIB: −0.74 dB, −0.69 dB, −0.25 dB; along MTD: ABC: −1.87 dB, −0.69 dB, −0.51 dB; AIB: −0.9 dB, −0.5 dB, −0.34 dB) are also different, because there is a major trabecular direction in cancellous bone. The results along MTD further indicate the variability in ABC and AIB induced by the trabecular microstructure in cancellous bone and illustrate that ABC and AIB are sensitive to the trabecular microstructure.

It is the same as the experimental results in the previous study [40]; ABC and AIB both had negative correlations with the porosity. The ABC and AIB are influenced by the trabecular microstructure deterioration direction; the slops of linear fitting between them (ABC and AIB) and porosity indicate that, for the linear evaluation of BV/TV or porosity using ABC or AIB to the age-related and spaceflight-induced bone loss, the deterioration direction maybe should be considered.

In a previous study it was suggested that the ultrasound backscatter was affected by the anisotropic microstructure [19]. And to the same reason induced bone loss, to the trabecular microstructures are various, but the deterioration direction is the same. In the present study, the effect of deterioration direction is investigated. And it is significant that when the deterioration direction is considered, the liner assessment of BV/TV by the ABC or AIB is more accurate.

This study investigated the relationship between ultrasonic backscatter and cancellous bone microstructure deterioration and indicated that the ultrasonic backscatter was affected by cancellous bone microstructure deterioration direction, and we just discussed the parameters of ABC and AIB. Therefore, the study of the trabecular microstructure effect on the different ultrasonic backscatter parameters without the porosity should be elaborated upon in the future.

5. Conclusion

The variabilities of ABC and AIB induced by different direction deteriorations of trabecular microstructure were investigated. ABC and AIB showed significantly negative correlations with the porosity of the cancellous bone. ABC and AIB were sensitive to the trabecular microstructure; they were confined to erosion procedures from which three different direction cancellous bone microstructure changes were revealed. The ABC and AIB are affected by the trabecular microstructure deterioration direction. When using ABC and AIB accurately in the evaluation of cancellous bone mass on different reason related bone loss, the effect of the trabecular microstructure deterioration direction maybe should be considered.

Conflicts of Interest

The authors state that they have no conflicts of interest.

Acknowledgments

This work was supported by the NSFC (11525416, 11504057, 11604054, and 11327405) and the International Scientific and Technological Cooperation Project of Shanghai (17510710700).

References

L. Vico, V. E. Novikov, J. M. Very, and C. Alexandre, “Bone histomorphometric comparison of rat tibial metaphysis after 7-day tail suspension vs. 7-day spaceflight,” Aviation, Space, and Environmental Medicine, vol. 62, no. 1, pp. 26–31, 1991.

D. Grimm, J. Grosse, M. Wehland et al., “The impact of microgravity on bone in humans,” Bone, vol. 87, pp. 44–56, 2016.

S. Khosla and B. L. Riggs, “Pathophysiology of age-related bone loss and osteoporosis,” Endocrinology and Metabolism Clinics of North America, vol. 34, no. 4, pp. 1015–1030, 2005.

R. Y. Lau and X. Guo, “A Review on Current Osteoporosis Research: With Special Focus on Disuse Bone Loss,” Journal of Osteoporosis, vol. 2011, Article ID 293808, 6 pages, 2011.

M. P. Nagaraja and D. Risin, “The current state of bone loss research: data from spaceflight and microgravity simulators,” Journal of Cellular Biochemistry, vol. 114, no. 5, pp. 1001–1008, 2013.

J. D. Sibonga, “Spaceflight-induced bone loss: Is there an Osteoporosis Risk?” Current Osteoporosis Reports, vol. 11, no. 2, pp. 92–98, 2013.

B. K. Hoffmeister, A. P. Holt, and S. C. Kaste, “Effect of the cortex on ultrasonic backscatter measurements of cancellous bone,” Physics in Medicine and Biology, vol. 56, no. 19, article no. 006, pp. 6243–6255, 2011.

B. K. Hoffmeister, D. P. Johnson, J. A. Janeski et al., “Ultrasonic characterization of human cancellous bone in vitro using three different apparent backscatter parameters in the frequency range 0.6-15.0 MHz,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 55, no. 7, pp. 1442–1452, 2008.

K. H. Chong, B. K. Poh, N. A. Jamil, N. A. Kamaruddin, and P. Deurenberg, “Radial quantitative ultrasound and dual energy X-ray absorptiometry: Intermethod agreement for bone status assessment in children,” BioMed Research International, vol. 2015, Article ID 232876, 7 pages, 2015.

D.-A. Ta, G.-H. Zhou, W.-Q. Wang, and J.-G. Yu, “Measurement of Spectral Maximum Shift of Ultrasonic Backscatter Signals in Cancellous Bone,” in Proceedings of the IEEE Engineering in Medicine and Biology 27th Annual Conference, pp. 2703–2706, Shanghai, China, January 2006.

D. Ta, W. Wang, K. Huang, Y. Wang, and L. H. Le, “Analysis of frequency dependence of ultrasonic backscatter coefficient in cancellous bone,” The Journal of the Acoustical Society of America, vol. 124, no. 6, pp. 4083–4090, 2009.

Y. Li, D. Liu, K. Xu, D. Ta, L. H. Le, and W. Wang, “Transverse and Oblique Long Bone Fracture Evaluation by Low Order Ultrasonic Guided Waves: A Simulation Study,” BioMed Research International, vol. 2017, Article ID 3083141, 10 pages, 2017.

B. K. Hoffmeister, A. R. Wilson, M. J. Gilbert, and M. E. Sellers, “A backscatter difference technique for ultrasonic bone assessment,” The Journal of the Acoustical Society of America, vol. 132, no. 6, pp. 4069–4076, 2012.

C. C. Liu, D. A. Ta, and W. Q. Wang, “Effect of trabecular bone material properties on ultrasonic backscattering signals,” Chinese Journal of Acoustics, vol. 38, no. 1, pp. 73–82, 2014.

K. A. Wear, “Ultrasonic scattering from cancellous bone: A review,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 55, no. 7, pp. 1432–1441, 2008.

A. Hosokawa, “Numerical analysis of variability in ultrasound propagation properties induced by trabecular microstructure in cancellous bone,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 56, no. 4, pp. 738–747, 2009.

C. Apostolopoulos and D. Deligianni, “Influence of microarchitecture on cancellous bone ultrasonic backscatter,” in Proceedings of the 8th International Workshop on Mathematical Methods in Scattering Theory and Biomedical Engineering, pp. 211–217, September 2007.

B. K. Hoffmeister, C. I. Jones III, G. J. Caldwell, and S. C. Kaste, “Ultrasonic characterization of cancellous bone using apparent integrated backscatter,” Physics in Medicine and Biology, vol. 51, no. 11, article no. 002, pp. 2715–2727, 2006.

C. Liu, D. Ta, F. Fujita et al., “The relationship between ultrasonic backscatter and trabecular anisotropic microstructure in cancellous bone,” Journal of Applied Physics, vol. 115, no. 6, Article ID 064906, 2014.

C. Liu, T. Tang, F. Xu et al., “Signal of Interest Selection Standard for Ultrasonic Backscatter in Cancellous Bone Evaluation,” Ultrasound in Medicine & Biology, vol. 41, no. 10, pp. 2714–2721, 2015.

C. Liu, H. Han, D. Ta, and W. Wang, “Effect of selected signals of interest on ultrasonic backscattering measurement in cancellous bones,” Science China Physics, Mechanics & Astronomy, vol. 56, no. 7, pp. 1310–1316, 2013.

B. K. Hoffmeister, J. A. McPherson, M. R. Smathers, P. L. Spinolo, and M. E. Sellers, “Ultrasonic backscatter from cancellous bone: The apparent backscatter transfer function,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 62, no. 12, pp. 2115–2125, 2015.

K. A. Wear, “Frequency dependence of ultrasonic backscatter from human trabecular bone: Theory and experiment,” The Journal of the Acoustical Society of America, vol. 106, no. 6, pp. 3659–2664, 1999.

K. A. Wear, A. P. Stuber, and J. C. Reynolds, “Relationships of ultrasonic backscatter with ultrasonic attenuation, sound speed and bone mineral density in human calcaneus,” Ultrasound in Medicine & Biology, vol. 26, no. 8, pp. 1311–1316, 2000.

R. Smolíková, M. P. M. P. Wachowiak, and J. M. J. M. Zurada, “An information-theoretic approach to estimating ultrasound backscatter characteristics.,” Computers in Biology and Medicine, vol. 34, no. 4, pp. 355–370, 2004.

J. Wolff, The Law of Bone Remodelling, Springer, Berlin, Germany, 1986.

S. Ramtani, J. M. Garcia, and M. Doblare, “Computer simulation of an adaptive damage-bone remodeling law applied to three unit-bone bars structure,” Computers in Biology and Medicine, vol. 34, no. 4, pp. 259–273, 2004.

H. M. Frost, “Wolff's Law and bone's structural adaptations to mechanical usage: an overview for clinicians,” The Angle Orthodontist, vol. 64, no. 3, pp. 175–188, 1994.

A. Hosokawa, “Effect of porosity distribution in the propagation direction on ultrasound waves through cancellous bone,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 57, no. 6, pp. 1320–1328, 2010.

G. Luo, J. J. Kaufman, A. Chiabrera et al., “Computational methods for ultrasonic bone assessment,” Ultrasound in Medicine & Biology, vol. 25, no. 5, pp. 823–830, 1999.

A. Hosokawa, “Development of a numerical cancellous bone model for finite-difference time-domain simulations of ultrasound propagation,” IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, vol. 55, no. 6, pp. 1219–1233, 2008.

R. P. Gilbert, P. Guyenne, and J. Li, “Numerical investigation of ultrasonic attenuation through 2D trabecular bone structures reconstructed from CT scans and random realizations,” Computers in Biology and Medicine, vol. 45, no. 1, pp. 143–156, 2014.

D. Sreenivasan, P. T. Tu, M. Dickinson et al., “Computer modelling integrated with micro-CT and material testing provides additional insight to evaluate bone treatments: Application to a beta-glycan derived whey protein mice model,” Computers in Biology and Medicine, vol. 68, pp. 9–20, 2016.

J. Litniewski, “Statistics of envelope of high-frequency ultrasonic backscatter from trabecular bone: Simulation study,” Archives of Acoustics, vol. 35, no. 3, pp. 349–360, 2010.

L. J. Gibson, “The mechanical behaviour of cancellous bone,” Journal of Biomechanics, vol. 18, no. 5, pp. 317–328, 1985.

G. Haïat, F. Padilla, F. Peyrin, and P. Laugier, “Variation of ultrasonic parameters with microstructure and material properties of trabecular bone: A 3D model simulation,” Journal of Bone and Mineral Research, vol. 22, no. 5, pp. 665–674, 2007.

J. Virieux, “SH- wave propagation in heterogeneous media: velocity- stress finite-difference method.,” Geophysics, vol. 49, no. 11, pp. 1933–1942, 1984.

S. Chaffaí, F. Peyrin, S. Nuzzo, R. Porcher, G. Berger, and P. Laugier, “Ultrasonic characterization of human cancellous bone using transmission and backscatter measurements: Relationships to density and microstructure,” Bone, vol. 30, no. 1, pp. 229–237, 2002.

O. Riekkinen, M. A. Hakulinen, J. Töyräs, and J. S. Jurvelin, “Spatial variation of acoustic properties is related with mechanical properties of trabecular bone,” Physics in Medicine and Biology, vol. 52, no. 23, pp. 6961–6968, 2007.

Y.-Q. Jiang, C.-C. Liu, R.-Y. Li et al., “Analysis of apparent integrated backscatter coefficient and backscattered spectral centroid shift in calcaneus invivo for the ultrasonic evaluation of osteoporosis,” Ultrasound in Medicine & Biology, vol. 40, no. 6, pp. 1307–1317, 2014.