Search NASA⌕ Search

SEARCH · Search NASA

Results for “random forest”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Symmetry in polarimetric remote sensing

Relationships among polarimetric backscattering coefficients are derived from the viewpoint of symmetry groups. For both reciprocal and non-reciprocal media, symmetry encountered in remote sensing due to reflection, rotation, azimuthal, and centrical symmetry groups is considered. The derived properties are general and valid to all scattering mechanisms, including volume and surface scatterings and their interactions, in a given symmetrical configuration. The scattering coefficients calculated from theoretical models for layer random media and rough surfaces are shown to obey the symmetry relations. Use of symmetry properties in remote sensing of structural and environmental responses of scattering media is also discussed. Orientations of spheroidal scatterers described by spherical, uniform, planophile, plagiothile, erectophile, and extremophile distributions are considered to derive their polarimetric backscattering characteristics. These distributions can be identified from the observed scattering coefficients by comparison with theoretical symmetry calculations. A new parameter is then defined to study scattering structures in geophysical media. Observations from polarimetric data acquired by the Jet Propulsion Laboratory airborne synthetic aperture radar over forests, sea ice, and sea surface are presented. Experimental evidences of the symmetry relationships are shown and their use in polarimetric remote sensing is illustrated. For forests, the coniferous forest in Mt. Shasta area (California) and mixed forest near Presque Isle (Maine) exhibit characteristics of the centrical symmetry at C-band. For sea ice in the Beaufort Sea, multi-year sea ice has a cross-polarized ratio e close to e(sub 0), calculated from symmetry, due to the randomness in the scattering structure. First-year sea ice has e much smaller than e(sub 0) due to the preferential alignment of the columnar structure of the ice. From polarimetric data of a sea surface in the Bering Sea, it is observed that e and e(sub 0) are increasing with incident angle and e is greater than e(sub 0) at L-band because of the directional feature of sea surface waves. Symmetry properties of geophysical media can also be used to calibrate polarimetric radars.

Nghiem, S. V.↗
NASA NTRS2011

Forest Biomass Mapping From Lidar and Radar Synergies

The use of lidar and radar instruments to measure forest structure attributes such as height and biomass at global scales is being considered for a future Earth Observation satellite mission, DESDynI (Deformation, Ecosystem Structure, and Dynamics of Ice). Large footprint lidar makes a direct measurement of the heights of scatterers in the illuminated footprint and can yield accurate information about the vertical profile of the canopy within lidar footprint samples. Synthetic Aperture Radar (SAR) is known to sense the canopy volume, especially at longer wavelengths and provides image data. Methods for biomass mapping by a combination of lidar sampling and radar mapping need to be developed. In this study, several issues in this respect were investigated using aircraft borne lidar and SAR data in Howland, Maine, USA. The stepwise regression selected the height indices rh50 and rh75 of the Laser Vegetation Imaging Sensor (LVIS) data for predicting field measured biomass with a R(exp 2) of 0.71 and RMSE of 31.33 Mg/ha. The above-ground biomass map generated from this regression model was considered to represent the true biomass of the area and used as a reference map since no better biomass map exists for the area. Random samples were taken from the biomass map and the correlation between the sampled biomass and co-located SAR signature was studied. The best models were used to extend the biomass from lidar samples into all forested areas in the study area, which mimics a procedure that could be used for the future DESDYnI Mission. It was found that depending on the data types used (quad-pol or dual-pol) the SAR data can predict the lidar biomass samples with R2 of 0.63-0.71, RMSE of 32.0-28.2 Mg/ha up to biomass levels of 200-250 Mg/ha. The mean biomass of the study area calculated from the biomass maps generated by lidar- SAR synergy 63 was within 10% of the reference biomass map derived from LVIS data. The results from this study are preliminary, but do show the potential of the combined use of lidar samples and radar imagery for forest biomass mapping. Various issues regarding lidar/radar data synergies for biomass mapping are discussed in the paper.

Sun, Guoqing↗
NASA NTRS1992

