Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial dependence”

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 937 records · Page 52

The NASA Soil Moisture Active Passive (SMAP) Mission - Science and Data Product Development Status

The Soil Moisture Active Passive (SMAP) mission, planned for launch in late 2014, has the objective of frequent, global mapping of near-surface soil moisture and its freeze-thaw state. The SMAP measurement system utilizes an L-band radar and radiometer sharing a rotating 6-meter mesh reflector antenna. The instruments will operate on a spacecraft in a 685 km polar orbit with 6am/6pm nodal crossings, viewing the surface at a constant 40-degree incidence angle with a 1000-km swath width, providing 3-day global coverage. Data from the instruments will yield global maps of soil moisture and freeze/thaw state at 10 km and 3 km resolutions, respectively, every two to three days. The 10-km soil moisture product will be generated using a combined radar and radiometer retrieval algorithm. SMAP will also provide a radiometer-only soil moisture product at 40-km spatial resolution and a radar-only soil moisture product at 3-km resolution. The relative accuracies of these products will vary regionally and will depend on surface characteristics such as vegetation water content, vegetation type, surface roughness, and landscape heterogeneity. The SMAP soil moisture and freeze/thaw measurements will enable significantly improved estimates of the fluxes of water, energy and carbon between the land and atmosphere. Soil moisture and freeze/thaw controls of these fluxes are key factors in the performance of models used for weather and climate predictions and for quantifYing the global carbon balance. Soil moisture measurements are also of importance in modeling and predicting extreme events such as floods and droughts. The algorithms and data products for SMAP are being developed in the SMAP Science Data System (SDS) Testbed. In the Testbed algorithms are developed and evaluated using simulated SMAP observations as well as observational data from current airborne and spaceborne L-band sensors including data from the SMOS and Aquarius missions. We report here on the development status of the SMAP data products. The Testbed simulations are designed to capture various sources of errors in the products including environment effects, instrument effects (nonideal aspects of the measurement system), and retrieval algorithm errors. The SMAP project has developed a Calibration and Validation (Cal/Val) Plan that is designed to support algorithm development (pre-launch) and data product validation (post-launch). A key component of the Cal/Val Plan is the identification, characterization, and instrumentation of sites that can be used to calibrate and validate the sensor data (Level l) and derived geophysical products (Level 2 and higher).

Nloku, E.↗

Monitoring Drought Conditions in the Navajo Nation Using NASA Earth Observations

The Navajo Nation, a 65,700 sq km Native American territory located in the southwestern United States, has been increasingly impacted by severe drought events and changes in climate. These events are coupled with a lack of domestic water infrastructure and economic resources, leaving approximately one-third of the population without access to potable water in their homes. Current methods of monitoring drought are dependent on state-based monthly Standardized Precipitation Index value maps calculated by the Western Regional Climate Center. However, these maps do not provide the spatial resolution needed to illustrate differences in drought severity across the vast Nation. To better understand and monitor drought events and drought regime changes in the Navajo Nation, this project created a geodatabase of historical climate information specific to the area, and a decision support tool to calculate average Standardized Precipitation Index values for user-specified areas. The tool and geodatabase use Tropical Rainfall Monitoring Mission (TRMM) and Global Precipitation Monitor (GPM) observed precipitation data and Parameter-elevation Relationships on Independent Slopes Model modeled historical precipitation data, as well as NASA's modeled Land Data Assimilation Systems deep soil moisture, evaporation, and transpiration data products. The geodatabase and decision support tool will allow resource managers in the Navajo Nation to utilize current and future NASA Earth observation data for increased decision-making capacity regarding future climate change impact on water resources.

Ly, Vickie↗

Overview of the Digitization Workflow Post Image Acquisition of Apollo Lunar and Antarctic Meteorite Samples Using Agisoft Photoscan for the NASA 3D Astromaterials Virtual Samples Collection

