ALBERT

Description: In many structures the final stress states are dependent on the sequence of construction or the stress states at various stages of construction are of interest. Such problems can be analyzed using finite element programs that have the capability of adding (birthing) elements to simulate the progress of construction. However, the usual procedure of assembling elements may lead to numerical instabilities or stress states that are unrealistic. Both problems are demonstrated in the analysis of a structure using the program ADINA. A technique which combines application of a preload with element birthing to overcome these problems is described and illustrated.

Description: A technique utilizing finite element analysis, liquid impact kinematics, and momentum theory is described and compared to single-drop impact test data performed on various configurations of coated ceramic material. The method correlates well with test data and is useful in predicting the single-drop impact damage velocity threshold for low-density, coated ceramic materials.

Description: An algorithm is outlined for simulating the contact surface erosion for impact problems. The algorithm dynamically relocates the contact surface as projectile and target materials exceed their failure criterion. Example computations of axisymmetric and oblique impacts are compared with experimental data.

Description: A study of the dynamic characteristics of a coupled translational-rotational system is given. The formulation of the problem considers the soil-structure interaction effects by utilizing the impedance functions at the foundation of a structure. Due to the fact that the coefficient matrix in the characteristic equation is frequency dependent in nature, iterations have to be performed to find the nature frequencies of the system. Examples and discussions are presented. Comparisons of the analytical results from various approaches are also given.

Description: Two analytical techniques applicable to large deflection dynamic response calculations for pressure loaded composite sandwich panels are demonstrated. One technique utilizes finite element modeling with a single equivalent layer representing the face sheets and core. The other technique utilizes the modal analysis computer code DEPROP which was recently modified to include transverse shear deformation in a core layer. The example problem consists of a simply supported rectangular sandwich panel. Included are comparisons of linear and nonlinear static response calculations, in addition to dynamic response calculations.

Description: The design of thin shell structures with respect to elastoplastic buckling requires an extended analysis of the influence of initial imperfections. For conservative design, the most critical defect should be assumed with the maximum allowable magnitude. This defect is closely related to the initial postbuckling behavior. An algorithm is given for the quasi-static analysis of the postbuckling behavior of structures that exhibit multiple buckling points. the algorithm based upon an energy criterion allows the computation of the critical perturbation which will be employed for the definition of the critical defect. For computational efficiency, the algorithm uses the reduced basis technique with automatic update of the modal basis. The method is applied to the axisymmetric buckling of cylindrical shells under axial compression, and conclusions are given for future research.

Description: A general approach is required for describing matter of behavior when the failure is likely to involve growth and coalescence of a large number of fractures. Failures of this kind appear frequently in rapid dynamic processes, particularly in the formation of spall fragments. An approach to formulating constitutive relations that accounts for the opening, shear and growth of an ensemble of cracks is discussed. The approach accounts for plastic flow accompanying fragmentation. The resulting constitutive relations were incorporated into a Lagrangian computer program. A theoretical approach to coalescence is described. The simplest formulation uses a linear Liouville equation, with crack growth limited by the mean free path of cracks, assumed constant. This approach allows for an anisotropic distribution of cracks. An alternative approach in which the decrease of the mean free path with increasing crack size is accounted for, but the crack distribution is assumed isotropic is described. A reduction of the governing Liouville equation to an ordinary differential equation of third order is possible, and the result can be used to determine how mean free path decreases with increasing crack size.

Description: The design of a composite panel requires some way of finding the minimum thickness laminate which will withstand the load requirements without failure. The mathematical complexity of this problem dictates the use of nonlinear optimization techniques. Specialized laminate optimization programs were developed which are compact and efficient enough to run on microcomputers. Only stresses at a point and inplane loads and deflections are considered. The programs are simple to use and require no knowledge of optimization. Techniques are developed which find minimum thickness laminates with either ply ratios or ply angles as design variables. A method is presented for finding the optimum orientation for the axis of symmetry of an orthotropic laminate. The orthotropic laminate program uses an approximate failure theory, which speed up computations dramatically.

Description: The evaluation of the total probability of a plastic collapse failure P sub f for a highly redundant structure of random interdependent plastic moments acted on by random interdepedent loads is a difficult and computationally very costly process. The evaluation of reasonable bounds to this probability requires the use of second moment algebra which involves man statistical parameters. A computer program which selects the best strategy for minimizing the interval between upper and lower bounds of P sub f is now in its final stage of development. The relative importance of various uncertainties involved in the computational process on the resulting bounds of P sub f, sensitivity is analyzed. Response sensitivities for both mode and system reliability of an ideal plastic portal frame are shown.

Description: A finite element modeling technique which utilizes a triangular element called TSHEL with 45 degrees of freedom and seven point integration was tested for analysis of thin plate and shell structures. The element formulation is based on the degenerate solid shell concept and the mixed formulation with assumed independent inplane and transverse shear strains. The effectiveness of the present modeling technique which features combined use of elements with kinematic modes and those without kinematic modes to eliminate both locking and spurious kinematic modes at the global structural levels are shown.

Description: The stress derivative technique for control of keystoning deformation in under-integrated finte elements is based on expansion of the stress in a Taylor series about the element center and retention of additional terms beyond the constant stress term. It has the advantage over other control techniques that keystoning resistance is provided by actual rather than artificial material properties. Application of this technique to the quadrilateral ring elements used for modelling solids of revolution subjected to axisymmetric loads is described. In a cylindrical coordinate system additional terms appear in the formulation which must be dealt with in arriving at a workable keystoning control scheme.

Description: An existing program is currently being adapted to perform finite element analysis by distributing substructures over a network of four Apple IIe microcomputers connected to a shared disk. In this network, one microcomputer controls the entire process while the others perform the analysis on each substructure in parallel. This substructure analysis is used in an iterative, fully stressed, structural resizing procedure. This procedure allows experimentatation with resizing in which all analyses are not completed during a single iteration. This research gives some insight on how to configure multidiscriplinary analysis and optimization procedures for decomposable engineering systems using either high performance engineering workstations or a parallel processor supercomputer. In addition, the operational experience gained facilitates the implementation of analysis programs on these new computers when they become available in an engineering environment.

Description: Described are the experiences gained from solving for the dynamic response of two simple structures on an experimental Multiple Instruction Multiple Data (MIMD) computer called the finite element machine. Introduced are MIMD computing concepts, describing how the concurrent algorithmic techniques implemented and giving results for the two example problems. The results show computational speedups of up to 7.83 using eight of the finite element machine processors and indicate that significant computational speedups are possible for large order structural computations.

Description: The FLEX/32 Multicomputer is a generic environment for cooperating multiple processors. The FLEX/32 supports a number of different processors, making it heterogeneous in terms of the instruction sets it supports, and homogeneous in its ability to provide consistent storage and input/output facilities to its differing processors. These facilities are accessed through standard 32-bit VMEbus connections. The FLEX/32 supports the full UNIX System V Operating System and languages associated with it, plus the extended ConCurrent C and Concurrent FORTRAN 77 languages that allow programming of concurrent software at a high level. Direct programming support at all levels is provided by the environment hardware for concurrent software execution and optimization, including hardware support for shared resource access arbitration, conditional critical region arbitration, and interprocessor messages.

Description: Research aimed at faster, more cost effective parallel machines and algorithms for improving designer productivity with finite element computations is discussed. A set of 8 boards, containing 4 nearest neighbor connected arrays of commercially available floating point chips and substantial memory, are inserted into a commercially available machine. One-tenth Mflop (64 bit operation) processors provide an 89% efficiency when solving the equations arising in a finite element problem for a single variable regular grid of size 40 by 40 by 40. This is approximately 15 to 20 times faster than a much more expensive machine such as a VAX 11/780 used in double precision. The efficiency falls off as faster or more processors are envisaged because communication times become dominant. A novel successive overrelaxation algorithm which uses cyclic reduction in order to permit data transfer and computation to overlap in time is proposed.

