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.

Development and characterization of a blood mimicking fluid for high intensity focused ultrasound

A blood mimicking fluid (BMF) has been developed for the acoustic and thermal characterizations of high intensity focused ultrasound (HIFU) ablation devices. The BMF is based on a degassed and de-ionized water solution dispersed with low density polyethylene microspheres, nylon particles, gellan gum, and glycerol. A broad range of physical parameters, including attenuation coefficient, speed of sound, viscosity, thermal conductivity, and diffusivity, were characterized as a function of temperature (20-70 degrees C). The nonlinear parameter B/A and backscatter coefficient were also measured at room temperature. Importantly, the attenuation coefficient is linearly proportional to the frequency (2-8 MHz) with a slope of about 0.2 dB cm(-1) MHz(-1) in the 20-70 degrees C range as in the case of human blood. Furthermore, sound speed and bloodlike backscattering indicate the usefulness of the BMF for ultrasound flow imaging and ultrasound-guided HIFU applications. Most of the other temperature-dependent physical parameters are also close to the reported values in human blood. These properties make it a unique HIFU research tool for developing standardized exposimetry techniques, validating numerical models, and determining the safety and efficacy of HIFU ablation devices.

Comparison of the Faran Cylinder Model and the Weak Scattering Model for predicting the frequency dependence of backscatter from human cancellous femur in vitro

This letter presents the first side-by-side comparison of the Faran Cylinder Model and the Weak Scattering Model for predicting backscatter from human femur. Both models are applied to the same dataset of frequency-dependent backscatter coefficients from 26 human femur cancellous bone samples in vitro. The Faran Cylinder Model predicts a slightly slower rate of increase of backscatter with frequency than the Weak Scattering Model, but both models are in reasonable agreement with the data and with each other, given the uncertainty in the measurements.

Interference between wave modes may contribute to the apparent negative dispersion observed in cancellous bone

Previous work has shown that ultrasonic waves propagating through cancellous bone often exhibit a linear-with-frequency attenuation coefficient, but a decrease in phase velocity with frequency (negative dispersion) that is inconsistent with the causality-imposed Kramers-Kronig relations. In the current study, interfering wave modes similar to those observed in bone are shown to potentially contribute to the observed negative dispersion. Biot theory, the modified Biot-Attenborogh model, and experimental results are used to aid in simulating multiple-mode wave propagation through cancellous bone. Simulations entail constructing individual wave modes exhibiting a positive dispersion using plausible velocities and amplitudes, and then summing the individual modes to create mixed-mode output wave forms. Results of the simulations indicate that mixed-mode wave forms can exhibit negative dispersion when analyzed conventionally under the assumption that only one wave is present, even when the individual interfering waves exhibit positive dispersions in accordance with the Kramers-Kronig relations. Furthermore, negative dispersion is observed when little or no visual evidence of interference exists in the time-domain data. Understanding the mechanisms responsible for the observed negative dispersion could aid in determining the true material properties of cancellous bone, as opposed to the apparent properties measured using conventional data analysis 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.

A method for improved standardization of in vivo calcaneal time-domain speed-of-sound measurements

Although calcaneal speed of sound (SOS) is an effective predictor of osteoporotic fracture risk, clinical SOS measurements exhibit a high degree of inter-system variability. Calcaneal SOS is usually computed from time-of-flight measurements of broadband ultrasound pulses that propagate through the foot. In order to minimize the effects of multi-path interference, many investigators measure time-of-flight from markers near the leading edge of the pulse. The calcaneus is a highly attenuating, highly inhomogeneous bone that distorts propagating ultrasound pulses via frequency-dependent attenuation, reverberation, dispersion, multiple scattering, and refraction. This pulse distortion can produce errors in leading-edge transit-time marker-based SOS measurements. In this paper, an equation to predict dependence of time-domain SOS measurements on system parameters (center frequency and bandwidth), transit-time marker location, and bone properties (attenuation coefficient and thickness) is validated with through-transmission measurements in a bone-mimicking phantom and in 73 women in vivo, using a clinical bone sonometer. In order to test the utility of the formula for suppressing system dependence of SOS measurements, a wideband laboratory data acquisition system was used to make a second set of through-transmission measurements on the phantom. The compensation formula reduced system-dependent leading-edge transit-time marker-based SOS measurements in the phantom from 41 m/s to 5 m/s and reduced average transit-time marker-related SOS variability in 73 women from 40 m/s to 10 m/s. The compensation formula can be used to improve standardization in bone sonometry.

