Relationships among calcaneal backscatter, attenuation, sound speed, hip bone mineral density, and age in normal adult women

The Biomechanics of Musculoskeletal Tissues during Activities of Daily Living: Dynamic Assessment Using Quantitative Transmission-Mode Ultrasound Techniques

The measurement of musculoskeletal tissue properties and loading patterns during physical activity is important for understanding the adaptation mechanisms of tissues such as bone, tendon, and muscle tissues, particularly with injury and repair. Although the properties and loading of these connective tissues have been quantified using direct measurement techniques, these methods are highly invasive and often prevent or interfere with normal activity patterns. Indirect biomechanical methods, such as estimates based on electromyography, ultrasound, and inverse dynamics, are used more widely but are known to yield different parameter values than direct measurements. Through a series of literature searches of electronic databases, including Pubmed, Embase, Web of Science, and IEEE Explore, this paper reviews current methods used for the in vivo measurement of human musculoskeletal tissue and describes the operating principals, application, and emerging research findings gained from the use of quantitative transmission-mode ultrasound measurement techniques to non-invasively characterize human bone, tendon, and muscle properties at rest and during activities of daily living. In contrast to standard ultrasound imaging approaches, these techniques assess the interaction between ultrasound compression waves and connective tissues to provide quantifiable parameters associated with the structure, instantaneous elastic modulus, and density of tissues. By taking advantage of the physical relationship between the axial velocity of ultrasound compression waves and the instantaneous modulus of the propagation material, these techniques can also be used to estimate the in vivo loading environment of relatively superficial soft connective tissues during sports and activities of daily living. This paper highlights key findings from clinical studies in which quantitative transmission-mode ultrasound has been used to measure the properties and loading of bone, tendon, and muscle tissue during common physical activities in healthy and pathological populations.

The Protocol of Ultrasonic Backscatter Measurements of Musculoskeletal Properties

This study aims to introduce the protocol for ultrasonic backscatter measurements of musculoskeletal properties based on a novel ultrasonic backscatter bone diagnostic (UBBD) instrument. Dual-energy X-ray absorptiometry (DXA) can be adopted to measure bone mineral density (BMD) in the hip, spine, legs and the whole body. The muscle and fat mass in the legs and the whole body can be also calculated by DXA body composition analysis. Based on the proposed protocol for backscatter measurements by UBBD, ultrasonic backscatter signals can be measured in vivo, deriving three backscatter parameters [apparent integral backscatter (AIB), backscatter signal peak amplitude (BSPA) and the corresponding arrival time (BSPT)]. AIB may provide important diagnostic information about bone properties. BSPA and BSPT may be important indicators of muscle and fat properties. The standardized backscatter measurement protocol of the UBBD instrument may have the potential to evaluate musculoskeletal characteristics, providing help for promoting the application of the backscatter technique in the clinical diagnosis of musculoskeletal disorders (MSDs), such as osteoporosis and muscular atrophy.

Effect of transducer position on ultrasonic backscatter measurements of cancellous bone

Ultrasonic backscatter techniques are being developed to detect changes in bone caused by osteoporosis and other diseases. Backscatter measurements performed at peripheral skeletal sites such as the heel may place the interrogated region of bone tissue in the acoustic near field of the transducer. The purpose of this study is to investigate how measurements in the near field affect backscatter parameters used for ultrasonic bone assessment. Ultrasonic measurements were performed in a water tank using a planar 2.25 MHz transducer. Signals were acquired for five transducer-specimen distances: N/4, N/2, 3 N/4, N, and 5 N/4, where N is the near-field distance, a location that represents the transition from the near field to far field. Five backscatter parameters previously identified as potentially useful for ultrasonic bone assessment purposes were measured: apparent integrated backscatter, frequency slope of apparent backscatter (FSAB), frequency intercept of apparent backscatter, normalized mean of the backscatter difference, and backscatter amplitude decay constant. All five parameters depended on transducer-specimen distance to varying degrees with FSAB exhibiting the greatest dependence on distance. These results suggest that laboratory studies of bone should evaluate the performance of backscatter parameters using transducer-specimen distances that may be encountered clinically including distances where the ultrasonically interrogated region is in the near field of the transducer.

