Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “synthetic grid data”

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 91 records · Page 5

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms.

Knapmeyer, Martin↗

Multi-Sensor Cloud and Aerosol Retrieval Simulator and Remote Sensing from Model Parameters : Aerosols - Part 2

The Multi-sensor Cloud Retrieval Simulator (MCRS) produces a simulated radiance product from any high-resolution general circulation model with interactive aerosol as if a specific sensor such as the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of the atmospheric column and land ocean surface at a specific location. Previously the MCRS code only included contributions from atmosphere and clouds in its radiance calculations and did not incorporate properties of aerosols. In this paper we added a new aerosol properties module to the MCRS code that allows users to insert a mixture of up to 15 different aerosol species in any of 36 vertical layers. This new MCRS code is now known as MCARS (Multi-sensor Cloud and Aerosol Retrieval Simulator). Inclusion of an aerosol module into MCARS not only allows for extensive, tightly controlled testing of various aspects of satellite operational cloud and aerosol properties retrieval algorithms, but also provides a platform for comparing cloud and aerosol models against satellite measurements. This kind of two-way platform can improve the efficacy of model parameterizations of measured satellite radiances, allowing the assessment of model skill consistently with the retrieval algorithm. The MCARS code provides dynamic controls for appearance of cloud and aerosol layers. Thereby detailed quantitative studies of the impacts of various atmospheric components can be controlled. In this paper we illustrate the operation of MCARS by deriving simulated radiances from various data field output by the Goddard Earth Observing System version 5 (GEOS-5) model. The model aerosol fields are prepared for translation to simulated radiance using the same model sub grid variability parameterizations as are used for cloud and atmospheric properties profiles, namely the ICA technique. After MCARS computes modeled sensor radiances equivalent to their observed counterparts, these radiances are presented as input to operational remote-sensing algorithms. Specifically, the MCARS-computed radiances are input into the processing chain used to produce the MODIS Data Collection 6 aerosol product (MOYD04). TheMOYD04 product is of course normally produced from MOYD021KM MODIS Level-1B radiance product directly acquired by the MODIS instrument. MCARS matches the format and metadata of a MOYD021KM product. The resulting MCARS output can be directly provided to MODAPS (MODIS Adaptive Processing System) as input to various operational atmospheric retrieval algorithms. Thus the operational algorithms can be tested directly without needing to make any software changes to accommodate an alternative input source. We show direct application of this synthetic product in analysis of the performance of the MOD04 operational algorithm. We use biomass-burning case studies over Amazonia employed in a recent Working Group on Numerical Experimentation (WGNE)-sponsored study of aerosol impacts on numerical weather prediction (Freitas et al., 2015). We demonstrate that a known low bias in retrieved MODIS aerosol optical depth appears to be due to a disconnect between actual column relative humidity and the value assumed by the MODIS aerosol product.

aerosol retrieval↗

A Comparative Study of Simulated and Measured Main Landing Gear Noise for Large Civil Transports

Computational results for the NASA 26%-scale model of a six-wheel main landing gear with and without a toboggan-shaped noise reduction fairing are presented. The model is a high-fidelity representation of a Boeing 777-200 aircraft main landing gear. A lattice Boltzmann method was used to simulate the unsteady flow around the model in isolation. The computations were conducted in free-air at a Mach number of 0.17, matching a recent acoustic test of the same gear model in the Virginia Tech Stability Wind Tunnel in its anechoic configuration. Results obtained on a set of grids with successively finer spatial resolution demonstrate the challenge in resolving/capturing the flow field for the smaller components of the gear and their associated interactions, and the resulting effects on the high-frequency segment of the farfield noise spectrum. Farfield noise spectra were computed based on an FWH integral approach, with simulated pressures on the model solid surfaces or flow-field data extracted on a set of permeable surfaces enclosing the model as input. Comparison of these spectra with microphone array measurements obtained in the tunnel indicated that, for the present complex gear model, the permeable surfaces provide a more accurate representation of farfield noise, suggesting that volumetric effects are not negligible. The present study also demonstrates that good agreement between simulated and measured farfield noise can be achieved if consistent post-processing is applied to both physical and synthetic pressure records at array microphone locations.

