Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “height estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

The Role of Orograph and Parallax Corrections on High Resolution Geostationary Satellite Rainfall Estimates for Flash Flood Applications

The current generation of geosynchronous satellites exhibits considerably improved capabilities in the area of resolution, gridding accuracy, and sampling frequency as compared to their predecessors. These improvements have made it possible to accurately observe the life cycle of small scale, short-live phenomenon like rapidly developing thunderstorms, at a very high spatial and temporal resolutions. While the gain in the improved resolution is not significant for synoptic scale cloud systems, it plays a major role on the computation of precipitation values for mesoscale and stonn scale systems. Two of the important factor on the accurate precision of precipitation from satellite imagery are the position of the cloud tops as viewed by the satellite and the influence of orographic effects on the distribution of precipitation. The first problem has to do with the fact that the accurate estimation of precipitation from data collected by a satellite in geosynchronous orbit requires the knowledge of the exact position of the cloud tops with respect to the ground below. This is not a problem when a cloud is located directly below the satellite; at large viewing angles the geographic coordinates on satellite images are dependent on cloud heights and distance from the sub-satellite point. The latitude and longitude coordinates for high convective cloud tops are displaced away from the sub-satellite point and may be shifted by as much as 20 Km from the sea level coordinates. The second problem has to do with the variations in rainfall distribution with elevation. Ground observations have shown that precipitation amounts tend to increase with height and that the slope of the hill or mountain that is facing the prevailing wind normally receives greater rainfall then do the lee slopes. The purpose of the study is to show the recent developments at the Office of Research and Applications (ORA) at the National Oceanic and Atmospheric Administration (NOAA/NESDIS) in Camp Springs, MD, USA, to adjust any satellite rainfall estimation technique and account for orographic and parallax corrections. Description and examples of the procedure applied to the current NOAA/NESDIS experimental satellite rainfall estimation technique for flash flood applications will be presented at the conference.

Vicente, Gilberto A.↗
NASA NTRS2019

Simultaenous Retrieval of Surface Roughness Parameters from Combined Active-Passive SMAP Observations

Soil roughness strongly influences processes like erosion, infiltration, moisture and evaporation of soils as well as growth of agricultural plants. An approach to soil roughness based on active-passive microwave covariation is proposed in order to simultaneously retrieve the vertical RMS height (s) and horizontal correlation length (l) of soil surfaces from simultaneously measured radar and radiometer microwave signatures. The approach is based on a retrieval algorithm for active-passive covariation including the improved Integral Equation Method (I2EM). The algorithm is tested with the global active-passive microwave observations of the SMAP mission. The developed roughness retrieval algorithm shows independence of permittivity for > 10 [-] due to the covariation formalism. Results reveal that s and l can be estimated simultaneously by the proposed approach since surface patterns of non-vegetated areas become evident on global scale. In regions with sandy deserts, like the Sahara or the outback in Australia, determined and confirm rather smooth to semi-rough surface roughness patterns with small vertical RMS heights and corresponding higher horizontal correlation lengths.

correlation length↗
NASA NTRS2009

Studies of a Space Lidar Instrument for the ASCENDS Mission

We have developed a lidar technique for measuring the tropospheric C02 concentrations as a candidate for NASA's planned ASCENDS mission. Our technique uses two pulsed laser transmitters allowing simultaneous measurement of a C02 absorption line in the 1570 nm band, 02 extinction in the Oxygen A-band and surface height and backscatter. The lidar measures the energy and time of flight of the laser echoes reflected from the atmosphere and surface. The lasers are stepped in wavelength across the C02 line and an 02 line region during the measurement. The receiver uses a telescope and photon counting detectors, and measures the background light and energies of the laser echoes from the surface along with scattering from any aerosols in the path. The gas extinction and column densities for the C02 and 02 gases are estimated from the ratio of the on- and off- line signals via the DIAL technique. Pulsed laser signals and a time gated receiver are used to isolate the laser echo signals from the surface, to reject laser photons scattered in the atmosphere, and measure the surface height and scattering profile in the path. We have recently completed a second design study for the space instrument. For the study, we selected a nominal sun-synchronous orbit with an altitude of 400 km and equator crossing time of 1:30 pm, and a receiver telescope with 1.5 m diameter.