Description: Voyager 1 images show 14 volcanic centers wholly or partly within the Kane Patera quadrangle of Io, which are divided into four major classes: (1) shield with parallel flows; (2) shield with early radial fan shapd flows; (3) shield with radial fan shaped flows, surfaces of flows textured with longitudinal ridges; and (4) depression surrounded by plateau-forming scarp-bounded, untextured deposits. The interpretation attempted here hinges largely on the ability to distinguish lava flows from pyroclastic flows by remote sensing.

Description: Structural and tectonic interpretations of planetary surfaces rely strongly on visual determination of regional structural grain. This grain can be very complex and confusing, and sorting out of discrete trends in time and space is of utmost importance. This study is a test of these techniques applied to a well known area having several discrete structural grains. In the Bighorn Basin region of Wyoming, a largely overlooked N10E structural grain has been verified with detailed structural analysis and indicates a significant change in stress orientation at the end of the Laramide orogeny.

Description: Tectoism in the Valles Marineris appears to have been accompanied by volcanism. The proposed volcanic features, though probably contemporaneous with the gigantic ones in the Tharsis area, are composed of small, mafic and, possibly, somewhat larger felsic flows. The size of these features is similar to that of volcanic flows on the Earth.

Description: Experiments were carried out in a steel pressure device using controlled amounts of water and thermite melt to examine the mechanical energy released on explosive mixing following the initial contact of the two materials. An experimental design was used to allow the direct calculation of the mechanical energy by the dynamic lift of the device as recorded both optically and physically. A large number of experiments were run to accurately determine the optimum mixture of water and melt for the conversion of thermal to mechanical energy. The maximum efficiency observed was about 12% at a water/thermite mass ratio of 0.50. These experiments are the basis for the development of models of hydroexplosions and melt fragmentation. Particles collected from the experimental products are similar in size and shape to pyroclasts produced by much larger hydrovolcanic explosions. Melt rupture at optimum ratios produces very fine particles whereas rupture at high or low water/melt ratios produces large melt fragments. Grain surface textures in the experimental products are also related to the water/melt ratio and the mechanism of explosive mixing. It is thus possible to have qualitative information about the nature of the explosion from the sizes and shapes of the fragments produced.

Description: The design, manufacture, and testing of an isolation system placed under a bench which simulates a space platform during testing in a hard vacuum are described. A low natural frequency isolation system is used to decouple ground vibrations from the bench. The materials used are vacuum compatible and do not introduce unwanted matter which would contaminate the payload optics and/or the testing environment. The system accommodates payloads of varying weights and envelopes.

Description: The potential cost and performance advantages of welding was understood but ignored by solar panel manufacturers in the U.S. Although NASA, DOD and COMSAT have supported welding development efforts, soldering remains the only U.S. space qualified method for interconnecting solar cells. The reason is that no U.S. satellite prime contractor found it necessary, due to mission requirements, to abandon the space proven soldering process. It appears that the proposed NASA space station program will provide an array requirement, a 10 year operation in a low Earth orbital environment, that mandates welding. The status of welding technology in the U.S. is assessed.

Description: Significant new geologic information has been revealed by comparing 1:5 million scale geologic maps of the equatorial zone quadrangles of Mercury (H-6, H-7 and H-8) to Earth-based elevation profiles and surface reflectivity maps of Mercury obtained in the early 1970's at the Arecibo (PR) and Goldstone (CA) radar facilities. These data consist of 23 Goldstone images and profiles of polarized return data at 12.5-cm wavelength and one Arecibo profile. Radar data with 150-m vertical accuracy and 10- to 20-km horizontal resolution are available for areas between latitudes 13 N. and 11 S. In general, these data sets show excellent correlation between: (1) relative elevation and roughness differences that are reflected by mapped geologic contacts; (2) mapped ridges and scarps that display distinctive radar signatures; and (3) position and morphology of crater-and-basin topographic elements. Inferences can also be drawn about topographic and geologic terrain beyond the area imaged by Mariner 10 cameras.

Description: On Mars, the association of gullied escarpments and chaotic terrain is evidence for failure and scarp retreat of poorly consolidated materials. Some martian gullies have no surface outlets and may have drained through subterranean channels. Similar features, though on a much smaller scale, can be seen in alluvium along terrestrial river banks in semiarid regions, such as the Rio Puerco Valley of central New Mexico. Many of the escarpments along the Rio Puerco are developing through formation of collapse gullies, which drain through soil pipes. Gully development can be monitored on aerial photographs taken in 1935, 1962, and 1980. A regression model was developed to quantify gully evolution over a known time span. Soil pipes and their associated collapse gullies make recognizable signatures on the air photos. The areal extent of this signature can be normalized to the scarp length of each pipe-gully system, which makes comparisons between systems possible.

Description: Geomorphic studies of impact structures in central Australia are being used to understand the complexities of fluvial dissection in the heavily cratered terrains of Mars. At Henbury, Northern Territory, approximately 12 small meteorite craters have interacted with a semiarid drainage system. The detailed mapping of the geologic and structural features at Henbury allowed this study to concentrate on degradational landforms. The breaching of crater rims by gullies was facilitated by the northward movement of sheetwash along an extensive pediment surface extending from the Bacon Range. South-facing crater rims have been preferentially breached because gullies on those sides were able to tap the largest amounts of runoff. At crater 6 a probable rim-gully system has captured the headward reaches of a pre-impact stream channel. The interactive history of impacts and drainage development is critical to understanding the relationships in the heavily cratered uplands of Mars. Whereas Henbury craters are younger than 4700 yrs. B.P., the Gosses Bluff structure formed about 130 million years ago. The bluff is essentially an etched central peak composed of resistant sandstone units. Fluvial erosion of this structure is also discussed.

Description: Eastern Acidalia Planitia contains a wide variety of terrain types on which are thousands of subkilometer volcanoes. Apparent morphometric variations were previously reported as a function of terrain type for the cones in the Cydonia area and extended to the rest of Acidalia for which high resolution Viking imagery exist. Crater counts are included for the six types of plains identified, density distributions of subkilometer cones found on each type of terrain, and orphometric data by morphological subclass as a function of terrain for more than 1400 cones.

Description: Preliminary mapping shows East Butte to be a single, large cumulo-dome composed dominantly of rhyolite which can be classified into three main groups based on color and structure. The rhyolite of East Butte is aphanitic with phenocrysts of sanidine and quartz which vary from 1 to 5 mm in length. Vesicular reddish black inclusions of basalt up to 10 cm in length, found in all varieties of the East Butte rhyolites are believed to have originated from fragmentation of the basalt walls of the conduit by rhyolitic magma as it was emplaced. Most of the inclusions contain plagioclase phenocrysts. These phenocrysts measure up to 1 to 2 cm in length and have a typical euhedral, tabular habit. A 250-m diameter depression which has the appearance of a crater is located at the top of East Butte. Evidence supporting the fact that the depression is a crater is displayed by three small (3 to 5 m in height) mounds of massive rhyolite which border the depression.

Description: Some of the geological relationships observed in the Mauna Loa sulfur flow may apply in considering volcanic processes on Io. Given the presence of sulfur/sulfur compounds in the eruption plumes and on the surface of Io, it is likely that extensive secondary deposits of sulfur exist, some of which may be of fumarolic origin and analogous to the Mauna Loa deposit. Given the likelihood of silicate volcanism of Io based on the inferred material properties of some flows, and the attendant high temperatures for silicate volcanism, it is likely that the secondary surface deposits of sulfur would have been mobilized without being heated to the high viscosity stage. Mobilized sulfur flows on Io may flow long distances as a result of: (1) low viscosities in the melting range; (2) sustained effusion resulting from continued heating source area; (3) continued remobilization within the flow as a consequence of surges from the source; and (4) extension via lava tubes, or similar conduits through which there is little heat loss. Sulfur flows may form a relatively thin veneer over silicate flows and other surface units, given their fluidity and low mobilization temperature. Active splashing and splattering may spread sulfur over a wider area contributing the bright blooms observed in association with some Ionian flows.