Ultrasonic Through-Transmission Measurements of Human Musculoskeletal and Fat Properties

The study described here was aimed at investigating the feasibility of using the ultrasonic through-transmission technique to estimate human musculoskeletal and fat properties. Five hundred eighty-two volunteers were assessed by dual-energy X-ray absorptiometry (DXA) and ultrasonic transmission techniques. Bone mineral density (BMD), muscle and fat mass were measured for both legs and the whole body. Hip BMD and spine BMD were also measured. Ultrasonic transmission measurements were performed on the heel, and the measured parameters were broadband ultrasound attenuation (BUA), speed of sound (SOS), ultrasonic stiffness index (SI), T-score and Z-score, which were significantly correlated with all measured BMDs. The optimal correlation was observed between SI and left-leg BMD (p < 0.001) before and after adjustment for age, sex and body mass index (BMI). The linear and partial correlation analyses revealed that BUA and SOS were closely associated with muscle and fat mass, respectively. Multiple regressions revealed that muscle and fat mass significantly contributed to the prediction of transmission parameters, explaining up to 17.83% (p < 0.001) variance independently of BMD. The results suggest that the ultrasonic through-transmission technique could help in the clinical diagnosis of skeletal and muscular system diseases.

In Vivo Comparison of Backscatter Techniques for Ultrasonic Bone Assessment at the Femoral Neck

Ultrasonic techniques are being developed to detect changes in cancellous bone caused by osteoporosis. The goal of this study was to test the relative in vivo performance of eight backscatter parameters developed over the last several years for ultrasonic bone assessment: apparent integrated backscatter (AIB), frequency slope of apparent backscatter (FSAB), frequency intercept of apparent backscatter (FIAB), normalized mean of the backscatter difference (nMBD), normalized slope of the backscatter difference (nSBD), normalized intercept of the backscatter difference (nIBD), normalized backscatter amplitude ratio (nBAR) and backscatter amplitude decay constant (BADC). Backscatter measurements were performed on the left and right femoral necks of 80 adult volunteers (age = 25 ± 11 y) using an imaging system equipped with a convex array transducer. For comparison, additional ultrasonic measurements were performed at the left and right heel using a commercially available heel-bone ultrasonometer that measured the stiffness index. Six of the eight backscatter parameters (all but nSBD and nIBD) exhibited similar and highly significant (p < 0.000001) left-right correlations (0.51 ≤ R ≤ 0.68), indicating sensitivity to naturally occurring variations in bone tissue. Left-right correlations for the stiffness index measured at the heel (R = 0.75) were not significantly better than those produced by AIB, FSAB and FIAB. The short-term precisions of AIB, nMBD, nBAR and BADC (7.8%-11.7%) were comparable to that of the stiffness index measured with the heel-bone ultrasonometer (7.5%).

Ultrasonic Backscatter Measurements of Human Cortical and Trabecular Bone Densities in a Head-Down Bed-Rest Study

This study aims to investigate the feasibility of quantitative ultrasonic backscatter in evaluating human cortical and trabecular bone densities in vivo based on a head-down-tilt bed rest study, with 36 participants tested through 90 d of bed rest and 180 d of recovery. Backscatter measurements were performed using an ultrasonic backscatter bone diagnostic instrument. Backscatter parameters were calculated with a dynamic signal-of-interest method, which was proposed to ensure the same ultrasonic interrogated volume in cortical and trabecular bones. The backscatter parameters exhibited significant correlations with site-matched bone densities provided by high-resolution peripheral quantitative computed tomography (0.33 < |R| < 0.72, p < 0.05). Some bone densities and backscatter parameters exhibited significant changes after the 90-d bed rest. The proposed method can be used to characterize bone densities, and the portable ultrasonic backscatter bone diagnostic device might be used to non-invasively reveal mean bone loss (across a group of people) after long-term bed rest and microgravity conditions of spaceflight missions.

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.

