SEARCH · Search NASA

Results for “Localization Error”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Understanding and Utilizing PBL Height Data from Multiple Observing Systems in the GEOS System

The accuracy of PBL height simulation is a key issue in many applications including forecasting near surface meteorology and air quality, however, it is a very challenging problem due to the lack of not only comprehensive, global Planetary Boundary Layer (PBL) observations but also a strategy and infrastructure to utilize PBL height data from a variety of sensors. Following the designation of PBL as an incubation class observable in the 2017 Decadal Survey, the PBL Incubation Study Team Report [14] made clear that “a future global PBL observing system requires modeling and data assimilation as essential components.” There is an urgent need for global modeling development in order to utilize Program of Record (POR) observations, assess their impacts, and identify gaps to be filled by future PBL missions. Our overall objective is to develop PBL data assimilation capabilities in the NASA Global Earth Observing System (GEOS), focusing on PBL height from multiple observing systems, to support the assessment and use of future PBL observations. The NASA GEOS system is composed of the GEOS global atmospheric general circulation model (AGCM) and the atmospheric data assimilation system (ADAS). The PBL parameterizations include the “Lock” K-profile scheme driven by surface and cloud-top buoyancy fluxes ([4]), and the “Louis” local scheme for stable conditions based on the Richardson number ([5]). Above the mixed layer defined by the Lock surface plume, shallow cumulus convection is represented by the mass flux scheme of [9]. Additional parameterizations are summarized in [1]. The ADAS employs the hybrid 4D Ensemble- Variational (EnVar) configuration ([15]), with the ensemble providing flow-dependent background error covariance information. The resultant analysis increments are fed back to the forecast model through the 4D incremental analysis update (IAU) approach ([11]). In this study, PBL height data are being or have been generated from radiosondes, GNSS RO, satellite (CATS, CALIPSO and ICESat-2) and ground-based (MPLNET) lidars, and wind profiler. Investigations have been conducted to specify quality marks for PBL height retrievals for the data assimilation purpose. These PBL height data have different strengths and weaknesses ([2], [3], [6], [7], [8], [10]), and the satellite PBL height data provide better global coverage and complement in-situ PBL height data. Radiosondes offer high accuracy and in situ measurement of temperature and humidity profiles, but with poor spatio-temporal sampling. The in-situ observing systems like MPLNET and wind profiler provide long history of PBL height records at each station. The GNSS RO based PBL height is retrieved based on the sharp gradients in refractivity profile that represent the fine vertical structure of temperature and moisture changes above the PBL. However, not all RO refractivity profiles reach the surface depending on location and regime, and RO refractivity retrievals can be negatively biased below 2km. The PBL height data from satellite lidars provide high resolution along track PBL height retrievals, but over land they are affected by previous day convective PBL aerosol and strongly associated with mixing layer and retrievals cannot be made below thick, attenuating clouds. A successful assimilation of PBL height data requires a thorough understanding of the observing method and the retrieval algorithm for each observing system in order to use the PBL height data from multiple observing systems properly. Due to the sensitivity of PBL height data to the observing method and choice of algorithm, it is important to use a model definition appropriate for each observation type to compute differences between PBL height data and model PBL height (OmFs). The GEOS model currently includes two PBL height definitions suitable for direct comparison with observed PBL height, and additional definitions are being added in this study. Evaluation of different model PBL height definitions is underway. Meanwhile, efforts have been made in the GEOS data assimilation system to develop PBL height data assimilation capability. PBL height data can be assimilated using two different approaches. The traditional approach is to construct an observation operator and its tangent linear and adjoint, which link control variables to PBL height data from each observing system. This observation operator can be very complicated, e.g., the lidar-based PBL height observation operator includes the backscatter lidar forward observation operator, the algorithm to derive PBL height from attenuated total backscatter, interpolation, and calculations handling the mismatch between observed and model scales. The other approach is to augment PBL height to the control variable vector, and it is adopted in this study. The latter approach was also used in previous studies, e.g., the assimilation of PBL height data from radiosonde and aircraft in the Real Time Mesoscale Analysis (RTMA) system for a dispersion modelling study ([13]); the PBL height assimilation study using lidar PBL height data at Greensburg, Kansas for a field campaign ([12]). The PBL height assimilation from multiple observing systems in this study allows us to take advantage of the diverse PBL height data that provide much better global coverage collectively under different meteorological conditions and with different temporal and spatial scales. As all the PBL heights are tightly coupled with the PBL thermodynamic variables, the strong correlations, which are provided by the 4D ensemble forecast, enable PBL height data from various sources to interact and combine coherently and provide additional information for PBL temperature and moisture fields. The results of comparisons among PBL height data from different sources and the evaluation of the model PBL height definitions with the PBL height data will be presented, and the PBL height data synergy strategies and preliminary results will also be discussed at the conference.

