Search NASA⌕ Search

SEARCH · Search NASA

Results for “local approximate GPs”

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.

A localized ensemble of approximate Gaussian processes for fast sequential emulation

More attention has been given to the computational cost associated with the fitting of an emulator. Substantially less attention is given to the computational cost of using that emulator for prediction. This is primarily because the cost of fitting an emulator is usually far greater than that of obtaining a single prediction, and predictions can often be obtained in parallel. In many settings, especially those requiring Markov Chain Monte Carlo, predictions may arrive sequentially and parallelization is not possible. In this case, using an emulator procedure which can produce accurate predictions efficiently can lead to substantial time savings in practice. In this paper, we propose a global model approximate Gaussian process framework via extension of a popular local approximate Gaussian process (laGP) framework. Our proposed emulator can be viewed as a treed Gaussian process where the leaf nodes are laGP models, and the tree structure is learned greedily as a function of the prediction stream. The suggested method (called leapGP) has interpretable tuning parameters which control the time‐memory trade‐off. One reasonable choice of settings leads to an emulator with a training cost and makes predictions rapidly with an asymptotic amortized cost of .

97 MATHEMATICS AND COMPUTING↗

Airborne gravity and other geophysical techniques for understanding the lithosphere beneath the West Antarctic Ice Sheet

As part of a program entitled Corridor Aerogeophysics of the Southeastern Ross Transect Zone (CASERTZ), an aerogeophysical platform was developed to study the interaction of geological and glaciological processes in West Antarctica. A de Havilland Twin Otter was equipped with an ice-penetrating radar, a proton precession magnetometer, an airborne gravity system, and a laser altimeter. The 60-MHz ice-penetrating radar can recover sub-ice topography with an accuracy of about 10 m through 3 km of comparatively warm West Antarctic ice, while the laser altimeter profiling of the ice surface is accurate to approximately 1 m. The magnetic field observations are accurate to several nT, and the gravity measurements are accurate to better than 3 mGal. The aircraft is navigated by a local radio transponder network, while differential positioning techniques based on the Global Positioning System (GPS) satellites are used for recovering high-resolution horizontal and vertical positions. Attitude information from an inertial navigation system is used to correct the laser altimetry and a digital pressure transducer is used to recover vertical positions and accelerations in the absence of satellite positioning. Continuous base-station observations are made for the differential GPS positioning and the removal of ionospheric noise from the airborne magnetometer measurements.

Bell, Robin E.↗

Total ionospheric electron content calibration using SERIES GPS satellite data

The current status of the Deep Space Network advanced systems research into ionospheric calibration techniques, based on Global Positioning System (GPS) data is described. A GPS-based calibration system is planned to replace the currently used Faraday rotation method by 1989. The SERIES receiver system determines the differential group delay of signals transmitted at two different carrier frequencies. This differential delay includes an ionospheric component and a GPS transmitter offset. The transmitter offsets are different for each GPS satellite. Tests were conducted to assess the effect of the offsets on the ionospheric calibration accuracy. From the obtained data, the total electron content and GPS transmitter offsets are calculated by a least squares estimation method employing a local model of total ionospheric electron content. The end product is an estimation of the total ionospheric content for an arbitrary line-of-sight direction. For the presented polynomial fitting technique, the systematic error due to mismodeling is estimated to be approximately 6 x 10 to the 16th power el/sq m, while the formal error is approximately 2 x 10 to the 16th power el/sq m. The final goal is an error of 3 x 10 to the 16th power el/sq m (approximately 0.7 ns at 2.3 GHz).

Lanyi, G.↗

Ground and Airborne Methane Measurements using Optical Parametric Amplifiers

