Guided ultrasonic waves in long bones: modelling, experiment and in vivo application

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.

Ultrasonic Guided Waves in Bone: A Decade of Advancement in Review

The use of guided wave ultrasonography as a means to assess cortical bone quality has been a significant practice in bone quantitative ultrasound for more than 20 years. In this article, the key developments within the technology of ultrasonic guided waves (UGW) in long bones during the past decade are documented. The covered topics include data acquisition configurations available for measuring bone guided waveforms, signal processing techniques applied to bone UGW, numerical modeling of ultrasonic wave propagation in cortical long bones, formulation of inverse approaches to extract bone properties from observed ultrasonic signals, and clinical studies to establish the technology's application and efficacy. The review concludes by highlighting specific challenging problems and future research directions. In general, the primary purpose of this work is to provide a comprehensive overview of bone guided-wave ultrasound, especially for newcomers to this scientific field.

Study on ultrasonic wave propagation in equine leg bone for screening bucked shin

For simple, safe, portable, and inexpensive evaluation suitable for leg bone diseases of racehorses in the field, an ultrasonic measurement technique was applied to evaluate wave velocities. A digital model of the third metacarpal bone with the bucked shin was fabricated using high-resolution peripheral quantitative computerized tomography data of a racehorse. This model was anisotropic and heterogeneous, and was constructed using the measured ultrasonic wave velocities in the bone. With this model, ultrasonic wave propagation along the bone axis was simulated using the elastic finite-difference time-domain method. We found two main waves with different propagation velocities. The fast-waves showed a wave velocity close to the longitudinal wave in the axial direction. However, the apparent velocities changed dramatically owing to bone surface irregularities (changes of the shape) in the area of bucked shin. The slow-waves showed a wave velocity close to the shear wave, which was unaffected by the bone surface irregularities. The simple comparison of different wave behaviors may be a suitable parameter for the initial in vivo screening of bucked shin in the legs of racehorses, which can be performed in the field.

Recycled windshield glass as new material for producing ultrasonic phantoms of cortical bone-healing stages

The quantitative ultrasound technique was used to evaluate bone-mimicking phantoms; however, these phantoms do not mimic the intermediate stages of cortical bone healing. We propose using windshield glass as an original material to produce phantoms that mimic the characteristics of three different stages of cortical bone healing. This material was processed via a route that included breaking, grinding, compacting, drying, and sintering in four temperature groups: 625 °C, 645 °C, 657 °C, and 663 °C. The parameters evaluated were the ultrasonic longitudinal phase velocity (cL), corrected (αc) ultrasonic attenuation coefficient, and bulk density (ρs). The results showed that the mean values ofcL,αc,andρsvaried from 2, 398 to 4, 406 m·s^-1, 3 to 10 dB·cm^-1, and 1, 563 to 2, 089 kg·m^-3, respectively. The phantoms exhibited properties comparable with the three stages of cortical bone healing and can be employed in diagnostic and therapeutic studies using ultrasound.

Identification of long-range ultrasonic guided wave characteristics in cortical bone by modelling

The propagation of ultrasonic guided waves in cortical bone has potential to inform medical caregivers about the condition of the bone structure. However, as waveguides, human long bones such as the tibia are complex in terms of their material behavior and their geometric features. They exhibit anisotropic elasticity and internal damping. For the first time, wave propagation is modelled in the irregular hollow tibial cross-section, which varies along its long axis. Semi-analytical, frequency domain, and time domain finite element analyses providing complimentary information about long-range wave propagation characteristics in such a waveguide are applied to the mid-diaphyseal region of a human tibia. Simulating the guided waves generated by a contact transducer, the signals received in axial transmission indicate the consistent presence of low phase velocity non-dispersive propagating modes. The guided waves capable of traveling long distances have strong potential for diagnosis of fracture healing.

Experimental and Computational Investigation of Guided Waves in a Human Skull

We investigate guided (Lamb) waves in a human cadaver skull through experiments and computational simulations. Ultrasonic wedge transducers and scanning laser Doppler vibrometry are used respectively to excite and measure Lamb waves propagating in the cranial bone of a degassed skull. Measurements are performed over a section of the parietal bone and temporal bone spanning the squamous suture. The experimental data are analyzed for the identification of wave modes and the characterization of dispersion properties. In the parietal bone, for instance, the A₀ wave mode is excited between 200 and 600 kHz, and higher-order Lamb waves are excited from 1 to 1.8 MHz. From the experimental dispersion curves and average thickness extracted from the skull computed tomography scan, we estimate average isotropic material properties that capture the essential dispersion characteristics using a semi-analytical finite-element model. We also explore the leaky and non-leaky wave behavior of the degassed skull with water loading in the cranial cavity. Successful excitation of leaky Lamb waves is confirmed (for higher-order wave modes with phase velocity faster than the speed of sound in water) from 500 kHz to 1.5 MHz, which may find applications in imaging and therapeutics at the brain periphery or skull-brain interface (e.g., for metastases). The non-leaky A₀ Lamb wave mode propagates between 200 and 600 kHz, with or without fluid loading, for potential use in skull-related diagnostics and imaging (e.g., for sutures).

Shear-Actuation and Vibrometer Reception of Penetrating Ultrasonic Guided Wave Modes in Human Tibia

The hollow long bones of the human appendicular skeleton are known to support the propagation of ultrasonic guided waves, whose potential for diagnosing bone health is being investigated. In this study, ultrasonic guided waves propagating in the diaphysis of human tibia are characterized experimentally and numerically in the frequency range around 200 kHz. The experiment involves a unique combination of omni-directional shear transducer-based excitation and detection using a 1D laser Doppler vibrometer. The cluster of phase velocities obtained from a linear array of time-history data using space-time Fourier transform is found to be in the non-dispersive low-phase velocity region of the dispersion curves obtained for a tibial cross-section. Time-domain finite element analysis revealed that the displacement components normal to the surface are significant, even though the loading is from a shear transducer. Furthermore, semi-analytical finite element analysis revealed that the wave structures of the wave modes contained within the cluster of low-phase velocity modes are consistent with the displacement profiles obtained from the time-domain analysis. The experimental results show that the low-phase velocity mode cluster has sufficient intensity to propagate axially at least 85 mm in the mid-diaphyseal region.

Ex Vivo Assessment of Cortical Bone Properties Using Low-Frequency Ultrasonic Guided Waves