Description: Researchers define damage-tolerant structural systems as those systems which not only have adeqate intact strength to withstand initial failure but also adequate residual strength to minimize the possibility of, and hence the consequences of, further failure. The incorporation of damage tolerance cannot be done in total isolation of the function being required of the system and the costs associated with obtaining improved damage tolerance. The approach, therefore, is to formulate multiple-objective, multi-level decision support problems (DSP), the solutions of which represent a compromise between higher costs and higher damage tolerance. Mulitple-objective decision support problems are easily solved in the linear domain. These formulations, however, include both linear and nonlinear constraints and goals, which in the past, have not been considered due to the resulting complexity. Here, researchers: (1) present a complete discussion and description of decision support problems; (2) identify what further research needs to be done in order to obtain information that is required but not known for solving problems using these models; and (3) identify what needs to be done to implement this prototype method in practice.

Description: Ocean-going vehicles and platforms are among the largest structures in the world and are subjected to relatively harsh conditions of motions and loads. Some of them, such as semi-submersible platforms, are a relatively new type of structure and hence there is no formal, well evolved and established structural design code as there is for more traditional structures. More recently, efforts have also been made to develop a design method of this type for ships and other ocean structures. One of the many advantages of a rationally based design method is versatility; it can be used for structures that have widely differing purposes, measures of merit, shapes and sizes. The purpose is to describe a rationally based design method that has been developed within the field of ocean structures, in order that persons dealing with other types of structure can judge whether and to what extent its various features may be useful for those other types. Also, even though some features may not be applicable they might stimulate some useful ideas.

Description: Perhaps the greatest chance for exhumation, or burial, of a landscape by terrestrial processes exists near the boundaries of the climatic belts. In the Southern Hemisphere, there is comparatively little land area within Budel's zone of extra-tropical valley formation, which contains most of the examples of exhumed topography in the Northern Hemisphere. The only examples of resurrected landforms that occur within Budel's tropical zone are located near the boundary of this zone, where climate may have changed during the Pleistocene. The ages of exhumed landforms sampled are not equally distributed through geologic time. Most of the exhumed features were created either during the Precambrian or the Tertiary periods which are commonly cited as episodes of significant landform development.

Description: The Dry Valleys of Antarctica are examples of polar deserts which are rare geological features on the Earth. Such deserts typically have high salinities associated with their closed-basin waters and on many surficial materials throughout them. In order to examine the possible sources for the salts observed in association with the soils in the Dry Valleys. The chloride and bromide concentrations of the water leachates from 58 soils and core samples were measured. The Cl/Br ratio for seawater is 289 and ratios measured for most of the 58 soils studied (greater than 85% of the soils studied) was larger than the seawater ratio (ratios typically were greater than 1000 and ranged up to 50,000). The enrichment in Cl relative to Br is strong evidence that the alts present within the soils were derived from seawater during ordinary evaporation processes, and not from the deposition of Cl and Br from aerosols or from rock weathering as has often been suggested.

Description: Field studies of terrestrial landforms and the processes that shape them provide new directions to the study of planetary features. Investigations discussed address principally mudflow phenomena and drainage development. At the Valley of 10,000 Smokes (Katmai, AK) and Mount St. Helens, WA, studies of the development of erosional landforms (in particular, drainage) on fresh, new surfaces permitted analysis of the result of competition between geomorphic processes. Of specific interest is the development of stream pattern as a function of the competition between perennial seepage overland flow (from glacial or groundwater sources), ephemeral overland flow (from pluvial or seasonal melt sources), and ephemeral/perennial groundwater sapping, as a function of time since initial resurfacing, material properties, and seasonal/annual environmental conditions.

Description: Desert pavement is a general term describing a surface that typically consists of a thin layer of cm-sized rock fragments set on top of a layer of finer material in which no fragments are found. An understanding of desert pavement is important to planetary geology because they may play a major role in the formation and visibility of various aeolian features such as wind streaks, which are important on Mars and may be important on Venus. A field study was conducted in Amboy, California to determine the formation mechanism of desert pavements. The probable sequence of events for the formation and evolution of a typical desert pavement surface, based on this experiment and the work of others, is as follows. Starting with a layer of surface material consisting of both fine particles and rock fragments, aeolian deflation will rapidly erode the surface until an armored lag is developed, after which aeolian processes become less important. The concentration of fragments then slowly increases as new fragments are brought to the surface from the subsurface and as fragments move downslope by sheet wash. Sheet wash would be responsible for removing very fine particles from the surface and for moving the fragments relative to one another, forming interlocks.

Description: The use of microcomputers in the design of a cable catenary large space antenna system is discussed. The development of a system design capability, data base utilization, systems integration, program structure and logic, and integrated graphics output are discussed.

Description: An exceptionally elegant method for structural optimization with constraints on the static response presented by Shield and Prager is discussed. Their derivation of the optimality condition was facilitated by a reformulation of the structural elasticity equations in terms of what was then a new variational principle, the principle of stationary mutual potential energy. Their optimality condition relates the design variable to an appropriately defined mutual strain energy. An alternative but related approach, based upon the principle of stationary mutual complementary energy, presented by N. C. Haung, is also discussed. The simplicity of these principles lies in the facts that the energy functionals are stationary at the solution to the field equations and that their stationary value is proportional to the quantity to be optimized.

Description: Field studies of terrestrial landforms and the processes that shape them provide new directions to the study of planetary features. These studies, conducted in Iceland and in Antarctica, investigated physical and chemical weathering mechanisms and rates, eolitan processes, mudflow phenomena, drainage development, and catastrophic fluvial and volcanic phenomena. Continuing investigations in Iceland fall in three main catagories: (1) catastrophic floods of the Jokulsa a Fjollum, (2) lahars associated with explosive volcanic eruptions of Askja caldera, and (3) rates of eolian abrasion in cold, volcanic deserts. The ice-free valleys of Antarctica, in particular those in South Victoria Land, have much is common with the surface of Mars. In addition to providing independent support for the application of the Iceland findings to consideration of the martian erosional system, the Antarctic observations also provide analogies to other martian phenomena. For example, a family of sand dunes in Victoria Valley are stabilized by the incorporation of snow as beds.

Description: Development of drainage networks by erosion by emergent groundwater (sapping) is being modelled by a combination of laboratory experiments and theoretical modelling. Miniature drainage networks formed in fine-grained sediments share many morphologic characteristics of Martian and terrestrial networks suspected to be formed by sapping processes. A larger and better instrumented sapping box was constructed to further explore the processes of sapping and the morphology of resulting networks. The experiments to be conducted in the sapping box will investigate the roles of several factors in controlling network morphology. The mechanics of sapping of fine-grained sediments were investigated in experiments in a two-dimensional sapping chamber and through development of a theoretical model. Results of extensive tests on sapping erosion of fine-grained, cohesionless sediment were analyzed with a theoretical model of the mechanics of sapping erosion and transport.

Description: The 18 May 1980 eruption of Mount St. Helens was initiated by a massive rockslide-debris avalanche which completely transformed the upper 25 km of the North Fork Toutle River valley. The debris was generated by one of the largest gravitational mass movements ever recorded on Earth. Moving at an average velocity of 35 m/s, the debris avalanche buried approximately 60 sq km of terrain to an average depth of 45 m with unconsolidated, poorly sorted volcaniclastic material, all within a period of 10 minutes. Where exposed and unaltered by subsequent lahars and pyroclastic flows, the new terrain surface was characterized predominantly by hummocks, closed depressions, and the absence of an identifiable channel network. Following emplacement of the debris avalanche, a complex interrelationship of fluvial and mass wasting processes immediately began operating to return the impacted area to an equilibrium status through the removal of material (potential energy) and re-establishment of graded conditions. In an attempt to chronicle the morphologic evolution of this unique environmental setting, a systematic series of interpretative maps of several selected areas was produced. These maps, which document the rate and character of active geomorphic processes, are discussed.

