SEARCH · Engineering Papers
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.
Utilization of remote sensor data by Virginia state agencies including aerial photography, radar detection, and infrared detectors
A system for deriving direct digital indications of frequency and phase difference between two incoming pulse trains adaptable for collision avoidance systems or the like. A pair of radar beams 152 and 152A are directed toward a target 153 and corresponding beams 154 and 154A returning therefrom are detected. A digital difference circuit 110 forms a pulse train 66 from the Doppler shift frequencies of each beam pair having a repetition rate functionally related to the difference in magnitude of the shift frequencies. Pulses from the pulse train are counted as a function of time. Visual indications thereof on display 144 are correlative to target position relative to beams 152 and 152A.
Abstract A new methodology is developed to construct a global (60°S–60°N) long‐term (2000–2019) high‐resolution (∼10‐km h) mesoscale convective system (MCS) database by tracking MCS jointly using geostationary satellite infrared brightness temperature ( T b ) and precipitation feature (PF) characteristics from the Integrated Multi‐satellitE Retrievals for GPM precipitation data sets. Independent validation shows that the satellite‐based MCS data set is able to reproduce important MCS statistics derived from ground‐based radar network observations in the United States and China. We show that by carefully considering key PF characteristics in addition to T b signatures, the new method significantly improves upon previous T b ‐only methods in detecting MCSs in the midlatitudes for all seasons. Results show that MCSs account for over 50% of annual total rainfall across most of the tropical belt and in selected regions of the midlatitudes, with a strong seasonality over many regions of the globe. The tracking database allows Lagrangian aspects such as MCS lifetime and translational speed and direction to be analyzed. The longest‐lived MCSs preferentially occur over the subtropical oceans. The land MCSs have higher cloud‐tops associated with more intense convection, and oceanic MCSs have much higher rainfall production. While MCSs are observed in many regions of the globe, there are fundamental differences in their dynamic and thermodynamic structures that warrant a better understanding of processes that control their evolution. This global database provides significant opportunities for observational and modeling studies of MCSs, their characteristics, and roles in regional and global water and energy cycles, as well as their hydrologic and other impacts.
Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.
Method for estimation of atmospheric turbulence encountered by airplanes utilizes wealth of data captured by multichannel digital flight-data recorders and air-traffic-control radar. Developed as part of continuing effort to understand how airplanes respond to such potentially hazardous phenomena as: clear-air turbulence generated by destabilized wind-shear layers above mountains and thunderstorms, and microbursts (intense downdrafts striking ground), associated with thunderstorms. Reconstructed wind fields used to predict and avoid future hazards.
The planetary ephemerides approximation for radar astronomy is discussed, and, in particular, the effect of this approximation on the performance of the programmable local oscillator (PLO) used in Goldstone Solar System Radar is presented. Four different approaches are considered and it is shown that the Gram polynomials outperform the commonly used technique based on Chebyshev polynomials. These methods are used to analyze the mean square, the phase error, and the frequency tracking error in the presence of the worst case Doppler shift that one may encounter within the solar system. It is shown that in the worst case the phase error is under one degree and the frequency tracking error less than one hertz when the frequency to the PLO is updated every millisecond.
This dissertation focuses on the signal processing problems associated with the detection of hazardous windshears using airborne Doppler radar when weak weather returns are in the presence of strong clutter returns. In light of the frequent inadequacy of spectral-processing oriented clutter suppression methods, we model a clutter signal as multiple sinusoids plus Gaussian noise, and propose adaptive filtering approaches that better capture the temporal characteristics of the signal process. This idea leads to two research topics in signal processing: (1) signal modeling and parameter estimation, and (2) adaptive filtering in this particular signal environment. A high-resolution, low SNR threshold maximum likelihood (ML) frequency estimation and signal modeling algorithm is devised and proves capable of delineating both the spectral and temporal nature of the clutter return. Furthermore, the Least Mean Square (LMS) -based adaptive filter's performance for the proposed signal model is investigated, and promising simulation results have testified to its potential for clutter rejection leading to more accurate estimation of windspeed thus obtaining a better assessment of the windshear hazard.
Significant efforts have been made to decompose polarimetric radar data into several simple scattering components. The components which are selected because of their physical significance can be used to classify SAR (Synthetic Aperture Radar) image data. If particular components can be related to forest parameters, inversion procedures may be developed to estimate these parameters from the scattering components. Several methods have been used to decompose an averaged Stoke's matrix or covariance matrix into three components representing odd (surface), even (double-bounce) and diffuse (volume) scatterings. With these decomposition techniques, phenomena, such as canopy-ground interactions, randomness of orientation, and size of scatters can be examined from SAR data. In this study we applied the method recently reported by van Zyl (1992) to decompose averaged backscattering covariance matrices extracted from JPL SAR images over forest stands in Maine, USA. These stands are mostly mixed stands of coniferous and deciduous trees. Biomass data have been derived from field measurements of DBH and tree density using allometric equations. The interpretation of the decompositions and relationships with measured stand biomass are presented in this paper.
This paper represents the first attempt to use TRMM rainfall information to estimate the four dimensional latent heating structure over the global tropics for February 1998. The mean latent heating profiles over six oceanic regions (TOGA COARE IFA, Central Pacific, S. Pacific Convergence Zone, East Pacific, Indian Ocean and Atlantic Ocean) and three continental regions (S. America, Central Africa and Australia) are estimated and studied. The heating profiles obtained from the results of diagnostic budget studies over a broad range of geographic locations are used to provide comparisons and indirect validation for the heating algorithm estimated heating profiles. Three different latent heating algorithms, the Goddard Convective-Stratiform (CSH) heating, the Goddard Profiling (GPROF) heating, and the Hydrometeor heating (HH) are used and their results are intercompared. The horizontal distribution or patterns of latent heat release from the three different heating retrieval methods are quite similar. They all can identify the areas of major convective activity (i.e., a well defined ITCZ in the Pacific, a distinct SPCZ) in the global tropics. The magnitude of their estimated latent heating release is also not in bad agreement with each other and with those determined from diagnostic budget studies. However, the major difference among these three heating retrieval algorithms is the altitude of the maximum heating level. The CSH algorithm estimated heating profiles only show one maximum heating level, and the level varies between convective activity from various geographic locations. These features are in good agreement with diagnostic budget studies. By contrast, two maximum heating levels were found using the GPROF heating and HH algorithms. The latent heating profiles estimated from all three methods can not show cooling between active convective events. We also examined the impact of different TMI (Multi-channel Passive Microwave Sensor) and PR (Precipitation Radar) rainfall information on latent heating structures.
A characteristic feature of rainfall statistics is that they depend on the space and time scales over which rain data are averaged. A previously developed spectral model of rain statistics that is designed to capture this property, predicts power law scaling behavior for the second moment statistics of area-averaged rain rate on the averaging length scale L as L right arrow 0. In the present work a more efficient method of estimating the model parameters is presented, and used to fit the model to the statistics of area-averaged rain rate derived from gridded radar precipitation data from TOGA COARE. Statistical properties of the data and the model predictions are compared over a wide range of averaging scales. An extension of the spectral model scaling relations to describe the dependence of the average fraction of grid boxes within an area containing nonzero rain (the "rainy area fraction") on the grid scale L is also explored.
Stormwater runoff plumes, municipal wastewater plumes, and natural hydrocarbon seeps are important pollution hazards for the heavily populated Southern California Bight (SCB). Due to their small size, dynamic and episodic nature, these hazards are difficult to sample adequately using traditional in situ oceanographic methods. Complex coastal circulation and persistent cloud cover can further complicate detection and monitoring of these hazards. We use imagery from space-borne synthetic aperture radar (SAR), complemented by field measurements, to examine these hazards in the SCB. The hazards are detectable in SAR imagery because they deposit surfactants on the sea surface, smoothing capillary and small gravity waves to produce areas of reduced backscatter compared with the surrounding ocean. We suggest that high-resolution SAR, which obtains useful data regardless of darkness or cloud cover, could be an important observational tool for assessment and monitoring of coastal marine pollution hazards in the SCB and other urbanized coastal regions.
Shallow subsurface ice deposits on Mars and the Moon have been identified as potential targets for in situ resource utilization by future crewed missions to the surfaces of these bodies. On Mars, such deposits have been identified through a variety of orbital measurements at the poles, mid-latitudes, and even preserved by volcanic ash near the Tharsis Montes. At some locations, massive ice deposits have been documented to be 100 m in thickness at depths as shallow as 1 m. On the Moon, the presence of thick, massive ice has not been confirmed; however, orbital observations indicate that the upper 1-2 m of regolith may contain a few weight percent in thermally favorable locations. Eruptions of the Askja Volcano in March, 1875 and November, 1961 deposited abundant pyroclasts, which blanketed and insulated fresh snowfall that later densified into massive ice and is preserved today. The pyroclasts consist of a buff colored pumice from the eruption in 1875 and basaltic lapilli and ash from the eruption in 1961. The largely unvegetated, unconsolidated nature of the pyroclasts and their stratigraphic relationships with shallow subsurface ice make them potentially analogous to some ice deposits within regolith and pyroclasts at the Moon and Mars. Ground penetrating radar(GPR) can be used by future human or robotic missions to identify water ice similarly preserved in the shallow subsurface of these bodies. Our team conducted over 66 GPR surveys inside and surrounding the Askja caldera in August 2019 at 200 and 400 MHz to map both the volcanic deposits and the subsurface ice deposits. We also used a hammer drill augur to take boreholes in order to confirm subsurface stratigraphy and the presence of ice down to 1-1.5 m and aerial surveys using a Mavic 2 Pro quadcopter at each site for additional context. We observed shallow ice deposits within pore spaces of the pyroclasts at depths of ~15-30 cm and pure ice deposits at varying depths (0.6-1 m) as thick as 2-3 m. We also plan to use the data to characterize the frequency-dependent attenuation of the radar signal as it travels through volcanoclastic material and ice in the subsurface. Our investigation of Askja as a planetary analog will provide insight into analytical methods that can be used to investigate subsurface water ice from surface operations at other terrestrial bodies.
The capability to assess the mission performance of novel navigation technologies from a systems-level approach and to provide quantitative results is crucial. With growing interest and innovations from NASA’s commercial partners, data-driven results that quantify the technological impact will guide research developments and facilitate stakeholder decision-making while requiring less time and resources. This paper presents a method to evaluate various sensor combinations and their impact on the overall system performance using an existing six degrees-of-freedom, physics-based engineering simulation for a government reference lunar lander as a testbed. Selected navigation technologies over various technology readiness levels, including Inertial Measurement Units, Navigation Doppler Lidar, and radar altimeter-radar velocimeter, were studied in this paper. Results using this testbed to perform sensitivity studies and to provide quantitative assessments are reported. One key finding is that there are diminishing returns for reducing sensor errors. The most influential sensor parameter to landing success are vehicle configuration and mission dependent. Results from this method could be used by stakeholders to make data-driven systems-based decisions and by technology developers as guidance for the specific parameter improvements that will have the most significant impact on mission success. This capability could be further used to assess alternative scenarios should one type of technology become unavailable or to re-assess initial assumptions and identify appropriate requirements from a system perspective.
Within the intermountain west, monitoring fuel loads is a major concern for wildland fire management efforts. To address this concern, we partnered with the U.S. Forest Service to inform the agency which forested areas should be prioritized for prescribed burning and fuel reduction near human communities in the Bridger-Teton National Forest, Wyoming. We created burn maps, fuel load maps, and a tutorial document to identify forest impact trends and provide the partner with the tools to replicate project methods for expansion to other wildfire crisis strategy sites. These end products were made using two NASA Earth observations: Landsat 8 Operational Land Imager and Shuttle Radar Topography Mission. Based on our random forest analysis, our maps identified 998 acres within the Wildland Urban Interface that are predicted to have high fuel loading and high burn severity within the Bridger-Teton National Forest. Forested areas closer to heavily populated areas such as Jackson, Kelly, Moran, New Forks Lake, and Star Valley Ranch should be prioritized for fuel reduction. However, our random forest model analysis was limited to using vegetation and topographical indices with no field data for model validation. Therefore, future studies should use field data for model validation to improve model accuracy and additionally incorporate Global Ecosystem Dynamics Investigation data into models to create better predictions of forested areas with high fuel load and high burn severity.
High-resolution near-surface imaging has important applications in civil engineering, infrastructure inspection, military threat detection, geological characterization, and lunar and planetary exploration. Zero-offset, singlechannel ground penetrating radar (GPR) imaging is an established technique for near-surface target imaging and sensing but often suffers from low spatial resolution and imaging artifacts, especially of deep structures. In response, we formulate the GPR imaging as a dual-sparsity optimization problem, and develop a super-resolution electromagnetic imaging method based on a fast iterative shrinkage-thresholding algorithm. We develop our GPR imaging method in the framework of electromagnetic exploding-reflectors simulation theory, therefore the imaging method is computationally efficient. In this work, we demonstrate through synthetic and field data examples that our method can produce sharper, more reliable images with fewer artifacts compared with single-pass reverse-time migration GPR method, thus leading to improved near-surface interpretation and object identification.
The Federal Aviation Administration's (FAA) Terminal Doppler Weather Radar (TDWR) system was primarily designed to address the operational needs of pilots in the avoidance of low-altitude wind shears upon takeoff and landing at airports. One of the primary methods of wind-shear detection for the TDWR system is the gust-front detection algorithm. The algorithm is designed to detect gust fronts that produce a wind-shear hazard and/or sustained wind shifts. It serves the hazard warning function by providing an estimate of the wind-speed gain for aircraft penetrating the gust front. The gust-front detection and wind-shift algorithms together serve a planning function by providing forecasted gust-front locations and estimates of the horizontal wind vector behind the front, respectively. This information is used by air traffic managers to determine arrival and departure runway configurations and aircraft movements to minimize the impact of wind shifts on airport capacity. This paper describes the gust-front detection and wind-shift algorithms to be fielded in the initial TDWR systems. Results of a quantitative performance evaluation using Doppler radar data collected during TDWR operational demonstrations at the Denver, Kansas City, and Orlando airports are presented. The algorithms were found to be operationally useful by the FAA airport controllers and supervisors.
Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.
Double-parameter expansions, leading to an instantaneous multiscale analysis, are used to obtain time-frequency graphs from time-series data for interpretation of scatterometer response to sea surface geometry. This method is used to define, as well as possible, one-dimensional signal input (elevation and slope). However, some important fundamental parameters, such as the asymmetric development of short scales in crests and troughs, the long-wave amplitude modulation, the different zero-crossings, and so forth, are easily extracted from the wavelet images.