Konig, Benedikt↗

Compactly‐Supported Nonstationary Kernels for Computing Exact Gaussian Processes on Big Data

The Gaussian process (GP) is a widely used method for analyzing large-scale data sets, including spatio-temporal measurements of nonlinear processes that are now commonplace in the environmental sciences. Traditional implementations of GPs involve stationary kernels (also termed covariance functions) that limit their flexibility, and exact methods for inference that prevent application to data sets with more than about 10,000 points. Modern approaches to address stationarity assumptions generally fail to accommodate large data sets, while all attempts to address scalability focus on approximating the Gaussian likelihood, which can involve subjectivity and lead to inaccuracies. In this work, we explicitly derive an alternative kernel that can discover and encode both sparsity and nonstationarity. We embed the kernel within a fully Bayesian GP model and leverage high-performance computing resources to enable the analysis of massive data sets. We demonstrate the favorable performance of our novel kernel relative to existing exact and approximate GP methods across a variety of synthetic data examples. Furthermore, we conduct space–time prediction based on more than 1 million measurements of daily maximum temperature and verify that our results outperform state-of-the-art methods in the Earth sciences. More broadly, having access to exact GPs that use ultra-scalable, sparsity-discovering, nonstationary kernels allows GP methods to truly compete with a wide variety of machine learning methods.

Gaussian processes↗

Life Cycle Cost Analysis of Prestressed Concrete Poles Subjected to Wind, Surges, and Waves

Prestressed concrete (PC) poles are becoming popular choices to support coastal power transmission systems. However, the existing literature does not offer a detailed analysis of the effectiveness of PC poles in terms of long-term vulnerabilities and the direct and indirect costs. This is due to (1) lack of fragility models for PC poles and (2) lack of probabilistic wind, storm surge, and wave models in coastal settings. In this study, we address these gaps through a series of Monte Carlo simulations to estimate fragility of PC poles as a function of age and hazard (wind, surges, and waves) intensity, and the development of a probabilistic hazard model based on 10,000 years of synthetic tropical cyclone data. The probabilistic hazard model is used in conjunction with high-resolution hydrodynamic models to generate realizations of coastal wind, storm surges, and waves for the Louisiana and Mississippi coasts. A comprehensive life cycle cost analysis for a service life of 70 years considering direct and indirect losses is conducted to compare the performance of a transmission line located in Pascagoula, Mississippi, when wood poles are replaced by PC poles. Results showed that aging has a minor effect on the reliability of PC poles, highlighting the advantages of replacing wood poles with PC poles, especially in coastal areas. In addition, PC poles are significantly more cost-effective compared with wood poles over their life cycle, leading to an estimated saving of $11.55 million (68.17% reduction). The results of this study provide key insight to inform decision-making processes to keep the coastal power grids resilient and cost-effective against future storm hazards.

natural disasters↗

Achieving Accuracy Requirements for Forest Biomass Mapping: A Data Fusion Method for Estimating Forest Biomass and LiDAR Sampling Error with Spaceborne Data

