The spatio-temporal strain response of oedematous and nonoedematous tissue to sustained compression in vivo

Pitting Is Not Only a Measure of Oedema Presence: Using High-Frequency Ultrasound to Guide Pitting Test Standardisation for Assessment of Lymphoedema

The pitting qualities of lymphoedema tissue change with disease progression. However, little is known about the underlying tissue response to the pitting test or the tissue characteristics that enhance or resist indentation. The pitting test is currently unstandardised, and the influence of test technique on pitting outcomes is unknown. Understanding how tissue reacts to applied pressure will build evidence for the standardisation of the pitting test. Ninety pitting test sites from fifteen women with unilateral breast cancer-related lymphoedema were evaluated using high-frequency ultrasound (HFUS), bioelectrical impedance spectroscopy (BIS), and limb volume measures. Three sites on each lymphoedema and non-lymphoedema arm were subject to a 60-s (s) staged pitting test, with changes in tissue features captured with ultrasound imaging before, throughout, and after the pitting test. Pitting qualities of tissues varied greatly, with lymphoedema sites pitting more frequently (p < 0.001) with greater depth (p < 0.001) and requiring a longer recovery time (p = 0.002) than contralateral unaffected tissue. Pitting is not solely attributable to oedema volume. Non-structural and structural characteristics of dermal and subcutaneous layers also influence tissue responses to sustained pressure. To enhance the validity and reliability of pitting assessment, a 60 s staged pitting test with an observation of tissue recovery is recommended for lymphoedema presentations.

Acoustic Radiation Force Impulse Elastography Assessment of Lymphoedema Tissue: An Insight into Tissue Stiffness

Palpation remains essential for evaluating lymphoedema to detect subtle changes that may indicate progression. As palpation sense is not quantifiable, this study investigates the utility of ultrasound elastography to quantify stiffness of lymphoedema tissue and explore the influence of the pitting test on tissue stiffness. Fifteen women with unilateral arm lymphoedema were scanned using a Siemens S3000 Acuson ultrasound (Siemens, Germany) with 18 MHz and 9 MHz linear transducers to assess tissue structure and tissue stiffness with Acoustic Radiation Force Impulse elastography. Ninety sites were assessed, three on each of the lymphoedema-affected and contralateral unaffected arms. A subgroup of seven lymphoedema-affected sites included additional elastography imaging after a 60-s pitting test. Dermal tissue stiffness was greater than subcutaneous tissue stiffness regardless of the presence of pathology (p < 0.001). Lymphoedematous tissue exhibited a higher dermal to subcutaneous tissue stiffness ratio than contralateral sites (p = 0.005). Subgroup analysis indicated that the pitting test reduces dermal tissue stiffness (p = 0.018) and may alter the stiffness of the subcutaneous tissue layer. Elastography demonstrates potential as a complement to lymphoedema palpation assessment. The novel pre-test and post-pitting elastography imaging protocol yielded information representative of lymphoedema tissue characteristics that could not be ascertained from pre-test elastography images alone.

Estimation of Mechanical and Transport Parameters in Cancers Using Short Time Poroelastography

Mechanical and transport properties of cancers such as Young's modulus (YM), Poisson's ratio (PR), and vascular permeability (VP) have great clinical importance in cancer diagnosis, prognosis, and treatment. However, non-invasive estimation of these parameters in vivo is challenged by many practical factors. Elasticity imaging methods, such as "poroelastography", require prolonged data acquisition, which can limit their clinical applicability. In this paper, we investigate a new method to perform poroelastography experiments, which results in shorter temporal acquisition windows. This method is referred to as "short-time poroelastography" (STPE). Finite element (FE) and ultrasound simulations demonstrate that, using STPE, it is possible to accurately estimate YM, PR (within 10% error) using windows of observation (WoOs) of length as short as 1 underlying strain Time Constant (TC). The error was found to be almost negligible (< 3%) when using WoOs longer than 2 strain TCs. In the case of VP estimation, WoOs of at least 2 strain TCs are required to obtain an error < 8% (in simulations). The stricter requirement for the estimation of VP with respect to YM and PR is due its reliance on the transient strain behavior while YM and PR depend on the steady state strain values only. In vivo experimental data are used as a proof-of-principle of the potential applicability of the proposed methodology in vivo. The use of STPE may provide a means to efficiently perform poroelastography experiments without compromising the accuracy of the estimated tissue properties.

