Comb Detection for Measuring Shear Wave Propagation.

Plane wave compounding (PWC) is widely used to measure the propagation of shear waves. Implementing PWC on most commercial ultrasound scanners is challenging because all channel (>128) data must be processed or transferred to the host computing unit in real time. Comb detection transmits multiple focused beams simultaneously and results in a reduced number of receive lines to be processed in parallel. These comb beams are scanned laterally to acquire receive lines at different lateral positions in order to obtain data over a large region of interest (ROI). One of the potential issues with using multiple simultaneously transmitted beams is the issue of crosstalk between the beams. Crosstalk is analyzed through simulated beam patterns, simulated B-mode images, and motion data from shear wave elastography (SWE) experiments. Using a Hamming window on transmit and receive can suppress crosstalk to 1.2% root-mean-square error (RMSE, normalized RMSE to the peak magnitude of the reference signal) for shear wave motion signals. Four comb beams with three laterally scanned locations cover almost the entire field of view (FOV) and achieve the same frame rate as PWC with three angles. Phantom and in vivo studies demonstrate comparable motion data of comb detection to PWC in terms of motion signal quality and measured phase velocity. In addition, comb detection provides motion with lower noise and stronger signals than PWC, which is believed to be due to the advantages of transmitting focused beams rather than plane waves (PWs).

Non-invasive Measurement of the Viscoelasticity of the Optic Nerve and Sclera for Assessing Papilledema: A Pilot Clinical Study.

OBJECTIVE: The purpose of this study was to evaluate our novel ultrasound vibro-elastography (UVE) technique for assessing patients with papilledema by non-invasively measuring shear wave speed (SWS), elasticity and viscosity properties of the optic nerve and sclera. METHODS: Shear wave speeds were measured at three frequencies-100, 150 and 200 Hz-on the optic nerve and sclera tissues for assessing patients with papilledema resulting from idiopathic intracranial hypertension (IIH). The method was evaluated in six papilledema patients and six controls on two separate locations for each participant (i.e., optic nerve and posterior sclera). SWSs of the optic nerve and sclera were analyzed by using a 2-D speed map technique within a circular region of interest (ROI) (i.e., the diameter of the ROI was 1.5 mm × 3.0 mm at the optic nerve and sclera, respectively). Elasticity and viscosity were then analyzed using the wave speed dispersion over the three frequencies. RESULTS: We measured values of SWS at both locations, optic nerve and sclera, of the right eye and left eye at three different frequencies in IIH patients and controls. The SWS (mean ± standard deviation [m/s]) of the right eye was significantly higher at the sclera in IIH patients compared with controls (i.e., patients vs. controls: 5.91 ± 0.54 vs. 3.86 ± 0.56, p < 0.0001 at 100 Hz), but there was no significant difference at the optic nerve (i.e., patients vs. controls: 3.62 ± 0.39 vs. 3.36 ± 0.35, p = 0.1100 at 100Hz). We observed increased elasticity (kPa) in IIH patients, indicating there are significant differences in elasticity between patients and controls at the optic nerve and sclera (i.e., right eye [patients vs. controls]: 14.42 ± 6.59 vs. 6.5 ± 5.71, p = 0.0065 [optic nerve]; 33.04 ± 10.62 vs. 9.16 ± 7.15, p < 0.0001 [sclera]). Viscosity was also (Pa·s) higher in the sclera and optic nerve of the left eye (i.e., left eye [patient vs. control]: 8.89 ± 4.37 vs. 7.27 ± 5.01, p = 0.3790 (optic nerve); 16.05 ± 10.79 vs. 8.49 ± 6.09, p < 0.0194 [sclera]). CONCLUSION: This research illustrates the feasibility of using our UVE system to evaluate stiffness of different tissues in the eye non-invasively. It suggests that the viscoelasticity of the posterior sclera is higher than that of the optic nerve. We found that the posterior sclera is stiffer than the optic nerve in patients with papilledema resulting from IIH, making UVE a potential non-invasive technique for assessing papilledema.

Factors Associated With Ultrasound Color Doppler Twinkling by Breast Biopsy Markers: In Vitro and Ex Vivo Evaluation of 35 Commercially Available Markers.

