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.
This report presents the Cryogenic Limb Array Etalon Spectrometer (CLAES) product improvement by use of the GSFC Data Assimilation System (DAS). The first task is to plug line of sight gradients derived from the CTM for 2/20/92 into the forward model of our retrieval software (RSW) in order to assess the impact on the retrieved quantities. The reporting period covers 12 May 2000 - 21 December 2000.
An important component of NASA's Program for Arctic Regional Climate Assessment (PARCA) is a mass balance investigation of the Greenland Ice Sheet. The mass balance is calculated by taking the difference between the areally Integrated snow accumulation and the net ice discharge of the ice sheet. Uncertainties in this calculation Include the snow accumulation rate, which has traditionally been determined by interpolating data from ice core samples taken from isolated spots across the ice sheet. The sparse data associated with ice cores juxtaposed against the high spatial and temporal resolution provided by remote sensing , has motivated scientists to investigate relationships between accumulation rate and microwave observations as an option for obtaining spatially contiguous estimates. The objective of this PARCA continuation proposal was to complete an estimate of surface accumulation rate on the Greenland Ice Sheet derived from C-band radar backscatter data compiled in the ERS-1 SAR mosaic of data acquired during, September-November, 1992. An empirical equation, based on elevation and latitude, is used to determine the mean annual temperature. We examine the influence of accumulation rate, and mean annual temperature on C-band radar backscatter using a forward model, which incorporates snow metamorphosis and radar backscatter components. Our model is run over a range of accumulation and temperature conditions. Based on the model results, we generate a look-up table, which uniquely maps the measured radar backscatter, and mean annual temperature to accumulation rate. Our results compare favorably with in situ accumulation rate measurements falling within our study area.
Our research focused on Kelvin waves in the tropical lower stratosphere, and resulted in three papers published or submitted to the Journal of Geophysical Research. The first of these, published in 2002, used temperature data from the Microwave Limb Sounder to examine the amplitude, frequency, phase, and spatial structure of leading modes of Kelvin waves. The second and third, submitted late in 2002 and currently in revision, described the response to Kelvin waves of various trace constituents measured by MLS and CLAES (methane, nitrous oxide, CFC-12, and ozone in the second paper; water vapor in the third paper). Water vapor is a special case because the vertical structure induced by Kelvin waves is convolved with water vapor's seasonally varying vertical profile induced by seasonal variations in temperature at the tropical tropopause. Forward modeling indicated that the vertical resolution of MLS was indeed adequate to capture this complicated structure, yet it was not visible in the MLS data, though the Kelvin wave signature was clear on certain UARS levels. The effects of Kelvin waves on the tropical tropopause and on stratosphere- troposphere exchange cannot be quantified from UARS data because of poor vertical resolution and sensitivity in that region. It is recommended that this analysis be repeated using data from the new MLS and HIRDLS instruments aboard Aura, and that priority be given to fine-scale retrievals of temperature, water vapor, and ozone in the tropical tropopause region.
The session "Lunar Remote Sensing: Fire, Ice, and Regolith" included the following:Compositional and Structural Study of the Aristarchus Plateau from Integrated UV-VIS-NIR Spectral Data; Clementine 2.7-?m Data: Mapping the Mare and Searching for Water; On the Search for Water at the Lunar Poles: Results of Forward Modeling of Permanently Shaded Areas and Lunar Prospector Measurements; Searching the Moon for Aluminous Mare Basalts Using Compositional Remote-Sensing Constraints I: Finding the Regions of Interest; Semi-automated Extraction of Contours from Lunar Topographic Maps; Basalts in Mare Humorum and S.E. Procellarum; The Hansteen and Helmet Volcanic Dome Regions on the Moon: Stratigraphy and Ages; Derivation of Elemental Abundance Maps at 15-km Spatial Resolution from the Merging of Clementine Optical and Lunar Prospector Geochemical Data; Remote Sensing and Geologic Studies of the Balmer Region of the Moon; Lava Flows in Mare Nubium and Mare Cognitum: A Geological History Based on Analysis of Multispectral Data; Development of Ground-based Lunar VIS/NEAR IR Spectral Imager; A BRDF Measurement Apparatus for Lab-based Samples; A New Source of High Resolution Lunar Images: Amateur Astronomers! ; Leakage of Gamma Rays and Neutrons from Thick Targets Bombarded by Energetic Protons; Progress on Reviving Lunar Orbiter: Scanning, Archiving, and Cartographic Processing at USGS; Modeling Lateral and Vertical Mixing by Impact Cratering with Applications for the Moon; Optical Maturity Study of Stuart#s Crater Candidate Impact; Evidence for Three Basins Beneath Oceanus Procellarum; and Ellipses of the South Pole-Aitken Basin: Implications for Basin Formation.
Hansteen Alpha is one of the few remaining locations on the Moon thought to be formed by highlands volcanism. Hansteen Alpha is a triangular shaped feature located in the southern portion of Oceanus Procellarum (12 degrees W, 50 degrees S) and its size is approximately 25 km on each side. As described by Hawke et al., there is clear evidence that: 1) Hansteen Alpha was emplaced by extrusive volcanic processes; and 2) it was formed by a viscous lava that should be enriched in Th. However, in the study of Hawke et al. using available Lunar Prospector (LP) Th data, it was concluded that the Hansteen Alpha region was not greatly enriched in Th as would be expected for a highly evolved, viscous lava. It was further concluded based on other compositional data that the magma that formed Hansteen Alpha did not correspond to any known rock type. Here we revisit the interpretation of Th abundances at Hansteen Alpha for a couple of reasons. First, the size of Hansteen Alpha is smaller than the spatial resolution of the LP Gamma-ray Spectrometer (LP-GRS) from which the Th abundances were derived. Therefore, the LP-GRS pixels covering Hansteen Alpha may not truly represent the Th abundance of the Hansteen Alpha feature. Second, recent work has led to a much greater understanding of the Th spatial distribution for small-area features on the lunar surface. In particular, using forward modeling techniques, we have developed the ability to obtain information about Th abundances for features that are at or smaller than the FWHM spatial resolution (approximately [80 square kilometers]) of the LP-GRS data.
Precipitation assimilation poses a special challenge in that the forward model for rain in a global forecast system is based on parameterized physics, which can have large systematic errors that must be rectified to use precipitation data effectively within a standard statistical analysis framework. We examine some key issues in precipitation assimilation and describe several exploratory studies in assimilating rainfall and latent heating information in NASA's global data assimilation systems using the forecast model as a weak constraint. We present results from two research activities. The first is the assimilation of surface rainfall data using a time-continuous variational assimilation based on a column model of the full moist physics. The second is the assimilation of convective and stratiform latent heating retrievals from microwave sensors using a variational technique with physical parameters in the moist physics schemes as a control variable. We will show the impact of assimilating these data on analyses and forecasts. Among the lessons learned are (1) that the time-continuous application of moisture/temperature tendency corrections to mitigate model deficiencies offers an effective strategy for assimilating precipitation information, and (2) that the model prognostic variables must be allowed to directly respond to an improved rain and latent heating field within an analysis cycle to reap the full benefit of assimilating precipitation information. of microwave radiances versus retrieval information in raining areas, and initial efforts in developing ensemble techniques such as Kalman filter/smoother for precipitation assimilation. Looking to the future, we discuss new research directions including the assimilation
A method of retrieving cloud microphysical properties using combined observations from both cloud radar and lidar is introduced. This retrieval makes use of an improvement to the traditional optimal estimation retrieval method, whereby a series of corrections are applied to the state vector during the search for an iterative solution. This allows faster convergence to a solution and is less processor intensive. The method is first applied to a synthetic cloud t o demonstrate its validity, and it is shown that the retrieval reliably reproduces vertical profiles of ice water content. The retrieval method is then applied to radar and lidar observations from the CRYSTAL-FACE experiment, and vertical profiles of ice crystal diameter, number concentration, and ice water content are retrieved for a cirrus cloud layers observed one day of that experiment. The validity of the relationship between visible extinction coefficient and radar reflectivity was examined. While synthetic tests showed such a functional relationship, the measured data only partially supported such a conclusion. This is due to errors in the forward model (as explained above) as well as errors in the data sets, including possible mismatch between lidar and radar profiles or errors in the optical depth. Empirical relationships between number concentrations and mean particle diameter were also examined. The results indicate that a distinct and robust relationship exists between these retrieved quantities and it is argued that such a relationship is more than an artifact of the retrieval process offering insight into the nature of the microphysical processes taking place in cirrus.
In operational weather forecasting, the assimilation of brightness temperatures from satellite sounders, instead of assimilation of 1D-retrievals has become increasingly common practice over the last two decades. Compared to these systems, assimilation of trace gases is still at a relatively early stage of development, and efforts to directly assimilate radiances instead of retrieved products have just begun a few years ago, partially because it requires much more computation power due to the employment of a radiative transport forward model (FM). This paper will focus on a method to assimilate SBUV/2 radiances (albedos) into the Global Earth Observation System Ozone Data Assimilation Scheme (GEOS-03DAS). While SBUV-type instruments cannot compete with newer sensors in terms of spectral and horizontal resolution, they feature a continuous data record back to 1978, which makes them very valuable for trend studies. Assimilation can help spreading their ground coverage over the whole globe, as has been previously demonstrated with the GEOS-03DAS using SBUV Version 6 ozone profiles. Now, the DAS has been updated to use the newly released SBUV Version 8 data. We will compare pre]lmlnarv results of SBUV radiance assimilation with the assimilation of retrieved ozone profiles, discuss methods to deal with the increased computational load, and try to assess the error characteristics and future potential of the new approach.
Ozone profiles are derived from backscattered radiances in the ultraviolet spectra (290-340 nm) measured by the nadir-viewing Global Ozone Monitoring Experiment using optimal estimation. Tropospheric O3 is directly retrieved with the tropopause as one of the retrieval levels. To optimize the retrieval and improve the fitting precision needed for tropospheric O3, we perform extensive wavelength and radiometric calibrations and improve forward model inputs. Retrieved O3 profiles and tropospheric O3 agree well with coincident ozonesonde measurements, and the integrated total O3 agrees very well with Earth Probe TOMS and Dobson/Brewer total O3. The global distribution of tropospheric O3 clearly shows the influences of biomass burning, convection, and air pollution, and is generally consistent with our current understanding.
Torque and plunge force during the initial plunge phase in Friction Stir Welding were measured for a 0.5 inch diameter pin entering a 2219 aluminum alloy plate. Weld structures were preserved for metallographic observation by making emergency stops at various plunge depths. The plunging pin tool is seen to be surrounded by a very fine grained layer of recrystallized metal extending substantially below the bottom of the pin, implying a shear interface in the metal below and not at the tool-metal interface. Torque and plunge force during the initial plunge phase in Friction Stir Welding are calculated from a straight forward model based on a concept to plastic flow in the vicinity of the plunging tool compatible with structural observations. The concept: a disk of weld metal seized to and rotating with the bottom of the pin is squeezed out laterally by the plunge force and extruded upwards in a hollow cylinder around the tool. As the shear surface separating rotating disk from stationary weld metal engulfs fresh metal, the fresh metal is subjected to severe shear deformation, which results in its recrystallization. Encouraging agreement between computations and measured torque and plunge force is obtained.
Recent improvements in the spatial and temporal retrieval of land-basecd cryospheric change information south of 42.5 detree allows fairly robust construction of forward model predictions of the time-rate of change in gravity.
Space-based measurements of reflected sunlight in the near-infrared (NIR) region promise to yield accurate and precise observations of the global distribution of atmospheric CO2. The Orbiting Carbon Observatory (OCO) is a future NASA mission, which will use this technique to measure the column-averaged dry air mole fraction of CO2 (XCO2) with the precision and accuracy needed to quantify CO2 sources and sinks on regional scales (approx.1000 x 1000 sq km and to characterize their variability on seasonal timescales. Here, we have used the OCO retrieval algorithm to retrieve XCO2 and surface pressure from space-based Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY) measurements and from coincident ground-based Fourier transform spectrometer (FTS) measurements of the O2 A band at 0.76 mm and the 1.58 mm CO2 band for Park Falls,Wisconsin. Even after accounting for a systematic error in our representation of the O2 absorption cross sections, we still obtained a positive bias between SCIAMACHY and FTS XCO2 retrievals of approx.3.5%. Additionally, the retrieved surface pressures from SCIAMACHY systematically underestimate measurements of a calibrated pressure sensor at the FTS site. These findings lead us to speculate about inadequacies in the forward model of our retrieval algorithm. By assuming a 1% intensity offset in the O2 A band region for the SCIAMACHY XCO2 retrieval, we significantly improved the spectral fit and achieved better consistency between SCIAMACHY and FTS XCO2 retrievals. We compared the seasonal cycle of XCO2 at Park Falls from SCIAMACHY and FTS retrievals with calculations of the Model of Atmospheric Transport and Chemistry/Carnegie-Ames-Stanford Approach (MATCH/CASA) and found a good qualitative agreement but with MATCH/CASA underestimating the measured seasonal amplitude. Furthermore, since SCIAMACHY observations are similar in viewing geometry and spectral range to those of OCO, this study represents an important test of the OCO retrieval algorithm and validation concept using NIR spectra measured from space. Finally, we argue that significant improvements in precision and accuracy could be obtained from a dedicated CO2 instrument such as OCO, which has much higher spectral and spatial resolutions than SCIAMACHY. These measurements would then provide critical data for improving our understanding of the carbon cycle and carbon sources and sinks.
The Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) radiance profile dataset of 1978/79 was reconditioned and reprocessed to Version 6 (V6) profiles of temperature and species that are improved significantly over those from Version 5 (V5). The LIMS V6 dataset was archived for public use in 2002. Improvements for its ozone include: (1) a more accurate accounting for instrument and spacecraft motion effects in the radiances, (2) the use of better spectroscopic line parameters for its ozone forward model, (3) retrievals of all its scans, (4) more accurate and compatible temperature versus pressure profiles (or T(p)), which are needed for the registration of the ozone radiances and for the removal of temperature effects from them, and (5) a better accounting for interfering species in the lower stratosphere. The retrieved V6 ozone profiles extend from near cloud top altitudes to about 80 km and from 64S to 84N latitude with better sampling along the orbit than for the V5 dataset. Calculated estimates of the single-profile precision and accuracy are provided for the V6 ozone from this study. Precision estimates based on the data themselves are of order 3% or better from 1 to 30 hPa. Estimates of total systematic error for a single profile are hard to generalize because the separate sources of error may not all be of the same sign and they depend somewhat on the atmospheric state. It is estimated that the V6 zonal mean ozone distributions are accurate to within 9% to 7% from 50 hPa to 3 hPa, respectively. Effects of a temperature bias can be significant and may be present at 1 to 2 hPa though. There may be ozone biases of order 10% at those levels due to possible biases of up to +2 K, but there is no indication of a similar problem elsewhere in the stratosphere. Simulation studies show that the LIMS retrievals are also underestimating slightly the small amplitudes of the atmospheric temperature tides, which affect its retrieved day/night ozone differences. There are small biases in the middle to lower stratosphere for the ascending versus descending node LIMS ozone, due principally to not accounting for the asymmetric weighting of its radiances across the tangent layer. The estimates of total accuracy were assessed by comparing the daily zonal mean LIMS ozone distributions against those from the Nimbus 7 SBUV Version 8 (V8) dataset for the same period. Generally, the LIMS V6 ozone agrees well with SBUV, except perhaps in the tropical lower stratosphere where the LIMS ozone is less. Still, the accuracy for LIMS V6 ozone in the lower stratosphere is improved over that found for LIMS V5, as indicated by several LIMS comparisons with ECC ozonesonde profiles. The LIMS V6 ozone is considered especially suitable for detailed studies of large-scale stratospheric processes above the 100-hPa level. Comparison of diurnal, photochemical model calculations with the monthly-averaged, upper stratospheric ozone observed from LIMS V6 indicates only a slight ozone deficit for the model at about 2 hPa. However, that deficit exhibits little to no seasonal variation and is in good agreement with similar model comparisons for a seasonal time series of ozone obtained with ground-based microwave instruments. Because the LIMS V6 ozone in the lower stratosphere has improved accuracy and sampling versus that of V5, it should now be possible to conduct quantitative studies of ozone transport and chemistry for the northern hemisphere, polar stratospheric winter of 1978/79 a time period when the catalytic loss of ozone due to reactive chlorine should not have been a major factor for the Arctic region.
The SOVA algorithm was originally developed under the Resilient Systems and Operations Project of the Engineering for Complex Systems Program from NASA s Aerospace Technology Enterprise as a conceptual framework to support real-time autonomous system mission and contingency management. The algorithm and its software implementation were formulated for generic application to autonomous flight vehicle systems, and its efficacy was demonstrated by simulation within the problem domain of Unmanned Aerial Vehicle autonomous flight management. The approach itself is based upon the precept that autonomous decision making for a very complex system can be made tractable by distillation of the system state to a manageable set of strategic objectives (e.g. maintain power margin, maintain mission timeline, and et cetera), which if attended to, will result in a favorable outcome. From any given starting point, the attainability of the end-states resulting from a set of candidate decisions is assessed by propagating a system model forward in time while qualitatively mapping simulated states into margins on strategic objectives using fuzzy inference systems. The expected return value of each candidate decision is evaluated as the product of the assigned value of the end-state with the assessed attainability of the end-state. The candidate decision yielding the highest expected return value is selected for implementation; thus, the approach provides a software framework for intelligent autonomous risk management. The name adopted for the technique incorporates its essential elements: Strategic Objective Valuation and Attainability (SOVA). Maximum value of the approach is realized for systems where human intervention is unavailable in the timeframe within which critical control decisions must be made. The Far Ultraviolet Spectroscopic Explorer (FUSE) satellite, launched in 1999, has been collecting science data for eight years.[1] At its beginning of life, FUSE had six gyros in two IRUs and four reaction wheels. Over time through various failures, the satellite has been left with one reaction wheel on the vehicle skew axis and two gyros. To remain operational, a control scheme has been implemented using the magnetic torque rods and the remaining momentum wheel.[2] As a consequence, there are attitude regions where there is insufficient torque authority to overcome environmental disturbances (e.g. gravity gradient torques). The situation is further complicated by the fact that these attitude regions shift inertially with time as the spacecraft moves through earth s magnetic field during the course of its orbit. Under these conditions, the burden of planning targets and target-to-target slew maneuvers has increased significantly since the beginning of the mission.[3] Individual targets must be selected so that the magnetic field remains roughly aligned with the skew wheel axis to provide enough control authority to the other two orthogonal axes. If the field moves too far away from the skew axis, the lack of control authority allows environmental torques to pull the satellite away from the target and can potentially cause it to tumble. Slew maneuver planning must factor the stability of targets at the beginning and end, and the torque authority at all points along the slew. Due to the time varying magnetic field geometry relative to any two inertial targets, small modifications in slew maneuver timing can make large differences in the achievability of a maneuver.
Watersheds draining into the Gulf of Alaska (GoA) experience large seasonal and inter-annual variations of water in the form of rain, snow, and ice, but accurate constraints on these variations have been difficult to obtain. Over larger geographic regions, water variations can be inferred directly from the Gravity Recovery and Climate Experiment (GRACE) data. However, because GoA variations occur over such a small region, the inferred average value of water flux increases as the applied smoothing of the GRACE data decreases. We use this observed scaling together with scaling results obtained from forward models to infer a seasonal amplitude of 115 plus or minus 20 cubic kilometers of water and an average contribution to sea level rise over the two years of data of 0.31 plus or minus 0.09 millimeters per year. These results suggest that accelerated melting that began in the late 1990s, as inferred from altimetry, continues unabated.
The quality of the retrieved temperature-versus-pressure (or T(p)) profiles is described for the middle atmosphere for the publicly available Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) Version 1.07 (V1.07) data set. The primary sources of systematic error for the SABER results below about 70 km are (1) errors in the measured radiances, (2) biases in the forward model, and (3) uncertainties in the corrections for ozone and in the determination of the reference pressure for the retrieved profiles. Comparisons with other correlative data sets indicate that SABER T(p) is too high by 1-3 K in the lower stratosphere but then too low by 1 K near the stratopause and by 2 K in the middle mesosphere. There is little difference between the local thermodynamic equilibrium (LTE) algorithm results below about 70 km from V1.07 and V1.06, but there are substantial improvements/differences for the non-LTE results of V1.07 for the upper mesosphere and lower thermosphere (UMLT) region. In particular, the V1.07 algorithm uses monthly, diurnally averaged CO2 profiles versus latitude from the Whole Atmosphere Community Climate Model. This change has improved the consistency of the character of the tides in its kinetic temperature (T(sub k)). The T(sub k) profiles agree with UMLT values obtained from ground-based measurements of column-averaged OH and O2 emissions and of the Na lidar returns, at least within their mutual uncertainties. SABER T(sub k) values obtained near the mesopause with its daytime algorithm also agree well with the falling sphere climatology at high northern latitudes in summer. It is concluded that the SABER data set can be the basis for improved, diurnal-to-interannual-scale temperatures for the middle atmosphere and especially for its UMLT region.
A computer program performs level-2 processing of thermal-microwave-radiance data from observations of the limb of the Earth by the Earth Observing System (EOS) Microwave Limb Sounder (MLS). The purpose of the processing is to estimate the composition and temperature of the atmosphere versus altitude from .8 to .90 km. "Level-2" as used here is a specialists f term signifying both vertical profiles of geophysical parameters along the measurement track of the instrument and processing performed by this or other software to generate such profiles. Designed to be flexible, the program is controlled via a configuration file that defines all aspects of processing, including contents of state and measurement vectors, configurations of forward models, measurement and calibration data to be read, and the manner of inverting the models to obtain the desired estimates. The program can operate in a parallel form in which one instance of the program acts a master, coordinating the work of multiple slave instances on a cluster of computers, each slave operating on a portion of the data. Optionally, the configuration file can be made to instruct the software to produce files of simulated radiances based on state vectors formed from sets of geophysical data-product files taken as input.
The gravity wave (GW)-resolving capabilities of 118-GHz saturated thermal radiances acquired throughout the stratosphere by the Microwave Limb Sounder (MLS) on the Aura satellite are investigated and initial results presented. Because the saturated (optically thick) radiances resolve GW perturbations from a given altitude at different horizontal locations, variances are evaluated at 12 pressure altitudes between 21 and 51 km using the 40 saturated radiances found at the bottom of each limb scan. Forward modeling simulations show that these variances are controlled mostly by GWs with vertical wavelengths z 5 km and horizontal along-track wavelengths of y 100-200 km. The tilted cigar-shaped three-dimensional weighting functions yield highly selective responses to GWs of high intrinsic frequency that propagate toward the instrument. The latter property is used to infer the net meridional component of GW propagation by differencing the variances acquired from ascending (A) and descending (D) orbits. Because of improved vertical resolution and sensitivity, Aura MLS GW variances are 5?8 times larger than those from the Upper Atmosphere Research Satellite (UARS) MLS. Like UARS MLS variances, monthly-mean Aura MLS variances in January and July 2005 are enhanced when local background wind speeds are large, due largely to GW visibility effects. Zonal asymmetries in variance maps reveal enhanced GW activity at high latitudes due to forcing by flow over major mountain ranges and at tropical and subtropical latitudes due to enhanced deep convective generation as inferred from contemporaneous MLS cloud-ice data. At 21-28-km altitude (heights not measured by the UARS MLS), GW variance in the tropics is systematically enhanced and shows clear variations with the phase of the quasi-biennial oscillation, in general agreement with GW temperature variances derived from radiosonde, rocketsonde, and limb-scan vertical profiles.