Estimation of Vascular Permeability in Irregularly Shaped Cancers Using Ultrasound Poroelastography

OBJECTIVE

Vascular permeability (VP) is a mechanical parameter which plays an important role in cancer initiation, metastasis, and progression. To date, there are only a few non-invasive methods that can be used to image VP in solid tumors. Most of these methods require the use of contrast agents and are expensive, limiting widespread use.

METHODS

In this paper, we propose a new method to image VP in tumors, which is based on the use of ultrasound poroelastography. Estimation of VP by poroelastography requires knowledge of the Young's modulus (YM), Poisson's ratio (PR), and strain time constant (TC) in the tumors. In our method, we find the ellipse which best fits the tumor (regardless of its shape) using an eigen-system-based fitting technique and estimate the YM and PR using Eshelby's elliptic inclusion formulation. A Fourier method is used to estimate the axial strain TC, which does not require any initial guess and is highly robust to noise.

RESULTS

It is demonstrated that the proposed method can estimate VP in irregularly shaped tumors with an accuracy of above [Formula: see text] using ultrasound simulation data with signal-to-noise ratio of 20 dB or higher. In vivo feasibility of the proposed technique is demonstrated in an orthotopic mouse model of breast cancer.

CONCLUSION

The proposed imaging method can provide accurate and localized estimation of VP in cancers non-invasively and cost-effectively.

SIGNIFICANCE

Accurate and non-invasive assessment of VP can have a significant impact on diagnosis, prognosis, and treatment of cancers.

Estimation of mechanical parameters in cancers by empirical orthogonal function analysis of poroelastography data

Ultrasound poroelastography is a non-invasive imaging modality that has been shown to be capable of estimating mechanical parameters such as Young's modulus (YM), Poisson's ratio (PR) and vascular permeability (VP) in cancers. However, experimental poroelastographic data are inherently noisy because of the requirement of relatively long temporal data acquisitions often in hand-held mode conditions. In this paper, we propose a new method, which allows accurate estimation of YM and PR from denoised steady state axial and lateral strains by empirical orthogonal function (EOF) analysis of poroelastographic data. The method also allows estimation of VP from the time constant (TC) of the first expansion coefficient (EC) of the temporal axial strain, which has larger dynamic range and lower noise in comparison to the actual temporal axial strain curve. We validated our technique through finite element (FE) and ultrasound simulations and tested the in vivo feasibility in experimental data obtained from a cancer animal model. The percent relative errors (PRE) in the estimation of YM, PR and VP using the EOF analysis as applied to ultrasound simulation data were 3.27%, 3.10%, 14.22%, respectively (at SNR of 20 dB). Based on the high level of accuracy by EOF analysis, the proposed technique may become a useful signal processing technique for applications focusing on the estimation of the mechanical behavior of cancers.

Correlation Between Renal Cortical Stiffness and Histological Determinants by Point Shear-Wave Elastography in Patients With Kidney Transplantation

INTRODUCTION

Renal allograft biopsy is the gold standard for the detection of histological lesions of chronic allograft dysfunction. The identification of a noninvasive routine test would be desirable. Elastosonography is used to assess tissue stiffness according to viscosity, and no data are available on the use of point quantification shear-wave elastography (ElastPQ) for the evaluation of renal chronic lesions.

RESEARCH QUESTION

To evaluate the feasibility of ElastPQ to assess cortical allograft stiffness and to determine the correlation of clinical, biological, and pathological factors with the diagnostic accuracy of kidney stiffness values in patients with histological lesions.

DESIGN

Forty-two patients underwent kidney transplant biopsy and 10 valid measurements of ElastPQ, blindly performed by 2 operators. The interobserver reproducibility was assessed according to intraclass correlation coefficient. The ElastPQ measurements and the clinical data were compared using the Spearman correlation analysis.

RESULTS

97.6% reliable measurements were obtained using ElastPQ, with an excellent interobserver agreement. The kidney stiffness was significantly higher in the patients with a time since transplantation >12 months and was correlated with chronic lesions (interstitial fibrosis, tubular atrophy transplant glomerulopathy, and mesangial matrix), with the interstitial fibrosis/tubular atrophy, score and with the sum of the scores of the chronic lesions. Mesangial matrix increase is the only independent determinant of kidney stiffness.

DISCUSSION

