Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizontal interpolation”

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 37 records · Page 2

Forecast experiments with the NASA/GLA stratospheric/tropospheric data assimilation system

For the first time, a four-dimensional stratospheric/tropospheric data assimilation system with a top analysis level at 0.4 mb has been developed and used to produce physically consistent gridded analyses for the stratosphere as well as the troposphere for a period during the First GARP Global Experiment (FGGE) and Limb Infrared Monitor of the Stratosphere (LIMS) (November 1978-May 1979). The system consists of a two-dimensional optimum interpolation analysis with 18 mandatory pressure levels and a 19-level fourth order stratospheric/tropospheric general circulation model with a horizontal resolution of 4 (latitude) by 5 deg (longitude) and a top at 0.3 mb. The system allows the utilization of stratospheric data including LIMS, Tiros-N retrievals, rocketsondes and vertical temperature profile radiometer soundings in addition to the other FGGE level 2b data. These data are analyzed every six hours. In order to examine the quality of the analyzed data, forecast experiments starting from different analyses are performed for the period of the stratospheric sudden warming of late February 1979. The results indicate that by employing the present four-dimensional assimilation approach, the medium-range forecast skill for this event is improved.

Takano, Kenji↗

Receive Mode Analysis and Design of Microstrip Reflectarrays

Traditionally microstrip or printed reflectarrays are designed using the transmit mode technique. In this method, the size of each printed element is chosen so as to provide the required value of the reflection phase such that a collimated beam results along a given direction. The reflection phase of each printed element is approximated using an infinite array model. The infinite array model is an excellent engineering approximation for a large microstrip array since the size or orientation of elements exhibits a slow spatial variation. In this model, the reflection phase from a given printed element is approximated by that of an infinite array of elements of the same size and orientation when illuminated by a local plane wave. Thus the reflection phase is a function of the size (or orientation) of the element, the elevation and azimuth angles of incidence of a local plane wave, and polarization. Typically, one computes the reflection phase of the infinite array as a function of several parameters such as size/orientation, elevation and azimuth angles of incidence, and in some cases for vertical and horizontal polarization. The design requires the selection of the size/orientation of the printed element to realize the required phase by interpolating or curve fitting all the computed data. This is a substantially complicated problem, especially in applications requiring a computationally intensive commercial code to determine the reflection phase. In dual polarization applications requiring rectangular patches, one needs to determine the reflection phase as a function of five parameters (dimensions of the rectangular patch, elevation and azimuth angles of incidence, and polarization). This is an extremely complex problem. The new method employs the reciprocity principle and reaction concept, two well-known concepts in electromagnetics to derive the receive mode analysis and design techniques. In the "receive mode design" technique, the reflection phase is computed for a plane wave incident on the reflectarray from the direction of the beam peak. In antenna applications with a single collimated beam, this method is extremely simple since all printed elements see the same angles of incidence. Thus the number of parameters is reduced by two when compared to the transmit mode design. The reflection phase computation as a function of five parameters in the rectangular patch array discussed previously is reduced to a computational problem with three parameters in the receive mode. Furthermore, if the beam peak is in the broadside direction, the receive mode design is polarization independent and the reflection phase computation is a function of two parameters only. For a square patch array, it is a function of the size, one parameter only, thus making it extremely simple.

Rengarajan, Sembiam↗

Chronology, eruption duration, and atmospheric contribution of Apollinaris Patera, Mars

Geologic mapping from Viking image data of the Martian volcano Apollinaris Patera allowed identification of 6 major events that shaped its current morphology. Derivation of new topographic data allowed accurate estimates of the volume of erupted products from which estimates of an eruption duration are presented for the edifice and its corresponding atmospheric water contribution. Topographic data were acquired using stereophotogrammetric (Viking picnos 603A42, 639A92) techniques in both profiling and contouring modes (1 km contour). The profiling mode results in a more precise measurement than the contouring mode, but is limited in its areal coverage. The contour data are used in a more general sense to provide synoptic coverage for the volcano. These new stereophotogrammetric measurements constrain the topography to an accuracy of approximately 800 m vertically and approximately 1000 m horizontally. Conversion of the derived contour map to a raster based digital elevation model (DEM) was done by a growing contours interpolation. An ambiguity arises in this volume calculation due to an uncertainty regarding the actual base of the volcano relative to the pre-existing topography. Thus, the volume estimate was bracketed by using the 100 m and the 750 m elevations, which generally correspond to the lower and upper portions of the circumferential basal scarp, respectively. The resulting volumes are 103,000 cu km and 97,000 cu km; therefore the total volume of Apollinaris Patera is estimated to be approximately 100,000 cu km.

