Changes in the echoes from ultrasonic contrast agents with imaging parameters

Superharmonic imaging with chirp coded excitation: filtering spectrally overlapped harmonics

Superharmonic imaging improves the spatial resolution by using the higher order harmonics generated in tissue. The superharmonic component is formed by combining the third, fourth, and fifth harmonics, which have low energy content and therefore poor SNR. This study uses coded excitation to increase the excitation energy. The SNR improvement is achieved on the receiver side by performing pulse compression with harmonic matched filters. The use of coded signals also introduces new filtering capabilities that are not possible with pulsed excitation. This is especially important when using wideband signals. For narrowband signals, the spectral boundaries of the harmonics are clearly separated and thus easy to filter; however, the available imaging bandwidth is underused. Wideband excitation is preferable for harmonic imaging applications to preserve axial resolution, but it generates spectrally overlapping harmonics that are not possible to filter in time and frequency domains. After pulse compression, this overlap increases the range side lobes, which appear as imaging artifacts and reduce the Bmode image quality. In this study, the isolation of higher order harmonics was achieved in another domain by using the fan chirp transform (FChT). To show the effect of excitation bandwidth in superharmonic imaging, measurements were performed by using linear frequency modulated chirp excitation with varying bandwidths of 10% to 50%. Superharmonic imaging was performed on a wire phantom using a wideband chirp excitation. Results were presented with and without applying the FChT filtering technique by comparing the spatial resolution and side lobe levels. Wideband excitation signals achieved a better resolution as expected, however range side lobes as high as -23 dB were observed for the superharmonic component of chirp excitation with 50% fractional bandwidth. The proposed filtering technique achieved >50 dB range side lobe suppression and improved the image quality without affecting the axial resolution.

Quantitative contrast-enhanced ultrasound imaging: a review of sources of variability

Ultrasound provides a valuable tool for medical diagnosis offering real-time imaging with excellent spatial resolution and low cost. The advent of microbubble contrast agents has provided the additional ability to obtain essential quantitative information relating to tissue vascularity, tissue perfusion and even endothelial wall function. This technique has shown great promise for diagnosis and monitoring in a wide range of clinical conditions such as cardiovascular diseases and cancer, with considerable potential benefits in terms of patient care. A key challenge of this technique, however, is the existence of significant variations in the imaging results, and the lack of understanding regarding their origin. The aim of this paper is to review the potential sources of variability in the quantification of tissue perfusion based on microbubble contrast-enhanced ultrasound images. These are divided into the following three categories: (i) factors relating to the scanner setting, which include transmission power, transmission focal depth, dynamic range, signal gain and transmission frequency, (ii) factors relating to the patient, which include body physical differences, physiological interaction of body with bubbles, propagation and attenuation through tissue, and tissue motion, and (iii) factors relating to the microbubbles, which include the type of bubbles and their stability, preparation and injection and dosage. It has been shown that the factors in all the three categories can significantly affect the imaging results and contribute to the variations observed. How these factors influence quantitative imaging is explained and possible methods for reducing such variations are discussed.

A consistent tissue attenuation coefficient estimator using bubble harmonic echoes

The ultrasonic property of soft tissue can be quantified by its attenuation coefficient α. Traditionally the backscattering signal of tissue is used to estimate α. To improve precision, a large number of spatially independent samples of tissue echoes are required for averaging. In this paper, we propose a new estimation method, which makes use of microbubbles to provide temporally independent samples for averaging. It is easier for temporal sampling to maintain ergodicity and provide a large number of independent samples for statistical averaging. A stochastic model for the harmonic signals of an ideal bubble attenuated by tissue is derived based on Kuc's and Miller's works. An estimator of α is then presented. This estimator is consistent and could be biased because of the unknown squarelaw relation between the second and fundamental harmonics for non-ideal bubble oscillation. In experimental works, we design a simplified phantom for demonstrating the performance of the proposed estimator. It is shown that both first and second harmonics can estimate α consistently. However, the interference of the tissue backscattering signal may cause additional estimation error using the first harmonic.

Direct observations of ultrasound microbubble contrast agent interaction with the microvessel wall

