Measurement of dependence of backscatter coefficient from cylinders on frequency and diameter using focused transducers--with applications in trabecular bone

Mechanisms of Interaction of Ultrasound With Cancellous Bone: A Review

Ultrasound is now a clinically accepted modality in the management of osteoporosis. The most common commercial clinical devices assess fracture risk from measurements of attenuation and sound speed in cancellous bone. This review discusses fundamental mechanisms underlying the interaction between ultrasound and cancellous bone. Because of its two-phase structure (mineralized trabecular network embedded in soft tissue-marrow), its anisotropy, and its inhomogeneity, cancellous bone is more difficult to characterize than most soft tissues. Experimental data for the dependencies of attenuation, sound speed, dispersion, and scattering on ultrasound frequency, bone mineral density, composition, microstructure, and mechanical properties are presented. The relative roles of absorption, scattering, and phase cancellation in determining attenuation measurements in vitro and in vivo are delineated. Common speed of sound metrics, which entail measurements of transit times of pulse leading edges (to avoid multipath interference), are greatly influenced by attenuation, dispersion, and system properties, including center frequency and bandwidth. However, a theoretical model has been shown to be effective for correction for these confounding factors in vitro and in vivo. Theoretical and phantom models are presented to elucidate why cancellous bone exhibits negative dispersion, unlike soft tissue, which exhibits positive dispersion. Signal processing methods are presented for separating "fast" and "slow" waves (predicted by poroelasticity theory and supported in cancellous bone) even when the two waves overlap in time and frequency domains. Models to explain dependencies of scattering on frequency and mean trabecular thickness are presented and compared with measurements. Anisotropy, the effect of the fluid filler medium (marrow in vivo or water in vitro), phantoms, computational modeling of ultrasound propagation, acoustic microscopy, and nonlinear properties in cancellous bone are also discussed.

Characterization of a polymer, open-cell rigid foam that simulates the ultrasonic properties of cancellous bone

Materials that simulate the ultrasonic properties of tissues are used widely for clinical and research purposes. However, relatively few materials are known to simulate the ultrasonic properties of cancellous bone. The goal of the present study was to investigate the suitability of using a polymer, open-cell rigid foam (OCRF) produced by Sawbones^®. Measurements were performed on OCRF specimens with four different densities. Ultrasonic speed of sound and normalized broadband ultrasonic attenuation were measured with a 0.5 MHz transducer. Three backscatter parameters were measured with a 5 MHz transducer: apparent integrated backscatter, frequency slope of apparent backscatter, and normalized mean of the backscatter difference. X-ray micro-computed tomography was used to measure the microstructural characteristics of the OCRF specimens. The trabecular thickness and relative bone volume of the OCRF specimens were similar to those of human cancellous bone, but the trabecular separation was greater. In most cases, the ultrasonic properties of the OCRF specimens were similar to values reported in the literature for cancellous bone, including dependence on density. In addition, the OCRF specimens exhibited an ultrasonic anisotropy similar to that reported for cancellous bone.

A contactless ultrasonic surface wave approach to characterize distributed cracking damage in concrete

We describe an approach that utilizes ultrasonic surface wave backscatter measurements to characterize the volume content of relatively small distributed defects (microcrack networks) in concrete. A simplified weak scattering model is used to demonstrate that the scattered wave field projected in the direction of the surface wave propagation is relatively insensitive to scatterers that are smaller than the propagating wavelength, while the scattered field projected in the opposite direction is more sensitive to sub-wavelength scatterers. Distributed microcracks in the concrete serve as the small scatterers that interact with a propagating surface wave. Data from a finite element simulation were used to demonstrate the viability of the proposed approach, and also to optimize a testing configuration to collect data. Simulations were validated through experimental measurements of ultrasonic backscattered surface waves from test samples of concrete constructed with different concentrations of fiber filler (0.0, 0.3 and 0.6%) to mimic increasing microcrack volume density and then samples with actual cracking induced by controlled thermal cycles. A surface wave was induced in the concrete samples by a 50kHz ultrasonic source operating 10mm above the surface at an angle of incidence of 9°. Silicon-based miniature MEMS acoustic sensors located a few millimeters above the concrete surface both behind and in front of the sender were used to detect leaky ultrasonic surface waves emanating from concrete. A normalized backscattered energy parameter was calculated from the signals. Statistically significant differences in the normalized backscattered energy were observed between concrete samples with varying levels of simulated and actual cracking damage volume.

