ALBERT

Description: The long term trends (least square linear regression with time) of ozone content at seven European, seven North American, three Japanese and two tropical stations during 21 years (1964 to 1984) are analyzed. In all regions negative trends are observed during the 1970s, but are partly compensated by limited periods of positive trends during the late 1960s and late 1970s. Solely the North American ozone data show negative trends in all 10 year periods. When the long term ozone trends are evaluated for each month of the year separately, a seasonal variation is revealed, which in Europe and North America has largest negative trends in late winter and spring. While in Europe the negative trends in winter/spring are partly compensated by positive trends in summer, in North America the summer values reach only zero, retaining the significant negative trend in annual mean values. In contrast to the antarctic ozone hole, the spring reduction of ozone in Europe and in North America is associated with stratospheric temperatures increasing in the analyzed period and therefore is consistent with the major natural ozone production and loss processes.

Description: Rocket payloads designed to measure small scale electron density irregularities and ion properties in the middle atmosphere were flown with each of the three main salvos of the MAC/Epsilon campaign conducted at the Andoya Rocket Range, Norway, during October to November 1987. Fixed bias, hemispheric nose tip probes measured small scale electron density irregularities, indicative of neutral air turbulence, during the rocket's ascent; and subsequently, parachute-borne Gerdien condensers measured the region's polar electrical conductivity, ion mobility and density. One rocket was launched during daylight (October 15, 1052:20 UT), and the other two launches occurred at night (October 21, 2134 UT: November 12, 0021:40 UT) under moderately disturbed conditions which enhanced the detection and measurement of turbulence structures. A preliminary analysis of the real time data displays indicates the presence of small scale electron density irregularities in the altitude range of 60 to 90 km. Ongoing data reduction will determine turbulence parameters and also the region's electrical properties below 90 km.

Description: A real time winds (RTW) system from Saskatoon operated with the Tromsoe M.F. (partial reflection) radar on a continuous basis, June to December 1987. Profiles with 3 km resolution were obtained every 5 minutes with weak ionization, and few geomagnetic disturbances limited the observations normally to 80 to 110 km. However, daily mean winds, tidal characteristics (24, 12 h) such as amplitudes, phases and wavelengths, and gravity wave characteristics (intensities, mean directions) are available throughout this interval, which includes MAC-SINE and Epsilon. This is particularly valuable in defining the background state for some experiments, e.g., rockets, and for comparison with related parameters from the lidar and other radars (EISCAT, SOUSY-VHF). Comparisons with dynamical parameters from Saskatoon (52 N) are made: the zonal circulation was weaker at Tromsoe, tidal amplitudes smaller, and summer 12 h tidal wavelengths shorter (approx. 80 km vs approx. 100 km). The fall transition for this tide occurred in September, earlier than observed elsewhere. Initial comparisons with other experimental systems are also made.

Description: Wind measurements near the mesopause level were provided in Badary (Eastern Siberia, USSR) by the D1 method using a low frequency (200 kHz) radio transmitter during 1976 to 1986. The seasonal variations of wind level are analyzed and compared with other midlatitude wind measurements and with some theoretical preditions. It is shown that tidal parameters depend not only on latitude but on the longitude as well. The annual variation of diurnal and semidurnal tide (zonal and meridional) have maxima in summer and minima in winter. There are significant differences between annual variations of phases for diurnal and semidiurnal tides. These differences are especially distinct for summer and autumn. There is no systematic seasonal variation of the phase of the semidiurnal tide (average value similar to 6 hours).

Description: In the winter of 1983 to 1984, the research institutes of the Soviet Union took an active part in the accomplishment of the project Winter in Northern Europe (MAP/WINE) of the Middle Atmosphere Program. Different methods were used to measure temperature, direction and velocity of wind, turbulence, electron concentration in the lower ionosphere, and radio wave absorption. The study of the stratopheric warmings and the related changes in the mesosphere and lower ionosphere was considered of special importance. The analysis of the obtained data has shown, in particular, that during the stratospheric warmings the western wind in winter time becomes weaker and even reverses. At the same time period the electron concentration and the radio wave absorption in the lower ionosphere are often reduced. It is also observed that the high absorption zones move from west to east. These results confirm the concept about the role of the cyclonic circumpolar vortex in the transport of the auroral air to temperate latitudes and about the appearance of conditions for the winter anomalous radio wave absorption.

Description: Electric fields and currents in the middle atmosphere are important elements of the modern picture of this region. Balloon instruments, reaching the level of the stratosphere, were used extensively for the experimental work. The research has shown good progress, both in the MAP period and in the years before and after. The knowledge was increased about, e.g., the upper atmosphere potential, the electric properties of the medium itself and about the coupling with magnetospheric (ionospheric) fields and currents. Also various measurements have brought about a discussion of the possible existence of hitherto unknown sources. Throughout the MAP period the work on a possible definition of an electric index has continued.

Description: The results of the restoration of the optical parameters of the El Chichon aerosol cloud at Lambda 2.2 micron are presented. The eight channel teleradiometer FAZA for visible and near infrared spectral regions is described in some detail. Mathematical aspects of data processing for the ill-posed problems are discussed. The results of the restoration of the volume emission rate of the molecular oxygen on Lambda 1.27 micron in the middle atmosphere and of the aerosol on Lambda 0.84 micron are presented.

Description: MAP/GLOBUS 1983 was a project for the study of stratospheric trace gases and dynamics. A respective field campaign was performed in September/October 1983 in Western Europe. A large number of measurements were taken by instruments based on the ground, on airplane, balloons, and satellite. The structure of the campaign is described, and a survey of the results are given.

Description: Results of wind measurements at the midlatitude ionospheric D region are presented. The wind regime of the lower thermosphere is rather sensitive to stratospheric temperature variations, especially to sudden stratospheric warmings. The longitudinal effect in D region dynamics was revealed on the basis of simultaneous wind measurements at some points located practically at the same latitude but in different climatic regions. The distance differences are observed in the statistical distributions of wind parameters, during winter the average zonal wind speed over East Siberia was about twice that over Central Europe, and the semidiurnal zonal tide is weaker over East Siberia. The data on the seasonal reconstruction of circulation and the response of the D region wind field to the stratospheric warmings depend on the intensity and locations of stratospheric disturbances in relation to the observatory. These experimental facts are interpreted as a meteorological control of the D region and as a dependence of the lower thermosphere dynamics on the conditions of dissipation of internal waves propagating from the troposphere and stratosphere.

Description: The equatorial region has attracted many MAD studies mainly based on data of limited locations and resolutions. Established at NIEO are: (1) Climatology of the equatorial middle atmosphere (all of the mean zonal flow, the meridional and/or east-west circulations and the planetary/gravity waves are described based on massive, reliable data statistics); (2) Troposphere-stratosphere coupling at the equator (the candidate location of NIEO is just at the stratospheric fountain area where the tracers and waves are pumped up into the middle atmosphere); and (3) Mesosphere-thermosphere coupling at the equator; thermospheric superrotation, which may be caused either by ion drag or by tidal breaking, is examined in detail by observations covering a wide altitude range from the mesosphere through the thermosphere.

Description: Incoherent scatter radar (ISR) has become the most powerful means of studying the ionosphere from the ground. Many of the ideas and methods underlying the troposphere and stratosphere (ST) radars have been taken over from ISR. Whereas the theory of refractive index fluctuations in the lower atmosphere, depending as it does on turbulence, is poorly understood, the theory of the refractivity fluctuations in the ionosphere, which depend on thermal fluctuations, is known in great detail. The underlying theory is one of the most successful theories in plasma physics, and allows for many detailed investigations of a number of parameters such as electron density, electron temperature, ion temperature, electron mean velocity, and ion mean velocity as well as parameters pertaining to composition, neutral density and others. Here, the author reviews the fundamental processes involved in the scattering from a plasma undergoing thermal or near thermal fluctuations in density. The fundamental scattering properties of the plasma to the physical parameters characterizing them from first principles. He does not discuss the observation process itself, as the observational principles are quite similar whether they are applied to a neutral gas or a fluctuating plasma.

Description: A review of the history of radar techniques which have been applied to atmospheric observation is given. The author starts with ionosphere observation with the ionosonde, symbolizing as it does the earliest history of radar observation, and proceeds to later developments in radar observation such as the use of partial reflection, meteor, and incoherent scatter radars. Mesosphere stratosphere troposphere (MST) radars are discussed in terms of lower atmosphere observation.

Description: Several tests were designed to determine the correct error variances for the GEM-T1 gravitational solution which was derived exclusively from satellite tracking data. The basic method employs both wholly independent and dependent subset data solutions and produces a full field coefficient by coefficient estimate of the model uncertainties. The GEM-T1 errors were further analyzed using a method based upon eigenvalue-eigenvector analysis which calibrates the entire covariance matrix. Dependent satellite and independent altimetric and surface gravity data sets, as well as independent satellite deep resonance information, confirm essentially the same error assessment.

Description: A natural extension of the recent satellite derived potential coefficient models is the development of high degree (maximum 180 or 360) expansions. Such expansions are based on the combination of the satellite derived models with terrestrial gravity data and satellite altimeter data. Such models are useful for more precise geoid undulation computations, for simulation studies involving different typed of future missions (e.g., gradiometry), and as reference fields for different types of gravimetric computations. The attention is to the effect of the terrain, ellipsoidal terms, and weighting. The basic methods used for the high degree solutions are reviewed. Various correction terms are described and recent models are discussed and compared.

