Search NASA⌕ Search

SEARCH · Search NASA

Results for “height function”

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 991 records · Page 55

Biomechanical Modeling Analysis of Loads Configuration for Squat Exercise

INTRODUCTION: Long duration space travel will expose astronauts to extended periods of reduced gravity. Since gravity is not present to assist loading, astronauts will use resistive and aerobic exercise regimes for the duration of the space flight to minimize loss of bone density, muscle mass and aerobic capacity that occurs during exposure to a reduced gravity environment. Unlike the International Space Station (ISS), the area available for an exercise device in the next generation of spacecraft for travel to the Moon or to Mars is limited and therefore compact resistance exercise device prototypes are being developed. The Advanced Resistive Exercise Device (ARED) currently on the ISS is being used as a benchmark for the functional performance of these new devices. Biomechanical data collection and computational modeling aid the device design process by quantifying the joint torques and the musculoskeletal forces that occur during exercises performed on the prototype devices. METHODS The computational models currently under development utilize the OpenSim [1] software platform, consisting of open source code for musculoskeletal modeling, using biomechanical input data from test subjects for estimation of muscle and joint loads. The OpenSim Full Body Model [2] is used for all analyses. The model incorporates simplified wrap surfaces, a new knee model and updated lower body muscle parameters derived from cadaver measurements and magnetic resonance imaging of young adults. The upper body uses torque actuators at the lumbar and extremity joints. The test subjects who volunteer for this study are instrumented with reflective markers for motion capture data collection while performing squat exercising on the Hybrid Ultimate Lifting Kit (HULK) prototype device (ZIN Technologies, Middleburg Heights, OH). Ground reaction force data is collected with force plates under the feet, and device loading is recorded through load cells internal to the HULK. Test variables include the applied device load and the dual cable long bar or single cable T-bar interface between the test subject and the device. Data is also obtained using free weights with the identical loading for a comparison to the resistively loaded exercise device trials. The data drives the OpenSim biomechanical model, which has been scaled to match the anthropometrics of the test subject, to calculate the body loads. RESULTS Lower body kinematics, joint moments, joint forces and muscle forces are obtained from the OpenSim biomechanical analysis of the squat exercises under different loading conditions. Preliminary results from the model for the loading conditions will be presented as will hypotheses developed for follow on work.

Human Factors Engineering↗

An error covariance model for sea surface topography and velocity derived from TOPEX/POSEIDON altimetry

In order to facilitate the use of satellite-derived sea surface topography and velocity oceanographic models, methodology is presented for deriving the total error covariance and its geographic distribution from TOPEX/POSEIDON measurements. The model is formulated using a parametric model fit to the altimeter range observations. The topography and velocity modeled with spherical harmonic expansions whose coefficients are found through optimal adjustment to the altimeter range residuals using Bayesian statistics. All other parameters, including the orbit, geoid, surface models, and range corrections are provided as unadjusted parameters. The maximum likelihood estimates and errors are derived from the probability density function of the altimeter range residuals conditioned with a priori information. Estimates of model errors for the unadjusted parameters are obtained from the TOPEX/POSEIDON postlaunch verification results and the error covariances for the orbit and the geoid, except for the ocean tides. The error in the ocean tides is modeled, first, as the difference between two global tide models and, second, as the correction to the present tide model, the correction derived from the TOPEX/POSEIDON data. A formal error covariance propagation scheme is used to derive the total error. Our global total error estimate for the TOPEX/POSEIDON topography relative to the geoid for one 10-day period is found tio be 11 cm RMS. When the error in the geoid is removed, thereby providing an estimate of the time dependent error, the uncertainty in the topography is 3.5 cm root mean square (RMS). This level of accuracy is consistent with direct comparisons of TOPEX/POSEIDON altimeter heights with tide gauge measurements at 28 stations. In addition, the error correlation length scales are derived globally in both east-west and north-south directions, which should prove useful for data assimilation. The largest error correlation length scales are found in the tropics. Errors in the velocity field are smallest in midlatitude regions. For both variables the largest errors caused by uncertainty in the geoid. More accurate representations of the geoid await a dedicated geopotential satellite mission. Substantial improvements in the accuracy of ocean tide models are expected in the very near future from research with TOPEX/POSEIDON data.