Comparison of Two Inexpensive Rapid Prototyping Methods for Manufacturing Filament Target Ultrasound Phantoms

Current use of 3-D printers to manufacture ultrasound phantoms is limited to relatively expensive photopolymer jetting printers. The present work investigates the feasibility of using two common and inexpensive 3-D printer technologies, fused deposition modeling (FDM) and digital light processing (DLP), to print custom filament target phantoms. Acoustic characteristics obtained from printed solid blocks indicated that the printing materials-acrylonitrile butadiene styrene and polylactic acid for FDM and a photopolymer for DLP printing-were appropriate for use as scatterers. A regular grid of filaments was printed to study printing accuracy. As a proof of concept of the phantom manufacturing process, a complex pattern of filament targets was placed in de-ionized water to create a phantom, which was then imaged using an ultrasound imager. The pattern was clearly identifiable, although multiple reflections were observed, which underscores the importance of future work to enhance printing resolution. This goal is deemed possible using improvement of the DLP printing setup.

Nonlinear attenuation and dispersion in human calcaneus in vitro: statistical validation and relationships to microarchitecture

Through-transmission measurements were performed on 30 human calcaneus samples in vitro. Nonlinear attenuation and dispersion measurements were investigated by estimating 95% confidence intervals of coefficients of polynomial expansions of log magnitude and phase of transmission coefficients. Bone mineral density (BMD) was measured with dual x-ray absorptiometry. Microarchitecture was measured with microcomputed tomography. Statistically significant nonlinear attenuation and nonzero dispersion were confirmed for a clinical bandwidth of 300-750 kHz in 40%-43% of bone samples. The mean linear coefficient for attenuation was 10.3 dB/cm MHz [95% confidence interval (CI): 9.0-11.6 dB/cm MHz]. The mean quadratic coefficient for attenuation was 1.6 dB/cm MHz(2) (95% CI: 0.4-2.8 dB/cm MHz(2)). Nonlinear attenuation provided little information regarding BMD or microarchitecture. The quadratic coefficient for phase (which is related to dispersion) showed moderate correlations with BMD (r = -0.65; 95% CI: -0.82 to -0.36), bone surface-to-volume ratio (r = 0.47; 95% CI: 0.12-0.72) and trabecular thickness (r = -0.40; 95% CI: -0.67 to -0.03). Dispersion was proportional to bone volume fraction raised to an exponent of 2.1 ± 0.2, which is similar to the value for parallel nylon-wire phantoms (2.4 ± 0.2) and supports a multiple-scattering model for dispersion.

Scattering by single physically large and weak scatterers in the beam of a single-element transducer

Quantitative ultrasound techniques are generally applied to characterize media whose scattering sites are considered to be small compared to a wavelength. In this study, the backscattered response of single weakly scattering spheres and cylinders with diameters comparable to the beam width of a 2.25 MHz single-element transducer were simulated and measured in the transducer focal plane to investigate the impact of physically large scatterers. The responses from large single spherical scatterers at the focus were found to closely match the plane-wave response. The responses from large cylindrical scatterers at the focus were found to differ from the plane-wave response by a factor of f(-1). Normalized spectra from simulations and measurements were in close agreement: the fall-off of the responses as a function of lateral position agreed to within 2 dB for spherical scatterers and to within 3.5 dB for cylindrical scatterers. In both measurement and simulation, single scatterer diameter estimates were biased by less than 3% for a more highly focused transducer compared to estimates for a more weakly focused transducer. The results suggest that quantitative ultrasound techniques may produce physically meaningful size estimates for media whose response is dominated by scatterers comparable in size to the transducer beam.

Comparison of phase velocity in trabecular bone mimicking-phantoms by time domain numerical (EFIT) and analytical multiple scattering approaches

