Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random variables”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

Spectrally Simplified Approach for Leveraging Legacy Geostationary Oceanic Observations

The use of multispectral geostationary satellites to study aquatic ecosystems improves the temporal frequency of observations and mitigates cloud obstruction, but no operational capability presently exists for the coastal and inland waters of the United States. The Advanced Baseline Imager (ABI) on the current iteration of the Geostationary Operational Environmental Satellites, termed the R Series (GOES-R), however, provides sub-hourly imagery and the opportunity to overcome this deficit and to leverage a large repository of existing GOES-R aquatic observations. The fulfillment of this opportunity is assessed herein using a spectrally simplified, two-channel aquatic algorithm consistent with ABI wave bands to estimate the diffuse attenuation coefficient for photosynthetically available radiation, K(d)(PAR). First, an in situ ABI dataset was synthesized using a globally representative dataset of above- and in-water radiometric data products. Values of K(d)(PAR) were estimated by fitting the ratio of the shortest and longest visible wave bands from the in situ ABI dataset to coincident, in situ K(d)(PAR) data products. The algorithm was evaluated based on an iterative cross-validation analysis in which 80% of the dataset was randomly partitioned for fitting and the remaining 20% was used for validation. The iteration producing the median coefficient of determination (R2) value (0.88) resulted in a root mean square difference of 0.319 m−1, or 8.5% of the range in the validation dataset. Second, coincident mid-day images of central and southern California from ABI and from the Moderate Resolution Imaging Spectroradiometer (MODIS) were compared using Google Earth Engine (GEE). GEE default ABI reflectance values were adjusted based on a near infrared signal. Matchups between the ABI and MODIS imagery indicated similar spatial variability (R2 = 0.60) between ABI adjusted blue-to-red reflectance ratio values and MODIS default diffuse attenuation coefficient for spectral downward irradiance at 490 nm, K(d)(490), values. This work demonstrates that if an operational capability to provide- ABI aquatic data products was realized, the spectral configuration of ABI would potentially support a sub-hourly, visible aquatic data product that is applicable to water-mass tracing and physical oceanography research.

Advanced Baseline Imager↗

Validation of a Remote Sensing Based Index of Forest Disturbance Using Streamwater Nitrogen Data

Vegetation disturbances are known to alter the functioning of forested ecosystems by contributing to export ('leakage') of dissolved nitrogen (N), typically nitrate-N, from watersheds that can contribute to acidification of acid-sensitive streams, leaching of base cations, and eutrophication of downstream receiving waters. Yet, at a landscape scale, direct evaluation of how disturbance is linked to spatial variability in N leakage is complicated by the fact that disturbances operate at different spatial scales, over different timescales, and at different intensities. In this paper we explore whether data from synoptic streamwater surveys conducted in an Appalachian oak-dominated forested river basin in western MD (USA) can be used to test and validate a scalable, synthetic, and integrative forest disturbance index (FDI) derived from Landsat imagery. In particular, we found support for the hypothesis that the interannual variation in spring baseflow total dissolved nitrogen (TDN) and nitrate-N concentrations measured at 35 randomly selected stream stations varied as a linear function of the change in FDI computed for the corresponding set of subwatersheds. Our results demonstrate that the combined effects of forest disturbances can be detected using synoptic water quality data. It appears that careful timing of the synoptic baseflow sampling under comparable phenological and hydrometeorological conditions increased our ability to identify a forest disturbance signal.

Eshleman, Keith N.↗

Quantifying Post-fire Vegetation Recovery on the Colorado Front Range

