Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Waveform Averaging”

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.

174 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Optical communication with semiconductor laser diodes

Slot timing recovery in a direct detection optical PPM communication system can be achieved by processing the photodetector waveform with a nonlinear device whose output forms the input to a phase lock group. The choice of a simple transition detector as the nonlinearity is shown to give satisfactory synchronization performance. The rms phase error of the recovered slot clock and the effect of slot timing jitter on the bit error probability were directly measured. The experimental system consisted of an AlGaAs laser diode (lambda = 834 nm) and a silicon avalanche photodiode (APD) photodetector and used Q=4 PPM signaling operated at a source data rate of 25 megabits/second. The mathematical model developed to characterize system performance is shown to be in good agreement with actual performance measurements. The use of the recovered slot clock in the receiver resulted in no degradation in receiver sensitivity compared to a system with perfect slot timing. The system achieved a bit error probability of 10 to the minus 6 power at received signal energies corresponding to an average of less than 60 detected photons per information bit.

Davidson, F.↗
DOE OSTI2022

Resource Allocation for Single Carrier Massive MIMO Systems

Resource allocation in orthogonal frequency division multiplexing (OFDM) systems is performed through allocating blocks of subcarriers to each user. Even though OFDM is the primary waveform for 5G NR systems, research reports have noted that single carrier modulation (SCM) offers several advantages over OFDM in massive multiple input multiple output (MIMO) systems, making it a preferred candidate for some future applications such as massive machine type communications (mMTC). This paper presents a method for SCM resource allocation and the relevant information recovery algorithms at the receiver. Our emphasis is on cyclic prefixed SCM, where highly flexible and efficient frequency domain detection algorithms enable the operation of many simultaneous users in a massive MIMO uplink scenario. The proposed resource allocation method allows the number of users to exceed the number of antennas at the base station (BS). Each single carrier transmission is partitioned into L interleaved streams, and each user is allocated a number of such streams. One major benefit of SCM is that each data symbol is spread over the entire bandwidth. As such, the receiver performance is dictated by the average channel gain across the transmission band rather than the channel gain at a given frequency bin or a small group of frequencies. In the proposed setup, each stream may be thought of as a resource block in SCM, analogous to resource blocks in OFDM. Hence, in the context of this paper, the terms resource blocks and streams may be used interchangeably.

5G and Beyond Communications↗
DOE OSTI2021

Dynamic Catalysis Fundamentals: I. Fast calculation of limit cycles in dynamic catalysis

Dynamic catalysis—the forced oscillation of catalytic reaction coordinate potential energy surfaces (PES)—has recently emerged as a promising method for the acceleration of heterogeneously-catalyzed reactions. Theoretical study of enhancement of rates and supra-equilibrium product yield via dynamic catalysis has, to-date, been severely limited by onerous computational demands of forward integration of stiff, coupled ordinary differential equations (ODEs) that are necessary to quantitatively describe periodic cycling between PESs. Here, we establish a new approach that reduces, by ≳108×, the computational cost of finding the time-averaged rate at dynamic steady state (i.e. the limit cycle for linear and nonlinear systems of kinetic equations). Our developments are motivated by and conceived from physical and mathematical insight drawn from examination of a simple, didactic case study for which closed-form solutions of rate enhancement are derived in explicit terms of periods of oscillation and elementary step rate constants. Generalization of such closed-form solutions to more complex catalytic systems is achieved by introducing a periodic boundary condition requiring the dynamic steady state solution to have the same periodicity as the kinetic oscillations and solving the corresponding differential equations by linear algebra or Newton-Raphson-based approaches. The methodology is well-suited to extension to non-linear systems for which we detail the potential for multiple solutions or solutions with different periodicities. For linear and non-linear systems alike, the acute decrement in computational expense enables rapid optimization of oscillation waveforms and, consequently, accelerates understanding of the key catalyst properties that enable maximization of reaction rates, conversions, and selectivities during dynamic catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
NASA NTRS2019

The Distribution, Composition, and Particle Properties of Mars Mesospheric Aerosols: an Analysis of CRISM Visible/near-IR Limb Spectra with Context from Near-Coincident MCS and MARCI Observations