The synergistic use of active and passive remote sensing (i.e., data fusion) demonstrates the ability of spaceborne light detection and ranging (LiDAR), synthetic aperture radar (SAR) and multispectral imagery for achieving the accuracy requirements of a global forest biomass mapping mission. This data fusion approach also provides a means to extend 3D information from discrete spaceborne LiDAR measurements of forest structure across scales much larger than that of the LiDAR footprint. For estimating biomass, these measurements mix a number of errors including those associated with LiDAR footprint sampling over regional - global extents. A general framework for mapping above ground live forest biomass (AGB) with a data fusion approach is presented and verified using data from NASA field campaigns near Howland, ME, USA, to assess AGB and LiDAR sampling errors across a regionally representative landscape. We combined SAR and Landsat-derived optical (passive optical) image data to identify forest patches, and used image and simulated spaceborne LiDAR data to compute AGB and estimate LiDAR sampling error for forest patches and 100m, 250m, 500m, and 1km grid cells. Forest patches were delineated with Landsat-derived data and airborne SAR imagery, and simulated spaceborne LiDAR (SSL) data were derived from orbit and cloud cover simulations and airborne data from NASA's Laser Vegetation Imaging Sensor (L VIS). At both the patch and grid scales, we evaluated differences in AGB estimation and sampling error from the combined use of LiDAR with both SAR and passive optical and with either SAR or passive optical alone. This data fusion approach demonstrates that incorporating forest patches into the AGB mapping framework can provide sub-grid forest information for coarser grid-level AGB reporting, and that combining simulated spaceborne LiDAR with SAR and passive optical data are most useful for estimating AGB when measurements from LiDAR are limited because they minimized forest AGB sampling errors by 15 - 38%. Furthermore, spaceborne global scale accuracy requirements were achieved. At least 80% of the grid cells at 100m, 250m, 500m, and 1km grid levels met AGB density accuracy requirements using a combination of passive optical and SAR along with machine learning methods to predict vegetation structure metrics for forested areas without LiDAR samples. Finally, using either passive optical or SAR, accuracy requirements were met at the 500m and 250m grid level, respectively.

Montesano, P. M.↗

Rectification of terrain induced distortions in radar imagery

This paper describes a technique to generate geocoded synthetic aperture radar (SAR) imagery corrected for terrain induced geometric distortions. This algorithm transforms the raw slant range image, generated by the signal processor, into a map registered product, resampled to either Universal Transverse Mercator (UTM) or Polar Stereographic projections, and corrected for foreshortening. The technique utilizes the space platform trajectory information in conjunction with a digital elevation map (DEM) of the target area to generate an ortho-radar map with near-autonomous operation. The current procedure requires only two to three tie-points to compensate for the platform position uncertainty that results in translational error between the image and the DEM. This approach is unique in that it does not require generation of a simulated radar image from the DEM or a grid of tie-points to characterize the image-to-map distortions. Rather, it models the inherent distortions based on knowledge of the radar data collection characteristics, the signal Doppler parameters, and the local terrain height to automatically predict the registration transformation. This algorithm has been implemented on a minicomputer system equipped with an array processor and a large random-access memory to optimize the throughput.

Kwok, Ronald↗

Forecasting of Storm-Surge Floods Using ADCIRC and Optimized DEMs

Increasing the accuracy of storm-surge flood forecasts is essential for improving preparedness for hurricanes and other severe storms and, in particular, for optimizing evacuation scenarios. An interactive database, developed by WorldWinds, Inc., contains atlases of storm-surge flood levels for the Louisiana/Mississippi gulf coast region. These atlases were developed to improve forecasting of flooding along the coastline and estuaries and in adjacent inland areas. Storm-surge heights depend on a complex interaction of several factors, including: storm size, central minimum pressure, forward speed of motion, bottom topography near the point of landfall, astronomical tides, and, most importantly, maximum wind speed. The information in the atlases was generated in over 100 computational simulations, partly by use of a parallel-processing version of the ADvanced CIRCulation (ADCIRC) model. ADCIRC is a nonlinear computational model of hydrodynamics, developed by the U.S. Army Corps of Engineers and the US Navy, as a family of two- and three-dimensional finite-element-based codes. It affords a capability for simulating tidal circulation and storm-surge propagation over very large computational domains, while simultaneously providing high-resolution output in areas of complex shoreline and bathymetry. The ADCIRC finite-element grid for this project covered the Gulf of Mexico and contiguous basins, extending into the deep Atlantic Ocean with progressively higher resolution approaching the study area. The advantage of using ADCIRC over other storm-surge models, such as SLOSH, is that input conditions can include all or part of wind stress, tides, wave stress, and river discharge, which serve to make the model output more accurate. To keep the computational load manageable, this work was conducted using only the wind stress, calculated by using historical data from Hurricane Camille, as the input condition for the model. Hurricane storm-surge simulations were performed on an eight-node Linux computer cluster. Each node contained dual 2-GHz processors, 2GB of memory, and a 40GB hard drive. The digital elevation model (DEM) for this region was specified using a combination of Navy data (over water), NOAA data (for the coastline), and optimized Interferometric Synthetic Aperture Radar data (over land). This high-resolution topographical data of the Mississippi coastal region provided the ADCIRC model with improved input with which to calculate improved storm-surge forecasts.