Description: Linear estimation theory, along with a new technique to compute relative data weights, was applied to the determination of the Earth's geopotential field and other geophysical model parameters using a combination of satellite ground-based tracking data, satellite altimetry data, and the surface gravimetry data. The relative data weights for the inhomogeneous data sets are estimated simultaneously with the gravity field and other geophysical and orbit parameters in a least squares approach to produce the University of Texas gravity field models. New techniques to perform calibration of the formal covariance matrix for the geopotential solution were developed to obtain a reliable gravity field error estimate. Different techniques, which include orbit residual analysis, surface gravity anomaly residual analysis, subset gravity solution comparisons and consider covariance analysis, were applied to investigate the reliability of the calibration.

Description: Various aspects of solar-terrestrial science are briefly discussed. The evolution of the concept of the Solar-Terrestrial Observatory (STO) is discussed. Several remote sensing concepts were proposed. The capability to image global atmospheric, ionospheric, and magnetospheric phenomena was thought to hold considerable promise for understanding global processes. All groups of the workshop expressed a desire to develop a means for focusing research activities on data acquired from past space missions. International cooperation, new vantage points for viewing the Earth and the Sun, Explorer-class satellites, Space Shuttle missions, and other astronaut-associated space platform missions are discussed.

Description: The general strategy to advance knowledge of the ionospheric component of the solar terrestrial system should consist of a three pronged attack on the problem. Ionospheric models should be refined by utilization of existing and new data bases. The data generated in the future should emphasize spatial and temporal gradients and their relation to other events in the solar terrestrial system. In parallel with the improvement in modeling, it will be necessary to initiate a program of advanced instrument development. In particular, emphasis should be placed on the area of improved imaging techniques. The third general activity to be supported should be active experiments related to a better understanding of the basic physics of interactions occurring in the ionospheric environment. These strategies are briefly discussed.

Description: By the early 1990s, magnetospheric physics will have progressed primarily through observations made from Explorer-class spacecraft, sounding rockets, ground based facilities, and shuttle based experiments. The global geospace science (GGS) element of the International Solar Terrestrial Physics program, when combined with contributions to the ESA Cluster mission and ground based and computer modeling programs, will form the basis for a major U.S. initiative in magnetospheric physics. The scientific objectives of the GGS program involve the study of energy transport throughout geospace. The Cluster mission will investigate turbulence and boundary phenomena in geospace, particularly at high latitudes on the dayside and in the region of the neutral sheet at geocentric distances of about 20 earth radii on the night side of the earth. The current state of knowledge is reviewed and the goals of these missions are briefly discussed.

Description: Most of the understanding of the thermosphere resulted from the analysis of data accrued through the Atmosphere Explorer satellites, the Dynamics Explorer 2 satellite, and observations from rockets, balloons, and ground based instruments. However, new questions were posed by the data that have not yet been answered. The mesosphere and lower thermosphere have been less thoroughly studied because of the difficulty of accessibility on a global scale, and many rather fundamental characteristics of these regions are not well understood. A wide variety of measurement platforms can be used to implement various parts of a measurement strategy, but the major thrusts of the International Solar Terrestrial Physics Program would require Explorer-class missions. A remote sensing mission to explore the mesosphere and lower thermosphere and one and two Explorer-type spacecraft to enable a mission into the thermosphere itself would provide the essential components of a productive program of exploration of this important region of the upper atomsphere. Theoretical mission options are explored.

Description: In MAP observations, it was found that: (1) gravity waves in selected or filtered portions of data are fit for monochromatic structures, whereas (2) those in fully continuous and resolved observations take universal continuous spectra. It is possible to explain (2) by dispersion of quasi-monochromatic (or slowly varying) wave packets observed locally as (1), since the medium atmosphere is unsteady and nonuniform. Complete verification of the wave-mean flow interactions by tracking individual wave packets seems hopeless, because the wave induced flow cannot be distinguished from the basic flow independent of the waves. Instead, the primitive picture is looked at before MAP, that is, the atmosphere is just like an entertainment stage illuminated by cocktail lights of quasi-monochromatic gravity waves. The wave parameters are regarded as functions of time and spatial coordinates. The observational evidences (1) and (2) suggest that the wave parameter field is rather homogeneous, which can be explained by interference of quasi-monochromatic wave packets.

Description: Recent studies have revealed that vertical wave number spectra of wind velocity and temperture fluctuations in the troposphere and the lower stratosphere are fairly well explained by a saturated gravity wave spectrum. But N(2) (N:Brunt-Vaisala (BV) frequency) spectra seem to be better for testing the scaling of the vertical wave number spectra in layers with different stratifications, beause its energy density is proportional only to the background value of N(2), while that for temperature depends on both the BV frequency and the potential temperature. From temperature profiles observed in June to August 1987 over the MU Observatory, Japan, by using a radiosonde with 30 m height resolution, N(2) spectra are determined in the 2 to 8.5 km (troposphere) and 18.5 to 25 km (lower stratosphere) ranges. Although individual spectra show fairly large day-by-day variability, the slope of the median of 34 spectra agrees reasonably with the theoretical value of -1 in the wave number range of 6 x 10(-4) similar to 3 x 10(-3) (c/m). The ratio of the spectral energy between these two height regions is about equal to the ratio of N(2), consistent with the prediction of saturated gravity wave theory.

Description: In recent measurements of mesospheric winds with the MU radar at Shigaraki (34.9 N, 136.1 E), Japan, a couple of wind velocity data was obtained showing that a monochromatic inertia-gravity wave was propagating vertically in the mesosphere. The data show that the velocity amplitude of the wave did not increase exponentially above a height and further show that the wave form had broken down at a level. Evaluating the changes of the atmospheric stability in the wave field from the measured wind data and from the estimated wave parameters, the breakdown of the wave was shown to be connected with the occurrence of convective instability in the wave field.

Description: The Chung Li Radar (24.91 N; 121.24 E) has been operating since 1986. A five beam observational configuration was used on a regular basis to study the various dynamics processes in the atmosphere-lower stratosphere height region. Due to its geographical location, the annual Typhoon and Mei-Yu seasons provide good opportunities to study the various interesting dynamic processes such as instabilities, generation of gravity waves, wave mean field interaction, etc. Three dimensional air motions due to these fronts are presented. Special cases of gravity wave generation, propagation and their effects on the turbulent layers are discussed.

Description: On the nights of October 21 and 28, 1987 (UT), two Nike Orion payloads (NASA 31.066 and 31.067) were launched from Andoya, Norway, as part of the MAC/Epsilon campaign, to study auroral energetics and their effect on the middle atmosphere. Each payload carried instrumentation to measure relativistic electrons from 0.1 to 1.0 MeV in 12 differential channels, and Bremsstrahlung X-rays from greater than 5 to greater than 80 KeV in 5 integral channels. In addition, instrumentation to measure ion densities and electric fields were also included on these and, in the case of 31.066, on other near simultaneous payloads. The first flight, 31.066, was launched under pre-magnetic midnight conditions during relatively stable auroral conditions. Flight 31.067 was launched during post-breakup conditions at which time pulsations of approx. 100 seconds duration were evident. The measured radiations including their spectral characteristics are compared for these two events, to appraise their effect on the electrodynamic properties of the middle atmosphere as determined by other rocket-borne measurements.

Description: The EISCAT incoherent scatter radar facility in Tromsoe, Norway was operated during the MAC/SINE campaign for 78 hours in the period 10 June to 17 July 1987, and during the MAC/Epsilon campaign for 90 hours in the period 15 October to 5 November 1987. The VHF (224 MHz) radar operations during MAC/SINE yielded most interesting observations of strong coherent echoes from the mesopause region. Characteristic data of these polar mesospheric summer echoes are presented. The UHF (933 MHz) radar operations during MAC/Epsilon were done with 18 deg off zenith beam and allows the deduction of meridonal and horizontal wind components as well as radial velocity spectra in addition to the usual electron density profiles in the D and lower E regions. Some results from the VHF and UHF radars indicating the presence of gravity waves are examined.

Description: The continuing series of horizontal wind measurements by the spaced-antenna real time winds (RTW) method was supplemented by a phase coherent system for two years. Vertical motions are inferred from the complex autocorrelation functions, and an RTW system provides 5 min samples from 60 to 110 km. Comparisons with full interferometric 3-D velocity measurements confirm the validity of this approach. Following comparisons and corrections with the horizontal winds, mean summer and winter (24 h) days of vertical motions are shown. Tidal fluctuations are evident. In summer the motions are downward, consistent with data from Poker Flat, and the suggestion of Coy et al. (1986) that these represent Eulerian motions. The expected upward Lagrangian motion then results from adding up upward Stokes' drift. The winter motions are more complex, and are discussed in the context of gravity wave fluxes and possible meridional cells. The divergence of the vertical flux of zonal momentum is also calculated and found to be similar to the coriolis torque due to the meridional winds.