Description: Pressure ridges are surface features on basaltic lava flows and, as with other surface features, they may be related to the emplacement of a flow and the rheological properties of the lava. Since many ridges are of sufficient size to be detected on high resolution orbital images, an understanding of pressure ridges could provide a means for interpreting volcanic flows on other terrestrial planets. Some proposed formation mechanisms are reviewed and three different types of pressure ridges are identified on the basis of morphology. Type 1 ridges are the most common and are associated with multiple flow unit pahoehoe in which the ridges are embayed by secondary toe fed lava. They tend to be restricted to wider sections or margins of the flow and to be oriented longitudinal to flow direction; however, oblique or transverse orientation is not uncommon. Bulbous squeeze ups are common within cracks and may reflect relative timing of crack formation. The interior structure of type 1 ridges consists of an upper slab section which generally contains columnar joints and a lower massive section with an irregular surface. This basic distinction may mark the thickness of the surface crust when ridge formation was initiated. Type 2 ridges occur in association with type 1 and are very similar with the exception of the secondary squeeze out material. Instead of only filling cracks, the secondary material on these ridges originated from underneath a thin crust and flowed as toes or channels from the top and sides of the ridge. Type 3 ridges have much steeper sides (almost vertical at the top) than the other types. Medial cracks are very wide and the crack walls are convex upward. No squeeze ups are present. The main difference between type 3 and the others may be reflection of viscosity.

Description: A model for the radiative cooling of thermally well mixed lava flows is presented and the relationship between effusion rate and length and area is analyzed. If radiative cooling is the prime mode of heat loss for a lava flow, one should expect to see a stronger correlation between the effusion rate and the plan area of the flow, than between effusion rate and just flow length. Different flows on a single volcano with differing initial temperatures, volatile content, and gross compositions should yield different areas for a given effusion rate. Likewise, a range of slopes for the relationship between effusion rate and flow area should result from comparisons between different volcanoes. As a test of these ideas, available data on the effusion rates, lengths, and areas of Hawaiian and Etnean flow is studied. It was found that: (1) the effusion rate/area correlation was statistically more significant than the correlation between effusion rate and length for four out of the five eruption episodes which met the necessary criteria of more than three individual flows with area, length, and effusion rate independently measured; (2) that there exists a minimum length and area for a given effusion rate, reflecting competition between overall characteristic proportionality between effusion rate and flow length, width, and area.

Description: Between 1217 and 1620 hours (PDT), on May 18, 1980, the magmatic eruption column of Mount St. Helens formed an ash fountain and pyroclastic flows dominated the eruption process over tephra ejection. Eurption-rate pulsations generally increased to a maximum at 1600 to 1700 hrs. After 1620 hrs, the eruption assumed an open-vent discharge with strong, vertical ejection of tephra. Relative eruption rates (relative mass flux rates) of the pyroclastic flows were determined by correlating sequential photographs and SLAR images, obtained during the eruption, with stratigraphy and surface morphology of the deposits.

Description: Reexamination of Upheaval Dome in the Canyonlands National Park, Utah, shows that the structure of this remarkable feature conforms with that expected for a deeply eroded astrobleme. The structure is definitely not compatible with an origin due simply to plastic flowage of salt and other rocks in the underlying Paradox Formation. The most strongly deformed rocks are bounded by a series of circumferential listric faults. The convergent displacement of the rocks corresponds to the deformation that results from collapse of a transient cavity produced by high speed impact. From considerations of the probable depth of exposure of the impact structure and upward extrapolation of the listric faults, the final collapsed crater is estimated to be about 8 to 10 km in diameter; the impacting body was on the order of 0.5 km in diameter.

Description: The effects of space environments on damping materials and damping designs on flexible structures were investigated. The following items were examined: damping of flexible spacecraft appendages; composite loss factor (n sub s) vs. time in high vacuum for damped test beams and damping of flexible structures. The STEP experiments show inherent damping of flexible structures in space effective possible damping design configurations for space structures, effects of passively damped components on the system loss factor of flexible structures and the effect of space environment on properties of damping materials.

Description: A brief synopsis of the low Earth orbit (LEO) satellite environment is presented including neutral and ionic species. Two ground based atomic and molecular beam instruments are described which are capable of simulating the interaction of spacecraft surfaces with the LEO environment and detecting the results of these interactions. The first detects mass spectrometrically low level fluxes of reactively and nonreactively surface scattered species as a function of scattering angle and velocity while the second ultrahigh velocity (UHV) molecular beam, laser induced fluorescence apparatus is capable of measuring chemiluminescence produced by either gas phase or gas-surface interactions. A number of proposed experiments are described.

Description: A brief review of atmospheric composition in low Earth orbit is presented. The flux of ambient atomic oxygen incident on a surface orbiting in this environment is described. Estimates are presented of the fluence of atomic oxygen to which satellite surfaces in various orbits are exposed.

Description: Interplanetary dust particles (micrometeoroids) are expected to form well-defined craters upon impacting exposed material in space. Studying the frequency and features of these craters will provide data on the mass-flux distribution of micrometeoroids and, to a lesser extent, on the velocity magnitude and direction. This experiment will study impact craters produced by micrometeoroids on selected materials (metals and glasses in the form of thick targets) to obtain valuable technological and scientific data. Specifically, the studies will focus on determining micrometeoroid composition and mass-flux distribution. Analyses will also be made on the distribution of impact velocity vectors.

Description: Three-dimensional finite element modeling techniques are used to synthesize geodetic and seismological results for 1979 Imperial Valley earthquake. The strategy pursued consists of two principal steps. In the first step, the seismologically-derived coseismic fault slip is taken as a function of position in the fault plane and is applied directly to the three-dimensional dislocation model. In the second step, a physical model of stresses and constitutive parameters is perturbed so as to reproduce the observed fault slip. Hence, the principal features of the coseismic slip pattern are explained by a stress-driven fault model in which: (1) a spatially unresolved asperity is found equivalent to a stress drop of 18 MPa averaged over an area of 15 sq km, and (2) driving stress is essentially absent on the fault segment overlapping the 1940 earthquake rupture zone.

Description: Atmospheric interactions involving the nitrogen species are varied and complex. These interactions include photochemical reactions, initiated by the absorption of solar photons and chemical kinetic reactions, which involve both homogeneous (gas-to-gas reactions) and heterogeneous (gas-to-particle) reactions. Another important atmospheric interaction is the production of nitrogen oxides by atmospheric lightning. The nitrogen cycle strongly couples the biosphere and atmosphere. Many nitrogen species are produced by biogenic processes. Once in the atmosphere nitrogen oxides are photochemically and chemically transformed to nitrates, which are returned to the biosphere via precipitation, dry deposition and aerosols to close the biosphere-atmosphere nitrogen cycle. The sources, sinks and photochemistry/chemistry of the nitrogen species; atmospheric nitrogen species; souces and sinks of nitrous oxide; sources; sinks and photochemistry/chemistry of ammonia; seasonal variation of the vertical distribution of ammonia in the troposphere; surface and atmospheric sources of the nitrogen species, and seasonal variation of ground level ammonia are summarized.

Description: If all biotic sinks of atmospheric CO2 reported were added a value of about 0.4 Gt C/yr would be found. For each category, a very high (non-conservative) estimate was used. This still does not provide a sufficient basis for achieving a balance between the sources and sinks of atmospheric CO2. The bulk of the discrepancy lies in a combination of errors in the major terms, the greatest being in a combination of errors in the major terms, the greatest being in the net biotic release and ocean uptake segments, but smaller errors or biases may exist in calculations of the rate of atmospheric CO2 increase and total fossil fuel use as well. The reason why biotic sinks are not capable of balancing the CO2 increase via nutrient-matching in the short-term is apparent from a comparison of the stoichiometry of the sources and sinks. The burning of fossil fuels and forest biomass releases much more CO2-carbon than is sequestered as organic carbon.

Description: The initial attempts to model the atmospheric CO2 distribution, including couplings to the ocean and biosphere as sources and sinks of atmospheric CO2, encourage the notion that this approach will lead to useful quantitative constraints on CO2 fluxes. Realization of this objective will require: (1) continued improvement in the realism of the global transport modeling; (2) extended timeline of atmospheric CO2 monitoring, which improved precision and improved definition of the uncertainties in the measured CO2 amounts; and (3) given an accurate knowledge of model capabilities and limitations and given a good understanding of CO2 observations and their limitations, there is a need for good ideas concerning what quantitative information on the carbon cycle can be inferred from global modeling.

