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.
The applicability of the various electromagnetic models of scattering from the ocean are reviewed. These models include the small perturbation method, the geometric optics solution, the composite model, and the exact integral equation solution. The restrictions on the electromagnetic models are discussed.
Topics addressed included: determination of sources of propagation and dissipation of atmospheric disturbances and waves; measurement parameters; comparison of different method and instruments; suitable combinations of instruments; information exchange; and training courses.
Simulated radar signatures within the melting layer in stratiform rain, namely the radar bright band, are checked by means of comparisons with simultaneous measurements of the bright band made by the EDOP (X-band) and CRS (W-band) airborne Doppler radars during the CRYSTAL-FACE campaign in 2002. A stratified-sphere model, allowing the fractional water content to vary along the radius of the particle, is used to compute the scattering properties of individual melting snowflakes. Using the effective dielectric constants computed by the conjugate gradient-fast Fourier transform (CGFFT) numerical method for X and W bands, and expressing the fractional water content of melting particle as an exponential function in particle radius, it is found that at X band the simulated radar bright-band profiles are in an excellent agreement with the measured profiles. It is also found that the simulated W-band profiles usually resemble the shapes of the measured bright-band profiles even though persistent offsets between them are present. These offsets, however, can be explained by the attenuation caused by cloud water and water vapor at W band. This is confirmed by the comparisons of the radar profiles made in the rain regions where the un-attenuated W-band reflectivity profiles can be estimated through the X- and W band Doppler velocity measurements. The bright-band model described in this paper has the potential to be used effectively for both radar and radiometer algorithms relevant to the TRMM and GPM satellite missions.
Satellite altimetric data routinely map sea surface topography by measuring the ocean return signal. One source of altimeter measurement contamination occurs when the radar ocean backscatter becomes unusually large, a situation termed a Sigma-0 bloom. Past research suggests Sigma-0 blooms are associated with weak wind and natural surface slick conditions where capillary waves at the air–sea interface are suppressed. To date, no explicit connection between these conditions and Sigma-0 bloom presence has been provided. Using a series of simplified equations, our reduced model determines capillary wave heights from estimates of planktonic carbon concentrations and regional wind speed. Our results suggest that the radar signal reflection increases as capillary wave height decreases. This relationship depends on surfactant concentration, surfactant composition, and wind speed. Model sensitivity analysis indicates that the interface reflectivity depends on biological activity and wind speed. Our proposed simplified model provides a method to identify potential Sigma-0 bloom regions. We conclude that because of the demonstrated impact of biological surfactants on ocean roughness, it is necessary to consider the biological activity, i.e., phytoplankton bloom events, when interpreting signals from radar altimetry and when developing ocean hydrology models.
Plastic or composite pipelines have been the bane of the utility locating industry because conduits are neither conductive nor magnetic, which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipeline where the resistive cover allows for adequate penetration of radar energy. However, GPR cannot be used in areas where the soil cover is conductive due to significant clay and/or salt content. This study takes a comprehensive look at near-surface geophysical methods that potentially are useful for locating buried plastic/composite pipelines, either singly or in combination with other geophysical methods. Specifically, this modeling study uses computational numerical methods to forward model the response of GPR, resistivity, seismic, and gravity gradiometry methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft to 4 ft; (3) various degrees in contrast in physical properties (i.e., permittivity, elasticity, electrical resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar frequencies, seismic source frequency, electrode spacing).
Polynomial functions of time are used to specify the components of the quaternion which represents the nominal attitude of the Venus Radar mapper spacecraft during mapping. The following constraints must be satisfied in order to obtain acceptable synthetic array radar data: the nominal attitude function must have a large dynamic range, the sensor orientation must be known very accurately, the attitude reference function must use as little memory as possible, and the spacecraft must operate autonomously. Fitting polynomials to the components of the desired quaternion function is a straightforward method for providing a very dynamic nominal attitude using a minimum amount of on-board computer resources. Although the attitude from the polynomials may not be exactly the one requested by the radar designers, the polynomial coefficients are known, so they do not contribute to the attitude uncertainty. Frequent coefficient updates are not required, so the spacecraft can operate autonomously.
This study analyzes turbulent energy fluxes in the Arctic atmospheric boundary layer (ABL) using measurements with a small uncrewed aircraft system (sUAS). Turbulent fluxes constitute a major part of the atmospheric energy budget and influence the surface heat balance by distributing energy vertically in the atmosphere. However, only few in situ measurements of the vertical profile of turbulent fluxes in the Arctic ABL exist. The study presents a method to derive turbulent heat fluxes from DataHawk2 sUAS turbulence measurements, based on the flux gradient method with a parameterization of the turbulent exchange coefficient. This parameterization is derived from high-resolution horizontal wind speed measurements in combination with formulations for the turbulent Prandtl number and anisotropy depending on stability. Measurements were taken during the MOSAiC (Multidisciplinary drifting Observatory for the Study of Arctic Climate) expedition in the Arctic sea ice during the melt season of 2020. For three example cases from this campaign, vertical profiles of turbulence parameters and turbulent heat fluxes are presented and compared to balloon-borne, radar, and near-surface measurements. The combination of all measurements draws a consistent picture of ABL conditions and demonstrates the unique potential of the presented method for studying turbulent exchange processes in the vertical ABL profile with sUAS measurements.
The analysis technique and a part of the results obtained from CAT radar echoes from higher troposphere and lower stratosphere are presented. First, the effect of processing distortion caused by the periodogram method using FFT algorithm on the slowly fading ground clutter echo is discussed. It is shown that an extremely narrow clutter spectrum can spill over the entire frequency range if the data are truncated at a tie sorter than their correlation time affecting largely the estimation of the CAT spectrum contribution, especially when the latter is a few tens of dB weaker than the former. A nonlinear least squares fitting procedure is used to parameterize the observed power spectrum in terms of CAT echo power, Doppler shift, spectral width, and the parameters which specify the shape of the clutter component.
In 1971 asteroid Vesta was observed in a region of the sky in which it had never been observed before. Its photometric lightcurve had two distinct maxima. Those observations have been the only strong evidence to support a rotation period of about 10 hr 41 min. Lightcurves made in 1982, when Vesta was at the same aspect as 1971, do not show two different maxima. It is concluded that there was a systematic error in the 1971 observations. At this time a definitive statement cannot be made about the true period of Vesta, although the 5 hr 20 min period does appear more plausible. Radar echoes in 1988 and 1992 should resolve the problem. The shorter rotation period was assumed and the photometric astrometry method applied. The sidereal period is 5 hr 20 min 31.68 sec (0.2225889 + or 0.0000002 days), the rotation is prograde, and the coordinates of the north pole are 103 deg longitude and + 43 deg latitude with an uncertainty of about 6 deg.
A computer scheme for extracting singularities has been developed and applied to single Doppler radar wind data. The scheme is planned for use in real-time wind and singularity analysis and forecasting. The method, known as Doppler Operational Variational Extraction of Singularities is outlined, focusing on the principle of local symmetry. Results are presented from the application of the scheme to a storm-generated gust front in Oklahoma on May 28, 1987.
A technique that uses the radar return from natural targets and at least one trihedral corner reflector to calibrate compressed polarimetric radar data is described. Calibration for relative amplitude and phase, absolute amplitude, and system crosstalk is addressed. The crosstalk calibration method does not require any external calibration targets. Because compressed data are used, the transmitting and receiving systems as must be modeled as reciprocal. Even though the inferred transmit and receive matrices are not simply related to the physical systems, the true Stokes maxtrix for each pixel can be determined accurately. The method is illustrated by estimating the crosstalk parameters of the NASA/JPL aircraft for different types of terrain and for two frequencies. For a C-band system, the crosstalk is less than -20 dB for all ranges in the images. The crosstalk of the L-band system is a function of range and may be as poor as -10 dB in the near range, leading to a noticeable distortion of the polarization signatures.
A hybrid Finite Element Method (FEM)/Method of Moments (MoM) technique in conjunction with the Asymptotic Waveform Evaluation (AWE) technique is applied to obtain radar cross section (RCS) of a cavity-backed aperture in an infinite ground plane over a frequency range. The hybrid FEM/MoM technique when applied to the cavity-backed aperture results in an integro-differential equation with electric field as the unknown variable, the electric field obtained from the solution of the integro-differential equation is expanded in Taylor series. The coefficients of the Taylor series are obtained using the frequency derivatives of the integro-differential equation formed by the hybrid FEM/MoM technique. The series is then matched via the Pade approximation to a rational polynomial, which can be used to extrapolate the electric field over a frequency range. The RCS of the cavity-backed aperture is calculated using the electric field at different frequencies. Numerical results for a rectangular cavity, a circular cavity, and a material filled cavity are presented over a frequency range. Good agreement between AWE and the exact solution over the frequency range is obtained.
Changes in the shape of the Earth's surface can be routinely measured with precisions better than centimeters. Processes below the surface often drive these changes and as a result, investigators require models with inversion methods to characterize the sources. Simplex inverts any combination of GPS (global positioning system), UAVSAR (uninhabited aerial vehicle synthetic aperture radar), and InSAR (interferometric synthetic aperture radar) data simultaneously for elastic response from fault and fluid motions. It can be used to solve for multiple faults and parameters, all of which can be specified or allowed to vary. The software can be used to study long-term tectonic motions and the faults responsible for those motions, or can be used to invert for co-seismic slip from earthquakes. Solutions involving estimation of fault motion and changes in fluid reservoirs such as magma or water are possible. Any arbitrary number of faults or parameters can be considered. Simplex specifically solves for any of location, geometry, fault slip, and expansion/contraction of a single or multiple faults. It inverts GPS and InSAR data for elastic dislocations in a half-space. Slip parameters include strike slip, dip slip, and tensile dislocations. It includes a map interface for both setting up the models and viewing the results. Results, including faults, and observed, computed, and residual displacements, are output in text format, a map interface, and can be exported to KML. The software interfaces with the QuakeTables database allowing a user to select existing fault parameters or data. Simplex can be accessed through the QuakeSim portal graphical user interface or run from a UNIX command line.
Estimation of path attenuation is a critical part of retrieving precipitation parameters using measurements from the Dual-Frequency Precipitation Radar (DPR) on board the Global Precipitation Measurement Mission (GPM) satellite. In this paper, we describe the latest implementation of the Surface Reference Technique that uses surface scattering properties to infer path attenuation through the precipitation. Both single- and dual-frequency versions of this method are available and while the dual-frequency version appears to be more accurate at moderate rain rates, the single-frequency approach at Ku-band is needed when the Ka-band data are not available. Despite improvements afforded by the dual-frequency version of the method, other methods such as the Hitschfeld-Bordan and standard dual-frequency approaches offer advantages particularly at lighter rain rates and at near-nadir incidence angles over land. Weighted averages of the results from these methods appear to offer the best estimate of path attenuation presently available.
Compact polarimetry has been shown to be an interesting alternative mode to full polarimetry when global coverage and revisit time are key issues. It consists on transmitting a single polarization, while receiving on two. Several critical points have been identified, one being the Faraday rotation (FR) correction and the other the calibration. When a low frequency electromagnetic wave travels through the ionosphere, it undergoes a rotation of the polarization plane about the radar line of sight for a linearly polarized wave, and a simple phase shift for a circularly polarized wave. In a low frequency radar, the only possible choice of the transmit polarization is the circular one, in order to guaranty that the scattering element on the ground is illuminated with a constant polarization independently of the ionosphere state. This will allow meaningful time series analysis, interferometry as long as the Faraday rotation effect is corrected for the return path. In full-polarimetric (FP) mode, two techniques allow to estimate the FR: Freeman method using linearly polarized data, and Bickel and Bates theory based on the transformation of the measured scattering matrix to a circular basis. In CP mode, an alternate procedure is presented which relies on the bare surface scattering properties. These bare surfaces are selected by the conformity coefficient, invariant with FR. This coefficient is compared to other published classifications to show its potential in distinguishing three different scattering types: surface, doublebounce and volume. The performances of the bare surfaces selection and FR estimation are evaluated on PALSAR and airborne data. Once the bare surfaces are selected and Faraday angle estimated over them, the correction can be applied over the whole scene. The algorithm is compared with both FP techniques. In the last part of the paper, the calibration of a CP system from the point of view of classical matrix transformation methods in polarimetry is proposed.
Over the past two decades, detailed computational simulations of the intricate three dimensional structures of ice-phase crystals and aggregates of those crystals have been developed. The microwave single-scattering properties of these simulated ice particles have been computed and used to improve quantitative estimates of snow rates and to better deRne the vertical structure of snow water contents in deep convective systems, as derived from satellite-borne passive microwave and/or radar remote sensing measurements. The same icephase particles have more recently been used as the starting point for simulations of melting precipitation using computational melting methods. In the current study, a heuristic melting method, as well as a physically-based melting procedure based on smoothed-particle hydrodynamics, are applied to ice particle models to describe the full evolution of the particles from dry snow to liquid drops. The discrete dipole approximation is utilized to calculate the single-scattering properties of the mixed-phase particles throughout the melting process. Then, the properties of the particles are “mapped” into simpliRed microphysical simulations of particle size spectra in the melting layers of stratiform, precipitating clouds. The bulk single-scattering properties of the melting layers and the sensitivity of their properties to modeling assumptions are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.
Range resolution of a conventional pulsed Doppler radar is determined by the scattering volume defined by the transmitted pulse shape. To increase the resolution, the length of the pulse must be reduced. Reducing the pulse length also reduces the transmitted power and hense the signal to noise ratio unless the peak power capability of the transmitter is greatly increased. Improved range resolution may also be attained through the use of various pulse coding methods, but such methods are sometimes difficult to implement from a hardware standpoint. The frequency-hopping (F-H) technique described increases the range resolution of pulse Doppler MST (mesosphere stratosphere troposphere) radar without the need for extensive modifications to the radar transmitter. This technique consists of sending a repeated sequence of pulses, each pulse in the sequence being transmitted at a unique radio frequency that is under the control of a microcomputer. This technique is discussed along with other radar parameters.
Soil moisture is a key parameter in understanding the global water cycle and in predicting natural hazards. Polarimetric radar measurements have been used for estimating soil moisture of bare surfaces. In order to estimate soil moisture accurately, the surface roughness effect must be compensated properly. In addition, these algorithms will not produce accurate results for vegetated surfaces. It is difficult to retrieve soil moisture of a vegetated surface since the radar backscattering cross section is sensitive to the vegetation structure and environmental conditions such as the ground slope. Therefore, it is necessary to develop a method to estimate the effect of the surface roughness and vegetation reliably. One way to remove the roughness effect and the vegetation contamination is to take advantage of the temporal variation of soil moisture. In order to understand the global hydrologic cycle, it is desirable to measure soil moisture with one- to two-days revisit. Using these frequent measurements, a time series approach can be implemented to improve the soil moisture retrieval accuracy.