Robinson, M. S.↗

Determination of stratospheric temperature and height gradients from Nimbus 3 radiation data

To improve the specification of stratospheric, horizontal, temperature and geopotential height fields needed for high-flying aircraft, we derived a technique to estimate data between satellite tracks using interpolated infrared interferometer spectrometer 15-micron radiation data from Nimbus 3. The interpolation is based on the observed gradients of the medium resolution infrared 15-micron radiances between subsatellite tracks. The technique was verified with radiosonde data taken within 6 hr of the satellite data. Comparison indicates that the technique developed here produces analyses that are in general agreement with those from radiosonde data. In addition, this technique provides details over areas of sparse data not shown by conventional techniques.

Nicholas, G. W.↗

Design considerations for a real-time ocular counterroll instrument

A real-time algorithm for measuring three-dimensional movement of the human eye, especially torsional movement, is presented. As its input, the system uses images of the eyeball taken at video rate. The amount of horizontal and vertical movement is extracted using a pupil tracking technique. The torsional movement is then measured by computing the discrete cross-correlation function between the circular samples of successive images of the iris patterns and searching for the position of the peak of the function. A local least square interpolation around the peak of the cross-correlation function is used to produce nearly unbiased estimates of torsion angle with accuracy of about 3-4 arcmin. Accuracies of better than 0.03 deg are achievable in torsional measurement with SNR higher than 36 dB. Horizontal and vertical rotations of up to + or - 13 deg can occur simultaneously with torsion without introducing any appreciable error in the counterrolling measurement process.

Hatamian, M.↗

Objective analysis of tidal fields in the Atlantic and Indian Oceans

An objective analysis technique has been developed to extrapolate tidal amplitudes and phases over entire ocean basins using existing gauge data and the altimetric measurements which are now beginning to be provided by satellite oceanography. The technique was previously tested in the Lake Superior basin. The method has now been developed and applied in the Atlantic-Indian ocean basins using a 6 deg x 6 deg grid to test its essential features. The functions used in the interpolation are the eigenfunctions of the velocity potential (Proudman functions) which are computed numerically from a knowledge of the basin's bottom topography, the horizontal plan form and the necessary boundary conditions. These functions are characteristic of the particular basin. The gravitational normal modes of the basin are computed as part of the investigation, they are used to obtain the theoretical forced solutions for the tidal constituents, the latter provide the simulated data for the testing of the method and serve as a guide in choosing the most energetic modes for the objective analysis. The results of the objective analysis of the M2 and K1 tidal constituents indicate the possibility of recovering the tidal signal with a degree of accuracy well within the error bounds of present day satellite techniques.

Sanchez, B. V.↗

Intercomparison of haloe and ER-2 aircraft H2O and CH4 observations collected during the second Airborne Arctic Stratospheric Experiment (AASE-II)

HALOE observations of H2O and CH4 are compared with in situ techniques aboard the ER-2 aircraft during the northern winter of 1991/92, in particular for the dates 911208, 920108, 920217, 920222, and 920320 when the spatial coincidences are close, within +/- 1 deg latitude and +/- 12 deg longitude. The results reveal the limitations of comparing high resolution in situ aircraft data with a remote sounding limb scanner. Of the five comparison dates, three had HALOE/ER-2 coincidences which occurred near the edge of the Arctic vortex; the vertical variability in the HALOE results and the horizontal variability in the ER-2 observations on these days show that the vortex edge is not a region where exact agreement can be expected except by chance. On a 4th comparison date, there was substantial overlap away from the vortex edge, although for some species the aircraft data show considerable variability near the coincidence point. On the 5th comparison date, the ER-2 had no overlap in altitude with the lowest HALOE observations; however, a short linear interpolation over about 1 km altitude results in smooth composite profiles. Generally speaking the agreement between HALOE and the ER-2 at overlap altitudes is about 12 percent in the case of water vapor, which shows low horizontal and vertical variability. The agreement for methane, with limited data having altitude overlap, is better than 6 percent.

Tuck, A. F.↗

Determination of stratospheric temperature and height gradients from nimbus 3 radiation data