The corrected Waterman-Truell model and the Elastodynamic Finite Integration Technique were used to analyze the ultrasonic wave dispersion in trabecular bones mimicking phantoms. A simple two-phase model of the trabecular bone is assumed; the trabeculae structure and the bone marrow. The phase velocity for frequencies within the range from 400kHz to 800kHz were computed for different scatterer arrays varying their dimensions and number. The theoretical and numerical results were compared to experimental published data, obtained from a mimicking phantom composed by a periodic array of nylon shreds (trabeculae array) immersed in a water tank. Our results showed an excellent consistency when compared to experimental data. The negative dispersions of -8.48m/s/MHz and -9.16m/s/MHz were computed by the multiple scattering method and the numerical approach, respectively, where the latter is closer to the experimental dispersion of -12.09m/s/MHz. Similar result has been reported in the literature, where the dispersion predicted by the Generalized Self-Consistent Method [J. Acoust. Soc. Am. 124 (2008) 4047] is -9.96m/s/MHz.

A review of tissue substitutes for ultrasound imaging

The characterization and calibration of ultrasound imaging systems requires tissue-mimicking phantoms with known acoustic properties, dimensions and internal features. Tissue phantoms are available commercially for a range of medical applications. However, commercial phantoms may not be suitable in ultrasound system design or for evaluation of novel imaging techniques. It is often desirable to have the ability to tailor acoustic properties and phantom configurations for specific applications. A multitude of tissue-mimicking materials and phantoms are described in the literature that have been created using a variety of materials and preparation techniques and that have modeled a range of biological systems. This paper reviews ultrasound tissue-mimicking materials and phantom fabrication techniques that have been developed over the past four decades, and describes the benefits and disadvantages of the processes. Both soft tissue and hard tissue substitutes are explored.

Frequency dependence of average phase shift from human calcaneus in vitro

If dispersion in a medium is weak and approximately linear with frequency (over the experimental band of frequencies), then it can be shown that the constant term in a polynomial representation of phase shift as a function of frequency can produce errors in measurements of phase-velocity differences in through-transmission, substitution experiments. A method for suppressing the effects of the constant phase shift in the context of the single-wave-model was tested on measurements from 30 cancellous human calcaneus samples in vitro. Without adjustment for constant phase shifts, the estimated phase velocity at 500 kHz was 1516+/-6 m/s (mean+/-standard error), and the estimated dispersion was -24+/-4 m/s MHz (mean+/-standard error). With adjustment for constant phase shifts, the estimated mean velocity decreased by 4-9 m/s, and the estimated magnitude of mean dispersion decreased by 50%-100%. The average correlation coefficient between the measured attenuation coefficient and frequency was 0.997+/-0.0026 (mean+/-standard deviation), suggesting that the signal for each sample was dominated by one wave. A single-wave, linearly dispersive model conformed to measured complex transfer functions from the 30 cancellous-bone samples with an average root-mean-square error of 1.9%+/-1.0%.

Influence of the filling fluid on frequency-dependent velocity and attenuation in cancellous bones between 0.35 and 2.5 MHz

The paper is focused on experiments on human cancellous bones filled with different fluids with the goal of evaluating their contribution to velocity dispersion, absorption, and scattering mechanisms. The specimens were measured first filled with marrow and subsequently, after marrow removal, with water and alcohol. No significant influence of the fluids was evidenced on the attenuation coefficient. Given the absence of impact of viscosity of the saturating fluid, the authors hypothesized that the source of attenuation is associated with viscoelastic absorption in the solid trabeculae and with scattering. Alteration of scattering obtained by changing the acoustic impedance mismatch between the fluid (alcohol vs water) and the trabeculae was reflected neither in the attenuation nor in its slope. This led the authors to suggest that longitudinal-to-shear scattering together with absorption in the solid phase are candidates as main sources for the attenuation. The differences in velocity values indicate that the elastic properties of the fluid are main determinants of the phase velocity. This finding is particularly significant in the context of /in vivo/ measurements, because it demonstrates that the subject-dependent properties of marrow may partly explain the inter-subject variability of speed of sound values.

Simulation of elastic wave scattering in cells and tissues at the microscopic level

