Negative dispersion in bone: the role of interference in measurements of the apparent phase velocity of two temporally overlapping signals

Measurement of Ultrasound Parameters of Bovine Cancellous Bone as a Function of Frequency for a Range of Porosities via Through-Transmission Ultrasonic Spectroscopy

The relationship between ultrasonic parameters (attenuation coefficients and velocity) and bone porosity in bovine cancellous bone is explored to understand the possibility of fracture risk diagnosis associated with osteoporosis by applying ultrasound. In vitro measurements of ultrasonic parameters on twenty-one bovine cancellous bone samples from tibia were conducted, using ultrasonic spectroscopy in the through-transmission mode. Transducers of three different center frequencies were used to cover a wide diagnostic frequency range between 1.0–7.8 MHz. The nonlinear relationship of porosity and normalized attenuation coefficient (nATTN) and normalized broadband attenuation coefficient (nBUA) were well described by a third-order polynomial fit, whereas porosity and the phase velocity (UV) were found to be negatively correlated with the linear correlation coefficients of −0.93, −0.89 and −0.83 at 2.25, 5.00 and 7.50 MHz, respectively. The results imply that the ultrasound parameters attain maximum values for the bone sample with the lowest porosity, and then decrease for samples with greater porosity for the range of porosities in our samples for all frequencies. Spatial variation in the ultrasound parameters was found to be caused by non-uniform pore size distribution, which was examined at five different locations within the same bone specimen. However, it did not affect the relationship of ultrasound parameters and porosity at these frequencies.

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.

A μCT-based investigation of the influence of tissue modulus variation, anisotropy and inhomogeneity on ultrasound propagation in trabecular bone

Ultrasound propagation is widely used in the diagnosis of osteoporosis by providing information on bone mechanical quality. When it loses calcium, the tissue properties will first decrease. However, limited research about the influence of tissue properties on ultrasound propagation have been done due to the cumbersome experiment. The goal of this study was to explore the relationships between tissue modulus (Es) and speed of sound (SOS) through numerical simulations, and to study the influence of Es on the acoustical behavior in characterizing the local structural anisotropy and inhomogeneity. In this work, three-dimensional finite element (FE) simulations were performed on a cubic high-resolution (15μm) bovine trabecular bone sample (4×4×4mm(3), BV/TV=0.18) mapped from micro-computed tomography. Ultrasound excitations of 50kHz, 500kHz and 2MHz were applied in three orthogonal axes and the first arriving signal (FAS) was collected to quantify wave velocity. In this study, a strong power law relationship between Es and SOS was measured with estimated exponential index β=2.08-3.44 for proximal-distal (PD), anterior-posterior (AP) and medial-lateral (ML), respectively (all R(2)>0.95). For various Es, a positive dispersion of sound speed with respect to sound frequency was observed and the velocity dispersion magnitude (VDM) was measured. Also, with Es=15GPa in three orientations, the SOS in PD axis is 2009±120m/s, faster than that of AP (1762±106m/s) and ML (1798±132m/s) (f=2MHz) directions. Besides, the standard deviation of SOS increases with the sound frequency and the Es in all directions except for that at 50kHz. For the mechanical properties, the apparent modulus with certain Es was highest in the longitudinal direction compared with the transverse directions. It indicates that the tissue modulus combining with anisotropy and inhomogeneity has great influence on ultrasound propagation. Simulation results agree well with theoretical and experimental results.

Solid volume fraction estimation of bone:marrow replica models using ultrasound transit time spectroscopy