Valenti, Elizabeth↗

An automated system for mosaicking spaceborne SAR imagery

An automated system has been developed for mosaicking spaceborne synthetic aperture radar (SAR) imagery. The system is capable of producing multiframe mosaics for large-scale mapping by combining images in both the along-track direction and adjacent cross-track swaths from ascending and descending passes. The system requires no operator interaction and is capable of achieving high registration accuracy. The output product is a geocoded mosaic on a standard map grid such as UTM or polar stereographic. The procedure described in detail in this paper consists essentially of remapping the individual image frames into these standard grids, frame-to-frame image registration and radiometric smoothing of the seams. These procedures are directly applicable to both the Magellan Venus Mapper and a scanning SAR design such as Radarsat, Eos SAR in addition to merging image frames from traditional SAR systems such as SEASAT and SIR-B. With minor modifications, it may also be applied to spaceborne optical sensor data to generate large-scale mosaics efficiently and with a high degree of accuracy. The system has been tested with SEASAT, SIR-B and Landsat TM data. Examples presented in this paper include a 38-frame mosaic of the Yukon River basin in central Alaska, a 33-frame mosaic of southern California and a three-frame terrain-corrected geocoded mosaic of the Wind River basin in Wyoming.

Kwok, Ronald↗

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗

Cloud Structure and Cloud-Radiation Interaction: From Smallest Observable Scales (Less than 10 m) to GCM Grid Scales (More than 100 km)

In this presentation we review the fractal nature of internal cloud structure from cm- to km-scales as captured by in-situ probes during long horizontal penetrations by aircraft. We uncover the non-Poissonian spatial distribution of cloud droplets at submeter scales and confirm scale-invariant behavior for large scales. Based on these structural characteristics, we generate simple fractal cloud models that reproduce statistical scaling properties of real clouds. These stochastic models represent a link between nonlinear science, in general, and cloud-radiation interaction, in particular. Next we run three-dimensional radiative transfer computations on these synthetic fractal clouds and compare the structure of the resulting radiation fields with the known structure of the cloud model and with satellite images of real clouds. The different behaviors observed for small and large-scale variabilities will be discussed in detail. We find that while the large-scale fluctuations of the resulting radiation fields resemble those in the original scale-invariant cloud structure, the radiation at small scales is much smoother than its cloud liquid water counterpart. This violates scale-invariance and produces a scale-break at 0.2-0.5 km that is clearly observed in high-resolution satellite data such as from Landsat. Finally, we show how radiative transfer Green function theory in the photon diffusion limit explains (and predicts) the above phenomena of "radiative smoothing."

Marshak, Alexander↗

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics↗

Development of a Two Dimensional Synthetic Aperture Radiometer at L-Band

A radiometer that uses aperture synthesis in two dimensions is being built as part of research under NASA's Instrument Incubator Program. The instrument development team consists of engineers at the Goddard Space Flight Center, the University of Massachusetts and Quadrant Engineering. This will be an aircraft instrument operating at L-band which builds on the heritage of Electronically Steered Thinned Array Radiometer (ESTAR). This instrument is a next step in the development of aperture synthesis (STAR technology) to meet the goal of a future mission to monitor soil moisture globally from space. The instrument will be designed to fly on the NASA P-3 aircraft in a nadir pointing mode, although other options are possible. The antenna will consist of an array of modules in a rectangular grid. Each module will be comprised of a printed circuit dual-polarized patch and integrated receiver. The distribution of modules within the rectangular array will be adjustable so that several different imaging configurations (e.g. '+', 'Y', 'T') can be employed. The integrated receiver will provide amplification and conversion to infrared (IF). The IF signal will be routed to a processor where the required correlations performed. The I and Q channels will be created digitally and the correlations will be done digitally in this processor. The digitization will be done with sufficient bits to study the effects of quantization on radiometer performance. A computer/controller will store the data for conversion to an image and will also perform temperature control and other data interfacing and housekeeping tasks. The design of critical components has been completed and hardware is being assembled to test the individual elements. It is expected that a complete two-channel correlator will be tested by the end of 2000 and that the complete instrument will be ready for flight tests the following summer (2001).