ElastPQ is a noninvasive, reproducible, and sensitive diagnostic tool able to detect moderate/severe chronic lesions. Its routine use during follow-up can identify patients eligible for biopsy, which remains the gold standard exam for detecting chronic allograft dysfunction.

Assessment of Mechanical Properties of Tissue in Breast Cancer-Related Lymphedema Using Ultrasound Elastography

Breast cancer-related lymphedema is a consequence of a malfunctioning lymphatic drainage system resulting from surgery or some other form of treatment. In the initial stages, minor and reversible increases in the fluid volume of the arm are evident. As the stages progress over time, the underlying pathophysiology dramatically changes with an irreversible increase in arm volume most likely due to a chronic local inflammation leading to adipose tissue hypertrophy and fibrosis. Clinicians have subjective ways to stage the degree and severity such as the pitting test which entails manually comparing the elasticity of the affected and unaffected arms. Several imaging modalities can be used but ultrasound appears to be the most preferred because it is affordable, safe, and portable. Unfortunately, ultrasonography is not typically used for staging lymphedema, because the appearance of the affected and unaffected arms is similar in B-mode ultrasound images. However, novel ultrasound techniques have emerged, such as elastography, which may be able to identify changes in mechanical properties of the tissue related to detection and staging of lymphedema. This paper presents a novel technique to compare the mechanical properties of the affected and unaffected arms using quasi-static ultrasound elastography to provide an objective alternative to the current subjective assessment. Elastography is based on time delay estimation (TDE) from ultrasound images to infer displacement and mechanical properties of the tissue. We further introduce a novel method for TDE by incorporating higher order derivatives of the ultrasound data into a cost function and propose a novel optimization approach to efficiently minimize the cost function. This method works reliably with our challenging patient data. We collected radio frequency ultrasound data from both arms of seven patients with stage 2 lymphedema, at six different locations in each arm. The ratio of strain in skin, subcutaneous fat, and skeletal muscle divided by strain in the standoff gel pad was calculated in the unaffected and affected arms. The p -values using a Wilcoxon sign-rank test for the skin, subcutaneous fat, and skeletal muscle were 1.24×10^-5 , 1.77×10^-8 , and 8.11×10^-7 respectively, showing differences between the unaffected and affected arms with a very high level of significance.

Quantitative characterization of mechanically indented in vivo human skin in adults and infants using optical coherence tomography

Influenced by both the intrinsic viscoelasticity of the tissue constituents and the time-evolved redistribution of fluid within the tissue, the biomechanical response of skin can reflect not only localized pathology but also systemic physiology of an individual. While clinical diagnosis of skin pathologies typically relies on visual inspection and manual palpation, a more objective and quantitative approach for tissue characterization is highly desirable. Optical coherence tomography (OCT) is an interferometry-based imaging modality that enables in vivo assessment of cross-sectional tissue morphology with micron-scale resolution, which surpasses those of most standard clinical imaging tools, such as ultrasound imaging and magnetic resonance imaging. This pilot study investigates the feasibility of characterizing the biomechanical response of in vivo human skin using OCT. OCT-based quantitative metrics were developed and demonstrated on the human subject data, where a significant difference between deformed and nondeformed skin was revealed. Additionally, the quantified postindentation recovery results revealed differences between aged (adult) and young (infant) skin. These suggest that OCT has the potential to quantitatively assess the mechanically perturbed skin as well as distinguish different physiological conditions of the skin, such as changes with age or disease.

Design and Testing of a Single-Element Ultrasound Viscoelastography System for Point-of-Care Edema Quantification

Management of fluid overload in patients with end-stage renal disease represents a unique challenge to clinical practice because of the lack of accurate and objective measurement methods. Currently, peripheral edema is subjectively assessed by palpation of the patient's extremities, ostensibly a qualitative indication of tissue viscoelastic properties. New robust quantitative estimates of tissue fluid content would allow clinicians to better guide treatment, minimizing reactive treatment decision making. Ultrasound viscoelastography (UVE) can be used to estimate strain in viscoelastic tissue, deriving material properties that can help guide treatment. We are developing and testing a simple, low-cost UVE system using a single-element imaging transducer that is simpler and less computationally demanding than array-based systems. This benchtop validation study tested the feasibility of using the UVE system by measuring the mechanical properties of a tissue-mimicking material under large strains. We generated depth-dependent creep curves and viscoelastic parameter maps of time constants and elastic moduli for the Kelvin model of viscoelasticity. During testing, the UVE system performed well, with mean UVE-measured strain matching standard mechanical testing with maximum absolute errors ≤4%. Motion tracking revealed high correlation and signal-to-noise ratios, indicating that the system is reliable.