The acceptance of broadband ultrasound attenuation (BUA) for the assessment of osteoporosis suffers from a limited understanding of both ultrasound wave propagation through cancellous bone and its exact dependence upon the material and structural properties. It has recently been proposed that ultrasound wave propagation in cancellous bone may be described by a concept of parallel sonic rays; the transit time of each ray defined by the proportion of bone and marrow propagated. A Transit Time Spectrum (TTS) describes the proportion of sonic rays having a particular transit time, effectively describing the lateral inhomogeneity of transit times over the surface aperture of the receive ultrasound transducer. The aim of this study was to test the hypothesis that the solid volume fraction (SVF) of simplified bone:marrow replica models may be reliably estimated from the corresponding ultrasound transit time spectrum. Transit time spectra were derived via digital deconvolution of the experimentally measured input and output ultrasonic signals, and compared to predicted TTS based on the parallel sonic ray concept, demonstrating agreement in both position and amplitude of spectral peaks. Solid volume fraction was calculated from the TTS; agreement between true (geometric calculation) with predicted (computer simulation) and experimentally-derived values were R(2)=99.9% and R(2)=97.3% respectively. It is therefore envisaged that ultrasound transit time spectroscopy (UTTS) offers the potential to reliably estimate bone mineral density and hence the established T-score parameter for clinical osteoporosis assessment.

Conventional, Bayesian, and Modified Prony's methods for characterizing fast and slow waves in equine cancellous bone

Conventional, Bayesian, and the modified least-squares Prony's plus curve-fitting (MLSP + CF) methods were applied to data acquired using 1 MHz center frequency, broadband transducers on a single equine cancellous bone specimen that was systematically shortened from 11.8 mm down to 0.5 mm for a total of 24 sample thicknesses. Due to overlapping fast and slow waves, conventional analysis methods were restricted to data from sample thicknesses ranging from 11.8 mm to 6.0 mm. In contrast, Bayesian and MLSP + CF methods successfully separated fast and slow waves and provided reliable estimates of the ultrasonic properties of fast and slow waves for sample thicknesses ranging from 11.8 mm down to 3.5 mm. Comparisons of the three methods were carried out for phase velocity at the center frequency and the slope of the attenuation coefficient for the fast and slow waves. Good agreement among the three methods was also observed for average signal loss at the center frequency. The Bayesian and MLSP + CF approaches were able to separate the fast and slow waves and provide good estimates of the fast and slow wave properties even when the two wave modes overlapped in both time and frequency domains making conventional analysis methods unreliable.

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.

Application of an effective medium theory for modeling ultrasound wave propagation in healing long bones

Quantitative ultrasound has recently drawn significant interest in the monitoring of the bone healing process. Several research groups have studied ultrasound propagation in healing bones numerically, assuming callus to be a homogeneous and isotropic medium, thus neglecting the multiple scattering phenomena that occur due to the porous nature of callus. In this study, we model ultrasound wave propagation in healing long bones using an iterative effective medium approximation (IEMA), which has been shown to be significantly accurate for highly concentrated elastic mixtures. First, the effectiveness of IEMA in bone characterization is examined: (a) by comparing the theoretical phase velocities with experimental measurements in cancellous bone mimicking phantoms, and (b) by simulating wave propagation in complex healing bone geometries by using IEMA. The original material properties of cortical bone and callus were derived using serial scanning acoustic microscopy (SAM) images from previous animal studies. Guided wave analysis is performed for different healing stages and the results clearly indicate that IEMA predictions could provide supplementary information for bone assessment during the healing process. This methodology could potentially be applied in numerical studies dealing with wave propagation in composite media such as healing or osteoporotic bones in order to reduce the simulation time and simplify the study of complicated geometries with a significant porous nature.

Experimental and computer simulation validation of ultrasound phase interference created by lateral inhomogeneity of transit time in replica bone: marrow composite models

The measurement of broadband ultrasound attenuation in cancellous bone for the assessment of osteoporosis follows a parabolic-type dependence with bone volume fraction, having minima values corresponding to both entire bone and entire marrow. Langton has recently proposed that the primary attenuation mechanism is phase interference due to variations in propagation transit time through the test sample as detected over the phase-sensitive surface of the receive ultrasound transducer. This fundamentally simple concept assumes that the propagation may be considered as an array of parallel ‘sonic rays’. The transit time of each ray is defined by the proportion of bone and marrow propagated, being a minimum (tmin) solely through bone and a maximum (tmax) solely through marrow, from which a transit time spectrum, may be defined describing the proportion of sonic rays having a particular transit time.

