SEARCH · Engineering Papers

Results for “distributed 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 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Decorrelation distance of snow in the Colorado River Basin

The problem of estimating areal averages from point measurement has been extensively studied by mining engineers and hydrologists. Its application to satellite measurements has recently been introduced. The semivariaogram has been used in many geostatistical applications to estimate spatial structures of observed properties, such as mineral distributions. An examination is made of snow variations in Colorado from daily snow data collected in 11 SNOTEL stations. The associated semivariogram is estimated. The objective is to estimate the spatial structure of the snow field so that the point data can be used for comparison with, and validation for, satellite measurements.

Chang, A. T. C.↗

NASA NTRS2012

A Probabilistic Mass Estimation Algorithm for a Novel 7- Channel Capacitive Sample Verification Sensor

A document describes an algorithm created to estimate the mass placed on a sample verification sensor (SVS) designed for lunar or planetary robotic sample return missions. A novel SVS measures the capacitance between a rigid bottom plate and an elastic top membrane in seven locations. As additional sample material (soil and/or small rocks) is placed on the top membrane, the deformation of the membrane increases the capacitance. The mass estimation algorithm addresses both the calibration of each SVS channel, and also addresses how to combine the capacitances read from each of the seven channels into a single mass estimate. The probabilistic approach combines the channels according to the variance observed during the training phase, and provides not only the mass estimate, but also a value for the certainty of the estimate. SVS capacitance data is collected for known masses under a wide variety of possible loading scenarios, though in all cases, the distribution of sample within the canister is expected to be approximately uniform. A capacitance-vs-mass curve is fitted to this data, and is subsequently used to determine the mass estimate for the single channel s capacitance reading during the measurement phase. This results in seven different mass estimates, one for each SVS channel. Moreover, the variance of the calibration data is used to place a Gaussian probability distribution function (pdf) around this mass estimate. To blend these seven estimates, the seven pdfs are combined into a single Gaussian distribution function, providing the final mean and variance of the estimate. This blending technique essentially takes the final estimate as an average of the estimates of the seven channels, weighted by the inverse of the channel s variance.

Wolf, Michael↗

NASA NTRS2021

SHARPEN: A Scheme to Restore the Distribution of Averaged Precipitation Fields

A key strategy in obtaining complete global coverage of high-resolution precipitation is to combine observations from multiple fields, such as the intermittent passive microwave observations, precipitation propagated in time using motion vectors, and geosynchronous infrared observations. These separate precipitation fields can be combined through weighted averaging, which produces estimates that are generally superior to the individual parent fields. However, the process of averaging changes the distribution of the precipitation values, leading to an increase in precipitating area and a decrease in the values of high precipitation rates, a phenomenon observed in IMERG. To mitigate this issue, we introduce a new scheme called SHARPEN (Scheme for Histogram Adjustment with Ranked Precipitation Estimates in the Neighborhood), which recovers the distribution of the averaged precipitation field based on the idea of quantile mapping applied to the local environment. When implemented in IMERG, precipitation estimates from SHARPEN exhibit a distribution that resembles that of the original instantaneous observations, with matching precipitating area and peak precipitation rates. Case studies demonstrate its improved ability in bridging between the parent precipitation fields. Evaluation against ground observations reveals a distinct improvement in precipitation detection skill, but also a slightly reduced correlation likely because of a sharper precipitation field. The increased computational demand of SHARPEN can be mitigated by striding over multiple grid boxes, which has only marginal impacts on the accuracy of the estimates. SHARPEN can be applied to any precipitation algorithm that produces an average from multiple input precipitation fields and is being considered for implementation in IMERG V07.

Jackson Tan↗

NASA NTRS2006

Estimation Of Organ Doses From Solar Particle Events For Future Space Exploration Missions