Description: Recently an interesting advancement in the study of nonmigrating tides has occurred. There have been two distinct approaches in works on this subject. One is based on mechanistic models as considering nonuniform global distribution of water vapor or heating only on land, solving a set of linear equations. It is found that insolation absorption of the nonuniformly distributed water vapor produces only weak nonmigrating tides in the lower thermosphere; the planetary boundary layer heating on land can explain the enhanced tides on land and those with short vertical wavelengths in the stratosphere. The other approach is novel and uses simulation on the general circulation model (GCM). This realistic model can reproduce tides globally and in many details. The enhancement of two nonmigrating modes as eastward traveling modes with a wave number 3 and westward traveling modes with a wave number 5 is in surprisingly good agreement with observation at sea level, at 700 mb and even at 300 mb.

Description: Continuous radar measurements of the wind oscillations caused by the solar diurnal and semidiurnal tides in the 80 to 110 km region of the atmosphere at the geographically conjugate stations of Adelaide (35 S, 138 E) and Kyoto (35 N, 136 E) are compared for the period 1983 to 1985. At the solstices it is found for both the 24- and 12-hr tides that the NS and EW oscillations tend to be in-phase and out-of-phase, respectively, behavior which indicates strong tidal asymmetries. The asymmetries in the 12-hr tide are consistent with a strong contribution from the (2,3) mode while the asymmetries in the 24-hr tide are ascribed to the effects of mean winds and dissipation acting to distort the (1,1) mode as it propagates up through the middle atmosphere.

Description: Middle atmosphere temperatures were measured during the MAP/WINE campaign by various ground-based techniques, by rocket instruments, and by satellites. Respective data were analyzed for atmospheric thermal mean state as well as for long and short period variations. A brief survey of the results is given. Monthly mean temperatures agree well with the new CIRA model. Long period (planetary) waves frequently exhibit peculiar vertical amplitude and phase structures, resembling those of standing waves. Short period oscillations tend to begin breaking well below the stratosphere.

Description: During MAP/WINE small scale structure and turbulence in the mesosphere and lower thermosphere was studied in situ by rocket-borne instruments as well as from the ground by remote sensing techniques. The eight salvoes launched during the campaign resulted in a wealth of information on the dynamical structure of these regions. The experimental results are reviewed and their interpretation is discussed in terms of gravity waves and turbulence. It is shown that eddy diffusion coefficients and turbulent energy dissipation rates may be derived from the in situ measurements in a consistent manner. The observations are also shown to be consistent with the hypothesis that turbulence can be created by a process of gravity wave saturation.

Description: The project Winter in Northern Europe (WINE) of the international Middle Atmosphere Program (MAP) comprised a multinational study of the structure, dynamics and composition of the middle atmosphere in winter at high latitudes. Coordinated field measurements were performed during the winter 1983 to 1984 by a large number of ground-based, air-borne, rocket-borne and satellite-borne instruments. Many of the individual experiments were performed in the European sector of the high latitude and polar atmosphere. Studies of the stratosphere, were, in addition, expanded to hemispheric scales by the use of data obtained from remotely sensing satellites. Beyond its direct scientific results, which are reviewed, MAP/WINE has stimulated quite a number of follow-on experiments and projects which address the aeronomy of the middle atmosphere at high and polar latitudes.

Description: Long term observations of relativistic electrons in the earth's outer magnetosphere show a strong solar cycle dependence with a prominent intensity maximum during the approach to solar minimum. This population therefore closely corresponds to the presence of high speed solar wind streams emanating from solar coronal holes. Using a numerical code, the precipitating electron energy deposition in the earth's upper and middle atmosphere were calculated. Observed events (typically persisting several days) would have maximum effect in the 40 to 60 km altitude range with peak energy depositions greater than 110 keV/cu cm-s. It is suggested that this electron population could play an important long term role in modulating lower D region ionization and middle atmospheric ozone chemistry. Methods are described of observing middle atmospheric and lower ionospheric effects of the electrons including balloon, riometer, and space-based ozone sensor systems. A particularly promising approach may involve the monitoring of global Schumann resonance modes which are sensitive to global changes in the properties of the earth-ionosphere cavity. Present work indicates that Schumann resonance properties are moderately correlated with the flux of precipitating relativistic electrons thus offering the possibility of continuously monitoring this aspect of magnetosphere-atmosphere coupling.

Description: While the Langmuir probe controlled by rocket propagation experiments by the University of Illinois at midlatitude revealed the existence of a permanent D region turning point (DTP), similar measurements over the Thumba equatorial station did not clearly bring out the above daytime feature. Moreover, the calibration constant (ratio of electron density to the current drawn by the Langmuir probe) increased with height (in the 70 to 100 km region) in the case of the midlatitude observations whereas the recent measurements over Thumba showed a decrease up to about 90 km followed by an increase above 90 km. Secondly, there is the problem of reconciling the station oriented observations from the COSPAR family with the ground based radio propagation measurements from the URSI family. Thirdly, new information on Winter in Northern Europe (WINE) and in USSR is available by asking for its incorporation into any global model such as the IRI. The results of investigation of the above aspects are presented.

Description: The reduction is reported of stratospheric ion conductivities in the altitude range of 20 to 27 km attributable to the aerosols injected into the stratosphere by the eruption of volcano Nevado Del Ruiz on November 13, 1985. Three balloon experiments were conducted from Hyderabad, India (17.5 N, 78.6 E) carrying a Langmuir probe payload for measuring stratospheric ion conductivities. The first flight took place about 9 months before the volcanic eruption, the second 3 weeks after the eruption and the third about a year later. Lidar observations from Japan, Hawaii and Europe reported detection of aerosol layers in the 18 to 25 km altitude range attributable to the Nevado Del Ruiz volcanic eruption. A comparison of the conductivity profiles shows that the reduction of ion conductivities is: 57.3 percent at 20 km and 31 percent at 25 km. A year after the eruption, conductivities at all heights tended to recover.

Description: Under MAP, two rockets were launched from Thumba (8.5 N, 76.8 E) around 1030 hrs Lt with identical payloads on 7 and 10 March 1986 for D region studies. Positive ion densities were measured by spherical probe and Gerdien condenser and electron densities were measured by Langmuir probe and propagation experiments. In both flights a valley in ionization height profile was noticed around 83 km. The density of ionization at this altitude was about 4 x 10(2) cu cm. A detailed positive ion-chemical scheme was used to reproduce the measured ionization height profiles. The density of NO needed to reproduce the valley in ionization at 83 km came around 5 x 10(5) cu cm. A photochemical treatment without diffusion process was found inadequate to explain this value of NO. Calculations showed that the value of vertical eddy diffusion needed to reproduce the value of NO was around 10(6)sq cm/s. Interestingly, the same value of eddy diffusion coefficient was obtained when derived in the manner described by Thrane and his coworkers using only the positive ion current data of spherical probes.

Description: The middle atmosphere (MA) has long been known as an absorber of radio waves, and as a region of nonlinear interactions among waves. The region of highest transverse conductivity near the top of the MA provides a common return for global thunderstorm, auroral Birkeland, and ionospheric dynamo currents, with possibilities for coupling among them. Their associated fields and other transverse fields map to lower altitudes depending on scale size. Evidence now exists for motion-driven aerosol generators, and for charge trapped at the base of magnetic field lines, both capable of producing large MA electric fields. Ionospheric Maxwell currents (curl H) parallel to the magnetic field appear to map to lower altitudes, with rapidly time-varying components appearing as displacement currents in the stratosphere. Lightning couples a (primarily ELF and ULF) current transient to the ionosphere and magnetosphere whose wave shape is largely dependent on the MA conductivity profile. Electrical energy is of direct significance mainly in the upper MA, but electrodynamic transport of minor constituents such as smoke particles or CN may be important at other altitudes.

Description: The method for the determination of the optical density of tropospheric and mesospheric aerosols and for the estimation of the scattering phase function in the forward directions (0 to 30 deg) is proposed. The method is based on measurements of the brightness of the twilight horizon with the high resolution limb camera having five separated fields of view. In some detail the features of the camera are described and the mathematical aspects of remote sensing data inversion are discussed.

Description: A pulsed ruby lidar has been in regular operation at the tropical station Trivandrum (8 deg 33 sec N, 76 deg 57 sec E) since October 1986. The lidar data were analyzed to obtain monthly mean aerosol extinction at lower stratospheric altitudes. The monthly mean variation of aerosol extinction shows a peak in March to April with a trough in the winter months. This behavior is compared with that of temperature at the same altitude, obtained from balloonsonde measurements at the same station. It is found that there is a negative correlation between the two with high extinction values corresponding to low temperature values and vice versa. This is attributed to the stratospheric aerosol microphysical processes. The association of stratospheric aerosol extinction with the tropopause altitude and temperature was studied along with the implications of these results.

Description: Considerable progress was made in understanding gravity waves and their effects in the middle atmosphere during the MAP and MAC periods. During this time, gravity waves were recognized to play a central role in controlling the large scale circulation and the thermal and constituent structure of this region through wave transports of energy and momentum, a significant induced meridional circulation, and through the action of wave induced turbulence. Both theoretical and observational studies also have contributed to the understanding of the gravity wave spectrum, its temporal and spatial variability, and the processes responsible for wave saturation. As a result, the propagation, interactions, and detailed effects of such motions in the middle atmosphere are beginning to be understood. An overview is provided.