We report on an initial airborne demonstration of atmospheric methane column measurements at 1.65 micrometers using a widely tunable, seeded optical parametric amplifier (OPA) lidar and a photon counting detector. Methane is an important greenhouse gas and accurate knowledge of its sources and sinks is needed for climate modeling. Our lidar system uses 20 pulses at increasing wavelengths and integrated path differential absorption (IPDA) to map a methane line at 1650.9 nanometers. The wavelengths are generated by using a Nd:YAG pump laser at 1064.5 nanometers and distributed feedback diode laser at 1650.9 nanometers and a periodically-poled lithium niobate (PPLN) crystal. The pulse width was 3 nanoseconds and the pulse repetition rate was 6.28 KHz. The outgoing energy was approximately 13 microJoules/pulse. A commercial 20 nanometer diameter fiber-coupled telescope with a photon counting detector operated in analog mode with a 0.8 nanometer bandpass filter was used as the lidar receiver. The lidar system was integrated on NASA's DC-8 flying laboratory, based at Dryden Airborne operations Facility (DAOF) in Palmdale CA. Three flights were performed in the central valley of California. Each flight lasted about 2.5 hours and it consisted of several flight segments at constant altitudes at approximately 3, 4.5, 6, 7.6, 9.1, 10.6 km (l0, 15, 20, 25, 30, 35 kft). An in-situ cavity ring down spectrometer made by Picarro Inc. was flown along with the lidar instrument provided us with the "truth" i.e. the local CH4, CO2 and H2O concentrations at the constant flight altitude segments. Using the aircraft's altitude, GPS, and meteorological data we calculated the theoretical differential optical depth of the methane absorption at increasing altitudes. Our results showed good agreement between the experimentally derived optical depth measurements from the lidar instrument and theoretical calculations as the flight altitude was increased from 3 to 10.6 kilometers, assuming a constant methane mixing ratio of 1.8 parts per million. The in-situ spectrometer did not show any significant deviations from the ambient concentrations. Further analysis using meteorological data from the Global Modeling and Assimilation Office (http://gmao.gsfc.nasa.gov/) to derive the theoretical optical depth also showed good agreement with the experimentally derived values. The OPA lidar system with slight modifications has also been used to measure CO2, water vapor, and CO in the near and mid-infrared spectral regions on the ground.

Riris, Haris↗

Ionospheric calibration for single frequency altimeter measurements

This study is a preliminary analysis of the effectiveness (in terms of altimeter calibration accuracy) of various ionosphere models and the Global Positioning System (GPS) to calibrate single frequency altimeter height measurements for ionospheric path delay. In particular, the research focused on ingesting GPS Total Electron Content (TEC) data into the physical Parameterized Real-Time Ionospheric Specification Model (PRISM), which estimates the composition of the ionosphere using independent empirical and physical models and has the capability of adjusting to additional ionospheric measurements. Two types of GPS data were used to adjust the PRISM model: GPS receiver station data mapped from line-of-sight observations to the vertical at the point of interest and a grid map (generated at the Jet Propulsion Laboratory) of GPS derived TEC in a sun-fixed longitude frame. The adjusted PRISM TEC values, as well as predictions by the International Reference Ionosphere (IRI-90), a climatological (monthly mean) model of the ionosphere, were compared to TOPEX dual-frequency TEC measurements (considered as truth) for a number of TOPEX sub-satellite tracks. For a 13.6 GHz altimeter, a Total Electron Content (TEC) of 1 TECU 10(exp 16) electrons/sq m corresponds to approximately 0.218 centimeters of range delay. A maximum expected TEC (at solar maximum or during solar storms) of 10(exp 18) electrons/sq m will create 22 centimeters of range delay. Compared with the TOPEX data, the PRISM predictions were generally accurate within the TECU when the sub-satellite track of interest passed within 300 to 400 km of the GPS TEC data or when the track passed through a night-time ionosphere. If neither was the case, in particular if the track passed through a local noon ionosphere, the PRISM values differed by more than 10 TECU and by as much as 40 TECU. The IRI-90 model, with no current ability to unseat GPS data, predicted TEC to a slightly higher error of 12 TECU. The performance of PRISM is very promising for predicting TEC and will prove useful for calibrating single frequency altimeter height measurements for ionospheric path delay. When adjusted to the GPS line-of-sight data the PRISM URSI empirical model predicted TEC over a day's period to within a global error of 8.60 TECU rms during a nighttime ionosphere and 9.74 TECU rms during the day. When adjusted to the GPS derived TEC grid, the PRISM parametrized model predicted TEC to within an error of 8.47 TECU rms for a nighttime ionosphere and 12.83 TECU rms during the day. However, the grid cannot be considered globally due to the lack of sufficient numbers of GPS stations and large latitude gaps in GPS data. It is the opinion of the authors that using the PRISM model and adjusting to the global sun-fixed TEC grid regenerated with a localized weighted interpolation offers the best possibility of meeting the 10 TECU global rms (or 2 cm at 13.6 GHz) ionosphere range correction accuracy requirement of TOPEX/Poseidon and should be the subject of further study. However, it is clear that the anticipated requirement of 34 TECU global rms for TOPEX/Poseidon Follow-On (corresponding to the TOPEX/Poseidon performance) can not be met with any realizable combination of existing models and data assimilation schemes.

Schreiner, William S.↗

Atmospheric pressure loading effects on Global Positioning System coordinate determinations

Earth deformation signals caused by atmospheric pressure loading are detected in vertical position estimates at Global Positioning System (GPS) stations. Surface displacements due to changes in atmospheric pressure account for up to 24% of the total variance in the GPS height estimates. The detected loading signals are larger at higher latitudes where pressure variations are greatest; the largest effect is observed at Fairbanks, Alaska (latitude 65 deg), with a signal root mean square (RMS) of 5 mm. Out of 19 continuously operating GPS sites (with a mean of 281 daily solutions per site), 18 show a positive correlation between the GPS vertical estimates and the modeled loading displacements. Accounting for loading reduces the variance of the vertical station positions on 12 of the 19 sites investigated. Removing the modeled pressure loading from GPS determinations of baseline length for baselines longer than 6000 km reduces the variance on 73 of the 117 baselines investigated. The slight increase in variance for some of the sites and baselines is consistent with expected statistical fluctuations. The results from most stations are consistent with approximately 65% of the modeled pressure load being found in the GPS vertical position measurements. Removing an annual signal from both the measured heights and the modeled load time series leaves this value unchanged. The source of the remaining discrepancy between the modeled and observed loading signal may be the result of (1) anisotropic effects in the Earth's loading response, (2) errors in GPS estimates of tropospheric delay, (3) errors in the surface pressure data, or (4) annual signals in the time series of loading and station heights. In addition, we find that using site dependent coefficients, determined by fitting local pressure to the modeled radial displacements, reduces the variance of the measured station heights as well as or better than using the global convolution sum.

Vandam, Tonie M.↗

Three-Dimensional Planetary Surface Tracking Based on a Simple Ultra-Wideband Impulse-Radio Infrastructure

Several prototype ultra-wideband (UWB) impulse-radio (IR) tracking systems are currently under development at NASA Johnson Space Center (JSC). These systems are being studied for use in tracking of Lunar/Mars rovers and astronauts during early exploration missions when satellite navigation systems (such as GPS) are not available. To date, the systems that have been designed and tested are intended only for two-dimensional location and tracking, but these designs can all be extended to three-dimensional tracking with only minor modifications and increases in complexity. In this presentation, we will briefly review the design and performance of two of the current 2-D systems: one designed specifically for short-range, extremely high-precision tracking (approximately 1-2 cm resolution) and the other designed specifically for much longer range tracking with less stringent precision requirements (1-2 m resolution). We will then discuss a new multi-purpose system design based on a simple UWB-IR architecture that can be deployed easily on a planetary surface to support arbitrary three-dimensional localization and tracking applications. We will discuss utilization of this system as an infrastructure to provide both short-range and long-range tracking and analyze the localization performance of the system in several different configurations. We will give theoretical performance bounds for some canonical system configurations and compare these performance bounds with both numerical simulations of the system as well as actual experimental system performance evaluations.

Barton, Richard J.↗

Central Valley Water Resources: Improving California Groundwater Assessments using GRACE and InSAR Datasets for Water Resource Management

California’s Central Valley is one of the most productive agricultural areas in the world, producing approximately $20 billion in crops annually. The recent California droughts of 2007-2010 and 2011-2017 resulted in increased groundwater pumping in the Central Valley to adequately irrigate farmland. Overdrafting of the Central Valley aquifer results in groundwater depletion, land subsidence, and permanent loss of groundwater storage. In 2014, depletion of groundwater led the state of California to enact the Sustainable Groundwater Management Act (SGMA), requiring critically overdrafted, high, and medium priority sub-basins to reach sustainable levels of groundwater pumping and recharge by 2042. SGMA allows local Groundwater Sustainability Agencies the authority to create Groundwater Sustainability Plans at the sub-basin level. To assist California’s Department of Water Resources, this project quantified groundwater change and land subsidence in Central Valley sub-basins with sparse or unreliable well and GPS data. This was done using NASA’s Gravity Recovery and Climate Experiment (GRACE), GRACE Follow-On (GRACE-FO), and interferograms derived from Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Advanced Land Observing Satellite 2 (ALOS-2) Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2). Time series of the GRACE and InSAR data were compared with well and GPS data in data-dense sub-basins to determine the feasibility of these datasets for groundwater storage and subsidence monitoring. We found thatGRACE and InSAR data are effective tools for determining groundwater change and land subsidence and can be used on their own to monitor sub-basins in the absence of well and GPS data.