Many thousands of contrast ultrasound studies have been conducted in clinics around the world. In addition, the microbubbles employed in these examinations are being widely investigated to deliver drugs and genes. Here, for the first time, the oscillation of these microbubbles in small vessels is directly observed and shown to be substantially different than that predicted by previous models and imaged within large fluid volumes. Using pulsed ultrasound with a center frequency of 1 MHz and peak rarefactional pressure of 0.8 or 2.0 MPa, microbubble expansion was significantly reduced when microbubbles were constrained within small vessels in the rat cecum (p<0.05). A model for microbubble oscillation within compliant vessels is presented that accurately predicts oscillation and vessel displacement within small vessels. As a result of the decreased oscillation in small vessels, a large resting microbubble diameter resulting from agent fusion or a high mechanical index was required to bring the agent shell into contact with the endothelium. Also, contact with the endothelium was observed during asymmetrical collapse, not during expansion. These results will be used to improve the design of drug delivery techniques using microbubbles.

Radiation-force assisted targeting facilitates ultrasonic molecular imaging

Ultrasonic molecular imaging employs contrast agents, such as microbubbles, nanoparticles, or liposomes, coated with ligands specific for receptors expressed on cells at sites of angiogenesis, inflammation, or thrombus. Concentration of these highly echogenic contrast agents at a target site enhances the ultrasound signal received from that site, promoting ultrasonic detection and analysis of disease states. In this article, we show that acoustic radiation force can be used to displace targeted contrast agents to a vessel wall, greatly increasing the number of agents binding to available surface receptors. We provide a theoretical evaluation of the magnitude of acoustic radiation force and show that it is possible to displace micron-sized agents physiologically relevant distances. Following this, we show in a series of experiments that acoustic radiation force can enhance the binding of targeted agents: The number of biotinylated microbubbles adherent to a synthetic vessel coated with avidin increases as much as 20-fold when acoustic radiation force is applied; the adhesion of contrast agents targeted to alpha(v)beta3 expressed on human umbilical vein endothelial cells increases 27-fold within a mimetic vessel when radiation force is applied; and finally, the image signal-to-noise ratio in a phantom vessel increases up to 25 dB using a combination of radiation force and a targeted contrast agent, over use of a targeted contrast agent alone.

Ultrasonic analysis of peptide- and antibody-targeted microbubble contrast agents for molecular imaging of alphavbeta3-expressing cells

The goal of targeted ultrasound contrast agents is to significantly and selectively enhance the detection of a targeted vascular site. In this manuscript, three distinct contrast agents targeted to the alphavbeta3 integrin are examined. The alphavbeta3 integrin has been shown to be highly expressed on metastatic tumors and endothelial cells during neovascularization, and its expression has been shown to correlate with tumor grade. Specific adhesion of these contrast agents to alphavbeta3-expressing cell monolayers is demonstrated in vitro, and compared with that of nontargeted agents. Acoustic studies illustrate a backscatter amplitude increase from monolayers exposed to the targeted contrast agents of up to 13-fold (22 dB) relative to enhancement due to control bubbles. A linear dependence between the echo amplitude and bubble concentration was observed for bound agents. The decorrelation of the echo from adherent targeted agents is observed over successive pulses as a function of acoustic pressure and bubble density. Frequency-domain analysis demonstrates that adherent targeted bubbles exhibit high-amplitude narrowband echo components, in contrast to the primarily wideband response from free microbubbles. Results suggest that adherent targeted contrast agents are differentiable from free-floating microbubbles, that targeted contrast agents provide higher sensitivity in the detection of angiogenesis, and that conventional ultrasound imaging techniques such as signal subtraction or decorrelation detection can be used to detect integrin-expressing vasculature with sufficient signal-to-noise.

Radiation-Force Assisted Targeting Facilitates Ultrasonic Molecular Imaging

Ultrasonic molecular imaging employs contrast agents, such as microbubbles, nanoparticles, or liposomes, coated with ligands specific for receptors expressed on cells at sites of angiogenesis, inflammation, or thrombus. Concentration of these highly echogenic contrast agents at a target site enhances the ultrasound signal received from that site, promoting ultrasonic detection and analysis of disease states. In this article, we show that acoustic radiation force can be used to displace targeted contrast agents to a vessel wall, greatly increasing the number of agents binding to available surface receptors. We provide a theoretical evaluation of the magnitude of acoustic radiation force and show that it is possible to displace micron-sized agents physiologically relevant distances. Following this, we show in a series of experiments that acoustic radiation force can enhance the binding of targeted agents: The number of biotinylated microbubbles adherent to a synthetic vessel coated with avidin increases as much as 20-fold when acoustic radiation force is applied; the adhesion of contrast agents targeted to alpha(v)beta3 expressed on human umbilical vein endothelial cells increases 27-fold within a mimetic vessel when radiation force is applied; and finally, the image signal-to-noise ratio in a phantom vessel increases up to 25 dB using a combination of radiation force and a targeted contrast agent, over use of a targeted contrast agent alone.