Description: The Middle Atmosphere Program (MAP) has provided a focus for considerable research on atmospherical radiative, chemical, and dynamical processes and the mutual coupling among these processes. In particular, major advances have occurred in the understanding of constituent transport as a result of near-global measurements obtained during MAP from several satellite based instruments (e.g., LIMS, SAMS, SAGE, and SSU among others). Using selected portions of these data, the development is reviewed of progress in understanding transport processes with special emphasis on dynamically active periods. Examples are presented which demonstrate coupling between chemistry and dynamics. In addition to the constituent data, the use is reviewed of Ertel's potential vorticity, inferred from satellite temperature data, as a diagnostic for interpreting transport phenomena. Finally, the use is briefly illustrated of 3-D model simulations, in conjunction with the satellite data, for providing additional insight into fundamental transport mechanisms.

Description: Two major international campaigns dedicated to the study of middle atmosphere dynamics in high latitudes were successfully completed in 1987. MAC/SINE (Middle Atmosphere Cooperation/Summer in Northern Europe) was carried out during the period 7 June to 19 July, whereas MAC/Epsilon took place in the period 12 October to 15 November. In both campaigns a large number of ground based and rocket techniques were used in a concerted effort to map the dynamical structure of the middle atmosphere over Northern Europe. Although the analysis of the observations has only just started, it is clear that a large and unique data set was obtained, which is believed will provide new insight into the dynamical processes in this interesting region of the atmosphere. A brief overview of the campaigns, their scientific aims, organization and structure is presented.

Description: The monthly mean zonal winds and semidiurnal tides at 80 and 90 km, in January and July, at Christchurch (44 S) for the period 1978 to 1986 are presented. There are significant trends but evidence for solar control of the mean zonal wind and the semidiurnal tide is not conclusive.

Description: The geomagnetic control of ionospheric D region dynamics was revealed and confirmed on the basis of radiophysical wind measurements (1978 to 1983) over East Siberia. The monthly mean parameters of the wind system are different for quiet and disturbed conditions. There is an increase in stability of the meridional wind with increasing level of geomagnetic activity. The influence of geomagnetic storms on the measured wind is considered on the basis of 31 events. There are effects on the phase of the semidiurnal tidal wind, but variations of amplitude are weak. The effect of the geomagnetic storm depends on the intensity and is more clear-cut for the A sub p is greater than 100.

Description: 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.

Description: For the MAP/WINE winter temperature and wind measurements of rockets were combined with SSU radiances (Stratospheric Sounder Unit onboard the NOAA satellites) and stratopause heights from the Solar Mesosphere Explorer (SME) to get a retrieved data set including all available information. By means of this data set a hemispheric geopotential height, temperature and geostrophic wind fields eddy transports for wave mean flow interaction and potential vorticity for the interpretation of nonlinear wave breaking could be computed. Wave reflection at critical lines was investigated with respect of stratospheric warmings. The meridional gradient of the potential vorticity and focusing of wave activity is compared with derived data from satellite observations during other winters.

Description: During the MAP/WINE campaign in winter 1983 to 1984 several instrumental techniques, such as meteorological rockets, sounding rockets, MST radar and incoherent scatter radar, were applied to measure wind velocities in the middle atmosphere. Profiles of mean, tidal and fluctuating wind velocities were obtained up to 90 to 100 km altitude. These are compared with profiles from models, measurements at other locations and at other times as well as satellite derived data. The results are discussed in terms of ageostropic winds, planetary waves, tidal modes and the possibility of a saturated gravity wave spectrum in the mesosphere.

Description: The middle atmosphere is weakly ionized, collision dominated plasma. Rocket-borne dc probes have proven effective in measuring conductivity and electron and ion concentrations in this plasma, and in some cases electron temperatures and ion composition were observed. During MAP, dc probes were also used successfully to study turbulence in the middle atmosphere, using the plasma as a passive tracer for the dynamical processes in the nonionized gas. Basic principles are reviewed of the dc-probe technique and examples are presented of results obtained during MAP.

Description: The recent revival and strong motivation for research in middle atmospheric electrodynamics can be attributed, in large part, to the discovery of large (V/m) electric fields within the lower mesosphere during the decade prior to MAP. Subsequent rocket soundings appeared to verify the preliminary findings. During the MAP era, more sophisticated techniques have been employed to obtain measurements which respond positively to criticisms of earlier results, and which provide more insight regarding the character of the fields. The occurrence of mesospheric V/m electric fields now seems to require the presence of aerosols, of local winds and related dynamics, and of an atmospheric electrical conductivity less than 10(-10)S/m. Furthermore, new theoretical ideas describing the origin of the V/m fields are consistent with the measurements. The current status of results regarding V/m fields in the middle atmosphere is reviewed in light of the more widely accepted electric field structure for this region from rocket, balloon and modeling results.

Description: Prior to MAP, virtually nothing was known about gravity wave spectra in the atmosphere. The development of observational techniques has played a major role in these studies. Radar and lidar have been particularly important since they can measure atmospheric parameters continuously over large height ranges. Some advances made are: (1) The observed fluctuations and power spectra in the free atmosphere are mostly if not entirely due to a superposition of gravity waves, which can be modeled by the Garrett Munk (GM) model; (2) There is no evidence that 2-D turbulence makes a significant contribution to the observed fluctuations. In any case, the agreement between observations and the GM model shows that the 2DT contribution must be relatively small; (3) Spectra versus vertical wave number are saturated at large wave number, with theory and observations indicating that t approximately equals 3; and (4) Vertical velocity fluctuations and spectra measured near rough terrain are strongly contaminated by mountain waves. But over very flat terrain the spectra are dominated by gravity waves at periods shorter than about 6 hours and apparently by synoptic scale velocities at periods longer than 6 hours. Thus it may be possible to study synoptic scale vertical velocities using radars located in very flat terrain.

Description: The aim of the MASH project is to study the dynamics of the middle atmosphere in the Southern Hemisphere, emphasizing inter-hemispheric differences. Both observational data and data from simulations with numerical models are being used. It is intended that MASH will be complemented by parallel studies on the transport and photochemistry of trace species in the Southern Hemisphere. Impetus for such studies has come from the unexpected finding of a springtime ozone hole over Antarctica. A summary of recent progress with the MASH project is given. Data from polar orbiting satellites are used to discuss the large scale circulation found in the Southern Hemisphere at extratropical latitudes. Comparisons are made with that of the Northern Hemisphere. Particular attention is paid to the springtime final warming, the most spectacular large scale phenomenon in the statosphere of the Southern Hemisphere. The circulation before and after this event has to be taken into account in theories for the formation and subsequent disappearance of the ozone hole.

Description: Evidence is cited from these studies and those of others showing the different nature of the yearly variations of the middle atmospheres of the Northern and Southern Hemispheres. The Northern Hemisphere middle atmosphere is shown to be characterized by annual variations in planetary wave number one amplitude and the accompanying ozone transports. The Southern Hemisphere middle atmosphere is shown to be characterized by semiannual variations in the amplitude of planetary wave number one and the accompanying ozone transports. The amplitude of wave number two in both hemispheres appears to vary annually. Examination is made of the nature of the planetary wave forcing in both hemispheres as well as the planetary wave propagation characteristics in both hemispheres in an attempt to better understand this.

Description: The understanding of turbulence in the middle atmosphere has improved considerably during the MAP period. For a theoretical viewpoint, several advances were made including understanding the ways in which turbulence is generated, and the differences between the rates of diffusion of momentum and heat. Experimentally, a proper understanding of how radars can be used to measure turbulence has emerged, and turbulent energy dissipation rates in the middle atmosphere were measured with MF, HF, and VHF radars. New rocket techniques were developed which have enabled detailed studies of the fine structure of turbulence to be made. While some discrepancies between techniques still exist, these will undoubtedly be resolved soon, and these different techniques are already providing a great improvement in the understanding of turbulence on a global scale.

Description: Numerical integrations of the structure function of tropospheric inhomogeneities have been performed to assess the impact of water vapor fluctuations on connected element interferometry (CEI). The expectation value of the RMS troposphere error for a differential spacecraft-quasar observation is derived by integrating the spatial refractivity structure function along raypaths to both the spacecraft and quasar from two spatially separated sites. Correlations between the tropospheric conditions at the two sites, which can become significant for short baseline observations, are fully accounted for in this calculation. Temporal effects are treated by assuming a frozen-flow model in which a fixed spatial distribution blows over both sites. Two nominal observation scenarios are considered, along with variations to study the dependence of the resultant differential troposphere errors on baseline length, observation time, source separation angle, and elevation. Consecutive differential observations are found to be almost completely uncorrelated, implying that averaging many repeated differential observations can quickly reduce the troposphere error.

Description: Comparison of segment development of a number of arcs has shown that consistent relationships between segmentation, volcanism and variable stresses exists. Researchers successfully modeled these relationships using the conceptual model of lithospheric buckling of Yamaoka et al. (1986; 1987). Lithosphere buckling (deformation) provides the needed mechanism to explain segmentation phenomenon; offsets in volcanic fronts, distribution of calderas within segments, variable segment stresses and the chemical diversity seen between segment boundary and segment interior magmas.