The aim of this study was to test the hypothesis that there is a dependence of phase interference upon the lateral inhomogeneity of transit time by comparing experimental measurements and computer simulation predictions of ultrasound propagation through a range of relatively simplistic solid:liquid models.

From qualitative and quantitative comparison of the experimental and computer simulation results, there is an extremely high degree of agreement of 94.2%–99.0% between the two approaches. This combined experimental and computer simulation study has successfully demonstrated that lateral inhomogeneity of transit time has significant potential for phase interference to occur if a phase-sensitive receive ultrasound transducer is implemented as in most commercial ultrasound bone analysis devices.

Phase velocity and attenuation predictions of waves in cancellous bone using an iterative effective medium approximation

The quantitative determination of wave dispersion and attenuation in bone is an open research area as the factors responsible for ultrasound absorption and scattering in composite biological tissues have not been completely explained. In this study, we use the iterative effective medium approximation (IEMA) proposed in [1] so as to calculate phase velocity and attenuation in media with properties similar to those of cancellous bones. Calculations are performed for a frequency range of 0.4-0.8 MHz and for different inclusions' volume concentrations and sizes. Our numerical results are compared with previous experimental findings so as to assess the effectiveness of IEMA. It was made clear that attenuation and phase velocity estimations could provide supplementary information for cancellous bone characterization.

Cancellous bone fast and slow waves obtained with Bayesian probability theory correlate with porosity from computed tomography

A Bayesian probability theory approach for separating overlapping ultrasonic fast and slow waves in cancellous bone has been previously introduced. The goals of this study were to investigate whether the fast and slow waves obtained from Bayesian separation of an apparently single mode signal individually correlate with porosity and to isolate the fast and slow waves from medial-lateral insonification of the calcaneus. The Bayesian technique was applied to trabecular bone data from eight human calcanei insonified in the medial-lateral direction. The phase velocity, slope of attenuation (nBUA), and amplitude were determined for both the fast and slow waves. The porosity was assessed by micro-computed tomography (microCT) and ranged from 78.7% to 94.1%. The method successfully separated the fast and slow waves from medial-lateral insonification of the calcaneus. The phase velocity for both the fast and slow wave modes showed an inverse correlation with porosity (R(2) = 0.73 and R(2) = 0.86, respectively). The slope of attenuation for both wave modes also had a negative correlation with porosity (fast wave: R(2) = 0.73, slow wave: R(2) = 0.53). The fast wave amplitude decreased with increasing porosity (R(2) = 0.66). Conversely, the slow wave amplitude modestly increased with increasing porosity (R(2) = 0.39).

Determining attenuation properties of interfering fast and slow ultrasonic waves in cancellous bone

Previous studies have shown that interference between fast waves and slow waves can lead to observed negative dispersion in cancellous bone. In this study, the effects of overlapping fast and slow waves on measurements of the apparent attenuation as a function of propagation distance are investigated along with methods of analysis used to determine the attenuation properties. Two methods are applied to simulated data that were generated based on experimentally acquired signals taken from a bovine specimen. The first method uses a time-domain approach that was dictated by constraints imposed by the partial overlap of fast and slow waves. The second method uses a frequency-domain log-spectral subtraction technique on the separated fast and slow waves. Applying the time-domain analysis to the broadband data yields apparent attenuation behavior that is larger in the early stages of propagation and decreases as the wave travels deeper. In contrast, performing frequency-domain analysis on the separated fast waves and slow waves results in attenuation coefficients that are independent of propagation distance. Results suggest that features arising from the analysis of overlapping two-mode data may represent an alternate explanation for the previously reported apparent dependence on propagation distance of the attenuation coefficient of cancellous bone.

Quantitative photoacoustic imaging: correcting for heterogeneous light fluence distributions using diffuse optical tomography