Review of ultrasound image guidance in external beam radiotherapy part II: intra-fraction motion management and novel applications

Imaging has become an essential tool in modern radiotherapy (RT), being used to plan dose delivery prior to treatment and verify target position before and during treatment. Ultrasound (US) imaging is cost-effective in providing excellent contrast at high resolution for depicting soft tissue targets apart from those shielded by the lungs or cranium. As a result, it is increasingly used in RT setup verification for the measurement of inter-fraction motion, the subject of Part I of this review (Fontanarosa et al 2015 Phys. Med. Biol. 60 R77-114). The combination of rapid imaging and zero ionising radiation dose makes US highly suitable for estimating intra-fraction motion. The current paper (Part II of the review) covers this topic. The basic technology for US motion estimation, and its current clinical application to the prostate, is described here, along with recent developments in robust motion-estimation algorithms, and three dimensional (3D) imaging. Together, these are likely to drive an increase in the number of future clinical studies and the range of cancer sites in which US motion management is applied. Also reviewed are selections of existing and proposed novel applications of US imaging to RT. These are driven by exciting developments in structural, functional and molecular US imaging and analytical techniques such as backscatter tissue analysis, elastography, photoacoustography, contrast-specific imaging, dynamic contrast analysis, microvascular and super-resolution imaging, and targeted microbubbles. Such techniques show promise for predicting and measuring the outcome of RT, quantifying normal tissue toxicity, improving tumour definition and defining a biological target volume that describes radiation sensitive regions of the tumour. US offers easy, low cost and efficient integration of these techniques into the RT workflow. US contrast technology also has potential to be used actively to assist RT by manipulating the tumour cell environment and by improving the delivery of radiosensitising agents. Finally, US imaging offers various ways to measure dose in 3D. If technical problems can be overcome, these hold potential for wide-dissemination of cost-effective pre-treatment dose verification and in vivo dose monitoring methods. It is concluded that US imaging could eventually contribute to all aspects of the RT workflow.

Investigation of In Vivo skin stiffness anisotropy in breast cancer related lymphoedema

There is a limited range of suitable measurement techniques for detecting and assessing breast cancer related lymphoedema (BCRL). This study investigated the suitability of using skin stiffness measurements, with a particular focus on the variation in stiffness with measurement direction (known as anisotropy). In addition to comparing affected tissue with the unaffected tissue on the corresponding site on the opposite limb, volunteers without BCRL were tested to establish the normal variability in stiffness anisotropy between these two corresponding regions of skin on each opposite limb. Multi-directional stiffness was measured with an Extensometer, within the higher stiffness region that skin typically displays at high applied strains, using a previously established protocol developed by the authors. Healthy volunteers showed no significant difference in anisotropy between regions of skin on opposite limbs (mean decrease of 4.7 +/-2.5% between non-dominant and dominant arms), whereas BCRL sufferers showed a significant difference between limbs (mean decrease of 51.0+/-16.3% between unaffected and affected arms). A large difference in anisotropy was apparent even for those with recent onset of the condition, indicating that the technique may have potential to be useful for early detection. This difference also appeared to increase with duration since onset. Therefore, measurement of stiffness anisotropy has potential value for the clinical assessment and diagnosis of skin conditions such as BCRL. The promising results justify a larger study with a larger number of participants.

Effect of temporal acquisition parameters on image quality of strain time constant elastography

Ultrasound methods to image the time constant (TC) of elastographic tissue parameters have been recently developed. Elastographic TC images from creep or stress relaxation tests have been shown to provide information on the viscoelastic and poroelastic behavior of tissues. However, the effect of temporal ultrasonic acquisition parameters and input noise on the image quality of the resultant strain TC elastograms has not been fully investigated yet. Understanding such effects could have important implications for clinical applications of these novel techniques. This work reports a simulation study aimed at investigating the effects of varying windows of observation, acquisition frame rate, and strain signal-to-noise ratio (SNR) on the image quality of elastographic TC estimates. A pilot experimental study was used to corroborate the simulation results in specific testing conditions. The results of this work suggest that the total acquisition time necessary for accurate strain TC estimates has a linear dependence to the underlying strain TC (as estimated from the theoretical strain-vs.-time curve). The results also indicate that it might be possible to make accurate estimates of the elastographic TC (within 10% error) using windows of observation as small as 20% of the underlying TC, provided sufficiently fast acquisition rates (>100 Hz for typical acquisition depths). The limited experimental data reported in this study statistically confirm the simulation trends, proving that the proposed model can be used as upper bound guidance for the correct execution of the experiments.