The early diagnosis of osteoporosis through bone quality assessment is a major public health challenge. Research in axial transmission using ultrasonic guided waves has shown the method to be sensitive to the geometrical and mechanical properties of the cortical layer in long bones. However, because of the asymmetric nature of cortical bone, the introduction of a more elaborate numerical model than the analytical plate and cylinder models, as well as its inversion, continues to be challenging. The aim of this article is, therefore, to implement a bone-like geometry using semianalytical finite-element (SAFE) modeling to perform the inverse characterization of ex vivo radii at low frequencies (< 60 kHz). Five cadaveric radiuses were taken from donors aged between 53 and 88 and tested using a typical axial transmission configuration at the middle of the diaphysis. The dispersion curves of the propagating modes were measured experimentally and then systematically compared with the solutions obtained with the SAFE method. For each sample, four parameters were estimated using a parameter identification procedure: 1) the bulk density; 2) the thickness; 3) the outer diameter; and 4) a shape factor (SF). The results showed a moderate agreement between the predicted bulk density and the average voxel value calculated from X-ray computed tomography images. Furthermore, a good agreement was observed between the geometrical parameters (thickness, outer diameter, and SF) and the reference images.

Artificial neural network to estimate micro-architectural properties of cortical bone using ultrasonic attenuation: A 2-D numerical study

The goal of this study is to estimate micro-architectural parameters of cortical porosity such as pore diameter (φ), pore density (ρ) and porosity (ν) of cortical bone from ultrasound frequency dependent attenuation using an artificial neural network (ANN). First, heterogeneous structures with controlled pore diameters and pore densities (mono-disperse) were generated, to mimic simplified structure of cortical bone. Then, more realistic structures were obtained from high resolution CT scans of human cortical bone. 2-D finite-difference time-domain simulations were conducted to calculate the frequency-dependent attenuation in the 1-8 MHz range. An ANN was then trained with the ultrasonic attenuation at different frequencies as the input feature vectors while the output was set as the micro-architectural parameters (pore diameter, pore density and porosity). The ANN is composed of three fully connected dense layers with 24, 12 and 6 neurons, connected to the output layer. The dataset was trained over 6000 epochs with a batch size of 16. The trained ANN exhibits the ability to predict the micro-architectural parameters with high accuracy and low losses. ANN approaches could potentially be used as a tool to help inform physics-based modelling of ultrasound propagation in complex media such as cortical bone. This will lead to the solution of inverse-problems to retrieve bone micro-architectural parameters from ultrasound measurements for the non-invasive diagnosis and monitoring osteoporosis.

The Ability of Ultrasonic Backscatter Parametric Imaging to Characterize Bovine Trabecular Bone

The ultrasonic backscatter technique holds the promise of characterizing bone density and microstructure. This paper conducts ultrasonic backscatter parametric imaging based on measurements of apparent integrated backscatter (AIB), spectral centroid shift (SCS), frequency slope of apparent backscatter (FSAB), and frequency intercept of apparent backscatter (FIAB) for representing trabecular bone mass and microstructure. We scanned 33 bovine trabecular bone samples using a 7.5 MHz focused transducer in a 20 mm × 20 mm region of interest (ROI) with a step interval of 0.05 mm. Images based on the ultrasonic backscatter parameters (i.e., AIB, SCS, FSAB, and FIAB) were constructed to compare with photographic images of the specimens as well as two-dimensional (2D) μ-CT images from approximately the same depth and location of the specimen. Similar structures and trabecular alignments can be observed among these images. Statistical analyses demonstrated that the means and standard deviations of the ultrasonic backscatter parameters exhibited significant correlations with bone density (|R| = 0.45-0.78, p < 0.01) and bone microstructure (|R| = 0.44-0.87, p < 0.001). Some bovine trabecular bone microstructure parameters were independently associated with the ultrasonic backscatter parameters (ΔR² = 4.18%-44.45%, p < 0.05) after adjustment for bone apparent density (BAD). The results show that ultrasonic backscatter parametric imaging can provide a direct view of the trabecular microstructure and can reflect information about the density and microstructure of trabecular bone.

Ultrasound-Based Estimates of Cortical Bone Thickness and Porosity Are Associated With Nontraumatic Fractures in Postmenopausal Women: A Pilot Study

Recent ultrasound (US) axial transmission techniques exploit the multimode waveguide response of long bones to yield estimates of cortical bone structure characteristics. This pilot cross-sectional study aimed to evaluate the performance at the one-third distal radius of a bidirectional axial transmission technique (BDAT) to discriminate between fractured and nonfractured postmenopausal women. Cortical thickness (Ct.Th) and porosity (Ct.Po) estimates were obtained for 201 postmenopausal women: 109 were nonfractured (62.6 ± 7.8 years), 92 with one or more nontraumatic fractures (68.8 ± 9.2 years), 17 with hip fractures (66.1 ± 10.3 years), 32 with vertebral fractures (72.4 ± 7.9 years), and 17 with wrist fractures (67.8 ± 9.6 years). The areal bone mineral density (aBMD) was obtained using DXA at the femur and spine. Femoral aBMD correlated weakly, but significantly with Ct.Th (R = 0.23, p < 0.001) and Ct.Po (R = -0.15, p < 0.05). Femoral aBMD and both US parameters were significantly different between the subgroup of all nontraumatic fractures combined and the control group (p < 0.05). The main findings were that (1) Ct.Po was discriminant for all nontraumatic fractures combined (OR = 1.39; area under the receiver operating characteristic curve [AUC] equal to 0.71), for vertebral (OR = 1.96; AUC = 0.84) and wrist fractures (OR = 1.80; AUC = 0.71), whereas Ct.Th was discriminant for hip fractures only (OR = 2.01; AUC = 0.72); there was a significant association (2) between increased Ct.Po and vertebral and wrist fractures when these fractures were not associated with any measured aBMD variables; (3) between increased Ct.Po and all nontraumatic fractures combined independently of aBMD neck; and (4) between decreased Ct.Th and hip fractures independently of aBMD femur. BDAT variables showed comparable performance to that of aBMD neck with all types of fractures (OR = 1.48; AUC = 0.72) and that of aBMD femur with hip fractures (OR = 2.21; AUC = 0.70). If these results are confirmed in prospective studies, cortical BDAT measurements may be considered useful for assessing fracture risk in postmenopausal women. © 2019 American Society for Bone and Mineral Research.

Influence of guided waves in bone on pulse-inversion contrast-enhanced ultrasound

PURPOSE

Guided waves generated from bone cortex inevitably act on microbubbles flowing through skeletal muscle capillaries in contrast-enhanced ultrasound (CEUS) and might influence the image quality. However, the action mechanism underlying the guided waves influence is still unknown, especially under contrast pulse-inversion transmission mode. This study aimed to clarify the influence of guided waves on pulse-inversion CEUS, which was investigated via in vitro infusion experiments.

METHOD

Tibia guided waves were detected at pulse-inversion transmission and then characterized by using a short-time Fourier transform energy distribution. Using results at normal incidence as a baseline, the influence of guided wave dispersion on the contrast and resolution of pulse-inversion CEUS was investigated at an oblique incidence through continuous microbubbles infusion experiments in a vessel-tibia flow phantom.

RESULTS