Improving broadband ultrasound attenuation assessment in cancellous bone by mitigating the influence of cortical bone: Phantom and in-vitro study

The purpose of this work was to present a new approach that allows the influence of cortical bone on noninvasive measurement of broadband ultrasound attenuation (BUA) to be corrected. The method, implemented here at 1 MHz makes use of backscattered signal and once refined and clinically confirmed, it would offer an alternative to ionizing radiation based methods, such as DEXA (Dual-energy X-ray absorptiometry), quantitative computed tomography (QCT), radiographic absorptiometry (RA) or single X-ray absorptiometry (SXA), which are clinically approved for assessment of progress of osteoporosis. In addition, as the method employs reflected waves, it might substantially enhance the applicability of BUA - from being suitable to peripheral bones only it would extend this applicability to include such embedded bones as hip and femoral neck. The proposed approach allows the cortical layer parameters used for correction and the corrected value and parameter of the cancellous bone (BUA) to be determined simultaneously from the single (pulse-echo) bone backscattered wave; to the best of the authors' knowledge such approach was not previously reported. The validity of the method was tested using acoustic data obtained from a custom-designed bone-mimicking phantom and a calf femur. The relative error of the attenuation coefficient assessment was determined to be 3.9% and 4.7% for the bone phantom and calf bone specimens, respectively. When the cortical shell influence was not taken into account the corresponding errors were considerably higher 8.3% (artificial bone) and 9.2% (calf femur). As indicated above, once clinically proven, the use of this BUA measurement technique in reflection mode would augment diagnostic power of the attending physician by permitting to include bones, which are not accessible for transmission mode evaluation, e.g. hip, spine, humerus and femoral neck.

Relationships among ultrasonic and mechanical properties of cancellous bone in human calcaneus in vitro

Clinical bone sonometers applied at the calcaneus measure broadband ultrasound attenuation and speed of sound. However, the relation of ultrasound measurements to bone strength is not well-characterized. Addressing this issue, we assessed the extent to which ultrasonic measurements convey in vitro mechanical properties in 25 human calcaneal cancellous bone specimens (approximately 2×4×2cm). Normalized broadband ultrasound attenuation, speed of sound, and broadband ultrasound backscatter were measured with 500kHz transducers. To assess mechanical properties, non-linear finite element analysis, based on micro-computed tomography images (34-micron cubic voxel), was used to estimate apparent elastic modulus, overall specimen stiffness, and apparent yield stress, with models typically having approximately 25-30 million elements. We found that ultrasound parameters were correlated with mechanical properties with R=0.70-0.82 (p<0.001). Multiple regression analysis indicated that ultrasound measurements provide additional information regarding mechanical properties beyond that provided by bone quantity alone (p≤0.05). Adding ultrasound variables to linear regression models based on bone quantity improved adjusted squared correlation coefficients from 0.65 to 0.77 (stiffness), 0.76 to 0.81 (apparent modulus), and 0.67 to 0.73 (yield stress). These results indicate that ultrasound can provide complementary (to bone quantity) information regarding mechanical behavior of cancellous bone.

Ultrasound Speed of Sound Measurements in Trabecular Bone Using the Echographic Response of a Metallic Pin

Bone quality is an important parameter in spine surgery, but its clinical assessment remains difficult. The aim of the work described here was to demonstrate in vitro the feasibility of employing quantitative ultrasound to retrieve bone mechanical properties using an echographic technique taking advantage of the presence of a metallic pin inserted in bone tissue. A metallic pin was inserted in bone tissue perpendicular to the transducer axis. The echographic response of the bone sample was determined, and the echo of the pin inserted in bone tissue and water were compared to determine speed of sound, which was compared with bone volume fraction. A 2-D finite-element model was developed to assess the effect of positioning errors. There was a significant correlation between speed of sound and bone volume fraction (R(2) = 0.6). The numerical results indicate the relative robustness of the measurement method, which could be useful to estimate bone quality intra-operatively.

Correlations of linear and nonlinear ultrasound parameters with density and microarchitectural parameters in trabecular bone