Remote sensing of Earth terrain

Research findings are summarized for projects dealing with the following: application of theoretical models to active and passive remote sensing of saline ice; radiative transfer theory for polarimetric remote sensing of pine forest; scattering of electromagnetic waves from a dense medium consisting of correlated Mie scatterers with size distribution and applications to dry snow; variance of phase fluctuations of waves propagating through a random medium; theoretical modeling for passive microwave remote sensing of earth terrain; polarimetric signatures of a canopy of dielectric cylinders based on first and second order vector radiative transfer theory; branching model for vegetation; polarimetric passive remote sensing of periodic surfaces; composite volume and surface scattering model; and radar image classification.

Kong, J. A.↗
NASA NTRS1993

Remote sensing of Earth terrain

Progress report on remote sensing of Earth terrain covering the period from Jan. to June 1993 is presented. Areas of research include: radiative transfer model for active and passive remote sensing of vegetation canopy; polarimetric thermal emission from rough ocean surfaces; polarimetric passive remote sensing of ocean wind vectors; polarimetric thermal emission from periodic water surfaces; layer model with tandom spheriodal scatterers for remote sensing of vegetation canopy; application of theoretical models to active and passive remote sensing of saline ice; radiative transfer theory for polarimetric remote sensing of pine forest; scattering of electromagnetic waves from a dense medium consisting of correlated mie scatterers with size distributions and applications to dry snow; variance of phase fluctuations of waves propagating through a random medium; polarimetric signatures of a canopy of dielectric cylinders based on first and second order vector radiative transfer theory; branching model for vegetation; polarimetric passive remote sensing of periodic surfaces; composite volume and surface scattering model; and radar image classification.

Kong, J. A.↗
NASA NTRS1999

Lidar Altimeter Measurements of Canopy Structure: Methods and Validation for Closed Canopy, Broadleaf Forests

Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.

Harding, D. J.↗
NASA NTRS1997

Evolution of Structure in the Intergalactic Medium and the Nature of the LY-Alpha Forest

We have performed a detailed statistical study of the evolution of structure in a photoionized intergalactic medium (IGM) using analytical simulations to extend the calculation into the mildly nonlinear density regime found to prevail at z = 3. Our work is based on a simple fundamental conjecture: that the probability distribution function of the density of baryonic diffuse matter in the universe is described by a lognormal (LN) random field. The LN distribution has several attractive features and follows plausibly from the assumption of initial linear Gaussian density and velocity fluctuations at arbitrarily early times. Starting with a suitably normalized power spectrum of primordial fluctuations in a universe dominated by cold dark matter (CDM), we compute the behavior of the baryonic matter, which moves slowly toward minima in the dark matter potential on scales larger than the Jeans length. We have computed two models that succeed in matching observations. One is a nonstandard CDM model with OMEGA = 1, h = 0.5, and GAMMA = 0.3, and the other is a low-density flat model with a cosmological constant (LCDM), with OMEGA = 0.4, OMEGA(sub LAMBDA) = 0.6, and h = 0.65. In both models, the variance of the density distribution function grows with time, reaching unity at about z = 4, where the simulation yields spectra that closely resemble the Ly-alpha forest absorption seen in the spectra of high-z quasars. The calculations also successfully predict the observed properties of the Ly-alpha forest clouds and their evolution from z = 4 down to at least z = 2, assuming a constant intensity for the metagalactic UV background over this redshift range. However, in our model the forest is not due to discrete clouds, but rather to fluctuations in a continuous intergalactic medium. At z = 3; typical clouds with measured neutral hydrogen column densities N(sub H I) = 10(exp 13.3), 10(exp 13.5), and 10(exp 11.5) /sq cm correspond to fluctuations with mean total densities approximately 10, 1, and 0.1 times the universal mean baryon density. Perhaps surprisingly, fluctuations whose amplitudes are less than or equal to the mean density still appear as "clouds" because in our model more than 70% of the volume of the IGM at z = 3 is filled with gas at densities below the mean value.