Description: The experiments performed with the inflatable falling sphere technique, for middle atmosphere studies, are reported. It is shown to be a potentially high accurate and independent source of temperature measurement and an intrinsic method for establishing accuracy of other atmospheric measurement techniques. Theoretical derivation, simulations, and actual measurements show that the sphere's temperature data are accurate. It is demonstrated that retrieved temperatures from falling spheres are not significantly affected by linear bias in density caused by uncertainties in sphere mass, volume, or cross sectional area. Case studies illustrate the sphere's capability to produce accurate temperatures. Comparisons with Datasonde temperature measurements obtained close in time and space are in agreement below 60 km.

Description: The simplest model for the Matachewan-Hearst Dike (MHD) magmas is assimilation-fractional crystallization (AFC), presumably occurring at the base of the crust during underplating. Subduction zone enriched mantle sources are not required. Trace elements suggest that the mantle sources for the MHD were depleted, but possessed a degree of heterogeneity. Rates of assimilation were approximately 0.5 (= Ma/Mc); the contaminant mass was less than 20 percent. The contaminant was dominated by tonalites-randodiorites, similar to xenoliths and rocks in the Kapuskasing Structural Zone (KSZ). Assimilation of partial melts of light-rare earth and garnet-bearing basaltic precursors may have produced some the MHD magmas. Apparently, previous underplating-AFC processes had already produced a thick crust. The silicic granitoid assimilant for the MHD magmas was probably produced by earlier processing of underplated mafic crust (4, 5, 10, 21 and 30). Calculations suggest that the derived silicic rocks possess negative Ta and Ti anomalies even though they were not the product of subduction.

Description: The European ARISTOTELES mission aims at the determination of the Earth's gravity field at short wavelength with a global coverage. Gravity gradient measurements will be achieved during six months by the GRADIO instrument onboard a dedicated satellite in a near dawn-dusk sun-synchronous orbit at an altitude of 200 km. The objective is an accuracy of better than 5 mgals for gravity anomalies, at ground level for blocks of 1 x 1 deg. According to present knowledge of the potential, the recovery of higher spherical harmonics (degree and order greater than 30) is of main importance. This leads to focus on the variations of the measured components T(sub ij) of the gravity gradient tensor, at frequencies greater than 5 x 10(exp -3) Hz. The resolution, required for the gradiometer is 10(exp -2) Eotvos (i.e., 10(exp -11)/s squared) with an averaging time of 4 s.

Description: Under contract of the European Space Agency a system study for a spaceborne gravity field recovery mission was performed, covering as a secondary mission objective geodetic point positioning in the cm range as well. It was demonstrated that under the given programmatic constraints including dual launch and a very tight development schedule, a six months gravity field mission in a 200 km near polar, dawn-dusk orbit is adequate to determine gravity anomalies to better than 5 mgal with a spatial resolution of 100 x 100 km half wavelength. This will enable scientists to determine improved spherical harmonic coefficients of the Earth gravity field equation to the order and degree of 180 or better.

Description: The idea of satellite-to-satellite tracking in the high-low mode has received renewed attention in light of the uncertain future of NASA's proposed low-low mission, Geopotential Research Mission (GRM). The principal disadvantage with a high-low system is the increased time interval required to obtain global coverage since the intersatellite visibility is often obscured by Earth. The U.S. Air Force has begun to investigate high-low satellite-to-satellite tracking between the Global Positioning System (GPS) of satellites (high component) and NASA's Space Transportation System (STS), the shuttle (low component). Because the GPS satellites form, or will form, a constellation enabling continuous three-dimensional tracking of a low-altitude orbiter, there will be no data gaps due to lack of intervisibility. Furthermore, all three components of the gravitation vector are estimable at altitude, a given grid of which gives a stronger estimate of gravity on Earth's surface than a similar grid of line-of-sight gravitation components. The proposed Air Force mission is STAGE (Shuttle-GPS Tracking for Anomalous Gravitation Estimation) and is designed for local gravity field determinations since the shuttle will likely not achieve polar orbits. The motivation for STAGE was the feasibility to obtain reasonable accuracies with absolutely minimal cost. Instead of simulating drag-free orbits, STAGE uses direct measurements of the nongravitational forces obtained by an inertial package onboard the shuttle. The sort of accuracies that would be achievable from STAGE vis-a-vis other satellite tracking missions such as GRM and European Space Agency's POPSAT-GRM are analyzed.

Description: A new, geometrical, first order, nonresonant, frozen orbit theory was developed based on Orlov's uniformly rotating plane of constant inclination. Perturbation spectra generated from a 90th order subset of OSU86F are shown for the ill-fated 1984 JHU/APL SAGE proposal for a pair of TRANSIT satellites at 400 km altitude with a 93.5 deg inclination.

Description: The NASA Ocean Topography Experiment satellite TOPEX will carry a microwave altimeter accurate to a few centimeters for the measurement of ocean height. The capability can be fully exploited only if TOPEX altitude can be independently determined to 15 cm or better. This in turn requires an accurate gravity model. The gravity will be tuned with selected nine 10-day arcs of laser ranging, which will be the baseline tracking data type, collected in the first six months of TOPEX flight. TOPEX will also carry onboard an experimental Global Positioning System (GPS) flight receiver capable of simultaneously observing six GPS satellites above its horizon to demonstrate the capability of GPS carrier phase and P-code pseudorange for precise determination of the TOPEX orbit. It was found that subdecimeter orbit accuracy can be achieved with a mere two-hour arc of GPS tracking data, provided that simultaneous measurements are also made at six of more ground tracking sites. The precision GPS data from TOPEX are also valuable for refining the gravity model. An efficient technique is presented for gravity tuning using GPS measurements. Unlike conventional global gravity tuning, this technique solves for far fewer gravity parameters in each filter run. These gravity parameters yield local gravity anomalies which can later be combined with the solutions over other parts of the earth to generate a global gravity map. No supercomputing power will be needed for such combining. The approaches used in this study are described and preliminary results of a covariance analysis presented.

Description: The dominant error source in estimating the orbital position of a satellite from ground based tracking data is the modeling of the Earth's gravity field. The resulting orbit error due to gravity field model errors are predominantly long wavelength in nature. This results in an orbit error signature that is strongly correlated over distances on the size of ocean basins. Anderle and Hoskin (1977) have shown that the orbit error along a given ground track also is correlated to some degree with the orbit error along adjacent ground tracks. This cross track correlation is verified here and is found to be significant out to nearly 1000 kilometers in the case of TOPEX/POSEIDON when using the GEM-T1 gravity model. Finally, it was determined that even the orbit error at points where ascending and descending ground traces cross is somewhat correlated. The implication of these various correlations is that the orbit error due to gravity error is geographically correlated. Such correlations have direct implications when using altimetry to recover oceanographic signals.

Description: Just as the Earth's surface deforms tectonically, so too does the gravity field evolve with time. Now that precise geodesy is yielding observations of these deformations it is important that concomitant, temporal changes in the gravity field be monitored. Although these temporal changes are minute they are observable: changes in the J2 component of the gravity field were inferred from satellite (LAGEOS) tracking data; changes in other components of the gravity field would likely be detected by Geopotential Research Mission (GRM), a proposed but unapproved NASA gravity field mission. Satellite gradiometers were also proposed for high-precision gravity field mapping. Using simple models of geodynamic processes such as viscous postglacial rebound of the solid Earth, great subduction zone earthquakes and seasonal glacial mass fluctuations, we predict temporal changes in gravity gradients at spacecraft altitudes. It was found that these proposed gravity gradient satellite missions should have sensitivities equal to or better than 10(exp -4) E in order to reliably detect these changes. It was also found that satellite altimetry yields little promise of useful detection of time variations in gravity.

Description: Recently, gravity measurements were made on a tall, very stable television transmitting tower in order to detect a non-Newtonian gravitational force. This experiment required the upward continuation of gravity from the Earth's surface to points as high as only 600 m above ground. The upward continuation was based on a set of gravity anomalies in the vicinity of the tower whose data distribution exhibits essential circular symmetry and appropriate radial attenuation. Two methods were applied to perform the upward continuation - least-squares solution of a local harmonic expansion and least-squares collocation. Both methods yield comparable results, and have estimated accuracies on the order of 50 microGal or better (1 microGal = 10(exp -8) m/sq s). This order of accuracy is commensurate with the tower gravity measurments (which have an estimated accuracy of 20 microGal), and enabled a definitive detection of non-Newtonian gravity. As expected, such precise upward continuations require very dense data near the tower. Less expected was the requirement of data (though sparse) up to 220 km away from the tower (in the case that only an ellipsoidal reference gravity is applied).

Description: Integrated geodesy is a method in which a wide variety of surveying measurements are modeled in terms of geometric positions and the earth's geopotential. Using heterogeneous data, both geometric and gravimetric quantities are simultaneously estimated by a least-squares procedure. Heretofore, geodetic leveling differences were reduced into pseudo-observables using assumed values of gravity prior to their inclusion into integrated geodesy least-squares adjustments. The errors in estimates of geometric and gravimetric quantities obtained from integrated geodesy adjustments of geodetic leveling difference, potential differences and Helmert height differences are compared.

