Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampled value”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

The Origin of Carbon-bearing Volatiles in Surprise Valley Hot Springs in the Great Basin: Carbon Isotope and Water Chemistry Characterizations

There are numerous hydrothermal fields within the Great Basin of North America, some of which have been exploited for geothermal resources. With methane and other carbon-bearing compounds being observed, in some cases with high concentrations, however, their origins and formation conditions remain unknown. Thus, studying hydrothermal springs in this area provides us an opportunity to expand our knowledge of subsurface (bio)chemical processes that generate organic compounds in hydrothermal systems, and aid in future development and exploration of potential energy resources as well. While isotope measurement has long been used for recognition of their origins, there are several secondary processes that may generate variations in isotopic compositions: oxidation, re-equilibration of methane and other alkanes with CO2, mixing with compounds of other sources, etc. Therefore, in addition to isotopic analysis, other evidence, including water chemistry and rock compositions, are necessary to identify volatile compounds of different sources. Surprise Valley Hot Springs (SVHS, 41 deg 32'N, 120 deg 5'W), located in a typical basin and range province valley in northeastern California, is a terrestrial hydrothermal spring system of the Great Basin. Previous geophysical studies indicated the presence of clay-rich volcanic and sedimentary rocks of Tertiary age beneath the lava flows in late Tertiary and Quaternary. Water and gas samples were collected for a variety of chemical and isotope composition analyses, including in-situ pH, alkalinity, conductivity, oxidation reduction potential (ORP), major and trace elements, and C and H isotope measurements. Fluids issuing from SVHS can be classified as Na-(Cl)-SO4 type, with the major cation and anion being Na+ and SO4(2-), respectively. Thermodynamic calculation using ORP and major element data indicated that sulfate is the most dominant sulfur species, which is consistent with anion analysis results. Aquifer temperatures at depth estimated by both dissolved SiO2 and Na-K-Ca geothermometers are in the range of 125.0 to 135.4 C, and higher than the values measured at orifices (77.3 to 90.0 C). CO2 and homologs of straight chain alkanes (C1-C5) were identified in gas samples. Carbon isotope values of alkanes increase with carbon numbers. The C-13 fractionation between CO2 and dissolved inorganic carbon suggests they are out of carbon isotope equilibrium. The hypothesis regarding the formation of carbon-bearing compounds in SVHS may involve two processes: 1) Under high heat flow conditions which are caused by regional faulting and crustal extension, original high molecular weight organic compounds (kerogens) in clay-rich rocks decomposed to generate methane and other alkane homologs. 2) The SVHS area is associated with outflow structures, and distant from the heat source. Anaerobic oxidation of methane (AOM) with sulfate at shallow depth (< 90 C) is suggested as being responsible for the generation of CO2 in SVHS.

Fu, Qi↗
NASA NTRS2001

Dynamic Fatigue of ULE Glass

Ultra Low Expansion (ULE) glass is used in a number of applications which require a low thermal expansion coefficient. One such application is telescope mirror elements. An allowable stress can be calculated for this material based upon modulus of rupture data; however, this does not take into account the problem of delayed failure. Delayed failure, due to stress corrosion can significantly shorten the lifetime of a glass article. Knowledge of the factors governing the rate of subcritical flaw growth in a given environment enables the development of relations between lifetime, applied stress and failure probability for the material under study. Dynamic fatigue is one method of obtaining the necessary information to develop these relationships. In this study, the dynamic fatigue method was used to construct time-to-failure diagrams for both 230/270 ground and optically polished samples. The grinding and polishing process reduces the surface flaw size and subsurface damage, and relieves residual stress by removing materials with successively smaller grinding media. This resulted in an increase in the strength of the optic during the grinding and polishing sequence. There was also an increase in the lifetime due to grinding and polishing. It was found that using the fatigue parameters determined from the 230/270 grit surface are not significantly different from the optically polished values. Although the lower bound of the polished samples is more conservative, neither time-to-failure curves lie beyond the upper or lower bound of the confidence limits. Therefore, designers preferring conservative limits could use samples without residual stress present (polished samples) to determine the fatigue parameters and inert Weibull parameters from samples with the service condition surface, to determine time-to-failure of the optical element.