The effect of phase cancellation on estimates of broadband ultrasound attenuation and backscatter coefficient in human calcaneus in vitro

Broadband ultrasound attenuation (BUA) is a clinically proven indicator of osteoporotic fracture risk. BUA measurements are typically performed in throughtransmission with single-element phase sensitive (PS) receivers and therefore can be compromised by phase cancellation artifact. Phase-insensitive (PI) receivers suppress phase cancellation artifact. To study the effect of phase cancellation on BUA measurements, through-transmission measurements were performed on 16 human calcaneus samples in vitro using a two-dimensional receiver array that enabled PS and PI BUA estimation. The means plus or minus standard deviations for BUA measurements were 22.1 +/- 15.8 dB/MHz (PS) and 17.6 +/- 7.2 dB/MHz (PI), suggesting that, on the average, approximately 20% of PS BUA values in vitro can be attributed to phase cancellation artifact. Therefore, although cortical plates are often regarded as the primary source of phase cancellation artifact, the heterogeneity of cancellous bone in the calcaneal interior may also be a significant source. Backscatter coefficient estimates in human calcaneus that are based on PS attenuation compensation overestimate 1) average magnitude of backscatter coefficient at 500 kHz by a factor of about 1.6 +/- 0.3 and 2) average exponent (n) of frequency dependence by about 0.34 +/- 0.12 (where backscatter coefficient is fit to a power law form proportional to frequency to the nth power).

Introduction to the special issue on diagnostic and therapeutic applications of ultrasound in bone--Part I

The 13 papers and one item of correspondence in this special issue focus on diagnostic and therapeutic applications of ultrasound in bone. In Part I manuscripts discuss instrumentation, numerical modeling, applications, and guided waves.

The dependence of time-domain speed-of-sound measurements on center frequency, bandwidth, and transit-time marker in human calcaneus in vitro

Time-domain speed-of-sound (SOS) measurements in calcaneus are effective predictors of osteoporotic fracture risk. High attenuation and dispersion in bone, however, produce severe distortion of transmitted pulses that leads to ambiguity of time-domain SOS measurements. An equation to predict the effects of system parameters (center frequency and bandwidth), algorithm parameters (pulse arrival-time marker), and bone properties (attenuation coefficient and thickness) on time-domain SOS estimates is derived for media with attenuation that varies linearly with frequency. The equation is validated using data from a bone-mimicking phantom and from 30 human calcaneus samples in vitro. The data suggest that the effects of dispersion are small compared with the effects of frequency-dependent attenuation. The equation can be used to retroactively compensate data. System-related variations in SOS are shown to decrease as the pulse-arrival-time marker is moved toward the pulse center. Therefore, compared with other time-domain measures of SOS, group velocity exhibits the minimum system dependence.

The effect of phase cancellation on estimates of calcaneal broadband ultrasound attenuation in vivo

Broadband ultrasonic attenuation (BUA) is a clinically-accepted measurement for prediction of osteoporotic fracture risk. Typical clinical BUA measurements are performed with phase-sensitive receivers and, therefore, can be affected by phase cancellation. In order to separate the effects of conventional attenuation (absorption plus scattering) from phase cancellation, BUA was measured on phantoms with acrylic wedge phase aberrators and on 73 women using both phase sensitive (PS) and phase insensitive (PI) reception. A clinical bone sonometer with a two-dimensional (2-D) receiver array was used. PI BUA measurements on phantoms with acrylic wedge phase aberrators were found to be far more resistant to phase cancellation than PS BUA measurements. In data from 73 women, means and standard deviations for BUA measurements were 81.4 +/- 21.4 dB/MHz (PS) and 67.2 +/- 9.7 dB/MHz (PI). The magnitude of the discrepancy between PS BUA and PI BUA tended to increase with bone mineral density (BMD).

Group velocity, phase velocity, and dispersion in human calcaneus in vivo

Commercial bone sonometers measure broadband ultrasonic attenuation and/or speed of sound (SOS) in order to assess bone status. Phase velocity, which is usually measured in frequency domain, is a fundamental material property of bone that is related to SOS, which is usually measured in time domain. Four previous in vitro studies indicate that phase velocity in human cancellous bone decreases with frequency (i.e., negative dispersion). In order to investigate frequency-dependent phase velocity in vivo, through-transmission measurements were performed in 73 women using a GE Lunar Achilles Insight commercial bone sonometer. Average phase velocity at 500 kHz was 1489 +/- 55 m/s (mean +/- standard deviation). Average dispersion rate was -59 +/- 52 m/sMHz. Group velocity was usually lower than phase velocity, as is expected for negatively dispersive media. Using a stratified model to represent cancellous bone, the reductions in phase velocity and dispersion rate in vivo as opposed to in vitro can be explained by (1) the presence of marrow instead of water as a fluid filler, and (2) the decreased porosity of bones of living (compared with deceased) subjects.