Water Resources↗

Central Valley Water Resources: Improving California Groundwater Assessments using GRACE and InSAR Datasets for Water Resource Management

California’s Central Valley is one of the most productive agricultural areas in the world, producing approximately $20 billion in crops annually. The recent California droughts of 2007-2010 and 2012-2019 resulted in increased groundwater pumping in the Central Valley to adequately irrigate farmland. Overdrafting of the Central Valley aquifer results in groundwater depletion, land subsidence, and permanent loss of groundwater storage. In 2014, depletion of groundwater led the state of California to enact the Sustainable Groundwater Management Act (SGMA),requiring critically overdrafted, high, and medium priority sub-basins to reach sustainable levels of groundwater pumping and recharge by 2042. SGMA allows local Groundwater Sustainability Agencies (GSAs) the authority to create Groundwater Sustainability Plans (GSPs) at the sub-basin level. To assist California’s Department of Water Resources (DWR), this project quantified groundwater change and land subsidence in Central Valley sub-basins with sparse or unreliable well and Geographic Positioning Systems (GPS) data. This was done using NASA’s Gravity Recovery and Climate Experiment (GRACE), GRACE FollowOn (GRACE-FO), and interferograms derived from Sentinel-1 C-band Synthetic Aperture Radar (C-SAR) and Advanced Land Observing Satellite 2 (ALOS-2)Phased Array L-band Synthetic Aperture Radar 2 (PALSAR-2). Time series of the GRACE and InSAR data were compared with well and GPS data in data-dense sub-basins to determine the feasibility of these datasets for groundwater storage and subsidence monitoring. We found that GRACE and InSAR data are effective tools for determining groundwater change and land subsidence and can be used on their own to monitor sub-basins in the absence of well and GPS data