Description: For the northern part of the Federal Republic of Germany, new quasigeoid solutions were computed by least squares collocation and fast Fourier transformation (FFT) techniques using point and mean gravity data, a digital terrain model, and a global geopotential model. As severe accuracy limitations for precise regional quasigeoid determination come from global model uncertainties, different geopotential models were investigated by combining them with gravimetric data and comparing the quasigeoid heights with GPS and leveling. Optimum results were obtained by a global model tailored to gravity data in Europe. Collocation and FFT results based on this model agree well. The comparison with GPS and leveling yields rms discrepancies of + or - 2 cm over approximately 400 km range.

Description: The definition of geoidal undulations is given and after summarizing the methods of determination of the geoid, computed geoidal undulations by some of the methods for several points in Kenya are compared to the results obtained by the satellite gravimetric solutions. Results from astrogeodetic levelling and satellite altimetry show some reasonable agreement with the satellite gravimetric geoids while results by Doppler satellite positioning indicate that good agreement can be obtained if the orthometric heights for the points are adjusted to a uniform system.

Description: Using Stokes' formula, a gravimetric geoid was calculated for Canada. The input data are as follows: 15 x 15' block averages were used for Canada and the USA and 1 x 1 deg block averages and satellite model (GEM-T1) provided values for the remaining part of the Earth. The geoid was calculated at 6398 points covering the area within the points rho(sub i)(phi sub i; lambda sub i) (lambda is + west): rho sub 1(40,125); rho sub 2(75,184); rho sub 3(75,10); and rho sub 4(40,60). The computed geoid refers to the GRS1980 and reaches a local minimum of -47.3 meters around the western part of Hudson Bay. A contour map of the geoid is shown.

Description: Model computations show that changes of sampling interval introduce only 0.3 cm changes, whereas zero padding provides an improvement of more than 5 cm in the fast Fourier transformation (FFT) generated geoid. For the Global Positioning System (GPS) survey of Franklin County, Ohio, the parameters selected as a result of model computations, allow large reduction in local data requirements while still retaining the cm accuracy when tapering and padding is applied. The results are shown in tables.

Description: Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Description: Fast Fourier transform (FFT) methods provide a fast and efficient means of processing large amounts of gravity or geoid data in local gravity field modelling. The FFT methods, however, has a number of theoretical and practical limitations, especially the use of flat-earth approximation, and the requirements for gridded data. In spite of this the method often yields excellent results in practice when compared to other more rigorous (and computationally expensive) methods, such as least-squares collocation. The good performance of the FFT methods illustrate that the theoretical approximations are offset by the capability of taking into account more data in larger areas, especially important for geoid predictions. For best results good data gridding algorithms are essential. In practice truncated collocation approaches may be used. For large areas at high latitudes the gridding must be done using suitable map projections such as UTM, to avoid trivial errors caused by the meridian convergence. The FFT methods are compared to ground truth data in New Mexico (xi, eta from delta g), Scandinavia (N from delta g, the geoid fits to 15 cm over 2000 km), and areas of the Atlantic (delta g from satellite altimetry using Wiener filtering). In all cases the FFT methods yields results comparable or superior to other methods.

Description: Different types of gravity prediction methods for local and regional gravity evaluation are developed, tested, and compared. Four different test areas were particularly selected in view of different prediction requirements. Also different parts of the spectrum of the gravity field were considered.

Description: One of the main purposes of geodesy is to determine the gravity field of the Earth in the space outside its physical surface. This purpose can be pursued without any particular knowledge of the internal density even if the exact shape of the physical surface of the Earth is not known, though this seems to entangle the two domains, as it was in the old Stoke's theory before the appearance of Molodensky's approach. Nevertheless, even when large, dense and homogeneous data sets are available, it was always recognized that subtracting from the gravity field the effect of the outer layer of the masses (topographic effect) yields a much smoother field. This is obviously more important when a sparse data set is bad so that any smoothing of the gravity field helps in interpolating between the data without raising the modeling error, this approach is generally followed because it has become very cheap in terms of computing time since the appearance of spectral techniques. The mathematical description of the Inverse Gravimetric Problem (IGP) is dominated mainly by two principles, which in loose terms can be formulated as follows: the knowledge of the external gravity field determines mainly the lateral variations of the density; and the deeper the density anomaly giving rise to a gravity anomaly, the more improperly posed is the problem of recovering the former from the latter. The statistical relation between rho and n (and its inverse) is also investigated in its general form, proving that degree cross-covariances have to be introduced to describe the behavior of rho. The problem of the simultaneous estimate of a spherical anomalous potential and of the external, topographic masses is addressed criticizing the choice of the mixed collection approach.

Description: The most widely known modern method for estimating gravity field values from observed data is least-squares collocation. Its advantages are that it can make estimates at arbitrary locations based on irregularly spaced observations, and that it makes use of statistical information about errors in the input data while providing corresponding information about the quality of the output estimates. Disadvantages of collocation include the necessity of inverting square matrices of dimension equal to the number of data values and the need to assume covariance models for the gravity field and the data errors. Fourier methods are an important alternative to collocation; having the advantage of greater computational efficiency, but requiring data estimates to be on a regular grid and not using or providing statistical accuracy information. The GEOFAST algorithm is an implementation of collocation that achieves high computational efficiency by transforming the estimation equations into the frequency domain where an accurate approximation may be made to reduce the workload. The forward and inverse Fast Fourier Transforms (FFTs) are utilized. The accuracy and computational efficiency of the GEOFAST algorithm is demonstrated using two sets of synthetic gravity data: marine gravity for an ocean trench region including wavelengths longer than 200 km; and local land gravity containing wavelengths as short as 5 km. These results are discussed along with issues such as the advantages of first removing reference field models before carrying out the estimation algorithm.

Description: Masses associated with the topography, bathymetry, and its isostatic compensation are a dominant source of gravity field variations, especially at shorter wavelengths. On global scales the topographic/isostatic effects are also significant, except for the lowest harmonics. In practice, though, global effects need not be taken into account as such effects are included in the coefficients of the geopotential reference fields. On local scales, the short-wavelength gravity variations due to the topography may, in rugged terrain, be an order of magnitude larger than other effects. In such cases, explicit or implicit terrain reduction procedures are mandatory in order to obtain good prediction results. Such effects may be computed by space-domain integration or by fast Fourier transformation (FFT) methods. Numerical examples are given for areas of the Canadian Rockies. In principle, good knowledge of the topographic densities is required to produce the smoothest residual field. Densities may be determined from sample measurements or by gravimetric means, but both are somewhat troublesome methods in practice. The use of a standard density, e.g., 2.67 g/cu cm, may often yield satisfactory results and may be put within a consistent theoretical framework. The independence of density assumptions is the key point of the classical Molodensky approach to the geodetic boundary value problem. The Molodensky solutions take into account that land gravity field observations are done on a non-level surface. Molodensky's problem may be solved by integral expansions or more effective FFT methods, but the solution should not be intermixed with the use of terrain reductions. The methods are actually complimentary and may both be required in order to obtain the smoothest possible signal, least prone to aliasing and other effects coming from sparse data coverage, typical of rugged topography.

Description: Growing interest in high precision studies of the Earth's gravitational field warrant a critical review of precision requirements to yield useful results. Several problems are now under consideration. All of these problems involve, more or less, the precise value of the vertical gradients of gravity. The principle conclusion from this review is that the essential absence of Free Air Vertical Gravity Gradient control and actual values of gravimeter calibrations require serious attention. Large errors in high topography on official published gravity maps also cannot be ignored.

Description: A balloon-borne experiment, flown at 30 km altitude over New Mexico, was used to test dynamic differential Global Positioning System (GPS) tracking in support of gravimetry at high-altitudes. The experiment package contained a gravimeter (Vibrating String Accelerometer), a full complement of inertial instruments, a TI-4100 GPS receiver and a radar transponder. The flight was supported by two GPS receivers on the ground near the flight path. From the 8 hour flight, about a forty minute period was selected for analysis. Differential GPS phase measurements were used to estimate changes in position over the sample time interval, or average velocity. In addition to average velocity, differential positions and numerical averages of acceleration were obtained in three components. Gravitational acceleration was estimated by correcting for accelerations due to translational motion, ignoring all rotational effects.

Description: One of the six absolute gravity instruments developed and built by the Joint Institute for Laboratory Astrophysics (JILA) between 1982 and 1985 was tested under a variety of environmental conditions between May 1987 and 1988. Of the 30 sites visited during this period, 10 were occupied more than once. These reobservations indicate repeatability between 1 and 4 microgals.

Description: The Rapid Geodetic Survey System (RGSS) is a system employing a high-accuracy gimbaled inertial platform. It provides a cost-effective capability for accurate direct measurement of the change in position, elevation, gravity intensity and deflection of the vertical from an initial point. The RGSS is an adaptation of the production version of the U.S. Army Position and Azimuth Determining System (PADS). Several hardware and software enhancements to improve the performance of the system, primarily for gravity vector survey, have occurred over the last few years. The basic principles for the control of error in the survey measurements due to noise and systematic error are discussed below. Actual acceptance test results for the RGSS which indicate an inherent capability of the system to measure change in the deflection of the vertical to a few-tenths of an arcsecond over survey periods of one to two hours using careful survey techniques are also presented. Finally a simple method to extend the capability of the system for longer duration surveys is indicated.