Radiation protection practices define the effective dose as a weighted sum of equivalent dose over major organ sites for radiation cancer risks. Since a crew personnel dosimeter does not make direct measurement of the effective dose, it has been estimated with skin-dose measurements and radiation transport codes for ISS and STS missions. If sufficient protection is not provided near solar maximum, the radiation risk can be significant due to exposure to sporadic solar particle events (SPEs) as well as to the continuous galactic cosmic radiation (GCR) on future exploratory-class and long-duration missions. For accurate estimates of overall fatal cancer risks from SPEs, the specific doses at various blood forming organs (BFOs) were considered, because proton fluences and doses vary considerably across marrow regions. Previous estimates of BFO doses from SPEs have used an average body-shielding distribution for the bone marrow based on the computerized anatomical man model (CAM). With the development of an 82-point body-shielding distribution at BFOs, the mean and variance of SPE doses in the major active marrow regions (head and neck, chest, abdomen, pelvis and thighs) will be presented. Consideration of the detailed distribution of bone marrow sites is one of many requirements to improve the estimation of effective doses for radiation cancer risks.

Kim, Myung-Hee↗

NASA NTRS2007

Tropical Cyclones Feed More Heavy Rain in a Warmer Climate

The possible linkage of tropical cyclones (TC) to global warming is a hotly debated scientific topic, with immense societal impacts. Most of the debate has been focused on the issue of uncertainty in the use of non-research quality data for long-term trend analyses, especially with regard to TC intensity provided by TC forecasting centers. On the other hand, it is well known that TCs are associated with heavy rain during the processes of genesis and intensification, and that there are growing evidences that rainfall characteristics (not total rainfall) are most likely to be affected by global warming. Yet, satellite rainfall data have not been exploited in any recent studies of linkage between tropical cyclones (TC) and global warming. This is mostly due to the large uncertainties associated with detection of long-term trend in satellite rainfall estimates over the ocean. This problem, as we demonstrate in this paper, can be alleviated by examining rainfall distribution, rather than rainfall total. This paper is the first to use research-quality, satellite-derived rainfall from TRMM and GPCP over the tropical oceans to estimate shift in rainfall distribution during the TC season, and its relationships with TCs, and sea surface temperature (SST) in the two major ocean basins, the northern Atlantic and the northern Pacific for 1979-2005. From the rainfall distribution, we derive the TC contributions to rainfall in various extreme rainfall categories as a function to time. Our results show a definitive trend indicating that TCs are contributing increasingly to heavier rain events, i.e., intense TC's are more frequent in the last 27 years. The TC contribution to top 5% heavy rain has nearly doubled in the last two decades in the North Atlantic, and has increased by about 10% in the North Pacific. The different rate of increase in TC contribution to heavy rain may be related to the different rates of different rate of expansion of the warm pool (SST >2S0 C) area in the two oceans.

Lau, K.-M.↗

NASA NTRS2014

Estimation of Crop Gross Primary Production (GPP): I. Impact of MODIS Observation Footprint and Impact of Vegetation BRDF Characteristics