Interlaboratory comparison of ultrasonic backscatter coefficient measurements from 2 to 9 MHz

OBJECTIVE

As are the attenuation coefficient and sound speed, the backscatter coefficient is a fundamental ultrasonic property that has been used to characterize many tissues. Unfortunately, there is currently far less standardization for the ultrasonic backscatter measurement than for the other two, as evidenced by a previous American Institute of Ultrasound in Medicine (AIUM)-sponsored interlaboratory comparison of ultrasonic backscatter, attenuation, and speed measurements (J Ultrasound Med 1999; 18:615-631). To explore reasons for these disparities, the AIUM Endowment for Education and Research recently supported this second interlaboratory comparison, which extends the upper limit of the frequency range from 7 to 9 MHz.

METHODS

Eleven laboratories were provided with standard test objects designed and manufactured at the University of Wisconsin (Madison, WI). Each laboratory was asked to perform ultrasonic measurements of sound speed, attenuation coefficients, and backscatter coefficients. Each laboratory was blinded to the values of the ultrasonic properties of the test objects at the time the measurements were performed.

RESULTS

Eight of the 11 laboratories submitted results. The range of variation of absolute magnitude of backscatter coefficient measurements was about 2 orders of magnitude. If the results of 1 outlier laboratory are excluded, then the range is reduced to about 1 order of magnitude. Agreement regarding frequency dependence of backscatter was better than reported in the previous interlaboratory comparison. For example, when scatterers were small compared with the ultrasonic wavelength, experimental frequency-dependent backscatter coefficient data obtained by the participating laboratories were usually consistent with the expected Rayleigh scattering behavior (proportional to frequency to the fourth power).

CONCLUSIONS

Greater standardization of backscatter measurement methods is needed. Measurements of frequency dependence of backscatter are more consistent than measurements of absolute magnitude.

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.

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

A theory for the elastic scattering response from a cylinder insonified by a plane wave was previously derived by Faran. In the present paper, the empirical relationship between Faran's theory and measurements of backscatter coefficient from cylindrical targets using focused transducers is investigated. Experimental measurements of dependence of backscatter coefficient on frequency and diameter for nylon wires are reported. It is found that, under certain conditions (including weak, incoherent scattering), backscatter coefficient measurements from collections of cylindrical scatterers may be meaningfully compared with Faran's model predictions. At low frequencies, the theory and experimental measurements exhibit similar dependences on frequency and diameter, provided that the scatterers are not too densely packed. At higher frequencies, the fine structure of Faran's predictions becomes difficult to reproduce experimentally with a focused transducer. Implications regarding applications to characterization of trabecular bone are discussed.

The dependence of ultrasonic backscatter on trabecular thickness in human calcaneus: theoretical and experimental results

Trabecular thickness within cancellous bone is an important determinant of osteoporotic fracture risk. Noninvasive assessment of trabecular thickness potentially could yield useful diagnostic information. Faran's theory of elastic scattering from a cylindrical object immersed in a fluid has been used to predict the dependence of ultrasonic backscatter on trabecular thickness. The theory predicts that, in the range of morphological and material properties expected for trabecular bone, the backscatter coefficient at 500 kHz should be approximately proportional to trabecular thickness to the power of 2.9. Experimental measurements of backscatter coefficient were performed on 43 human calcaneus samples in vitro. Mean trabecular thicknesses on the 43 samples were assessed using micro computed tomography (CT). A power law fit to the data showed that the backscatter coefficient empirically varied as trabecular thickness to the 2.8 power. The 95% confidence interval for this exponent was 1.7 to 3.9. The square of the correlation coefficient for the linear regression to the log transformed data was 0.40. This suggests that 40% of variations in backscatter may be attributed to variations in trabecular thickness. These results reinforce previous studies that offered validation for the Faran cylinder model for prediction of scattering properties of cancellous bone, and provide added evidence for the potential diagnostic utility of the backscatter measurement.

The effect of trabecular material properties on the frequency dependence of backscatter from cancellous bone