Tsaoussi, Lucia S.↗

Light Scattered from Polished Optical Surfaces: Wings of the Point Spread Function

Random figure errors from the polishing process plus particles on the main mirrors in a telescope cause an extended point spread function (PSF) declining approximately as the inverse square of the sine of the angle from a star from about 100 micro-rad to a right angle. The decline in at least one case, and probably in general, proceeds as the inverse cube at smaller angles where the usual focal plane aperture radius is chosen. The photometric error due to misalignment by one Airy ring spacing with an aperture of n rings depends on the net variance in the figure. It is approximately 60/(n+1)(3) when using the data of Kormendy (1973). A typical value is 6 x 10 to the -5th power per ring of misalignment with n = 100 rings. The encircled power may be modulated on a time scale of hours by parts per thousand in a wavelength dependent manner due to relative humidity effects on mirror dust. The scattering according to an inverse power law is due to a random walk in aberration height caused by a multitude of facets and slope errors left by the polishing process. A deviation from such a law at grazing emergence may permit monitoring the dust effects.

Kenknight, C. E.↗

A Microwave Radiometric Method to Obtain the Average Path Profile of Atmospheric Temperature and Humidity Structure Parameters and Its Application to Optical Propagation System Assessment

The values of the key atmospheric propagation parameters Ct2, Cq2, and Ctq are highly dependent upon the vertical height within the atmosphere thus making it necessary to specify profiles of these values along the atmospheric propagation path. The remote sensing method suggested and described in this work makes use of a rapidly integrating microwave profiling radiometer to capture profiles of temperature and humidity through the atmosphere. The integration times of currently available profiling radiometers are such that they are approaching the temporal intervals over which one can possibly make meaningful assessments of these key atmospheric parameters. Since these parameters are fundamental to all propagation conditions, they can be used to obtain Cn2 profiles for any frequency, including those for an optical propagation path. In this case the important performance parameters of the prevailing isoplanatic angle and Greenwood frequency can be obtained. The integration times are such that Kolmogorov turbulence theory and the Taylor frozen-flow hypothesis must be transcended. Appropriate modifications to these classical approaches are derived from first principles and an expression for the structure functions are obtained. The theory is then applied to an experimental scenario and shows very good results.

Atmospheric turbulence↗

Residual Strength Analyses of Monolithic Structures