The 3D Virtual Astromaterials Samples (3DVAS) collection is a multi-year funded project to create a digital database of sixty Apollo Lunar and Antarctic Meteorite samples following non-destructive documentation conservation protocols. After initial image processing, the photos are evaluated and processed using unique structure-from-motion photogrammetric techniques in a high performance modelling software designed to create a 3D model from 2D images: Agisoft Photoscan Pro. Agisoft Photoscan Pro uses image processing algorithms and techniques originating in computer vision to resolve 3D models for accurate and detailed visualization of a subject. The software provides a stepwise process that is tailored per model based on spatial and specular reflectance properties, for example. The process includes: photo alignment, creation of a dense point cloud, mesh, and finally texture. Photo alignment is dependent on model properties. The 3DVAS process requires a special rotation platform with calibrated photogrammetric targets, specific distance rotation protocols, and a contrasting background for alignment and scale accuracy. As a result of the photographic process, alignment will complete with two mirrored hemispheres that, in a sense, represent the 2D images overlapping to create a 3D model. Each dense point cloud is analyzed with provided statistical measures in a gradual selection process to eliminate outliers. The point cloud is reduced to include only data valuable to the final model. When a precise dense point cloud is achieved, a mesh and texture are applied. Each model is scaled with scale bar accuracies within 100 microns. Each sample has its own intimate process for modelling; there is no standard for the parameters required in the final creation of a high resolution model. By processing multiple samples, a skill is gained in practice to allow a close definition of the original sample and will result in the most detailed version of the sample shell. This process completes one-fifth of the 3DVAS protocol for providing accurate digital documentation. Each model shell is merged with X-ray Computed Tomography data to create a full volumetric sample. All 3DVAS data will be served on NASA's Astromaterials Acquisition and Curation website with an early subset of data available in 2019 and the 3D Virtual Astromaterials Samples Collection launch in 2020.

Thomas, Andi B.↗

Projection Moire Interferometry for Rotorcraft Applications: Deformation Measurements of Active Twist Rotor Blades

Projection Moire Interferometry (PMI) has been used during wind tunnel tests to obtain azimuthally dependent blade bending and twist measurements for a 4-bladed Active Twist Rotor (ATR) system in simulated forward flight. The ATR concept offers a means to reduce rotor vibratory loads and noise by using piezoelectric active fiber composite actuators embedded in the blade structure to twist each blade as they rotate throughout the rotor azimuth. The twist imparted on the blades for blade control causes significant changes in blade loading, resulting in complex blade deformation consisting of coupled bending and twist. Measurement of this blade deformation is critical in understanding the overall behavior of the ATR system and the physical mechanisms causing the reduction in rotor loads and noise. PMI is a non-contacting, video-based optical measurement technique capable of obtaining spatially continuous structural deformation measurements over the entire object surface within the PMI system field-of-view. When applied to rotorcraft testing, PMI can be used to measure the azimuth-dependent blade bending and twist along the full span of the rotor blade. This paper presents the PMI technique as applied to rotorcraft testing, and provides results obtained during the ATR tests demonstrating the PMI system performance. PMI measurements acquired at select blade actuation conditions generating minimum and maximum rotor loads are provided to explore the interrelationship between rotor loads, blade bending, and twist.

Fleming, Gary A.↗

The Spatial Variation of Polar Rain Electrons and its Cause

It is generally accepted that field aligned electrons in the solar wind can follow field lines connected to Earth and precipitate in the polar ionosphere where they are known as polar rain. Few-hundred eV, field-aligned electrons of the solar wind "strahl" carry the interplanetary heat flux moving out from the sun and these electrons precipitate in either the northern or southern hemisphere depending on the magnetic field direction. These electrons produce enhanced polar rain in one hemisphere or the other although weaker polar rain is usually produced in the opposite hemisphere by whatever electrons are moving in the opposite direction. Although much evidence exists for this simple free entry mechanism, it has also long been known that there are spatial variations in the energies and intensities of the precipitating electrons. The present work compares electron distribution functions measured by the ACE spacecraft in the solar wind with those measured by the DMSP spacecraft at 800 km altitude over the polar cap. It is found that shifting the DMSP distribution functions in energy by amounts ranging from 10's to a few hundred eV produces quite good agreement with simultaneous ACE measurements. Over most of the polar cap this DMSP energy shift must be positive to achieve this agreement, suggesting the electrons have been decelerated by a field aligned potential as they move from the solar wind to low altitudes. The largest shifts occur on the nightside and on the dawn or dusk side, with the latter depending on the plasma convection pattern which is controlled by the orientation of the IMF. Nearer the cusp the shift is smaller or even negative. Since more massive tailward flowing magnetosheath ions are unable io follow the field lines into the magnetotail like the electrons, a field aligned potential is expected to develop to exclude low energy electrons and prevent an excessive charge imbalance. Such a potential would also produce the deceleration of those electrons that reach low altitudes. This improved understanding of polar rain should increase the utility of polar rain measurements as a diagnostic of the magnetosphere magnetic field configuration.

Fairfield, D. H.↗

The preferred-mode coherent structure in the near field of an axisymmetric jet with and without excitation