Bi, Hongguang↗
DOE OSTI2024

Quantifying the Variability of “Fixed-Width” Buffers on Harvested Lands in Western Oregon and Washington

In contemporary forest management, buffers of unharvested trees are left along streams to protect riparian and aquatic ecosystems. Buffer regulations often focus on specific minimum width requirements, which aid in straightforward regulation and application, but minimum widths also suggest buffered edges are uniform and contain little variability. Conceptual papers suggest that alternative buffer configurations may offer greater flexibility in landscape-level protection, increase forest complexity, and enhance aquatic and riparian biodiversity. However, before considering alternatives to fixed-width buffers, it is necessary to quantify the inherent variability in current buffer practices present on the landscape. In this study, we used aerial imagery to quantify variability of buffer widths in two hundred randomly selected recently harvested units on managed land in Oregon and Washington with both fish-bearing and non-fish-bearing sections of stream. Wider buffers on larger streams had a greater magnitude of variability, but when normalized by stream size, variation was greater in smaller streams, and overall, variation ranged from 25% to 50% of the mean width. Despite local variability, buffer widths rarely fell below 9.14 m (30 ft). The variation quantified here provides an initial measure of variability to inform future management, given emerging interest in variable retention buffers.

58 GEOSCIENCES↗
NASA NTRS1999

Unlocking the Keys to Vortex/Flame Interactions in Turbulent Gas-Jet Diffusion Flames--Dynamic Behavior Explored on the Space Shuttle

Most combustion processes in industrial applications (e.g., furnaces and engines) and in nature (e.g., forest fires) are turbulent. A better understanding of turbulent combustion could lead to improved combustor design, with enhanced efficiency and reduced emissions. Despite its importance, turbulent combustion is poorly understood because of its complexity. The rapidly changing and random behavior of such flames currently prevents detailed analysis, whether experimentally or computationally. However, it is possible to learn about the fundamental behavior of turbulent flames by exploring the controlled interaction of steady laminar flames and artificially induced flow vortices. These interactions are an inherent part of turbulent flames, and understanding them is essential to the characterization of turbulent combustion. Well-controlled and defined experiments of vortex interaction with laminar flames are not possible in normal gravity because of the interference of buoyancy- (i.e., gravity) induced vortices. Therefore, a joint microgravity study was established by researchers from the Science and Technology Development Corp. and the NASA Lewis Research Center. The experimental study culminated in the conduct of the Turbulent Gas-Jet Diffusion Flames (TGDF) Experiment on the STS-87 space shuttle mission in November 1997. The fully automated hardware, shown in photo, was designed and built at Lewis. During the mission, the experiment was housed in a Get Away Special (GAS) canister in the cargo bay.

Stocker, Dennis P.↗
NASA NTRS1993

Polarimetric Remote Sensing of Geophysical Medium Structures

Polarimetric remote sensing of structures in geophysical media is studied in this paper based on their symmetry properties. Orientations of spheroidal scatterers described by spherical, uniform, planophile, plagiothile, erectophile, and extremophile distributions are considered to derive their polarimetric backscattering characteristics. These distributions can be identified from the observed scattering coefficients by comparison with theoretical symmetry calculations. A new parameter is defined to study scattering structures in geophysical media. Experimental observations from polarimetric data acquired by the Jet Propulsion Laboratory airborne synthetic aperture radar over forests, sea ice, and sea surface are presented to illustrate the use of symmetry properties. For forests, the coniferous forest in Mount Shasta area and mixed forests neir Presque Isle show evidence of the centrical symmetry at C band. In sea ice from the Beaufort Sea, multiyear sea ice has a cross-polarized ratio e close to e(sub 0), calculated from symmetry, due to the randomness in the scattering structure. For first-year sea ice, e is much smaller than e(sub 0) as a result of preferential alignment of the columnar structure of the ice. From polarimetric data of a sea surface in the Bering sea, it is observed that e and e(sub 0) are increasing with incident angle and e is greater than e(sub 0) at L band because of the directional feature of sea surface waves. Use of symmetry properties of geophysical media for polarimetric radar calibration is also suggested.