The specificity of molecular and functional photoacoustic (PA) images depends on the accuracy of the photoacoustic absorption spectroscopy. The PA signal is proportional to the product of the optical absorption coefficient and local light fluence; quantitative PA measurements of the optical absorption coefficient therefore require an accurate estimation of optical fluence. Light-modeling aided by diffuse optical tomography (DOT) can be used to map the required fluence and to reduce errors in traditional PA spectroscopic analysis. As a proof-of-concept, we designed a tissue-mimicking phantom to demonstrate how fluence-related artifacts in PA images can lead to misrepresentations of tissue properties. To correct for these inaccuracies, the internal fluence in the tissue phantom was estimated by using DOT to reconstruct spatial distributions of the absorption and reduced scattering coefficients of multiple targets within the phantom. The derived fluence map, which only consisted of low spatial frequency components, was used to correct PA images of the phantom. Once calibrated to a known absorber, this method reduced errors in estimated absorption coefficients from 33% to 6%. These results experimentally demonstrate that combining DOT with PA imaging can significantly reduce fluence-related errors in PA images, while producing quantitatively accurate, high-resolution images of the optical absorption coefficient.

Measurements of ultrasonic phase velocities and attenuation of slow waves in cellular aluminum foams as cancellous bone-mimicking phantoms

The water-saturated aluminum foams with an open network of interconnected ligaments were investigated by ultrasonic transmission technique for the suitability as cancellous bone-mimicking phantoms. The phase velocities and attenuation of nine samples covering three pores per inch (5, 10, and 20 PPI) and three aluminum volume fractions (5, 8, and 12% AVF) were measured over a frequency range of 0.7-1.3 MHz. The ligament thickness and pore sizes of the phantoms and low-density human cancellous bones are similar. A strong slow wave and a weak fast wave are observed for all samples while the latter is not visible without significant amplification (100x). This study reports the characteristics of slow wave, whose speeds are less than the sound speed of the saturating water and decrease mildly with AVF and PPI with an average 1469 m/s. Seven out of nine samples show positive dispersion and the rest show minor negative dispersion. Attenuation increases with AVF, PPI, and frequency except for the 20 PPI samples, which exhibit non-increasing attenuation level with fluctuations due to scattering. The phase velocities agree with Biot's porous medium theory. The RMSE is 16.0 m/s (1%) at n = 1.5. Below and above this value, the RMSE decreases mildly and rises sharply, respectively.

Multi-frequency characterization of the speed of sound and attenuation coefficient for longitudinal transmission of freshly excised human skulls