Accurate estimation of gross primary production (GPP) is essential for carbon cycle and climate change studies. Three AmeriFlux crop sites of maize and soybean were selected for this study. Two of the sites were irrigated and the other one was rainfed. The normalized difference vegetation index (NDVI), the enhanced vegetation index (EVI), the green band chlorophyll index (CIgreen), and the green band wide dynamic range vegetation index (WDRVIgreen) were computed from the moderate resolution imaging spectroradiometer (MODIS) surface reflectance data. We examined the impacts of the MODIS observation footprint and the vegetation bidirectional reflectance distribution function (BRDF) on crop daily GPP estimation with the four spectral vegetation indices (VIs - NDVI, EVI, WDRVIgreen and CIgreen) where GPP was predicted with two linear models, with and without offset: GPP = a × VI × PAR and GPP = a × VI × PAR + b. Model performance was evaluated with coefficient of determination (R2), root mean square error (RMSE), and coefficient of variation (CV). The MODIS data were filtered into four categories and four experiments were conducted to assess the impacts. The first experiment included all observations. The second experiment only included observations with view zenith angle (VZA) ≤ 35◦ to constrain growth of the footprint size,which achieved a better grid cell match with the agricultural fields. The third experiment included only forward scatter observations with VZA ≤ 35◦. The fourth experiment included only backscatter observations with VZA ≤ 35◦. Overall, the EVI yielded the most consistently strong relationships to daily GPP under all examined conditions. The model GPP = a × VI × PAR + b had better performance than the model GPP = a × VI × PAR, and the offset was significant for most cases. Better performance was obtained for the irrigated field than its counterpart rainfed field. Comparison of experiment 2 vs. experiment 1 was used to examine the observation footprint impact whereas comparison of experiment 4 vs. experiment 3 was used to examine the BRDF impact. Changes in R2, RMSE,CV and changes in model coefficients "a" and "b" (experiment 2 vs. experiment 1; and experiment 4 vs. experiment 3) were indicators of the impacts. The second experiment produced better performance than the first experiment, increasing R2 (↑0.13) and reducing RMSE (↓0.68 g C m−2 d−1) and CV (↓9%). For each VI, the slope of GPP = a × VI × PAR in the second experiment for each crop type changed little while the slope and intercept of GPP = a × VI × PAR + b varied field by field. The CIgreen was least affected by the MODIS observation footprint in estimating crop daily GPP (R2, ↑0.08; RMSE, ↓0.42 g C m−2 d−1; and CV, ↓7%). Footprint most affected the NDVI (R2, ↑0.15; CV, ↓10%) and the EVI (RMSE, ↓0.84 g C m−2 d−1). The vegetation BRDF impact also caused variation of model performance and change of model coefficients. Significantly different slopes were obtained for forward vs. backscatter observations, especially for the CIgreen and the NDVI. Both the footprint impact and the BRDF impact varied with crop types, irrigation options, model options and VI options.

Chlorophyll↗

NASA NTRS2015

Construction of Polarimetric Radar-Based Reference Rain Maps for the Iowa Flood Studies Campaign

The Global Precipitation Measurement (GPM) Mission Iowa Flood Studies (IFloodS) campaign was conducted in central and northeastern Iowa during the months of April-June, 2013. Specific science objectives for IFloodS included quantification of uncertainties in satellite and ground-based estimates of precipitation, 4-D characterization of precipitation physical processes and associated parameters (e.g., size distributions, water contents, types, structure etc.), assessment of the impact of precipitation estimation uncertainty and physical processes on hydrologic predictive skill, and refinement of field observations and data analysis approaches as they pertain to future GPM integrated hydrologic validation and related field studies. In addition to field campaign archival of raw and processed satellite data (including precipitation products), key ground-based platforms such as the NASA NPOL S-band and D3R Ka/Ku-band dual-polarimetric radars, University of Iowa X-band dual-polarimetric radars, a large network of paired rain gauge platforms, and a large network of 2D Video and Parsivel disdrometers were deployed. In something of a canonical approach, the radar (NPOL in particular), gauge and disdrometer observational assets were deployed to create a consistent high-quality distributed (time and space sampling) radar-based ground "reference" rainfall dataset, with known uncertainties, that could be used for assessing the satellite-based precipitation products at a range of space/time scales. Subsequently, the impact of uncertainties in the satellite products could be evaluated relative to the ground-benchmark in coupled weather, land-surface and distributed hydrologic modeling frameworks as related to flood prediction. Relative to establishing the ground-based "benchmark", numerous avenues were pursued in the making and verification of IFloodS "reference" dual-polarimetric radar-based rain maps, and this study documents the process and results as they pertain specifically to efforts using the NPOL radar dataset. The initial portions of the "process" involved dual-polarimetric quality control procedures which employed standard phase and correlation-based approaches to removal of clutter and non-meteorological echo. Calculation of a scale-adaptive KDP was accomplished using the method of Wang and Chandrasekar (2009; J. Atmos. Oceanic Tech.). A dual-polarimetric blockage algorithm based on Lang et al. (2009; J. Atmos. Oceanic Tech.) was then implemented to correct radar reflectivity and differential reflectivity at low elevation angles. Next, hydrometeor identification algorithms were run to identify liquid and ice hydrometeors. After the quality control and data preparation steps were completed several different dual-polarimetric rain estimation algorithms were employed to estimate rainfall rates using rainfall scans collected approximately every two to three minutes throughout the campaign. These algorithms included a polarimetrically-tuned Z-R algorithm that adjusts for drop oscillations (via Bringi et al., 2004, J. Atmos. Oceanic Tech.), and several different hybrid polarimetric variable approaches, including one that made use of parameters tuned to IFloodS 2D Video Disdrometer measurements. Finally, a hybrid scan algorithm was designed to merge the rain rate estimates from multiple low level elevation angle scans (where blockages could not be appropriately corrected) in order to create individual low-level rain maps. Individual rain maps at each time step were subsequently accumulated over multiple time scales for comparison to gauge network data. The comparison results and overall error character depended strongly on rain event type, polarimetric estimator applied, and range from the radar. We will present the outcome of these comparisons and their impact on constructing composited "reference" rainfall maps at select time and space scales.