Forest composition and structure in the Colorado Front Range has been altered by changing wildfire regimes, particularly in response to increased fire severity. Subsequent reductions in post-fire tree regeneration can result in chronic impacts to ecological function and water quality. This project partnered with the U.S. Forest Service to estimate tree canopy recovery from 15 to 21 years after four Colorado Front Range fires that burned between 1996 and 2002— the Bobcat, Buffalo Creek, Hayman, and High Meadows fires—using Landsat 5 Thematic Mapper (TM), Landsat 7 Enhanced Thematic Mapper (ETM+), and Landsat 8 Operational Land Imager (OLI). First, we evaluated relationships between recovery metrics derived from spectral vegetation indices (i.e. post-fire vegetation, normalized burn ratio, percent unrecovered vegetation) and field-collected counts of post-fire tree seedling regeneration. Then, we used high-resolution imagery paired with Landsat timeseries spectral and geomorphometric data to map forest canopy cover 15-21 years post-fire. Finally, we assessed ecological drivers of recovery including topography, climate, soils, and fire severity using Random Forest. While we found that basic vegetation index metrics failed to distinguish field-measured tree regeneration, our data showed that both aspect and burn severity were important drivers of vegetation recovery. North-facing aspects tended to have higher pre-fire Normalized Difference Vegetation Index (NDVI) and thus higher burn severities. However, areas that burned at high severity experienced slower vegetation recovery, regardless of aspect suggesting that high severity fire can mute favorable site conditions on north-facing aspects (i.e. moister, more developed soils). While this work propels our understanding of variables that influence vegetative recovery, vegetation type, forest conversion, and watershed function, more work should focus on what vegetation indices identifying areas of recovery indicate across different forest types and burn severities.

NASA DEVELOP↗

Real-Time Adaptive Drag Minimization Wind Tunnel Investigation of a Flexible Wing with Variable Camber Continuous Trailing Edge Flap System

This paper reports the results of a recently completed real-time adaptive drag minimization wind tunnel investigation of a highly flexible wing wind tunnel model equipped with the Variable Camber Continuous Trailing Flap (VCCTEF) technology at the University of Washington Aeronautical Laboratory (UWAL). The wind tunnel investigation is funded by NASA SBIR Phase II contract with Scientific Systems Company, Inc. (SSCI) and University of Washington (UW) as a subcontractor. The wind tunnel model is a sub-scale Common Research Model (CRM) wing constructed of foam core and fiberglass skin and is aeroelastically scaled to achieve a wing tip deflection of 10% of the wing semi-span which represents a typical wing tip deflection for a modern transport such as Boeing 787. The jig-shape twist of the CRM wing is optimized using a CART3D aero-structural model to achieve the minimum induced drag for the design cruise lift coefficient of 0.5. The wing is equipped with two chord wise cambered segments for each of the six span wise flap sections for a total of 12 individual flap segments that comprise the VCCTEF system. Each of the 12 flap segments is actively controlled by an electric servo-actuator. The real-time adaptive drag optimization strategy includes an on-board aerodynamic model identification, a model excitation, and a real-time drag optimization. The on-board aerodynamic model is constructed parametrically as a function of the angle of attack and flap positions to model the lift and drag coefficients of the wing. The lift coefficient models include a linear model and a second-order model. The drag coefficient models include a quadratic model and a higher-order up to 6th-order model to accurately model the drag coefficient at high angles of attack. The onboard aerodynamic model identification includes a recursive least-squares (RLS) algorithm and a batch least-squares (BLS) algorithm designed to estimate the model parameters. The model excitation method is designed to sample the input set that comprises the angle of attack and the flap positions. Three model excitation methods are developed: random excitation method, sweep method, and iterative angle-of-attack seeking method. The real-time drag optimization includes a generic algorithm developed by SSCI and several optimization methods developed by NASA which include a second-order gradient Newton-Raphson optimization method, an iterative gradient optimization method, a pseudo-inverse optimization method, an analytical optimization method, and an iterative refinement optimization method. The first wind tunnel test entry took place in September 2017. This test revealed major hardware issues and required further redesign of the flap servo mechanisms. The second test entry took place in April 2018. However, the test was not successful due to the issues with the onboard aerodynamic model identification RLS algorithm which incorrectly identified model parameters. This test also provides an experimental comparison study between the VCCTEF and a variable camber discrete trailing edge flap (VCDTEF) without the elastomer transition mechanisms. The experimental result confirms the benefit of the VCCTEF which produces lower drag by 5% than the VCDTEF. The third and final test entry took place in June 2018 after the issues with the RLS algorithm have been identified and corrected. Additional improvements were implemented. These include the BLS algorithm, the iterative angle-of-attack seeking method, the iterative gradient optimization method, and the pseudo-inverse optimization method. The test objectives were successfully demonstrated as the real-time drag optimization identifies several optimal solutions at off-design lift coefficients. The iterative gradient optimization method is found to achieve up to 4.7% drag reduction for the off-design lift coefficient of 0.7. The pseudo-inverse optimization method which does not require the drag coefficient model is found to be quite effective in reducing drag. Up to 9.4% drag reduction for the off-design lift coefficient of 0.7 is achieved with the pseudo-inverse optimization method. The wind tunnel investigation demonstrates the potential of real-time drag optimization technology. Several new capabilities are developed that could enable future adaptive wing technologies for flexible wings equipped with drag control devices such as the VCCTEF.