Abshire, J. B.↗
NASA NTRS1981

Vertical extensions of galactic spiral arms

The vertical structure of the large scale galactic shock along spiral arms is studied in order to understand the extended features or the wings of spiral arms observed by Kepner (1970) and recently checked by Soukup (1979). Numerical calculations of two-dimensional gas flow were performed for two models of the gas. The results show that the shock front, standing perfectly perpendicular to the galactic plane, extends impressively above the scale height of the interstellar gas. In the extreme case of isothermal gas, the shock solution exists at any distance from the plane. In a realistic model in which the kinetic temperature of the gas is stratified, increasing with height, the shock can extend up to 700 pc. For both models, although the velocity component perpendicular to the galactic plane is small, the solutions corresponding to layered one-dimensional flow contain no shocks at such heights. In order to study the geometry of the shock front and the flow pattern near it, a simplified model is adopted, in which the Coriolis force is neglected. However, estimates show that including it will enhance the shock strength. The compression at the shock front makes the gas observable at high latitudes.

Soukup, J. E.↗
NASA NTRS1985

The extent of the local hi halo

Forty-five high-latitude, OB stars have been observed in the Ly alpha and 21 cm lines of HI in an effort to map out the vertical distribution and extent of the local HI halo. The 25 stars for which a reliable HI colum density can be obtained from Ly alpha lie between 60 and 3100 pc from the plane. The principal result is that the total column density of HI at z 1 kpc is, on the average, 5 + or - 3 x 10 the 19th power/sq cm, or 15% of the total sub HI. At relatively low z the data toward some stars suggest a low effective scale height and fairly high average foreground density, while toward others the effective scale height is large and the average density is low. This can be understood as the result of irregularities in the interstellar medium. A model with half of the HI mass in clouds having radii of a few pc and a Gaussian vertical distribution with sigma sub 2 = 135 pc, and half of the mass in an exponential component with a scale height of 500 pc, gives a satisfactory fit to the data. The technique of comparing Ly alpha and 21 cm column densities is also used to discuss the problem of estimating the distance to several possibly subluminous stars.

Lockman, F. J.↗
NASA NTRS1986

The extent of the local H I halo

Forty-five high-latitude, OB stars have been observed in the Ly alpha and 21 cm lines of HI in an effort to map out the vertical distribution and extent of the local HI halo. The 25 stars for which a reliable HI column density can be obtained from Ly alpha lie between 60 and 3100 pc from the plane. The principal result is that the total column density of HI at an absolute value of z greater than 1 kpc is, on the average, 5 + or - 3 x 10 to the 19th power/sq cm, or 15 percent of the total sub HI. At relatively low z the data toward some stars suggest a low effective scale height and fairly high average foreground density, while toward others the effective scale height is large and the average density is low. This can be understood as the result of irregularities in the interstellar medium. A model with half of the HI mass in clouds having radii of a few pc and a Gaussian vertical distribution with sigma sub 2 = 135 pc, and half of the mass in an exponential component with a scale height of 500 pc, gives a satisfactory fit to the data. The technique of comparing Ly alpha and 21 cm column densities is also used to discuss the problem of estimating the distance to several possibly subluminous stars.

Lockman, F. J.↗
NASA NTRS1994

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

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

Tsaoussi, Lucia S.↗
NASA NTRS2016

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

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

airborne measurements↗
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 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↗
NASA NTRSDate not supplied

Adaptive Pressure Profile Method to Locate the Isolator Shock Train Leading Edge Given Limited Pressure Information

To maximize the performance of high-speed air-breathing engines, such as dual-mode scramjets, the streamwise location of the shock train leading edge (STLE) is ideally placed as far upstream in the isolator as possible while avoiding engine unstart. Thus, it is of interest to quantify and control the STLE location as the vehicle travels along its flight trajectory. The STLE location is typically quantified using wall static pressure measurements but there are often restrictions on the number and placement of transducers, thus reducing the accuracy and overall capability of STLE detection methods. In this work, the Adaptive Pressure Profile (APP) method is introduced to address such sparsity concerns. This method is data driven and does not heavily rely on prior information about the flow regime or engine model. Instead, the \mname method uses real-time pressure measurements from a small number of transducers to adaptively learn the isolator pressure profile. This adaptively-learned profile is fit to the pressure data at each time instance to estimate STLE location. The \mname method produces accurate estimates even when (1) the STLE location is not bounded by two or more transducers or (2) when the STLE location is between two transducers that are situated several duct heights apart. Data from two direct-connect isolator models are used to evaluate the accuracy of the \mname method and demonstrate its robustness for different back-pressure scenarios and transducer configurations.