The Compact Reconnaissance Imaging Spectral Mapper (CRISM) onboard the Mars Reconnaissance Orbiter (MRO) obtains pole-to-pole observations (i.e., full MRO orbits) of vertical profiles for visible/near-IR spectra (λ=0.4–4.0 μm), which are ideally suited to identifying the composition and particle sizes of Mars ice and dust aerosols over 50–100 km altitudes in the Mars mesosphere. Within the coverage limitations of the CRISM limb data set, a distinct compositional dichotomy is found in Mars mesospheric ice aerosols. CO2 ice clouds appear during the aphelion period of Mars orbit (Solar Longitudes, Ls∼0–160°) at low latitudes (∼20S–10N) over specific longitude regions (Meridiani, Valles Marineris) and at typical altitudes of 55–75 km. Apart from faint water ice hazes below 55 km, mesospheric H2O ice clouds are primarily restricted to the perihelion orbital range (Ls∼160 – 350°) at northern and southern mid-to-low latitudes with less apparent longitudinal dependences. Mars mesospheric CO2 clouds are presented in CRISM spectra with a surprisingly large range of particle sizes (cross section weighted radii, Reff=0.3 to 2.2 μm). The smaller particle sizes (Reff ≤1 μm) appear concentrated near the spatial (latitude and altitude) boundaries of their global occurrences. CRISM spectra of mesospheric CO2 clouds also show evidence of iridescence, indicating very narrow particle size distributions (effective variance, Veff∼0.03) and so very abrupt CO2 cloud nucleation. Furthermore, these clouds are sometimes accompanied by altitude coincident peaks in 1.27 μm O2 dayglow, which indicates very dry, cold regions of formation. Mesospheric water ice clouds generally exhibit small particle sizes (Reff=0.1–0.3 μm), although larger particle sizes (Reff=0.4–0.7 μm) appear infrequently. On average, water ice cloud particle sizes decrease with altitude over 50–80 km in the perihelion mesosphere. Water ice mass appears similar in clouds over a large range of observed cloud particle sizes, with particle number densities increasing to ∼10 cm−3 for Reff=0.2 μm. Near coincident Mars Climate Sounder (MCS) temperature and aerosol profile measurements for a subset of CRISM mesospheric aerosol measurements indicate near saturation (H2O and CO2) conditions for ice clouds and distinct mesospheric temperature increases associated with mesospheric dust loading. Dayside (3 pm) mesospheric CO2 clouds with larger particle sizes (Reff ≥0.5 μm) scatter surface infrared emission in MCS limb infrared radiances, as well as solar irradiance in the MCS solar band channel. Scattering of surface infrared emission is most strikingly presented in nighttime (3 am) MCS observations at 55–60 km altitudes, indicating extensive mesospheric nighttime CO2 clouds with considerably larger particle sizes (Reff∼7 μm). Mesospheric CO2 ice clouds present cirrus-like waveforms over extensive latitude and longitude regions (10°×10°), as revealed in coincident Mars Color Imager (MARCI) nadir imaging. Solar tides, gravity waves, and the large orbital variation of the extended thermal structure of the Mars atmosphere influence all of these behaviors. Mesospheric dust aerosols appear infrequently over the non-global (planet encircling) dust storm era of the CRISM limb data set (2009–2016), and exhibit smaller particle sizes (Reff=0.2–0.7 μm) relative to dust in the lower atmosphere. One isolated case of an aphelion (Ls=96°) mesospheric dust layer with large dust particle sizes (Reff ∼2 μm) over Syria Planum may reflect high altitude, non-local transport of dust over elevated regions.

Clancy, Todd R.↗
NASA NTRS2019

The distribution, composition, and particle properties of Mars mesospheric aerosols: An analysis of CRISM visible/near-IR limb spectra with context from near-coincident MCS and MARCI observations

