# How is it possible for an Ultrasound device to correctly interpret a negative density change in tissue?

I understand the principles of Ultrasound Imaging, and the mathematics behind sonar velocity, impedance, and reflection. I also understand that an Ultrasound device recieves an echo produced by impedance (~ density) changes and can calculate the distance of the echo source and it's density.

Recap: Ultrasoinc wave travels along medium, encounters section of different density/impedance andis partially reflected. The echo is picked up by the Piezoelectric elements, and, depending on Amplitude and time offset, depth and density are estimated.

So, if a wave travels from a less dense to a denser material, everything is clear to me. But if a wave travels from a denser to a less dense material, AFAIK the same echo response is sent back (Since the degree of reflection is quadratic to the fraction of Impedances ). How can a negative density gradient be correctly interpreted?

I'm really interested in this and I'd be really thankful for an explanation or a pointer to literature that explains this. I did search, and I was unable to locate a source that does. I understand that this could have to do with the transducer focusing the waves, but I cannot tell how that should contribute. Thank you for your time!

A remark up front: I am not a specialist on ultrasound. However, I am a Geophysicist and we use more or less the same principle to image the Earth's interior, just a very different frequency range compared to your problem at hand; in exploration seismics typically 5-100Hz compared to frequencies in the Megahertz range for ultrasound imaging.

A transducer emits a pulse into the medium (body). This pulse travels with a certain speed and at each "interface" (which can either be a change of density or a change of velocity) a certain part of the incoming energy gets reflected, and the other part gets transmitted through the interface into the next layer. In seismics, these are normally rock formations; in ultrasound imaging, these are different parts of your body.

The reflection (transmission) coefficient, i.e. the percentage of energy that gets reflected (transmitted) at interfaces, is determined by the acoustic impedance contrast. For pulses hitting an interface at normal incidence, the reflection coefficient is given by

$$R = \frac{z_1-z_0}{z_1+z_0}$$

where z1 is the impedance (density times velocity) in the transmission layer while z0 is the impedance in the incoming/reflection layer. The transmission coefficient is given by

$$T = \frac{2z_0}{z_1+z_0}$$

and obviously $T=1-R$ (conservation of energy; in general, I don't consider absorption here). In ultrasound imaging, as far as I know, a constant velocity is assumed for the entire body and this velocity is roughly around 1500m/s (the speed of sound in water). With this assumption, you can then determine the density contrast of two different parts of the body from the amplitude of the reflected and recorded pulse. We can't make the same constant velocity assumption for the Earth, that's why reflection seismics is a bit more complicated and first of all we have to determine some estimate for the velocity within the Earth. If we assume a constant velocity, the reflection coefficient simply becomes

$$R = \frac{\rho_1-\rho_0}{\rho_1+\rho_0}$$

where rho is the density of the medium on the two sides of the interface.

Now, if we go from a medium with small density to a medium with larger density, rho1 > rho0 and consequently the reflection coefficient R will be greater than zero. On the other hand, if we go from a medium with large density to a medium with a smaller density, rho1 < rho0. This means, the numerator in above formula will be negative, and the reflection coefficient will be negative. What does that mean?

It means, if the transducer emits a positive pulse (i.e. a waveform that goes from a zero amplitude to a positive amplitude first), we will upon return receive another positive pulse if the reflection took place at an interface where the density became larger. If the density became smaller at the interface, we will actually receive a negative pulse back, i.e. a waveform that goes from a zero amplitude to a negative amplitude first). If we display, for instance, pulses with a positive amplitude in black and pulses with a negative amplitude in white on a grey background, we can distinguish from what kind of density contrast the reflection actually came. Here is an example from the world of seismics: on the left is the impedance contrast, then comes the reflection coefficient (as you can see, positive at the first interface, negative at the second interface), then comes the wavelet (the pulse emitted by the seismic source), and then the seismic trace that we record (positive amplitudes coloured in blue, negative amplitudes coloured in red).

I hope this answers your question. The explanation is somewhat simplified, there are of course lots of pitfalls and other things to consider in practice (for instance, waves might not hit interfaces at normal incidence - then the behaviour of the reflection and transmission coefficients become much more complicated; there will be frequency-dependent absorption, etc, etc).