LeVine, D. M.↗

Software and Algorithms for Biomedical Image Data Processing and Visualization

A new software equipped with novel image processing algorithms and graphical-user-interface (GUI) tools has been designed for automated analysis and processing of large amounts of biomedical image data. The software, called PlaqTrak, has been specifically used for analysis of plaque on teeth of patients. New algorithms have been developed and implemented to segment teeth of interest from surrounding gum, and a real-time image-based morphing procedure is used to automatically overlay a grid onto each segmented tooth. Pattern recognition methods are used to classify plaque from surrounding gum and enamel, while ignoring glare effects due to the reflection of camera light and ambient light from enamel regions. The PlaqTrak system integrates these components into a single software suite with an easy-to-use GUI (see Figure 1) that allows users to do an end-to-end run of a patient s record, including tooth segmentation of all teeth, grid morphing of each segmented tooth, and plaque classification of each tooth image. The automated and accurate processing of the captured images to segment each tooth [see Figure 2(a)] and then detect plaque on a tooth-by-tooth basis is a critical component of the PlaqTrak system to do clinical trials and analysis with minimal human intervention. These features offer distinct advantages over other competing systems that analyze groups of teeth or synthetic teeth. PlaqTrak divides each segmented tooth into eight regions using an advanced graphics morphing procedure [see results on a chipped tooth in Figure 2(b)], and a pattern recognition classifier is then used to locate plaque [red regions in Figure 2(d)] and enamel regions. The morphing allows analysis within regions of teeth, thereby facilitating detailed statistical analysis such as the amount of plaque present on the biting surfaces on teeth. This software system is applicable to a host of biomedical applications, such as cell analysis and life detection, or robotic applications, such as product inspection or assembly of parts in space and industry.

Talukder, Ashit↗

Development of a Two Dimensional Synthetic Aperture Radiometer at L-Band

A radiometer that uses aperture synthesis in two dimensions is being built as part of research under NASA's Instrument Incubator Program. The instrument development team consists of engineers at the Goddard Space Flight Center, the University of Massachusetts and Quadrant Engineering. This will be an aircraft instrument operating at L-band which builds on the heritage of ESTAR. The choice of L-band was made because the problem of achieving adequate resolution in space is most critical at this wavelength and because a polarimetric, conical scanning airborne radiometer for future experiments to validate soil moisture and ocean salinity retrieval algorithms is not currently available. The instrument will be designed to fly on the NASA P-3 aircraft in a nadir pointing mode, although other options are possible. The antenna will consist of an array of modules arranged in a rectangular grid. Each module will be comprised of a printed circuit dual-polarized patch and integrated receiver. The distribution of modules within the rectangular array will be adjustable so that several different imaging configurations (e.g. "+","Y", "T") can be employed. The integrated receiver will provide amplification and conversion to IF. The IF signal will be routed to a processor where the required correlations performed. The I and Q channels will be created digitally and the correlations will be done digitally in this processor. The digitization will be done with sufficient bits to study the effects of quantization on radiometer performance. A computer/controller will store the data for conversion to an image and will also perform temperature control and other data interfacing and housekeeping tasks. The instrument is currently in the bread boarding phase of development. A design of the critical components has been completed and hardware is being assembled to test the individual elements. It is expected that a complete 2-channel correlator will be tested by the summer of 2000 and that the complete instrument will be ready for flight tests the following summer (2001).

LeVine, D. M.↗

Nature and origin of interstellar diamond from the Allende CV3 meteorite