Description: A science that chooses the globe as it fundamental biogeophysical unit forces extraordinary conceptual difficulties. The roles of energy flow, matter cycles, carbon cycle, air pollution, global effects, air water interactions are discussed.

Description: Measurements of atmospheric aerosol backscatter coefficients, using a coherent CO2 lidar at 9.25- and 10.6-micron wavelengths, are described. Vertical profiles of the volume backscatter coefficient beta have been measured to a 10-km altitude over the Pasadena, CA, region. These measurements indicate a wide range of variability in beta both in and above the local boundary layer. Certain profiles also indicate a significant enhancement in beta at the 9.25-micron wavelength compared with beta at the 10.6-micron wavelength, which possibly indicates a major contribution to the volume backscatter from ammonium sulfate aerosol particles.

Description: While in the lobes of the distant magnetotail, ISEE-3 encountered regions of compressed magnetic field at a rate of several per day. The duration of these events was 5 to 20 minutes and they were observed 10 to 30 minutes following the onset of substorm activity near the earth. During each event, the lobe magnetic field tilted first northward and then southward with the inflection point near the time of peak field strength. Following the compression events, the lobe field weakened and retained a southward component for 20 to 40 minutes. It is suggested that these traveling compression regions are the lobe signatures of plasmoids moving rapidly down the tail in the plasma sheet. Comparison of ISEE-3 compression event times with substorm onset times yielded propagation speeds of 350 to 750 km/s.

Description: Ionospheric rocket sounding data for eight nighttime auroral events are used to characterize relativistic electron showers and their effects on atmospheric ozone. The rockets were launched from the Poker Flat Research Range in Alaska and from Andoya, Norway over the period 1976-82. Energetic fluxes were always detected but were of insufficient magnitude to produce significant changes in stratospheric ozone. However, middle atmospheric energy sources were found to be dominated by relativistic electrons and X-ray bremmstrahlung, the latter from 40-55 km and the former from 55-60 km altitudes. The ionizing radiation is concluded to be a significant factor in mesospheric ion conductivity, mobility, electric field structure and analytical models for the ion-neutral chemistry.

Description: The dynamics of the polar thermosphere are examined by using observations made from the Dynamics Explorer 2 satellite. The results used in this study were obtained primarily from the Fabry-Perot interferometer (FPI) and the wind and temperature spectrometer (WATS) during the time period from September 1981 through January 1982. Two primary geophysical conditions were examined: these were the southern summer and the northern winter polar regions. The results support the conclusion that above 60 degrees of latitude the neutral winds are strongly controlled by ion/neutral frictional momentum transfer resulting from magnetospheric convection. This implies that the natural coordinate system within which to display the neutral winds in the high polar thermosphere is magnetic. The collected observations of this study were used to assess the validity of two of the large thermospheric general circulation models. The result of this assessment was that the models reasonably represent the vector winds at high altitudes but do not, at present, accurately simulate the thermodynamics of that regime.

Description: A geometrically nonlinear finite element analysis of cohesive failure in typical joints is presented. Cracked-lap-shear joints were chosen for analysis. Results obtained from linear and nonlinear analysis show that nonlinear effects, due to large rotations, significantly affect the calculated mode 1, crack opening, and mode 2, inplane shear, strain-energy-release rates. The ratio of the mode 1 to mode 2 strain-energy-release rates (G1/G2) was found to be strongly affected by the adhesive modulus and the adherend thickness. The ratios between 0.2 and 0.8 can be obtained by varying adherend thickness and using either a single or double cracked-lap-shear specimen configuration. Debond growth rate data, together with the analysis, indicate that mode 1 strain-energy-release rate governs debond growth. Results from the present analysis agree well with experimentally measured joint opening displacements. Previously announced in STAR as N83-13497

Description: Line strengths of N2O and its isotopic derivatives in the 1120-1440/cm region were measured at low pressure and high resolution (0.0054/cm). The band strength, rotationless dipole moment matrix elements, and F factor coefficients were considered. First-order nondegenerate perturbation theory was employed to derive explicit expressions for the rotationless dipole moment matrix elements and F factor coefficients. This made it possible to obtain general expressions for the F factor. The derived expressions were also applicable to CO2 bands.

Description: The atmospheric and oceanic mass transport associated with the southern oscillation/El Nino will inevitably induce an interannual variation in the length of day. An empirical correlation study is conducted by comparing the Southern Oscillation Index time series and the interannual length-of-day variation (obtained by removing the long-period and short-period variations from a BIH length-of-day series) for the period 1957-1983. The two series have an encouraging qualitative correlation, in particular with respect to El Nino events; and the linear correlation coefficient is found to be 0.55. It is believed that much, if not most, interannual length-of-day variation is caused by the southern oscillation, and the true correlation is considerably higher than its apparent value considering the fact that the Southern Oscillation Index is merely an indicator derived from two local atmospheric measurements.

Description: Infrared solar absorption spectra recorded at 0.02-per cm resolution during a balloon flight from Alamogordo, NM (33 deg N), on March 23, 1981, have been analyzed for the possible presence of absorption by formic acid (HCOOH). An absorption feature at 1105 per cm has been tentatively identified in upper tropospheric spectra as due to the nu-6 band Q branch. A preliminary analysis indicates a concentration of about 0.6 ppbv and 0.4 ppbv near 8 and 10 km, respectively.

Description: The UV spectrometer polarimeter on the Solar Maximum Mission has been utilized to measure mesospheric O3 altitude profiles by the technique of solar occultation. Sunset data are presented for 1980, during the fall equinoctal period within + or - 20 deg of the geographic equator. Mean O3 concentrations are (40, 16, 5.5, and 1.5) x 10 to the 9th/cu cm at 50, 55, 60, and 65 km, respectively. Some profiles exhibit altitude structure which is wavelike. The mean O3 profile is fit best with the results of a time-dependent model if the assumed water-vapor mixing ratio employed varies from 6 ppm at 65 km.

Description: A comprehensive formula is worked out for the photochemical time constant of minor constituents in the middle atmosphere. The formula is particularly useful for evaluating the time constants for families of some chemical species that are strongly coupled by rapid exchange reactions. Height profiles of photochemical time constants are calculated for individual species and their families based on the chemical reaction rate constants recommended in the recent WMO and JPL reports. Potential exchange reactions among various family members are discussed, and the effects of the choice of family membership on the time constant are evaluated.

Description: Sixty-four observations of the orbital node made by the Lageos satellite over a five year time interval reveal an acceleration of (-8.1 + or - 1.8) x 10 to the -8 power arcseconds day/2 due to a source which is not presently modeled in the GEODYN orbit determination computer program. This acceleration cannot be explained by the ocean tide with 18.6 year period, assuming it to be an equilibrium tide. Instead it seems to be due to postglacial rebound, which changes the J(2) coefficient in the spherical harmonic expansion of the earth's gravitational field at the rate of (-8.2 + or - 18) x 10 to the -19th power/s; this in turn accelerates the node. This rate does not agree with the -32 x 10 to the -19th power/s predicted by Wu and Peltier's (1982) L2 model, which has upper and lower mantle effective viscosities of 10 to the 21st and 22nd powers Pa's, respectively. It does agree well with their L1 model, which gives about 10 x 10 to the 19th power/s. Since the effective viscosity is 10 to the 21st power Pa s throughout the entire mantle in the L1 model, the results support the contentions that the efective viscosity is near 10 to the 21st power Pa s everyhere in the mantle, and this relatively low value for the effective viscosity may have permitted several degrees of polar wander due to glaciation during the Quaternary Ice Age. Previously announced in STAR as N84-13705

Description: ISEE-3 magnetic-field measurements in the region of the geomagnetic tail from -80 to -220 earth radii are reported and discussed. A well-ordered field structure is found, comprising two 7-8-nT lobes separated by a plasma sheet, an embedded neutral sheet with significant By fields, and an intermittent plasma-sheet boundary layer with 5-nT-amplitude (peak-to-peak) electromagnetic waves. The plasma-sheet Bz distribution changes from principally northern orientation near the earth to an approximately equal north-south distribution at 200-220 earth radii. These findings are considered to be in general agreement with magnetic-reconnection models of the magnetosphere, with reconnection either throughout the region observed (in tearing-mode or plasmoid-formation models) or at a constant (about 220-earth-radii) or variable (40-80 to 220-earth-radii) X line (in X-line models).