Visualization and phase-locked hot-wire measurements have been used to investigate the preferred-mode coherent structure in the near field of an axisymmetric jet for different Reynolds numbers and for the two limiting asymptotic states of the exit boundary layer: laminar and fully developed turbulent. Distributions of phase-averaged longitudinal and lateral velocities, coherent vorticity, coherent and incoherent turbulence Reynolds stresses, and coherent strain rate and production have been obtained over the spatial extent of the coherent structure for a Strouhal number of 0.3. It is found that these properties are independent of the initial conditions, but weakly dependent on the Reynolds number up to about 100,000. The naturally occurring coherent structure in the unexcited jet, educed via an optimized conditional sampling scheme, shows no clear dependence on the jet Reynolds number up to 800,000.

Hussain, A. K. M. F.↗

A simulation of rotor-stator interaction using the Euler equations and patched grids

An unsteady Euler code to study rotor-stator interaction problem was developed. The code uses patched grids that move relative to each other to simulate the motion of the rotor airfoils with respect to the stator airfoils. The Osher integration scheme is used in conjunction with an implicit relaxation approach. The scheme is second order accurate in space and time, and is also TVD in each spatial direction. The numerical results were found to be periodic in time, thus demonstrating the capability of the integration and zonal schemes in simulating periodic time dependent flow. The pressure contours obtained are almost oscillation free because of the TVD nature of the scheme. A new procedure was developed to simulate flows about bodies that move relative to each other. This capability should prove to be very useful in the areas of rotor-stator interaction, propeller-nacelle interaction, and helicopter rotor-fuselage interaction.

Rai, M. M.↗

Deformations and strains in a thick adherend lap joint

Displacement fields in a thick adherend lap joint were measured by high-sensitivity moire interferometry. Contour maps of in-plane U and V displacements were obtained across adhesive and adherend surfaces. Loads ranged from a modest load to a near-failure load. Quantitative results are given for displacements and strains in the adhesive and along the adhesive/adherend boundary lines. The results show nearly constant shear strain in the adhesive, nonlinear strains as a function of load or average shear stress, and viscoelastic or time-dependent response. Longitudinal normal strains in the adhesive are nearly two orders of magnitude less than the shear strains. With its subwavelength displacement resolution and high spatial resolution, moire interferometry is especially well suited for deformation studies of adhesive joints.

Post, D.↗

X-ray constraints on large scale clustering

HEAO-1 A2 all-sky survey data are analyzed to investigate the autocorrelation function of the surface brightness fluctuations of the extra-galactic 2-10-keV X-ray background on angular scales ranging from 3 to 27 deg. No significant signal has been detected on any scale. The derived upper limits set important bounds on correlated emission on scales of 10 Mpc or greater. X-ray data are found to be compatible with optical estimates of the two-point spatial correlation function of rich clusters of galaxies and with recent predictions for their cosmological evolution. Constraints on clustering of active galactic nuclei and on its dependence on cosmic time are analyzed in the light of results of recent optical studies.

De Zotti, G.↗

Radio interferometer observations of solar wind turbulence from the orbit of Helios to the solar corona

We report observations of Very Long Baseline Interferometer (VLBI) phase scintillations due to turbulence in the solar wind. The observations were made at 5.00 and 8.42 GHz with the Very Long Baseline Array (VLBA) on three dates in 1991 July and August. We observed the sources 0851 + 202 and 0735 + 178 at solar elongations ranging from 2.66 deg to 13.29 deg; the closest approach of the line of sight to the Sun ranged from 10 to 49.8 solar radii. We have also included previously unpublished 5 GHz VLBI phase scintillation measurements from 1989. These measurements probe solar wind density fluctuations on spatial scales from about 200 km to 2000 km. Our measurements are in quite good agreement with the Coles & Harmon model for the radio phase structure function, which was largely determined from observations on both shorter and larger spatial scales. Departures from the Coles & Harmon functions are attributable to day-to-day variations in the solar wind conditions. Phase scintillations at the greatest solar elongations are in very good agreement with extrapolated estimates from direct measurements made with the Helios spacecraft at slightly larger heliocentric distances. Thus there is a consistency between the in-situ spacecraft and radio sensing measurements of density turbulence. All of the VLBI data are consistent with a Kolmogorov spectrum for the density fluctuations, although at the closest elongations there may be excess power at small spatial scales. An advantage of interferometric techniques over other radio propagation measurements is that they provide a measure of the anisotropy of the irregularities. Our observations at closest approach (10 solar radii) show weak evidence for anisotropic, field-aligned density irregularities with an axial ratio of order 2. This degree of anisotropy would appear to be less than that measured at similar solar elongations but on smaller spatial scales by Armstrong and colleagues. Finally, a combination of the radio propagation data and in situ Helios measurements is used to determine the heliocentric distance dependence of the normalization coefficient of the density power spectrum, C(sup 2)(sub N). Modeling C(sup2)(sub N) varies as (R/solar radii)(sup -Kappa), we find Kappa = 3.72 +/- 0.30. Subject to a number of assumptions, this result is consistent with the conclusion of D. A. Roberts that turbulence within 0.28 astronomical units adheres to a WKB formula for turbulence amplitude as a function of heliocentric distance.