Ultrasonic Analysis of Peptide- and Antibody-Targeted Microbubble Contrast Agents for Molecular Imaging of α v β 3 -Expressing Cells

The goal of targeted ultrasound contrast agents is to significantly and selectively enhance the detection of a targeted vascular site. In this manuscript, three distinct contrast agents targeted to the alphavbeta3 integrin are examined. The alphavbeta3 integrin has been shown to be highly expressed on metastatic tumors and endothelial cells during neovascularization, and its expression has been shown to correlate with tumor grade. Specific adhesion of these contrast agents to alphavbeta3-expressing cell monolayers is demonstrated in vitro, and compared with that of nontargeted agents. Acoustic studies illustrate a backscatter amplitude increase from monolayers exposed to the targeted contrast agents of up to 13-fold (22 dB) relative to enhancement due to control bubbles. A linear dependence between the echo amplitude and bubble concentration was observed for bound agents. The decorrelation of the echo from adherent targeted agents is observed over successive pulses as a function of acoustic pressure and bubble density. Frequency-domain analysis demonstrates that adherent targeted bubbles exhibit high-amplitude narrowband echo components, in contrast to the primarily wideband response from free microbubbles. Results suggest that adherent targeted contrast agents are differentiable from free-floating microbubbles, that targeted contrast agents provide higher sensitivity in the detection of angiogenesis, and that conventional ultrasound imaging techniques such as signal subtraction or decorrelation detection can be used to detect integrin-expressing vasculature with sufficient signal-to-noise.

Anin vitrostudy of a microbubble contrast agent using a clinical ultrasound imaging system

Optimal insonation settings for contrast imaging are yet to be specified, mainly due to the lack of good understanding of the behaviour of the microbubbles. A satisfactory model that explains the behaviour of individual contrast agent scatterers has not yet been reported in the literature. An in vitro system based on a commercial scanner (ATL HDI3000) has been developed to investigate the backscatter of such agents. Suspensions of Definity were introduced in an anechoic tank. The frequency of transmitted ultrasound varied from 1 to 5 MHz, pulse period from 2 to 10 periods and peak negative acoustic pressure from 0.08 to 1.7 MPa. The backscatter at the fundamental and second harmonic frequency windows from the agent was normalized in terms of the corresponding components of backscatter from a blood mimicking fluid suspension. The agent provided a dominant resonance effect at 1.6 MHz transmit frequency. Second harmonic normalized backscatter averaged around 9 dB higher than the fundamental. The normalized fundamental backscatter intensity was linear with peak negative pressure. The second harmonic at resonance peaked at 0.5 MPa suggestive of bubble disruption above such pressure. The system proved capable of illustrating the ultrasonic behaviour of Definity in vitro, and the investigation suggested particular insonation conditions for optimal image enhancement using Definity.

Volume scattering of distributed microbubbles and its influence on blood flow estimation

In recent years, microbubble contrast agents have become a potential adjunct in Doppler ultrasound diagnosis. In this paper, we show that volume scattering makes the effective band in Doppler spectrum shift downward after injection of microbubbles. Because the insonified volume comprises a collection of distributed microbubbles, the statistical properties such as the autocorrelation function and ensemble average power spectrum of the echoes from a collection of distributed microbubbles were derived first. It can be observed that, beyond a critical frequency, the theoretical volume backscattering cross section derived from the ensemble average power spectrum of microbubbles decreases with frequency. On the contrary, the volume backscattering cross section of red cells increases with frequency. Using two-dimensional (2-D) Fourier transform, the variation in Doppler spectrum caused by different volume backscattering cross section can be demonstrated, and the consequential downward shifts of the estimated Doppler parameters (e.g., the mean and maximum Doppler shifts, and the variance of Doppler power spectrum) after microbubble injection are shown. In addition, it can be observed that the variation gets larger as the transmitted bandwidth increases. And, the variations in Doppler parameters estimated with experimental data are presented to verify the theoretical deviations.

The ultrasonic weak short-pulse responses of microbubbles based on a two-frequency approximation