Water Resources↗

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.↗

The Sun Radio Interferometer Space Experiment (SunRISE) Mission

The Sun Radio Interferometer Space Experiment (SunRISE) will provide an entirely new view on particle acceleration and transport in the inner heliosphere by creating the first low radio frequency interferometer in space to localize heliospheric radio emissions. By imaging and determining the location of decametric-hectometric (DH, 0.1 MHz–23 MHz) solar radio bursts, SunRISE will provide key information on particle acceleration mechanisms associated with coronal mass ejections (CMEs) and the magnetic field topology from active regions into interplanetary space. The SunRISE Observatory will consist of six space vehicles in a passive formation, in orbits designed to keep them within approximately 10 km of each other, and flying in a supersynchronous geosynchronous Earth orbit (GEO). Each space vehicle consists of a Solar DH-GNSS payload and a 6U form factor spacecraft. The SunRISE Observatory together with significant ground-based processing, will enable imaging of the Sun in a portion of the spectrum that is blocked by the ionosphere and cannot be observed from Earth. Key aspects that enable this mission are that only position knowledge of the space vehicles is required, not active control, and that the architecture involves a modest amount of on-board processing coupled with significant ground-based processing for navigation, position determination, and science operations. Mission-enabling advances in software-defined radios, GPS navigation and timing, and small spacecraft technologies, developed and flown on the DARPA High Frequency Research (DHFR) and the Community Initiative for Continuing Earth Radio Occultation (CICERO) have made this mission affordable and low-risk. The SunRISE mission will exploit the multiple spacecraft per aperture (MSPA) capability of NASA’s Deep Space Network (DSN), for more efficient data transfers of larger data volumes, and utilize commercial access to space, in which the SunRISE space vehicles will be carried to their target orbit as secondary payloads in conjunction with a larger host spacecraft intended for GEO.