Radar↗

NASA NTRS2014

Data Analysis & Statistical Methods for Command File Errors

This paper explains current work on modeling for managing the risk of command file errors. It is focused on analyzing actual data from a JPL spaceflight mission to build models for evaluating and predicting error rates as a function of several key variables. We constructed a rich dataset by considering the number of errors, the number of files radiated, including the number commands and blocks in each file, as well as subjective estimates of workload and operational novelty. We have assessed these data using different curve fitting and distribution fitting techniques, such as multiple regression analysis, and maximum likelihood estimation to see how much of the variability in the error rates can be explained with these. We have also used goodness of fit testing strategies and principal component analysis to further assess our data. Finally, we constructed a model of expected error rates based on the what these statistics bore out as critical drivers to the error rate. This model allows project management to evaluate the error rate against a theoretically expected rate as well as anticipate future error rates.

Correlation Analysis↗

NASA NTRS2003

Statistical Properties of Maximum Likelihood Estimators of Power Law Spectra Information

A simple power law model consisting of a single spectral index, sigma(sub 2), is believed to be an adequate description of the galactic cosmic-ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at the knee energy, E(sub k), to a steeper spectral index sigma(sub 2) greater than sigma(sub 1) above E(sub k). The maximum likelihood (ML) procedure was developed for estimating the single parameter sigma(sub 1) of a simple power law energy spectrum and generalized to estimate the three spectral parameters of the broken power law energy spectrum from simulated detector responses and real cosmic-ray data. The statistical properties of the ML estimator were investigated and shown to have the three desirable properties: (Pl) consistency (asymptotically unbiased), (P2) efficiency (asymptotically attains the Cramer-Rao minimum variance bound), and (P3) asymptotically normally distributed, under a wide range of potential detector response functions. Attainment of these properties necessarily implies that the ML estimation procedure provides the best unbiased estimator possible. While simulation studies can easily determine if a given estimation procedure provides an unbiased estimate of the spectra information, and whether or not the estimator is approximately normally distributed, attainment of the Cramer-Rao bound (CRB) can only be ascertained by calculating the CRB for an assumed energy spectrum- detector response function combination, which can be quite formidable in practice. However, the effort in calculating the CRB is very worthwhile because it provides the necessary means to compare the efficiency of competing estimation techniques and, furthermore, provides a stopping rule in the search for the best unbiased estimator. Consequently, the CRB for both the simple and broken power law energy spectra are derived herein and the conditions under which they are stained in practice are investigated.

Howell, L. W., Jr.↗

NASA NTRS2002

Statistical Properties of Maximum Likelihood Estimators of Power Law Spectra Information