Spangler, Steven R.↗

Jet Noise Physics and Modeling Using First-principles Simulations

An extensive analysis of our jet DNS database has provided for the first time the complex correlations that are the core of many statistical jet noise models, including MGBK. We have also for the first time explicitly computed the noise from different components of a commonly used noise source as proposed in many modeling approaches. Key findings are: (1) While two-point (space and time) velocity statistics are well-fitted by decaying exponentials, even for our low-Reynolds-number jet, spatially integrated fourth-order space/retarded-time correlations, which constitute the noise "source" in MGBK, are instead well-fitted by Gaussians. The width of these Gaussians depends (by a factor of 2) on which components are considered. This is counter to current modeling practice, (2) A standard decomposition of the Lighthill source is shown by direct evaluation to be somewhat artificial since the noise from these nominally separate components is in fact highly correlated. We anticipate that the same will be the case for the Lilley source, and (3) The far-field sound is computed in a way that explicitly includes all quadrupole cancellations, yet evaluating the Lighthill integral for only a small part of the jet yields a far-field noise far louder than that from the whole jet due to missing nonquadrupole cancellations. Details of this study are discussed in a draft of a paper included as appendix A.

Freund, Jonathan B.↗

Low- to Mid-Latitude X-Ray Emission from Jupiter

The Chandra X-ray Observatory (CXO) observed Jupiter during the period 2003 February 24-26 for approx.40 hours (4 Jupiter rotations), using both the spectroscopy array of the Advanced CCD Imaging Spectrometer (ACIS-S) and the imaging array of the High-Resolution Camera (HRC-I). Two ACIS-S exposures, each approx.8.5 hr long, were separated by an HRC-I exposure of approx.20 hr. The low- to mid-latitude non-auroral disk X-ray emission is much more spatially uniform than the auroral emission. However, the low- to mid-latitude X-ray count rate shows a small but statistically significant hour angle dependence, and is higher in regions of relatively low surface magnetic field strength, confirming ROSAT results. In addition, the spectrum from the low surface field region shows an enhancement in the energy band 1.14- 1.38 keV, perhaps partly due to line emission from that region. Correlation of surface magnetic field strength with count rate is not found for the 2000 December HRC-I data, at a time when solar activity was high. The low- to mid-latitude disk X-ray count rate observed by the HRC-I in the 2003 February observation is about 50% of that observed in 2000 December, roughly consistent with a decrease in the solar activity index (F10.7 cm flux) by a similar amount over the same time period. The low- to mid-latitude X-ray emission does not show any oscillations similar to the -45 minute oscillations sometimes seen from the northern auroral zone. The temporal variation in Jupiter's non-auroral X-ray emission exhibits similarities to variations in solar X-ray flux observed by GOES and TIMED/SEE. The two ACIS-S 0.3-2 keV low- to mid-latitude X-ray spectra are harder than the auroral spectrum, and are different from each other at energies above 0.7 keV, showing variability in Jupiter s non-auroral X-ray emission on a time scale of a day. The 0.3-2.0 keV X-ray power emitted at low- to mid-latitudes is 0.21 GW and 0.39 GW for the first and second ACIS-S exposures, respectively. We suggest that X-ray emission from Jupiter's disk may be largely generated by solar X-rays resonantly and fluorescently scattered in its upper atmosphere, especially at times of high incident solar X-ray flux. However, the correlation of higher count rate with low surface magnetic-field strength indicates the presence of some secondary component, possibly ion precipitation from radiation belts closer to the planet than elsewhere at low- to mid-latitudes.

Bhardwaj, Anil↗

Low- to Middle-Latitude X-Ray Emission from Jupiter