The Compact Reconnaissance Imaging Spectral Mapper (CRISM) onboard the Mars Reconnaissance Orbiter (MRO) obtains pole-to-pole observations (i.e., full MRO orbits) of vertical profiles for visible/near-IR spectra (λ= 0.4–4.0 μm), which are ideally suited to identifying the composition and particle sizes of Mars ice and dust aerosols over 50–100 km altitudes in the Mars mesosphere. Within the coverage limitations of the CRISM limb data set, a distinct compositional dichotomy is found in Mars mesospheric ice aerosols. CO2 ice clouds appear during the aphelion period of Mars orbit (Solar Longitudes, L(s )∼ 0–160°) at low latitudes (∼20S–10N) over specific longitude regions (Meridiani, Valles Marineris) and at typical altitudes of 55–75 km. Apart from faint water ice hazes below 55 km, mesospheric H2O ice clouds are primarily restricted to the perihelion orbital range (L(s)∼160 – 350°) at northern and southern mid-to-low latitudes with less apparent longitudinal dependences. Mars mesospheric CO2 clouds are presented in CRISM spectra with a surprisingly large range of particle sizes (cross section weighted radii, R(eff) = 0.3 to 2.2 μm). The smaller particle sizes (R(eff) ≤1 μm) appear concentrated near the spatial (latitude and altitude) boundaries of their global occurrences. CRISM spectra of mesospheric CO2 clouds also show evidence of iridescence, indicating very narrow particle size distributions (effective variance, V(eff) ∼ 0.03) and so very abrupt CO2 cloud nucleation. Furthermore, these clouds are sometimes accompanied by altitude coincident peaks in 1.27 μm O2 dayglow, which indicates very dry, cold regions of formation. Mesospheric water ice clouds generally exhibit small particle sizes (R(eff) = 0.1–0.3 μm), although larger particle sizes (R(eff) = 0.4–0.7 μm) appear infrequently. On average, water ice cloud particle sizes decrease with altitude over 50–80 km in the perihelion mesosphere. Water ice mass appears similar in clouds over a large range of observed cloud particle sizes, with particle number densities increasing to ∼10 per cu. cm for R(eff) = 0.2 μm. Near coincident Mars Climate Sounder (MCS) temperature and aerosol profile measurements for a subset of CRISM mesospheric aerosol measurements indicate near saturation (H2O and CO2) conditions for ice clouds and distinct mesospheric temperature increases associated with mesospheric dust loading. Dayside (3 pm) mesospheric CO2 clouds with larger particle sizes (R(eff) ≥0.5 μm) scatter surface infrared emission in MCS limb infrared radiances, as well as solar irradiance in the MCS solar band channel. Scattering of surface infrared emission is most strikingly presented in nighttime (3 am) MCS observations at 55–60 km altitudes, indicating extensive mesospheric nighttime CO2 clouds with considerably larger particle sizes (R(eff)∼7 μm). Mesospheric CO2 ice clouds present cirrus-like waveforms over extensive latitude and longitude regions (10°×10°), as revealed in coincident Mars Color Imager (MARCI) nadir imaging. Solar tides, gravity waves, and the large orbital variation of the extended thermal structure of the Mars atmosphere influence all of these behaviors. Mesospheric dust aerosols appear infrequently over the non-global (planet encircling) dust storm era of the CRISM limb data set (2009–2016), and exhibit smaller particle sizes (R(eff) = 0.2–0.7 μm) relative to dust in the lower atmosphere. One isolated case of an aphelion (L(s) = 96°) mesospheric dust layer with large dust particle sizes (R(eff) ∼2 μm) over Syria Planum may reflect high altitude, non-local transport of dust over elevated regions.

R. Todd Clancy↗
NASA NTRS2001

Wave Normal and Poynting Vector Calculations using the Cassini Radio and Plasma Wave Instrument

Wave normal and Poynting vector measurements from the Cassini radio and plasma wave instrument (RPWS) are used to examine the propagation characteristics of various plasma waves during the Earth flyby on August 18, 1999. Using the five-channel waveform receiver (WFR), the wave normal vector is determined using the Means method for a lightning-induced whistler, equatorial chorus, and a series of low-frequency emissions observed while Cassini was in the magnetosheath. The Poynting vector for these emissions is also calculated from the five components measured by the WFR. The propagation characteristics of the lightning-induced whistler were found to be consistent with the whistler wave mode of propagation, with propagation antiparallel to the magnetic field (southward) at Cassini. The sferic associated with this whistler was observed by both Cassini and the Stanford VLF group at the Palmer Station in Antarctica. Analysis of the arrival direction of the sferic at the Palmer Station suggests that the lightning stroke is in the same sector as Cassini. Chorus was observed very close (within a few degrees) to the magnetic equator during the flyby. The chorus was found to propagate primarily away from the magnetic equator and was observed to change direction as Cassini crossed the magnetic equator. This suggests that the source region of the chorus is very near the magnetic equator. The low-frequency emission in the magnetosheath has many of the characteristics of lion roars. The average value of the angle between the wave normal vector and the local magnetic field was found to be 16 degrees, and the emissions ranged in frequency from 0. 19 to 0.75 f(sub ce), where f(sub ce) is the electron cyclotron frequency. The wave normal vectors of these waves were primarily in one direction for each individual burst (either parallel or antiparallel to the local field) but varied in direction throughout the magnetosheath. This suggests that the sources of the emissions are far from the spacecraft and that there are multiple source regions.

Hospodarsky, G. B.↗
NASA NTRS2012