Y. Zhu↗

NASA NTRS2006

Identifications of Four Integral Sources in the Galactic Plane via CHANDRA Localizations

Hard X-ray imaging of the Galactic plane by the INTEGRAL satellite is uncovering large numbers of 20-100 keV "IGR" sources. We present results from Chandra, INTEGRAL, optical, and IR observations of four IGR sources: three sources in the Norma region of the Galaxy(1GR J16195-4945,IGR J16207-5129, and IGR J16167-4957) and one that is closer to the Galactic center (IGR 5171 95-4100). In all four cases, one relatively bright Chandra source is seen in the INTEGRAL error circle, and these are likely to be the soft X-ray counterparts of the IGR sources. They have hard 0.3-10 keV spectra with power-law photon indices of Gamma = 0.5-1.1. While many previously studied IGR sources show high column densities (NH approx. 10(exp 23)-10(exp 24)/sq cm), only IGR J16195-4945 has a column density that could be as high as 10(exp 23)/sq cm. Using optical and IR sky survey catalogs and our own photometry, we have obtained identifications for all four sources. The J-band magnitudes are in the range 14.9-10.4, and we have used the optical/IR spectral energy distributions (SEDs) to constrain the nature of the sources. Blackbody components with temperature lower limits of >9400 K for IGR J16195-4945 and >18,000 K for IGR J16207-5129 indicate that these are very likely high-mass X-ray binaries (HMXBs). However, for IGR 516167-4957 and IGR J17195-4100, low extinction and the SEDs indicate later spectral types for the putative companions, suggesting that these are not HMXBs.

Tomsick, John A.↗

NASA NTRS1975

Atmospheric response to variations in sea-surface temperature

A two-week prediction experiment was performed with the GISS atmospheric model on a global data set beginning 20 December 1972 to test the sensitivity of the model to sea-surface temperature (SST) variations. Use of observed SST's in place of climatological monthly mean sea temperatures for surface flux calculations in the model was found to have a marked local effect on predicted precipitation over the ocean, with enhanced convection over warm SST anomalies. However, use of observed SST's did not lead to any detectable general improvement in forecast skill. The influence of the SST anomalies on daily predicted fields of pressure and geopotential was small up to about one week compared with the growth of prediction error, and no greater over a two-week period than that resulting from random errors in the initial meteorological state. The 14-day average fields of sea-level pressure and 500 mb height predicted by the model were similarly insensitive to the SST anomalies.

Spar, J.↗

NASA NTRS1978

Lunar occultation of Saturn. IV - Astrometric results from observations of the satellites

The method of determining local lunar limb slopes, and the consequent time scale needed for diameter studies, from accurate occultation timings at two nearby telescopes is described. Results for photoelectric observations made at Mauna Kea Observatory during the occultation of Saturn's satellites on March 30, 1974, are discussed. Analysis of all observations of occultations of Saturn's satellites during 1974 indicates possible errors in the ephemerides of Saturn and its satellites.

Dunham, D. W.↗

NASA NTRS1984

Mean excitation energies for stopping powers in various materials using local plasma oscillator strengths

The basic model of Lindhard and Scharff, known as the local plasma model, is used to study the effects on stopping power of the chemical and physical state of the medium. Unlike previous work with the local plasma model, in which individual electron shifts in the plasma frequency were estimated empirically, he Pines correction derived for a degenerate Fermi gas is shown herein to provide a reasonable estimate, even on the atomic scale. Thus, the model is moved to a complete theoretical base requiring no empirical adjustments, as characteristic of past applications. The principal remaining error is in the overestimation of the low-energy absorption properties that are characteristic of the plasma model in the region of the atomic discrete spectrum, although higher-energy phenomena are accurately represented, and even excitation-to-ionization ratios are given to fair accuracy. Mean excitation energies for covalent-bonded gases and solids, for ionic gases and crystals, and for metals are calculated using first-order models of the bonded states.

Wilson, J. W.↗

NASA NTRS1985