The Chandra X-ray Observatory (CXO) observed Jupiter during the period 24-26 February 2003 for approx. 40 hours (4 Jupiter rotations), using both the spectroscopy array of the Advanced CCD Imaging Spectrometer (ACIS-S) and the imaging array of the High-Resolution Camera (HRC-I). Two ACIS-S exposures, each -8.5 hours long, were separated by an HRC-I exposure of approx. 20 hours. The low- to middle-latitude nonauroral disk X-ray emission is much more spatially uniform than the auroral emission. However, the low- to middle-latitude X-ray count rate shows a small but statistically significant hour angle dependence and depends on surface magnetic field strength. In addition, the X-ray spectra from regions corresponding to 3-5 gauss and 5-7 gauss surface fields show significant differences in the energy band 1.26-1.38 keV, perhaps partly due to line emission occurring in the 3-5 gauss region but not the 5-7 gauss region. A similar correlation of surface magnetic field strength with count rate is found for the 18 December 2000 HRC-I data, at a time when solar activity was high. The low- to middle-latitude disk X-ray count rate observed by the HRC-I in the February 2003 observation is about 50% of that observed in December 2000, roughly consistent with a decrease in the solar activity index (F10.7 cm flux) by a similar amount over the same time period. The low- to middle-latitude X-ray emission does not show any oscillations similar to the approx. 45 min oscillations sometimes seen from the northern auroral zone. The temporal variation in Jupiter's nonauroral X-ray emission exhibits similarities to variations in solar X-ray flux observed by GOES and TIMED/SEE. The two ACIS-S 0.3-2.0 keV low- to middle-latitude X-ray spectra are harder than the auroral spectrum and are different from each other at energies above 0.7 keV, showing variability in Jupiter's nonauroral X-ray emission on a timescale of a day. The 0.3-2.0 keV X-ray power emitted at low to middle latitudes is 0.21 GW and 0.39 GW for the first and second ACIS-S exposures, respectively. We suggest that X-ray emission from Jupiter's disk may be largely generated by the scattering and fluorescence of solar X rays in its upper atmosphere, especially at times of high incident solar X-ray flux. However, the dependence of count rate on surface magnetic-field strength may indicate the presence of some secondary component, possibly ion precipitation from radiation belts close to the planet.

Bhardwaj, Anil↗

A Critical Look at Deriving Monthly Aerosol Optical Depth from Satellite Data

Satellite-derived aerosol data sets, such as those provided by NASA's Moderate Resolution Imaging Spectroradiometer (MODIS) instruments, are greatly improving our understanding of global aerosol optical depth (AOD). Yet, there are sampling issues. MODIS specific orbital geometry, convolved with the need to avoid bright surfaces (glint, desert, clouds, etc.), means that AOD can be under- or over-sampled in places. When deriving downstream products, such as daily or monthly gridded AOD, one must consider the spatial and temporal density of the measurements relative to the gradients of the true AOD. Additionally, retrieval confidence criteria should be considered. Averaged products are highly dependent on choices made for data aggregation and weighting, and sampling errors can be further propagated when deriving regional or global mean AOD. Different choices for aggregation and weighting result in estimates of regional and global means varying by 30% or more. The impacts of a particular averaging algorithm vary by region and surface type and can be shown to represent different tolerance for clouds and retrieval confidence.

Levy, R. C.↗

GERMcode: A Stochastic Model for Space Radiation Risk Assessment

A new computer model, the GCR Event-based Risk Model code (GERMcode), was developed to describe biophysical events from high-energy protons and high charge and energy (HZE) particles that have been studied at the NASA Space Radiation Laboratory (NSRL) for the purpose of simulating space radiation biological effects. In the GERMcode, the biophysical description of the passage of HZE particles in tissue and shielding materials is made with a stochastic approach that includes both particle track structure and nuclear interactions. The GERMcode accounts for the major nuclear interaction processes of importance for describing heavy ion beams, including nuclear fragmentation, elastic scattering, and knockout-cascade processes by using the quantum multiple scattering fragmentation (QMSFRG) model. The QMSFRG model has been shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections. For NSRL applications, the GERMcode evaluates a set of biophysical properties, such as the Poisson distribution of particles or delta-ray hits for a given cellular area and particle dose, the radial dose on tissue, and the frequency distribution of energy deposition in a DNA volume. By utilizing the ProE/Fishbowl ray-tracing analysis, the GERMcode will be used as a bi-directional radiation transport model for future spacecraft shielding analysis in support of Mars mission risk assessments. Recent radiobiological experiments suggest the need for new approaches to risk assessment that include time-dependent biological events due to the signaling times for activation and relaxation of biological processes in cells and tissue. Thus, the tracking of the temporal and spatial distribution of events in tissue is a major goal of the GERMcode in support of the simulation of biological processes important in GCR risk assessments. In order to validate our approach, basic radiobiological responses such as cell survival curves, mutation, chromosomal aberrations, and representative mouse tumor induction curves are implemented into the GERMcode. Extension of these descriptions to other endpoints related to non-targeted effects and biochemical pathway responses will be discussed.