BACKGROUND. Targeted axillary lymph node dissection after neoadjuvant systemic therapy (NST) for breast cancer depends on identifying marked metastatic lymph nodes. However, ultrasound visualization of biopsy markers is challenging. OBJECTIVE. The purpose of our study was to identify biopsy markers that show actionable twinkling in cadaveric breast and to assess the association of actionable twinkling with markers' surface roughness. METHODS. Commercial breast biopsy markers were evaluated for twinkling artifact in various experimental conditions relating to scanning medium (solid gel phantom, ultrasound coupling gel, cadaveric breast), transducer (ML6-15, 9L, C1-6), and embedding material (present vs absent). Markers were assigned twinkling scores from 0 (confident in no twinkling) to 4 (confident in exuberant twinkling); a score of 3 or greater represented actionable twinkling (sufficient confidence to rely solely on twinkling for target localization). Markers were hierarchically advanced to evaluation with increasingly complex media if showing at least minimal twinkling for a given medium. A 3D coherence optical profiler measured marker surface roughness. Mixed-effects proportional odds regression models assessed associations between twinkling scores and transducer and embedding material; Wilcoxon rank sum test evaluated associations between actionable twinkling and surface roughness. RESULTS. Thirty-five markers (21 with embedding material) were evaluated. Ten markers without embedding material advanced to evaluation in cadaveric breast. Higher twinkling scores were associated with presence of embedding material (odds ratio [OR] = 5.05 in solid gel phantom, 9.84 in coupling gel) and transducer (using the C1-6 transducer as reference; 9L transducer: OR = 0.36, 0.83, and 0.04 in solid gel phantom, ultrasound coupling gel, and cadaveric breast; ML6-15 transducer: OR = 0.07, 0.18, and 0.00 respectively; post hoc p between 9L and ML6-15: p < .001, p = .02, and p = .04). In cadaveric breast, three markers (Cork, Professional Q, MRI [Flex]) exhibited actionable twinkling for two or more transducers; surface roughness was significantly higher for markers with than without actionable twinkling for C1-6 (median values: 0.97 vs 0.35, p = .02) and 9L (1.75 vs 0.36; p = .002) transducers. CONCLUSION. Certain breast biopsy markers exhibited actionable twinkling in cadaveric breast. Twinkling was observed with greater confidence for the C1-6 and 9L transducers than the ML6-15 transducer. Actionable twinkling was associated with higher marker surface roughness. CLINICAL IMPACT. Use of twinkling for marker detection could impact preoperative or intraoperative localization after NST.

Ultrasonographic Detection and Surgical Retrieval of a Nonmetallic Twinkle Marker in Breast Cancer: Pilot Study.

Purpose To evaluate the short-term safety of a nonmetallic twinkle marker and compare its conspicuity at color Doppler US with that of standard breast biopsy clips and radioactive seeds by using B-mode US in axillary lymph nodes. Materials and Methods This prospective study (November 2020-July 2021) of participants with node-positive breast cancer who completed chemotherapy involved placing a twinkle marker at the time of preoperative radioactive seed localization. A five-point scoring system (1 = easiest, 5 = most difficult) was used to rate the ease of identifying the clip, seed, and twinkle marker on postlocalization sonograms, mammograms, specimen radiographs, and gross pathologic specimens. Descriptive statistics were used. Results Eight women (mean age, 57 years ± 16 [SD]) were enrolled. The median scores for US conspicuity of each device were 3.9 (range, 3.7-5.0) for the radioactive seed, 2.4 (range, 1.0-5.0) for the clip, and 2.0 (range, 1.0-4.3) for the twinkle marker. In six of eight participants, the twinkle marker was the most identifiable at US. The seeds, clips, and twinkle markers were scored "very easy" to identify on seven of eight postlocalization mammograms. The surgeon retrieved all eight twinkle markers 1-3 days after localization. In all 16 interpretations, the seeds, clips, and twinkle markers were rated as very easy to identify on specimen radiographs. The clip was the most difficult device to identify at pathologic examination in all participants, and the twinkle marker was the easiest to identify in seven of eight participants. Conclusion This pilot study demonstrates that the safety and ease of US detection of a twinkling tissue marker may be comparable to a biopsy clip. Keywords: Ultrasonography, US-Doppler, Breast, Localization, Surgery Clinical trial registration no. NCT04674852 © RSNA, 2022.

Measurement of wave propagation through a tube using dual transducers for elastography in arteries.