A simple power law model consisting of a single spectral index, a is believed to be an adequate description of the galactic cosmic-ray (GCR) proton flux at energies below 10(exp 13) eV, with a transition at the knee energy, E(sub k), to a steeper spectral index alpha(sub 2) greater than alpha(sub 1) above E(sub k). The Maximum likelihood (ML) procedure was developed for estimating the single parameter alpha(sub 1) of a simple power law energy spectrum and generalized to estimate the three spectral parameters of the broken power law energy spectrum from simulated detector responses and real cosmic-ray data. The statistical properties of the ML estimator were investigated and shown to have the three desirable properties: (P1) consistency (asymptotically unbiased). (P2) efficiency asymptotically attains the Cramer-Rao minimum variance bound), and (P3) asymptotically normally distributed, under a wide range of potential detector response functions. Attainment of these properties necessarily implies that the ML estimation procedure provides the best unbiased estimator possible. While simulation studies can easily determine if a given estimation procedure provides an unbiased estimate of the spectra information, and whether or not the estimator is approximately normally distributed, attainment of the Cramer-Rao bound (CRB) can only he ascertained by calculating the CRB for an assumed energy spectrum-detector response function combination, which can be quite formidable in practice. However. the effort in calculating the CRB is very worthwhile because it provides the necessary means to compare the efficiency of competing estimation techniques and, furthermore, provides a stopping rule in the search for the best unbiased estimator. Consequently, the CRB for both the simple and broken power law energy spectra are derived herein and the conditions under which they are attained in practice are investigated. The ML technique is then extended to estimate spectra information from an arbitrary number of astrophysics data sets produced by vastly different science instruments. This theory and its successful implementation will facilitate the interpretation of spectral information from multiple astrophysics missions and thereby permit the derivation of superior spectral parameter estimates based on the combination of data sets.

Howell, L. W.↗

NASA NTRS1975

High-energy galactic gamma radiation from cosmic rays concentrated in spiral arms

A model for the emission of high-energy (exceeding 100 Mev) gamma-rays from the galactic disk has been developed and compared with recent SAS-2 observations. In the calculation, it is assumed that (1) the high energy galactic gamma-rays result primarily from the interaction of cosmic rays with galactic matter, (2) the cosmic-ray density is proportional to the matter density on the scale of galactic arms, and (3) the matter in the Galaxy is distributed in a spiral pattern consistent with density-wave theory and experimental data on the matter distribution that is available, including the 21-cm H I line emission, continuum emission from H II regions, and data currently being used to estimate the H2 density. The calculated galactic-longitude distribution of gamma rays is in good agreement with the SAS-2 observations in relative shape and absolute flux. As a corollary, the nonuniform cosmic-ray distribution of this model tends to support the galactic origin of the fraction of cosmic rays which is important in the production of high-energy photons. Modifications of the basic model show that the gamma-ray flux is relatively sensitive to large variations of the assumed distribution of molecular hydrogen in the Galaxy.

Bignami, G. F.↗

NASA NTRS1994

Optical studies of interstellar material in low density regions of the Galaxy

We analyze high-resolution Na I and Ca II interstellar absorption line data obtained in an earlier spectroscopic survey of 57 stars along extended sight lines through the Galactic disk and halo. We find that the Na I lines trace a diffuse cloudy medium and the CA II lines trace both the cloudy medium and a more extended (intercloud) medium. High latitude and interarm sight lines that do not cross spiral arms have clouds that are more diffuse on average than those along sight lines that cross spiral arms. Spiral structure may play an important role in determinating the average absorption properties along extended sight lines and/or interesting physical differences may exist between sight lines that cross spiral arms and those that do not. These might include a harder radiation field and/or higher electron tempertures along the high latitude and 'clean' interarm sight lines. On average, 10% of the Ca II column density occurs at velocities forbidden by the Galactic rotation law by more than 10 km/s. In contrast, only a small precentage of the Na I column density occurs at these velocites. The Ca II to Na I ratio increases by a factor of 15 over forbidden velocities from 0 to 50 km/s and rises rapidly thereafter. A two component model of the Ca II column density per unit velocity over the range l = 325 deg to 360 deg indicates that two distinct distributions exists, one with sigma = 8 km/s and one with sigma = 21 km/s. As much as 60% of the Ca II column density at forbidden velocities may be associated with the faster distribution, which we attribute to warm intercloud material. We estimate expontential scale heights of 0.4-0.5 kpc for the neutral gas traced by the E(B-V), Na I, and H I distributions along the low density sight lines, and we find that Ca II has a larger scale height of 0.8 kpc.