Tucker, Dennis S.↗
NASA NTRS1992

Aerodynamic roughness measured in the field and simulated in a wind tunnel

This study evaluates how well values of aerodynamic surface roughness, z sub 0, measured over scale models in wind tunnels correlate with values of z sub 0 measured at full scale in the field. A field experiment was conducted in which values of z sub 0 and u* (wind friction speed) were measured over three arrays of non-erodible roughness elements on a dry lake bed. Wind profiles were measured by ten anemometers on a 15 m mast under thermally neutral atmospheric conditions. Values of z sub 0 increased from .00014 m (dry lake bed only) to .026 m with increasing roughness element density. The three roughness element arrays were simulated at 1/10 and 1/20 scale in an open-circuit atmospheric boundary-layer wind tunnel. Velocities were measured with a boundary-layer pitot-tube rake from the same relative position within the scale model arrays as the anemometers were relative to the field arrays. Each array at each scale was sampled three times at five freestream velocities. Average values of z sub 0 for each model array at each scale were compared with full-scale values of z sub 0 obtained in the field. The field vs. wind tunnel correspondence of z sub 0 is found to be z sub 0 field = 0.2661 x (z sub(0 model) x scale(exp -1))exp .8159.

Sullivan, Robert↗
NASA NTRS2000

The Preparation Conditions of Chromium Doped ZnSe and Their Effect on the Infrared Luminescence Properties

We report the investigation by photoluminescence lifetime measurements of the near-IR emissions from a series of chromium-doped ZnSe samples, correlated to their preparation conditions. The samples were polycrystalline or single crystals prepared by post growth diffusion doping or single crystals doped during growth by the physical vapor transport method. Room temperature lifetime values between 6 and 8 micro seconds were measured for samples with Cr2+ from low 10(exp 17) to high 10(exp 18) / cubic cm range. Lifetime data taken down to 78 K was found to be rather temperature independent, reconfirming previous reports indicating a quantum yield of the corresponding emission of close to 100% at room temperature. A strong decrease in the room temperature lifetime was found for chromium concentrations higher than 10(exp 19) / cubic CM.

Burger, A.↗
NASA NTRS1990

On the determination of iridium in diverse geological samples employing HPGe-coincidence/NaI(Tl)-anticoincidence spectrometry

This paper reports the Ir content of a variety of geological samples determined by the high-purity Ge-coincidence/NaI(Tl)-anticoincidence gamma-ray spectrometry (henceforth referred to as coincidence/anticoincidence technique) and by the conventional INAA. The advantages of this technique are: (1) the Ir content of the samples is obtained (ppm to a fraction of ppb ranges) not only by the 468.1 keV peak as in the conventional INAA but also by the 784.6 keV and 920.9 keV sum peaks, which gives more confidence in the values obtained; and (2) it is well suited for the samples with high Compton background for which it is difficult to measure the Ir content by the conventional INAA technique. The practical sensitivity of this technique depends on the sample matrix. Under present experimental conditions, it varied from 0.1 ng for Mn nodules and 0.004 ng for Libyan Desert Glass. Iridium values obtained on small (about 1 microg) olivine grains demonstrate the potential application of this new technique to microsamples. The principle and methodology of this new technique as well as its advantages and disadvantages over the conventional INAA are discussed.

Murali, A. V.↗
NASA NTRS1978

ALC/50/ values for some polymeric materials

Apparent lethal concentrations for 50 per cent of the test animals within a 30-min exposure period (ALC/50/) were determined for seventeen samples of polymeric materials, using the screening test method. The materials evaluated included resin-glass composites, film composites, and miscellaneous resins. ALC(50) values, based on weight of original sample charged, ranged from 24 to 110 mg/l. Modified phenolic resins seemed to exhibit less toxicity than the baseline epoxy resins. Among the film composites evaluated, only flame modified polyvinyl fluoride appeared to exhibit less toxicity than the baseline polyvinyl fluoride film.

Hilado, C. J.↗
NASA NTRS2025

Erosion Results of the MISSE 9-15 Polymers and Composites Experiment 1-4 (PCE 1-4)