Description: Using plasma electron and magnetic field measurements from ISEE 3, 220 earth radii from earth, it is found that the magnetotail at that distance is a coherent structure that evidently waves about through distances comparable to its own lateral scale size. For about one-third of the time it was inside the magnetotail, ISEE 3 was in the plasma sheet. During quiet times the plasma sheet is apparently quite thin, but in response to geomagnetic activity it expands, becoming filled with hot plasma flowing tailward at speeds sometimes exceeding 1000 km/sec, and forces the magnetotail cross-section itself to expand. The plasma sheet's expansion is delayed typically by about 30 minutes from the onset of the associated geomagnetic activity (often a clearly identified isolated substorm). The magnetic field in the newly-expanded plasma sheet usually exhibits a few-minute steep northward excursion followed by a more prolonged (and often steep) southward excursion. These are believed to be the signatures of arrival of a plasmoid formed and released near the earth at the onset of the corresponding geomagnetic activity. The discreteness of these plasma releases through the magnetotail and their close association with onsets of geomagnetic activity at earth suggest that they are consequences of spontaneous release, probably by magnetic reconnection, of energy and plasma earlier stored in the magnetotail.

Description: The effects of solar activity on stratospheric ozone were studied between 1970 and 1976 by the NIMBUS 4 backscattered ultraviolet experiment. Results show that, after correcting for instrument sensitivity, globally averaged ozone decreased by 3-4 percent above 2 mbar to less than 1 percent at 10 mbar, as solar activity decreased from solar maximum to solar minimum. This systematic ozone decrease (at all pressure levels) and the seemingly periodic oscillation (generally a 2 year period) seem to be associated with conventional indices of solar activity, which suggests a solar UV-ozone relationship. However, since the ozone depletion, especially at 40 km, is characteristic of atmospheric fluorocarbon injection effects, the solar cycle ozone relationship should be qualified: it may exist if the solar UV flux varies by 15-20 percent in the shorter-wavelength region (less than 200 nm).

Description: The Cold Arctic Mesopause Program (CAMP) was conducted at ESRANGE, Sweden, in July/August 1982. During the time period of several weeks, the temperature was monitored by ground-based OH emission spectrometers and by satellite radiance measurements. Rocket launchings occurred on the nights of 3/4 and 11/12 August. On 3/4 August, seven rocket payloads were launched during a period of noctilucent cloud sighting over ESRANGE. The presence of the NLC was confirmed by several rocket-borne photometer profiles. The temperature measurements showed that the temperature profiles in the stratosphere and lower mesosphere were near the expected values of high latitude summer models. A large amplitude wave structure with three temperature minima of 139K, 114K and 111K were observed at altitudes between 83 and 94 km. The temperature minimum at 83 km was the location of the observed NLC. The temperature minima caused by the growth of the gravity wave amplitude in the highly stable mesosphere provide the regions for the growth of particles by nucleation to optical scattering size, as well as regions where the nuclei for condensation can be formed through ion chemistry paths.

Description: Ocean floor basalts studied from the Galapagos Ridge, FAMOUS area, Cayman Trough and Kilauea east rift contain 20-200 ppm carbon and 0.3-2.8 ppn nitrogen as sums of the vesicle-filling gases CO2 and N2 and dissolved species. The wide range of carbon contents found is due partly to the different extent of outgassing of vesicle-filling gases and partly to depth dependency of dissolved CO2 in the basalts. Sulfate commonly exists with sulfide in these basalts, and the sulfate/sulfide ratio increases with increasing water content, perhaps reflecting the higher oxidation potential in basalt melt of the higher water content.

Description: Studies have been conducted to develop efficient techniques to simulate crack extension and to examine various local and global fracture criteria. Of the considered criteria, the crack-tip-opening angle (CTOA) or displacement (CTOD) at a specified distance from the crack tip was shown to be most suited for modeling stable crack growth and instability during the fracture process. The results obtained in a number of studies show the necessity for studying different crack configurations when assessing the validity of any fracture criteria. One of the objectives of the present investigation is related to a critical evaluation of the CTOD growth criterion using an elastic-plastic finite element analysis under monotonic loading to failure. The analysis was found to predict three stages of crack growth behavior under monotonic loading to failure. Calculated CTOD values agreed well with experimental values for crack growth initiation.

Description: ISEE-3 measurements indicate that a broad mantle-like boundary layer plasma often exists within the distant geomagnetic tail lobes at all latitudes, directly adjacent to the tail magnetopause. The presence of this boundary layer at large tail distances indicates that plasma from the magnetosheath often crosses the magnetopause locally along much of the length of the tail, and is evidence that the tail is 'open'.

Description: Line parameters for the nu(5) and 2nu(9) bands and associated hot bands of HNO3 have been calculated and compared with laboratory spectra, and the results are presented. Spectral intervals near 870/cm for which best agreement was obtained are used to quantitatively analyze HNO3 absorption features in 0.02/cm resolution stratospheric solar absorption spectra.

Description: Ozone data were obtained in the free troposphere at altitudes of 5 to 7.5 km and at latitudes from 44 degrees N to 46 degrees S during a 3-week period from October 19 to November 4, 1982. Flight trajectories included the continental U.S., Central America, and the western coast of South America. Maximum observed ozone was 110 ppb at about 41 degrees S latitude and at an altitude of 7.3 km (msl). Analysis of the data as a function of latitude showed a region of elevated ozone in the southern latitudes extending from about 2 degrees to 30 degrees in which latitudinal ozone averages (2 degree bands) peaked at about 80 ppb. Ozone concentrations, measured at northern latitudes, were lower, averaging about 35 ppb. A discussion of the data, including comparison with other ozone data sets, is presented.

Description: Total vertical column amounts of NO have been derived from infrared solar absorption spectra recorded near sunrise and sunset with the 0.01 per cm resolution Fourier transform interferometer at the national Solar Observatory on Kitt peak (elevation 2095 m, latitude 31.9 degrees N) on February 23, 1981. The results show an increase in NO concentration in the morning, late afternoon values about 40 percent higher than in the morning, and a decrease in NO concentration prior to sunset. The measured diurnal changes in the total vertical column amount are compared with values obtained from time-dependent photochemical calculations.

Description: Measurements of SO2 emission rates and concentrations and of particle distribution, size, shape, and composition were made in quiescent volcanic plumes emitted into the troposphere from Poas and Arenal volcanos, Costa Rica, and Colima volcano, Mexico. SO2 emission rates were 700 + or - 180 metric tons per day (t/d) for Poas, 210 + or - 30 t/d for Arenal, and 320 + or - 50 t/d for Colima. The concentrations of SO2 calculated from the COSPEC/lidar data were 5-380 ppb. Concentrations of SO2 measured directly by flame photometry were 10-250 ppb. Particles collected in the plumes with a quartz crystal microbalance impactor were mostly less than 3 microns in diameter and consisted of droplets of dilute sulfur-bearing solutions and minor amounts of layer silicate particles coated with a sulfur-bearing film or crust. Total particle concentrations were 4.7 micrograms per cu.m for Poas and 18.8 micrograms per cu.m for Colima. Comparison of concentrations of SO2 in the plumes with gas samples collected at fumaroles on the ground suggests that the plumes are diluted by the atmosphere by factors of up to 100,000.

Description: A two-dimensional photochemical model based on diabatic circulation has been used to simulate the behavior of N2O, CFCl3(F-11) and CF2Cl2(F-12) The circulation is based on estimates of net heating from the ground to 60 km. Eddy diffusion has been reduced with respect to other model studies with Kzz = 2000 sq cm/s everywhere bove 100 mbar. Resulting tracer profiles show reasonable agreement with measured profiles in the tropics and fall off much more sharply with altitude than those produced by models using larger values of Kzz. The agreement obtained is at least as good as that obtained with adjustable, eddy diffusion parameters. The diabatic circulation treatment is more closely related to related to real physical processes and thus more easily interpreted. Diffusive mixing appears to be more important in determining the details of the tracer distributions than the basic morphology.