Finite-element fracture simulation methodology predicts the residual strength of damaged aircraft structures. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic- plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intention of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels.In this regard, a series of fracture tests were conducted on both flat and curved aluminum alloy integrally-stiffened panels. These flat panels were subjected to uniaxial tension and during the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (wc) using three-dimensional code (ZIP3D) and the plane-strain core height (hJ using two-dimensional code (STAGS). These values were then used in the STAGS analysis to predict the fracture behavior of the integrally-stiffened panels. The analyses modeled stable tearing, buckling, and crack branching at the integral stiffener using different values of critical CTOA for different material thicknesses and orientation. Comparisons were made between measured and predicted load-crack extension, out-of-plane displacements and local deformations around the crack tip region. Simultaneously, three-dimensional capabilities to model crack branching and to monitor stable crack growth of multiple cracks in a large thick integrally-stiffened flat panels were implemented in three-dimensional finite element code (ZIP3D) and tested by analyzing the integrally-stiffened panels tested at Alcoa. The residual strength of the panels predicted from STAGS and ZP3D code compared very well with experimental data. In recent times, STAGS software has been updated with new features and now one can have combinations of solid and shell elements in the residual strength analysis of integrally-stiffened panels.

Forth, Scott↗

Ames Optimized TCA Configuration

Configuration design at Ames was carried out with the SYN87-SB (single block) Euler code using a 193 x 49 x 65 C-H grid. The Euler solver is coupled to the constrained (NPSOL) and the unconstrained (QNMDIF) optimization packages. Since the single block grid is able to model only wing-body configurations, the nacelle/diverter effects were included in the optimization process by SYN87's option to superimpose the nacelle/diverter interference pressures on the wing. These interference pressures were calculated using the AIRPLANE code. AIRPLANE is an Euler solver that uses a unstructured tetrahedral mesh and is capable of computations about arbitrary complete configurations. In addition, the buoyancy effects of the nacelle/diverters were also included in the design process by imposing the pressure field obtained during the design process onto the triangulated surfaces of the nacelle/diverter mesh generated by AIRPLANE. The interference pressures and nacelle buoyancy effects are added to the final forces after each flow field calculation. Full details of the (recently enhanced) ghost nacelle capability are given in a related talk. The pseudo nacelle corrections were greatly improved during this design cycle. During the Ref H and Cycle 1 design activities, the nacelles were only translated and pitched. In the cycle 2 design effort the nacelles can translate vertically, and pitch to accommodate the changes in the lower surface geometry. The diverter heights (between their leading and trailing edges) were modified during design as the shape of the lower wing changed, with the drag of the diverter changing accordingly. Both adjoint and finite difference gradients were used during optimization. The adjoint-based gradients were found to give good direction in the design space for configurations near the starting point, but as the design approached a minimum, the finite difference gradients were found to be more accurate. Use of finite difference gradients was limited by the CPU time limit available on the Cray machines. A typical optimization run using finite difference gradients can use only 30 to 40 design variables and one optimization iteration within the 8 hour queue limit for the chosen grid size and convergence level. The efficiency afforded by the adjoint method allowed for 50-120 design variables and 5-10 optimization iterations in the 8 hour queue. Geometric perturbations to the wing and fuselage were made using the Hicks/Henne (HH) shape functions. The HH functions were distributed uniformly along the chords of the wing defining sections and lofted linearly. During single-surface design, constraints on thickness and volume at selected wing stations were imposed. Both fuselage camber and cross-sectional area distributions were permitted to change during design. The major disadvantage to the use of these functions is the inherent surface waviness produced by repeated use of such functions. Many smoothing operations were required following optimization runs to produce a configuration with reasonable smoothness. Wagner functions were also used on the wing sections but were never used on the fuselage. The Wagner functions are a family of increasingly oscillatory functions that have also been used extensively in airfoil design. The leading and trailing edge regions of the wing were designed by use of polynomial and monomial functions respectively. Twist was attempted but was abandoned because of little performance improvement available from changing the baseline twist.

Cliff, Susan E.↗

Topography of the Moon from the Clementine Lidar

Range measurements from the lidar instrument carried aboard the Clementine spacecraft have been used to produce an accurate global topographic model of the Moon. This paper discusses the function of the lidar; the acquisition, processing, and filtering of observations to produce a global topographic model; and the determination of parameters that define the fundamental shape of the Moon. Our topographic model: a 72nd degree and order spherical harmonic expansion of lunar radii, is designated Goddard Lunar Topography Model 2 (GLTM 2). This topographic field has an absolute vertical accuracy of approximately 100 m and a spatial resolution of 2.5 deg. The field shows that the Moon can be described as a sphere with maximum positive and negative deviations of approx. 8 km, both occurring on the farside, in the areas of the Korolev and South Pole-Aitken (S.P.-Aitken) basins. The amplitude spectrum of the topography shows more power at longer wavelengths as compared to previous models, owing to more complete sampling of the surface, particularly the farside. A comparison of elevations derived from the Clementine lidar to control point elevations from the Apollo laser altimeters indicates that measured relative topographic heights generally agree to within approx. 200 in over the maria. While the major axis of the lunar gravity field is aligned in the Earth-Moon direction, the major axis of topography is displaced from this line by approximately 10 deg to the cast and intersects the farside 24 deg north of the equator. The magnitude of impact basin topography is greater than the lunar flattening (approx. 2 km) and equatorial ellipticity (approx. 800 m), which imposes a significant challenge to interpreting the lunar figure. The floors of mare basins are shown to lie close to an equipotential surface, while the floors of unflooded large basins, except for S.P.-Aitken, lie above this equipotential. The radii of basin floors are thus consistent with a hydrostatic mechanism for the absence of significant farside maria except for S.P.-Aitken, whose depth and lack of mare require significant internal compositional and/or thermal heterogeneity. A macroscale surface roughness map shows that roughness at length scales of 10(exp 1) - 10(exp 2) km correlates with elevation and surface age.

Smith, David E.↗

Residual Strength Characterization of a Curved Integrally-Stiffened Panel

Over the years, Finite-element fracture simulation methodology has been very well established at NASA Langley to predict the residual strength of damaged aircraft structures. This methodology has been experimentally verified at NASA Langley for structures ranging from laboratory coupons up to full-scale built-up structural components with single and multiple-site damage cracking. The methodology uses the critical crack-tip-opening-angle (CTOA) fracture criterion to characterize the fracture behavior of the material. The CTOA fracture criterion assumes that stable crack growth occurs when the crack-tip angle reaches a constant critical value. The use of the CTOA criterion requires an elastic-plastic, finite-element analysis. The critical CTOA value is determined by simulating fracture behavior in laboratory specimens, such as a compact specimen, to obtain the angle that best fits the observed test behavior. The critical CTOA value appears to be independent of loading, crack length, and in-plane dimensions. However, it is a function of material thickness and local crack-front constraint. Modeling the local constraint requires either a three-dimensional analysis or a two-dimensional analysis with an approximation to account for the constraint effects. In recent times as the aircraft industry is leaning towards monolithic structures with the intension of reducing part count and manufacturing cost, there has been a consistent effort at NASA Langley to extend critical CTOA based numerical methodology in the analysis of integrally-stiffened panels. In this regard, a series of fracture tests were conducted on curved aluminum-alloy integrally-stiffened panels. These curved panels were subjected to uniaxial tension and pressure loading. During the test, applied load-crack extension, out-of-plane displacements and local deformations around the crack tip region were measured. Compact and middle-crack tension specimens were tested to determine the critical angle (psi(sub c) using three-dimensional code (ZIP3D) and the plane-strain core height (h(sub c)) using two-dimensional code (STAGS).

Seshadri, B. R.↗

Full-scale Wind-tunnel and Flight Tests of a Fairchild 22 Airplane Equipped with External-airfoil Flaps

Wind-tunnel and flight tests have been made of a Fairchild 22 airplane equipped with a wing having external-airfoil flaps that also perform the function of ailerons. Lift, drag, and pitching-moment coefficients of the airplane with several flap settings, and the rolling- and yawing-moment coefficients with the flaps deflected as ailerons were measured in the full-scale tunnel with the horizontal tail surfaces and propeller removed. The effect of the flaps on the low speed and on the take-off and landing characteristics, the effectiveness of flaps when used as ailerons, and the forces required to operate them as ailerons were determined in flight. The wind-tunnel tests showed that the flaps increased the maximum lift coefficient of the airplane from 1.51 with the flap in the minimum drag position to 2.12 with the flap in the minimum drag position to 2.12 with the flap deflected 30 degrees. In the flight tests the minimum speed decreased from 46.8 miles per hour with the flaps up to 41.3 miles per hour with the flaps deflected. The required take-off run to attain a height of 50 feet was reduced from 820 to 750 feet and the landing run from a height of 50 feet was reduced from 930 to 480 feet. The flaps for this installation gave lateral control that was not entirely satisfactory. Their rolling action was good but the adverse yaw resulting from their use was greater than is considerable, and the stick forces required to operate them increased too rapidly with speed.

Reed, Warren D↗

How Did I Get Here?

What country was responsible for much of the early work in laser altimetry, who was the first scientist to look at a vegetation profile generated by a laser, who first related laser ranges to forest structure? The purpose of this paper is to look back 3+ decades ago when forestry lidar was in its formative years and to give credit to those who first thought to employ a laser to measure the Earth's surface. As with much scientific research, advances were often made independently and concurrently in a number of countries. Functioning lasers were first demonstrated in 1960 in the USA and in 1961 in the USSR, but research into the use of lasers as forest measurement tools did not begin for another IS years. Initially, with respect to Earth resources, lasers were employed to measure sea ice surface roughness, to make near-shore bathymetric measurements, to penetrate forests to make detailed topographic measurements, and to fluoresce oceanic phytoplankton for surface current studies. Some of these early studies noted that forest profiles were evident but in fact added noise to topographic retrievals. Trees became the signal rather than noise in 1976 when researchers in the Canadian Forest Service went to an IUFRO conference in Oslo and reported on efforts to better estimate timber volume in tropical forests. They employed large-scale airphotos (1:500 up to 1:4000) to measure the top of the tropical canopy, an avionic radar to measure the location of the ground beneath the canopy, and a barometric sensor to record aircraft/radar height above a datum (e.g., sea level). Their radar did not work well in dense vegetation and mentioned in passing that a laser altimeter had worked well in Canadian forests. But in a report 2 years later they reported that, on a study area in Costa Rica, the radar ground line was much improved and no mention is made of a lidar altimeter. In the USSR in 1977, Russian researchers felled a birch and a spruce, aimed a helium-neon, 0.63 micron laser with a spot size of approx 25mm, at the horizontal trees, produced a profilograph, compared it to tape measurements, and concluded that, with increased power, such a laser could be mounted on an aircraft to remotely measure forest canopies. In 1979, they mounted their He-Ne laser on an AN-2 biplane and acquired their fIrSt airborne profiles. Scientists with TRANARG, a mapping/surveying company in Caracas, Venezuela, reported on 1976 flights employing a helium-neon lidar to collect over 11,000 km of lidar profiles spaced 1.5 km apart in order to construct a topographic map to help site a new reservoir. They noted 35-40 m median canopy heights with emergents up to 55 m in their profiles. These studies and others that utilize these height and canopy density metrics for forest mensuration are reviewed in this glance backwards at the history of forestry lidar.

Nelson, Ross↗

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). V. First Results on the Magnetic Field Orientation of Galaxies

We present the analysis of the magnetic field ( B -field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA of 14 nearby (<20 Mpc) galaxies with resolved (5"-18" ; 90 pc–1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 μm . We compute the magnetic pitch-angle (𝜓 B ) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter ( ζ ) to estimate the disordered-to-ordered ratio of spiral B -fields. We find FIR and radio wavelengths to not generally trace the same B -field morphology in galaxies. The 𝜓 B profiles tend to be more ordered across all galactocentric radii in radio ( ζ 6cm = 0.93 + 0.03) than in FIR ( ζ 154μm = 0.84 + 0.14). For spiral galaxies, FIR B -fields are 2%–75% more turbulent than the radio B -fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B -fields along the galactic outflows than radio polarization. Our results suggest that the B -fields associated with dense, dusty, turbulent star-forming regions (those traced at FIR) are less ordered than warmer, less dense regions (those traced at radio) of the interstellar medium. The FIR B -fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of a few hundred parsecs in the disk of spiral galaxies than the radio B -fields.

magnetic field (B-field)↗

Extragalactic Magnetism with SOFIA (SALSA Legacy Program) - V. First Results on the Magnetic Field Orientation of Galaxies

We present the analysis of the magnetic field (B-field) structure of galaxies measured with far-infrared (FIR) and radio (3 and 6 cm) polarimetric observations. We use the first data release of the Survey of extragALactic magnetiSm with SOFIA (SALSA) of 14 nearby (< 20 Mpc) galaxies with resolved (5′′ − 18′′; 90 pc–1 kpc) imaging polarimetric observations using SOFIA/HAWC+ from 53 to 214 μm. We compute the magnetic pitch angle (Ψ B ) profiles as a function of the galactocentric radius. We introduce a new magnetic alignment parameter (ζ) to estimate the disordered-to-ordered ratio of spiral B -fields. We find FIR and radio wavelengths to not generally trace the same B -field morphology in galaxies. The Ψ B profiles tend to be more ordered across all galactocentric radii in radio (ζ 6cm = 0.93 ± 0.03) than in FIR (ζ 154μm = 0.84 ± 0.14). For spiral galaxies, FIR B-fields are 2−75% more turbulent than the radio B -fields. For starburst galaxies, we find that FIR polarization is a better tracer of the B -fields along the galactic outflows than radio polarization. Our results suggest that the B -fields associated with dense, dusty, turbulent star-forming regions, those traced at FIR, are less ordered than warmer, less-dense regions, those traced at radio, of the interstellar medium. The FIR B -fields seem to be more sensitive to the activity of the star-forming regions and molecular clouds within a vertical height of few hundred pc in the disk of spiral galaxies than the radio B -fields.

SALSA↗

Implications of Large Scale Shifts in Tropospheric NO(x) Levels in the Remote Tropical Pacific

A major observation recorded during NASA's western Pacific Exploratory Mission (PEM-West B) was the large shift in tropical NO levels as a function of geographical location. High-altitude NO levels exceeding 100 pptv were observed during portions of tropical flights 5-8, while values almost never exceeded 20 pptv during tropical flights 9 and 10. The geographical regions encompassing these two flight groupings are here labeled "high" and "low" NOx regimes. A comparison of these two regimes, based on back trajectories and chemical tracers, suggests that air parcels in both were strongly influenced by deep convection. The low NO(x) regime appears to have been predominantly impacted by marine convection, whereas the high NO(x) regime shows evidence of having been more influenced by deep convection over a continental land mass. DMSP satellite observations point strongly toward lightning as the major source of NOx in the latter regime. Photochemical ozone formation in the high NOx regime exceeded that for low NO(x) by factors of 2 to 6, whereas O3 destruction in the low NO(x)regime exceeded that for high NOx by factors of up to 3. Taking the tropopause height to be 17 km, estimates of the net photochemical effect on the O3 column revealed that the high NO(x) regime led to a small net production. By contrast, the low NOx regime was shown to destroy O3 at the rate of 3.4 % per day. One proposed mechanism for off-setting this projected large deficit would involve the transport of 03 rich midlatitude air into the tropics. Alternatively, it is suggested that O3 within the tropics may be overall near self-sustaining with respect to photochemical activity. This scenario would require that some tropical regions, unsampled at the time of PEM-B, display significant net column O3 production, leading to an overall balanced budget for the "greater" tropical Pacific basin. Details concerning the chemical nature of such regimes are discussed.

Crawford, J. H.↗

Regional Variability in Convection and Rain Retrievals from the TRMM Microwave Imager (TMI)

Precipitation Radar (PR) on board the TRMM satellite shows that the average height of 30 dBz in convective rain areas of the tropics varies significantly from one region to the other. When the convection is weak this height is shallow and when convection is strong this height extends deeper into the troposphere. The brightness temperature (Tb) measured by the microwave radiometer by itself does not reflect this nature of convection satisfactorily. Radiative transfer simulations of Tbs reveal that this could be due to the variations in the vertical distribution of optically active water and ice hydrometeors and their density, shape, and size. These variations are not coupled uniquely to the strength of the convective updrafts, and as a result the Tbs do not reflect properly the convective strength indicated by PR. Because of this deficiency in the Tbs the rain rate deduced from them differs from that of PR. For this reason, to improve the estimation of rain rate we have developed an empirical method. In this method a parameter based on the areal extent of the Tbs that exceed a certain magnitude is included along with the Tbs. Rain rate deduced with this approach is better correlated with that of PR when compared to the current Version 5 operational algorithm. The percentage of rain volume as a function of rain rate, for a given region of 5deg lat. X 5deg long. over a period of three months, deduced from this method, is also in better agreement with that of the PR.

Prabhakara, C.↗

GPS Altimetry

The advent of satellite altimetry has greatly improved our ability to observe global ocean circulation. However, the swath of a single, nadir-viewing satellite altimeter is only a few km and the track spacing is several hundred km to resolve the two-dimensional structure of ocean eddies. Our goal is to increase spatial and temporal coverage by monitoring Global Positioning System (GPS) signals reflected from the ocean. A constellation of spacecraft would each carry a GPS receiver capable of recording 8 reflections simultaneously. The reflections are well distributed in azimuth and elevation and can be tracked continuously while the satellite is in view, and another is then acquired, as illustrated below. The diagram depicts a new approach at altimetry measurements where ocean surface reflected GPS signals are simultaneously tracked and processed in a GPS flight receiver in space. The reflected GPS signals from the ocean must be compared precisely with the direct GPS signals in order to infer the characteristics of the ocean from the combined data set. Understanding the features and accuracy of GPS altimetry measurement is crucial to establishing its suitability for oceanography. Preliminary work has enabled us to theoretically model the signal output of the correlator for a variety of system parameters such as wind speed (sea roughness), receiver height, incidence angle, receiver range and Doppler filter bandwidth and antenna gain. Expected signal-to-noise ratio has been estimated from which we have inferred, to a first approximation, the basic receiver gain requirements for a space-based altimeter and the expected range raw error. In 1998, work on a different task led to the extraction of the first reflected GPS signal observed from a spaceborne receiver during the 1995 Space Transportation System-68 (STS-68) Shuttle Radar Laboratory-2 (SRL-2) high resolution synthetic aperture radar mission. Good comparisons with our signal models have been obtained. Having established that only modest signal-to-noise ratios are obtainable unless very high gain antennas are used, we expect that single measurements might not provide an estimate of sea state parameters as accurate as that obtainable with traditional remote sensing instruments, if costs are limited. Therefore spatial and temporal averaging of many measurements is required. Since the receiving satellite tracks do not repeat, measurements in a given area will be used to refine the local solution as a function of time, and will define the spatial resolution. In order to make this instrument viable in space, we need to detect and process many scattered signals. In 1998 we made a first step at understanding what configurations of antenna gain and orientations capture the largest number of viable signals to be used in the subsequent spatial/temporal averaging process. Additional information is contained in the original.

Lichten, Stephen M.↗

Towards an Optimal Estimation Retrieval of Cirrus Cloud Optical and Microphysical Properties Using Hyperspectral Shortwave Instruments and A Fast Radiative Transfer Algorithm

Cirrus cloud retrieval products (here, cloud optical depth, effective particle size, and cloud top height) are important inputs into numerical weather and climate models. Uncertainties in such retrieval products, as a matter of course, propagate downstream, impacting model calculations. Improvements in high-quality global cirrus cloud optical and microphysical data products from satellite observations are needed to understand and reduce retrieval uncertainties. Hyperspectral instruments produce high-resolution and information-dense radiance spectra, thus offering the opportunity to reduce uncertainties in retrieval products. An optimal estimation-based cirrus cloud optical and microphysical product retrieval is under development for the NASA Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) and the Earth Surface Mineral Dust Source Investigation (EMIT) instruments, and the forthcoming Climate Absolute Radiance and Refractivity Observatory Pathfinder (CLARREO-Pathfinder). The spectral coverage of all three instruments includes the ultraviolet, visible, and near-infrared. In this study a very fast radiative transfer model, the Principal Component-based Radiative Transfer Model in the solar spectral region (PCRTM-Solar), is used for forward modeling computations. This will reduce the computational burden in the forward radiance and Jacobian calculations as the cost function is minimized, enabling the entire AVIRIS, EMIT, and CLARREO-Pathfinder spectral range to be used in the optimal estimation framework. Using the entire spectrum will maximize the information content, resulting in a more robust and more accurate retrieval. As a first step, the retrieval is being designed for single layer ice clouds over open ocean water. Preliminary results will be shown.

Jeffrey Mast↗

A Method for Improving Hotspot Directional Signatures in BRDF Models Used for MODIS

The semi-empirical, kernel-driven, linear RossThick-LiSparseReciprocal (RTLSR) Bidirectional Reflectance Distribution Function (BRDF) model is used to generate the routine MODIS BRDFAlbedo product due to its global applicability and the underlying physics. A challenge of this model in regard to surface reflectance anisotropy effects comes from its underestimation of the directional reflectance signatures near the Sun illumination direction; also known as the hotspot effect. In this study, a method has been developed for improving the ability of the RTLSR model to simulate the magnitude and width of the hotspot effect. The method corrects the volumetric scattering component of the RTLSR model using an exponential approximation of a physical hotspot kernel, which recreates the hotspot magnitude and width using two free parameters (C(sub 1) and C(sub 2), respectively). The approach allows one to reconstruct, with reasonable accuracy, the hotspot effect by adjusting or using the prior values of these two hotspot variables. Our results demonstrate that: (1) significant improvements in capturing hotspot effect can be made to this method by using the inverted hotspot parameters; (2) the reciprocal nature allow this method to be more adaptive for simulating the hotspot height and width with high accuracy, especially in cases where hotspot signatures are available; and (3) while the new approach is consistent with the heritage RTLSR model inversion used to estimate intrinsic narrowband and broadband albedos, it presents some differences for vegetation clumping index (CI) retrievals. With the hotspot-related model parameters determined a priori, this method offers improved performance for various ecological remote sensing applications; including the estimation of canopy structure parameters.

airborne measurements↗

Solar Radiation Pressure Binning for the Geosynchronous Orbit

Orbital maintenance parameters for individual satellites or groups of satellites have traditionally been set by examining orbital parameters alone, such as through apogee and perigee height binning; this approach ignored the other factors that governed an individual satellite's susceptibility to non-conservative forces. In the atmospheric drag regime, this problem has been addressed by the introduction of the "energy dissipation rate," a quantity that represents the amount of energy being removed from the orbit; such an approach is able to consider both atmospheric density and satellite frontal area characteristics and thus serve as a mechanism for binning satellites of similar behavior. The geo-synchronous orbit (of broader definition than the geostationary orbit -- here taken to be from 1300 to 1800 minutes in orbital period) is not affected by drag; rather, its principal non-conservative force is that of solar radiation pressure -- the momentum imparted to the satellite by solar radiometric energy. While this perturbation is solved for as part of the orbit determination update, no binning or division scheme, analogous to the drag regime, has been developed for the geo-synchronous orbit. The present analysis has begun such an effort by examining the behavior of geosynchronous rocket bodies and non-stabilized payloads as a function of solar radiation pressure susceptibility. A preliminary examination of binning techniques used in the drag regime gives initial guidance regarding the criteria for useful bin divisions. Applying these criteria to the object type, solar radiation pressure, and resultant state vector accuracy for the analyzed dataset, a single division of "large" satellites into two bins for the purposes of setting related sensor tasking and orbit determination (OD) controls is suggested. When an accompanying analysis of high area-to-mass objects is complete, a full set of binning recommendations for the geosynchronous orbit will be available.

Hejduk, M. D.↗