Apparatus for and Method of Remotely Sensing Sub-Surface Water Temperatures

Ser. No. 064,375 filed June 22, 1987, "METHOD OF REMOTELY DETECTING SUBMARINES USING A LASER."

1 Field of the Invention

The present invention relates to the measurement of sub-surface ocean temperature profiles and more particularly to an improved method of and apparatus for remotely measuring such temperature profiles from an aircraft or the like.

2. Description of the Prior Art

Profiles of ocean water temperature as a function of depth are important basic scientific data used by the oceanographic community for several purposes including determination of the effect on climate. A technique employed in the past for measuring such profiles involves the use of many temperature sensing elements spaced on a cable and towed by a ship. This is costly, time consuming and generally unsuited to high spatial and temporal resolution coverage of large ocean areas. Furthermore, the technique cannot be used for mapping ocean temperature profiles, that is, temperatures over wide areas of water, because the measuring time is too long compared to the time over which sea temperatures vary.

Another technique is described in Patent No. 4,123,160 in which a laser beam is used to illuminate the water and observation is made of the Raman scatter from the monomer and hydrogen bonded polymeric forms of water, the ratio of which is a function of temperature. This technique is vulnerable to interference from high background illumination, such as sunlight, because of the relatively wide optical bandwidth of the Raman scattering. Furthermore, in this technique there is differential absorption over the Raman band as light transits the water column. Depolarization effects of the water column also limit the effectiveness of the technique when polarization spectroscopy is employed.

Still another laser remote sensing method has been used in limited experiments, see "SPEED OF SOUND AND TEMPERATURE IN THE OCEAN BY BRILLOUIN SCATTERING" by Hirschberg, et al., Applied Optics, Aug. 1984, pages 2624-2628, inclusive. This method relies on the wavelength shift associated with Brillouin scattering from the water. This shift, however, is small so that extremely high resolution is required in optical measurement of the wavelength shift. Typically a Fabry-Perot interferometer is used to resolve the Brillouin shift. However, an interferometer requires a well collimated light source which generally is incompatible with remote sensing applications, that is, because of spreading, light must be collected from a much larger field of view than is possible with an interferometer.

This invention is directed to the measurement of sub-surface ocean temperatures while avoiding these disadvantages.

A general object of the invention is the provision of a technique of remotely and rapidly measuring ocean temperature profiles without interference from high background illumination, such as sunlight.

Another object is the provision of such measuring apparatus having relatively high overall efficiency.

Still another object is the provision of a method of and apparatus for remotely measuring ocean temperature profiles without the need for a precision interferometer.

A further object is the provision of such a technique and apparatus having relatively high sensitivity.

These and other objects of the invention are achieved by using the "self-heterodyne" of the wavelength-shifted Brillouin scatter with the unshifted Rayleigh scatter mainly from impurities in the water. The "self-heterodyne" action allows mixing of these signals from each volume element of the illuminated water column independently and measuring the frequency of the heterodyne signal. Since this frequency is directly related to the water temperature and since time is directly proportional to depth, the resultant time-temperature pattern is equivalent to a temperature-depth profile.

For a better understanding of the invention as well as other objects and further features thereof, reference is made to the following description which is to be read in conjunction with the accompanying drawings wherein:

FIG. 1 is a schematic view of an aircraft over a body of water and involved in the practice of the invention.

FIG. 2 is a schematic block diagram of the transmitter and receiver components comprising the apparatus embodying the invention.

FIG. 3 is a diagram illustrating frequency products resulting from mixing Rayleigh and Brillouin backscatter components in accordance with the practice of the invention.

FIG. 4 is a set of curves depicting the relation of ocean temperature resolution and depth for different monitoring altitudes.

It should be noted that the same reference numerals identify identical components in each of the figures.

For a better understanding of the subject invention, reference is made to the following description and appended claims in conjunction with the above-described drawings.

Referring now to the drawings, FIG. 1 depicts an aircraft 10 in flight at a predetermined altitude over a body 12 of water such as an ocean having a transmitter 13 for directing a laser beam 14 into the water and a receiver 15 for receiving a return signal 16 from the water. Transmitter 13 generates a pulsed output beam 14 having a wavelength selected to penetrate the water 12, and return signal 16 contains Brillouin and Rayleigh backscatter components. In accordance with this invention, receiver 15 processes the Brillouin and Rayleigh signals to produce data indicative of the temperature at various depths to provide temperature-depth information.