Previous experimental studies indicate that backscatter coefficient for human calcaneal trabecular bone varies approximately as frequency cubed. This frequency dependence has been shown to be consistent with a model in which trabeculae are thought of as long thin cylinders composed of a substance with the same material properties as hydroxyapatite. The true material properties of human trabecular bone are not known however. Based on reported measurements of material properties of many bones and bone-like substances, it is possible that the density and longitudinal sound speed of trabecular bone material are far lower than the hydroxyapatite model would suggest. In this letter, it is shown that the frequency dependence of backscatter is still expected to be approximately cubic for wide ranges for density and longitudinal sound speed (encompassing the conceivable ranges for trabecular bone).

Autocorrelation and cepstral methods for measurement of tibial cortical thickness

Cortical thickness of the tibia is related to stress fracture risk and overall skeletal status. Two methods are proposed for estimating tibial cortical thickness based on power spectra of ultrasonic echoes containing reflections from front and back surfaces of the cortex. The locations of the peaks in the autocorrelation function and the cepstrum are related to cortical thickness. Data were acquired on plastic plates in order to validate the methodology. These data indicate high correlations between estimated and true thickness with correlation coefficients r = 0.99, (95% confidence interval: 0.993-1.00) for the autocorrelation method and r = 0.99 (95% CI: 0.996-1.00) for the cepstral method. Data on six tibia samples in vitro indicate correlation coefficients of r = 0.92 (95% CI: 0.72-1.00) for the autocorrelation method and r = 0.85 (95% CI: 0.62-0.94) for the cepstral method. Estimates of precisions of the two methods were 0.3 +/- 0.1 mm (autocorrelation method) and 0.5 +/- 0.2 mm (cepstral method). One measurement in a human volunteer in vivo demonstrated clinical feasibility of the measurement and good agreement with cortical thickness assessed using peripheral quantitative computed tomography (QCT). This technology offers the promise of an inexpensive, fast, portable, simple, nonionizing technique for assessing skeletal status.

Characterization of trabecular bone using the backscattered spectral centroid shift

Ultrasonic attenuation in bone in vivo is generally measured using a through-transmission method at the calcaneus. Although attenuation in calcaneus has been demonstrated to be a useful predictor for osteoporotic fracture risk, measurements at other clinically important sites, such as hip and spine, could potentially contain additional useful diagnostic information. Through-transmission measurements may not be feasible at these sites due to complex bone shapes and the increased amount of intervening soft tissue. Centroid shift from the backscattered signal is an index of attenuation slope and has been used previously to characterize soft tissues. In this paper, centroid shift from signals backscattered from 30 trabecular bone samples in vitro were measured. Attenuation slope also was measured using a through-transmission method. The correlation coefficient between centroid shift and attenuation slope was -0.71. The 95% confidence interval was (-0.86, -0.47). These results suggest that the backscattered spectral centroid shift may contain useful diagnostic information potentially applicable to hip and spine.

A Gaussian framework for modeling effects of frequency-dependent attenuation, frequency-dependent scattering, and gating

For a wide range of applications in medical ultrasound, power spectra of received signals are approximately Gaussian. It has been established previously that an ultrasound beam with a Gaussian spectrum propagating through a medium with linear attenuation remains Gaussian. In this paper, Gaussian transformations are derived to model the effects of scattering (according to a power law, as is commonly applicable in soft tissues, especially over limited frequency ranges) and gating (with a Hamming window, a commonly used gate function). These approximations are shown to be quite accurate even for relatively broad band systems with fractional bandwidths approaching 100%. The theory is validated by experiments in phantoms consisting of glass particles suspended in agar.

Fundamental precision limitations for measurements of frequency dependence of backscatter: applications in tissue-mimicking phantoms and trabecular bone

Various models for ultrasonic scattering from trabecular bone have been proposed. They may be evaluated to a certain extent by comparison with experimental measurements. In order to appreciate limitations of these comparisons, it is important to understand measurement precision. In this article, an approach proposed by Lizzi and co-workers is adapted to model precision of estimates of frequency-dependent backscatter for scattering targets (such as trabecular bone) that contain many scatterers per resolution cell. This approach predicts uncertainties in backscatter due to the random nature of the interference of echoes from individual scatterers as they are summed at the receiver. The model is validated in experiments on a soft-tissue-mimicking phantom and on 24 human calcaneus samples interrogated in vitro. It is found that while random interference effects only partially explain measured variations in the magnitude of backscatter, they are virtually entirely responsible for observed variations in the frequency dependence (exponent of a power law fit) of backscatter.