Radiative Transfer Through Clouds and Its Applications in Support of the Ice, Cloud, and Land Elevation Satellite-2 (ICESat-2) Mission

The Greenland and Antarctic ice sheets, which contain enough ice to raise sea level by about 7 and 60 m, respectively, are losing mass at an increasing rate. To acquire continuous information of the cryosphere, after the Ice, Cloud, and land Elevation Satellite (ICESat) (2003-2010), NASA is actively planning for the ICESat-2 mission. Both ICESat and ICESat-2 are space-borne lidar altimetry systems. The systems measure the time of flight of the arriving photons that are reflected by the surface to deduce the elevation of the underlying terrain. As one of NASA's top priority missions, ICESat-2 is scheduled to launch in 2016. One of the major science goals of ICESat-2 is to quantify the ice sheet mass balance to determine its contributions to the sea level change and its impacts on ocean circulation (Abdalati et al. 2010). Compared to ICESat, which operates at 40 Hz and records the reflected laser energy as a waveform, the significantly improved ICESat-2 lidar employs a 532 nm micro-pulse photon counting system that operates at a high frequency of 10kHz with single photon detectability (Yang et al. 2012). To achieve its science goals, ICESat-2 requires the ability of detecting the elevation change with an accuracy of 0.2 cm/year over the entire ice sheet. Since every photon emitted by the lidar system will travel through the atmosphere, clouds can certainly affect the flight time of the arriving photons. Forward scattering by cloud particles increases the photon path length, thus resulting in biases in ice sheet elevation measurements known as atmospheric path delay (Duta et al. 2001, Yang et al. 2010, 2011). To ensure the accuracy of ICESat-2 surface elevation measurements, it is critical to understand how clouds would affect the travel time of arriving photons. In this talk, we will first present a framework that simulates the behavior of a space-borne 532 mn micro-pulse photon counting lidar in cloudy and clear atmospheres. To investigate the process of laser propagation through clouds, a 3-D Monte Carlo radiative transfer model is used to simulate the photon path distribution and the Poisson distribution is adopted for the number of photon returns. Since the photon counting system only registers the time of the first arriving photon within the detector "dead time", the retrieved average surface elevation tends to bias towards higher values. This is known as the first photon bias. With the scenarios simulated here, the first photon bias for clear sky is about 6.5 cm. Clouds affect surface altimetry in two ways: (1) cloud attenuation lowers the average number of arriving photons and hence reduces the first photon bias; (2) cloud forward scattering increases the photon path length and makes the surface appear further away from the satellite. Compared to clear sky, the average surface elevation detected by the photon counting system for cloudy sky with optical depth 1.0 is 4.0 to 6.0 cm lower for the simulations conducted. The effect of surface roughness on the accuracy of elevation retrievals will also discussed.

Yang, Yuekui↗
NASA NTRS1999

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

Accurate measurements of surface heights and atmospheric backscatter have been demonstrated with the SLA, MOLA and LITE space lidar. Recent MOLA measurements of the Mars surface have 40 cm resolution and have reduced the global uncertainty in Mars topography from a few km to approx. 10 m. GLAS is a next generation lidar being developed as part of NASA's Icesat Mission for Earth orbit . The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, determine the height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS will fly on a small dedicated spacecraft in a polar orbit at 598 km altitude with an inclination of 94 degrees. GLAS is scheduled to launch in summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will measure the vertical distance to the ice sheet from orbit with 1064 nm pulses from a Nd:Yag laser at 40 Hz. Each 5 nsec wide laser pulse is used for a single range measurement. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a I m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser, Errors in knowledge of the laser beam pointing angle can bias height measurements of sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam centroid to about 8 urad is required to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith whose measurements are combined with a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser pulse is measured with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the star camera and gyroscope permits the angular offsets of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring atmospheric backscatter profiles at both 1064 and 532 nm. The 1064 nm measurements use an analog detector and profile the height and vertical structure of thicker clouds. Measurements at 532 nm use new highly sensitive photon counting detectors, and measure the height distributions of very thin clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗
DOE OSTI2020

Detection and Perception of Sound by Eagles and Surrogate Raptors