Frequency-dispersive property of tibia guided waves was observed at phases 0° and 180°, which improved the contrast of CEUS and reduced its resolution. Pulse-inversion CEUS balanced the contrast enhancement and resolution degeneration induced by guided waves. By contrast, contrast-to-tissue ratio of pulse-inversion CEUS increased by up to 109.1 ± 13.2% (P < 0.05) due to guided waves and its resolution was up to 0.9 ± 0.1 times that of baseline.

CONCLUSIONS

Alterations of contrast and resolution in pulse-inversion CEUS induced by guided waves might provide an additional assessment for the capillary perfusion in the skeletal muscle near the bone cortex.

Can low-frequency guided waves at the tibia paired with machine learning differentiate between healthy and osteopenic/osteoporotic subjects? A pilot study

PURPOSE

Axial transmission quantitative acoustics (ax-QA) has shown to be a promising tool for assessing bone health and properties in a safe, inexpensive, and portable manner. This study investigated the efficacy of low-frequency ax-QA measured at the tibia, paired with a support vector machine (SVM) approach for combining multiple acoustic indicators, to diagnose osteoporosis as defined by bone mineral density.

METHODS

This pilot study measured 41 female subjects using ax-QA (flexural mode, 3 kHz) at the tibia and using dual X-ray absorptiometry (DXA) at the lumbar spine, femoral neck, and distal radius. For each location, a threshold classifier and SVM were trained to differentiate between healthy and non-healthy subjects based on the phase velocity at different frequencies. Receiver Operating Characteristics and area under curve values (AUC) were used to assess the classifiers' performances for various thresholds and class-weights.

RESULTS

The SVM outperformed the threshold classifier for all three bone locations at low false positive rates. While differentiation between healthy and non-healthy bone states was poor for the spine (AUC: 0.56 ± 0.04), good to moderate performances were observed for the radius (AUC: 0.83 ± 0.03) and hip (AUC: 0.71 ± 0.04).

CONCLUSIONS

Low-frequency ax-QA has demonstrated potential for complementing DXA in screening for osteoporosis at the radius and hip. Through further addition of acoustic indicators ax-QA could provide a diagnostic alternative in third-world countries, and bring bone health screening and monitoring into the hands of clinicians and general health practitioners everywhere.

A Quantitative Approach for the Bone-implant Osseointegration Assessment Based on Ultrasonic Elastic Guided Waves

Quantitative and reliable monitoring of osseointegration will help further evaluate the integrity of the orthopaedic construct to promote novel prosthesis design and allow early mobilisation. Quantitative assessment of the degree or the lack of osseointegration is important for the clinical management with the introduction of prosthetic implants to amputees. Acousto-ultrasonic wave propagation has been used in structural health monitoring as well as human health monitoring but so far has not extended to osseointegrated implants or prostheses. This paper presents an ultrasonic guided wave approach to assess the osseointegration of a novel implant. This study explores the potential of integrating structural health monitoring concepts into a new osseointegrated implant. The aim is to demonstrate the extension of acousto-ultrasonic techniques, which have been widely reported for the structural health monitoring of engineering structures, to assess the state of osseointegration of a bone and implant. To illustrate this potential, this paper will report on the experimental findings which investigated the unification of an aluminium implant and bone-like geometry surrogate. The core of the test specimen is filled with silicone and wrapped with plasticine to simulate the highly damped cancellous bone and soft tissue, respectively. To simulate the osseointegration process, a 2-h adhesive epoxy is used to bond the surrogate implant and a bone-like structure. A series of piezoelectric elements are bonded onto the surrogate implant to serve as actuators and sensors. The actuating piezoelectric element on an extramedullary strut is excited with a 1 MHz pulse signal. The reception of the ultrasonic wave by the sensing elements located on the adjacent and furthest struts is used to assess the integration of this implant to the parent bone structure. The study shows an Osseointegration Index can be formulated by using engineering and acousto-ultrasonic methods to measure the unification of a bone and implant. This also highlights a potential quantitative evaluation technique regardless of bone-implant geometry and soft tissue damping.

Bone cortical thickness and porosity assessment using ultrasound guided waves: An ex vivo validation study

Several studies showed the ability of the cortex of long bones such as the radius and tibia to guide mechanical waves. Such experimental evidence has given rise to the emergence of a category of quantitative ultrasound techniques, referred to as the axial transmission, specifically developed to measure the propagation of ultrasound guided waves in the cortical shell along the axis of long bones. An ultrasound axial transmission technique, with an automated approach to quantify cortical thickness and porosity is described. The guided modes propagating in the cortex are recorded with a 1-MHz custom made linear transducer array. Measurement of the dispersion curves is achieved using a two-dimensional spatio-temporal Fourier transform combined with singular value decomposition. Automatic parameters identification is obtained through the solution of an inverse problem in which the dispersion curves are predicted with a two-dimensional transverse isotropic free plate model. Thirty-one radii and fifteen tibiae harvested from human cadavers underwent axial transmission measurements. Estimates of cortical thickness and porosity were obtained on 40 samples out of 46. The reproducibility, given by the root mean square error of the standard deviation of estimates, was 0.11 mm for thickness and 1.9% for porosity. To assess accuracy, site-matched micro-computed tomography images of the bone specimens imaged at 9 μm voxel size served as the gold standard. Agreement between micro-computed tomography and axial transmission for quantification of thickness and porosity at the radius and tibia ranged from R²=0.63 for porosity (root mean square error RMSE=1.8%) to 0.89 for thickness (RMSE=0.3 mm). Despite an overall good agreement for porosity, the method performs less well for porosities lower than 10%. The heterogeneity and general complexity of cortical bone structure, which are not fully accounted for by our model, are suspected to weaken the model approximation. This study presents the first validation study for assessing cortical thickness and porosity using the axial transmission technique. The automatic signal processing minimizes operator-dependent errors for parameters determination. Recovering the waveguide characteristics, that is to say cortical thickness and porosity, could provide reliable information about skeletal status and future fracture risk.

Fatigue evaluation of long cortical bone using ultrasonic guided waves

Bone fatigue fracture is a progressive disease due to stress concentration. This study aims to evaluate the long bone fatigue damage using the ultrasonic guided waves. Two-dimensional finite-difference time-domain method was employed to simulate the ultrasonic guided wave propagation in the long bone under different elastic modulus. The experiment was conducted on a 3.8 mm-thick bovine bone plate. The phase velocities of two fundamental guided modes, A1 and S1, were measured by using the axial transmission technique. Simulation shows that the phase velocities of guided modes A1 and S1 decrease with the increasing of the fatigue damage. After 20,000 cycles of fatigue loading on the bone plate, the average phase velocities of A1 and S1 modes were 6.6% and 5.3% respectively, lower than those of the intact bone. The study suggests that ultrasonic guided waves can be potentially used to evaluate the fatigue damage in long bones.