Objective.Measuring waves induced with acoustic radiation force (ARF) in arteries has been studied over the last decade. To date, it remains a challenge to quantitatively assess the local arterial biomechanical properties. The cylindrical shape and waveguide behavior of waves propagating in the arterial wall pose complexities to determining the mechanical properties of the artery.Approach. In this paper, an artery-mimicking tube in water is examined utilizing three-dimensional measurements. The cross-section of the tube is measured while a transducer is translated over 41 different positions along the length of the tube. Motion in the radial direction is calculated using two components of motion which are measured from the two orthogonal views of the cross-section. This enables more accurate estimation of motion along the circumference of tube.Main results. The results provide more information to categorize the motion in tube wall into two types of responses: a transient response and a steady state response. The transient response is caused by ARF application and the waves travel along the length of the tube for a relatively short period of time. This corresponds to the axial and circumferential propagating waves. The two circumferential waves travel along the circumference of tube in CW (clockwise) and CCW (counter-clockwise) direction and result in a standing wave. By using a directional filter, the two waves were successfully separated, and their propagation was more clearly visualized. As a steady state response, a circumferential mode is generated showing a symmetric motion (i.e. the proximal and distal walls move in the opposite direction) following the transient response.Significance.This study provides a more comprehensive understanding of the waves produced in an artery-mimicking tube with ARF application, which will provide opportunities for improving measurement of arterial mechanical properties.

The influence of acoustic radiation force beam shape and location on wave spectral content for arterial dispersion ultrasound vibrometry.

Objective. Arterial dispersion ultrasound vibrometry (ADUV) relies on the use of guided waves in arterial geometries for shear wave elastography measurements. Both the generation of waves through the use of acoustic radiation force (ARF) and the techniques employed to infer the speed of the resulting wave motion affect the spectral content and accuracy of the measurement. In particular, the effects of the shape and location of the ARF beam in ADUV have not been widely studied. In this work, we investigated how such variations of the ARF beam affect the induced motion and the measurements in the dispersive modes that are excited.Approach.The study includes an experimental evaluation on an arterial phantom and anin vivovalidation of the observed trends, observing the two walls of the waveguide, simultaneously, when subjected to variations in the ARF beam extension (F/N) and focus location.Main results.Relying on the theory of guided waves in cylindrical shells, the shape of the beam controls the selection and nature of the induced modes, while the location affects the measured dispersion curves (i.e. variation of phase velocity with frequency or wavenumber, multiple modes) across the waveguide walls.Significance.This investigation is important to understand the spectral content variations in ADUV measurements and to maximize inversion accuracy by tuning the ARF beam settings in clinical applications.

Twinkling-guided ultrasound detection of polymethyl methacrylate as a potential breast biopsy marker: a comparative investigation.

Since its first description 25 years ago, color Doppler twinkling has been a compelling ultrasound feature in diagnosing urinary stones. While the fundamental cause of twinkling remains elusive, the distinctive twinkling signature is diagnostically valuable in clinical practice. It can be inferred that if an entity twinkles, it empirically has certain physical features. This work investigates a manipulable polymeric material, polymethyl methacrylate (PMMA), which twinkles and has measurable surface roughness and porosity that likely contribute to twinkling. Comparative investigation of these structural properties and of the twinkling signatures of breast biopsy markers made from PMMA and selected commercially available markers showed how twinkling can improve ultrasound detection of devices intentionally designed to twinkle. While this specific application of detecting breast biopsy markers by twinkling may provide a way to approach an unmet need in the care of patients with breast cancer, this work ultimately provides a platform from which the keys to unlocking the fundamental physics of twinkling can be rigorously explored.

A New Plane Wave Compounding Scheme Using Phase Compensation for Motion Detection.

Plane wave (PW) transmission has enabled multiple new applications, such as shear wave elastography, ultrafast Doppler imaging, and functional ultrasound imaging. PW compounding (PWC), which coherently sums the echo signals from multiple PW transmits with different angles, is widely used to improve B-mode image quality. When the motion between two speckle images is estimated, PWC suffers from an inherent displacement estimation error. This is derived theoretically and experimentally demonstrated in this work. We show that the phase difference between the acquired data with PW emissions with different angles is related to this error. When the absolute value of the phase difference is larger than π /2, the displacement estimation error occurs. A new scheme, named initial-phase-compensated PWC (IPCPWC), is proposed, which compensates the phase of echo signals from each PW transmit and maintains the absolute value of the phase difference smaller than π /2. The increased signal-to-noise ratio and reduced jitter of IPCPWC in motion data are demonstrated using tissue mimicking phantoms compared with PWC.

Toward improved accuracy in shear wave elastography of arteries through controlling the arterial response to ultrasound perturbation in-silico and in phantoms.