The ultrasonic short-pulse responses of microbubbles are of interest in cavitation, transient responses, and contrast imaging. We extend the two-frequency analytic solutions of Newhouse and Shankar [J. Acoust. Soc. Am. 75, 1473-1477 (1984)] to approximate the short-pulse responses of microbubbles in a low-amplitude field. Based on their results, there is an expected component near dc in the spectrum of bubble echoes excited by a short pulse. Here this component is named the low-frequency response, and its theoretical properties are verified experimentally. Including the fundamental and second-harmonic components, the weak short-pulse responses of microbubbles include three types of response. Our work has determined the constraint conditions under which this approximated solution can be used to analyze these short-pulse responses. This paper also provides the amplitude and spectral properties of these responses. The low-frequency response has a special bandwidth-dependent property and has potential applications in imaging and bubble sizing.

The behaviour of individual contrast agent microbubbles

In recent years, our knowledge of the behaviour of ultrasonic microbubble contrast agents has improved substantially through in vitro experiments. However, there has been a tendency to use high concentrations of contrast agents in suspension, so that ultrasonic backscatter data are generated by a cloud of microbubbles. Such experiments involve a variety of assumptions with validity that is open to question. In addition, high concentrations of microbubbles cannot be used to understand the behaviour of individual microbubble scatterers. This paper proposes a technique that minimises the number of assumptions that need to be made to interpret in vitro experimental data. The basis of the technique is a dilute suspension of microbubbles that makes single scattering events distinguishable. A commercial scanner was used to collect radio frequency (RF) data from suspensions of two different contrast agents, Quantison and Definity. Backscatter data were collected over a range of acoustic pressures. It was found that Definity provided a constant number of scattering events per unit volume of suspension for almost all applied acoustic pressures. Quantison demonstrated an increasing number of scattering events per unit volume with increasing acoustic pressure. Below 0.6 MPa, Quantison scatterers were not individually detectable and provided levels of backscatter similar to those of a blood-mimicking fluid, which suggests that Quantison microbubbles had almost linear scattering behaviour. At acoustic pressures greater than 0.6 MPa, both agents appeared to provide echoes from free bubbles. The change in the number of scatterers per unit volume with acoustic pressure cannot be demonstrated using high concentrations of contrast agent.

Understanding the limitations of ultrasonic backscatter measurements from microbubble populations

Despite over ten years of in vitro investigations of ultrasound contrast agents, the level of understanding of their behaviour in ultrasound fields is limited. Several problems associated with these investigations, particular to the nature of contrast agents, are discussed. Using a commercial scanner the RF normalized backscatter of two different contrast agents (Definity and Quantison) was measured at different suspension concentrations and acoustic pressures. Both contrast agents scattered ultrasound nonlinearly and the backscatter showed a dependence on acoustic pressure. In order to assess the average behaviour of the agents across the range of acoustic pressures and microbubble concentrations the experimental data were fitted to a theoretically acceptable model using nonlinear regression analysis. The analysis showed that both the backscatter and the attenuation of the Quantison suspensions displayed a higher order of dependence on acoustic pressure than the Definity suspensions. It was also discovered that Quantison microbubbles did not demonstrate uniform behaviour across the acoustic pressure range. At lower acoustic pressures the behaviour could not follow a model similar to that which predicted the behaviour at higher acoustic pressures, which was mainly due to the fact that free bubbles were released in a fashion dependent on acoustic pressure. The fact that two different populations of scatterers exist in the same suspensions makes the assessment of the behaviour of the particular agent impossible with the high concentrations that are commonly used. Very low concentration suspensions whereby single scattering events can be monitored should be more useful. In conclusion, the approach of using high microbubble concentrations in order to investigate the properties of ultrasonic contrast agents is limited in that the results of such studies cannot be used to understand the behaviour of single microbubbles.

In vitro acoustic characterisation of four intravenous ultrasonic contrast agents at 30 MHz

The acoustic properties of four ultrasonic contrast agents (Optison, Definity, SonoVue and Sonazoid) were studied at 30 MHz using a Boston Scientific ClearView Ultra intravascular ultrasound (US) scanner modified to allow access to the unprocessed US data. A range of contrast agent concentrations were studied using either saline or glucose as the diluent of choice. Mean backscatter power was measured over regions-of-interest (ROI) at distances of 1, 1.5, 2, 3, 4 and 5 mm from the centre of the intravascular probe and normalised to the US data collected from a standard glass reflector. For all of the agents, the mean backscatter power at 30 MHz varied in a linear manner with concentration between 0.01 million microbubbles/mL and 1 million microbubbles/mL. Furthermore, for two of the agents, mean backscatter enhancement was detectable at concentrations as low as 2 microbubbles/sample volume.