For medical applications of ultrasound inside the brain, it is necessary to understand the relationship between the apparent density of skull bone and its corresponding speed of sound and attenuation coefficient. Although there have been previous studies exploring this phenomenon, there is still a need to extend the measurements to cover more of the clinically relevant frequency range. The results of measurements of the longitudinal speed of sound and attenuation coefficient are presented for specimens of human calvaria. The study was performed for the frequencies of 0.27, 0.836, 1.402, 1.965 and 2.525 MHz. Specimens were obtained from fresh cadavers through a protocol with the Division of Anatomy of the University of Toronto. The protocol was approved by the Research Ethics Board of Sunnybrook Health Sciences Centre. The specimens were mounted in polycarbonate supports that were marked for stereoscopic positioning. Computer tomography (CT) scans of the skulls mounted on their supports were performed, and a three-dimensional skull surface was reconstructed. This surface was used to guide a positioning system to ensure the normal sound incidence of an acoustic signal. This signal was produced by a focused device with a diameter of 5 cm and a focal length of 10 cm. Measurements of delay in time of flight were carried out using a needle hydrophone. Measurements of effective transmitted energy were carried out using a radiation force method with a 10 µg resolution scale. Preliminary functions of speed of sound and attenuation coefficient, both of which are related to apparent density, were established using a multi-layer propagation model that takes into account speed of sound, density and thickness of the layer. An optimization process was executed from a large set of random functions and the best functions were chosen for those ones that closest reproduced the experimental observations. The final functions were obtained after a second pass of the optimization process was executed, but this time using a finite-difference time-difference solution of the Westervelt equation, which is more precise than the multi-layer model but much more time consuming for computation. For six of seven specimens, measurements were carried out on five locations on the calvaria, and for the other specimen three measurements were made. In total, measurements were carried out on 33 locations. Results indicated the presence of dispersion effects and that these effects are different according to the type of bone in the skull (cortical and trabecular). Additionally, both the speed of sound and attenuation showed dependence on the skull density that varied with the frequency. Using the optimal functions and the information of density from the CT scans, the average values (±s.d.) of the speed of sound for cortical bone were estimated to be 2384(± 130), 2471(± 90), 2504(± 120), 2327(± 90) and 2053(± 40) m s(-1) for the frequencies of 270, 836, 1402, 1965 and 2526 kHz, respectively. For trabecular bone, and in the same order of frequency values, the speeds of sound were 2140(± 130), 2300(± 100), 2219(± 200), 2133(± 130) and 1937(± 40) m s(-1), respectively. The average values of the attenuation coefficient for cortical bone were 33(± 9), 240(± 9) and 307(± 30) Np m(-1) for the frequencies of 270, 836, and 1402, respectively. For trabecular bone, and in the same order of frequency values, the average values of the attenuation coefficient were 34(± 13), 216(± 16) and 375(± 30) Np m(-1), respectively. For frequencies of 1.965 and 2.525 MHz, no measurable radiation force was detected with the setup used.

Inverse problems in cancellous bone: estimation of the ultrasonic properties of fast and slow waves using Bayesian probability theory

Quantitative ultrasonic characterization of cancellous bone can be complicated by artifacts introduced by analyzing acquired data consisting of two propagating waves (a fast wave and a slow wave) as if only one wave were present. Recovering the ultrasonic properties of overlapping fast and slow waves could therefore lead to enhancement of bone quality assessment. The current study uses Bayesian probability theory to estimate phase velocity and normalized broadband ultrasonic attenuation (nBUA) parameters in a model of fast and slow wave propagation. Calculations are carried out using Markov chain Monte Carlo with simulated annealing to approximate the marginal posterior probability densities for parameters in the model. The technique is applied to simulated data, to data acquired on two phantoms capable of generating two waves in acquired signals, and to data acquired on a human femur condyle specimen. The models are in good agreement with both the simulated and experimental data, and the values of the estimated ultrasonic parameters fall within expected ranges.

Cancellous bone analysis with modified least squares Prony's method and chirp filter: phantom experiments and simulation

The presence of two longitudinal waves in porous media is predicted by Biot's theory and has been confirmed experimentally in cancellous bone. When cancellous bone samples are interrogated in through-transmission, these two waves can overlap in time. Previously, the Modified Least-Squares Prony's (MLSP) method was validated for estimation of amplitudes, attenuation coefficients, and phase velocities of fast and slow waves, but tended to overestimate phase velocities by up to about 5%. In the present paper, a pre-processing chirp filter to mitigate the phase velocity bias is derived. The MLSP/chirp filter (MLSPCF) method was tested for decomposition of a 500 kHz-center-frequency signal containing two overlapping components: one passing through a low-density-polyethylene plate (fast wave) and another passing through a cancellous-bone-mimicking phantom material (slow wave). The chirp filter reduced phase velocity bias from 100 m/s (5.1%) to 69 m/s (3.5%) (fast wave) and from 29 m/s (1.9%) to 10 m/s (0.7%) (slow wave). Similar improvements were found for 1) measurements in polycarbonate (fast wave) and a cancellous-bone-mimicking phantom (slow wave), and 2) a simulation based on parameters mimicking bovine cancellous bone. The MLSPCF method did not offer consistent improvement in estimates of attenuation coefficient or amplitude.