Nghiem, S. V.↗
NASA NTRS2016

Modeling Forest Biomass and Growth: Coupling Long-Term Inventory and Lidar Data

Combining spatially-explicit long-term forest inventory and remotely sensed information from Light Detection and Ranging (LiDAR) datasets through statistical models can be a powerful tool for predicting and mapping above-ground biomass (AGB) at a range of geographic scales. We present and examine a novel modeling approach to improve prediction of AGB and estimate AGB growth using LiDAR data. The proposed model accommodates temporal misalignment between field measurements and remotely sensed data-a problem pervasive in such settings-by including multiple time-indexed measurements at plot locations to estimate AGB growth. We pursue a Bayesian modeling framework that allows for appropriately complex parameter associations and uncertainty propagation through to prediction. Specifically, we identify a space-varying coefficients model to predict and map AGB and its associated growth simultaneously. The proposed model is assessed using LiDAR data acquired from NASA Goddard's LiDAR, Hyper-spectral & Thermal imager and field inventory data from the Penobscot Experimental Forest in Bradley, Maine. The proposed model outperformed the time-invariant counterpart models in predictive performance as indicated by a substantial reduction in root mean squared error. The proposed model adequately accounts for temporal misalignment through the estimation of forest AGB growth and accommodates residual spatial dependence. Results from this analysis suggest that future AGB models informed using remotely sensed data, such as LiDAR, may be improved by adapting traditional modeling frameworks to account for temporal misalignment and spatial dependence using random effects.

Babcock, Chad↗
NASA NTRS1987

Continental land cover assessment using Landsat MSS data

A statistical procedure to assess level-II continental resources using Landsat MSS digital data is presented. The statistical procedure involves a two-stage cluster sample within a stratified random sample. The utility of this procedure is assessed by using it to estimate the areal extent of the conifer and hardwood resources of the continental U.S. National estimates of conifer and hardwood derived using this sampling procedure were within 3 percent of U.S. Forest Service (USFS) figures. According to the Landsat-based study, 11 percent of the country is conifer forest and 12 percent is hardwood. The corresponding USFS figures are 13 and 15 percent, respectively. Comparison of the MSS classification products and airphotos showed that the conifer cover class was correctly identified 74 percent of the time and hardwood 80 percent of the time. The average classification accuracy countrywide for the four cover types considered (conifer, hardwood, water, and 'other') is 74 percent, and the overall accuracy is 85 percent. The statistical procedure provides a method of incorporating Landsat MSS digital data as a second state for level-II continental resource assessment. Alternate data sources, e.g., satellite and aircraft photographic imagery, may also be used in conjunction with this statistical model.

Nelson, Ross↗
DOE OSTI2025

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2012

Geospatial Method for Computing Supplemental Multi-Decadal U.S. Coastal Land-Use and Land-Cover Classification Products, Using Landsat Data and C-CAP Products

This paper discusses the development and implementation of a geospatial data processing method and multi-decadal Landsat time series for computing general coastal U.S. land-use and land-cover (LULC) classifications and change products consisting of seven classes (water, barren, upland herbaceous, non-woody wetland, woody upland, woody wetland, and urban). Use of this approach extends the observational period of the NOAA-generated Coastal Change and Analysis Program (C-CAP) products by almost two decades, assuming the availability of one cloud free Landsat scene from any season for each targeted year. The Mobile Bay region in Alabama was used as a study area to develop, demonstrate, and validate the method that was applied to derive LULC products for nine dates at approximate five year intervals across a 34-year time span, using single dates of data for each classification in which forests were either leaf-on, leaf-off, or mixed senescent conditions. Classifications were computed and refined using decision rules in conjunction with unsupervised classification of Landsat data and C-CAP value-added products. Each classification's overall accuracy was assessed by comparing stratified random locations to available reference data, including higher spatial resolution satellite and aerial imagery, field survey data, and raw Landsat RGBs. Overall classification accuracies ranged from 83 to 91% with overall Kappa statistics ranging from 0.78 to 0.89. The accuracies are comparable to those from similar, generalized LULC products derived from C-CAP data. The Landsat MSS-based LULC product accuracies are similar to those from Landsat TM or ETM+ data. Accurate classifications were computed for all nine dates, yielding effective results regardless of season. This classification method yielded products that were used to compute LULC change products via additive GIS overlay techniques.