Polymers and other oxidizable materials on the exterior of spacecraft in the low Earth orbit (LEO) space environment can be eroded from reaction with atomic oxygen (AO). Therefore, in order to design durable spacecraft it is important to know the extent of erosion that will occur during a mission. This can be determined by knowing the LEO AO erosion yield, E y (volume loss per incident oxygen atom), of materials susceptible to AO reaction. In addition, recent flight experiments have shown that the AO E y can vary with the AO fluence and/or solar exposure. Therefore, obtaining AO E y data for materials flown on various spaceflight missions is important. NASA Glenn Research Center has flown numerous experiments as part of the Materials International Space Station Experiment (MISSE) missions on the exterior of the International Space Station (ISS) to characterize the LEO E y of polymers, composites, protective coatings, and other spacecraft materials. Recently, four Glenn experiments with 365 flight (F) samples were flown on ISS’s MISSE-Flight Facility (MISSE-FF). These experiments are the Polymers and Composites Experiment-1 (PCE-1) flown as part of the MISSE-9 mission, the PCE-2 flown as part of the MISSE-10 mission, the PCE-3 flown as part of the MISSE-12 and MISSE-15 missions, and the PCE-4 flown as part of the MISSE-13 mission. Although each experiment had numerous sample objectives, the primary objective was to determine the LEO AO E y of various spacecraft materials as a function of solar irradiation and AO fluence. This paper provides a summary of the erosion data for the PCE 1-4 AO E y samples. The AO E y for 150 samples flown in either the LEO ram, wake, zenith or nadir directions are provided. The AO ram fluence varied from 2.97×10 20 atoms/cm 2 after 0.89 years of direct space exposure (with relatively high levels of Si contamination) on MISSE 12 to 3.93×10 20 atoms/cm 2 after 1.17 years of direct space exposure on MISSE-10. The ram AO E y values for uncoated polymers range from 3.81×10 –25 cm 3 /atom for polytetrafluoroethylene (M9R-C20 F) exposed to an AO fluence of 3.44×10 20 atoms/cm 2 on MISSE-9 to 4.43×10 –23 cm 3 /atom for AO etched low density polyimide aerogel (M12R-C21 F) exposed to an AO fluence of 2.97×10 21 atoms/cm 2 on MISSE-12. Because of the low AO fluence and relatively high Si contamination, a number of PCE-3 wake samples experienced mass gain. Thus, AO E y values are not provided for these samples. Although there are calculated AO E y values for the zenith, wake and nadir samples, the ram AO E y for a particular material is a more reliable value in terms of AO exposure because the zenith, wake and nadir directions were exposed to either no or very little AO fluence and thus other space environmental factors (i.e. vacuum, thermal extremes and thermal cycling, and/or various types of radiation) are responsible for the mass loss.

Erosion yield↗
NASA NTRSDate not supplied

Improving Computational Efficiency of Prognostics Algorithms in Resource-Constrained Settings