The scattering of longitudinal and shear waves from spherical, nucleated cells and three-dimensional tissues with simple and hierarchical microstructures was numerically modeled at the microscopic level using an iterative multipole approach. The cells were modeled with a concentric core-shell (nucleus-cytoplasm) structure embedded in an extracellular matrix. Using vector multipole expansions and boundary conditions, scattering solutions were derived for single cells with either solid or fluid properties for each of the cell components. Tissues were modeled as structured packings of cells. Multiple scattering between cells was simulated using addition theorems to translate the multipole fields from cell to cell in an iterative process. Backscattering simulations of single cells indicated that changes in the shear properties and nuclear diameter had the greatest effect on the frequency spectra. Simulated wave field images and high-frequency spectra (15-75 MHz) from tissues containing 1211-2137 cells exhibited up to 20% enhancement of the field amplitudes at the plasma membrane, significant changes in spectral features due to neoplastic and other microstructural alterations, and a detection threshold of approximately 8.5% infiltration of tumor cells into normal tissue. These findings suggest that histology-based simulations may provide insight into fundamental ultrasound-tissue interactions and help in the development of new medical technologies.

The dependencies of phase velocity and dispersion on volume fraction in cancellous-bone-mimicking phantoms

Frequency-dependent phase velocity was measured in eight cancellous-bone-mimicking phantoms consisting of suspensions of randomly oriented nylon filaments (simulating trabeculae) in a soft-tissue-mimicking medium (simulating marrow). Trabecular thicknesses ranged from 152 to 356 mum. Volume fractions of nylon filament material ranged from 0% to 10%. Phase velocity varied approximately linearly with frequency over the range from 300 to 700 kHz. The increase in phase velocity (compared with phase velocity in a phantom containing no filaments) at 500 kHz was approximately proportional to volume fraction occupied by nylon filaments. The derivative of phase velocity with respect to frequency was negative and exhibited nonlinear, monotonically decreasing dependence on volume fraction. The dependencies of phase velocity and its derivative on volume fraction in these phantoms were similar to those reported in previous studies on (1) human cancellous bone and (2) phantoms consisting of parallel nylon wires immersed in water.

Analysis of frequency dependence of ultrasonic backscatter coefficient in cancellous bone

The ultrasonic scattering mechanism in cancellous bone is investigated theoretically and a model describing the frequency dependence of ultrasonic scattering from cancellous bone is presented. The ultrasonic backscatter coefficient (BSC) of bovine tibiae, human calcanei in vitro and in vivo, were measured and discussed. The data of BSC were also fitted by polynomial. The results demonstrate that BSC is a nonlinear function of frequency and increases with frequency. A good agreement was obtained between BSC values from theory and experiment. Also, the high correlation coefficient between BSC and bone mineral density was obtained, r=0.85+/-0.07 (mean+/-SD) (n=15, p<0.001). Based on the values of BSC, the status of cancellous bone and the degree of osteoporotic fracture risk may be assessed.

Velocity dispersion in trabecular bone: influence of multiple scattering and of absorption

Speed of sound measurements are widely used clinically to assess bone strength. Trabecular bone is an attenuating composite material in which negative values of velocity dispersion have been measured, this behavior remaining poorly explained physically. The aim of this work is to describe the ultrasonic propagation in trabecular bone modeled by infinite cylinders immersed in a saturating matrix, and to derive the physical determinants of velocity dispersion. A homogenization model accounting for the coupling of multiple scattering and absorption phenomena allows the computation of phase velocity and of dispersion while varying bone properties. The present model is adapted from the generalized self-consistent method developed in the work of Yang and Mal [(1994). "Multiple-scattering of elastic waves in a fiber-reinforced composite," J. Mech. Phys. Solids 42, 1945-1968]. It predicts negative values of velocity dispersion, in agreement with experimental results obtained in phantoms mimicking trabecular bone. In trabecular bone, mostly negative and also positive values of velocity dispersion are predicted by the model, which span within the range of values measured experimentally. Scattering effects are responsible for the negative values of dispersion, whereas the frequency dependence of the absorption coefficient in bone marrow and/or in the trabeculae results in an increase in dispersion, which may then become positive.

Ultrasonic attenuation in parallel-nylon-wire cancellous-bone-mimicking phantoms

Attenuation coefficients between 1.5 and 3.5 MHz were measured on four parallel-nylon-wire arrays (simulating cancellous bone) with four different wire diameters (150, 200, 250, and 300 microm). Interwire spacing was 800 microm for all four parallel-nylon-wire arrays. The measured frequency dependencies of attenuation were consistent with theoretical predications based on Faran's theory, which considers the component of attenuation due to scattering of longitudinal waves.