Robin L Hunt↗
NASA NTRSDate not supplied

Adaptive Pressure Profile Method to Locate the Isolator Shock Train Leading Edge Given Limited Pressure Information

To maximize the performance of high-speed air-breathing engines, such as dual-mode scramjets, the streamwise location of the shock train leading edge (STLE) is ideally placed as far upstream in the isolator as possible while avoiding engine unstart. Thus, it is of interest to quantify and control the STLE location as the vehicle travels along its flight trajectory. The STLE location is typically quantified using wall static pressure measurements but there are often restrictions on the number and placement of transducers, thus reducing the accuracy and overall capability of STLE detection methods. In this work, the Adaptive Pressure Profile (APP) method is introduced to address such sparsity concerns. This method is data driven and does not heavily rely on prior information about the flow regime or engine model. Instead, the \mname method uses real-time pressure measurements from a small number of transducers to adaptively learn the isolator pressure profile. This adaptively-learned profile is fit to the pressure data at each time instance to estimate STLE location. The \mname method produces accurate estimates even when (1) the STLE location is not bounded by two or more transducers or (2) when the STLE location is between two transducers that are situated several duct heights apart. Data from two direct-connect isolator models are used to evaluate the accuracy of the \mname method and demonstrate its robustness for different back-pressure scenarios and transducer configurations.

Robin Hunt↗
NASA NTRS1985

Simulations of SAR wave spectra using high spectral resolution estimates from the SCR and ROWS instruments

A numerical model for predicting the synthetic aperture radar (SAR) image of a moving ocean surface is described, and results are presented for two SIR-B data sets collected off the coast of Chile. Wave height spectra measured by the NASA radar ocean wave spectrometer (ROWS) and surface contour radar (SCR) were used as inputs to this model, and results are compared with actual SIR-B image spectra from orbits 91 and 106.

Lyzenga, D.↗
NASA NTRS1989

Probabilities and statistics for backscatter estimates obtained by a scatterometer

Methods for the recovery of winds near the surface of the ocean from measurements of the normalized radar backscattering cross section must recognize and make use of the statistics (i.e., the sampling variability) of the backscatter measurements. Radar backscatter values from a scatterometer are random variables with expected values given by a model. A model relates backscatter to properties of the waves on the ocean, which are in turn generated by the winds in the atmospheric marine boundary layer. The effective wind speed and direction at a known height for a neutrally stratified atmosphere are the values to be recovered from the model. The probability density function for the backscatter values is a normal probability distribution with the notable feature that the variance is a known function of the expected value. The sources of signal variability, the effects of this variability on the wind speed estimation, and criteria for the acceptance or rejection of models are discussed. A modified maximum likelihood method for estimating wind vectors is described. Ways to make corrections for the kinds of errors found for the Seasat SASS model function are described, and applications to a new scatterometer are given.

Pierson, Willard J., Jr.↗
NASA NTRS1979

Base line estimation using single passes of laser data

The laser data of the GEOS 3 satellite passes observed by four stations at Greenbelt (Maryland), Bermuda, Grand Turk Island (Bahamas) and Patrick Air Force Base (Florida), were employed to determine precise interstation base lines and relative heights in short orbital arcs of no more than 12-min duration. No more than five arcs of data are required to define the interstation base lines to 30-cm precision. Base lines running parallel to the orbital motion can be defined to submeter precision from a single short arc of data. Combining arcs of different orbital geometry in a common adjustment of two or more stations relative to the base station helps to compensate for weak base line definition in any single arc. This technique can be used for tracking such spacecraft as Lageos, a high-altitude retroreflector-carrying satellite designed for precise laser ranging studies.

Dunn, P. J.↗
NASA NTRS1984