In engineering and aerospace applications, it is vital to operational success to have insight into the expected performance and health of physical systems. The field of prognostics and health management provides quantitative methods for monitoring, predicting, and managing system health. Prognostics algorithms can be employed to assess the current state of a system, propagate the system throughout time, and predict potential anomalies or failures that may occur. While they can provide accurate prediction results, effective prognostics algorithms can be challenging to use in resource-constrained settings due to computational limitations and high computational latency, leading to obsolete predictions. Thus, computationally efficient and accurate algorithms are necessary for future remaining useful life predictions. In this work, we implement new algorithmic approaches for prediction, quantitatively compare them via a battery degradation use-case, and provide recommendations of potential improvements to a prognostics framework. One approach to prediction is through sampling, whereby the current state of a physical system is sampled many times and each sample is propagated forward until failure is reached, resulting in a distribution of failure values. To improve the efficiency of this process, we implemented five new algorithmic approaches to prediction, including three distinct sampling methods (standard Monte Carlo, Quasi-Monte Carlo, and Latin Hypercube Sampling), a variable time step algorithm, and a variable sample size algorithm. To compare the algorithms, we employ a variety of metrics designed specifically to analyze both computational efficiency and model accuracy. Our metrics include accuracy to compare the average predicted value to ground truth, mean absolute deviation to illustrate dispersion, specific percentile error to describe accuracy within a user-defined risk tolerance, and code run-time. To quantitatively analyze our results, we employ a use-case of degradation of a Lithium-ion battery. We use an electrochemistry-based model to describe the current health state of the battery, and implement our prediction algorithms to propagate forward in time until end-of-discharge (EOD) is reached. Notably, through this work it was found that none of our sampling approaches had a significant impact on computational efficiency or model accuracy in predicting EOD of the battery. We find that while the sampling methods are unique, the distributions they generate are similar, ultimately producing final predictions that are nearly identical. In exploring the effect of the time step within the prediction algorithm, we found that prediction accuracy was highly dependent on the time step used, and that implementing a variable time step within a particular prediction may provide an increase in computational efficiency while also maintaining prediction accuracy. Finally, implementing a variable sample size also affected prediction, and our results show that tuning both the magnitude and timing of the sample size adjustment can result in improved computation speed and maintained prediction accuracy. Taken together, our findings highlight the challenge of performing prognostics in resource-constrained settings, and illustrate the potential of developing new prediction algorithms to improve computational efficiency.

prognostics↗
NASA NTRS2011

Analysis of Mars Analogue Soil Samples Using Solid-Phase Microextraction, Organic Solvent Extraction and Gas Chromatography/Mass Spectrometry

Polycyclic aromatic hydrocarbons (PAHs) are robust and abundant molecules in extraterrestrial environments. They are found ubiquitously in the interstellar medium and have been identified in extracts of meteorites collected on Earth. PAHs are important target molecules for planetary exploration missions that investigate the organic inventory of planets, moons and small bodies. This study is part of an interdisciplinary preparation phase to search for organic molecules and life on Mars. We have investigated PAH compounds in desert soils to determine their composition, distribution and stability. Soil samples (Mars analogue soils) were collected at desert areas of Utah in the vicinity of the Mars Desert Research Station (MDRS), in the Arequipa region in Peru and from the Jutland region of Denmark. The aim of this study was to optimize the solid-phase microextraction (SPME) method for fast screening and determination of PAHs in soil samples. This method minimizes sample handling and preserves the chemical integrity of the sample. Complementary liquid extraction was used to obtain information on five- and six-ring PAH compounds. The measured concentrations of PAHs are, in general, very low, ranging from 1 to 60 ng g(sup -1). The texture of soils is mostly sandy loam with few samples being 100% silt. Collected soils are moderately basic with pH values of 8-9 except for the Salten Skov soil, which is slightly acidic. Although the diverse and variable microbial populations of the samples at the sample sites might have affected the levels and variety of PAHs detected, SPME appears to be a rapid, viable field sampling technique with implications for use on planetary missions.

Orzechowska, G. E.↗
NASA NTRS2015

Complete Tem-Tomography: 3D Structure of Gems Cluster

GEMS (glass with embedded metal and sulfide) grains in interplanetary dust particles (IDPs) are considered to be one of the ubiquitous and fundamental building blocks of solids in the Solar System. They have been considered to be interstellar silicate dust that survived various metamorphism or alteration processes in the protoplanetary disk but the elemental and isotopic composition measurements suggest that most of them have been formed in the protoplanetary disk as condensates from high temperature gas. This formation model is also supported by the formation of GEMS-like grains with respect to the size, mineral assemblage, texture and infrared spectrum by condensation experiments from mean GEMS composition materials. Previous GEMS studies were performed only with 2D observation by transmission electron microscopy (TEM) or scanning TEM (STEM). However, the 3D shape and structure of GEMS grains and the spatial distribution of Fe/FeS's has critical information about their formation and origin. Recently, the 3D structure of GEMS grains in ultrathin sections of cluster IDPs was revealed by electron tomography using a TEM/STEM (JEM-2100F, JEOL). However, CT images of thin sections mounted on Cu grids acquired by conventional TEM-tomography are limited to low tilt angles (e. g., less than absolute value of 75 deg. In fact, previous 3D TEM observations of GEMS were affected by some artifacts related to the limited tilt range in the TEM used. Complete tomographic images should be acquired by rotating the sample tilt angle over a range of more than absolute value of 80 deg otherwise the CT images lose their correct structures. In order to constrain the origin and formation process of GEMS grains more clearly, we performed complete electron tomography for GEMS grains. Here we report the sample preparation method we have developed for this study, and the preliminary results.