Description: The mid-latitude upper stratospheric ozone profiles obtained by the solar backscatter ultraviolet instrument on the Nimbus 7 satellite show a clear annual cycle both in the absolute ozone amounts between 0.98 and 15.6 mbar and in the magnitude of disturbances that reveal themselves as longitudinal structure. At the lowest pressure analyzed a winter maximum in ozone exists, but as one progresses downward in altitude a shift in the temporal phase of the annual cycle occurs in the vicinity of 3 to 4 mbar. Comparison of the observed behavior with the predictions of a one-dimensional photochemical model shows a systematic tendency for calculated ozone amount to be 20-27 percent below the data for pressures less than 7.8 mbar. The chemical model successfully predicts the change in phase of the annual cycle, although at a pressure greater than observed. Diagnosis of model results shows the observed shift to be closely coupled to the magnitude of the ozone column density near 3-4 mbar. The wavelength-dependent attenuation of the solar radiation field by ozone alters the relative magnitude of the molecular oxygen and ozone dissociation rates, leading to a change in the temporal phase of the annual cycle.

Description: One-dimensional radiative and photochemical models are used to determine how much CO2 must have been present to maintain a temperate early climate and to examine the consequences that are implied for the controls on atmospheric oxidation state. It is shown that CO2 concentrations of the order of 1000 PAL are required to keep the average surface temperature close to the present value, if albedo changes and heating by reduced greenhouse gases were relatively unimportant. The oxidation state of such a high-CO2, prebiotic atmosphere should have been largely determined by the balance between the H2O2 rainout rate and the rate at which hydrogen escaped to space, with only a weak dependence on the volcanic outgassing rate or on other speculative sources of H2. The implied upper limit on the ground-level O2 mixing ratio is approximately 10 to the -11th and is subject to less uncertainty than the results of previous models.

Description: A time-dependent, one-dimensional photochemical model of the troposphere is used to describe the vertical distribution of atmospheric trace constituents for summer-time conditions at midlatitudes in the Northern Hemisphere. The model incorporates a planetary boundary layer (PBL) parametrization and a detailed chemical mechanism that includes the photochemistry of important nonmethane hydrocarbon species formed during the oxidation process. One result of the parametrized PBL is that the concentrations of some trace species in the free troposphere are 20-30 percent higher than when mixing processes are described by a vertical eddy diffusion coefficient which is held constant with respect to height and time. The lifetime of the oxides of nitrogen against photochemical conversion to nitric acid during summertime conditions is on the order of six hours. This lifetime is short enough to deplete most of the NO(x) in the PBL so that other reactive nitrogen species are more abundant than NO(x) throughout the free troposphere.

Description: A petrologic model of the northern Mississippi Embayment, derived from gravity, seismic and rift data, is evaluated by converting the model to a magnetization model which is compared with satellite magnetic anomaly models. A magnetization contrast of approximately -0.54 A/m, determined from the petrologic model of the embayment compares favorably to values of -0.62 A/m and -0.45 A/m from a Magsat United States Apparent Magnetization Contrast Map and a published POGO magnetization contrast model, respectively. The petrologic model suggests that the magnetic anomaly low associated with the Mississippi Embayment may be largely due to the intrusion under non-oxidizing conditions of low Curie temperature gabbroic material at the base of the crust of the embayment. Near-surface mafic plutons, bordering the Mississippi Valley Graben, appear from aeromagnetic data to have higher magnetizations than the deeper gabbroic material; however, it is impossible to ascertain if this is due to compositional differences or similar material at shallower (lower temperature) depths. These results indicate that variations in the Curie temperatures of intrusions accompanying rifting may account for a large part of the wide range of magnetic anomalies associated with presently inactive rifts with normal heat flow.

Description: A discussion is presented concerning the role of water in its liquid and gaseous states in the chemistry and photochemistry of the earth's troposphere, with attention to the formation of the hydroxyl radical, oxidation chains involving the gases methane, carbon monoxide and ammonia, the atmospheric chemistry of the sulfur, hydrogen, halogen and nitrogen species, and the function of lightning as a source of tropospheric species. The phenomena of 'rainout', 'washout', and the aqueous chemistry of cloud an rain droplets and of water-covered aerosols, are noted. A section is devoted to the past and anticipated impact of anthropogenic activities on the chemistry and composition of the earth's atmosphere.

Description: Attenuation of the load transmitted to a vehicle occupant by the structure, either by modifying structural assembly, changing geometry of its elements, or adding specific load limiting devices to help dissipate kinetic energy was studied. General aviation aircraft have fuselage subfloors of a built up structure which are generally very stiff perpendicular to the floor of the cabin. The subfloor structure, designed to crush at an appropriate force level, can be used to advantage in a crash by dissipating energy through plastic buckling of the floor beams and frames. Simple closed form solutions to predict the mean crushing force levels of subfloor designs is useful in engineering practice. The crushing process of thin walled, plate formed, open structures with particular emphasis on L and cruciform shapes was analyzed. Lower and upper bound solutions for the mean crushing strength of cruciforms are obtained by considering modes of deformation which account for both bending and extensional deformation. It is the importance of extensional deformation to the energy absorption process, representing at least one third of the dissipated energy is shown. Previously announced in STAR as N82-33765

Description: Using the 9 layer general circulation model developed at the Goddard Laboratory for Atmospheric Sciences (GLAS), several 4 to 5 day integrations were made to assess the impact that latent heating processes (supersaturation and moist convective) have on the model forecasts. In an earlier study by Shukla (1981) it was hypothesized that because of strong interaction between dynamics and moist convection, small initial errors grow very fast and make short range forecasting difficult. The purpose of this study was to examine if prescribed heating rates can improve the forecasts for a few days.

Description: The methods for predicting the strength of composite laminates with fastener holes are reviewed for the cases of unloaded as well as loaded holes. Three categories of methods are discussed: (1) the hole-boundary stress (HBS) methods which compare the peak stresses and the laminate strength; (2) the linear elastic fracture (LEFM) methods which assume cracks emanating from the hole and predict failure when the stress intensity factor for a crack equals the laminate fracture toughness; and (3) two different two-parameter (TP) methods: the average-stress (ASTP) method, which predicts failure when the average stress over a distance from the hole equals the laminate strength; and the point-stress (PSTP) Whitney-Nuismer (1974) approach, which predicts failure when the stress near the hole equals the material strength. The PSTP is relatively accurate and is by far the most widely used of all the prediction methods. However, for large ranges of variables and cases, the current strength predictions may not always be accurate.

Description: Optimality criteria methods take advantage of some concepts as those of statically determinate or indeterminate structures, and certain variational principles of structural dynamics, to develop efficient algorithms for the sizing of structures that are subjected to stiffness-related constraints. Some of the methods and iterative strategies developed over the last decade for calculations of the Lagrange multipliers in stressand displacement-limited problems, as well as for satisfying the appropriate optimality criterion, are discussed. The application of these methods are illustrated by solving problems with stress and displacement constraints.

Description: The stiffened composite structural panel analysis and sizing code designated 'PASCO' encompasses both the generality required for the exploitation of composite materials' design flexibility and an accurate buckling analysis for the detection of complex buckling modes. PASCO can accordingly design for buckling, frequency, material strength, and panel stiffness requirements. Attention is given to an additional thermal loading design capability. Design studies illustrate the importance of the multiple load condition capability when thermal loads are present.

Description: The structure of the stratospheric spectrum, measured with a resolution of 0.0033/cm unapodized in the 7-90/cm spectral interval is described. The spectra from measurements recorded on April 24, 1979 from a balloon-borne platform at about 39 km altitude launched from the National Scientific Balloon Facility of Palestine, TX are presented. The main spectroscopic constituents include a few strong and saturated lines due to the rotational spectrum of water vapor, to magnetic dipole transitions of oxygen, and to the rotational spectrum of ozone. It is shown that the lines of ozone prevent the positive identification of the contribution to the submillimeter spectrum of several minor constituents such as HF, HCl, HCN, HNO3, and N2O. The assignments which are possible in the spectra are revised.