Drag Optimization↗

Elemental Analyses of Hypervelocity Microparticle Impact Sites on Interplanetary Dust Experiment Sensor Surfaces

The Interplanetary Dust Experiment (IDE) had over 450 electrically active ultra-high purity metal-oxide-silicon impact detectors located on the six primary sides of the Long Duration Exposure Facility (LDEF). Hypervelocity microparticles (approximately 0.2 to approximately 100 micron diameter) that struck the active sensors with enough energy to breakdown the 0.4 or 1.0 micron thick SiO2 insulator layer separating the silicon base (the negative electrode), and the 1000 A thick surface layer of aluminum (the positive electrode) caused electrical discharges that were recorded for the first year of orbit. The high purity Al-SiO2-Si substrates allowed detection of trace (ppm) amounts of hypervelocity impactor residues. After sputtering through a layer of surface contamination, secondary ion mass spectrometry (SIMS) was used to create two-dimensional elemental ion intensity maps of microparticle sites on the IDE sensors. The element intensities in the central craters of the impacts were corrected for relative ion yields and instrumental conditions and then normalized to silicon. The results classification resulted from the particles' origins as 'manmade', 'natural', or 'indeterminate'. The last classification resulted from the presence of too little impactor residue, analytical interference from high background contamination, the lack of information on silicon and aluminum residues, or a combination of these circumstances. Several analytical 'blank' discharges were induced on flight sensors by pressing down on the sensor surface with a pure silicon shard. Analyses of these blank discharges showed that the discharge energy blasts away the layer of surface contamination. Only Si and Al were detected inside the discharge zones, including the central craters, of these features. Thus far, a total of 79 randomly selected microparticle impact sites from the six primary sides of the LDEF were analyzed: 36 from tray C-9 (Leading (ram), or east, side), 18 from tray C-3 (Trailing (wake), or west, side), 12 from tray B-12 (north side), 4 from tray D-6 (south side), 3 from tray H-11 (space end), and 6 from tray G-10 (earth end). Residue from manmade debris was identified in craters on all trays (aluminum oxide particle residues were not detectable on the Al/Si substrates). These results were consistent with the IDE impact record which showed highly variable long term microparticle impact flux rates on the west, space, and Earth sides of the LDEF which could not be ascribed to astronomical variability of micrometeorite density. The IDE record also showed episodic bursts of microparticle impacts on the east, north, and south sides of the satellite, denoting passage through orbital debris clouds or rings.

Simon, C. G.↗

Detection in fixed and random noise in foveal and parafoveal vision explained by template learning

Foveal and parafoveal contrast detection thresholds for Gabor and checkerboard targets were measured in white noise by means of a two-interval forced-choice paradigm. Two white-noise conditions were used: fixed and twin. In the fixed noise condition a single noise sample was presented in both intervals of all the trials. In the twin noise condition the same noise sample was used in the two intervals of a trial, but a new sample was generated for each trial. Fixed noise conditions usually resulted in lower thresholds than twin noise. Template learning models are presented that attribute this advantage of fixed over twin noise either to fixed memory templates' reducing uncertainty by incorporation of the noise or to the introduction, by the learning process itself, of more variability in the twin noise condition. Quantitative predictions of the template learning process show that it contributes to the accelerating nonlinear increase in performance with signal amplitude at low signal-to-noise ratios.