Neilson, Tim↗

Crustal Deformation and the Seismic Cycle Across the Kodiak Islands, Alaska

The Kodiak Islands are located approximately 120 to 250 km from the Alaska-Aleutian Trench - and are within the southern extent of the 1964 Prince William Sound (M(sub w) = 9.2) earthquake rupture zone. Here we report new campaign GPS results (1993-2001) from northern Kodiak Island. The rate and orientation of the horizontal velocities, relative to a fixed North America, range from 25.3 plus or minus 1.4 mm/yr at N32.9 deg. W plus or minus 2.5 to 8.5 plus or minus 1.0 mm/yr at N59.7 deg. W plus or minus 6.5 deg. In addition to the northern Kodiak data, we analyzed data from three southern Kodiak Island stations. The inland stations from both the northern and southern networks indicate a counterclockwise rotation of the velocity vectors. These results are consistent with the hypothesis that the difference between the Pacific-North American plate motion and the orientation of the down going slab would lead to 4-8 mm/yr of left-lateral slip above the unlocked, down-dip portion of the main thrust zone. The northern and southern Kodiak geodetic data are consistent with a model that includes the viscoelastic response to (1) a downgoing Pacific plate interface that is locked at shallow depths, (2) local coseismic slip in the 1964 earthquake, and (3) interseismic creep down dip from the seismogenic zone. Based on the pre-1964 and post-1944 earthquake history, as well as the pattern of interseismic earthquakes across the plate boundary zone, we hypothesize that in southern Kodiak some strain is released in moderate to large earthquakes between the occurrences of great earthquakes like the 1964 event.

Sauber, Jeanne↗

New Global Electron Density Observations from GPS-RO in the D- and E-Region Ionosphere