A review of physical phenomena associated with ultrasonic contrast agents and illustrative clinical applications

Successful clinical applications of contrast agents involve an understanding of the physical interaction of ultrasound (US) with contrast agents. This paper reviews the physical phenomena involved in these interactions and discusses the relevant theoretical background for modeling US-contrast agent interactions. Measurement techniques using US to obtain information regarding contrast agents are summarized. Illustrative clinical applications are given in the second part of the paper. Recent developments in nonlinear imaging techniques and transient techniques are reviewed. New methods, such as depletion perfusion measurement, and high-frequency applications are included.

An in vitro system for the study of ultrasound contrast agents using a commercial imaging system

An in vitro system for the investigation of the behaviour of contrast microbubbles in an ultrasound field, that provides a full diagnostic range of settings, is yet to be presented in the literature. The evaluation of a good compromise of such a system is presented in this paper. It is based on (a) an HD13000 ATL scanner (Bothell, WA, USA) externally controlled by a PC and (b) on the use of well-defined reference materials. The suspensions of the reference ultrasonic scattering material are placed in an anechoic tank. The pulse length ranges from 2 to 10 cycles, the acoustic pressure from 0.08 to 1.8 MPa, the transmit frequency from 1 to 4.3 MHz, and the receive frequency from 1 to 8 MHz. The collection of 256 samples of RF data, at an offset distance from the transducer face, was performed at 20 MHz digitization rate, which corresponds to approximately 1 cm depth in water. Two particle suspensions are also presented for use as reference scatterers for contrast studies: (a) a suspension of Orgasol (ELF Atochem, Paris, France) particles (approximately 5 microm mean diameter) and (b) a suspension of Eccosphere (New Metals & Chemicals Ltd, Essex, UK) particles (approximately 50 microm mean diameter). A preliminary experiment with the contrast agent Definity (DuPont Pharmaceutical Co, Waltham, MA) showed that the above two materials are suitable for use as a reference for contrast backscatter.

Techniques for perfusion imaging with microbubble contrast agents

The acoustic properties of ultrasound contrast agents vary widely with agent composition and insonation conditions. For contrast imaging, methods are required to match RF and Doppler processing to each combination of transmission parameters and agent and tissue properties. We propose a method that uses the measured or modeled echoes from agent and tissue to specify directly the characteristics of RF and Doppler filters for contrast imaging. The proposed method is sufficiently general to cover most common imaging techniques including harmonic greyscale, Doppler, and pulse inversion imaging. Using this method, sample filters were designed to detect myocardial perfusion with the contrast agent Optison (Mallinckrodt Medical, St. Louis, MO) under selected imaging conditions. Simplified power Doppler filtering, using a weighted sum of the Doppler samples, matched the performance of more complicated matrix filters. By coordinating the selection of RF and Doppler filters rather than designing these filters sequentially, agent-to-tissue contrast was increased by up to 3.9 dB. Under some conditions, fundamental RF filtering outperformed harmonic filtering for intermittent Doppler imaging.

Contrast agent stability: a continuous B-mode imaging approach

The stability of contrast agents in suspensions with various dissolved gas levels has not been reported in the literature. An in vitro investigation has been carried out that studied the combined effect of varying the acoustic pressure along with degassing the suspension environment. In this study, the contrast agents were introduced into suspensions with different oxygen concentration levels, and their relative performance was assessed in terms of decay rate of their backscatter echoes. The partial pressures of oxygen in those solutions ranged between 1.5 and 26 kPa. Two IV and one arterial contrast agents were used: Definity, Quantison, and Myomap. It was found that Quantison and Myomap released free bubbles at high acoustic pressure that also dissolved faster in degassed suspensions. The backscatter decay for Definity did not depend on the air content of the suspensions. The destruction of bubbles was dependent on acoustic pressure. Different backscatter performance was observed by different populations of bubbles of the last two agents. The physical quantity of "overall backscatter" (OB) was defined as the integral of the decay rate over time of the backscatter of the contrast suspensions, and improved significantly the understanding of the behaviour of the agents. A quantitative analysis of the backscatter properties of contrast agents using a continuous imaging approach was difficult to achieve. This is due to the fact that the backscatter in the field of view is representative of a bubble population affected by the ultrasound (US) field, but this bubble population is not representative of the contrast suspension in the whole tank. Single frame insonation is suggested to avoid the effects of decay due to the ultrasonic field, and to measure a tank-representative backscatter. The definition of OB was useful, however, in understanding the behaviour of the agents.