A Preliminary Assessment of LANDSAT-4 Thematic Mapper Data

The geometric errors for image to map rectification of one Thematic Mapper (TM) scene of an area near Windsor, Ontario were studied. The scene had been produced on computer compatible tape by NASA and contained radiometric and system corrections for geometric distortions. The geometric properties of TM photographic imagery permitted very good fitting to map detail using simple scaling techniques to localized areas and, using simple scaling, the overall geometry remained with 500 meters or 0.4 mm at the image scale of 1:1,141,600. An affine transformation, permitting differential scaling, slightly improves the fit to about 400 meters or 0.35 mm at image scale. The imagery shows promise of having the needed additional resolution and spectral discrimination to provide map revision information in urban-rural areas where the MSS sensor is considered inadequate. The late-season prairie image, however, did not hold such promise, and judgement must be reserved until images are acquired at other seasons in this particular geographic area.

Goodenough, D. G.↗

NASA NTRS1991

On error sources during airborne measurements of the ambient electric field

The principal sources of errors during airborne measurements of the ambient electric field and charge are addressed. Results of their analysis are presented for critical survey. It is demonstrated that the volume electric charge has to be accounted for during such measurements, that charge being generated at the airframe and wing surface by droplets of clouds and precipitation colliding with the aircraft. The local effect of that space charge depends on the flight regime (air speed, altitude, particle size, and cloud elevation). Such a dependence is displayed in the relation between the collector conductivity of the aircraft discharging circuit - on one hand, and the sum of all the residual conductivities contributing to aircraft discharge - on the other. Arguments are given in favor of variability in the aircraft electric capacitance. Techniques are suggested for measuring from factors to describe the aircraft charge.

Evteev, B. F.↗

NASA NTRS1990

A compressed marine data set for geomagnetic field modeling

Some 13 million scalar magnetic field data points that have been collected from the world's ocean areas reside in the collection of the National Geophysical Data Center. In order to derive a suitable data set for modeling the geomagnetic field of the earth, each ship track is divided into 220 km segments. The distribution of the reduced data in position, time and local time is discussed. The along-track filtering process described has proved to be an effective method of condensing large numbers of shipborne magnetic data into a manageable and meaningful data set for field modeling. This process also provides the benefits of smoothing short-wavelength crystal anomalies, discarding data recorded during magnetically noisy periods, and assigning reasonable error estimates to be utilized in the least squares modeling.

Langel, R. A.↗

NASA NTRS1990

Robust, angstrom level circularity profilometry

A noncontacting approach is presented which involves measuring the local circumferential curvature of the test piece by simultaneously measuring its circumferential slope at two slightly displaced locations. A pair of sensing beams is scanned along the circumference, and a profile of curvature is built, from which the circularity profile is deduced. The sensing of curvature makes the approach insensitive to all types of vibration and drift and runout errors in the relative rotation. The special qualities of the approach are summarized which make it well suited to measuring cylindrical optics and enable it to accommodate radii as small as twenty millimeters.

Glenn, Paul↗

NASA NTRS2011

Subsurface Emission Effects in AMSR-E Measurements: Implications for Land Surface Microwave Emissivity Retrieval

An analysis of land surface microwave emission time series shows that the characteristic diurnal signature associated with subsurface emission in sandy deserts carry over to arid and semi-arid region worldwide. Prior work found that diurnal variation of Special Sensor Microwave/Imager (SSM/I) brightness temperatures in deserts was small relative to International Satellite Cloud Climatology Project land surface temperature (LST) variation and that the difference varied with surface type and was largest in sand sea regions. Here we find more widespread subsurface emission effects in Advanced Microwave Scanning Radiometer-EOS (AMSR-E) measurements. The AMSR-E orbit has equator crossing times near 01:30 and 13 :30 local time, resulting in sampling when near-surface temperature gradients are likely to be large and amplifying the influence of emission depth on effective emitting temperature relative to other factors. AMSR-E measurements are also temporally coincident with Moderate Resolution Imaging Spectroradiometer (MODIS) LST measurements, eliminating time lag as a source of LST uncertainty and reducing LST errors due to undetected clouds. This paper presents monthly global emissivity and emission depth index retrievals for 2003 at 11, 19, 37, and 89 GHz from AMSR-E, MODIS, and SSM/I time series data. Retrieval model fit error, stability, self-consistency, and land surface modeling results provide evidence for the validity of the subsurface emission hypothesis and the retrieval approach. An analysis of emission depth index, emissivity, precipitation, and vegetation index seasonal trends in northern and southern Africa suggests that changes in the emission depth index may be tied to changes in land surface moisture and vegetation conditions