A novel retrieval technique is developed for electron density (N(sub e)) in the D- and E-region (80-120 km) using the high-quality 50-Hz GPS radio occultation (GPS-RO) phase measurements. The new algorithm assumes a slow, linear variation in the F-region background when the GPS-RO passes through the D- and E-region, and extracts the N(sub e) profiles at 80-130 km from the phase advance signal caused by N(sub e). Unlike the conventional Abel function, the new approach produces a sharp N(sub e) weighting function in the lower ionosphere, and the N(sub e) retrievals are in good agreement with the IRI (International Reference Ionosphere) model in terms of monthly maps, zonal means and diurnal variations. The daytime GPS-RO N(sub e) profiles can be well characterized by the alpha-Chapman function of three parameters (N(sub mE), h(sub mE) and H), showing that the bottom of E-region is deepening and sharpening towards the summer pole. At high latitudes the monthly GPS-RO N(sub e) maps at 80-120 km reveal clear enhancement in the auroral zones, more prominent at night, as a result of energetic electron precipitation (EEP) from the outer radiation belt. The D-/E-region auroral N(sub e) is strongly correlated with K(sub p) on a daily basis. The new N(sub e) data allow further comprehensive analyses of the sporadic E (E(sub s)) phenomena in connection with the background N(sub e) in the E-region. The layered (2-10 km) and fluctuated (less than 2 km) E(sub s) components, namely N(sub e_Layer) than N(sub e_Pert), are extracted with respect to the background N( sub e_Region) on a profile-by-profile basis. The N(sub e_Layer) component has a strong but highly-refined peak at approximately 105 km, with an amplitude smaller than N(sub e_Region) approximately by an order of magnitude. The N(sub e_Pert) component, which was studied extensively in the past, is approximately 2 orders of magnitude weaker than N(sub e_Layer). Both N(sub e_Layer) and N(sub e_Pert) are subject to significant diurnal and semidiurnal variations, showing downward progression with local time in amplitude. The 11-year solar cycle dominates the N(sub e) interannual variations, showing larger N(sub e_Region) and N(sub e_Layer) but smaller N(sub e_Pert) amplitudes in the solar maximum years. Enhanced Ne profiles are often observed in the polar winter, showing good correlation with solar proton events (SPEs) and geomagnetic activity. The new methodology offers great potential for retrieving low N(sub e) in the D-region, where radio propagation and communication blackouts can occur due to enhanced ionization. For space weather applications it is recommended for GPSRO operations to raise the top of high-rate data acquisition to approximately 140 km in the future.