Dispersion characteristics of the flexural wave assessed using low frequency (50-150kHz) point-contact transducers: A feasibility study on bone-mimicking phantoms

Guided waves-based techniques are currently under development for quantitative cortical bone assessment. However, the signal interpretation is challenging due to multiple mode overlapping. To overcome this limitation, dry point-contact transducers have been used at low frequencies for a selective excitation of the zeroth order anti-symmetric Lamb A0 mode, a mode whose dispersion characteristics can be used to infer the thickness of the waveguide. In this paper, our purpose was to extend the technique by combining a dry point-contact transducers approach to the SVD-enhanced 2-D Fourier transform in order to measure the dispersion characteristics of the flexural mode. The robustness of our approach is assessed on bone-mimicking phantoms covered or not with soft tissue-mimicking layer. Experiments were also performed on a bovine bone. Dispersion characteristics of measured modes were extracted using a SVD-based signal processing technique. The thickness was obtained by fitting a free plate model to experimental data. The results show that, in all studied cases, the estimated thickness values are in good agreement with the actual thickness values. From the results, we speculate that in vivo cortical thickness assessment by measuring the flexural wave using point-contact transducers is feasible. However, this assumption has to be confirmed by further in vivo studies.

Towards assessing cortical bone porosity using low-frequency quantitative acoustics: A phantom-based study

PURPOSE

Cortical porosity is a key characteristic governing the structural properties and mechanical behaviour of bone, and its quantification is therefore critical for understanding and monitoring the development of various bone pathologies such as osteoporosis. Axial transmission quantitative acoustics has shown to be a promising technique for assessing bone health in a fast, non-invasive, and radiation-free manner. One major hurdle in bringing this approach to clinical application is the entanglement of the effects of individual characteristics (e.g. geometry, porosity, anisotropy etc.) on the measured wave propagation. In order to address this entanglement problem, we therefore propose a systematic bottom-up approach, in which only one bone property is varied, before addressing interaction effects. This work therefore investigated the sensitivity of low-frequency quantitative acoustics to changes in porosity as well as individual pore characteristics using specifically designed cortical bone phantoms.

MATERIALS AND METHODS

14 bone phantoms were designed with varying pore size, axial-, and radial pore number, resulting in porosities (bone volume fraction) between 0% and 15%, similar to porosity values found in human cortical bone. All phantoms were manufactured using laser sintering, measured using axial-transmission acoustics and analysed using a full-wave approach. Experimental results were compared to theoretical predictions based on a modified Timoshenko theory.

RESULTS

A clear dependence of phase velocity on frequency and porosity produced by increasing pore size or radial pore number was demonstrated, with the velocity decreasing by between 2-5 m/s per percent of additional porosity, which corresponds to -0.5% to -1.0% of wave speed. While the change in phase velocity due to axial pore number was consistent with the results due to pore size and radial pore number, the relative uncertainties for the estimates were too high to draw any conclusions for this parameter.

CONCLUSIONS

This work has shown the capability of low-frequency quantitative acoustics to reflect changes in porosity and individual pore characteristics and demonstrated that additive manufacturing is an appropriate method that allows the influence of individual bone properties on the wave propagation to be systematically assessed. The results of this work opens perspectives for the efficient development of a multi-frequency, multi-mode approach to screen, diagnose, and monitor bone pathologies in individuals.

Semi-analytical finite-element modeling approach for guided wave assessment of mechanical degradation in bones

Numerical models based on the Semi Analytical Finite-Element method are used to study the characteristics of guided wave modes supported by bone-like multi-layered tubular structures. The method is first validated using previous literature and experimental studies on phantoms mimicking healthy and osteoporotic conditions of cortical bone, and later used to study a trilayer marrow–bone–tissue system at varying mechanical degradation levels. The results show that bone condition strongly affects the modal properties of axially propagating guided waves and indicates that L(0,3) and F(1,6) are suitable modes for assessing the mechanical condition of the bone. The work here reports suitable modal selection and their dispersion properties which would the aid in development of a transduction mechanism for mechanical assessment of bones.

Relationships of the group velocity of the time-reversed Lamb wave with bone properties in cortical bone in vitro

The present study aims to investigate the feasibility of using the time-reversed Lamb wave as a new method for noninvasive characterization of long cortical bones. The group velocity of the time-reversed Lamb wave launched by using the modified time reversal method was measured in 15 bovine tibiae, and their correlations with the bone properties of the tibia were examined. The group velocity of the time-reversed Lamb wave showed significant positive correlations with the bone properties (r=0.55-0.81). The best univariate predictor of the group velocity of the time-reversed Lamb wave was the cortical thickness, yielding an adjusted squared correlation coefficient (r²) of 0.64. These results imply that the group velocity of the time-reversed Lamb wave, in addition to the velocities of the first arriving signal and the slow guided wave, could potentially be used as a discriminator for osteoporosis.

Predicting bone strength with ultrasonic guided waves

Recent bone quantitative ultrasound approaches exploit the multimode waveguide response of long bones for assessing properties such as cortical thickness and stiffness. Clinical applications remain, however, challenging, as the impact of soft tissue on guided waves characteristics is not fully understood yet. In particular, it must be clarified whether soft tissue must be incorporated in waveguide models needed to infer reliable cortical bone properties. We hypothesize that an inverse procedure using a free plate model can be applied to retrieve the thickness and stiffness of cortical bone from experimental data. This approach is first validated on a series of laboratory-controlled measurements performed on assemblies of bone- and soft tissue mimicking phantoms and then on in vivo measurements. The accuracy of the estimates is evaluated by comparison with reference values. To further support our hypothesis, these estimates are subsequently inserted into a bilayer model to test its accuracy. Our results show that the free plate model allows retrieving reliable waveguide properties, despite the presence of soft tissue. They also suggest that the more sophisticated bilayer model, although it is more precise to predict experimental data in the forward problem, could turn out to be hardly manageable for solving the inverse problem.

Simulation study of axial ultrasonic wave propagation in heterogeneous bovine cortical bone

The effect of the heterogeneity of the long cortical bone is an important factor when applying the axial transmission technique. In this study, the axial longitudinal wave velocity distributions in specimens from the mid-shaft of a bovine femur were measured, in the MHz range. Bilinear interpolation and the piecewise cubic Hermite interpolating polynomial method were used to construct three-dimensional (3D) axial velocity models with a resolution of 40 μm. By assuming the uniaxial anisotropy of the bone and using the results of previous experimental studies [Yamato, Matsukawa, Yanagitani, Yamazaki, Mizukawa, and Nagano (2008b). Calcified Tissue Int. 82, 162-169; Nakatsuji, Yamamoto, Suga, Yanagitani, Matsukawa, Yamazaki, and Matsuyama (2011). Jpn. J. Appl. Phys. 50, 07HF18], the distributions of all elastic moduli were estimated to obtain a 3D heterogeneous bone model and a uniform model. In the heterogeneous model, moduli at the surface were smaller than those inside the model. The elastic finite-difference time-domain method was used to simulate axial ultrasonic wave propagation in these models. In the heterogeneous model, the wavefront of the first arriving signal (FAS) was dependent on the heterogeneity, and the FAS velocity depended on the measured position. These phenomena were not observed in the uniform model.