Dispersion-based inversion has been proposed as a viable direction for materials characterization of arteries, allowing clinicians to better study cardiovascular conditions using shear wave elastography. However, these methods rely ona prioriknowledge of the vibrational modes dominating the propagating waves induced by acoustic radiation force excitation: differences between anticipated and real modal content are known to yield errors in the inversion. We seek to improve the accuracy of this process by modeling the artery as a fluid-immersed cylindrical waveguide and building an analytical framework to prescribe radiation force excitations that will selectively excite certain waveguide modes using ultrasound acoustic radiation force. We show that all even-numbered waveguide modes can be eliminated from the arterial response to perturbation, and confirm the efficacy of this approach within silicotests that show that odd modes are preferentially excited. Finally, by analyzing data from phantom tests, we find a set of ultrasound focal parameters that demonstrate the viability of inducing the desired odd-mode response in experiments.

Using Ultrasound Color Doppler Twinkling to Identify Biopsy Markers in the Breast and Axilla.

In breast radiology, ultrasound detection of biopsy markers or clips for localization purposes is often challenging, especially in the axilla. The purpose of this research was to test the hypothesis that the surface roughness of biopsy clips would elicit a twinkling signature on color Doppler, making them more readily identifiable by ultrasound. Ultrasound color Doppler imaging of 12 biopsy markers was performed and consensus scoring of the degree of twinkling (0 [no twinkling] to 4 [exuberant twinkling]) was obtained for each of the markers. The surface roughness characteristics of the markers were measured using 3-D coherence scanning interferometry. The 3 markers scoring at least 3 for twinkling in vitro were cork, Q and Vision. Of these 3 markers, only the cork marker scored a 4 ex vivo and in cadaveric tissue. Surface roughness metrics demonstrated a positive estimated correlation with the twinkling scores (rho = 0.33, 95% CI = [-0.48 to 0.84]). Of the 12 markers tested, the markers that twinkled corresponded to surface roughness measured with non-contact 3-D optical imaging. Qualitatively, lower color scales and color frequencies optimized twinkling, but the most specific qualitative predictor of confidence in twinkling was insensitivity to changes in color scale and color frequency values.

Multimodal guided wave inversion for arterial stiffness: methodology and validation in phantoms.

Arterial stiffness is an important biomarker for many cardiovascular diseases. Shear wave elastography is a recent technique aimed at estimating local arterial stiffness using guided wave inversion (GWI), i.e. matching the computed and measured wave dispersion. This paper develops and validates a new GWI approach by synthesizing various recent observations and algorithms: (a) refinements to signal processing to obtain more accurate experimental dispersion curves; (b) an efficient forward model to compute theoretical dispersion curves for immersed, incompressible cylindrical waveguides; (c) an optimization framework based on the recent observation that the measured dispersion curve is multimodal, i.e. it matches for not one but two different wave modes in two different frequency ranges. The resulting inversion approach is validated using extensive experimental data from rubber tube phantoms, not only for modulus estimation but also to simultaneously estimate modulus and wall thickness. The observations indicate that the modulus estimates are best performed with the information on wall thickness. The approach, which takes less than half a minute to run, is shown to be accurate, with the modulus estimated with less than 4% error for 70% of the experiments.

Time-Aligned Plane Wave Compounding Methods for High-Frame-Rate Shear Wave Elastography: Experimental Validation and Performance Assessment on Tissue Phantoms.

Shear wave elastography (SWE) is an ultrasonic technique able to quantitatively assess the mechanical properties of tissues by combining acoustic radiation force and ultrafast imaging. While utilizing coherent plane wave compounding enhances echo and shear wave motion signal-to-noise ratio (SNR), it also reduces the effective pulse repetition frequency (PRFe), affecting the accuracy of the measurements of motion and, consequently, of material properties. It is important to maintain both high-motion SNR and PRFe, particularly for the characterization of (material and/or geometrical) dispersive tissues such as arteries. This work proposes a method for SWE measurements with high SNR, while maintaining a high PRFe, using conventional clinical ultrasound scanners. A time alignment process is applied after acquiring data from plane wave transmissions at different angles. The time alignment uses interpolation to obtain data points at higher frame rates, and the time-aligned data are compounded to increase the SNR. The method is used for SWE in tissue-mimicking phantoms of different stiffness and is compared with traditional plane wave compounding. Increases of 58% and 36% in spatial and temporal bandwidth compared with conventional plane wave compounding, respectively, can be achieved for SWE measurements of representative arterial stiffness values. Improvements in phase velocity accuracy and bandwidth in an arterial phantom are also described, to emphasize the beneficial advantage in dispersive cases.