Galantowicz, John F.↗

NASA NTRS2012

Axisymmetric Flow Properties for Magnetic Elements of Differing Strength

Aspects of the structure and dynamics of the flows in the Sun's surface shear layer remain uncertain and yet are critically important for understanding the observed magnetic behavior. In our previous studies of the axisymmetric transport of magnetic elements we found systematic changes in both the differential rotation and the meridional flow over the course of Solar Cycle 23. Here we examine how those flows depend upon the strength (and presumably anchoring depth) of the magnetic elements. Line of sight magnetograms obtained by the HMI instrument aboard SDO over the course of Carrington Rotation 2097 were mapped to heliographic coordinates and averaged over 12 minutes to remove the 5-min oscillations. Data masks were constructed based on the field strength of each mapped pixel to isolate magnetic elements of differing field strength. We used Local Correlation Tracking of the unmasked data (separated in time by 1- to 8-hours) to determine the longitudinal and latitudinal motions of the magnetic elements. We then calculated average flow velocities as functions of latitude and longitude from the central meridian for approx 600 image pairs over the 27-day rotation. Variations with longitude indicate and characterize systematic errors in the flow measurements associated with changes in the signal from disk center to limb. Removing these systematic errors reveals changes in the axisymmetric flow properties that reflect changes in flow properties with depth in the surface shear layer.

Rightmire-Upton, Lisa↗

NASA NTRS2006

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

NASA NTRS2004

SIM spectral characteristics and accuracy analysis

The Space Interferometry Mission (SIM) will perform global astrometry (full sky), local wide-angle (15 degree) and narrow-angle (1 degree) observations to search extra-solar planets, and can calibrate stellar and galactic evolution theories. The astrometric accuracy of the SIM mission depends on spectral characteristics of the optics, detectors and targets. This paper will discuss the photometric throughput of the SIM instrument, and analyze effects of wavefront errors, optical mismatches and control biases as a function of wavelength. The color dependence models of the instrument optics including mirrors, lenses, field-stop and beam-splitter are presented. The performances of different detectors with a variety of coatings are compared. A model of the SIM fringe spectrometer is created. For early and late types of stars, brightness dependency errors are analyzed for different combinations of optics and detectors. Visibility loss due to imperfect optics is investigated in detail. Based on the models of instrument and estimated visibilities, the astrometric accuracies for various kinds of stars are evaluated. It is important to emphasize that not only light sources, mirrors, lenses, field stop and detectors are all wavelength dependent, but also fringe visibility loss, wavefront error, optics control error, etc. are all a function of wavelengths. For the first time the estimate of SIM performance is based on spectral analysis of all factors above, rather than monochromatic approximations of detected fringes, or simply adopted constants. This paper summarizes the astrometric accuracies for a wide range of stars and various combinations of optical design and detector configurations. It has been verified that SIM has astrometric accuracy of about 4 μas for targets with different spectra.

Shao, Michael↗

NASA NTRS1989

Design Of Combined Stochastic Feedforward/Feedback Control

Methodology accommodates variety of control structures and design techniques. In methodology for combined stochastic feedforward/feedback control, main objectives of feedforward and feedback control laws seen clearly. Inclusion of error-integral feedback, dynamic compensation, rate-command control structure, and like integral element of methodology. Another advantage of methodology flexibility to develop variety of techniques for design of feedback control with arbitrary structures to obtain feedback controller: includes stochastic output feedback, multiconfiguration control, decentralized control, or frequency and classical control methods. Control modes of system include capture and tracking of localizer and glideslope, crab, decrab, and flare. By use of recommended incremental implementation, control laws simulated on digital computer and connected with nonlinear digital simulation of aircraft and its systems.

Halyo, Nesim↗

NASA NTRS1989

Radiometric calibration of Landsat Thematic Mapper Thermal Band