Ultrasound strain elastography in assessment of cortical mechanical behavior in acute renal vein occlusion: in vivo animal model

To assess the correlation of quantitative ultrasound strain parameters with the severity of cortical edema in renal vein occlusion, we prospectively performed ultrasound strain elastography on a canine acute renal vein occlusion model prior to and following 10, 20, and 40min of renal vein ligation. Strain and strain relaxation time representing the deformation and relaxation of the renal cortices and reference soft tissue were produced by the external compression with the ultrasound transducer and estimated using commercially available 2-D speckle tracking software. Cortical thickness was additionally measured. Repeated-measures analysis of variance was used to examine the difference in cortical thickness, strain ratio (mean cortical strain divided by mean reference tissue strain), and strain relaxation time ratio (cortical relaxation time divided by reference tissue relaxation time) prior to and after renal vein ligation. Pearson's correlation coefficient was applied to test the relationship between strain parameters and the time of the renal vein ligation. There was a strong positive correlation between the duration of renal vein ligation and strain (R(2)=0.97) and strain relaxation time (R(2)=0.98) ratios. Significant differences in strain and strain relaxation time ratios were found at all measured timepoints (all P≪.001). Cortical thickness, however, showed no significant difference between timepoints (P=.065). Our result suggest that strain and strain relaxation time ratios may be used as quantitative markers for the assessment of the renal cortical mechanical behavior in subclinical acute renal vein occlusion.

Ultrasound strain zero-crossing elasticity measurement in assessment of renal allograft cortical hardness: a preliminary observation

To determine whether ultrasound strain zero-crossing elasticity measurement can be used to discriminate moderate cortical fibrosis or inflammation in renal allografts, we prospectively assessed cortical hardness with quasi-static ultrasound elastography in 38 renal transplant patients who underwent kidney biopsy from January 2013 to June 2013. With the Banff score criteria for renal cortical fibrosis as gold standard, 38 subjects were divided into two groups: group 1 (n = 18) with ≤25% cortical fibrosis and group 2 (n = 20) with >26% cortical fibrosis. We then divided this population again into group 3 (n = 20) with ≤ 25% inflammation and group 4 (n = 18) with >26% inflammation based on the Banff score for renal parenchyma inflammation. To estimate renal cortical hardness in both population divisions, we propose an ultrasound strain relative zero-crossing elasticity measurement (ZC) method. In this technique, the relative return to baseline, that is zero strain, of strain in the renal cortex is compared with that of strain in reference soft tissue (between the abdominal wall and pelvic muscles). Using the ZC point on the reference strain decompression slope as standard, we determined when cortical strain crossed zero during decompression. ZC was negative when cortical strain did not return or returned after the reference, whereas ZC was positive when cortical strain returned ahead of the reference. Fisher's exact test was used to examine the significance of differences in ZC between groups 1 and 2 and between groups 3 and 4. The accuracy of ZC in determining moderate cortical fibrosis and moderate inflammation was examined by receiver operating characteristic analysis. The intra-class correlation coefficient and analysis of variance were used to test inter-rater reliability and reproducibility. ZC had good inter-observer agreement (ICC = 0.912) and reproducibility (p = 0.979). ZCs were negative in 18 of 18 cases in group 1 and positive in 19 of 20 cases in group 2 (p ≪ 0.001), and were positive in 18 of 20 cases in group 3 and negative in 17 of 18 cases in group 4 (p ≪ 0.001). The area under the receiver operating characteristic curve was 0.992 ± 0.010 for fibrosis and 0.988 ± 0.021 for inflammation. ZC had 100% sensitivity and 95% specificity when zero strain was used as the cutoff value to determine moderate cortical fibrosis and 94% sensitivity and 90% specificity for inflammation. ZC is a new strain marker that could be straightforward to interpret and perform, making it a potentially practical approach for monitoring progression of cortical fibrosis or inflammation in renal allografts.

Effect of permeability on the performance of elastographic imaging techniques