Description: To assess the resolution of the local gravity field from kinematic measurements, a state model for motion in the gravity field of the earth is formulated. The resulting set of equations can accommodate gravity gradients, specific force, acceleration, velocity and position as input data and can take into account approximation errors as well as sensor errors.

Description: At the turn of the century, only pendulum apparatuses and torsion balances were available for general exploration work. Both of these early techniques were cumbersome and time-consuming. It was no wonder that the development of the gravity meter was welcomed with a universal sigh of relief. By 1935 potential field measurements with gravity meters supplanted gradient measurements with torsion balances. Potential field measurements are generally characterized by three types: absolute - measurements are made in fundamental units, traceable to national standards of length and time at each observation site; relative with absolute scale - differences in gravity are measured in fundamental units traceable to national standards of length and time; and relative - differences in gravity are measured with arbitrary scale. Improvements in the design of gravity meters since their introduction has led to a significant reduction in size and greatly increased precision. As the precision increased, applications expanded to include the measurement of crustal motion, the search for non-Newtonian forces, archeology, and civil engineering. Apart from enhancements to the astatic gravity meter, few developments in hardware were achieved. One of these was the vibrating string gravity meter which was developed in the 1950s and was employed briefly for marine and borehole applications. Another is the cryogenic gravity meter which utilizes the stability of superconducting current to achieve a relative instrument with extremely low drift suitable for tidal and secular gravity measurements. An advance in performing measurements from a moving platform was achieved with the development of the straight-line gravity meter. The latter part of the century also saw the rebirth of gradient measurements which offers advantages for observations from a moving platform. Definitive testing of the Bell gradiometer was recently reported.

Description: Although the amount of data yielded by the initial airborne and surface tests was modest, it was sufficient to demonstrate that the full gravity gradient tensor was successfully measured from moving platforms both in the air and on the surface. The measurements were effectively continuous with spatial along-track resolution limited only by choice of integration lengths taken to reduce noise. The airborne data were less noisy (800 E squared/Hz typical) than were the Gravity Gradiometer Survey System (GGSS) measurements taken at the surface (5000 E squared/Hz typical). Single tracks of surface gravity disturbances recovered from airborne data were accurate to 3 to 4 mgal in each component of gravity when compared to 5 x 5 mean gravity anomalies over a 90 km track. Multitrack processing yielded 2 to 3 mgal when compared to 5 x 5 mean anomalies. Deflection of the vertical recovery over a distance of 150 km was about one arcsecond.

Description: New integral formulas to determine geopotential coefficients from terrestrial gravity and satellite altimetry data are given. The formulas are based on the integration of data over the non-spherical surface of the Earth. The effect of the topography to low degrees and orders of coefficients is estimated numerically. Formulas for the solution of the gravimetric boundary value problem are derived.

Description: Traditionally, the global gravity field was described by representations based on the spherical harmonics (SH) expansion of the geopotential. The SH expansion coefficients were determined by fitting the Earth's gravity data as measured by many different methods including the use of artificial satellites. As gravity data have accumulated with increasingly better accuracies, more of the higher order SH expansion coefficients were determined. The SH representation is useful for describing the gravity field exterior to the Earth but is theoretically invalid on the Earth's surface and in the Earth's interior. A new global Earth model (NEWGEM) (KIM, 1987 and 1988a) was recently proposed to provide a unified description of the Earth's gravity field inside, on, and outside the Earth's surface using the Earth's mass density profile as deduced from seismic studies, elevation and bathymetric information, and local and global gravity data. Using NEWGEM, it is possible to determine the constraints on the mass distribution of the Earth imposed by gravity, topography, and seismic data. NEWGEM is useful in investigating a variety of geophysical phenomena. It is currently being utilized to develop a geophysical interpretation of Kaula's rule. The zeroth order NEWGEM is being used to numerically integrate spherical harmonic expansion coefficients and simultaneously determine the contribution of each layer in the model to a given coefficient. The numerically determined SH expansion coefficients are also being used to test the validity of SH expansions at the surface of the Earth by comparing the resulting SH expansion gravity model with exact calculations of the gravity at the Earth's surface.

Description: An important goal of geodesy is to determine the anomalous potential and its derivatives outside of the earth. Representing the surface anomalies by a series of spherical harmonics is useful since it is then possible to do a term by term solution of Laplace's equation and upward continuation. The problem of finding such a spherical harmonic series for anomaly values given on an equiangular surface grid is addressed. (This is a first step toward the more complicated problem of finding a function such that locally averaged values fit a grid of mean anomalies.) Three approaches to this fitting problem are discussed and compared: the discrete Fourier technique, the discrete integral technique, and a new approach. The peculiar nature of the equiangular grid, with its increasing density of (noisy) data toward the poles, causes each method to exhibit a different type of difficulty. The new method is shown to be practical as well as precise since the numerical conditioning problems which appear can be successfully handled by such well-known techniques as a (simple) Kalman filter.

Description: For the determination of geopotential coefficients, data can be used from rather different sources, e.g., satellite tracking, gravimetry, or altimetry. As each data type is particularly sensitive to certain wavelengths of the spherical harmonic coefficients it is of essential importance how they are treated in a combination solution. For example the longer wavelengths are well described by the coefficients of a model derived by satellite tracking, while other observation types such as gravity anomalies, delta g, and geoid heights, N, from altimetry contain only poor information for these long wavelengths. Therefore, the lower coefficients of the satellite model should be treated as being superior in the combination. In the combination a new method is presented which turns out to be highly suitable for this purpose due to its great flexibility combined with robustness.

Description: The ability of three high-order geopotential models (OSU81, GPM2 and OSU86E) to recover the gravity anomaly field (delta g) in the Australian region was tested. The region was divided into 2 x 2 deg blocks, and the mean and rms of the residual gravity (delta g measured - delta g modeled) was found to estimate the fit of the model to the point gravity data. The results showed that OSU81 and GPM2 performed similarly, recovering the delta g with a mean value of less than plus or minus 5 mGal in 63 and 70 percent of the blocks, respectively. However, both these models achieved a fit of worse that was plus or minus 13 mGal in 6 to 7 percent of cases. These were in areas either on or near the coast, or in the Central Australian region, inferring that for a precise geoid slope determination in these regions, a detailed analysis of delta g in region is needed. On the other hand, OSU86E produced a very good result, having a mean fit of less than plus or minus 5 mGal in 80 percent of the blocks, and worse than plus or minus 13 mGal in only 1 percent of cases. The rms values for this model were also improved over the other two models, indicating that for applications requiring highest precision, the preferred model is OSU86E.

Description: Major new computations of terrestrial gravitational field models were performed by the Geodynamics Branch of Goddard Space Flight Center (GSFC). This development has incorporated the present state of the art results in satellite geodesy and have relied upon a more consistent set of reference constants than was heretofore utilized in GSFC's GEM models. The solutions are complete in spherical harmonic coefficients out to degree 50 for the gravity field parameters. These models include adjustment for a subset of 66 ocean tidal coefficients for the long wavelength components of 12 major ocean tides. This tidal adjustment was made in the presence of 550 other fixed ocean tidal terms representing 32 major and minor ocean tides and the Wahr frequency dependent solid earth tidal model. In addition 5-day averaged values for Earth rotation and polar motion were derived for the time period of 1980 onward. Two types of models were computed. These are satellite only models relying exclusively on tracking data and combination models which have incorporated satellite altimetry and surface gravity data. The satellite observational data base consists of over 1100 orbital arcs of data on 31 satellites. A large percentage of these observations were provided by third generation laser stations (less than 5 cm). A calibration of the model accuracy of the GEM-T2 satellite only solution indicated that it was a significant improvement over previous models based solely upon tracking data. The rms geoid error for this field is 110 cm to degree and order 36. This is a major advancement over GEM-T1 whose errors were estimated to be 160 cm. An error propagation using the covariances of the GEM-T2 model for the TOPEX radial orbit component indicates that the rms radial errors are expected to be 12 cm. The combination solution, PGS-3337, is a preliminary effort leading to the development of GEM-T3. PGS-3337 has incorporated global sets of surface gravity data and the Seasat altimetry to produce a model complete to (50,50). A solution for the dynamic ocean topography to degree and order 10 was included as part of this adjustment.

Description: All particles in the atmosphere are not spherical. Moreover, the scattering properties of randomly oriented nonspherical particles are not equivalent to those of spherical particles no matter how the term equivalent is defined. This is especially true for scattering in the backward direction and at the infrared wavelengths at which some atmospheric particles have strong absorption bands. Thus calculations based on Mie theory of infrared backscattering by dry or insoluble atmospheric particles are suspect. To support this assertion, it was noted that peaks in laboratory-measured infrared backscattering spectra show appreciable shifts compared with those calculated using Mie theory. One example is ammonium sulfate. Some success was had in modeling backscattering spectra of ammonium sulfate particles using a simple statistical theory called the continuous distribution of ellipsoids (CDE) theory. In this theory, the scattering properties of an ensemble are calculated. Recently a modified version of this theory was applied to measured spectra of scattering by kaolin particles. The particles were platelike, so the probability distribution of ellipsoidal shapes was chosen to reflect this. As with ammonium sulfate, the wavelength of measured peak backscattering is shifted longward of that predicted by Mie theory.

