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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

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

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

DOE OSTI2020

On the Relationship Between the Marine Cold Air Outbreak M Parameter and Low-Level Cloud Heights in the Midlatitudes

Focusing on conditions of subsidence when low clouds are present, ground-based observations in both the North Atlantic and the Southern Ocean reveal strong relationships between cloud boundary (base and top heights) and different measures of lower tropospheric instability. The difference in potential temperature between the surface and 800 hPa (a metric called M) provides a stronger relationship than measures of inversion strength such as the lower tropospheric stability and estimated inversion strength. This is because (1) inversion strength itself does not correlate well with cloud boundaries, and (2) M contains information that appears important for cloud boundaries. These include the surface forcing through the use of sea surface rather than near-surface air temperature and an upper level close to the real cloud top. These results expand upon previous work on the importance of M as a predictor of cloud morphology. However, important differences are found in low-cloud conditions for the North Atlantic as compared to the Southern Ocean (for a given value of M): stronger inversions, deeper boundary layers, and much larger sea level pressures. Therefore, the relationship between cloud boundaries and M differs between the two regions. A general circulation model provides similar relationships as observed between M and both cloud top height and temperature but tends to place clouds higher and at colder temperatures than observed for a given M. This might cause issues with the representation of precipitation, cloud cover and radiation in the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Multimodal sensor fusion for real-time standoff estimation in directed energy deposition

In Laser Powder-based Direct Energy Deposition (LP-DED) systems, achieving consistency, precision and quality of produced parts requires tight control over printing parameters. One of the critical parameters is the standoff distance. Maintaining an optimal standoff height is crucial for achieving correct laser power density and powder catchment efficiency, as both laser and powder streams are focused at this distance. Here, this study introduces a novel approach using multimodal sensor fusion to predict standoff height in real-time. The proposed system integrates two low-profile, cost-effective sensors: an RGB coaxial camera and a high frequency and high dynamic range microphone. By utilizing a simple fully connected neural network, trained on a limited dataset, data fusion of these sensors allowed for the real-time prediction of the standoff height. The results demonstrate high resolution and accuracy of the predictions across multiple geometries and a wide range of standoff heights. This approach offers a simple, and cost-effective solution for real-time standoff height monitoring and lays the groundwork for future integration into commercial LP-DED systems.

42 ENGINEERING↗

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

DOE OSTI2025

Gas and dust dynamics in γ Cephei-type disks

Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗

DOE OSTI2026

Microstructure Validation of Graph Theory Model-Derived Cooling Rates in the Wire Arc Additive Manufacturing of ER70S-6 Steel

Wire arc additive manufacturing (WAAM) enables high-rate fabrication of large metallic components, but spatial variations in thermal history can lead to microstructural heterogeneity that requires efficient process models to evaluate. This study evaluates whether cooling rates extracted from a graph theory model (GTM)-based thermal simulation are consistent with the microstructural evolution observed in an ER70S-6 WAAM wall. Thermal histories from the model were analyzed at selected build heights, and cooling rates were extracted from the final thermal excursion through the austenite phase field. Microstructures at corresponding locations were characterized using electron backscatter diffraction (EBSD) to quantify grain size distributions, and pearlite interlamellar spacing was used as an additional indicator of cooling behavior. The modeled cooling rates were highest near the substrate and generally decreased with build height, consistent with the observed reduction in the fine grain fraction and the progressive shift in the grain size distribution as build height increased. Pearlite spacing trends also supported the modeled cooling rate variation. These results indicate that GTM-derived thermal histories can be post-processed into metallurgically meaningful cooling rate estimates for WAAM steel builds and linked to dataset specific empirical grain size distribution relationships for process–thermal history–microstructure assessment.

36 MATERIALS SCIENCE↗

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↗

NASA NTRS2020

Simultaneous Retrieval of Surface Roughness Parameters for Bare Soils from Combined Active-Passive Microwave SMAP Observations

An active-passive microwave retrieval algorithm for simultaneous determination of soil surface roughness parameters [vertical root-mean-square (RMS) height (s) and horizontal correlation length (l)] is presented for bare soils. The algorithm is based on active-passive microwave covariation, including the improved Integral Equation Method (I²EM), and is tested with global soil moisture active passive (SMAP) observations. The estimated retrieval results for s and l are overall consistent with values in the literature, indicating the validity of the proposed algorithm. Sensitivity analyses showed that the developed roughness retrieval algorithm is independent of permittivity for εₛ > 10 [-]. Furthermore, the physical model basis of this approach (I²EM) allows the application of different autocorrelation functions (ACF), such as Gaussian and exponential ACFs. Global roughness retrieval results confirm bare areas in deserts such as Sahara or Gobi. However, the type of ACF used within roughness parameter estimation is important. Retrieval results for the Gaussian ACF describe a rougher surface than retrieval results for the exponential ACF. No correlations were found between roughness results and the amount of precipitation or the soil texture, which could be due to the coarse spatial resolution of the SMAP data. The extension of this approach to vegetated soils is planned as an add-on study.

correlation length↗