In the present study, correlations of linear and nonlinear ultrasound parameters (speed of sound, normalized broadband ultrasound attenuation, and nonlinear parameter B/A) with bone mineral density and microarchitectural parameters were investigated in 28 bovine femoral trabecular bone samples in vitro. All three ultrasound parameters exhibited relatively high correlation coefficients with the indexes of bone quantity (bone mineral density and bone volume fraction) and lower correlation coefficients with the remaining microarchitectural parameters. These results suggest that B/A, in addition to speed of sound and attenuation, may have potential as an index for the assessment of bone status and osteoporosis.

Influence of cortical endplate on speed of sound in bovine femoral trabecular bone in vitro

Speed of sound (SOS) was measured in 14 bovine femoral trabecular bone samples with and without the cortical endplates with various thicknesses of 1.00, 1.31, 1.47, 1.75, and 2.00 mm. The presence of the cortical endplates resulted in an increase in the mean SOS of 16 m/s (+0.9%) to 91 m/s (+5.3%). The mean SOS measured in the samples with and without the cortical endplates exhibited similar significant correlations with apparent bone density (r = 0.86-0.91). All the SOS measurements were also found to be highly correlated with each other (r = 0.89-0.99).

Correlations between ultrasonic guided wave velocities and bone properties in bovine tibia in vitro

Correlations between ultrasonic guided wave velocities and bone properties were investigated in bovine tibia in vitro. The velocities of the first arriving signal and the slow guided wave, termed V(FAS) and V(SGW), along the long axis of the tibia were measured at 200 kHz in 20 bovine tibiae using the axial transmission technique. V(FAS) yielded significant negative correlation coefficients of -0.54 to -0.66 with the bone properties. In contrast, V(SGW) yielded strong positive correlation coefficients of 0.68-0.84. The best univariate predictor of V(FAS) and V(SGW) was the cortical thickness yielding adjusted squared correlation coefficients of 0.41 and 0.69, respectively.

Multi-site bone ultrasound measurements in elderly women with and without previous hip fractures

About 75% of patients suffering from osteoporosis are not diagnosed. This study describes a multi-site bone ultrasound method for osteoporosis diagnostics. In comparison with axial dual energy X-ray absorptiometry (DXA), the ultrasound method showed good diagnostic performance and could discriminate fracture subjects among elderly females.

INTRODUCTION

Axial DXA, the gold standard diagnostic method for osteoporosis, predicts fractures only moderately. At present, no reliable diagnostic methods are available at the primary health care level. Here, a multi-site ultrasound method is proposed for osteoporosis diagnostics.

METHODS

Thirty elderly women were examined using the ultrasound backscatter measurements in proximal femur, proximal radius, proximal and distal tibia in vivo. First, we predicted the areal bone mineral density (BMD) at femoral neck by ultrasound measurements in tibia combined with specific subject characteristics (density index, DI) and, second, we tested the ability of ultrasound backscatter measurements at proximal femur to discriminate between individuals with previously fractured hips from those without fractures. Areal BMD was determined by axial DXA.

RESULTS

Combined ultrasound parameters, cortical thickness at distal and proximal tibia, with age and weight of the subject, provided a significant estimate of BMD(neck) (r = 0.86, p < 0.001, n = 30). When inserted into FRAX (World Health Organization fracture risk assessment tool), the DI indicated the same treatment proposal as the BMD(neck) with 86% sensitivity and 100% specificity. The receiver operating characteristic analyses, with a combination of ultrasound parameters and patient characteristics, discriminated fracture subjects from the controls similarly as the model combining BMD(neck) and patient characteristics.

CONCLUSIONS

For the first time, ultrasound backscatter measurements of proximal femur were conducted in vivo. The results indicate that ultrasound parameters, combined with patient characteristics, may provide a means for osteoporosis diagnostics.

Exvivo ultrasound attenuation coefficient for human cervical and uterine tissue from 5 to 10 MHz