NASA Center ARC↗

Component response to random vibratory motion of the carrier vehicle

Two physical models of component plus supporting substructure are considered. Each model consists of a rigid body attached to a moving base by means of linear springs and viscous dampers. The second model differs from the first in that its dampers are elastically supported. The first model receives the more extensive treatment. Base motion, assumed a random translational motion parallel to a fixed axis, is prescribed only to the extent that the power spectral density (PSD) of its acceleration is given; and, as given, its plot on log-log graph paper is a series of straight line segments, each segment having an extremity in common with the adjacent segment. Closed expressions are given for the mean squares of base acceleration, base velocity, and base displacement. The component is restricted to planar motion and allowed two degrees of freedom, one translational and one rotational. Integral expressions are given for the mean squares of component response variables, the transfer functions essential to mean square computation being available via the equations of motion. Closed expressions are given for mean squares of certain of the response variables for the case wherein the base acceleration PSD is constant.

Tuell, L. P.↗

Elemental analyses of hypervelocity microparticle impact sites on Interplanetary Dust Experiment sensor surfaces

The Interplanetary Dust Experiment (IDE) had over 450 electrically active ultra-high purity metal-oxide-silicon impact detectors located on the six primary sides of the Long Duration Exposure Facility (LDEF). Hypervelocity microparticles (approximately 0.2 to approximately 100 micron diameter) that struck the active sensors with enough energy to break down the 0.4 or 1.0 micron thick SIO2 insulator layer separating the silicon base (the negative electrode), and the 1000 A thick surface layer of aluminum (the positive electrode) caused electrical discharges that were recorded for the first year of orbit. The high purity Al-SiO2-Si substrates allowed detection of trace (ppm) amounts of hypervelocity impactor residues. After sputtering through a layer of surface contamination, secondary ion mass spectrometry (SIMS) was used to create two-dimensional elemental ion intensity maps of microparticle impact sites on the IDE sensors. The element intensities in the central craters of the impacts were corrected for relative ion yields and instrumental conditions and then normalized to silicon. The results were used to classify the particles' origins as 'manmade,' 'natural,' or 'indeterminate.' The last classification resulted from the presence of too little impactor residue, analytical interference from high background contamination, the lack of information on silicon and aluminum residues, or a combination of these circumstances. Several analytical 'blank' discharges were induced on flight sensors by pressing down on the sensor surface with a pure silicon shard. Analyses of these blank discharges showed that the discharge energy blasts away the layer of surface contamination. Only Si and Al were detected inside the discharge zones, including the central craters of these features. Thus far a total of 79 randomly selected microparticle impact sites from the six primary sides of the LDEF have been analyzed: 36 from tray C-9 (Leading (ram), or East, side), 18 from tray C-3 (Trailing (wake), or West, side), 12 from tray B-12 (North side), 4 from tray D-6 (South side), 3 from tray H-11 (Space end), and 6 from tray G-10 (Earth end). Residue from manmade debris was identified in craters on all trays. (Aluminum oxide particle residues were not detectable on the Al/Si substrates.) These results were consistent with the IDE impact record which showed highly variable long term microparticle impact flux rates on the West, Space and Earth sides of the LDEF which could not be ascribed to astronomical variability of micrometeorite density. The IDE record also showed episodic bursts of microparticle impacts on the East, North, and South sides of the satellite, denoting passage through orbital debris clouds or rings.

Simon, Charles G.↗

Sound production due to large-scale coherent structures