Matsuno, J.↗
NASA NTRS1989

Resistivity of pristine and intercalated graphite fiber epoxy composites

Laminar composites were fabricated from pristine and bromine intercalated Amoco P-55, P-75, and P-100 graphite fibers and Hysol-Grafil EAG101-1 film epoxy. The thickness and r.f. eddy current resistivity of several samples were measured at grid points and averaged point by point to obtain final values. Although the values obtained this way have high precision (less than 3 percent deviation), the resistivity values appear to be 20 to 90 percent higher than resistivities measured on high aspect ratio samples using multi-point techniques, and by those predicted by theory. The temperature dependence of the resistivity indicates that the fibers are neither damaged nor deintercalated by the composite fabrication process. The resistivity of the composites is a function of sample thickness (i.e., resin content). Composite resistivity is dominated by fiber resistivity, so lowering the resistivity of the fibers, either through increased graphitization or intercalation, results in a lower composite resistivity. A modification of the simple rule of mixtures model appears to predict the conductivity of high aspect ratio samples measured along a fiber direction, but a directional dependence appears which is not predicted by the theory. The resistivity of these materials is clearly more complex than that of homogeneous materials.

Gaier, James R.↗
NASA NTRS1991

Resistivity of pristine and intercalated graphite fiber epoxy composites

Laminar composites were fabricated from pristine and bromine intercalated Amoco P-55, P-75, and P-100 graphite fibers and Hysol-Grafil EAG101-1 film epoxy. The thickness and r.f. eddy current resistivity of several samples were measured at grid points and averaged point by point to obtain final values. Although the values obtained this way have high precision (less than 3 percent deviation), the resistivity values appear to be 20 to 90 percent higher than resistivities measured on high aspect ratio samples using multi-point techniques, and by those predicted by theory. The temperature dependence of the resistivity indicates that the fibers are neither damaged nor deintercalated by the composite fabrication process. The resistivity of the composites is a function of sample thickness (i.e., resin content). Composite resistivity is dominated by fiber resistivity, so lowering the resistivity of the fibers, either through increased graphitization or intercalation, results in a lower composite resistivity. A modification of the simple rule of mixtures model appears to predict the conductivity of high aspect ratio samples measured along a fiber direction, but a directional dependence appears which is not predicted by the theory. The resistivity of these materials is clearly more complex than that of homogeneous materials.

Gaier, James R.↗
NASA NTRS1995

Measurement of photoemission and secondary emission from laboratory dust grains

The overall goal of this project is experimentally determine the emission properties of dust grains in order to provide theorists and modelers with an accurate data base to use in codes that predict the charging of grains in various plasma environments encountered in the magnetospheres of the planets. In general these modelers use values which have been measured on planar, bulk samples of the materials in question. The large enhancements expected due to the small size of grains can have a dramatic impact upon the predictions and the ultimate utility of these predictions. The first experimental measurement of energy resolved profiles of the secondary electron emission coefficient, 6, of sub-micron diameter particles has been accomplished. Bismuth particles in the size range of .022 to .165 micrometers were generated in a moderate pressure vacuum oven (average size is a function of oven temperature and pressure) and introduced into a high vacuum chamber where they interacted with a high energy electron beam (0.4 to 20 keV). Large enhancements in emission were observed with a peak value, delta(sub max) = 4. 5 measured for the ensemble of particles with a mean size of .022 micrometers. This is in contrast to the published value, delta(sub max) = 1.2, for bulk bismuth. The observed profiles are in general agreement with recent theoretical predictions made by Chow et al. at UCSD.

Hazelton, Robert C.↗
NASA NTRS2004

Computation of Variability in the Average Thermal and Mechanical Properties of a Melt-infiltrated Sic/sic Composite