One overarching objective of this program of study was the accumulation of objective, scientifically valid information relating to auditory performance of bald and golden eagles that may be used to guide the development of acoustic alerting/deterrence technologies intended to discourage encroachment into wind energy air spaces. To that end, analyses aimed at the characterization of sensitivity to sound in bald and golden eagles, along with findings in the supra-threshold, dynamic frequency spaces related to response latencies and amplitudes, leads us to conclude that bald, and golden eagles navigate the same basic working auditory space, as in other known and thus far characterized members of the diurnal raptor family. Specifically, bald and golden eagles, along with other raptor species within the group, operate in an auditory space characterized by a frequency band at least four octaves wide and centered on 2 kHz, with an upper frequency limit between 6 and 10 kHz at 80 dB SPL and a lower frequency limit that almost certainly extends below 0.2 kHz. Consequently, we recommend that signal designers use these data as a guideline in efforts to design effective and efficient acoustic alerting/deterrent systems. It is important to note that signal energy broadcast outside of this frequency band at moderate levels will not contribute to the efficacy of a deterrent but will add an unnecessary fraction to the overall acoustic pollution budget. The importance of this consideration is heightened by contemporaneous concerns related to the transmission of noise broadcast by wind energy farms. In addition, based on analyses of data acquired from red-tailed hawks using the same experimental paradigm and data acquisition system, we conclude that auditory function in the red-tailed hawk is sufficiently like that observed in bald and golden eagles to permit its use as a surrogate species. Response waveforms, threshold-frequency curves, and input-output characteristics match those of eagles closely. It should be noted however, that differences in sensitivity and slightly extended high-frequency limits of hearing should be taken into account when extrapolating findings from one species to the others. Although the inclusion of behavioral tests of red-tailed hawks to acoustic stimuli was beyond the scope of this investigation, future efforts to assess response parameters like signal-type preference and habituation rate will further elucidate their suitability to serve as eagle surrogates in behavioral studies; nonetheless, the species in question are well matched with respect to basic auditory performance. A second essential objective of this program of study was the acoustic characterization of a subset of calls comprising the vocal repertoires of bald and golden eagles that may be used to supplement auditory performance findings in the effort to guide the development of acoustic alerting signals. With regard to that objective, the vocal repertoires of both bald and golden eagle species are rich and varied. While similar in spectrographic structure, distinctive differences are also clear. Generally, golden eagles produce some calls with shorter durations, and similar “sounding” calls exhibit distinctively different spectrographic patterns than those of bald eagles. Both species produce calls that contain a wide variety of nonlinear elements that operate to enhance the rich and varied nature of commonly observed vocal products. Comparison of the average power spectra of commonly observed bald and golden eagle calls with threshold-frequency curves leads to the conclusion that call energies fall within the frequency bounds of hearing. Further, the acoustic energy of calls considered in this report tend to fall into overlapping, but different frequency ranges of the acoustic sensitivity curve. This condition may encourage signal designers to vary the frequency content of acoustic deterrence signals in the field. Finally, preliminary observations relating to the tendencies and proclivities of bald eagles to attend to the acoustic landscape lead to the conclusion that eagles monitor their immediate sound environment assiduously. Individuals respond to a variety of natural and synthetic sound signals reliably and, perhaps most relevant in the context of the engineering of acoustic alerting/deterrence technologies, habituation to most sounds considered in this effort was minimal. These preliminary results, while calling for extended behavioral testing, are promising and set the stage for the exportation of behavioral studies into real world scenarios.

17 WIND ENERGY↗
NASA NTRS2000

The Geoscience Laser Altimeter System (GLAS) for the ICESAT Mission