Attenuation estimation and imaging in the cervix has been utilized to evaluate the onset of cervical ripening during pregnancy. This feature has also been utilized for the acoustic characterization of leiomyomas and myometrial tissue. In this paper, we present direct narrowband substitution measurement values of the variation in the ultrasonic attenuation coefficient in ex vivo human uterine and cervical tissue, in the 5-10 MHz frequency range. At 5 MHz, the attenuation coefficient values are similar for the different orientations of uterine tissue with values of 4.1-4.2 dB/cm, 5.1 dB/cm for the leiomyoma, and 6.3 dB/cm for the cervix. As the frequency increases, the attenuation coefficient values increase and are also spread out, with a value of approximately 12.6 dB/cm for the uterus (both parallel and perpendicular), 16.0 for the leiomyoma, and 26.8 dB/cm for the cervix at 10 MHz. The attenuation coefficient measured increases monotonically over the frequency range measured following a power law.

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.

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.

Predictions of the modified Biot-Attenborough model for the dependence of phase velocity on porosity in cancellous bone

The modified Biot-Attenborough (MBA) model for acoustic wave propagation in porous media has been found useful to predict wave properties in cancellous bone. The present study is aimed at applying the MBA model to predict the dependence of phase velocity on porosity in cancellous bone. The MBA model predicts a phase velocity that decreases nonlinearly with porosity. The optimum values for input parameters of the MBA model, such as compressional speed c(m) of solid bone and phase velocity parameter s(2), were determined by comparing the predictions with previously published measurements in human calcaneus and bovine cancellous bone. The value of the phase velocity parameter s(2)=1.23 was obtained by curve fitting to the experimental data for 53 human calcaneus samples only, assuming a compressional speed c(m)=2500 m/s of solid bone. The root-mean-square error (RMSE) of the curve fit was 15.3m/s. The optimized value of s(2) for all 75 cancellous bone samples including 22 bovine samples was 1.42 with a value of 55 m/s for the RMSE of the curve fit. The latter fit was obtained by using of a value of c(m)=3200 m/s. Although the MBA model relies on the empirical parameters determined from experimental data, it is expected that the model can be usefully employed as a practical tool in the field of clinical ultrasonic bone assessment.

Phase velocity and normalized broadband ultrasonic attenuation in Polyacetal cuboid bone-mimicking phantoms

The variations of phase velocity and normalized broadband ultrasonic attenuation (nBUA) with porosity were investigated in Polyacetal cuboid bone-mimicking phantoms with circular cylindrical pores running normal to the surface along the three orthogonal axes. The frequency-dependent phase velocity and attenuation coefficient in the phantoms with porosities from 0% to 65.9% were measured from 0.65 to 1.10 MHz. The results showed that the phase velocity at 880 kHz decreased linearly with porosity, whereas the nBUA increased linearly with porosity. This study provides a useful insight into the relationships between ultrasonic properties and porosity in bone at porosities lower than 70%.

Bayesian estimation of the underlying bone properties from mixed fast and slow mode ultrasonic signals

We recently proposed that the observed apparent negative dispersion in bone can arise from the interference between fast wave and slow wave modes, each exhibiting positive dispersion [Marutyan et al., J. Acoust. Soc. Am. 120, EL55-EL61 (2006)]. In the current study, we applied Bayesian probability theory to solve the inverse problem: extracting the underlying properties of bone. Simulated mixed mode signals were analyzed using Bayesian probability. The calculations were implemented using the Markov chain Monte Carlo with simulated annealing to draw samples from the marginal posterior probability for each parameter.

What is ultrasound?