Spruce, J. P.↗
NASA NTRS2013

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI.

Comparison↗
DOE OSTI2024

Maximum a posteriori Ly α estimator (MAPLE): band power and covariance estimation of the 3D Ly α forest power spectrum

We present a novel maximum a posteriori estimator to jointly estimate band powers and the covariance of the three-dimensional power spectrum (P3D) of Ly $\alpha$ forest flux fluctuations, called MAPLE. Our Wiener-filter based algorithm reconstructs a window-deconvolved P3D in the presence of complex survey geometries typical for Ly $\alpha$ surveys that are sparsely sampled transverse to and densely sampled along the line of sight. We demonstrate our method on idealized Gaussian random fields with two selection functions: (i) a sparse sampling of 30 background sources per square degree designed to emulate the current Dark Energy Spectroscopic Instrument; (ii) a dense sampling of 900 background sources per square degree emulating the upcoming Prime Focus Spectrograph Galaxy Evolution Survey. Our proof-of-principle shows promise, especially since the algorithm can be extended to marginalize jointly over nuisance parameters and contaminants, i.e. offsets introduced by continuum fitting. Our code is implemented in JAX and is publicly available on GitHub.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2007

Potential of VIIRS Time Series Data for Aiding the USDA Forest Service Early Warning System for Forest Health Threats: A Gypsy Moth Defoliation Case Study