The Laser In space Technology Experiment, Shuttle Laser Altimeter and the Mars Observer Laser Altimeter have demonstrated accurate measurements of atmospheric backscatter and Surface heights from space. The recent MOLA measurements of the Mars surface have 40 cm vertical resolution and have reduced the global uncertainty in Mars topography from a few km to about 5 m. The Geoscience Laser Altimeter System (GLAS) is a next generation lidar for Earth orbit being developed as part of NASA's Icesat Mission. The GLAS design combines a 10 cm precision surface lidar with a sensitive dual wavelength cloud and aerosol lidar. GLAS will precisely measure the heights of the Earth's polar ice sheets, establish a grid of accurate height profiles of the Earth's land topography, and profile the vertical backscatter of clouds and aerosols on a global scale. GLAS is being developed to fly on a small dedicated spacecraft in a polar orbit with a 590 630 km altitude at inclination of 94 degrees. GLAS is scheduled to launch in the summer 2001 and to operate continuously for a minimum of 3 years with a goal of 5 years. The primary mission for GLAS is to measure the seasonal and annual changes in the heights of the Greenland and Antarctic ice sheets. GLAS will continuously measure the vertical distance from orbit to the Earth's surface with 1064 nm pulses from a ND:YAG laser at a 40 Hz rate. Each 5 nsec wide laser pulse is used to produce a single range measurement, and the laser spots have 66 m diameter and about 170 m center-center spacings. When over land GLAS will profile the heights of the topography and vegetation. The GLAS receiver uses a 1 m diameter telescope and a Si APD detector. The detector signal is sampled by an all digital receiver which records each surface echo waveform with I nsec resolution and a stored echo record lengths of either 200, 400, or 600 samples. Analysis of the echo waveforms within the instrument permits discrimination between cloud and surface echoes. Ground based echo analysis permits precise ranging, determining the roughness or slopes of the surface as well as the vertical distributions of vegetation illuminated by the laser. Accurate knowledge of the laser beam's pointing angle is needed to prevent height biases when over sloped surfaces. For surfaces with 2 deg. slopes, knowledge of pointing angle of the beam's centroid to about 8 urad is needed to achieve 10 cm height accuracy. GLAS uses a stellar reference system (SRS) to determine the pointing angle of each laser firing relative to inertial space. The SRS uses a high precision star camera oriented toward local zenith and a gyroscope to determine the inertial orientation of the SRS optical bench. The far field pattern of each laser is measured pulse relative to the star camera with a laser reference system (LRS). Optically measuring each laser far field pattern relative to the orientation of the star camera and gyroscope permits the precise pointing angle of each laser pulse to be determined. GLAS will also determine the vertical distributions of clouds and aerosols by measuring the vertical profile of laser energy backscattered by the atmosphere at both 1064 and 532 nm. The 1064 nm measurements use the Si APD detector and profile the height and vertical structure of thicker clouds. The measurements at 532 nm use new highly sensitive photon counting, detectors, and measure the height distributions of very thin Clouds and aerosol layers. With averaging these can be used to determine the height of the planetary boundary layer. The instrument design and expected performance will be discussed.

Abshire, James B.↗
NASA NTRS1999

Elevation Change of the Southern Greenland Ice Sheet from Satellite Radar Altimeter Data

Long-term changes in the thickness of the polar ice sheets are important indicators of climate change. Understanding the contributions to the global water mass balance from the accumulation or ablation of grounded ice in Greenland and Antarctica is considered crucial for determining the source of the about 2 mm/yr sea-level rise in the last century. Though the Antarctic ice sheet is much larger than its northern counterpart, the Greenland ice sheet is more likely to undergo dramatic changes in response to a warming trend. This can be attributed to the warmer Greenland climate, as well as a potential for amplification of a global warming trend in the polar regions of the Northern Hemisphere. In collaboration with Drs. Curt Davis and Craig Kluever of the University of Missouri, we are using data from satellite radar altimeters to measure changes in the elevation of the Southern Greenland ice sheet from 1978 to the present. Difficulties with systematic altimeter measurement errors, particularly in intersatellite comparisons, beset earlier studies of the Greenland ice sheet thickness. We use altimeter data collected contemporaneously over the global ocean to establish a reference for correcting ice-sheet data. In addition, the waveform data from the ice-sheet radar returns are reprocessed to better determine the range from the satellite to the ice surface. At JPL, we are focusing our efforts principally on the reduction of orbit errors and range biases in the measurement systems on the various altimeter missions. Our approach emphasizes global characterization and reduction of the long-period orbit errors and range biases using altimeter data from NASA's Ocean Pathfinder program. Along-track sea-height residuals are sequentially filtered and backwards smoothed, and the radial orbit errors are modeled as sinusoids with a wavelength equal to one revolution of the satellite. The amplitudes of the sinusoids are treated as exponentially-correlated noise processes with a time-constant of six days. Measurement errors (e.g., altimeter range bias) are simultaneously recovered as constant parameters. The corrections derived from the global ocean analysis are then applied over the Greenland ice sheet. The orbit error and measurement bias corrections for different missions are developed in a single framework to enable robust linkage of ice-sheet measurements from 1978 to the present. In 1998, we completed our re-evaluation of the 1978 Seasat and 1985-1989 Geosat Exact Repeat Mission data. The estimates of ice thickness over Southern Greenland (south of 72N and above 2000 m) from 1978 to 1988 show large regional variations (+/-18 cm/yr), but yield an overall rate of +1.5 +/- 0.5 cm/yr (one standard error). Accounting for systematic errors, the estimate may not be significantly different from the null growth rate. The average elevation change from 1978 to 1988 is too small to assess whether the Greenland ice sheet is undergoing a long-term change.

Haines, Bruce J.↗