Thermal conductivity and tensile properties (elastic modulus, proportional limit strength, inplane tensile strength, and strain to failure) of a SiC/SiC composite were experimentally determined at 816 and 1204 C. Tests were performed at room temperature on the same material to obtain interlaminar shear and tensile strength. For each thermal and mechanical property 24 tests were conducted to capture the variation. A random sampling method was used to quantify the variability exhibited by the mean values of the thermal and mechanical properties. The minimum number of tests required to characterize the mean value of each thermal or mechanical property for the composite was determined by varying the group size of the random sample.

Kalluri, Sreeramesh↗
NASA NTRS2014

Technique for Performing Dielectric Property Measurements at Microwave Frequencies

A method, system, apparatus, and computer readable medium has been provided with the ability to obtain a complex permittivity dielectric or a complex permeability micron of a sample in a cavity. One or more complex-valued resonance frequencies f(sub m) of the cavity, wherein each f(sub m) is a measurement, are obtained. Maxwell's equations are solved exactly for dielectric, and/or micron, using the f(sub m) as known quantities, thereby obtaining the dielectric and/or micron of the sample.

Barmatz, Martin B.↗
NASA NTRS2017

Osiris-REx Spacecraft Current Status and Forward Plans

The NASA New Frontiers OSIRIS-REx spacecraft executed a flawless launch on September 8, 2016 to begin its 23-month journey to near-Earth asteroid (101955). The primary objective of the OSIRIS-REx mission is to collect and return to Earth a pristine sample of regolith from the asteroid surface. The sampling event will occur after a two-year period of remote sensing that will ensure a high probability of successful sampling of a region on the asteroid surface having high science value and within well-defined geological context. The OSIRIS-REx instrument payload includes three high-resolution cameras (OCAMS), a visible and near-infrared spectrometer (OVIRS), a thermal imaging spectrometer (OTES), an X-ray imaging spectrometer (REXIS), and a laser altimeter (OLA). As the spacecraft follows its nominal outbound-cruise trajectory, the propulsion, power, communications, and science instruments have undergone basic functional tests, with no major issues. Outbound cruise science investigations include a search for Earth Trojan asteroids as the spacecraft approaches the Sun-Earth L4 Lagrangian point in February 2017. Additional instrument checkouts and calibrations will be carried out during the Earth gravity assist maneuver in September 2017. During the Earth-moon flyby, visual and spectral images will be acquired to validate instrument command sequences planned for Bennu remote sensing. The asteroid Bennu remote sensing campaign will yield high resolution maps of the temperature and thermal inertia, distributions of major minerals and concentrations of organic matter across the asteroid surface. A high resolution 3d shape model including local surface slopes and a high-resolution gravity field will also be determined. Together, these data will be used to generate four separate maps that will be used to select the sampling site(s). The Safety map will identify hazardous and safe operational regions on the asteroid surface. The Deliverability map will quantify the accuracy with which the navigation team can deliver the spacecraft to and from specific sites on the asteroid surface. The Sampleability map quantifies the regolith properties, providing an estimation of how much material would be sampled at different points on the surface. The final Science Value map synthesizes the chemical, mineralogical, and geological, observations to identify the areas of the asteroid surface with the highest science value. Here, priority is given to organic, water-rich regions that have been minimally altered by surface processes. Asteroid surface samples will be acquired with a touch-and-go sample acquisition system (TAGSAM) that uses high purity pressurized N2 gas to mobilize regolith into a stainless steel canister. Although the mission requirement is to collect at least 60 g of material, tests of the TAGSAM routinely exceeded 300 g of simulant in micro-gravity tests. After acquiring the sample, the spacecraft will depart Bennu in 2021 to begin its return journey, with the sample return capsule landing at the Utah Test and Training Range on September 23, 2023. The OSIRIS-REx science team will carry out a series of detailed chemical, mineralogical, isotopic, and spectral studies that will be used to determine the origin and history of Bennu and to relate high spatial resolution sample studies to the global geological context from remote sensing. The outline of the sample analysis plan is described in a companion abstract.

Messenger, Scott↗
NASA NTRS1984

A supercluster catalog