Decomposition of two-component ultrasound pulses in cancellous bone using modified least squares prony method--phantom experiment and simulation

Porous media such as cancellous bone often support the simultaneous propagation of two compressional waves. When small bone samples are interrogated in through-transmission with broadband sources, these two waves often overlap in time. The modified least-squares Prony's (MLSP) method was tested for decomposing a 500 kHz-center-frequency signal containing two overlapping components: one passing through a polycarbonate plate (to produce the "fast" wave) and another passing through a cancellous-bone-mimicking phantom (to produce the "slow" wave). The MLSP method yielded estimates of attenuation slopes accurate to within 7% (polycarbonate plate) and 2% (cancellous bone phantom). The MLSP method yielded estimates of phase velocities accurate to within 1.5% (both media). The MLSP method was also tested on simulated data generated using attenuation slopes and phase velocities corresponding to bovine cancellous bone. Throughout broad ranges of signal-to-noise ratio (SNR), the MLSP method yielded estimates of attenuation slope that were accurate to within 1.0% and estimates of phase velocity that were accurate to within 4.3% (fast wave) and 1.3% (slow wave).

Simplified expressions of the subtracted Kramers-Kronig relations using the expanded forms applied to ultrasonic power-law systems

The Kramers-Kronig (KK) relations are a large class of integral transformations that exploit the broad principle of simple causality in order to link the physical properties of matter and materials. In applications to the complex-valued wavenumber for acoustic propagation, the method of subtractions is used to form convergent integral relations between the phase velocity and the attenuation coefficient. When the method of subtractions is applied in the usual manner, the integrands in the relations become unnecessarily complicated. In this work, an expanded form of the subtracted relations is presented, which is essentially a truncated Taylor series expansion of the Hilbert transforms. The implementation of the relations only requires the explicit evaluation of two simply expressed integrals involving the Hilbert transform kernel. These two integrals determine the values of the other terms in the subtracted relations, demonstrating the computational efficiency of the technique. The method is illustrated analytically through its application to power-law attenuation coefficients and its associated dispersion, which are observed in a wide variety of materials. This approach explicitly shows the central role of the Hilbert transform kernel in the KK relations, which can become obscured in other formulations.

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.

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.

Bone sonometry: reducing phase aberration to improve estimates of broadband ultrasonic attenuation

Previous studies suggest that phase cancellation at the receiving transducer can result in the overestimation of the frequency dependent ultrasonic attenuation of bone, a quantity that has been shown to correlate with bone mineral density and ultimately with osteoporotic fracture risk. Evidence supporting this interpretation is provided by phase insensitive processing of the data, which appear to reduce the apparent overestimates of attenuation. The present study was designed to clarify the components underlying phase aberration artifacts in such through-transmission measurements by conducting systematic studies of the simplest possible test objects capable of introducing phase aberration. Experimental results are presented for a Lexan phantom over the frequency range 300-700 kHz and a Plexiglas phantom over the 3-7 MHz range. Both phantoms were flat and parallel plates featuring a step discontinuity milled into one of their initially flat sides. The through-transmitted signals were received by a 0.6 mm diameter membrane hydrophone that was raster scanned over a grid coaxial with the transmitting transducer. Signals received by the pseudoarray were processed offline to emulate phase sensitive and phase insensitive receivers with different aperture diameters. The data processed phase sensitively were focused to demonstrate the results of planar, geometrical, and correlation-based aberration correction methods. Results are presented illustrating the relative roles of interference in the ultrasonic field and phase cancellation at the receiving transducer in producing phase aberration artifacts. It was found that artifacts due to phase cancellation or interference can only be minimized with phase insensitive summation techniques by choosing an appropriately large receiving aperture. Data also suggest the potentially confounding role of time-and frequency-domain artifacts on ultrasonic measurements and illustrate the advantages of two-dimensional receiving arrays in determining the slope of attenuation (nBUA) for the clinical assessment of osteoporosis.

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.