Parametric study of the peak negative acoustic pressure distribution within the image plane of a phased array transducer

A study is presented in which the influence of pulse and transducer characteristics on the distribution of the peak negative acoustic pressure within the image plane of a phased array transducer is shown. First, the influence of the center frequency, the bandwidth of the transmitted pulse, the position of the focus point, and the properties of the phased array are investigated by computer simulation. These simulations show that both the maximal amplitude and the homogeneity of the distribution of the peak negative pressure within a sector image depend on these parameters. Most remarkably, they show that the maximal peak pressure does not necessarily occur on the geometrical symmetry axis of the array transducer but can occur at the outermost lateral lines of a sector image. Second, this effect is demonstrated experimentally by hydrophone measurements of the sound field produced by a standard, clinically used, 2.5-MHz phased array transducer. The experimental results agree well with the theoretical model. The relevance of this study for contrast echography and the safety measurements of ultrasound is discussed.

Optical imaging of contrast agent microbubbles in an ultrasound field with a 100-MHz camera

Ultrasound (US) contrast agents, used in the field of medical diagnosis, contain small microbubbles of a mean diameter of about 3 microm. The acoustic behavior of these bubbles in US field has been subject to many investigations. In this study, we propose a method to visualize the behavior of the bubbles in a 0.5-MHz US field under a microscope with a frame rate of 4 MHz. For low acoustic pressures (peak negative pressure of 0.12 MPa), the radius-time curve as measured from the optical images is in agreement with the theory. For higher acoustic pressures (peak negative pressure of 0.6 MPa), the recorded radius is significantly larger than predicted by theory and sudden change in the bubbles shapes has been noticed. The proposed method enables the study and characterization of individual bubbles and their encapsulation. It is expected that this will open new areas for quality control, US contrast imaging and US-guided drug delivery.

Experimental and theoretical evaluation of microbubble behavior: effect of transmitted phase and bubble size

Ultrasound contrast agents provide new opportunities to image vascular volume and flow rate directly. To accomplish this goal, new pulse sequences can be developed to detect specifically the presence of a microbubble or group of microbubbles. We consider a new scheme to detect the presence of contrast agents in the body by examining the effect of transmitted phase on the received echoes from single bubbles. In this study, three tools are uniquely combined to aid in the understanding of the effects of transmission parameters and bubble radius on the received echo. These tools allow for optical measurement of radial oscillations of single bubbles during insonation, acoustical study of echoes from single contrast agent bubbles, and the comparison of these experimental observations with theoretical predictions. A modified Herring equation with shell terms is solved for the time-dependent bubble radius and wall velocity, and these outputs are used to formulate the predicted echo from a single encapsulated bubble. The model is validated by direct comparison of the predicted radial oscillations with those measured optically. The transient bubble response is evaluated with a transducer excitation consisting of one-cycle pulses with a center frequency of 2.4-MHz. The experimental and theoretical results are in good agreement and predict that the transmission of two pulses with opposite polarity will yield similar time domain echoes with the first significant portion of the echo generated when the rarefactional half-cycle reaches the bubble.

Ultrasonic characterization of the nonlinear properties of contrast microbubbles

The nonlinear properties of microbubble contrast agents have been used to create contrast-specific imaging modalities such as harmonic imaging and subharmonic imaging. Thus, a better understanding of the nonlinear performance of contrast microbubbles may enhance the diagnostic capabilities of medical ultrasound (US) imaging. The first and second harmonic, the 1/2 order subharmonic and the 3/2 order ultraharmonic components in spectra of scattered signals from Optison microbubbles insonified at 2 and 4 MHz have been investigated using an in vitro laboratory pulse-echo system. The development of these signal components over time is quite different for 2-MHz insonification compared to 4-MHz insonification. Scattered subharmonic and ultraharmonic signals are much more time-dependent than first and second harmonic echoes. The dependence of the first and second harmonic, subharmonic and ultraharmonic components on acoustic pressure for 2-MHz insonification is similar to that for 4-MHz insonification. The first and second harmonic components increase linearly with acoustic pressure (in double logarithmic scales) and the subharmonic and ultraharmonic amplitudes undergo rapid growths in the intermediate acoustic pressure range and much slower increases at both lower and higher acoustic pressures.