This paper is based on material presented at the start of a Health Protection Agency meeting on ultrasound and infrasound. In answering the question 'what is ultrasound?', it shows that the simple description of a wave which transports mechanical energy through the local vibration of particles at frequencies of 20 kHz or more, with no net transport of the particles themselves, can in every respect be misleading or even incorrect. To explain the complexities responsible for this, the description of ultrasound is first built up from the fundamental properties of these local particle vibrations. This progresses through an exposition of the characteristics of linear waves, in order to explain the propensity for, and properties of, the nonlinear propagation which occurs in many practical ultrasonic fields. Given the Health Protection environment which framed the original presentation, explanation and examples are given of how these complexities affect issues of practical importance. These issues include the measurement and description of fields and exposures, and the ability of ultrasound to affect tissue (through microstreaming, streaming, cavitation, heating, etc.). It is noted that there are two very distinct regimes, in terms of wave characteristics and potential for bioeffect. The first concerns the use of ultrasound in liquids/solids, for measurement or material processing. For biomedical applications (where these two processes are termed diagnosis and therapy, respectively), the issue of hazard has been studied in depth, although this has not been done to such a degree for industrial uses of ultrasound in liquids/solids (sonar, non-destructive testing, ultrasonic processing etc.). However, in the second regime, that of the use of ultrasound in air, although the waves in question tend to be of much lower intensities than those used in liquids/solids, there is a greater mismatch between the extent to which hazard has been studied, and the growth in commercial applications for airborne ultrasound.

Ultrasonic characterization of cancellous bone using apparent integrated backscatter

Apparent integrated backscatter (AIB) is a measure of the frequency-averaged (integrated) backscattered power contained in some portion of a backscattered ultrasonic signal. AIB has been used extensively to study soft tissues, but its usefulness as a tissue characterization technique for cancellous bone has not been demonstrated. To address this, we performed measurements on 17 specimens of cancellous bone over two different frequency ranges using a 1 MHz and 5 MHz broadband ultrasonic transducer. Specimens were obtained from bovine tibiae and prepared in the shape of cubes (15 mm side length) with faces oriented along transverse (anterior, posterior, medial and lateral) and longitudinal (superior and inferior) principal anatomic directions. A mechanical scanning system was used to acquire multiple backscatter signals from each direction for each cube. AIB demonstrated highly significant linear correlations with bone mineral density (BMD) for both the transverse (R2 = 0.817) and longitudinal (R2 = 0.488) directions using the 5 MHz transducer. In contrast, the correlations with density were much weaker for the 1 MHz transducer (R2 = 0.007 transverse, R2 = 0.228 longitudinal). In all cases where a significant correlation was observed, AIB was found to decrease with increasing BMD.

Kramers-Kronig analysis of attenuation and dispersion in trabecular bone

A restricted-bandwidth form of the Kramers-Kronig dispersion relations is applied to in vitro measurements of ultrasonic attenuation and dispersion properties of trabecular bone specimens from bovine tibia. The Kramers-Kronig analysis utilizes only experimentally measured properties and avoids extrapolation of ultrasonic properties beyond the known bandwidth. Compensation for the portions of the Kramers-Kronig integrals over the unknown bandwidth is partially achieved by the method of subtractions, where a subtraction frequency acts as an adjustable parameter. Good agreement is found between experimentally measured and Kramers-Kronig reconstructed dispersions. The restricted-bandwidth approach improves upon other forms of the Kramers-Kronig relations and may provide further insight into how ultrasound interacts with trabecular bone.

Singular spectrum analysis applied to backscattered ultrasound signals from in vitro human cancellous bone specimens

Mean scatterer spacing (MSS) holds particular promise for the detection of changes in quasiperiodic tissue microstructures such as may occur during development of disease in the liver, spleen, or bones. Many techniques that may be applied for MSS estimation (temporal and spectral autocorrelation, power spectrum and cepstrum, higher order statistics, and quadratic transformation) characterize signals that contain a mixture of periodic and nonperiodic contributions. In contrast, singular spectrum analysis (SSA), a method usually applied in nonlinear dynamics, first identifies components of signals corresponding to periodic structures and, second, identifies dominant periodicity. Thus, SSA may better separate periodic structures from nonperiodic structures and noise. Using an ultrasound echo simulation model, we previously demonstrated SSA's potential to identify MSS of structures in quasiperiodic scattering media. The current work aims to observe the behavior of MSS estimation by SSA using ultrasound measurements in phantom materials (two parallel, nylon-line phantoms and four foam phantoms of different densities). The SSA was able to estimate not only the nylon-line distances but also nylon-line thickness. The method also was sensitive to the average pore-size differences of the four sponges. The algorithms then were applied to characterize human cancellous bone microarchitectures. Using 1-MHz center-frequency, radio-frequency ultrasound signals, MSS was measured in 24 in vitro bone samples and ranged from 1.0 to 1.7 mm. The SSA MSS estimates correlate significantly to MSS measured independently from synchrotron microtomography, r2 = 0.68. Thus, application of SSA to backscattered ultrasound signals seems to be useful for providing information linked to tissue microarchitecture that is not evident from clinical images.