To improve the specification of stratospheric horizontal temperature and geopotential height fields from satellite radiation data, needed for high flying aircraft, a technique was derived to estimate data between satellite tracks using interpolated IRIS 15-micron data from Nimbus III. The interpolation is based on the observed gradients of the MRIR 15-micron radiances between subsatellite tracks. The technique was verified with radiosonde data taken within 6 hours of the satellite data. The sample varied from 1126 pairs at low levels to 383 pairs at 10 mb using northern hemisphere data for June 15 to July 20, 1969. The data were separated into five latitude bands. The Rms temperature differences were generally from 2 to 5 C for all levels above 300 mb. From 500 to 300 mb RMS differences vary from 4 to 9C except at high latitudes which show values near 3C. The RMS differences between radiosonde heights and those calculated hydrostatically from the surface were from 30 to 280 meters increasing from the surface to 10 mb. Integration starting at 100 mb reduced the RMS difference in the stratosphere to 20 to 120 meters from 70 to 10 mb. From a comparison with actual operational maps at 50 and 10 mb, it appears the techniques developed produce analyses in general agreement with those from radiosonde data. In addition, they are able to indicate details over areas of sparse data not shown by conventional techniques.

Nicholas, G. W.↗

A New Lunar Digital Elevation Model from the Lunar Orbiter Laser Altimeter and SELENE Terrain Camera

We present an improved lunar digital elevation model (DEM) covering latitudes within +/-60 deg, at a horizontal resolution of 512 pixels per degree ( approx.60 m at the equator) and a typical vertical accuracy approx.3 to 4 m. This DEM is constructed from approx.4.5 ×10(exp 9) geodetically-accurate topographic heights from the Lunar Orbiter Laser Altimeter (LOLA) onboard the Lunar Reconnaissance Orbiter, to which we co-registered 43,200 stereo-derived DEMs (each 1 deg×1 deg) from the SELENE Terrain Camera (TC) ( approx.10(exp 10) pixels total). After co-registration, approximately 90% of the TC DEMs show root-mean-square vertical residuals with the LOLA data of < 5 m compared to approx.50% prior to co-registration. We use the co-registered TC data to estimate and correct orbital and pointing geolocation errors from the LOLA altimetric profiles (typically amounting to < 10 m horizontally and < 1 m vertically). By combining both co-registered datasets, we obtain a near-global DEM with high geodetic accuracy, and without the need for surface interpolation. We evaluate the resulting LOLA + TC merged DEM (designated as "SLDEM2015") with particular attention to quantifying seams and crossover errors.

Barker, M. K.↗

A comparison between Nimbus 5 THIR and ITPR temperatures and derived winds with rawinsonde data obtained in the AVE 2 experiment

During the period of May 11 and 12, 1974, NASA conducted its second Atmospheric Variability Experiment (AVE II) over the eastern United States. In this time interval, two Nimbus 5 orbits crossed the AVE II area, providing a series of ITPR soundings as well as THIR data. Horizontal temperature mapping of the AVE II cloud field is examined using two grid print map scales. Implied cloud top heights are compared with maximum radar-echo top reports. In addition, shelter temperatures in areas of clear sky are compared with the surface temperatures as determined from 11.5 micrometer radiometer data of the THIR experiment. The ITPR sounding accuracy is evaluated using interpolated radiosonde temperatures at times nearly coincident with the ITPR soundings. It was found that mean differences between the two data sets were as small as 1.3 C near 500 mb and as large as 2.9 C near the tropopause. The differences between ITPR and radiosonde temperatures at constant pressure levels were sufficient to induce significant differences in the horizontal temperature gradient. Cross sections of geostrophic wind along the orbital tracks were developed using a thermal wind buildup based on the ITPR temperature data and the radiosonde temperature data. Differences between the radiosonde and ITPR geostrophic winds could be explained on the basis of differences in the ITPR and radiosonde temperature gradients.

Arnold, J. E.↗

Interpolation problems in meteor radar analysis, part 7.6A

Meteor echoes come from random points in the observation volume, and are irregularly spaced in time. This precludes the use of fast fourier transformations (FFT) techniques on the raw data to give the spectrum of waves that are present. One way around this obstacle is to restrict our interest to a particular class of waves, and fit a corresponding model to the raw data. It is assumed that there is no horizontal variation across the observation volume for tides, but in the vertical this is certainly not the case. If, in addition, we are interested in other types of waves which may be present and whose periods are unknown, then examining the raw line-of-sight velocities does not tell us how to modify the model, since the line-of-sight direction is not fixed. This is the motivation for interpolation. Interpolation takes a temporal series of line-of-sight velocities, and transforms it to a temporal series of wind velocities for each orthogonal direction. The velocities along a given direction can then be examined readily for any waves in addition to tides.

Tetenbaum, D.↗

Accuracy of diagnostic heat and moisture budgets using SESAME-79 field data as revealed by observing system simulation experiments