Three-dimensional shear wave elastography on conventional ultrasound scanners with external vibration.

Two-dimensional (2D) ultrasound shear wave elastography (SWE) has been widely used for soft tissue properties assessment. Given that shear waves propagate in three dimensions (3D), extending SWE from 2D to 3D is important for comprehensive and accurate stiffness measurement. However, implementation of 3D SWE on a conventional ultrasound scanner is challenging due to the low volume rate (tens of Hertz) associated with limited parallel receive capability of the scanner's hardware beamformer. Therefore, we developed an external mechanical vibration-based 3D SWE technique allowing robust 3D shear wave tracking and speed reconstruction for conventional scanners. The aliased shear wave signal detected with a sub-Nyquist sampling frequency was corrected by leveraging the cyclic nature of the sinusoidal shear wave generated by the external vibrator. Shear wave signals from different sub-volumes were aligned in temporal direction to correct time delays from sequential pulse-echo events, followed by 3D speed reconstruction using a 3D local frequency estimation algorithm. The technique was validated on liver fibrosis phantoms with different stiffness, showing good correlation (r = 0.99, p < 0.001) with values measured from a state-of-the-art SWE system (GE LOGIQ E9). The phantoms with different stiffnesses can be well-differentiated regardless of the external vibrator position, indicating the feasibility of the 3D SWE with regard to different shear wave propagation scenarios. Finally, shear wave speed calculated by the 3D method correlated well with magnetic resonance elastography performed on human liver (r = 0.93, p = 0.02), demonstrating the in vivo feasibility. The proposed technique relies on low volume rate imaging and can be implemented on the widely available clinical ultrasound scanners, facilitating its clinical translation to improve liver fibrosis evaluation.

In Vivo Open- and Closed-chest Measurements of Left-Ventricular Myocardial Viscoelasticity using Lamb wave Dispersion Ultrasound Vibrometry (LDUV): A Feasibility Study.

Diastolic dysfunction causes close to half of congestive heart failures and is associated with increased stiffness in left-ventricular myocardium. A clinical tool capable of measuring viscoelasticity of the myocardium could be beneficial in clinical settings. We used Lamb wave Dispersion Ultrasound Vibrometry (LDUV) for assessing the feasibility of making in vivo non-invasive measurements of myocardial elasticity and viscosity in pigs. In vivo open-chest measurements of myocardial elasticity and viscosity obtained using a Fourier space based analysis of Lamb wave dispersion are reported. The approach was used to perform ECG-gated transthoracic in vivo measurements of group velocity, elasticity and viscosity throughout a single heart cycle. Group velocity, elasticity and viscosity in the frequency range 50-500 Hz increased from diastole to systole, consistent with contraction and relaxation of the myocardium. Systolic group velocity, elasticity and viscosity were 5.0 m/s, 19.1 kPa, 6.8 Pa·s, respectively. In diastole, the measured group velocity, elasticity and viscosity were 1.5 m/s, 5.1 kPa and 3.2 Pa·s, respectively.

Probe Oscillation Shear Wave Elastography: Initial In Vivo Results in Liver.

Shear wave elastography methods are able to accurately measure tissue stiffness, allowing these techniques to monitor the progression of hepatic fibrosis. While many methods rely on acoustic radiation force to generate shear waves for 2-D imaging, probe oscillation shear wave elastography (PROSE) provides an alternative approach by generating shear waves through continuous vibration of the ultrasound probe while simultaneously detecting the resulting motion. The generated shear wave field in in vivo liver is complicated, and the amplitude and quality of these shear waves can be influenced by the placement of the vibrating probe. To address these challenges, a real-time shear wave visualization tool was implemented to provide instantaneous visual feedback to optimize probe placement. Even with the real-time display, it was not possible to fully suppress residual motion with established filtering methods. To solve this problem, the shear wave signal in each frame was decoupled from motion and other sources through the use of a parameter-free empirical mode decomposition before calculating shear wave speeds. This method was evaluated in a phantom as well as in in vivo livers from five volunteers. PROSE results in the phantom as well as in vivo liver correlated well with independent measurements using the commercial General Electric Logiq E9 scanner.

Application of Acoustoelasticity to Evaluate Nonlinear Modulus in Ex Vivo Kidneys.