Mechanisms for attenuation in cancellous-bone-mimicking phantoms

Broadband ultrasound attenuation (BUA) in cancellous bone is useful for prediction of osteoporotic fracture risk, but its causes are not well understood. To investigate attenuation mechanisms, 9 cancellous-bone-mimicking phantoms containing nylon filaments (simulating bone trabeculae) embedded within soft-tissue-mimicking fluid (simulating marrow) were interrogated. The measurements of frequency-dependent attenuation coefficient had 3 separable components: 1) a linear (with frequency) component attributable to absorption in the soft-tissue-mimicking fluid, 2) a quasilinear (with frequency) component, which may include absorption in and longitudinal-shear mode conversion by the nylon filaments, and 3) a nonlinear (with frequency) component, which may be attributable to longitudinal-longitudinal scattering by the nylon filaments. The slope of total linear (with frequency) attenuation coefficient (sum of components #1 and #2) versus frequency was found to increase linearly with volume fraction, consistent with reported measurements on cancellous bone. Backscatter coefficient measurements in the 9 phantoms supported the claim that the nonlinear (with frequency) component of attenuation coefficient (component #3) was closely associated with longitudinal-longitudinal scattering. This work represents the first experimental separation of these 3 components of attenuation in cancellous bone-mimicking phantoms.

Frequency dependence of backscatter from thin, oblique, finite-length cylinders measured with a focused transducer-with applications in cancellous bone

A model is presented for the echo from a thin, oblique, finite-length cylinder. The echo is calculated from the line integral of the transducer directivity pattern along the cylinder axis. The model was validated with broadband pulse-echo measurements from (1) a perpendicular (to the ultrasound beam) nylon wire as a function of lateral displacement from the beam center, (2) a tilted nylon wire as a function of the angle of inclination relative to the ultrasound beam, and (3) a quasi-parallel-nylon-wire phantom, which mimicked the scattering properties of cancellous bone. The transducer directivity pattern (as a function of position and frequency) was measured with a membrane hydrophone. The model predicts an approximately cubic frequency dependence of backscatter coefficient from the phantom, as has been observed experimentally in cancellous bone. The model also predicts the relationship between the cylinder length and the exponent of a power law fit to backscatter coefficient versus frequency, which is 4 for very short (compared to a wavelength) cylinders and asymptotically approaches 3 for very long cylinders.

Relationship between ultrasonic attenuation, size and axial strain parameters for ex vivo atherosclerotic carotid plaque

Many ultrasonic parameters, primarily related to attenuation and scatterer size, have been used to characterize the composition of atherosclerotic plaque tissue. In this study, we combine elastographic (axial strain ratio) and ultrasonic tissue characterization parameters, namely the attenuation coefficient and a scattering parameter associated with an "equivalent" scatterer size to delineate between fibrous, calcified, and lipidic plaque tissue. We present results obtained from 44 ex vivo atherosclerotic plaque specimens obtained after carotid endarterectomy on human patients. Our results in the frequency range 2.5 - 7.5 MHz indicate that softer plaques (with higher values of the strain ratio) are usually associated with larger equivalent scatterer size estimates (200 - 500 microm) and lower values of the attenuation coefficient slope (<1 dB/cm/MHz). On the other hand, stiffer plaques (with lower strain ratio values) are associated with smaller equivalent scatterer size estimates (100 - 200 microm) and higher values of the attenuation coefficient slope (1 - 3 dB/cm/MHz). These results indicate that ultrasonic tissue characterization and strain parameters have the potential to differentiate between different plaque types. These parameters can be estimated from radio-frequency data acquired under in vivo conditions and may help the clinician decide on appropriate interventional techniques.

Ultrasonic scattering from cancellous bone: a review

This paper reviews theory, measurements, and computer simulations of scattering from cancellous bone reported by many laboratories. Three theoretical models (binary mixture, Faran cylinder, and weak scattering) for scattering from cancellous bone have demonstrated some consistency with measurements of backscatter. Backscatter is moderately correlated with bone mineral density in human calcaneus in vitro (r(2) = 0.66 - 0.68). Backscatter varies approximately as frequency cubed and trabecular thickness cubed in human calcaneus and femur in vitro. Backscatter from human calcaneus and bovine tibia exhibits substantial anisotropy. So far, backscatter has demonstrated only modest clinical utility. Computer simulation models have helped to elucidate mechanisms underlying scattering from cancellous bones.

Characterization of dense bovine cancellous bone tissue microstructure by ultrasonic backscattering using weak scattering models

A weak scattering model was proposed for the ultrasonic frequency-dependent backscatter in dense bovine cancellous bone, using two autocorrelation functions to describe the medium: one with discrete homogeneities (spherical distribution of equal spheres) and another, which considers tissue as an inhomogeneous continuum (densely populated medium). The inverse problem to estimate trabecular thickness of bone tissue has been addressed. A combination of the two autocorrelation functions was required to closely approximate the backscatter from bovine bone with various microarchitecture, given that the shape of trabeculae ranges from a rodlike to a platelike shape. Because of the large variation in trabecular thickness, both at an intraspecimen and an interspecimen level, thickness distributions for individual trabeculae for each bone specimen were obtained, and dominant trabecular sizes were determined. Comparison of backscatter measurements to theoretical predictions indicated that there were more than one dominant trabecular sizes that scatter sound for most specimens. Linear regression, performed between dominant trabecular thickness and estimated correlation length, showed significant linear correlation (R(2)=0.81). Attenuation due to scattering by a continuous distribution of scatterers was predicted to be linear over a frequency range from 0.3 to 0.9 MHz, suggesting a possibility that scattering may be a significant source of attenuation.

Ultrasonic tissue characterization via 2-D spectrum analysis: theory and in vitro measurements

A theoretical model is described for application in ultrasonic tissue characterization using a calibrated 2-D spectrum analysis method. This model relates 2-D spectra computed from ultrasonic backscatter signals to intrinsic physical properties of tissue microstructures, e.g., size, shape, and acoustic impedance. The model is applicable to most clinical diagnostic ultrasound systems. Two experiments employing two types of tissue architectures, spherical and cylindrical scatterers, are conducted using ultrasound with center frequencies of 10 and 40 MHz, respectively. Measurements of a tissue-mimicking phantom with an internal suspension of microscopic glass beads are used to validate the theoretical model. Results from in vitro muscle fibers are presented to further elucidate the utility of 2-D spectrum analysis in ultrasonic tissue characterization.

Evaluation of ultrasonic scattering in human cancellous bone by using a binary mixture model

A weak scattering model based on small perturbations in a binary mixture is developed to estimate the ultrasonic scattering from human cancellous bone, which is modelled as a random isotropic continuum containing identical scatters. Ultrasonic scattering is determined by both velocity fluctuation and density fluctuation, when k(2)a(2) << 1 is satisfied. Two kinds of trabeculae thickness distributions, i.e. even distribution and Gauss distribution, are applied in the calculation of attenuation and backscattering. Frequency dependence of the backscatter coefficient is found to be Af(3.13) and Af(2.84) with the Gauss distribution and an even distribution, respectively. Both backscattering and attenuation change significantly against porosity for the case of high porosity. The predicted results are close to the measured ones from the literature. The errors of this theoretical model are also discussed in this paper.

Ultrasonic characterization of human trabecular bone microstructure

New quantitative ultrasound (QUS) techniques involving ultrasound backscattering have been introduced for the assessment of bone quality. QUS parameters are affected by the transducer characteristics, e.g. frequency range, wave and pulse length. Although frequency-dependent backscattering has been studied extensively, understanding of the ultrasound scattering phenomenon in trabecular bone is still limited. In the present study, the relationships between QUS parameters and the microstructure of human trabecular bone were investigated experimentally and by using numerical simulations. Speed of sound (SOS), normalized broadband ultrasound attenuation (nBUA), average attenuation, integrated reflection coefficient (IRC) and broadband ultrasound backscatter (BUB) were measured for 26 human trabecular bone cylinders. Subsequently, a high-resolution microCT system was used to determine the microstructural parameters. Moreover, based on the sample-specific microCT data, a numerical model for ultrasound propagation was developed for the simulation of experimental measurements. Experimentally, significant relationships between the QUS parameters and microstructural parameters were demonstrated. The relationships were dependent on the frequency, and the strongest association (r = 0.88) between SOS and structural parameters was observed at a centre frequency of 5 MHz. nBUA, average attenuation, IRC and BUB showed somewhat lower linear correlations with the structural properties at a centre frequency of 5 MHz, as compared to those determined at lower frequencies. Multiple regression analyses revealed that the variation of acoustic parameters could best be explained by parameters reflecting the amount of mineralized tissue. A principal component analysis demonstrated that the strongest determinants of BUB and IRC were related to the trabecular structure. However, other structural characteristics contributed significantly to the prediction of the acoustic parameters as well. The two-dimensional numerical model introduced in the present study demonstrated good agreement with the experimental measurements. However, further studies with the simulation model are warranted to systematically investigate the relation between the structural parameters and ultrasound scattering.

The dependencies of phase velocity and dispersion on trabecular thickness and spacing in trabecular bone-mimicking phantoms

Frequency-dependent phase velocity was measured in trabecular-bone-mimicking phantoms consisting of two-dimensional arrays of parallel nylon wires (simulating trabeculae) with thicknesses ranging from 152 to 305 microm and spacings ranging from 700 to 1000 microm. Phase velocity varied approximately linearly with frequency over the range from 400 to 750 kHz. Dispersion was characterized by the slope of a linear least-squares regression fit to phase velocity versus frequency data. The increase in phase velocity (compared with that in water) at 500 kHz was approximately proportional to the (1) square of trabecular thickness, (2) inverse square of trabecular spacing, and (3) volume fraction occupied by nylon wires. The first derivative of phase velocity with respect to frequency was negative and exhibited nonlinear, monotonically decreasing dependencies on trabecular thickness and volume fraction. The dependencies of phase velocity and its first derivative on volume fraction in the phantoms were consistent with those reported in trabecular bone.

Comparison of measurements of phase velocity in human calcaneus to Biot theory

Biot's theory for elastic propagation in porous media has previously been shown to be useful for modeling the dependence of phase velocity on porosity in bovine cancellous bone in vitro. In the present study, Biot's theory is applied to measurements of porosity-dependent phase velocity in 53 human calcanea in vitro. Porosity was measured using microcomputed tomography for some samples (n = 23) and estimated based on bone mineral densitometry for the remaining samples (n = 30). The phase velocity at 500 kHz was measured in a water tank using a through-transmission technique. Biot's theory performed well for the prediction of the dependence of sound speed on porosity. The trend was quasilinear, but both the theory and experiment show similar slight curvature. The root mean square error (RMSE) of predicted versus measured sound speed was 15.8 m/s.

Prediction of density and mechanical properties of human trabecular bone in vitro by using ultrasound transmission and backscattering measurements at 0.2-6.7 MHz frequency range

The ultrasound (US) backscattering method has been introduced as an alternative for the through-transmission measurement of sound attenuation and speed in diagnosis of osteoporosis. Both attenuation and backscattering depend strongly on the US frequency. In this study, 20 human trabecular bone samples were measured in transmission and pulse-echo geometry in vitro. The aim of the study was to find the most sensitive frequency range for the quantitative ultrasound (QUS) analyses. Normalized broadband US attenuation (nBUA), speed of sound (SOS), broadband US backscatter (BUB) and integrated reflection coefficient (IRC) were determined for each sample. The samples were spatially scanned with five pairs of US transducers covering a frequency range of 0.2-6.7 MHz. Furthermore, mechanical properties and density of the same samples were determined. At all frequencies, SOS, BUB and IRC showed statistically significant linear correlations with the mechanical properties or density of human trabecular bone (0.51 < r < 0.82, 0.54 < r < 0.81 and 0.70 < r < 0.85, respectively). In contrast to SOS, IRC and BUB, nBUA showed statistically significant correlations with mechanical parameters or density at the centre frequency of 1 MHz only. Our results suggest that frequencies up to 5 MHz can be useful in QUS analyses for the prediction of bone mechanical properties and density. Since the use of higher frequencies provides better axial and spatial resolution, improved structural analyses may be possible. While extensive attenuation of high frequencies in trabecular bone limits the clinically feasible frequency range, selection of optimal frequency range for in vivo QUS application should be carefully considered.