Description: During the first year two areas of work were emphasized, analysis of aerosol data to provide improved estimates of backscatter over the Pacific ocean and development of a global model for backscatter at different wavelengths and preparation and planning for the aircraft flights of the GLOVE program. Significant progress was made in each of these areas. The analytical work is a continuation of the GAMETAG analysis, and was directed toward the development of a backscatter model for the Pacific Oceanic free troposphere. This has included a further evaluation of relationships between the optical effects at different wavelengths, a comparison of modeled optical effects at differing temporal resolutions, and an investigation of the effects of sampling on the modeled results. An initial investigation of the predictability of 9 to 10 micrometers Beta values from other data sets as well as some preliminary comparisons of the modeling results with experimental data. The comparison of the optical effects was extended to make the comparisons for different sets of time bases between 1 and 20 minutes.

Description: The focus of this effort is the development of a global-scale model of aerosol backscatter for laser atmospheric wind sounder (LAWS) design and performance studies. Background parameters are derived from aerosol data sets with global-scale spatial and/or temporal coverage, using objective statistical decomposition and/or a priori stratification based on supplementary data. Backscatter coefficients at the LAWS design wavelength are derived from background aerosol physical, chemical, and optical data, or from direct backscatter measurements at other wavelengths, using background conversion factors. Direct measurements of aerosol backscatter at 10.6 microns from the Royal Signals and Radar Establishment (RSRE) and the Wave Propagation Laboratory (WPL) were selected. The RSRE backscatter data processing code were optimized under low backscatter conditions, performed detailed analyses of collocated intercomparisons between the two lidars, and assisted in the analysis of the long-term backscatter climatologies from the two lidars. Timely presentation of global backscattering experiment (GLOBE) research results to the global geophysical community is required.

Description: The use of near real-time optical techniques is emphasized for the measurement of mid-tropospheric aerosol over the Central Pacific. The primary focus is on measurement of the aerosol size distribution over the range of particle diameters from 0.15 to 5.0 microns that are essential for modeling CO2 backscatter values in support of the laser atmospheric wind sounder (LAWS) program. The measurement system employs a LAS-X (Laser Aerosol Spectrometer-PMS, Boulder, CO) with a custom 256 channel pulse height analyzer and software for detailed measurement and analysis of aerosol size distributions. A thermal preheater system (Thermo Optic Aerosol Descriminator (TOAD) conditions the aerosol in a manner that allows the discrimination of the size distribution of individual aerosol components such as sulfuric acid, sulfates and refractory species. This allows assessment of the relative contribution of each component to the BCO2 signal. This is necessary since the different components have different sources, exhibit independent variability and provide different BCO2 signals for a given mass and particle size. Field activities involve experiments designed to examine both temporal and spatial variability of these aerosol components from ground based and aircraft platforms.

Description: The research is comprised of the following tasks: use of simple analytical and numerical models of a coupled troposphere-stratosphere system to examine the effects of radiation and ozone on planetary wave dynamics and the tropospheric circulation; use of satellite data obtained from the Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) instrument and Solar Backscattered Ultraviolet (SBUV) experiment, in conjunction with National Meteorological Center (NMC) data, to determine the planetary wave vertical structures, dominant wave spectra, ozone spectra, and time variations in diabatic heating rate; and synthesis of the modeling and observational results to provide a better understanding of the effects that stratospheric processes have on tropospheric dynamics.

Description: The continuous low-g environment of the orbiting space shuttle provided a setting for conducting geophysical fluid model experiments with a completely consistent representation of sphericity and the resultant radial gravity found on astrogeophysical objects. This is possible because in zero gravity an experiment can be constructed that has its own radial buoyancy forces. The dielectric forces in a liquid, which are linearly dependent on fluid temperature, give rise to an effectively radial buoyancy force when a radial electrostatic field is applied. The Geophysical Fluid Flow Cell (GFFC) experiment is an implementation of this idea in which fluid is contained between two rotating hemispheres that are differentially heated and stressed with a large ac voltage. The GFFC flew on Spacelab 3 in May 1985. Data in the form of global Schlieren images of convective patterns were obtained for a large variety of configurations. These included situations of rapid rotation (large Taylor numbers), low rotation, large and small thermal forcing, and situations with applied meridional temperature gradients. The analysis and interpretation of the GFFC-85 data are being conducted. Improvements were developed to the GFFC instrument that will allow for real-time (TV) display of convection data and for near-real-time interactive experiments. These experiments, on the transition to global turbulence, the breakdown of rapidly rotating convective planforms and other phenomena, are scheduled to be carried out on the International Microgravity Laboratory (IML-1) aboard the shuttle in June 1990.

Description: During the 1987 to 1988 academic year, three projects were finished and plans were made to redirect and focus work in a proposal now being reviewed. The completed work involves study of waves on an equatorial beta-plane in shear flow, investigation of the influence of orography on the index cycle, and analysis of a model of cloud street development in a thermally-forced, sheared environment. The proposed work involves study of boundary layer circulations supporting stratocumulus decks and investigation of how the radiative effects of these clouds modulate larger-scale flows such as those associated with the index oscillation.

Description: A subset of data representing most of the participating model groups was discussed and reviewed. The results are presented for intercomparisons of that subset of data. Subsequent to the workshop, the various modeling groups had the opportunity to either update or add to the data sets discussed at the workshop. The overall revised data sets as of December 1988 are presented. Data from all the models (where available) are grouped according to: (1) photochemistry and radiation (photodissociation coefficients, and ultraviolet heating and infrared cooling); (2) transport (net radiative heating, tropospheric source tracer experiment, time-dependent source conserved tracer experiment, stratospheric source tracer experiment, and O3 column); (3) current atmosphere (integrated columns of trace gases, Cl(y) and NO(y), nitrogen gases, chlorine gases, O(x) and HO(x) gases, and source gases); and (4) perturbed circulations and temperatures. The data presented provide a detailed summary of the two-dimensional picture of the atmosphere as seen by current atmospheric models.

Description: The modeling groups are listed along with a brief description of the respective models.

Description: In the attempts to understand the processes which affect ozone in the stratosphere, and to predict future ozone levels, a variety of computer models were developed and employed. To oversimplify somewhat, the models range from 1-D models with intricately detailed photochemistry but only the crudest representation of transport, to 3-D general circulation models with intricately detailed dynamics and no photochemistry at all. Each of these has its appropriate uses. It was felt that there was much to be gained from intermediate models incorporating an extensive treatment of photochemistry with a dynamical framework which at least recognizes that atmospheric motions are advective as well as diffusive, and that both chemistry and dynamics are subject to latitudinal and seasonal variations. Thus the 2-D, zonally averaged models have begun to play a larger role in the last seven years, both in attempting to understand observed distributions of trace species and in attempting to assess the probable effects of antropogenic perturbations. There are many choices to be made in developing the model. It was the purpose of the 2-D model Intercomparison Workshop to permit discussion of the choices made and the behavior of the resulting models. The goal was not to identify a best set of choices, but rather to identify areas in which the models are sensitive to the choices made, and to develop a sense of where these models as a class do well or poorly in simulating the observed atmosphere.

Description: Several factors contribute to the errors encountered. With the exception of the line-by-line model, all of the models employ simplifying assumptions that place fundamental limits on their accuracy and range of validity. For example, all 2-D modeling groups use the diffusivity factor approximation. This approximation produces little error in tropospheric H2O and CO2 cooling rates, but can produce significant errors in CO2 and O3 cooling rates at the stratopause. All models suffer from fundamental uncertainties in shapes and strengths of spectral lines. Thermal flux algorithms being used in 2-D tracer tranport models produce cooling rates that differ by as much as 40 percent for the same input model atmosphere. Disagreements of this magnitude are important since the thermal cooling rates must be subtracted from the almost-equal solar heating rates to derive the net radiative heating rates and the 2-D model diabatic circulation. For much of the annual cycle, the net radiative heating rates are comparable in magnitude to the cooling rate differences described. Many of the models underestimate the cooling rates in the middle and lower stratosphere. The consequences of these errors for the net heating rates and the diabatic circulation will depend on their meridional structure, which was not tested here. Other models underestimate the cooling near 1 mbar. Suchs errors pose potential problems for future interactive ozone assessment studies, since they could produce artificially-high temperatures and increased O3 destruction at these levels. These concerns suggest that a great deal of work is needed to improve the performance of thermal cooling rate algorithms used in the 2-D tracer transport models.

Description: The general theory of MST radar observations of gravity wave spectra is developed. This effort extends the previous results to include anisotropy and Doppler effects for the spectra, as well as the consequences for the multibeam configuration. The relationships between the observed one- or two-dimensional spectra for the line-of-sight velocity in the gravity wave spectra are derived. Expressions for cross spectra, as well as covariances between velocities observed on different beams, are computed. Using these results, studies are carried out to show how model parameters for gravity wave spectra can be estimated from the observed quantities. Model parameters include the variance, power law indices, anisotropy parameters, Doppler parameters, mean scale sizes, etc. Cases with different numbers of beams are investigated.