Elastography is a well-established imaging modality. While a number of studies aimed at evaluating the performance of elastographic techniques are retrievable in the literature, very little information is available on the effects that the presence of an underlying permeability contrast in the tissue may have on the resulting elastograms. Permeability is a fundamental tissue parameter, which characterizes the ease with which fluid can move within a tissue. This parameter plays a central role both biomechanically in the description of the temporal behavior of fluid-filled tissues and clinically in the development of a number of diagnostic and therapeutic modalities. In this paper, we present a simulation study that investigates selected elastographic image quality factors in nonhomogeneous materials, modeled as poroelastic media with different geometries and permeability contrasts. The results of this study indicate that the presence of an underlying permeability contrast may create a new contrast mechanism in the spatial and temporal distributions of the axial strains and the effective Poisson's ratios experienced by the tissue and as imaged by the corresponding elastograms. The effect of permeability on the elastographic image quality factors analyzed in this study was found to be a nonsymmetric function of the underlying mechanical contrast between background and target, the geometry of the material and the boundary conditions.

Medical ultrasound: imaging of soft tissue strain and elasticity

After X-radiography, ultrasound is now the most common of all the medical imaging technologies. For millennia, manual palpation has been used to assist in diagnosis, but it is subjective and restricted to larger and more superficial structures. Following an introduction to the subject of elasticity, the elasticity of biological soft tissues is discussed and published data are presented. The basic physical principles of pulse-echo and Doppler ultrasonic techniques are explained. The history of ultrasonic imaging of soft tissue strain and elasticity is summarized, together with a brief critique of previously published reviews. The relevant techniques-low-frequency vibration, step, freehand and physiological displacement, and radiation force (displacement, impulse, shear wave and acoustic emission)-are described. Tissue-mimicking materials are indispensible for the assessment of these techniques and their characteristics are reported. Emerging clinical applications in breast disease, cardiology, dermatology, gastroenterology, gynaecology, minimally invasive surgery, musculoskeletal studies, radiotherapy, tissue engineering, urology and vascular disease are critically discussed. It is concluded that ultrasonic imaging of soft tissue strain and elasticity is now sufficiently well developed to have clinical utility. The potential for further research is examined and it is anticipated that the technology will become a powerful mainstream investigative tool.

What challenges must be overcome before ultrasound elasticity imaging is ready for the clinic?

Ultrasound elasticity imaging has been a research interest for the past 20 years with the goal of generating novel images of soft tissues based on their material properties (i.e., stiffness and viscosity). The motivation for such an imaging modality lies in the fact that many soft tissues can share similar ultrasonic echogenicities, but may have very different mechanical properties that can be used to clearly visualize normal anatomy and delineate diseased tissues and masses. Recently, elasticity imaging techniques have moved from the laboratory to the clinical setting, where clinicians are beginning to characterize tissue stiffness as a diagnostic metric and commercial implementations of ultrasonic elasticity imaging are beginning to appear on the market. This article provides a foundation for elasticity imaging, an overview of current research and commercial realizations of elasticity imaging technology and a perspective on the current successes, limitations and potential for improvement of these imaging technologies.

Performance analysis of a new real-time elastographic time constant estimator

New elastographic techniques such as poroelastography and viscoelasticity imaging aim at imaging the temporal mechanical behavior of tissues. These techniques usually involve the use of curve fitting methods being applied to noisy data to estimate new elastographic parameters. As of today, however, current elastographic implementations of poroelastography and viscoelasticity imaging methods are in general too slow and not optimized for clinical applications. Furthermore, image quality performance of these new elastographic techniques is still largely unknown due to a paucity of data and the lack of systematic studies that analyze their performance limitations. In this paper, we propose a new elastographic time constant (TC) estimator, which is based on the use of the least square error (LSE) curve-fitting method and the Levenberg-Marquardt (LM) optimization rule as applied to noisy elastographic data obtained from a material in a creep-type experiment. The algorithm is executed on a massively parallel general purpose graphics processing unit (GPGPU) to achieve real-time performance. The estimator's performance is analyzed using simulations. Experimental results obtained from poroelastic phantoms are presented as a proof of principle of the new estimator's technical applicability on real experimental data. The results of this study demonstrate that the newly proposed elastographic estimator can produce highly accurate and sensitive elastographic TC estimates in real-time and at high signal-to-noise ratios.