Description: Infrared solar spectra have been obtained during a series of flights with a balloon borne Fourier Transform Spectrometer system. Major emphasis during these flights was placed on obtaining data during sunset. The spectra thus contain information on the atmospheric transmission over long stratospheric paths. These spectra have been analyzed to obtain information on the constituents responsible for the observed absorptions and the distribution with altitude of several of these constituents of interest in stratospheric chemistry. These spectra have also been used to determine the spectroscopic parameters for several compounds such as O3 which are difficult to study in the laboratory. A description of the instrumentation used to obtain the data, samples of the spectra obtained, and details of the analysis are given.

Description: A detailed comparison is presented of the predicted eigenfrequencies of twisted rotating plates as obtained by using two different shape functions. Primarily, rotating twisted plates of two different aspect ratios and two different thickness ratios are considered. The effects of rotation are included by using a 'stress smoothing' technique when calculating the augmented stiffness matrix. In addition, the effects of Coriolis acceleration, contributions from membrane behavior, setting angle and sweep angle are considered. The effects of geometric nonlinearity are briefly discussed. Finally, results of a brief study of cambered plates are presented.

Description: Delamination growth in compressively loaded composite laminates was studied analytically and experimentally. The configuration used was a laminate with an across-the-width delamination. An approximate super-position stress analysis was developed to quantify the effects of various geometric, material, and load parameters on mode 1 and mode 2 strain energy release rates G sub 1 and G sub 2, respectively. Calculated values of G sub 1 and G sub 2 were then compared with measured cyclic delamination growth rates to determine the relative importance of G sub 1 and G sub 2. High growth rates were observed only when G sub 1 was large. However, slow growth was observed even when G sub 1 was negligibly small. This growth was apparently due to a large value of G sub 2.

Description: Recent theoretical and observational studies of middle atmosphere electrodynamics are reviewed. Attention is given to observations of large electric fields in the mesosphere and stratosphere which suggest magnitudes of about one volt per sq m. Recommendations are offered with respect to areas of future study, with emphasis on studies of the morphology of large electric fields, and their relationship with external influences such as magnetospheric electric fields and tropospheric thunderstorms.

Description: The steady state gravity wave model of Schoeberl et al. (1983) is extended to compute wave breaking by disturbances originating at the earth's surface. For winter and summer mean zonal wind profiles, no waves reach the mesosphere unless the absolute value of the zonally averaged perturbed zonal velocity minus c is greater than approximately 20 m/sec. Gravity waves with c = 0 can only reach the winter mesosphere if planetary scale waves are present in the troposphere and the lower stratosphere, to provide strong zonal wind channels for upward wave propagation. This results in nonzonal wave breaking in the mesosphere which could provide in situ forcing of planetary waves. Dissipation of gravity waves by molecular viscosity and conduction can provide significant deceleration and heating/cooling in the 85-105 km region.

Description: The three-satellite ERBE will be flown in the mid-1980s and will provide upgraded data on changes in the input and outflow of energy which drives climatic phenomena. The spacecraft will carry radiometers, self-calibration equipment, and be programmed with inversion and averaging algorithms. Radiation will be scanned at 0.2-5 and 5-50 microns. The ACRIM detector, similar to that on the SMM spacecraft, will measure the incoming solar radiation at the 600 km sun-synchronous orbit travelled by the three ERBE spacecraft. Data transmitted to the NASA-Langley Center will be analyzed with algorithms which extract radiant flux data from the measured radiances.

Description: Results from an intercomparison of methods to measure carbon monoxide (CO), nitric oxide (NO), and the hydroxyl radical (OH) are discussed. The intercomparison was conducted at Wallops Island, Virginia, in July 1983 and included a laser differential absorption and three grab sample/gas chromatograph methods for CO, a laser-induced fluorescence (LIF) and two chemiluminescence methods for NO, and two LIF methods and a radiocarbon tracer method for OH. The intercomparison was conducted as a field measurement program involving ambient measurements of CO (150-300 ppbv) and NO (10-180 pptv) from a common manifold with controlled injection of CO in incremental steps from 20 to 500 ppbv and NO in steps from 10 to 220 pptv. Only ambient measurements of OH were made. The agreement between the techniques was on the order of 14 percent for CO and 17 percent for NO. Hardware difficulties during the OH tests resulted in a data base with insufficient data and uncertanties too large to permit a meaningful intercomposition.

Description: The ozone density of the earth's mesosphere in the 1.0-0.1 mbar (48 to 70 km) region has been measured at sunlit latitudes for the period from December 1981 until the present by an ultraviolet spectrometer on the Solar Mesosphere Explorer satellite. Results for 1982 are reported. The ozone mixing ratios are found to be highly variable in time and place, with maxima occurring in the winter hemispheres. The results show complex time variations at all pressure levels, with annual and semiannual variations apparent at most pressures and latitudes. A relative maximum occurs in July at the equator.

Description: In this study it is shown that the upward moving electron beams observed by DE-1 in the dayside polar cusp at 2 to 5 R sub E drive an electrostatic electron acoustic mode, rather than the whistler mode near resonance, as previously concluded. The characteristics of this instability are compared with the properties of hiss in the polar cusp observed by the Dynamics Explorer (DE-1) satellite. The hiss frequencies and funnel shaped frequency-time spectra are consistent with the electron acoustic mode. However, because the hiss often has a fluctuating magnetic field component near the axis of the funnel shape, the hiss is probably composed of both electron acoustic and whistler mode waves.

Description: The stratospheric emission in the submillimeter spectral region has been measured from 38.2 km altitude with a resolution of 0.0033 kayser. In most cases this resolution makes possible the separation of contributions caused by the individual rotational transitions. An atlas of the observed features is given as a guide to both the measurements of minor stratospheric constituents concentration and the evaluation of stratospheric transparency for the measurement of nontelluric sources.

Description: A transfer function approach is taken in constructing a spectral model of the acoustic-gravity wave response in a multiconstituent thermosphere. The model is then applied to describing the thermospheric response to various sources around the globe. Zonal spherical harmonics serve to model the horizontal variations in propagating waves which, when integrated with respect to height, generate a transfer function for a vertical source distribution in the thermosphere. Four wave components are characterized as resonance phenomena and are associated with magnetic activity and ionospheric disturbances. The waves are either trapped or propagate, the latter becoming significant when possessing frequencies above 3 cycles/day. The energy input is distributed by thermospheric winds. The disturbances decay slowly, mainly due to heat conduction and diffusion. Gravity waves appear abruptly and are connected to a sudden switching on or off of a source. Turn off of a source coincides with a reversal of the local atmospheric circulation.

Description: ISEE 1 electron spectrometer observations are used to study the electron distribution function in the low latitude boundary layer (LLBL) of the magnetosphere in the sunward direction. The study included 43 examples of well-defined layers. Using three-dimensional, model-independent distribution functions of electrons observed in the boundary layer, beams of field-aligned electrons with typically 100 eV energy maximums, streaming in either one, or both, directions have been identified. The phase densities in the LLBL and at lower altitudes are similar, and the beams observed in the LLBL are identified with those observed at lower altitudes. They form an important component of the plasma electrons in the LLBL, which otherwise resembles a mixture of magnetosheath-like and magnetospheric electrons.

Description: A special three-dimensional element based on the total Lagrangian description of the motion of a layered anisotropic composite medium is developed, validated and employed to analyze laminated anisotropic composite shells. The element contains the following features: geometric nonlinearity, dynamic (transient) behavior and arbitrary lamination scheme and lamina properties. Numerical results of nonlinear bending, natural vibration, and transient response are presented to illustrate the capabilities of the element.

Description: Two airborne surveys of atmospheric carbon dioxide concentration have been conducted over the Gulf Stream off the east coast of Virginia and North Carolina on September 7-8, 1983. In situ CO2 data were acquired at an aircraft altitude of 300 m on trajectories that transcected the Gulf Stream near 36 deg N 73 deg W. Data show evidence of a CO2 concentration increase by 4 ppm to 15 ppm above the nominal atmospheric background value of 345 ppm. These enhanced values were associated with the physical location of the Gulf Stream prior to the passage of a weak cold front.