The largest available complete redshift sample of rich clusters of galaxies is used to construct a three-dimensional catalog of superclusters complete to redhsift of about 0.08 or less. The criterion of membership in the catalog is the density enhancement of clusters over the mean density; the catalog is presented for various values of this parameter. The dependence of the identified sample of superclusters on the overdensity is studied, and the sensitivity of the identified superclusters to various quantities and selection procedures is investigated. The validity of the identification procedure is established by comparisons with random catalogs.

Bahcall, N. A.↗
NASA NTRS2008

Relating a Spectral Index from MODIS and Tower-based Measurements to Ecosystem Light Use Efficiency for a Fluxnet-Canada Coniferous Forest

As part of the North American Carbon Program effort to quantify the terrestrial carbon budget of North America, we have been examining the possibility of retrieving ecosystem light use efficiency (LUE, the carbon sequestered per unit photosynthetically active radiation) directly from satellite observations. Our novel approach has been to compare LUE derived from tower fluxes with LUE estimated using spectral indices computed from MODIS satellite observations over forests in the Fluxnet-Canada Research Network, using the MODIS narrow ocean bands acquired over land. We matched carbon flux data collected around the time of the MODIS mid-day overpass for over one hundred relatively clear days in five years (2001-2006) from a mature Douglas fir forest in British Columbia. We also examined hyperspectral reflectance data collected diurnally from the tower in conjunction with the eddy correlation fluxes and meteorological measurements made throughout the 2006 growing season at this site. The tower-based flux data provided an opportunity to examine diurnal and seasonal LUE processes and their relationship to spectral indices at the scale of the forest stand. We evaluated LUE in conjunction with the Photochemical Reflectance Index (PRI), a normalized difference spectral index that uses 531 nm and a reference band to capture responses to high light induced stress afforded by the xanthophyll cycle. Canopy structure information, retrieved from airborne laser scanning radar (LiDAR) observations, was used to partition the forest canopy into sunlit and shaded fractions throughout the day, on numerous days during 2006. At each observation period throughout a day, the PRI was examined for the sunlit, shaded, and intermediate canopy segments defined by their instantaneous position relative to the solar principal plane (SPP). The sunlit sector was associated with the illumination "hotspot" (the reflectance backscatter maximum), the shaded sector with the "cold or dark spot" (the reflectance forward scatter minimum), while the intermediate, mixed sunlit/shade sector was located in the cross-plane to the SPP. The PRI indices clearly captured the differences in leaf groups, with sunlit foliage exhibiting the lowest values on sunny days throughout the 2006 season. When tower-based canopy-level LUE was recalculated to estimate foliage-based values (LUE(sub foilage) for the three foliage groups under their incident light environments, a strong linear relationship for PRI:LUE(sub foilage) was demonstrated (0.6 less than or equal to r(sup 2) less than or equal to 0.8, n=822, P<0.0001). The MODIS data represent relatively large areas when acquired at nadir (approx.1 sq km) or at variable off-nadir view angles (greater than or equal to 1 sq km) looking forward or aft. Nevertheless, a similar relationship between MODIS PRI and tower-based LUE was obtained from satellite observations (r(sup 2) = 0.76, n=105, P= 0.026) when the azimuth offsets from the SPP for off-nadir observations were considered. At this relatively high latitude of 50 degrees, the MODIS directional observations were offset from the SPP by approximately 50 degrees, but still represented backscatter or forward scatter sectors of the bidirectional reflectance distribution function (BRDF). The backscatter observations sampled the sunlit forest and provided lower PRI values, in general, than the forward scatter observations from the shaded forest. Since the hotspot and darkspot were not typically directly observed, the dynamic range for MODIS PRI was less than that observed in the SPP at the canopy level; therefore, MODIS PRI values were more similar to those observed in sifu in the BRDF cross-plane. While not ideal in terms of spatial resolution or optimal viewing configuration, the MODIS observations nevertheless provide a means to monitor forest under stress using narrow spectral band indices and off-nadir observations. This research has stimulated several spin-off studies for remote sensinf LUE, and demonstrates the importance of the connection between ecosystem structure and physiological function.

Middleton, Elizabeth M.↗