Observing system simulation experiments are used to investigate the accuracy of diagnostic heat and moisture budgets which employ the AVE-SESAME 1979 data. The time-including, four-dimensional data set of a mesoscale model is used to simulate rawinsonde observations from AVE-SESAME 1979. The 5 C/day (heat budget) and 2 g/kg per day (moisture budget) magnitudes of error obtained indicate difficulties in the diagnosis of the heating rate in weak convective systems. The influences exerted by observational frequency, objective analysis, observational density, vertical interpolation, and observational errors on the budget are also studied, and it is found that the temporal and spatial resolution of the SESAME regional network is marginal for diagnosing convective effects on a horizontal time scale of 550 x 550 km.

Kuo, Y.-H.↗

Assimilation of global cloud imagery from multiple satellites

Synoptic images of the global cloud field are created from measurements of IR radiation taken aboard four geostationary and two polar-orbiting platforms simultaneously observing the earth's cloud field. A series of spatial and temporal interpolations, together with data reliability criteria, are used to map synoptic measurements from the geostationary satellites and asynoptic measurements from the polar-orbiting satellites. The synoptic global cloud imagery (GCI) that results has a horizontal resolution of 0.7 deg in longitude by 0.35 deg in latitude and a temporal resolution of 3 h, providing an unprecedented view of the global cloud pattern. Each composite image represents a nearly instantaneous snapshot of the entire earth's cloud field. Collectively, the composite imagery resolve, on a global basis, most of the variability inherent to organized convection, including several harmonics of the diurnal cycle. Because of its homogeneous properties, the GCI lends itself easily to a variety of space-time analyses useful for studying global cloud behavior.

Tanaka, Ken↗

Frequency and Angular Variations of Land Surface Microwave Emissivities: Can we Estimate SSM/T and AMSU Emissivities from SSM/I Emissivities?

To retrieve temperature and humidity profiles from SSM/T and AMSU, it is important to quantify the contribution of the Earth surface emission. So far, no global estimates of the land surface emissivities are available at SSM/T and AMSU frequencies and scanning conditions. The land surface emissivities have been previously calculated for the globe from the SSM/I conical scanner between 19 and 85 GHz. To analyze the feasibility of deriving SSM/T and AMSU land surface emissivities from SSM/I emissivities, the spectral and angular variations of the emissivities are studied, with the help of ground-based measurements, models and satellite estimates. Up to 100 GHz, for snow and ice free areas, the SSM/T and AMSU emissivities can be derived with useful accuracy from the SSM/I emissivities- The emissivities can be linearly interpolated in frequency. Based on ground-based emissivity measurements of various surface types, a simple model is proposed to estimate SSM/T and AMSU emissivities for all zenith angles knowing only the emissivities for the vertical and horizontal polarizations at 53 deg zenith angle. The method is tested on the SSM/T-2 91.655 GHz channels. The mean difference between the SSM/T-2 and SSM/I-derived emissivities is less than or equal to 0.01 for all zenith angles with an r.m.s. difference of approx. = 0.02. Above 100 GHz, preliminary results are presented at 150 GHz, based on SSM/T-2 observations and are compared with the very few estimations available in the literature.

Prigent, Catherine↗

Earth Global Reference Atmospheric Model (GRAM) Overview and Comparison to Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2)

Earth-GRAM (Earth-Global Reference Atmospheric Model) Overview: Provides monthly statistics at any point in the atmosphere; Monthly, geographic, altitude variation; Current Version - Earth-GRAM 2016, C++; Output Includes - pressure, density, temperature, horizontal and vertical winds, speed of sound, atmospheric constituents; Used by engineering community to create atmospheric dispersions at a rapid runtime; Not a forecast model. MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2) Background: Developed by NASA Goddard Modelling and Assimilation Office; Horizontal Resolution: 0.625 degrees by 0.5 degrees longitude-by-latitude grid (NCEP (National Centers for Environmental Prediction) reanalysis I vs. 2.5 degrees by 2.5 degrees currently used in Earth-GRAM); Vertical resolution: 72 model layers or interpolated to 42 pressure levels to 0.1 hectopascals (hPa) (NCEP reanalysis I, surface to 10hPa at 17 pressure levels).

White, Patrick W.↗

High-Resolution Mesoscale Model Setup for the Eastern Range and Wallops Flight Facility