Multichannel processing for dispersion curves extraction of ultrasonic axial-transmission signals: Comparisons and case studies

Some pioneering studies have shown the clinical feasibility of long bones evaluation using ultrasonic guided waves. Such a strategy is typically designed to determine the dispersion information of the guided modes to infer the elastic and structural characteristics of cortical bone. However, there are still some challenges to extract multimode dispersion curves due to many practical limitations, e.g., high spectral density of modes, limited spectral resolution and poor signal-to-noise ratio. Recently, two representative signal processing methods have been proposed to improve the dispersion curves extraction. The first method is based on singular value decomposition (SVD) with advantages of multi-emitter and multi-receiver configuration for enhanced mode extraction; the second one uses linear Radon transform (LRT) with high-resolution imaging of the dispersion curves. To clarify the pros and cons, a face to face comparison was performed between the two methods. The results suggest that the LRT method is suitable to separate the guided modes at low frequency-thickness-product ( fh) range; for multimode signals in broadband fh range, the SVD-based method shows more robust performances for weak mode enhancement and noise filtering. Different methods are valuable to cover the entire fh range for processing ultrasonic axial transmission signals measured in long cortical bones.

Improving accuracy through density correction in guided wave tomography

The accurate quantification of wall loss caused by corrosion is critical to the reliable life estimation of pipes and pressure vessels. Traditional thickness gauging by scanning a probe is slow and requires access to all points on the surface; this is impractical in many cases as corrosion often occurs where access is restricted, such as beneath supports where water collects. Guided wave tomography presents a solution to this; by transmitting guided waves through the region of interest and exploiting their dispersive nature, it is possible to build up a map of thickness. While the best results have been seen when using the fundamental modes A0 and S0 at low frequency, the complex scattering of the waves causes errors within the reconstruction. It is demonstrated that these lead to an underestimate in wall loss for A0 but an overestimate for S0. Further analysis showed that this error was related to density variation, which was proportional to thickness. It was demonstrated how this could be corrected for in the reconstructions, in many cases resulting in the near-elimination of the error across a range of defects, and greatly improving the accuracy of life estimates from guided wave tomography.

Reliability of Phase Velocity Measurements of Flexural Acoustic Waves in the Human Tibia In-Vivo

Purpose

Axial-transmission acoustics have shown to be a promising technique to measure individual bone properties and detect bone pathologies. With the ultimate goal being the in-vivo application of such systems, quantification of the key aspects governing the reliability is crucial to bring this method towards clinical use.

Materials and Methods

This work presents a systematic reliability study quantifying the sources of variability and their magnitudes of in-vivo measurements using axial-transmission acoustics. 42 healthy subjects were measured by an experienced operator twice per week, over a four-month period, resulting in over 150000 wave measurements. In a complementary study to assess the influence of different operators performing the measurements, 10 novice operators were trained, and each measured 5 subjects on a single occasion, using the same measurement protocol as in the first part of the study.

Results

The estimated standard error for the measurement protocol used to collect the study data was ∼ 17 m/s (∼ 4% of the grand mean) and the index of dependability, as a measure of reliability, was Φ = 0.81. It was shown that the method is suitable for multi-operator use and that the reliability can be improved efficiently by additional measurements with device repositioning, while additional measurements without repositioning cannot improve the reliability substantially. Phase velocity values were found to be significantly higher in males than in females (p < 10⁻⁵) and an intra-class correlation coefficient of r = 0.70 was found between the legs of each subject.

Conclusions

The high reliability of this non-invasive approach and its intrinsic sensitivity to mechanical properties opens perspectives for the rapid and inexpensive clinical assessment of bone pathologies, as well as for monitoring programmes without any radiation exposure for the patient.

Computational Study of the Effect of Cortical Porosity on Ultrasound Wave Propagation in Healthy and Osteoporotic Long Bones

Computational studies on the evaluation of bone status in cases of pathologies have gained significant interest in recent years. This work presents a parametric and systematic numerical study on ultrasound propagation in cortical bone models to investigate the effect of changes in cortical porosity and the occurrence of large basic multicellular units, simply called non-refilled resorption lacunae (RL), on the velocity of the first arriving signal (FAS). Two-dimensional geometries of cortical bone are established for various microstructural models mimicking normal and pathological tissue states. Emphasis is given on the detection of RL formation which may provoke the thinning of the cortical cortex and the increase of porosity at a later stage of the disease. The central excitation frequencies 0.5 and 1 MHz are examined. The proposed configuration consists of one point source and multiple successive receivers in order to calculate the FAS velocity in small propagation paths (local velocity) and derive a variation profile along the cortical surface. It was shown that: (a) the local FAS velocity can capture porosity changes including the occurrence of RL with different number, size and depth of formation; and (b) the excitation frequency 0.5 MHz is more sensitive for the assessment of cortical microstructure.

Influence of Guided Waves in Tibia on Non-linear Scattering of Contrast Agents

The aim of this study was to elucidate the linear and non-linear responses of ultrasound contrast agent (UCA) to frequency-dispersive guided waves from the tibia cortex, particularly two individual modes, S0 (1.23 MHz) and A1 (2.06 MHz). The UCA responses to guided waves were illustrated through the Marmottant model derived from measured guided waves, and then verified by continuous infusion experiments in a vessel-tibia flow phantom. These UCA responses were further evaluated by the enhanced ratio of peak values and the resolutions of UCA backscattered echoes. Because of the individual modes S0 and A1 in the tibia, the peak values of the UCA backscattered echoes were enhanced by 83.57 ± 7.35% (p < 0.05) and 80.77 ± 6.60% (p < 0.01) in the UCA subharmonic frequency and subharmonic imaging, respectively. However, corresponding resolutions were 0.78 ± 0.07 (p < 0.05) and 0.72 ± 0.12 (p < 0.01) times those without guided wave disturbances, respectively. Even though the resolution was partly degenerated, the subharmonic detection sensitivity of UCA was improved by the guided waves. Thus, UCA responses to the double-frequency guided waves should be further explored to benefit the detection of capillary perfusion in tissue layers near the bone cortex, particularly for perfusion imaging in the free flaps and skeletal muscles.

Ultrasound propagation in cortical bone: Axial transmission and backscattering simulations