The acoustic pressure fluctuations due to large-scale finite amplitude disturbances in a free turbulent shear flow are calculated. The flow is decomposed into three component scales; the mean motion, the large-scale wave-like disturbance, and the small-scale random turbulence. The effect of the large-scale structure on the flow is isolated by applying both a spatial and phase average on the governing differential equations and by initially taking the small-scale turbulence to be in energetic equilibrium with the mean flow. The subsequent temporal evolution of the flow is computed from global energetic rate equations for the different component scales. Lighthill's theory is then applied to the region with the flowfield as the source and an observer located outside the flowfield in a region of uniform velocity. Since the time history of all flow variables is known, a minimum of simplifying assumptions for the Lighthill stress tensor is required, including no far-field approximations. A phase average is used to isolate the pressure fluctuations due to the large-scale structure, and also to isolate the dynamic process responsible. Variation of mean square pressure with distance from the source is computed to determine the acoustic far-field location and decay rate, and, in addition, spectra at various acoustic field locations are computed and analyzed. Also included are the effects of varying the growth and decay of the large-scale disturbance on the sound produced.

Gatski, T. B.↗

Wind-forced variability of the deep eastern north Pacific: Observations of seafloor pressure and abyssal currents

Data from an array of bottom pressure gauges and a string of current meters in the vicinity of 47 deg N, 139 deg W, are used to examine the deep-ocean variability forced by ocean surface wind stress curl from August 1987 to June 1988. Bottom geostrophic currents are computed from the pressure gauge array, and these correspond well to the long-period directly measured currents at 3000 m. The supratidal-period bottom pressure variations are coherent at 95% confidence with the wind stress curl in period bands of 3-4 days and 15-60 days but removed in distances of 400 and 700 km to the northwest and the southeast, respectively. A linear, two-layer hydrodynamic model is used to examine the theoretical forcing produced by random-phased surface wind fields for the conditions of the eastern north Pacific and the 15- to 60-day-period observed response is reproduced credibly. To model 3- to 15-day variations, more realistic models are required.

Niller, P. P.↗

Radiative transfer theory for inhomogeneous media with random extinction and scattering coefficients

The small-angle scattering approximation of the scalar radiative transfer equation (RTE) is examined for the case where the extinction and scattering coefficients have a component that is a deterministic function of position along the propagation path and a component that is a random function of position transverse to the propagation direction. It is found that the resulting stochastic RTE can be reduced to a system of two stochastic integrodifferential equations for the average and fluctuating components of the radiant intensity. Two transfer equations are obtained describing the average radiant intensity and the spatial correlation function of the intensity fluctuations. The average intensity equation is then solved and applied to a simple propagation scenario. An approximate solution is also derived for the equation giving the correlation function. The developed equations can be applied to problems involving short wavelength electromagnetic wave propagation through media possessing the variable characteristics of turbulence and turbidity, such as plasmas, the atmosphere, and the ocean.

Manning, Robert M.↗

A model for simulating random atmospheres as a function of latitude, season, and time

An empirical stochastic computer model was developed with the capability of generating random thermodynamic profiles of the atmosphere below an altitude of 99 km which are characteristic of any given season, latitude, and time of day. Samples of temperature, density, and pressure profiles generated by the model are statistically similar to measured profiles in a data base of over 6000 rocket and high-altitude atmospheric soundings; that is, means and standard deviations of modeled profiles and their vertical gradients are in close agreement with data. Model-generated samples can be used for Monte Carlo simulations of aircraft or spacecraft trajectories to predict or account for the effects on a vehicle's performance of atmospheric variability. Other potential uses for the model are in simulating pollutant dispersion patterns, variations in sound propagation, and other phenomena which are dependent on atmospheric properties, and in developing data-reduction software for satellite monitoring systems.

Campbell, J. W.↗

Invasion in the Niger Delta: Remote Sensing of Mangrove Conversion to Invasive Nypa fruticans from 2015-2020