Kim, Myung-Hee Y.↗

Using Principal Component Analysis (PCA) to Speed up Radiative Transfer (RT) Computations

Multiple scattering RT calculations time-consuming. Need a speed improvement of about 1000 (for OCO)! Solution: Make use of redundancies in spectra. Correlated-k (Lacis and Wang, Lacis and Oinas, Goody et al, Fu and Liou) Problem: Assume that spectral variation of atmospheric optical properties spatially correlated at all points along optical path. High accuracy (HI) and 2-stream (2S) calculations have high correlation. Single scattering (SS) computations highly scenario-dependent, but not time consuming. Perform SS and 2S calculations at every wavelength. Perform small number of HI computations. Need to compute correction factor B at every wavelength.

data sets↗

The Impact of ARM on Climate Modeling

Climate models are among humanity's most ambitious and elaborate creations. They are designed to simulate the interactions of the atmosphere, ocean, land surface, and cryosphere on time scales far beyond the limits of deterministic predictability, and including the effects of time-dependent external forcings. The processes involved include radiative transfer, fluid dynamics, microphysics, and some aspects of geochemistry, biology, and ecology. The models explicitly simulate processes on spatial scales ranging from the circumference of the Earth down to one hundred kilometers or smaller, and implicitly include the effects of processes on even smaller scales down to a micron or so. The atmospheric component of a climate model can be called an atmospheric global circulation model (AGCM). In an AGCM, calculations are done on a three-dimensional grid, which in some of today's climate models consists of several million grid cells. For each grid cell, about a dozen variables are time-stepped as the model integrates forward from its initial conditions. These so-called prognostic variables have special importance because they are the only things that a model remembers from one time step to the next; everything else is recreated on each time step by starting from the prognostic variables and the boundary conditions. The prognostic variables typically include information about the mass of dry air, the temperature, the wind components, water vapor, various condensed-water species, and at least a few chemical species such as ozone. A good way to understand how climate models work is to consider the lengthy and complex process used to develop one. Lets imagine that a new AGCM is to be created, starting from a blank piece of paper. The model may be intended for a particular class of applications, e.g., high-resolution simulations on time scales of a few decades. Before a single line of code is written, the conceptual foundation of the model must be designed through a creative envisioning that starts from the intended application and is based on current understanding of how the atmosphere works and the inventory of mathematical methods available.

atmospheric circulation↗

Multibaseline TanDEM-X Mangrove Height Estimation: The Selection of the Vertical Wavenumber

We generated a large-scale mangrove forest height map using multiple TanDEM-X (TDX) interferometric synthetic aperture radar (InSAR) acquisitions with various spatial baselines in order to improve the height estimation accuracy across a wide range of forest heights. The forest height inversion using InSAR data is strongly dependent upon the vertical wavenumber (i.e., perpendicular baseline). First, we investigated the role of the vertical wavenumber in forest height inversion from InSAR data using the sensitivity of the interferometric (volume) coherence to forest height. We used corrected but lower resolution and accuracy Shuttle Radar Topography Mission (SRTM) mangrove height maps as a priori information over Akanda and Pongara National Parks in Gabon to estimate lower and upper boundaries of the vertical wavenumber over test sites from the measured coherence-to-height sensitivity. Only TDX acquisitions within the boundaries of the vertical wavenumber were selected and combined for multibaseline mangrove height inversion. Mangrove forest height was obtained with multibaseline TDX acquisitions and was validated against the reference height derived from field measurement data providing improvements in multibaseline inversion over existing height estimates (i.e., SRTM height) and single-baseline inversions (multibaseline inversion: r(exp 2) = 0.98, root mean square error (RMSE) of 2.73 m; SRTM height: r(exp 2) = 0.86, RMSE = 7.21 m; single-baseline inversions: r(exp 2) = 0.08-0.97, RMSE = 3.86-11.10 m). As a result, to accurately estimate forest heights over a wide range (3-60 m), multibaseline InSAR acquisitions (at least three different baselines) are needed to exclude biases associated with the vertical wavenumber in forest height inversion.

vertical wavenumber.↗