This report details one of three experiments performed during FY 2007 for the NASA RPC (Rapid Prototyping Capability) at Stennis Space Center. This RPC experiment assesses the potential of VIIRS (Visible/Infrared Imager/Radiometer Suite) and MODIS (Moderate Resolution Imaging Spectroradiometer) data for detecting and monitoring forest defoliation from the non-native Eurasian gypsy moth (Lymantria dispar). The intent of the RPC experiment was to assess the degree to which VIIRS data can provide forest disturbance monitoring information as an input to a forest threat EWS (Early Warning System) as compared to the level of information that can be obtained from MODIS data. The USDA Forest Service (USFS) plans to use MODIS products for generating broad-scaled, regional monitoring products as input to an EWS for forest health threat assessment. NASA SSC is helping the USFS to evaluate and integrate currently available satellite remote sensing technologies and data products for the EWS, including the use of MODIS products for regional monitoring of forest disturbance. Gypsy moth defoliation of the mid-Appalachian highland region was selected as a case study. Gypsy moth is one of eight major forest insect threats listed in the Healthy Forest Restoration Act (HFRA) of 2003; the gypsy moth threatens eastern U.S. hardwood forests, which are also a concern highlighted in the HFRA of 2003. This region was selected for the project because extensive gypsy moth defoliation occurred there over multiple years during the MODIS operational period. This RPC experiment is relevant to several nationally important mapping applications, including agricultural efficiency, coastal management, ecological forecasting, disaster management, and carbon management. In this experiment, MODIS data and VIIRS data simulated from MODIS were assessed for their ability to contribute broad, regional geospatial information on gypsy moth defoliation. Landsat and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) data were used to assess the quality of gypsy moth defoliation mapping products derived from MODIS data and from simulated VIIRS data. The project focused on use of data from MODIS Terra as opposed to MODIS Aqua mainly because only MODIS Terra data was collected during 2000 and 2001-years with comparatively high amounts of gypsy moth defoliation within the study area. The project assessed the quality of VIIRS data simulation products. Hyperion data was employed to assess the quality of MODIS-based VIIRS simulation datasets using image correlation analysis techniques. The ART (Application Research Toolbox) software was used for data simulation. Correlation analysis between MODIS-simulated VIIRS data and Hyperion-simulated VIIRS data for red, NIR (near-infrared), and NDVI (Normalized Difference Vegetation Index) image data products collectively indicate that useful, effective VIIRS simulations can be produced using Hyperion and MODIS data sources. The r(exp 2) for red, NIR, and NDVI products were 0.56, 0.63, and 0.62, respectively, indicating a moderately high correlation between the 2 data sources. Temporal decorrelation from different data acquisition times and image misregistration may have lowered correlation results. The RPC experiment also generated MODIS-based time series data products using the TSPT (Time Series Product Tool) software. Time series of simulated VIIRS NDVI products were produced at approximately 400-meter resolution GSD (Ground Sampling Distance) at nadir for comparison to MODIS NDVI products at either 250- or 500-meter GSD. The project also computed MODIS (MOD02) NDMI (Normalized Difference Moisture Index) products at 500-meter GSD for comparison to NDVI-based products. For each year during 2000-2006, MODIS and VIIRS (simulated from MOD02) time series were computed during the peak gypsy moth defoliation time frame in the study area (approximately June 10 through July 27). Gypsy moth defoliation mapping products from simated VIIRS and MOD02 time series were produced using multiple methods, including image classification and change detection via image differencing. The latter enabled an automated defoliation detection product computed using percent change in maximum NDVI for a peak defoliation period during 2001 compared to maximum NDVI across the entire 2000-2006 time frame. Final gypsy moth defoliation mapping products were assessed for accuracy using randomly sampled locations found on available geospatial reference data (Landsat and ASTER data in conjunction with defoliation map data from the USFS). Extensive gypsy moth defoliation patches were evident on screen displays of multitemporal color composites derived from MODIS data and from simulated VIIRS vegetation index data. Such defoliation was particularly evident for 2001, although widespread denuded forests were also seen for 2000 and 2003. These visualizations were validated using aforementioned reference data. Defoliation patches were visible on displays of MODIS-based NDVI and NDMI data. The viewing of apparent defoliation patches on all of these products necessitated adoption of a specialized temporal data processing method (e.g., maximum NDVI during the peak defoliation time frame). The frequency of cloud cover necessitated this approach. Multitemporal simulated VIIRS and MODIS Terra data both produced effective general classifications of defoliated forest versus other land cover. For 2001, the MOD02-simulated VIIRS 400-meter NDVI classification produced a similar yet slightly lower overall accuracy (87.28 percent with 0.72 Kappa) than the MOD02 250-meter NDVI classification (88.44 percent with 0.75 Kappa). The MOD13 250-meter NDVI classification had a lower overall accuracy (79.13 percent) and a much lower Kappa (0.46). The report discusses accuracy assessment results in much more detail, comparing overall classification and individual class accuracy statistics for simulated VIIRS 400-meter NDVI, MOD02 250-meter NDVI, MOD02-500 meter NDVI, MOD13 250-meter NDVI, and MOD02 500-meter NDMI classifications. Automated defoliation detection products from simulated VIIRS and MOD02 data for 2001 also yielded similar, relatively high overall classification accuracy (85.55 percent for the VIIRS 400-meter NDVI versus 87.28 percent for the MOD02 250-meter NDVI). In contrast, the USFS aerial sketch map of gypsy moth defoliation showed a lower overall classification accuracy at 73.64 percent. The overall classification Kappa values were also similar for the VIIRS (approximately 0.67 Kappa) versus the MOD02 (approximately 0.72 Kappa) automated defoliation detection product, which were much higher than the values exhibited by the USFS sketch map product (overall Kappa of approximately 0.47). The report provides additional details on the accuracy of automated gypsy moth defoliation detection products compared with USFS sketch maps. The results suggest that VIIRS data can be effectively simulated from MODIS data and that VIIRS data will produce gypsy moth defoliation mapping products that are similar to MODIS-based products. The results of the RPC experiment indicate that VIIRS and MODIS data products have good potential for integration into the forest threat EWS. The accuracy assessment was performed only for 2001 because of time constraints and a relative scarcity of cloud-free Landsat and ASTER data for the peak defoliation period of the other years in the 2000-2006 time series. Additional work should be performed to assess the accuracy of gypsy moth defoliation detection products for additional years.The study area (mid-Appalachian highlands) and application (gypsy moth forest defoliation) are not necessarily representative of all forested regions and of all forest threat disturbance agents. Additional work should be performed on other inland and coastal regions as well as for other major forest threats.