Sembach, K. R.↗

NASA NTRS1989

Transmission overhaul and replacement predictions using Weibull and renewel theory

A method to estimate the frequency of transmission overhauls is presented. This method is based on the two-parameter Weibull statistical distribution for component life. A second method is presented to estimate the number of replacement components needed to support the transmission overhaul pattern. The second method is based on renewal theory. Confidence statistics are applied with both methods to improve the statistical estimate of sample behavior. A transmission example is also presented to illustrate the use of the methods. Transmission overhaul frequency and component replacement calculations are included in the example.

Savage, M.↗

NASA NTRS1989

Transmission overhaul and replacement predictions using Weibull and renewal theory

Savage, M.↗

NASA NTRS2005

Analysis of Upper Air, Ground and Remote Sensing Data for the ATLAS Field Campaign in San Juan, Puerto Rico

The Atlas San Juan Mission was conducted in February 2004 with the main objectives of observing the Urban Heat Island of San Juan, providing high resolution data of the land use for El Yunque Rain Forest and for calibrating remote sensors. The mission was coordinated with NASA staff members at Marsha& Stennis, Goddard, and Glenn. The Airborne Thermal and Land Applications Sensor (ATLAS) from NASA/Stennis, that operates in the visual and IR bands, was used as the main sensor and was flown over Puerto Rico in a Lear 23 jet plane. To support the data gathering effort by the ATLAS sensor, remote sensing observations and upper air soundings were conducted along with the deployment of a number of ground based weather stations and temperature sensors. This presentation focuses in the analysis of this complementary data for the Atlas San Juan Mission. Upper air data show that during the days of the mission the Caribbean mid and high atmospheres were relatively dry and highly stable reflecting positive surface lifted index, a necessary condition to conduct this suborbital campaign. Surface wind patterns at levels below 850mb were dominated by the easterly trades, while the jet stream at the edge of the troposphere dominated the westerly wind at levels above 500mb. The jet stream remained at high latitudes reducing the possibility of fronts. In consequence, only 8.4 mm of precipitation were reported during the entire mission. Observation of soundings located about 150 km apart reflected minimum variations of the boundary layer across the Island for levels below 850 meters and a uniform atmosphere for higher levels. The weather stations and the temperature sensors were placed at strategic locations to observe variations across the urban and rural landscapes. Time series plot of the stations' data show that heavily urbanized commercial areas have higher air temperatures than urban and suburban residential areas, and much higher temperatures than rural areas. Temperature differences [dT(U-R)] were obtained by subtracting the values of several stations h m a reference urban station, located m the commercial area of San Juan. These time series show that the UHI peaks during the morning between 10:00am and noon to an average of 4.5 C, a temporal pattern not previously observed in similar studies for continental cities. It is also observed a high variability of the UHI with the precipitation patterns even for short events. These results may be a reflection of a large land use density by low level buildings with an apparent absence of significant heat storage effects in the urban areas, and the importance of the surrounding soil and vegetation moisture in controlling the urban tropical climate. The ATLAS data was used to determine albedo and surface temperature patterns on a 10m scale for the study area. These data were used to calibrate the spatial distribution of the surface temperature when using remote sensing images from MODIS (Moderate Resolution Imaging Spectradiometer). Surface temperatures were estimated using the land surface temperature product MODII-L2 distributed by the Land Process Distributed Active Archive Center(LP DAAC). These results show the maximum, minimum and average temperatures in San Juan and in the entire Island at a resolution of 1 km. The information retrieved from MODIS for land surface temperatures reflected similar temporal and spatial variations as the weather stations and ATLAS measurements with a highest absolute offset of about 5 C due to the differences between surface and air temperatures.