GSP-R↗

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricityand river water management for industry and agriculture. Although the hydrological and geophysical impacts ofwater relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack ofin situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial waterstorage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satelliteobservations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The studysite is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia.Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that wasclearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period(from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200mmequivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, aforward model is developed to determine the increased water surface level corresponding to GRACE observations,estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannotcapture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the landsubsidence was calculated using the disk load model constructed based on the GRACE-derived lake level andLandsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the westerncoast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, whichcauses the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation ofthe surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidencefeature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE inconjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smallerspatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-basedalgorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysicalchanges from local to regional scales.

Tangdamrongsub, Natthachet↗

Quantifying Water Storage Change and Land Subsidence Induced by Reservoir Impoundment Using GRACE, Landsat, and GPS Data

The construction of hydropower dams is a common strategy to support a country's increasing need for electricity and river water management for industry and agriculture. Although the hydrological and geophysical impacts of water relocation are usually assessed prior to impoundment, their accuracy is generally limited due to the lack of in situ observations, especially in a remote area. This study presents a workflow to quantify the terrestrial water storage change (TWS) and land subsidence induced by a reservoir's water impoundment using multiple satellite observations (GRACE, Landsat), land surface models (CABLE, GLDAS, NCEP, ECMWF), and GPS data. The study site is the Bakun Dam, located in Sarawak, Malaysia, which is the largest hydropower dam in Southeast Asia. Commencing operation in late 2010, the dam induced a change of water mass and lake surface area that was clearly observed by GRACE and Landsat observations, respectively. During the 17-month impounding period (from August 2010 to December 2011), GRACE observed a dramatic increase of approximately 200 mm equivalent water height, while Landsat detected an increased lake extent of around 600 km2. In this paper, a forward model is developed to determine the increased water surface level corresponding to GRACE observations, estimated to be about 120 m. In contrast to GRACE, the TWS derived from land surface models cannot capture the increased TWS, due to the lack of reservoir routing algorithms in the models. In addition, the land subsidence was calculated using the disk load model constructed based on the GRACE-derived lake level and Landsat-derived lake extent; the result is validated with the GPS data from BIN1 station, located at the western coast of Borneo. The commencement stage of the Bakun Dam induces the large-scale land subsidence, which causes the GPS-BIN1 station to subside by ~9 mm, and move toward the Bakun Lake by ~4 mm. Computation of the surface displacements directly from GRACE spherical harmonic coefficient data fails to capture the subsidence feature, mainly due to the truncation error. Overall, this study demonstrates that evaluating GRACE in conjunction with Landsat, LSMs, and GPS data allows the exploitation of the gravity signal at a much smaller spatial scale than its intrinsic resolution. Benefiting from global coverage, the newly developed satellite-based algorithm is a valuable tool for assessing the impacts of reservoir operation on hydrological and geophysical changes from local to regional scales.

Natthachet Tangdamrongsub↗

A GPS-Based Pitot-Static Calibration Method Using Global Output-Error Optimization

Pressure-based airspeed and altitude measurements for aircraft typically require calibration of the installed system to account for pressure sensing errors such as those due to local flow field effects. In some cases, calibration is used to meet requirements such as those specified in Federal Aviation Regulation Part 25. Several methods are used for in-flight pitot-static calibration including tower fly-by, pacer aircraft, and trailing cone methods. In the 1990 s, the introduction of satellite-based positioning systems to the civilian market enabled new inflight calibration methods based on accurate ground speed measurements provided by Global Positioning Systems (GPS). Use of GPS for airspeed calibration has many advantages such as accuracy, ease of portability (e.g. hand-held) and the flexibility of operating in airspace without the limitations of test range boundaries or ground telemetry support. The current research was motivated by the need for a rapid and statistically accurate method for in-flight calibration of pitot-static systems for remotely piloted, dynamically-scaled research aircraft. Current calibration methods were deemed not practical for this application because of confined test range size and limited flight time available for each sortie. A method was developed that uses high data rate measurements of static and total pressure, and GPSbased ground speed measurements to compute the pressure errors over a range of airspeed. The novel application of this approach is the use of system identification methods that rapidly compute optimal pressure error models with defined confidence intervals in nearreal time. This method has been demonstrated in flight tests and has shown 2- bounds of approximately 0.2 kts with an order of magnitude reduction in test time over other methods. As part of this experiment, a unique database of wind measurements was acquired concurrently with the flight experiments, for the purpose of experimental validation of the optimization method. This paper describes the GPS-based pitot-static calibration method developed for the AirSTAR research test-bed operated as part of the Integrated Resilient Aircraft Controls (IRAC) project in the NASA Aviation Safety Program (AvSP). A description of the method will be provided and results from recent flight tests will be shown to illustrate the performance and advantages of this approach. Discussion of maneuver requirements and data reduction will be included as well as potential applications.

Foster, John V.↗

Vision Based Localization in Urban Environments

As part of DARPA's MARS2020 program, the Jet Propulsion Laboratory developed a vision-based system for localization in urban environments that requires neither GPS nor active sensors. System hardware consists of a pair of small FireWire cameras and a standard Pentium-based computer. The inputs to the software system consist of: 1) a crude grid-based map describing the positions of buildings, 2) an initial estimate of robot location and 3) the video streams produced by each camera. At each step during the traverse the system: captures new image data, finds image features hypothesized to lie on the outside of a building, computes the range to those features, determines an estimate of the robot's motion since the previous step and combines that data with the map to update a probabilistic representation of the robot's location. This probabilistic representation allows the system to simultaneously represent multiple possible locations, For our testing, we have derived the a priori map manually using non-orthorectified overhead imagery, although this process could be automated. The software system consists of two primary components. The first is the vision system which uses binocular stereo ranging together with a set of heuristics to identify features likely to be part of building exteriors and to compute an estimate of the robot's motion since the previous step. The resulting visual features and the associated range measurements are software component, a particle-filter based localization system. This system uses the map and the then fed to the second primary most recent results from the vision system to update the estimate of the robot's location. This report summarizes the design of both the hardware and software and will include the results of applying the system to the global localization of a robot over an approximately half-kilometer traverse across JPL'S Pasadena campus.

unmanned-ground vehicles↗