Mesoscale weather conditions can have an adverse effect on space launch, landing, ground processing, and weather advisories, watches, and warnings at the Eastern Range (ER) in Florida and Wallops Flight Facility (WFF) in Virginia. During summer, land-sea interactions across Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) lead to sea breeze front formation, which can spawn deep convection that can hinder operations and endanger personnel and resources. Many other weak locally-driven low-level boundaries and their interactions with the sea breeze front and each other can also initiate deep convection in the KSC/CCAFS area. These convective processes often last 60 minutes or less and pose a significant challenge to the local forecasters. Surface winds during the transition seasons (spring and fall) pose the most difficulties for the forecasters at WFF. They also encounter problems forecasting convective activity and temperature during those seasons. Therefore, accurate mesoscale model forecasts are needed to better forecast a variety of unique weather phenomena. Global and national scale models cannot properly resolve important local-scale weather features at each location due to their horizontal resolutions being much too coarse. Therefore, a properly tuned local data assimilation (DA) and forecast model at a high resolution is needed to provide improved capability. To accomplish this, a number of sensitivity tests were performed using the Weather Research and Forecasting (WRF) model in order to determine the best DA/model configuration for operational use at each of the space launch ranges to best predict winds, precipitation, and temperature. A set of Perl scripts to run the Gridpoint Statistical Interpolation (GSI)/WRF in real-time were provided by NASA's Short-term Prediction Research and Transition Center (SPoRT). The GSI can analyze many types of observational data including satellite, radar, and conventional data. The GSI/WRF scripts use a cycled GSI system similar to the operational North American Mesoscale (NAM) model. The scripts run a 12-hour pre-cycle in which data are assimilated from 12 hours prior up to the model initialization time. A number of different model configurations were tested for both the ER and WFF by varying the horizontal resolution on which the data assimilation was done. Three different grid configurations were run for the ER and two configurations were run for WFF for archive cases from 27 Aug 2013 through 10 Nov 2013. To quantify model performance, standard model output will be compared to the Meteorological Assimilation Data Ingest System (MADIS) data. The MADIS observation data will be compared to the WRF forecasts using the Model Evaluation Tools (MET) verification package. In addition, the National Centers for Environmental Prediction's Stage IV precipitation data will be used to validate the WRF precipitation forecasts. The author will summarize the relative skill of the various WRF configurations and how each configuration behaves relative to the others, as well as determine the best model configuration for each space launch range.

Watson, Leela R.↗

Assessment of Thematic Mapper Band-to-band Registration by the Block Correlation Method

The design of the Thematic Mapper (TM) multispectral radiometer makes it susceptible to band-to-band misregistration. To estimate band-to-band misregistration a block correlation method is employed. This method is chosen over other possible techniques (band differencing and flickering) because quantitative results are produced. The method correlates rectangular blocks of pixels from one band against blocks centered on identical pixels from a second band. The block pairs are shifted in pixel increments both vertically and horizontally with respect to each other and the correlation coefficient for each shift position is computed. The displacement corresponding to the maximum correlation is taken as the best estimate of registration error for each block pair. Subpixel shifts are estimated by a bi-quadratic interpolation of the correlation values surrounding the maximum correlation. To obtain statistical summaries for each band combination post processing of the block correlation results performed. The method results in estimates of registration error that are consistent with expectations.

Card, D. H.↗

Assimilation of Dual-Polarimetric Radar Observations with WRF GSI

Dual-polarimetric (dual-pol) radar typically transmits both horizontally and vertically polarized radio wave pulses. From the two different reflected power returns, more accurate estimate of liquid and solid cloud and precipitation can be provided. The upgrade of the traditional NWS WSR-88D radar to include dual-pol capabilities will soon be completed for the entire NEXRAD network. Therefore, the use of dual-pol radar network will have a broad impact in both research and operational communities. The assimilation of dual-pol radar data is especially challenging as few guidelines have been provided by previous research. It is our goal to examine how to best use dual-pol radar data to improve forecast of severe storm and forecast initialization. In recent years, the Development Testbed Center (DTC) has released the community Gridpoint Statistical Interpolation (GSI) DA system for the Weather Research and Forecasting (WRF) model. The community GSI system runs in independently environment, yet works functionally equivalent to operational centers. With collaboration with the NASA Short-term Prediction Research and Transition (SPoRT) Center, this study explores regional assimilation of the dual-pol radar variables from the WSR-88D radars for real case storms. Our presentation will highlight our recent effort on incorporating the horizontal reflectivity (ZH), differential reflectivity (ZDR), specific differential phase (KDP), and radial velocity (VR) data for initializing convective storms, with a significant focus being on an improved representation of hydrometeor fields. In addition, discussion will be provided on the development of enhanced assimilation procedures in the GSI system with respect to dual-pol variables. Beyond the dual-pol variable assimilation procedure developing within a GSI framework, highresolution (≤1 km) WRF model simulations and storm scale data assimilation experiments will be examined, emphasizing both model initialization and short-term forecast of precipitation fields and processes. Further details of the methodology of data assimilation, the impact of different dual-pol variables, the influence on precipitation forecast will be presented at the conference.

Li, Xuanli↗