Gonzalez, Jorge E.↗

NASA NTRS2018

Burden Versus Capability Statistical Analysis for Structural Probabilistic Risk Assessments

The purpose of the burden versus capability analysis is to analyze the ability of components to withstand the loads that they are subjected to. When designing a launch vehicle, there is always a trade-off for the strength of the components versus the weight of the vehicle. The vehicle needs to have some margin built in to the design, but this added margin should not add a significant amount of weight to the vehicle. When the material properties and limits are known, estimated loads can be used to ensure that the vehicle will survive launch loads. If the variation in the distributions can be quantified, the probability of failure can be estimated more accurately.

Green, Becky↗

NASA NTRS1991

An atlas of monthly mean distributions of GEOSAT sea surface height, SSMI surface wind speed, AVHRR/2 sea surface temperature, and ECMWF surface wind components during 1988

Monthly mean global distributions for 1988 are presented with a common color scale and geographical map. Distributions are included for sea surface height variation estimated from GEOSAT; surface wind speed estimated from the Special Sensor Microwave Imager on the Defense Meteorological Satellite Program spacecraft; sea surface temperature estimated from the Advanced Very High Resolution Radiometer on NOAA spacecrafts; and the Cartesian components of the 10m height wind vector computed by the European Center for Medium Range Weather Forecasting. Charts of monthly mean value, sampling distribution, and standard deviation value are displayed. Annual mean distributions are displayed.

Halpern, D.↗

NASA NTRS2021

Three-dimensional estimation of deciduous forest canopy structure and leaf area using multi-directional, leaf-on and leaf-off airborne lidar data

Airborne laser scanning (ALS) has been widely used to map gap probability and leaf area index (LAI) distribution at plot and landscape scales. As an indirect measurement, most ALS methods to estimate LAI combine waveform or point density information with supporting field measurements such as the leaf angle distribution, gap probability, or direct LAI measures. The development of a more independent estimation approach would facilitate more widespread use of existing ALS data to investigate patterns of forest structure and build realistic 3-D vegetation scenes to simulate remote sensing imagery and energy balance. Here, we develop a data processing workflow (named PVlad) using ALS point cloud apparent reflectance to estimate LAI and voxel-based leaf area density (LAD), aiming to reduce the need for associated field measurements such as the gap probability. The adaptation of the path volume (PV) concept derived from apparent reflectance integrates information from multi-directional ALS pulses, and quantifies the percentage exploration of each voxel for classification and occlusion correction, such that rigorous volumetric sampling approaches can be developed to derive LAI and LAD. The PVlad workflow was applied to discrete-return lidar data (Riegl VQ480i) acquired by NASA Goddard's LiDAR, Hyperspectral and Thermal Imager (G-LiHT) Airborne Imager during leaf-on (summer) and leaf-off (spring) conditions at the Smithsonian Environmental Research Center (SERC). The estimates of LAI and LAD captured structural differences between mature, logged, and intermediate-aged stands over eight deciduous forest plots. The derived LAI values were compared to field litter collection measurements, and the derived LAD vertical distribution was compared to the output of the VoxLAD model using terrestrial laser scan (TLS) field survey data. Using voxel sizes ranging from 0.5 m to 5 m, overall LAI estimation showed linear fitting coefficient bias and for 1 and 2 m voxel sizes, and vertical LAD distribution showed strong correlation with and for 0.5 and 1m voxel sizes. For every forest stand, upper-canopy LAD had a low variance for voxel sizes of ≤ . Application of PVlad to the G-LiHT and other similar ALS data archives enables the development of fine-resolution LAI map products, including voxelization of LAD for ecosystem science and radiative transfer simulations of remote sensing imagery or surface energy balance.

Tiangang Yin↗