Currently, dynamic elastography techniques estimate the linear elastic shear modulus of different body tissues. New methods that investigate other properties of soft tissues such as anisotropy, viscosity, and shear nonlinearity would provide more information about the structure and function of the tissue and might provide a better contrast than tissue stiffness and hence provide more effective diagnostic tools for some diseases. It has previously been shown that shear wave velocity in a medium changes due to an applied stress, a phenomenon called acoustoelasticity (AE). Applying a stress to compress a medium while measuring the shear wave velocity versus strain provides data with which the third-order nonlinear shear modulus can be estimated. To evaluate the feasibility of estimating , we evaluated ten ex vivo porcine kidneys embedded in 10% porcine gelatin to mimic the case of a transplanted kidney. Under assumptions of an elastic incompressible medium for AE measurements, the shear modulus was quantified at each compression level and the applied strain was assessed by measuring the change in the thickness of the kidney cortex. Finally, was calculated by applying the AE theory. Our results demonstrated that it is possible to estimate a nonlinear shear modulus by monitoring the changes in strain and due to kidney deformation. The magnitudes of are higher when the compression is performed progressively and when using a plate attached to the transducer. Nevertheless, the values obtained for are similar to those previously reported in the literature for breast tissue.

Acoustic Radiation Force-Induced Creep-Recovery (ARFICR): A Noninvasive Method to Characterize Tissue Viscoelasticity.

Ultrasound shear wave elastography is a promising noninvasive, low cost, and clinically viable tool for liver fibrosis staging. Current shear wave imaging technologies on clinical ultrasound scanners ignore shear wave dispersion and use a single group velocity measured over the shear wave bandwidth to estimate tissue elasticity. The center frequency and bandwidth of shear waves induced by acoustic radiation force depend on the ultrasound push beam (push duration, -number, etc.) and the viscoelasticity of the medium, and therefore are different across scanners from different vendors. As a result, scanners from different vendors may give different tissue elasticity measurements within the same patient. Various methods have been proposed to evaluate shear wave dispersion to better estimate tissue viscoelasticity. A rheological model such as the Kelvin-Voigt model is typically fitted to the shear wave dispersion to solve for the elasticity and viscosity of tissue. However, these rheological models impose strong assumptions about frequency dependence of elasticity and viscosity. Here, we propose a new method called Acoustic Radiation Force Induced Creep-Recovery (ARFICR) capable of quantifying rheological model-independent measurements of elasticity and viscosity for more robust tissue health assessment. In ARFICR, the creep-recovery time signal at the focus of the push beam is used to calculate the relative elasticity and viscosity (scaled by an unknown constant) over a wide frequency range. Shear waves generated during the ARFICR measurement are also detected and used to calculate the shear wave velocity at its center frequency, which is then used to calibrate the relative elasticity and viscosity to absolute elasticity and viscosity. In this paper, finite-element method simulations and experiments in tissue mimicking phantoms are used to validate and characterize the extent of viscoelastic quantification of ARFICR. The results suggest that ARFICR can measure tissue viscoelasticity reliably. Moreover, the results showed the strong frequency dependence of viscoelastic parameters in tissue mimicking phantoms and healthy liver.

An Ultrasound Surface Wave Technique for Assessing Skin and Lung Diseases.

Systemic sclerosis (SSc) is a multi-organ connective tissue disease characterized by immune dysregulation and organ fibrosis. Severe organ involvement, especially of the skin and lung, is the cause of morbidity and mortality in SSc. Interstitial lung disease (ILD) includes multiple lung disorders in which the lung tissue is fibrotic and stiffened. The purpose of this study was to translate ultrasound surface wave elastography (USWE) for assessing patients with SSc and/or ILD via measuring surface wave speeds of both skin and superficial lung tissue. Forty-one patients with both SSc and ILD and 30 healthy patients were enrolled in this study. An external harmonic vibration was used to generate the wave propagation on the skin or lung. Three excitation frequencies of 100, 150 and 200 Hz were used. An ultrasound probe was used to measure the wave propagation in the tissue non-invasively. Surface wave speeds were measured on the forearm and upper arm of both left and right arm, as well as the upper and lower lungs, through six intercostal spaces of patients and healthy patients. Viscoelasticity of the skin was calculated by the wave speed dispersion with frequency using the Voigt model. The magnitudes of surface wave speed and viscoelasticity of patients' skin were significantly higher than those of healthy patients (p <0.0001) for each location and each frequency. The surface wave speeds of patients' lung were significantly higher than those of healthy patients (p <0.0001) for each location and each frequency. USWE is a non-invasive and non-ionizing technique for measuring both skin and lung surface wave speed and may be useful for quantitative assessment of SSc and/or ILD.