Bone quality: where do we go from here?

As was true 10 years ago, tremendous interest surrounds the concept of "bone quality," as shown by the intense and growing research activity in the field. The urgency to advance knowledge in this area is motivated by the need to understand not only the causes of increased skeletal fragility with aging and disease, but also the mechanisms by which drugs reduce fracture risk. As reflected in the preceding articles, in the past decade collaborations between biologists, physicists, engineers and clinicians have led to new insights regarding the biological, material and structural features that contribute to skeletal fragility. Despite these new insights, important issues remain unresolved. Our challenges lie in identifying, describing and understanding the totality of features and characteristics that determine a bone's ability to resist fracture, and using this information to identify new therapeutic targets and develop better biomarkers and noninvasive imaging modalities. This summary is intended to integrate reports in the literature with presentations and discussions that occurred during the meeting, to highlight areas of consensus and those of continued controversy, and thus to identify critical areas for future research.

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).

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.

Temperature dependence of ultrasonic propagation speed and attenuation in canine tissue

Previously reported data on the temperature dependence of propagation speed in tissues generally span only temperature ranges up to 60 degrees C. However, with the emerging use of thermal ablative therapies, information on variation in this parameter over higher temperature ranges is needed. Measurements of the ultrasonic propagation speed and attenuation in tissue in vitro at discrete temperatures ranging from 25 to 95 degrees C was performed for canine liver, muscle, kidney and prostate using 3 and 5 MHz center frequencies. The objective was to produce information for calibrating temperature-monitoring algorithms during ablative therapy. Resulting curves of the propagation speed vs. temperature for these tissues can be divided into three regions. In the 25-40 degrees C range, the speed of sound increase rapidly with temperature. It increases moderately with temperature in the 40-70 degrees C range, and it then decreases with increasing temperature from 70-95 degrees C. Attenuation coefficient behavior with temperature is different for the various tissues. For liver, the attenuation coefficient is nearly constant with temperature. For kidney, attenuation increases approximately linearly with temperature, while for muscle and prostate tissue, curves of attenuation vs. temperature are flat in the 25-50 degrees C range, slowly rise at medium temperatures (50-70 degrees C), and level off at higher temperatures (70-90 degrees C). Measurements were also conducted on a distilled degassed water sample and the results closely follow values from the literature.

Ultrasonic wave speed measurement using the time-delay profile of rf-backscattered signals: simulation and experimental results

Conventional methods determine the ultrasonic wave speed measuring the medium path length propagated by a pulsed wave and the corresponding time-of-flight. In this work, the wave speed is determined without the need of the path length. A transmit transducer sends a pulsed wave into the medium (wave speed constant along the beam axis) and the backscattered signal is collected by a hydrophone placed at two distinct positions near the transmitted beam. The time-delay profile, between gated windows of the two rf-signals received by the hydrophone, is determined using a cross-correlation method. Also, a theoretical time-delay profile is determined considering the wave speed as a parameter. The estimated wave speed is obtained upon minimization of the rms error between theoretical and experimental time-delay profiles. A PZT conically focused transmitting transducer with center frequency of 3.3 MHz, focal depth of 30 mm, and beam full width (-3 dB) of 2 mm at the focus was used together with a PZT hydrophone (0.8 mm of aperture). The method was applied to three phantoms (wave speed of 1220, 1540, and 1720 m/s) and, in vitro, to fresh bovine liver sample, immersed in a temperature-controlled water bath. The results present a relative speed error less than 3% when compared with the sound speed obtained by a conventional method.