A more detailed block diagram of transmitter 13 and receiver 15 along with associated optics is shown in FIG. 2. Transmitter 13 comprises a laser 18 controlled by a timer 19 to produce output laser pulses 20. The output of laser 18 preferably is in the blue-green part of the spectrum (4500 to 5500 .ANG. in wavelength) and has a typical pulse repetition frequency of 100 Hz. A laser useful for this purpose is a frequency doubled Nd:YAG at a wavelength of 5300 .ANG.. Laser output 20 is directed by mirrors 21, 22 and 23 to the target area, in this application, the sea water 12 below aircraft 10.

Return signal 16 contains Brillouin and Rayleigh backscatter components resulting from the interaction of the laser pulses with the sea water. Mirrors 23, 24 and 25 direct the return light 16 to the receiver 15 comprising photodetector 26 having a cathode 26a, discriminator 27, analog-to-digital converter 28, analyzer unit 29 and a recorder 30. Phototector 26 preferably is a photomultiplier tube having a photocathode on which the return light 16 is incident. The PM tube converts the optical energy (photons) to electrons (current) while simultaneously amplifying the signal. The photocathode also performs the additional important function of mixing the Brillouin and Rayleigh components of signal 16 to provide a difference frequency, called the optical heterodyne frequency. The Brillouin backscatter component is frequency-shifted by the acoustic properties of the water whereas the Rayleigh backscatter component has a constant frequency and is analogous to the "local oscillator" in a superheterodyne receiver. The acoustic velocity of the water is a function of water temperature.

Referring now to FIG. 3, the Brillouin and Rayleigh components incident on cathode 26a of photodetector 26 are represented by curves 37; and 33 at frequencies f.sub.1 and f.sub.3, respectively, and by curve 34 at frequency f.sub.2. The signals at frequencies f.sub.1 and f.sub.3 are the frequency shifted Brillouin components, and the signal at f.sub.2 is the frequency unshifted Rayleigh component, the latter functioning as a "local oscillator" in an analogy to a superheterodyne receiver. The two difference frequencies between f.sub.3 and f.sub.2 and between f.sub.2 and f.sub.1 are the same frequency, f.sub.B and are additive in the baseband so that the process results in increased sensitivity. Moreover, the phase relationship of the Brillouin and Rayleigh components at photodetector 26 is correct for efficient heterodyning since both components are received from the same angle of view however wide and are produced by the same laser pulse stream.

The output (current) of photodetector 26 passes to discriminator 27, such as a delay line discriminator, which produces a voltage proportional to the baseband frequency f.sub.B. Converter 28 changes the output of discriminator 27 into digital form for analysis in analyzer 29. Analyzer 29 receives the digital frequency f.sub.B and, using the equations below, translates it into a signal or number representative of the water temperatures. By way of example, analyzer 29 may comprise a computer which compares the digital frequency input with a reference table to derive the equivalent temperature. The output of analyzer 29 is stored in recorder 30.

The accuracy of temperature measurement by analysis of Brillouin backscatter may be shown mathematically. The relation of signal-to-noise ratio (SNR) to the standard deviation of frequency measurement accuracy for a radar is given by the following expression. ##EQU1##

This standard deviation of the frequency also depends on a parameter T.sub.o which is usually taken to be the laser pulse width. In this case, however, the self broadening of the Brillouin shifted lines produces a Brillouin bandwidth of 480 MHz (See Hirschberg, J. G., et al, "SPEED OF SOUND AND TEMPERATURE IN THE OCEAN BY BRILLOUIN SCATTERING," Applied Optics, 23, 2624 (1984)). This requires that an effective T.sub.o of 2.1 ns be used in equation (1) rather than the laser pulse duration.

The Brillouin backscatter frequency f.sub.B is a function of the acoustic velocity v.sub.s and the laser optical frequency, f.sub.o, as follows

where n is the index of refraction and c is the velocity of light in free space. Differentiating equation (2) with respect to temperature yields

The sound velocity vs is given as a function of temperature by the following expression.

which, when differentiated with respect to temperature and evaluated at 10.degree. C., gives dv.sub.s /dT=4.82 m/s-.degree. C. Substituting this value into equation (3) and using 6.54.times.10.sup.14 Hz as the optical frequency for an assumed transmitter wavelength of 459 nm, the result obtained is

The temperature accuracy as a function of SNR can be obtained by dividing equation (1) by df.sub.b /dT. The resultant temperature accuracy expressed as a function of depth for a representative system is shown in FIG. 4.

While the invention has been described with reference to its preferred embodiments, it will be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the true spirit and scope of the invention. In addition, many modifications may be made to adapt a particular situation or material to the teaching of the invention without departing from its essential teachings.