Cortical bone is a heterogeneous, composite medium with a porosity from 5-10%. The characterization of cortical bone using ultrasonic techniques is a complicated procedure especially in numerical studies as several assumptions must be made to describe the concentration and size of pores. This study presents numerical simulations of ultrasound propagation in two-dimensional numerical models of cortical bone to investigate the effect of porosity on: a) the propagation of the first arriving signal (FAS) velocity using the axial transmission method, and b) the displacement and scattering amplitude in the backward direction. The excitation frequency 1 MHz was used and different receiving positions were examined to provide a variation profile of the examined parameters along cortical bone. Cortical porosity was simulated using ellipsoid scatterers and the concentrations of 0-10% were examined. The results indicate that the backscattering method is more appropriate for the evaluation of cortical porosity in comparison to the axial transmission method.

Propagation of time-reversed Lamb waves in bovine cortical bone in vitro

The present study aims to investigate the propagation of time-reversed Lamb waves in bovine cortical bone in vitro. The time-reversed Lamb waves were successfully launched at 200 kHz in 18 bovine tibiae through a time reversal process of Lamb waves. The group velocities of the time-reversed Lamb waves in the bovine tibiae were measured using the axial transmission technique. They showed a significant correlation with the cortical thickness and tended to follow the theoretical group velocity of the lowest order antisymmetrical Lamb wave fairly well, consistent with the behavior of the slow guided wave in long cortical bones.

Combined estimation of thickness and velocities using ultrasound guided waves: a pioneering study on in vitro cortical bone samples

This paper reports for the first time on inverse estimation of several bone properties from guided-wave measurements in human bone samples. Previously, related approaches have focused on ultrasonic estimation of a single bone property at a time. The method is based on two steps: the multi-Lamb mode response is analyzed using the singular value decomposition signal processing method recently introduced in the field, then an identification procedure is run to find thickness and anisotropic elastic properties of the considered specimen. Prior to the measurements on bone, the method is validated on cortical bone-mimicking phantoms. The repeatability and the trueness of the estimated parameters on bone-mimicking phantoms were found around a few percent. Estimation of cortical thickness on bone samples was in good agreement with cortical thickness derived from high-resolution peripheral quantitative computed tomography data analysis of the samples.

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.

Multi-frequency axial transmission bone ultrasonometer

The last decade has seen a surge in the development of axial transmission QUS (Quantitative UltraSound) technologies for the assessment of long bones using various modes of acoustic waves. The condition of cortical bones and the development of osteoporosis are determined by numerous mechanical, micro-structural, and geometrical or macro-structural bone properties like hardness, porosity and cortical thickness. Such complex manifestations of osteoporosis require the evaluation of multiple parameters with different sensitivities to the various properties of bone that are affected by the disease. This objective may be achieved by using a multi-frequency ultrasonic examination The ratio of the acoustic wavelength to the cortical thickness can be changed by varying the frequency of the ultrasonic pulse propagating through the long bone that results in the change in composition of the induced wave comprised of a set of numerous modes of guided, longitudinal, and surface acoustic waves. The multi-frequency axial transmission QUS method developed at Artann Laboratories (Trenton, NJ) is implemented in the Bone Ultrasonic Scanner (BUSS). In the current version of the BUSS, a train of ultrasonic pulses with 60, 100, 400, 800, and 1200 kHz frequencies is used. The developed technology was tested on a variety of bone phantoms simulating normal, osteopenic, and osteoporotic bones. The results of this study confirm the feasibility of the multi-frequency approach for the assessment of the processes leading to osteoporosis.

Excitation of ultrasonic Lamb waves using a phased array system with two array probes: phantom and in vitro bone studies

Long bones are good waveguides to support the propagation of ultrasonic guided waves. The low-order guided waves have been consistently observed in quantitative ultrasound bone studies. Selective excitation of these low-order guided modes requires oblique incidence of the ultrasound beam using a transducer-wedge system. It is generally assumed that an angle of incidence, θi, generates a specific phase velocity of interest, co, via Snell's law, θi=sin(-1)(vw/co) where vw is the velocity of the coupling medium. In this study, we investigated the excitation of guided waves within a 6.3-mm thick brass plate and a 6.5-mm thick bovine bone plate using an ultrasound phased array system with two 0.75-mm-pitch array probes. Arranging five elements as a group, the first group of a 16-element probe was used as a transmitter and a 64-element probe was a receiver array. The beam was steered for six angles (0°, 20°, 30°, 40°, 50°, and 60°) with a 1.6-MHz source signal. An adjoint Radon transform algorithm mapped the time-offset matrix into the frequency-phase velocity dispersion panels. The imaged Lamb plate modes were identified by the theoretical dispersion curves. The results show that the 0° excitation generated many modes with no modal discrimination and the oblique beam excited a spectrum of phase velocities spread asymmetrically about co. The width of the excitation region decreased as the steering angle increased, rendering modal selectivity at large angles. The phenomena were well predicted by the excitation function of the source influence theory. The low-order modes were better imaged at steering angle ⩾30° for both plates. The study has also demonstrated the feasibility of using the two-probe phased array system for future in vivo study.

Quantification of guided mode propagation in fractured long bones

Guided modes propagation in intact, fractured and healing long bone has drawn significant research interests. However, mode quantifications for the direct comparison are still necessary to address. The aim of the study is to analyze the mode interaction with a notch-fracture in the long bone and find quantitative ultrasound parameters sensitive to depth and width variation of the fracture. We analyzed the impacts of the partially and completely diaphyseal osteotomy on fundamental guided modes propagation using the two-dimension finite-difference time-domain (2D-FDTD) simulations. The long bones were built as three layer models by a cortical plate embedded between overlying soft tissue and inner-coated marrow. Narrowband low-frequency sinusoids (100 kHz) were employed to only excite two fundamental guided modes. The mode amplitude variations were investigated as functions of the gap-breakage width and depth. It is found that the transverse fractures have strong influences on the anti-symmetric mode A0 transmission and reflection, whereas amplitudes of the symmetric mode S0 are not sensitive to the fracture degree. The quantitative results consistently indicate that reflection energy and transmission coefficients of the S0 and A0 modes can be used to quantify the mode interaction in the fractured long bone and further to evaluate long bone fracture status. Future study is needed to investigate the physical experiments on realistic fractured long bone and to insure that the proposed ultrasound parameters can be used to quantitatively evaluate the long bone fracture in clinical application.

Association between low-frequency ultrasound and hip fractures -- comparison with DXA-based BMD

BACKGROUND

New methods for diagnosing osteoporosis and evaluating fracture risk are being developed. We aim to study the association between low-frequency (LF) axial transmission ultrasound and hip fracture risk in a population-based sample of older women.

METHODS

The study population consisted of 490 community-dwelling women (78-82 years). Ultrasound velocity (V(LF)) at mid-tibia was measured in 2006 using a low-frequency scanning axial transmission device. Bone mineral density (BMD) at proximal femur measured using dual-energy x-ray absorptiometry (DXA) was used as the reference method. The fracture history of the participants was collected from December 1997 until the end of 2010. Lifestyle-related risk factors and mobility were assessed at 1997.

