METHOD AND SYSTEM FOR THE DETECTION AND MEASUREMENT OF AIR PHENOMENA AND TRANSMITTER AND RECEIVER FOR USE IN THE SYSTEM

Abstract: A system and method for the detection and measurement of atmospheric air movement irregularities such as wind velocity vector, clear air turbulence, aircraft induced vorteces and turbulence, by means of electromagnetic waves. An air volume (10) under investigation is illuminated by a transmitter (1, 1A, 1B) with a beam of coherent electromagnetic energy (1E), and a resulting wave field (2E) is received and subsequently coherently demodulated and processed in process means to derive information on the existence of the said atmospheric irregularities and furthermore to give specific measurements of related parameters. At least one receiver (2, 2A, 2B) for said resulting wave field (2E), which is due to scattering in said air volume (10), is positioned at a bistatic location having a selected distance (2R) from said transmitter (1, 1A, 1B). Said electromagnetic energy (1E) is transmitted at microwave frequencies in a concentrated beam. Transmitter and receiver antennas (1A, 1B, 2A, 2B), respectively, are directed so as to cover said air volume (10) with a scattering angle ( THETA ) chosen to have such a (low) value that scattering from atmospheric irregularities in said air volume (10) has a higher power than noise sources in the system.

Description: Method and system for the detection and measurement of air phenomena and transmitter and receiver for use in the system. This invention relates to a system for the detection and measurement of velocity, turbulence, vorteces and similar irregularities or phenomena in air, including classification of such phenomena. These may comprise the wind velocity vector, clear air turbulence as well as aircraft induced vorteces and turbulence. Detection and measurement as contemplated here takes place by means of electromagnetic waves. The invention also comprises transmitter and receiver equipment for use in the system, as well as an associated method. The invention is primarily directed to radar-based technology for the measurement of wind, wind shear, aircraft wave vorteces and turbulence affecting air traffic safety. Thus, a simple and inexpensive bistatic radar is described theoretically and with practical examples. In general various and different fields of use of this invention will be possible, such as for meteorological purposes. A general background of interest in this connection can be found in the following publications: REFERENCES 1) Gjessing, Dag T, (1962) : Determination of permittivity variations in the tropo¬ sphere by scatter-propagation methods. IEE Monograph No 510E, April. 2) Gjessing, Dag T, (1969) : Atmospheric structure deduced from forward-scatter wave propagation experiments. Radio Science, Vol 4, No 12, pp 1195-1201, December. 3) Gjessing, Dag T, Anton G. Kjelaas and J. Nordø, (1969) : Spectral measurements and atmospheric stability. J of the Atmospheric Sciences, Vol 26, February. 4) Gjessing, Dag T, (1986) : Target adaptive matched illumination RADAR: Principles & applications. Peter Peregrinus Ltd, London. 5) Gjessing, Dag T and Jens Hjelmstad, (1989) : Artificial Perception of objects and scattering surfaces based on electromagnetic waves. Target Adaptive Matched illumination radar (TAMIR) . Radar 1989. Reference to publications listed above will be made in the following description. More particularly reference is made at this point to published European patent application No. 436.048A1, which is directed to a method and system for measuring atmospheric wind fields, i.e. to problems related to those with which the present invention is concerned. Thus, according to EP 436.048A1 an air volume under investigation is illuminated by a transmitter with a beam of coherent electromagnetic energy and a resulting, i.e. reflected wave field is received and subsequently coherently demodulated and processed to derive information on the atmospheric phenomena of interest. It is to be noted that the known method and system just referred to, is based on resulting waves due to coherent back-scattering (reflection) of transmitted VHF or UHF frequencies, which means that desired performance with respect to accuracy, resolution, flexibility and degrees of freedom in available functions and measurements, is not attained. Substantial improvements in these and other respects are obtained in the system according to the present invention, which basically uses microwave frequencies and measures coherent forward scatter waves. The novel and specific features according to this invention are set out more completely in the claims. The invention as well as resulting advantages will be explained more in detail in the following description, referring also to the drawings in which: Fig. 1 schematically illustrates the overall system according to the invention, with exemplary embodi- ments of transmitter and " receiver arrangements, Fig. * 2 shows examples of wavenumber spectra of interest in connection with wind-shear phenomena, Fig. 3 is a table of relationships characterizing the scattering from wind-shear phenomena, Fig. 4 is a diagram with curves illustrating the theore¬ tical relationship between the measured quantity and the slope of the refractive-index spectrum, Fig. 5 illustrates an interference pattern resulting from transmission of two closely spaced frequencies from two different apertures, Fig. 6 schematically shows an interferometric radar system with two apertures, Fig. 7 is a schematic illustration of the general multi- frequency bistatic radar system according to the invention, with indications of relevant variables in the system, Fig. 8 is a diagram showing as an example the effect of using 16 antennas and 16 frequencies in a common, monostatic radar system, Fig. 9 shows a three-dimensional diagram related to the example of Fig. 8, Fig. 10 is a diagram of turbulence spectra, and Fig. 11 illustrates an example of Doppler broadening. Introduction, Presentation of Problem, General Philosophy and Measuring System. Remote measurement of the wind field within a limited volume on the glidepath within which the wind field has major impact on the aircraft, is the primary topic of this descrip- tion. Several weather conditions are of importance. In this connection emphasis will be placed on four conditions. a) Little or no wind in the lowest part of the boundary layer monitored by in situ ground based instruments. Strong wind above certain height. The two atmospheric layers may also have marked differences in temperature. b) Local areas of strong down-and/or updrafts resulting from lee-wave occurrences or local convective plume phenomena. c) Strong anisotropic turbulence characterized as T 2 ( L1 ) κ-~. v d) Small scale aircraft wake vortex. It is the objective of this description to present a simple inexpensive system that will constitute a diagnostic tool based on remote assessments of air dynamics. The under¬ lying philosophy in regard to such a measuring system is the following: Although the ideal system probably will consist of a single radar covering a large geographical area capable of giving the wind vectors in a large 3-directional space, the air safety requirements may well be satisfied with a simple dedicated system which by remote sensing methods gives information about a set of parameters that is strongly related to, or is a direct consequence of, a given wind field adverse to air safety. - A multi-sensor-data-fusion system with a maximum of orthogonality (maximum degrees of freedom) should be employed. Learning with time, the system redundancy can be reduced thus reducing system cost. Potential Remote Sensed Parameters Describing Aircraft Wake Vortex, Wind Shear, Down Draft and Turbulence Conditions. When discussing such parameters pragmatically, the fundamentals of the observational system should first be closely defined. (See Fig. 1 to be explained below) . Secondly, one shall need a clear understanding regarding the prevailing atmospheric conditions. Since the degrees to which these conditions affect air traffic safety is dramati¬ cally different, it is important to design a system with a suff cient number of degrees of freedom so as to be able to make the appropriate classification of windshear phenomena and other athmospheric irregularities. Specifically, if the measurements are limited to measuring scattered power and Doppler shift which is the case with most radar concepts, the false alarm rate is expected to be considerable. As an example, a bird flying at a given speed through the scattering volume may well give rise to a scattered power and a Doppler shift which is similar to that caused by an aircraft wake vortex. The basic principles of the system according to the present invention shall now be described with reference to Fig. 1, showing the general architecture of the signal analysis and pattern recognition in the system. A trans¬ mitter generally denoted 1 with an antenna or antennas 1A and IB cooperates with a receiver generally denoted 2 having an antenna or antennas 2A. Thus, two or more antennas (aper¬ tures) can be provided at the transmitter 1 and/or the receiver 2. As indicated there is a certain distance 2R between the transmitter antenna 1A and the receiver antenna 2A.. From antenna 1A there is shown a beam of coherent electromagnetic energy, the axis of which is denoted IE. Correspondingly, receiver antenna 2A is adapted to receive a wave field as indicated at 2E, i.e. within a certain angle in space around the central axis indicated from antenna 2A. Accordingly, scattering of the transmitted electromagnetic energy IE from the common air volume indicated at 10, results in the received wave field 2E, which makes it possible to investigate atmospheric air movement phenomena occurring in the common air volume 10. Thus, the system illustrated in Fig. 1 utilizes the forward scattering mechanism and in this connection it is essential that the scattering angle θ is chosen to have such a low value that scattering from homogeneous isotropic turbu¬ lence in air volume 10 has a higher power than noise sources in the system. The scattering angle θ of course is related to the distance 2R separating the transmitter antenna 1A and the receiver antenna 2A, as well as the elevation of the central axes IE and 2E respectively of these antennas. In the case of detection or measurement relating to the glide¬ path for aircraft in connection with a runway, the arrange¬ ment of antennas 1A and 2A with their beams IE and 2E respec tively, is usually oriented transversally to the longitudina direction of the glidepath. Moreover, in practice and on the basis of suitable microwave frequencies employed, the scatte¬ ring angle θ preferably is chosen to be smaller than 20°. Typical specifications in the embodiment illustrated in Fig. 1 are as follows: Distance 2R — 7 km Transmitter power P τ = 5 WATTS Frequency F = 10 GHz (λ * - *** 3 cm) Antenna aperture d = 100 cm Transmitter antenna elevation α = 2° As indicated in Fig. 1, transmitter 1 comprises two transmitter circuits generally denoted 11 and 12 respec¬ tively, for frequences Fl and Fl + ΔF. The receiver 2 has corresponding receiver channels comprising coherent homodyne demodulater circuits 21 and 22 respectively. For this demodulation, a reference frequency F^ p , is transmitted by means of a separate channel or path in which antennas 13 and 23 are incorporated. As will be understood the required reference frequency can be provided by other means than illustrated as an example in Fig. 1. A preferred solution comprises local oscillators adapted to be locked to a common reference, such as a GPS satellite. Thus, according to the above, it is preferred in the system of the invention to employ a transmitter for two or more coherently related microwave frequences. Moreover, as will be understood by persons skilled in the art, it is an advantage that the transmitter and/or receiver antennas 1A and 2A are arranged to radiate a minimum of electromagnetic power in the direction of the other antenna, thus avoiding undue disturbance or masking of the forward scatter signal which is of interest. In order to cover an increased air volume or space the receiver and/or transmitter antennas 1A and 2A should preferably be adapted to be steered mechani¬ cally or electronically. This will make it possible to obtain a two or three dimensional map of atmospheric pheno¬ mena of interest in particular windshear features. Processor means or unit 25 shown in Fig. 1 in practice will be in the form of a computer with necessary programs for performing functions, processing signals and implementing algorithms to be described further below, examples of which are Delta-K-processing, algorithms for inversion, classifi¬ cation, generating images and windshear alerts. In this connection FFT processing and Kal an filtering are usually within the capabilities of processor unit 25. Block 26 is a display or data link unit. As indicated in Fig. 1 the wind assessment system is based on a highly sensitive coherent bistatic radar system retaining the transmitting frequency on the receiver site by direct transmission of a reference signal R j ^ p or by the use of stable reference oscillators at both ends. This coherent system makes it possible to measure ampli¬ tude and phase (Doppler broadening and Doppler shift) at the receiver. Thus, an important feature according to the invention consists therein that the coherent homodyne demodulators 21 and 22 give a Doppler shift in the plane through the trans¬ mitter and receiver antennas 1A and 2A and the scattering air volume 10. In order for the instantaneous bandwidth to be measured as a function of time, many frequencies are employed, as indicated with F- ^ and (Fη+ΔF) in Fig. 1. At the receiver 2, amplitude and phase are measured at all these frequences and the complex voltages thus obtained are multiplied (correlated) , this being one of the functions performed by processor unit 25.. The degree of correlation R(ΔF) is, as explained furthe below, a direct measure of bandwidth which in turn determine the slope ή of the refractive index spectrum §(K) - K "n . Similarly, by making use of spaced receiving antennas the spatial correlation properties of the scattered field ca be measured. From this, important information regarding the dynamics of the air can be obtained. Various aspects of this invention as briefly set out above, as well as other aspects, will be further explained i the following description. Four important classes of air motion being of most interest in this connection, are as follows: The measurable parameters in relation to these four classes, are: TRANSMISSION LOSS (SCATTERING CROSS SECTION) BANDWIDTH OF TRANSMISSION PATH - DOPPLER SPECTRUM (MOTION PATTERN) SPATIAL CORRELATION OF SCATTERED FIELD FREQUENCY DEPENDANCE REFRACTION (BENDING OF RAY) Thus, according to the principles of this invention, four different air dynamic classes can be characterized by 6 radio circuit parameters or signature domains. These relationships will now briefly be sketched with reference to earlier contributions according to references 1-5 listed above. Atmospheric Structure Deduced from Forward Scatter Wave Propagation. This topic has been subjected to detailed considerations with emphasis on communications applications as early as twenty to thirty years ago. This invention is in part based on the result of these early considerations for air safety application. First, consider the turbulence spectrum φ(K) - K "n . The spectrum parameter n characterizes the distribution of "eddy size". Turbulence spectrum φ (K) : It is well established through similar considerations that the 3-dimensional velocity spectrum is of the same form as that of the refraction index irregularity spectrum. We write this on the form φ(K) - K" n where n = 5/3 is a good estimate for high values of K in a well mixed isotropic atmosphere. Close to the ground along the glidepath, strong degrees of anisotropy and gradients of wind temperature and velocity has an effect on the turbulence parameter n. Furthermore, the large scale eddies are governed by the nearness of the ground leading to a cut-off of the spectrum at wavenumbers V- 2π where H is the height (e.g. glidepath) above ground. Fig. 2 illustrates this with a sketch of a vortex spectrum and a CAT spectrum included for comparison. A clear difference is seen between the resonant character of the vortex spectrum and the characteristic slope (K -5 - 3 ) of the CAT spectrum. It is seen that in addition to measuring the scattered power, the shape of the refrative index spectrum contributes to the assessment of air motion severity. Thus, the spatial wavenumber spectrum of refractive index is goverened by many dynamic factors. The measurable quantities in the table of Fig. 3 all contribute to an assessment of the turbulence parameter n, Fig. 3 showing some relationships characterizing the propagation medium for an over-the-horizon remote sensing circuit. Fig. 4 illustrates the table of Fig. 3, i.e. the theo¬ retical relationship between the measured quantity and the slope n of the refractive-index spectrum ( (K) ~κ- n ) The curves show the degree to which the spectrum slope affects the measured quantity. Length of the arrows shows the effect of refraction on the n deductions; in other words to what extent the various measurable radio quantities contribute to an assessment of the turbulence parameter n and to the refraction parameter a. We see, as an example, that measuring the circuit band¬ width gives us information about the turbulence parameter. Consequently focus is placed on this bandwidth parameter since this represents an inexpensive and simple method. With reference to this simple equation relating the turbulence spectrum slope n to the measurable quantity Δω (and to Fig. 4) , it is clear that this bandwidth parameter is a well pronounced descriptor for the turbulence parameter n. The turbulence parameter (slope of spectrum) can be deduced from the measurement of bandwidth, as will be seen from the following relationships regarding bandwidth Δω (see also Fig. 1) . m = SPECTRUM SLOPE (σ - φ(K)~K "m c = WAVE VELOCITY (3 X 10 8 M/sek.) 2R = DISTANCE between transmitter and receiver θ = SCATTERING ANGLE Furthermore, with reference to Fig. 2 there is a dramatic difference between the K' n turbulence spectra to be expected from CAT, respectively aircraft vortex. The wavenumber corresponding to the horizontal gradient (max φ(K)) for wind produced turbulence is in the order of 2π/400 (at 400 m height) whereas that for aircraft vortex is larger by an order of magnitude. In order to couple to the appropriate turbulence scale K, we can control the following parameters: Scattering angle θ - Carrier frequency F Beat frequency ΔF by measuring (E(K) and E(K+ΔK) and compute the covariance function between the two (E(K) E*(k+ΔK) . With this methods we couple to a scale Δ = c 2Δ or to a wavenumber K-- 2π _ 2πΔ AL (For details regarding the ΔK method, the reader is referred to ref. 4) . As an example illustrating this, consider a forward scatter system with the following specifications. - Path length 700 m - Height of common volume 400 m Scattering angle Carrier frequency F = 10 GHz. In order to couple to the CAT and height 400 m we require the following conditions to be fulfilled: ΔK = 2Δk sin θ/2 _2π_ =2 _2πΔ sinθ/2 400 C ΔF ** = 3.4 MHZ In order to obtain max coupling to the vortex structure on the other hand, we need a AK-- 2π 30m requiring a beat frequency ΔF = 45 MHz Having considered the spatial turbulence spectrum —* characterized by the simple ø(K) - K "n function where the parameter n is determined by a set of measuring techniques involving measurement of bandwidth, delay spectrum, spatial correlation etc, we shall now consider another domain, namel that of Doppler. Measurement of Air Motion by Doppler Considerations: Up to this point we have focused the attention on 3 dimensions of the problem, namely those of space. We have measured received field as a function of frequency (wave- number K) and deduced the spatial distribution of the scatte¬ ring elements We shall now concentrate on the fourth dimension, namely that of time or Doppler. The Doppler shift ω D which an EM wave is subjected to when scattered by scattering elements with velocity V is given by → —► ω D (K) = K * V As already shown, the scattering process selects spatial irregularities of scale L or wavenumber L given by K = k, - k_ where k1. is the wavenumber of the incident wave and ks„ the corresponding one for the scattered wave. Similarly, using a spatial interference pattern obtained by the use of several frequencies from one or more antennas we have ω D (ΔK) = ΔK * V Hence, in order to measure the angular distribution of the wind field (3 D motion of the scattering elements) , we shall have to vary the direction of the wavenumber ΔK so as to measure the vector velocity V. Consider now a bistatic system arranged transversally to the glidepath, where the beams are oriented in elevation and azimuth such that the scattering plane contains the glidepath. Hence the vector K = k ± - k s should be parallel to the glidepath. Using now two transmitting antennas spaced d meter apart with different frequencies, an interference pattern in the plane of the antennas will result as illustrated in Fig. 5. According to this example, when two closely spaced frequen¬ cies (3% spacing) are transmitted into the spatial volume, an interference pattern results. The direction of ΔK and its magnitude is dependent on ω, Δω, d and position R,θ in space. The Doppler shift caused by air motion is maximum when —>• the wind is normal to the interference line (V parallel to K), and obviously zero when parallel to the interference lines. Changing the spacing d between the antennas or between the frequencies Δω, has the same effect, as illustrated in Fig. 6. Thus, use is made of an interferometric radar system (ΔK system) with two apertures spaced distance d apart and having a frequency difference ΔF, a ΔK pattern with different orientation and spacing can be generated. It can be shown (as inferred from Fig. 6) that the -+ direction of the interference wavenumber ΔK is given by φ=dirΔΪ=tan' 1 (- --≤) ΔF R whereas the magnitude of ΔK (inverse of spacing between interference lines) is (this applies for narrow beams such that the scattering takes place close to the center line where φ is small,* see Fig. 6) . Consider as an example that we wish to measure the wind vector V in the glidepath plane. We make use of two antenna spaced d and two frequencies spaced ΔF so as to produce two sets of ΔK lines. One set is tilted 45° in the other direction. Hence F _ R L F~ d Under such conditions the spacing L between the interference lines (two sets of 90° difference in orientation) becomes Hence, for a wind in the direction of ΔK.,, the Doppler shift becomes _ π (ΔJ =-^Δ -V 2π £ D AK) =j2^V whereas the other orthogonal ΔK channel will give zero. In this context advantages are obtained according to the invention, in that the receiver has several receiving aper¬ tures and comprises a coherent homodyne demodulator that gives a directional Doppler shift in a plane orthogonal to the plane through the transmitter and receiver antennas and the scattering air volume. When two or more coherently related microwave frequencies are transmitted, the receiver preferably has a separate receiving aperture for each trans- mitted frequency and a signal processor for multiplying a homodyned signal at each aperture and frequency to generate a time series having a spectrum which gives the directional Doppler shift in any plane through the scattering air volume. Concluding this section on Doppler with the beam geometry in the glidepath plane we note that the particular ΔK system allows one to measure the wind velocity distri¬ bution in the plane of the glidepath. If there is a wind shear this leads to several Doppler lines the strength of which are determined by the relative size of the air strata of different speed, and by the refractivity variance within the strata. Up to this point we have constrained the discussion to the consideration of turbulence spectrum φ(K) and motion pattern P(ω) within a spatial region (air volume under inves¬ tigation) governed by the beam geometry, namely that illumi¬ nated by the transmitting beam and "seen" by the receiver. Through such considerations we analyze an area of some 100 m deep and 1000 m wide for a pathlength of 7 km. If the beam geometry is confined to a vertical plane, this means that the height interval covered around the glidepath is ± 50 m whereas the width with the glidepath as the centerline is ± 500 m. If spatial resolution in excess of this is a require¬ ment, we shall need a number of frequencies which is deter¬ mined by resolution in range and by the size of the common volume and we shall need a number of antennas which is deter¬ mined by the cross range resolution required and by the cross range size of common volume. This topic is covered in the following description referring to Figs. 7, 8 and 9. Resolution in Space, Imaging Capability. A Brief Summary of the Basic Foundation. Reference is first made to Fig. 7 showing the general multifrequency bistatic radar system according to this invention, with indications of magnitudes and variables to be discussed below. In the principle Fig. 7 corresponds to Fig. 1. Main blocks and functions illustrated in Fig. 7 are transmitter 71, receiver represented by homodyne mixer 72, complex multiplication 75 and processor 77, the two latter being in Fig. 1 incorporated in processor unit 25. Delta-k- processing as referred to in this description may take place in complex multiplication unit 75. How to achieve a 1- or 2- dimensional wavenumber representation of the scattering distribution in the common air volume 70, is readily apparent from Fig. 7. The third dimension is obtained by turning transmitter antennas 71A,B,C by 90°. We homodyne (complex conjugate multiplication) the received signal V R (K) e ωt - κ'x with that transmitted V T (K) e^ ωt thus obtaining V(K) e~- K i ' x . We then multiply each of the N carrier frequencies with each other obtaining N(N-l)/2 ΔF channels (Delta-k- processing) . Through this complex multiplication process we get <V*(K) V(K+ΔK)> K and obtain a phase angle e " ^ Δκ*x where The scattering scene (target or air volume 70) is then characterized in the Fourier domain as a hologram: <V*(K) V(K+ΔK)> K - φ(ΔK) ~ j σ(r) e - |UA r dr To image the scene we compute the inverse Fourier trans¬ form _ _ σ(r) - / φ(ΔK) e jΔκ dΔK The N carrier frequencies are spaced according to the Golo b ruler, and so are the receiving antennas 72A. In order to image a scene 70 of size X 0 Y 0 down to a resolution ΔXΔY, the following conditions shall have to be fulfilled for backscatter: ΔF. MAX 2ΔX R Range to scene d Distance between nearest antenna elements UN 2X Distance between extreme antenna elements cR MIN FΔY = d c = Wave velocity cR MAX FY = D ΔF = Frequency difference A numerical example of this onostatic 2-D imaging capa¬ bility through spaced frequencies and spaced antennas is shown in Figs. 8 and 9. According to Fig. 8, by employing a Golomb scheme for frequency spacing as well as for antenna spacing, the number of antennas/frequencies is reduced substantially (N = ~^ψ~ ) . In the current example use is made of 16 antennas and 16 frequencies giving a hologram (matrix) with 16(16-1) χ 16(16-1) χ4 QQQ elements (textels) . Fig. 9 illustrates performing a 2-D Fourier transformation of the two-dimensional ΔK spectrum (hologram) , whereby the spatial window function results. The resolution in range is 2AF M, AX whereas the range of interest with no ambiguities is c 2ΔM Λ£IW Similarly the corresponding cross range characteristic is Rc and Rc To get some feeling for the general dimensions involved with spaced frequencies and spaced antennas, let us present an example. Consider a scene in the air: - Size of scene Cross range dimension 1000 m Dimensions in range 4000 m Zλ M « AXy = F „, χ . =1. 5m - cR _ 'Mrm™= ■ ~15m FX„ With forward scattering ΔX and X 0 should be multiplied with sin θ/2 and the Y component should be multiplied by cos 2 θ/2 where θ is the scattering angle. The number of frequencies required is N r = N< 2 ^ = 9 AX The number of antenna elements required is The motion pattern is obtained through a Fourier trans¬ formation process in time φ(ΔK,ω) - j σ(r,t) e " ^ ωt dt. Hence we obtain a Doppler frequency ω associated with each point in the 2-D Fourier space (each textel) as ω D = ΔK • V In order to obtain the velocity pattern in physical space, i.e. velocity vector for a pixel, a sequence of condi¬ tional operations shall have to be performed involving changing ΔK in a systematic fashion so as to map the motion field in relation to the image structure. This requires sequential introduction of Doppler filters so as to resolve the ambiguities. Characterization of Air Motion (CAT, Wind Shear, Down Draft and Aircraft Wake Vortex) Based on Matched Illumination Concepts, an Example. Let us, in the form of an example specify a set of conditions which one can meet when flying down the glidepath. It is the purpose of this endeavour to form the basis of an assessment of the general technique in regard to distinguish between the various air dynamics factors affecting air safety. Distinction Based on Measurement of the Turbulence Spectrum φ(K) - K "n Through Measurement of Bandwidth: We illuminate the common volume in a bistatic arrange¬ ment with the scattering plane containing the glidepath such that K = k. - k is parallel to the glidepath. We space the frequencies F f to F M in such a way so as to obtain N(N-l) different frequency spacings ΔF according to the Golomb ruler. As an example we use 6 values of F around 10 GHz so as to obtain 15 different frequency spacings ranging from 0.5 MHz to 6 MHz. Based on the specification set given above, we can compute the φ(ΔK) spectrum (bandwidth and center frequency) by calculating φ(ΔK) - E(K) E*(K+ΔK) If we express the turbulence spectrum on the form φ(K) - P(K) ~ K "n and the bandwidth δK that determined by a 10 dB reduction in power, we have the following relationship: P(K a ) K, δK=K o (10 -l) Reference is made to Fig. 10 showing an approximate assessment of the turbulence ΔK spectrum for respectively clear air turbulence at height 400 m and 100 m, respectively, and also for the spatially exponentially damped wake vortex. In Fig. 10 the turbulence spectra φ(K) are plotted for various values of n and height of glidepath H. For compari¬ son, with normalized amplitude, the corresponding spectrum i plotted for the wake vortex assuming that the vortex pheno¬ menon with eddy size 30 m is damped to 1/e over a distance o 300 respectively 150 m. Note that there is a dramatic difference between the normalized CAT signature and that of the wake vortex. In practice there is also expected to be a 10 dB differ ence in spectral intensity between the two phenomena in that vortex is up by 10 dB relative to CAT. This concludes the → discussion on turbulence spectra φ(ΔK) . Measurement of Doppler: The Doppler shift f D is given by f D = -^ ΔK • V whereas the Doppler broadening is Coupling to the CAT at height 400 m we need a wavenumber κ = k i " k s = 400 * If the velocity variance is 2 m this leads to a Doppler broadening of 5 x 10 "3 Hz whereas the vortex at wavenumber 2π 3027? and an orbital velocity of 20 m/sec gives a dramatically higher Doppler broadening of 6 x 10 "1 Hz. This is illustrate in Fig. 11 (Motion Pattern Doppler) . The marked difference in scale size between the vortex and CAT phenomena (30 m and 400, respectively) coupled to the difference in orbital motion leads to a dramatic difference in Doppler broadening. Finally, confining the discussion to Doppler conside- ration, let us consider the case with wind shear. To couple to the air motion at height H a frequency ΔF given by ΔF= 2Jϊsinθ/2 should be used. If the wind at height H moves with a velocity V along the direction of the glidepath, the Doppler shift of the carrier frequency F will be or whereas the corresponding Doppler shift of the beat frequency is With the assumed wind shear of 10 m/s front wind above 1200 FT and 10 m/sec tail wind below this leads to Doppler shifts of f D (K) = 73 Hz whereas the ΔK Doppler shift is f D (ΔK) = 0.02 Hz In order to optimize the S/N ratio, the two F channels should be subjected to narrow filters prior to the V(F) V*(F+ΔF) multiplication. After the multiplication this multiplied signal should be subjected to further filtering and integration. Through this process, the S/N ratio is enhanced by the factor F ΔF It is to be noted at this point that two frequency channels and/or two antennas or apertures represent embodi¬ ments being within the scope of this invention. Narrow bandpass filtering in association with the processor means or unit in the receiver, is then contemplated. Finally, note that in order to evaluate the 4. wind phenomenon, namely that confined to vertical air motion, our beam geometry should be confined to a vertical plane and the common volume should include the glidepath. A down draft of 10 m/sec would then lead to a Doppler shift of 73 Hz on the 10 GHz carrier frequency and 0.02 Hz for the beat frequency at 3.3 Hz.

Claims: Claims 1. A system for the detection and measurement of atmos¬ pheric air movement irregularities such as wind velocity vector, clear air turbulence, aircraft induced vorteces and turbulence, by means of electromagnetic waves, whereby an air volume (10,70) under investigation is illuminated by a trans¬ mitter (1,1A,1B,71,71A,B,C) with a beam of coherent electro¬ magnetic energy (IE) , and a resulting wave field (2E) is received and subsequently coherently demodulated and processed in processor means (25) to derive information on the existence of the said atmospheric irregularities and furthermore to give specific measurements of related para¬ meters, characterized in that at least one receiver (2,2A,2B,72) for said resulting wave field (2E) , which is due to scattering in said air volume (10,70) , is positioned at a bistatic location having a sele¬ cted distance (2R) from said transmitter (1,1A,1B,71,71A,B,C) , that said electromagnetic energy (IE) is transmitted at microwave frequencies in a concentrated beam, from one or more locations, and that transmitter and receiver antennas (1A,1B,2A,2B,71A,B,C,72A) , respectively, are directed so as to cover said air volume (10,70) with a scattering angle (θ) chosen to have such a (low) value that scattering from athmospheric irregularities in said air volume (10,70) has a higher power than noise sources in the system. 2. A system according to claim 1, characterized in that said scattering angle (θ) is selected to be lower than 20°. 3. A system according to claim 1 or 2, characterized in that the transmitter and/or receiver antennas (1A,1B,2A,2B,71A,B,C,72A) are arranged to give a supressed radiated power (a null) in the direction of the other antennas (2A,2B,1A,1B,72A,71A,B,C) to reduce the direct leakage of electromagnetic energy between transmitter (1,71) and receiver (2,72). 4. A system according to claim 1, 2 or 3, characterized in that the receiver and/or transmitter antennas are adapted to be steered mechanically or electronically to provide a 2- or 3-dimensional map of said atmospheric irregularities, in particular windshear features. 5. A system according to any preceding claim, characterized in that the transmitter (1,71) and receiver (2,72) comprise local oscillators (11,12,21,22) adapted to be locked to a common reference, such as a GPS satellite. 6. A system according to any preceding claim, characterized in that the transmitter (1) comprises means (11,12) for transmitting two or more coherently related microwave frequencies (F.,, F.,+ΔF) . 7. A system according to claim 6, characterized in that separate, spatially displaced antennas (1A,1B) are provided for transmitting said coherently related microwave fre¬ quencies. 8. A system according to claim 6 or 7, characterized in that two receiver channels for said transmitted frequencies, comprise narrow bandpass filtering means (20A,20B) in associ¬ ation with said processor means (25) . 9. A system according to claim 6, 7 or 8 characterized in that the receiver has a separate receiving aperture (2A,2B) and a separate bandpass filter (20A,20B) for each transmitted frequency, and that said processor means (25) comprises a signal processor for multiplying a homodyned signal at each aperture and frequency to generate a time series having a spectrum which gives the directional Doppler shift in any plane through the scattering air volume (10) . 10. A system according to any preceding claim, characterized in that the or each receiver comprises a coherent ho odyne demodulator (21,22) for each channel that gives a directional Doppler shift. 11. A system according to any preceding claim, characterized in that the or each receiver has several receiving apertures and comprises a coherent homodyne demodulator for each aper¬ ture and/or frequency. 12. A system according to any preceding claim, characterized in that said receiver processor means (25,75) is adapted to perform Delta-K processing to achieve 1, 2 or 3-dimensional wavenumber representation of the scattering distribution in said air volume (10,70) . 13. A system according to any preceding claim, characterized in that said processor means (25) is adapted to generate a visual image (26) of said atmospheric irregularities through an inversion algorithm. 14. A method for the detection and measurement of atmos¬ pheric air movement irregularities such as wind velocity vector, clear air turbulence, aircraft induced vorteces and turbulence, by means of electromagnetic waves, whereby an air volume (10) under investigation is illuminated by a trans¬ mitter (1,1A,1B) with a beam of coherent electromagnetic energy (IE) , and a resulting wave field (2E) is received and subsequently coherently demodulated and processed in process means to derive information on the existence of the said atmospheric irregularities and furthermore to give specific measurements of related parameters, characterized in that said electromagnetic energy (IE) is transmitted at microwave frequencies in a concentrated beam from said trans¬ mitter (1,1A,IB) , that said resulting wave field (2E) , which is due to scatte¬ ring in said air volume (10) , is received by at least one receiver (2,2A,2B) positioned at a bistatic location having a selected distance (2R) from said transmitter (1,1A,1A), and that transmitter and receiver antennas (1A,1B,2A,2B) , respectively, are directed so as to cover said air volume (10) with a scattering angle (θ) chosen to have such a (low) value that scattering from athmospheric irregularities in said air volume (10) has a higher power than noise sources in the system. 15. A transmitter for use in a system according to claim 1, characterized by comprising means (11,12) for transmitting two or more coherently related microwave frequencies (F.,, F^ΔF) from the same or from different antennas (1A,1B). 16. A receiver for use in a system according to any one of claims 1-13, characterized by comprising processor means (25,75) adapted to perform Delta-K processing to achieve 1, 2 or 3-dimensional wavenumber representation of the scattering distribution in said air volume (10) . 17. A receiver according to claim 16 for use in a system as defined in claims 6, 7 or 8, characterized by comprising a separate receiver aperture (2A,2B) and a separate bandpass filter (20A,20B) for each transmitted frequency, and that said processor means (25) comprises a signal processor for multiplying a homodyned signal at each aperture and frequenc to generate a time series having a spectrum which gives the directional Doppler shift in any plane through the scatterin air volume (10) . 18. A receiver according to claim 16 or 17, characterized i that said processor means (25) is adapted to generate a visual image (26) of said atmospheric irregularities through an inversion algorithm.

Inventors: GJESSING DAG KAARE, HJELMSTAD JENS FREDRIK, GJESSING, DAG, KAARE, HJELMSTAD, JENS, FREDRIK
Applicants: SUSAR AS, SUSAR A/S
Application year: 1993
Main classification: G01S13
IPC codes: G01S13, G01W1