Data and experimental evidence which support the contention that the C delta diamonds may result from grain-grain collisions in supernova shocks in the interstellar medium are presented. Fragments of the Allende CV3 chondrite were acid-treated. A whitish powder was obtained. For the Analytical Electron Microscopy (AEM) a small drop of ethanol suspension was transferred onto holey carbon support films on 3 mm EM grids. The AEM was performed on transmission-thin fragments of the material which overlay holes in the film, to eliminate interference from the substrate. Electron Spectroscopy for Chemical Analysis (ESCA) was performed on a large aliquot of C. Diamond was identified by selected area electron diffraction. Scanning Transmission Electron Microscope / Energy Dispersive X-ray (STEM-EDS) microanalyses of the C delta diamond, using a light-element detector, show that oxygen and possibly nitrogen are the only impurities consistently present. ESCA spectra from bulk C delta material confirm the presence of N at a level of 0.35 percent or less. Under UV irradiation a yellow-red fluorescence is observed, consistent with that of natural diamonds containing substitutional N. Electron Energy Loss Spectra (EELS) were recorded at 2 eV resolution from the C delta diamond, high pressure synthetic diamond, a diamond film produced in a low pressure plasma by chemical vapor deposition (CVD) on a heated silicon substrate (Roy, 1987), graphite, and amorphous arc sputtered carbon. Comparison of the carbon K edge shape and fine structure shows the Allende C delta phase to be largely diamond, but with a significant pre-edge absorption feature indicative of transitions of C 1s electrons into pi asterisk orbitals which are absent in the purely sp(3)-bonded diamond but present in graphite and amorphous carbon.

Blake, David↗

Air quality impacts from the development of unconventional oil and gas well pads: Air toxics and other volatile organic compounds

Unconventional oil and natural gas development (UOGD) has expanded rapidly across the United States in recent decades and raised concerns about associated air quality impacts. While significant effort has been made to quantify methane emissions, relatively few observations have been made of Volatile Organic Compounds (VOCs), especially during drilling and completion of new wells. Extensive air monitoring during development of several large, multi-well pads in Broomfield, Colorado, in the Denver-Julesburg Basin, provides a novel opportunity to examine changes in local air toxics and other VOC concentrations during well drilling and completions and production. These operations offer an especially useful case to study as several management practices were implemented to reduce emissions (e.g., electrified, grid-powered drill rigs and closed loop fluid handling systems to reduce truck traffic and limit fluid handling on the pad). With simultaneous measurements of methane and 50 VOCs from October 2018 to December 2022 at as many as 19 sites near well pads, in adjacent neighborhoods, and at a more distant reference location, we identify impacts from each phase of well development and production. Use of weekly, time-integrated canisters, a Proton Transfer Reaction Mass Spectrometer (PTR-MS), continuous photoionization detectors (PID) to trigger canister collection upon detection of VOC-rich plumes, and an instrumented vehicle, provided a powerful suite of measurements to characterize both transient plumes and longer-term changes in air quality. Prior to the start of well development, VOC gradients were small across Broomfield. Once drilling commenced, concentrations of oil and gas (O&G) related VOCs, including alkanes and aromatics, increased around active well pads. Concentration increases were clearly apparent during certain operations, including drilling, coil tubing/millout operations, and production tubing installation. Emissions of C 8 –C 10 n-alkanes during drilling operations highlighted the importance of VOC emissions from synthetic drilling mud chosen to reduce odor impacts. More than 90 samples were collected of transient plumes. Using composition measurements, meteorological data, and information about well pad activities, these plumes were connected with specific UOGD operations including drilling, flowback, and production equipment maintenance. The chemical signatures of these plumes differed by operation type (e.g., C 8 –C 10 n-alkanes constituted a larger fraction of measured VOCs in drilling-related plumes). Concentrations of individual, oil and gas-related VOCs in these plumes were often several orders of magnitude higher than in background air, with maximum ethane and benzene concentrations of 79,600 and 819 ppbv, respectively. Because these plumes typically impact a monitoring site for just several minutes, they are easily missed by slower-responding instruments. Study measurements highlight future emission mitigation opportunities during UOGD operations, including better control of emissions from shakers that separate drill cuttings from drilling mud, production separator maintenance operations, and periodic emptying of sand cans during flowback operations.

54 ENVIRONMENTAL SCIENCES↗