RESULTS

During the total follow-up period (1997-2010), 130 women had one or more fractures, and 20 of them had a hip fracture. Low V(LF) (the lowest quartile) was associated with increased hip fracture risk when compared with V(LF) in the normal range (Odds ratio, OR = 3.3, 95% confidence interval (CI) 1.3-8.4). However, V(LF) was not related to fracture risk when all bone sites were considered. Osteoporotic femoral neck BMD was associated with higher risk of a hip fracture (OR = 4.1, 95% CI 1.6-10.5) and higher risk of any fracture (OR = 2.4, 95% CI 1.6-3.8) compared to the non-osteoporotic femoral neck BMD. Decreased VLF remained a significant risk factor for hip fracture when combined with lifestyle-related risk factors (OR = 3.3, 95% CI 1.2-9.0).

CONCLUSION

Low V(LF) was associated with hip fracture risk in older women even when combined with lifestyle-related risk factors. Further development of the method is needed to improve the measurement precision and to confirm the results.

A study on Rayleigh wave dispersion in bone according to Mindlin's Form II gradient elasticity

The classical elasticity cannot effectively describe bone's mechanical behavior since only homogeneous media and local stresses are assumed. Additionally, it cannot predict the dispersive nature of the Rayleigh wave which has been reported in experimental studies and was also demonstrated in a previous computational study by adopting Mindlin's Form II gradient elasticity. In this work Mindlin's theory is employed to analytically determine the dispersion of Rayleigh waves in a strain gradient elastic half-space. An isotropic semi-infinite space is considered with properties equal to those of bone and dynamic behavior suffering from microstructural effects. Microstructural effects are considered by incorporating four intrinsic parameters in the stress analysis. The results are presented in the form of group and phase velocity dispersion curves and compared with existing computational results and semi-analytical curves calculated for a simpler case of Rayleigh waves in dipolar gradient elastic half-spaces. Comparisons are also performed with the velocity of the first-order antisymmetric mode propagating in a dipolar plate so as to observe the Rayleigh asymptotic behavior. It is shown that Mindlin's Form II gradient elasticity can effectively describe the dispersive nature of Rayleigh waves. This study could be regarded as a step toward the ultrasonic characterization of bone.

Axial transmission method for long bone fracture evaluation by ultrasonic guided waves: simulation, phantom and in vitro experiments

Mode conversion occurs when the ultrasonic guided waves encounter fractures. The aim of this study was to investigate the feasibility of fracture assessment in long cortical bone using guided-mode conversion. Mode conversion behavior between the fundamental modes S0 and A0 was analyzed. The expressions proposed for modal velocity were used to identify the original and converted modes. Simulations and phantom experiments were performed using 1.0-mm-thick steel plates with a notch width of 0.5 mm and notch depths of 0.2, 0.4, 0.6 and 0.8 mm. Furthermore, in vitro experiments were carried out on nine ovine tibias with 1.0-mm-wide partial transverse gap break and cortical thickness varying from 2.10 to 3.88 mm. The study confirmed that mode conversion gradually becomes observable as fracture depth increases. Energy percentages of the converted modes correlated strongly with fracture depth, as illustrated by the frequency-sweeping experiments on steel phantoms (100-1100 kHz, r(2) = 0.97, p < 0.0069) and the fixed-frequency experiments on nine ovine tibias (250 kHz, r(2) = 0.97, p < 0.0056). The approaches described, including mode excitation, velocity expressions and energy percentage criteria, may also contribute to ultrasonic monitoring of long bone fracture healing.

Photo-acoustic excitation and optical detection of fundamental flexural guided wave in coated bone phantoms

Photo-acoustic (PA) imaging was combined with skeletal quantitative ultrasound (QUS) for assessment of human long bones. This approach permitted low-frequency excitation and detection of ultrasound so as to efficiently receive the thickness-sensitive fundamental flexural guided wave (FFGW) through a coating of soft tissue. The method was tested on seven axisymmetric bone phantoms, whose 1- to 5-mm wall thickness and 16-mm diameter mimicked those of the human radius. Phantoms were made of a composite material and coated with a 2.5- to 7.5-mm layer of soft material that mimicked soft tissue. Ultrasound was excited with a pulsed Nd:YAG laser at 1064-nm wavelength and received on the same side of the coated phantom with a heterodyne interferometer. The FFGW was detected at 30-kHz frequency. Fitting the FFGW phase velocity by the FLC(1,1) tube mode provided an accurate (9.5 ± 4.0%) wall thickness estimate. Ultrasonic in vivo characterization of cortical bone thickness may thus become possible.

Influence of porosity, pore size, and cortical thickness on the propagation of ultrasonic waves guided through the femoral neck cortex: a simulation study

The femoral neck is a common fracture site in elderly people. The cortical shell is thought to be the major contributor to the mechanical competence of the femoral neck, but its microstructural parameters are not sufficiently accessible under in vivo conditions with current X-ray-based methods. To systematically investigate the influences of pore size, porosity, and thickness of the femoral neck cortex on the propagation of ultrasound, we developed 96 different bone models (combining 6 different pore sizes with 4 different porosities and 4 different thicknesses) and simulated the ultrasound propagation using a finite-difference time-domain algorithm. The simulated single-element emitter and receiver array consisting of 16 elements (8 inferior and 8 superior) were placed at anterior and posterior sides of the bone, respectively (transverse transmission). From each simulation, we analyzed the waveform collected by each of the inferior receiver elements for the one with the shortest time of flight. The first arriving signal of this waveform, which is associated with the wave traveling through the cortical shell, was then evaluated for its three different waveform characteristics (TOF: time point of the first point of inflection of the received signal, Δt: difference between the time point at which the signal first crosses the zero baseline and TOF, and A: amplitude of the first extreme of the first arriving signal). From the analyses of these waveform characteristics, we were able to develop multivariate models to predict pore size, porosity, and cortical thickness, corresponding to the 96 different bone models, with remaining errors in the range of 50 μm for pore size, 1.5% for porosity, and 0.17 mm for cortical thickness.

Analysis of ultrasonic waves propagating in a bone plate over a water half-space with and without overlying soft tissue

Recent in vitro studies have shown that guided waves can characterize bone properties. However, for clinical applications to be viable, the soft-tissue layer should be considered. This study examined the effect of soft tissue on guided waves using a bovine bone plate over a water half-space and overlaid by a 4-mm gelatin-based soft-tissue mimic. The data (with and without soft tissue) clearly show a high-frequency, fast-propagating wave packet and a low-frequency, delayed phase group. The presence of soft tissue attenuates the signals significantly and increases mode density and number as predicted by theory. The data retain higher frequency content than the bone-plate data at large offsets. Using theoretical dispersion curves, the guided modes can be identified with mode 1 (similar to the A0 Lamb mode) minimally affected by the addition of soft tissue. There is infiltration of high-frequency, late-arriving energy within the low-velocity guided-wave regime. Results of travel-time calculation suggest that P-wave and PP-reflections/multiples within the soft tissue may be responsible for the high-frequency oscillations.