Thermomechanical evolution of oceanic fracture zones

To explain the topographic and gravitational signals of fracture zones (FZs), a realistic model is developed and compared with Seasat altimeter data across three major FZs. The model is based on the thermal, subsidence, and mechanical models of the oceanic lithosphere. The model is used to calculate the FZ topography and the evolution of the geoid signature across a FZ. When lithospheric flexure is included, the model explains the persistence of sharp FZ geoid steps on old seafloor as well as the asymmetry of the geoid steps across major FZs. The model does not explain the disappearance of the thermal part of the FZ geoid signature. The model profiles are used to access the various techniques for estimating the overall geoid offsets across FZs. It is shown that steps may be underestimated by as much as 50 percent. These predictions, if correct, will force reinterpretation of published geoid height versus age data derived from altimeter profiles across FZs.

Sandwell, D. T.↗
NASA NTRS2012

Analysis of Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear CO2 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed ~ 1 ppm random errors for 8-10 km altitudes and ~30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2lidar on the NASA DC-8 and added an 02lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected ~linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. We demonstrated measurements over the California Central Valley, to stratus cloud tops over the Pacific Ocean, over mountain regions with snow, and over several areas with broken clouds. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗
NASA NTRS2012

Pulsed Lidar Measurements of Atmospheric CO2 Column Absorption and Range During the ASCENDS 2009-2011 Airborne Campaigns

We have developed a pulsed lidar technique for measuring the tropospheric CO2 concentrations as a candidate for NASA's ASCENDS mission and have demonstrated the CO2 and O2 measurements from aircraft. Our technique uses two pulsed lasers allowing simultaneous measurement of a single CO2 absorption line near 1572 nm, O2 extinction in the Oxygen A-band, surface height and backscatter profile. The lasers are stepped in wavelength across the CO2 line and an O2 line doublet during the measurement. The column densities for the CO2 and O2 are estimated from the differential optical depths (DOD) of the scanned absorption lines via the IPDA technique. For the 2009 ASCENDS campaign we flew the CO2 lidar only on a Lear-25 aircraft, and measured the absorption line shapes of the CO2 line using 20 wavelength samples per scan. Measurements were made at stepped altitudes from 3 to 12.6 km over the Lamont OK, central Illinois, North Carolina, and over the Virginia Eastern Shore. Although the received signal energies were weaker than expected for ASCENDS, clear C02 line shapes were observed at all altitudes. Most flights had 5-6 altitude steps with 200-300 seconds of recorded measurements per step. We averaged every 10 seconds of measurements and used a cross-correlation approach to estimate the range to the scattering surface and the echo pulse energy at each wavelength. We then solved for the best-fit CO2 absorption line shape, and calculated the DOD of the fitted CO2 line, and computed its statistics at the various altitude steps. We compared them to CO2 optical depths calculated from spectroscopy based on HITRAN 2008 and the column number densities calculated from the airborne in-situ readings. The 2009 measurements have been analyzed in detail and they were similar on all flights. The results show clear CO2 line shape and absorption signals, which follow the expected changes with aircraft altitude from 3 to 13 km. They showed the expected nearly the linear dependence of DOD vs altitude. The measurements showed -1 ppm random errors for 8-10 km altitudes and -30 sec averaging times. For the 2010 ASCENDS campaigns we flew the CO2 lidar on the NASA DC-8 and added an O2 lidar channel. During July 2010 we made measurements of CO2 and O2 column absorption during longer flights over Railroad Valley NV, the Pacific Ocean and over Lamont OK. CO2 measurements were made with 30 steps/scan, 300 scans/sec and improved line resolution and receiver sensitivity. Analysis of the 2010 CO2 measurements shows the expected -linear change of DOD with altitude. For measurements at altitudes> 6 km the random errors were 0.3 ppm for 80 sec averaging times. For the summer 2011 ASCENDS campaigns we made further improvements to the lidar's CO2 line scan and receiver sensitivity. The seven flights in the 2011 Ascends campaign were flown over a wide variety of surface and cloud conditions in the US, which produced a wide variety of lidar signal conditions. Details of the lidar measurements and their analysis will be described in the presentation.

Abshire, J. B.↗
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↗