Engineering topics
Publications and source records attributed to Bowhill, S. A..
The Middle Atmosphere Program (MAP) took place from January 1, 1982 through December 31, 1986, and was followed by Middle Atmosphere Cooperation (MAC) through to the end of 1988. The inception and organization of the program are described, together with some of the salient features of its results.
A high resolution numerical technique is used to model the development of a periodically perturbed shear layer imbedded in an initially vertical gradient of a passive scalar. The technique follows the development of the vorticity through an initial linear growth state and well into the nonlinear development of Kelvin-Helmholtz billows, in the zero-viscosity, zero-diffusion limit. The resulting scalar distribution rapidly develops regions of extremely sharp scalar gradients, which wind around the periodically spaced vortical low gradient cores. Vertical cross sections through different parts of the billow structure are presented and compared with rocket measurements of electron density fine structure in the mesosphere. Gradient limits imposed by finite diffusion are calculated, and implications for atmospheric radar observations are discussed.
Various topics related to investigations of the middle atmosphere are discussed. Numerical weather prediction, performance characteristics of weather profiling radars, determination of gravity wave and turbulence parameters, case studies of gravity-wave propagation, turbulence and diffusion due to gravity waves, the climatology of gravity waves, mesosphere-stratosphere-troposphere radar, antenna arrays, and data management techniques are among the topics discussed.
A preliminary search has begun for evidence of tides in the 1-hr average line-of-sight mesospheric velocity data from the Urbana radar in the period 1978 to 1982, inclusive. Observations are restricted to the southeasterly component of those velocities. Since observations are only available for a fraction of a day due to the absence of night-time ionization in the Urbana mesosphere, it was decided to adopt an unusual procedure in the search; namely, to perform a Fourier analysis in the vertical direction and look for rotation in phase of vectors representing spatial frequency components. Propagating tidal modes would then show as vectors with a net rotation corresponding to their downward phase velocity. Five year monthly averages of hourly mean horizontal velocities inferred from the Urbana data are given. Consistent diurnal variation is seen for a number of months. These data were analyzed for vertical spatial periods of 3, 4.5, 6, 9, 12, and 24 km. When plotted as a function of time of day, many of the phasors tended to show a net rotation. The sense of rotation of a spinning vector in the complex plane can be determined objectively by computing the signed area swept out by the vector in saving from point to point. This calculation, indicated downward motion in 8 out of 12 months for the 9-km component and in 11 out of 12 months for the 24-km component. A comparison of the magnitudes of the 6 modes showed that the primary component was the 24-km component, a result confirmed by the tidal models of Forbes.
With a view to determining the role of severe weather in producing gravity waves, two tests were made. In the first, the wind speed measured at two nearby radiosonde stations, Peoria and Salem, was correlated with the stratosphere gravity-wave intensity at Urbana. Although the gravity-wave intensity fluctuated greatly from day to day, these is little if any correlation with the stratospheric wind speed. This suggests that orographic forcing is not a factor in generating gravity waves in Urbana. On the other hand, a clear correlation is found between cloud to heights exceeding 20,000 ft and an increased gravity-wave amplitude in the stratosphere.
Solar flares produce increases in coherent-scatter power from the mesosphere due to the increase in free electrons produced by X-ray photoionization. Thirteen such power enhancements were observed at Urbana. When such an enhancement occurs at an altitude containing a turbulence layer with constant strength, the relative enhancement of electon density is estimated from the enhancement in power. Such estimates of enchanced electron density are compared with estimates of the X-ray photoionization at that altitude, deduced from geostationary satellite measurements. It is found that possible types ion-chemical reaction scheme may be distinguished, and the nonflare ion-pair production function may be estimated. The type of ion-chemical scheme and the nonflare ion-production function are shown to depend on the solar zenith angle.
The Urbana radar, which operates at 40.92 MHz with a peak power of about 1.2 MW into a 100 x 120 m phased array antenna was improved. An accelerated data-acquisition system, a beam-steering system, and a transmit/receive switch were installed. With these changes, the radar is in regular operations for two hours every day around local noon gathering stratospheric and mesospheric data. Special campaigns are mounted in addition under severe weather conditions.
A simple technique for steering the beam of a multimodule phased array MST (mesosphere, stratosphere, troposphere) radar antenna is described. It is desirable to be able to point the antenna in multiple directions, so as to derive all components of the horizontal velocity. This was done on an experimental basis by adding parallel wire line to the feed and achieving a southward tilt of the antenna. It is proved possible to steer the beam through most of the available range without adversely affecting the VSWR seen from the transmitter. Calibrating the antenna direction can be accomplished by observing radio sources, though there are an inadequate number to cover all directions. For various assumed values of the aspect sensitivity in dB/deg, and the calculated antenna pattern, it is possible to calculate the effective pointing angle of the antenna, defined as that angle which would give an identical location for the centroid of the power spectrum if aspect sensitivity were absent. Using averaged apparent steering directions, eastward and northward winds were calculated for special radar runs simultaneous with 14 balloon launches at Peoria and results are presented.
Range resolution of a conventional pulsed Doppler radar is determined by the scattering volume defined by the transmitted pulse shape. To increase the resolution, the length of the pulse must be reduced. Reducing the pulse length also reduces the transmitted power and hense the signal to noise ratio unless the peak power capability of the transmitter is greatly increased. Improved range resolution may also be attained through the use of various pulse coding methods, but such methods are sometimes difficult to implement from a hardware standpoint. The frequency-hopping (F-H) technique described increases the range resolution of pulse Doppler MST (mesosphere stratosphere troposphere) radar without the need for extensive modifications to the radar transmitter. This technique consists of sending a repeated sequence of pulses, each pulse in the sequence being transmitted at a unique radio frequency that is under the control of a microcomputer. This technique is discussed along with other radar parameters.
A new data acquisition system was put into operation at Urbana in August 1984. It uses a standard Apple 2 microcomputer with 48 k RAM and a standard 5 1/4 inch floppy disk. Design criteria for the system is given. The system was implemented using fig-FORTH, a threaded interpretive language which permits easy interfacing to machine code. The throughput of this system is better by a factor of 6 than the PDP-15 minicomputer system previously used, and it has the real time display feature and provides the data in much more convenient form. The features which contribute to this improved performance is listed.
One of the most troublesome features in MSR (mesosphere stratosphere troposphere) velocity measurements is the determination of unwanted scatterer whose velocity is different from that of the surrounding atmosphere. Aircraft seen in the sidelobes of the antenna are the principal problem. Because coherent integration essentially eliminates echoes with line of sight velocities greater than 10 or 20 m/s, aircraft are seen only when their flight path is almost perpendicular to the line of sight. Then, they give large returns whose velocities may be positive or negative, and certainly different from that of the surrounding air. The glitches in the minute by minute velocity records are quite troublesome in that they may distort the statistics of the velocity. An objective way is therefore needed to remove sporadic points of this kind. For this purpose, the Tukey algorithm is appropriate and has some advantages over averaging. The Tukey algorithm, applied to a data array, uses for each data point the median of it and the two points surrounding it. If the three points form a monotonically increasing or decreasing sequence, the original point is copied with change. However, if the central point is remote from the other two, it is replaced by whichever of the two surrounding points is closest in value. The greatest effect of the Tukey algorithm is on data where the successive points are uncorrelated. Examples are given.
The coherent-scatter technique, as used with the Urbana radar, is able to measure relative changes in electron density at one altitude during the progress of a solar flare when that altitude contains a statistically steady turbulent layer. This work describes the analysis of Urbana coherent-scatter data from the times of 13 solar flares in the period from 1978 to 1983. Previous methods of measuring electron density changes in the D-region are summarized. Models of X-ray spectra, photoionization rates, and ion-recombination reaction schemes are reviewed. The coherent-scatter technique is briefly described, and a model is developed which relates changes in scattered power to changes in electron density. An analysis technique is developed using X-ray flux data from geostationary satellites and coherent scatter data from the Urbana radar which empirically distinguishes between proposed D-region ion-chemical schemes, and estimates the nonflare ion-pair production rate.
The Langmuir probe technique for measurement of electron concentration in the mesosphere is capable of excellent altitude resolution, of order 1 m. Measurements from nine daytime rocket flights carrying an electron density fine structure experiment frequently show small scale ionization structures in the altitude region 70 to 90 km. The irregularities are believed to be the result of turbulent advection of ions and electrons. The fine structure experiment flown by the University of Illinois is described and methods of analyzing the collected data is presented. Theories of homogeneous, isotropic turbulence are reviewed. Power spectra of the measured irregularities are calculated and compared to spectra predicted by turbulence theories.
A real time collection system was developed for the Urbana coherent scatter radar system. The new system, designed for use with a microcomputer, has several advantages over the old system implemented with a minicomputer. The software used to collect the data is described as well as the processing software used to analyze the data. In addition a magnetic tape format for coherent scatter data exchange is given.
Various topics relative to middle atmosphere research were discussed. meteorological and aeronomical requirements for mesosphere-stratosphere-troposphere (MST) radar networks, general circulation of the middle atmosphere, the interpretation of radar returns from clear air, spaced antenna and Doppler techniques for velocity measurement, and techniques for the study of gravity waves and turbulence are among the topics discussed.
The performance of an mesosphere-stratosphere-troposphere (MST) radar can be characterized by its system sensitivity and its range resolution. The former enables Doppler velocities to be determined even in the presence of very weak structures; the latter permits study of the fine structure within a turbulent region. Coding of transmitted signals has as its aim an increase in the effective radar sensitivity or range resolution without an increase in the peak transmitted power. This is accomplished by spreading the power in the frequency domain, giving better range resolution, without reducing the pulse width. Two basic techniques are used to accomplish this frequency dispersion: (1) using a type of pseudorandom code for the phase or amplitude within a single pulse, or within a finite sequence of pulses; (2) to code the frequency of the transmitted signal in some way. The various possibilities are discussed and are compared with the pulse-coding methods.
Real-time fast fourier transformer (FFT) processing of a MST radar data and cost-effective approaches to hardware FFT generation were studied. Previously devised hardware FFT configurations are described including the estimated number of chips used and the time required to perform a 1024-point FFT. The remaining entries in the table correspond to original designs, which presuppose the availability of a microcomputer and a modestly complicated hardware peripheral. These original designs, all of which implement a radix-4 FFT with twiddle factors, are assigned model numbers to make them easier to refer to.
As a first step in interchange of data from the Urbana mesosphere-stratosphere-troposphere (MST) radar, a sample tape has been prepared in 9-track 1600-bpi IBM format. It includes all Urbana data for April 1978 (the first month of operation of the radar). The 300-ft tape contains 260 h of typical mesospheric power and line-of-sight velocity data.