Invasive species are a leading threat to biodiversity worldwide. Nypa palm ( Nypa fruticans ) has emerged as the predominant invasive species in the Niger Delta region of Nigeria. While endemic mangroves have high rates of carbon sequestration, stabilize coastlines, and protect biodiversity, Nypa does not provide these services outside its native region of Southeast Asia. Oil exploration and urbanization in this region also exacerbates mangrove loss and Nypa spread. As Nypa is difficult to distinguish from endemic mangrove species in remotely sensed data, estimates of mangrove and ecosystem services losses in Nigeria are highly uncertain. Here, we analyze multisensor satellite data with machine learning to quantify the rapid expansion of Nypa from 2015-2020 in Nigeria. Using Landsat imagery and random forest classification, we quantify total potential Nypa extent in Nigeria in 2019. We then produced a Nypa extent map using iterative combinations of Sentinel-1 SAR, Sentinel-2 MSI, and ALOS PALSAR. Random forest classifications using SAR data from ALOS and Sentinel-1 were best suited for mapping Nypa extent with similar accuracies (78% and 75% respectively). Based on data availability and accuracy, we focused our change analysis on Sentinel-1 SAR. Our results show ~28,000 ha of mangroves were converted to Nypa in Nigeria by 2020 and covered a larger extent than endemic mangroves, compounding the effect of the existing degradation and deforestation in the region. We also compared forest height and complexity estimates from GEDI (Global Ecosystem Dynamics Investigation) LiDAR to further distinguish between endemic mangroves and Nypa in three dimensions. Nypa structural variability, measured by top-of-canopy height, vegetation cover, plant area index, and foliage height diversity, was lower than that of mangroves. At current rates of Nypa expansion, the entire area of study would be invaded by Nypa by 2028, with potentially detrimental consequences to the ecosystem services provided by mangroves.

GEE↗

In Vitro Selections with RNAs of Variable Length Converge on a Robust Catalytic Core

In vitro selection is a powerful tool that can be used to understand basic principles of molecular evolution. We used in vitro selection to understand how changes in length and the accumulation of point mutations enable the evolution of functional RNAs. Using RNA populations of various lengths, we performed a series of in vitro experiments to select for ribozymes with RNA ligase activity. We identified a core ribozyme structure that was robust to changes in RNA length, high levels of mutagenesis, and increased selection pressure. Elaboration on this core structure resulted in improved activity which we show is consistent with a larger trend among functional RNAs in which increasing motif size can lead to an exponential improvement in fitness. We conclude that elaboration on conserved core structures is a preferred mechanism in RNA evolution. This conclusion, drawn from selections of RNAs from random sequences, is consistent with proposed evolutionary histories of specific biological RNAs. More generally, our results indicate that modern RNA structures can be used to infer ancestral structures. Our observations also suggest a mechanism by which structural outcomes of early RNA evolution would be largely reproducible even though RNA fitness landscapes consist of disconnected clusters of functional sequences.

In vitro evolution↗

A Metric to Quantify Shared Visual Attention in Two-Person Teams

Introduction: Critical tasks in high-risk environments are often performed by teams, the members of which must work together efficiently. In some situations, the team members may have to work together to solve a particular problem, while in others it may be better for them to divide the work into separate tasks that can be completed in parallel. We hypothesize that these two team strategies can be differentiated on the basis of shared visual attention, measured by gaze tracking. 2) Methods: Gaze recordings were obtained for two-person flight crews flying a high-fidelity simulator (Gontar, Hoermann, 2014). Gaze was categorized with respect to 12 areas of interest (AOIs). We used these data to construct time series of 12 dimensional vectors, with each vector component representing one of the AOIs. At each time step, each vector component was set to 0, except for the one corresponding to the currently fixated AOI, which was set to 1. This time series could then be averaged in time, with the averaging window time (t) as a variable parameter. For example, when we average with a t of one minute, each vector component represents the proportion of time that the corresponding AOI was fixated within the corresponding one minute interval. We then computed the Pearson product-moment correlation coefficient between the gaze proportion vectors for each of the two crew members, at each point in time, resulting in a signal representing the time-varying correlation between gaze behaviors. We determined criteria for concluding correlated gaze behavior using two methods: first, a permutation test was applied to the subjects' data. When one crew member's gaze proportion vector is correlated with a random time sample from the other crewmember's data, a distribution of correlation values is obtained that differs markedly from the distribution obtained from temporally aligned samples. In addition to validating that the gaze tracker was functioning reasonably well, this also allows us to compute probabilities of coordinated behavior for each value of the correlation. As an alternative, we also tabulated distributions of correlation coefficients for synthetic data sets, in which the behavior was modeled as a first-order Markov process, and compared correlation distributions for identical processes with those for disparate processes, allowing us to choose criteria and estimate error rates. 3) Discussion: Our method of gaze correlation is able to measure shared visual attention, and can distinguish between activities involving different instruments. We plan to analyze whether pilots strategies of sharing visual attention can predict performance. Possible measurements of performance include expert ratings from instructors, fuel consumption, total task time, and failure rate. While developed for two-person crews, our approach can be applied to larger groups, using intra-class correlation coefficients instead of the Pearson product-moment correlation.