Spruce, Joseph P.↗
DOE OSTI2024

The 3D Lyman- α forest power spectrum from eBOSS DR16

We measure the three-dimensional power spectrum (P3D) of the transmitted flux in the Lyman-α (Ly α) forest using the complete extended Baryon Oscillation Spectroscopic Survey data release 16 (eBOSS DR16). This sample consists of ~205 000 quasar spectra in the redshift range 2 ≤ z ≤ 4 at an effective redshift z = 2.334. We propose a pair-count spectral estimator in configuration space, weighting each pair by exp( i k ∙ r), for wave vector k and pixel pair separation r, effectively measuring the anisotropic power spectrum without the need for fast Fourier transforms. This accounts for the window matrix in a tractable way, avoiding artefacts found in Fourier-transform based power spectrum estimators due to the sparse sampling transverse to the line of sight of Ly α skewers. We extensively test our pipeline on two sets of mocks: (i) idealized Gaussian random fields with a sparse sampling of Ly α skewers, and (ii) log-normal LyaCoLoRe mocks including realistic noise levels, the eBOSS survey geometry and contaminants. On eBOSS DR16 data, the Kaiser formula with a non-linear correction term obtained from hydrodynamic simulations yields a good fit to the power spectrum data in the range $(0.02 ≤ k ≤ 0.35)$ h Mpc -1 at the 1–2σ level with a covariance matrix derived from LyaCoLoRe mocks. We demonstrate a promising new approach for full-shape cosmological analyses of Ly α forest data from cosmological surveys such as eBOSS, the currently observing Dark Energy Spectroscopic Instrument and future surveys such as the Prime Focus Spectrograph, WEAVE-QSO, and 4MOST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Mock data sets for the Eboss and DESI Lyman-α forest surveys

We present a publicly-available code to generate sets of mock Lyman-α (Lyα) forest data that have realistic large-scale correlations including those due to the Baryonic Acoustic Oscillations (BAO). The primary purpose of these mocks is to test the analysis procedures of the Extended Baryon Oscillation Survey (eBOSS) and the Dark Energy Spectroscopy Instrument (DESI) surveys. The transmitted flux fraction, F(λ), of background quasars due to Lyα absorption in the intergalactic medium (IGM) is simulated using the Fluctuating Gunn-Petterson Approximation (FGPA) applied to Gaussian random fields produced through the use of fast Fourier transforms (FFT). The output includes the IGM-Lyα transmitted flux fraction along quasar lines of sight and a catalog of high-column-density systems appropriately placed at high-density regions of the IGM. This output serves as input to additional code that superimposes the IGM tranmission on realistic quasar spectra, adds absorption by high-column-density systems and metals, and simulates instrumental transmission and noise. Redshift space distortions (RSD) of the flux correlations are implemented by including the large-scale velocity-gradient field in the FGPA resulting in a correlation function of F(λ) that can be accurately predicted. One hundred realizations have been produced over the 14,000 deg 2 DESI survey footprint with 100 quasars per deg 2 . The analysis of these realizations shows that the correlations of F(λ) follows the prediction within the accuracy of eBOSS survey. Here, the most time-consuming part of the mock production occurs before application of the FGPA, and the existing pre-FGPA forests can be used to easily produce new mock sets with modified redshift-dependent bias parameters or observational conditions.

79 ASTRONOMY AND ASTROPHYSICS↗