Assessment of the fundamental flexural guided wave in cortical bone by an ultrasonic axial-transmission array transducer

The fundamental flexural guided wave (FFGW), as modeled, for example, by the A0 Lamb mode, is a clinically useful indicator of cortical bone thickness. In the work described in this article, we tested so-called multiridge-based analysis, based on the crazy climber algorithm and short-time Fourier transform, for assessment of the FFGW component recorded by a clinical array transducer featuring a limited number of elements. Methods included numerical finite-element simulations and experiments in bone phantoms and human radius specimens (n = 41). The proposed approach enabled extraction of the FFGW component and determination of its group velocity. This group velocity was in good agreement with theoretical predictions and possessed reasonable sensitivity to cortical width (r(2) = 0.51, p < 0.001) in the in vitro experiments. It is expected that the proposed approach enables related clinical application. Further work is still needed to analyze in more detail the challenges related to the impact of the overlying soft tissue.

Discrimination of fractures by low-frequency axial transmission ultrasound in postmenopausal females

SUMMARY

In this cross-sectional study, 95 postmenopausal women, with and without fracture history, were measured by low-frequency axial transmission ultrasound. The measured ultrasound velocity discriminated the fractured subjects from the nonfractured ones equally or better than peripheral quantitative computed tomography (pQCT) and dual energy x-ray absorptiometry (DXA). These results suggest that low-frequency ultrasound is suitable for bone fragility assessment.

INTRODUCTION

Quantitative low-frequency axial transmission ultrasound is a promising modality for assessing mineral density and geometrical properties of long bones such as radius and tibia. The aim of the current study was to evaluate the ability of low-frequency axial transmission ultrasound to discriminate fractures retrospectively in postmenopausal women.

METHODS

A cross-sectional study involved 95 female subjects aged 45-88 years, whose fracture information was gathered retrospectively. The fracture group was defined as subjects with one or more low-/moderate-energy fractures. The radius and tibial shaft were measured with a custom-made ultrasonometer to assess the velocity of the low-frequency first-arriving signal (V (LF)). Site-matched pQCT was used to measure volumetric cortical and subcortical bone mineral density (sBMD), and cortical thickness (CTh). Areal BMD (aBMD) was measured using DXA for the whole body (WB), lumbar spine, and hip.

RESULTS

The majority (19/32; 59 %) of the fractures were in the upper limb. V (LF) in the radius, but not in the tibia, discriminated fractures with an age- and BMI-adjusted odds ratio (OR) of 2.06 (95 % CI 1.21-3.50, p < 0.01). In the radius, CTh and cortical BMD (CBMD) significantly discriminated fractures, as did the total, cortical, and sBMD in the tibia (adjusted OR 1.35-2.15, p < 0.05). Sensitivity and specificity were similar among all the measurements (area under the receiver operating characteristic curve 0.74-0.81, p < 0.001).

CONCLUSIONS

Low-frequency axial transmission ultrasound in the radius was able to discriminate fractured subjects from the nonfractured ones. This suggests that low-frequency axial transmission ultrasound has the potential to assess bone fragility in postmenopausal women.

Study of factors affecting the magnitude and nature of ultrasound exposure with in vitro set-ups

Therapeutic ultrasound is a clinically applied method to improve fracture healing and holds great potential as a manipulator of biologic material relevant to tissue engineering approaches. Unfortunately, the cell stimulating property of ultrasound is not known, which inhibits the optimal use of this technique. Additionally, many in vitro studies in this field use ultrasound configurations that are vulnerable to errors during calibration and use. These errors arise from the structural simplicity and incomplete characterization of these configurations. In this study, pulse-echo ultrasound, laser Doppler vibrometry and Schlieren imaging were applied to noninvasively characterize common in vitro experimental configurations. Fine wire thermocouple measurements were conducted to characterize any possible temperature rise during the ultrasound exposures. The results quantified the frequency dependent sound transmission through culture wells and the standing wave effect within the cell volume. These effects can cause uncertainty of up to 700% in the actual ultrasound exposure experienced by the cell. A temperature rise of 2.7°C was measured from an ultrasound configuration commonly used in vitro ultrasound studies. Furthermore, wave mode conversion in culture wells was observed, emphasizing the complexity of these sonications. Similar type Lamb waves have been observed in bone in vivo. Thus, Lamb waves may be a mechanism for stimulating the cells.

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.

Comparative examination of human proximal tibiae in vitro by ultrasonic guided waves and pQCT

The velocity of ultrasonic guided waves in long bones is dependent upon two determinants of bone strength: the cortical thickness and the material properties. In this study, six human proximal tibiae in vitro were examined to test the efficacy of an ultrasonic method based on guided waves. Peripheral quantitative computed tomography (pQCT) was used as the comparative reference modality. The guided wave velocity (c(F)) was derived from two-dimensional (2-D) spatial-temporal waveform profiles formed by multiple ultrasonic signals acquired along the bones at 100 kHz frequency and passed wavelet processing. The ultrasonic profiles from the examined bones were ranged according to pQCT measurements of cortical thickness (CTh), and cortical bone mineral density (CBMD). Strong correlations between c(F) and CTh (r(S) = 0.83, p < 0.0001) and CBMD (r(S) = 0.88, p < 0.0001) resulted. The study confirmed a potential for guided waves to assess atrophic changes of the proximal tibia.

Analysis of superimposed ultrasonic guided waves in long bones by the joint approximate diagonalization of eigen-matrices algorithm

The parameters of ultrasonic guided waves (GWs) are very sensitive to mechanical and structural changes in long cortical bones. However, it is a challenge to obtain the group velocity and other parameters of GWs because of the presence of mixed multiple modes. This paper proposes a blind identification algorithm using the joint approximate diagonalization of eigen-matrices (JADE) and applies it to the separation of superimposed GWs in long bones. For the simulation case, the velocity of the single mode was calculated after separation. A strong agreement was obtained between the estimated velocity and the theoretical expectation. For the experiments in bovine long bones, by using the calculated velocity and a theoretical model, the cortical thickness (CTh) was obtained. For comparison with the JADE approach, an adaptive Gaussian chirplet time-frequency (ACGTF) method was also used to estimate the CTh. The results showed that the mean error of the CTh acquired by the JADE approach was 4.3%, which was smaller than that of the ACGTF method (13.6%). This suggested that the JADE algorithm may be used to separate the superimposed GWs and that the JADE algorithm could potentially be used to evaluate long bones.