Radiometric calibration of satellite-acquired data is essential for quantitative scientific studies, as well as for a variety of image-processing applications. This paper describes a multiyear, on-orbit radiometric calibration of the Landsat Thematic Mapper (TM) Band 6 conducted at DOE's Pacific Northwest Laboratory. Numerous Landsat TM scenes acquired and analyzed included day and night coverages at several geographical locations over several seasons. Concurrent with Landsat overpasses, thermal field and local meteorological (surface and radiosonde) measurements were collected. At-satellite (uncorrected) radiances and temperatures for water and nonwater land cover were compared to ground truth (GT) measurements after making adjustments for atmospheric (using LOWTRAN), mixed-pixel, and emissivity effects. Results indicate that, for both water and nonwater features, TM Band 6 average corrected temperature determinations using local radiosonde data to adjust for atmospheric effects, and using appropriate emissivities, are within 1.0 C of GT temperature values. Temperatures of water pixels derived from uncorrected TM Band 6 data varied roughly between 1 and 3 C of ground truth values for water temperatures ranging between 4 and 24 C. Moreover, corrections using nonlocal and noncoincident radiosonde data resulted in errors as large as 12 C. Corrections using the U.S. Standard Atmosphere gave temperature values within 1 to 2 C of GT. The average uncertainty for field instruments was + or - 0.2 C; average uncertainty for Landsat TM corrected temperature determinations was + or - 0.4 C.

Wukelic, G. E.↗

NASA NTRS1999

Global Velocities from VLBI

Precise geodetic Very Long Baseline Interferometry (VLBI) measurements have been made since 1979 at about 130 points on all major tectonic plates, including stable interiors and deformation zones. From the data set of about 2900 observing sessions and about 2.3 million observations, useful three-dimensional velocities can be derived for about 80 sites using an incremental least-squares adjustment of terrestrial, celestial, Earth rotation and site/session-specific parameters. The long history and high precision of the data yield formal errors for horizontal velocity as low as 0.1 mm/yr, but the limitation on the interpretation of individual site velocities is the tie to the terrestrial reference frame. Our studies indicate that the effect of converting precise relative VLBI velocities to individual site velocities is an error floor of about 0.4 mm/yr. Most VLBI horizontal velocities in stable plate interiors agree with the NUVEL-1A model, but there are significant departures in Africa and the Pacific. Vertical precision is worse by a factor of 2-3, and there are significant non-zero values that can be interpreted as post-glacial rebound, regional effects, and local disturbances.

Ma, Chopo↗

NASA NTRS2014

Implementing Badhwar-O'Neill Galactic Cosmic Ray Model for the Analysis of Space Radiation Exposure

For the analysis of radiation risks to astronauts and planning exploratory space missions, accurate energy spectrum of galactic cosmic radiation (GCR) is necessary. Characterization of the ionizing radiation environment is challenging because the interplanetary plasma and radiation fields are modulated by solar disturbances and the radiation doses received by astronauts in interplanetary space are likewise influenced. A model of the Badhwar‐O'Neill 2011 (BO11) GCR environment, which is represented by GCR deceleration potential theta, has been derived by utilizing all of the GCR measurements from balloons, satellites, and the newer NASA Advanced Composition Explorer (ACE). In the BO11 model, the solar modulation level is derived from the mean international sunspot numbers with time‐delay, which has been calibrated with actual flight instrument measurements to produce better GCR flux data fit during solar minima. GCR fluxes provided by the BO11 model were compared with various spacecraft measurements at 1 AU, and further comparisons were made for the tissue equivalent proportional counters measurements at low Earth orbits using the high‐charge and energy transport (HZETRN) code and various GCR models. For the comparison of the absorbed dose and dose equivalent calculations with the measurements by Radiation Assessment Detector (RAD) at Gale crater on Mars, the intensities and energies of GCR entering the heliosphere were calculated by using the BO11 model, which accounts for time‐dependent attenuation of the local interstellar spectrum of each element. The BO11 model, which has emphasized for the last 24 solar minima, showed in relatively good agreement with the RAD data for the first 200 sols, but it was resulted in to be less well during near the solar maximum of solar cycle 24 due to subtleties in the changing heliospheric conditions. By performing the error analysis of the BO11 model and the optimization in reducing overall uncertainty, the resultant BO13 model corrects the fit at solar maxima as well as being accurate at solar minima. The BO13 model is implemented to the NASA Space Cancer Risk model for the assessment of radiation risks. Overall cumulative probability distribution of solar modulation parameters represents the percentile rank of the average interplanetary GCR environment, and the probabilistic radiation risks can be assessed for various levels of GCR environment to support mission design and operational planning for future manned space exploration missions.

Kim, Myung-Hee Y.↗

NASA NTRS2019

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↗