Gontar, Patrick↗

Composite Loads Spectra for Select Space Propulsion Structural Components: Probabilistic Load Model Development

This effort is in support of the development of the expert system of computer codes to predict the loads on select structural components of a space propulsion engine. The development will be based primarily on the space shuttle main engine (SSME) test data base. Because of random variations of the many different sources of the loadings on the selected structural components and transients, a probabilistic approach to the problems was adopted. The goal of this task is to characterize all of the individual sources of loads at critical structural locations, such as the turbine blades, the transfer ducts, and liquid oxygen posts, using state-of-the-art probabilistic methods with varying levels of sophistication. The second phase of this work is the development of a composite load model based on a probabilistic synthesis of the individual load model previously developed. This model will be based on the stochastic combination of the load variables and not on the physical process for the combination of the individual loads to the composite load seen by the selected structural components.

Kurth, R.↗

Long term monitoring of methane in the atmosphere by multiplex gas chromatography

Methane is of interest in the study of the Earth's atmosphere because of its implication in the future global warming of the surface. This warming is produced by the absorption of infrared energy by trace gases. It has been estimated that in the next 40 to 50 years, methane could contribute 20 to 25 pct. as much atmospheric warming as that expected from carbon dioxide increases. Studies to examine sources, sinks, and cycles of methane will require analytical methods capable of continuous unattended measurement with temporal resolution of an hour or less for weeks at a time. Gas chromatography (GC) is one of the most practical methods available to conduct the analysis of air, but limitations in this technique still exist which can be alleviated with multiplex GC (MGC). MGC is a technique where many samples are pseudo-randomly introduced to the chromatograph without regard to the length of time required for an analysis. The resulting data must then be reduced using computational methods such as cross correlation. In the technique reported, a tube packed with silver oxide was used at the inlet of the GC column to create concentration pulses of methane in a sample stream of air. By using only one carrier, i.e., ambient air, an effective and accurate method to monitor the variations in concentration of methane in the atmosphere over long periods of time was developed. Methane in ambient air was monitored for an eight day period and an interesting temporal variability was found. This work has shown the utility of a relatively simple MGC for the analysis of a real environmental sample.

Valentin, Jose R.↗

Physically-based parameterization of spatially variable soil and vegetation using satellite multispectral data

A stochastic-geometric landsurface reflectance model is formulated and tested for the parameterization of spatially variable vegetation and soil at subpixel scales using satellite multispectral images without ground truth. Landscapes are conceptualized as 3-D Lambertian reflecting surfaces consisting of plant canopies, represented by solid geometric figures, superposed on a flat soil background. A computer simulation program is developed to investigate image characteristics at various spatial aggregations representative of satellite observational scales, or pixels. The evolution of the shape and structure of the red-infrared space, or scattergram, of typical semivegetated scenes is investigated by sequentially introducing model variables into the simulation. The analytical moments of the total pixel reflectance, including the mean, variance, spatial covariance, and cross-spectral covariance, are derived in terms of the moments of the individual fractional cover and reflectance components. The moments are applied to the solution of the inverse problem: The estimation of subpixel landscape properties on a pixel-by-pixel basis, given only one multispectral image and limited assumptions on the structure of the landscape. The landsurface reflectance model and inversion technique are tested using actual aerial radiometric data collected over regularly spaced pecan trees, and using both aerial and LANDSAT Thematic Mapper data obtained over discontinuous, randomly spaced conifer canopies in a natural forested watershed. Different amounts of solar backscattered diffuse radiation are assumed and the sensitivity of the estimated landsurface parameters to those amounts is examined.

Jasinski, Michael F.↗