SEARCH · Engineering Papers

Results for “spectral element method”

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 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2014

WE-G-18C-07: Accelerated Water/fat Separation in MRI for Radiotherapy Planning Using Multi-Band Imaging Techniques

Purpose: Dixon sequences are used to characterize disease processes, obtain good fat or water separation in cases where fat suppression fails and to obtain pseudo-CT datasets. Dixon's method uses at least two images acquired with different echo times and thus requires prolonged acquisition times. To overcome associated problems (e.g., for DCE/cine-MRI), we propose to use a method for water/fat separation based on spectrally selective RF pulses. Methods: Two alternating RF pulses were used, that imposes a fat selective phase cycling over the phase encoding lines, which results in a spatial shift for fat in the reconstructed image, identical to that in CAIPIRINHA. Associated aliasing artefacts were resolved using the encoding power of a multi-element receiver array, analogous to SENSE. In vivo measurements were performed on a 1.5T clinical MR-scanner in a healthy volunteer's legs, using a four channel receiver coil. Gradient echo images were acquired with TE/TR = 2.3/4.7ms, flip angle 20°, FOV 45×22.5cm{sup 2}, matrix 480×216, slice thickness 5mm. Dixon images were acquired with TE,1/TE,2/TR=2.2/4.6/7ms. All image reconstructions were done in Matlab using the ReconFrame toolbox (Gyrotools, Zurich, CH). Results: RF pulse alternation yields a fat image offset from the water image. Hence the water and fat images fold over, which is resolved using in-plane SENSE reconstruction. Using the proposed technique, we achieved excellent water/fat separation comparable to Dixon images, while acquiring images at only one echo time. Conclusion: The proposed technique yields both inphase water and fat images at arbitrary echo times and requires only one measurement, thereby shortening the acquisition time by a factor 2. In future work the technique may be extended to a multi-band water/fat separation sequence that is able to achieve single point water/fat separation in multiple slices at once and hence yields higher speed-up factors.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

SOLAR OPACITY CALCULATIONS USING THE SUPER-TRANSITION-ARRAY METHOD

A new opacity model has been developed based on the Super-Transition-Array (STA) method for the calculation of monochromatic opacities of plasmas in local thermodynamic equilibrium. The atomic code, named STAR (STA-Revised), is described and used to calculate spectral opacities for a solar model implementing the recent AGSS09 composition. Calculations are carried out throughout the solar radiative zone. The relative contributions of different chemical elements and atomic processes to the total Rosseland mean opacity are analyzed in detail. Monochromatic opacities and charge-state distributions are compared with the widely used Opacity Project (OP) code, for several elements near the radiation–convection interface. STAR Rosseland opacities for the solar mixture show a very good agreement with OP and the OPAL opacity code throughout the radiation zone. Finally, an explicit STA calculation was performed of the full AGSS09 photospheric mixture, including all heavy metals. It was shown that, due to their extremely low abundance, and despite being very good photon absorbers, the heavy elements do not affect the Rosseland opacity.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1991

Accuracy and speed in computing the Chebyshev collocation derivative

We studied several algorithms for computing the Chebyshev spectral derivative and compare their roundoff error. For a large number of collocation points, the elements of the Chebyshev differentiation matrix, if constructed in the usual way, are not computed accurately. A subtle cause is is found to account for the poor accuracy when computing the derivative by the matrix-vector multiplication method. Methods for accurately computing the elements of the matrix are presented, and we find that if the entities of the matrix are computed accurately, the roundoff error of the matrix-vector multiplication is as small as that of the transform-recursion algorithm. Results of CPU time usage are shown for several different algorithms for computing the derivative by the Chebyshev collocation method for a wide variety of two-dimensional grid sizes on both an IBM and a Cray 2 computer. We found that which algorithm is fastest on a particular machine depends not only on the grid size, but also on small details of the computer hardware as well. For most practical grid sizes used in computation, the even-odd decomposition algorithm is found to be faster than the transform-recursion method.

Don, Wai-Sun↗

NASA NTRS2003

Method for Balancing Detector Output to a Desired Level of Balance at a Frequency

A multi-gas sensor is provided which modulates a polarized light beam over a broadband of wavelengths between two alternating orthogonal polarization components. The two orthogonal polarization components of the polarization modulated beam are directed along two distinct optical paths. At least one optical path contains one or more spectral discrimination elements, with each spectral discrimination element having spectral absorption features of one or more gases of interest being measured. The two optical paths then intersect, and one orthogonal component of the intersected components is transmitted and the other orthogonal component is reflected. The combined polarization modulated beam is partitioned into one or more smaller spectral regions of interest where one or more gases of interest has an absorption band. The difference in intensity between the two orthogonal polarization components is then determined in each partitioned spectral region of interest as an indication of the spectral emission/absorption of the light beam by the gases of interest in the measurement path. The spectral emission/absorption is indicative of the concentration of the one or more gases of interest in the measurement path. More specifically, one embodiment of the present invention is a gas filter correlation radiometer which comprises a polarizer, a polarization modulator, a polarization beam splitter, a beam combiner, wavelength partitioning element, and detection element. The gases of interest are measured simultaneously and, further, can be measured independently or non-independently. Furthermore, optical or electronic element are provided to balance optical intensities between the two optical paths.

Sachse, Glenn W.↗

DOE OSTI2016

Brownian-like kinematics of ball milling for particulate structural modeling

Ball milling motion has been previously studied through computationally expensive, off-line experimental video processing and numerical simulations by the discrete element method. Here, this research establishes a more efficient formulation of the ball energetics and kinetics similar to the Brownian kinetic theory of statistical mechanics. Based on assumptions of thermomechanical equilibrium, negligible gravitational, aerodynamic and surface condition effects, and decoupled impact interaction among balls and with milled particulates, this model proposes mono-parametric spectral energy and velocity probability density functions akin to Maxwell-Boltzmann statistics, along with uniformly distributed impact directionality. The model predictions are calibrated and validated by comparison with published experimental measurements and computationally derived spectra. This descriptive Brownian-like motion model enables effective simulation of contact and impact, material deformation and micro-joining of ball milled bimetallic powders. A comprehensive simulation of the evolving internal fractal microstructure of the processed particulates is implemented at real-time computation speed, and its predictions are compared with experimental micrographs of ball milled Ni—Al particulates.

42 ENGINEERING↗

DOE OSTI2019

Real time time-dependent density functional theory using higher order finite-element methods

We present a computationally efficient approach to solve the time-dependent Kohn-Sham equations in real time using higher order finite-element spatial discretization, applicable to both pseudopotential and all-electron calculations. To this end, we develop an a priori mesh adaption technique, based on the semidiscrete (discrete in space but continuous in time) error estimate on the time-dependent Kohn-Sham orbitals, to construct an efficient finite-element discretization. Subsequently, we obtain the full-discrete error estimate to guide our choice of the time step. We employ spectral finite elements along with special reduced order quadrature to render the overlap matrix diagonal, thereby simplifying the inversion of the overlap matrix that features in the evaluation of the discrete time-evolution operator. We use the second-order Magnus operator as the time-evolution operator, wherein the action of the discrete Magnus operator, expressed as exponential of a matrix, on the Kohn-Sham orbitals is obtained efficiently through an adaptive Lanczos iteration. We observe close to optimal rates of convergence of the dipole moment with respect to spatial and temporal discretization, for both pseudopotential and all-electron calculations. We demonstrate a staggering 100-fold reduction in the computational time afforded by higher order finite elements over linear finite elements, for both pseudopotential and all-electron calculations. Further, for similar level of accuracy, we obtain significant computational savings by our approach as compared to state-of-the-art finite-difference methods. Here, we also demonstrate the competence of higher order finite elements for all-electron benchmark systems. Lastly, we observe good parallel scalability of the proposed method on many hundreds of processors.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Low-density laboratory spectra near the He ii λ 304 line

Aims. Our aim is to interpret the EUV spectra of the solar corona, one hopes for laboratory data of specific chemical elements obtained under coronal conditions. Methods. EUV spectra of He, C, N, O, F, Ne, S, Ar, Fe, and Ni in a 40 Å wide wavelength interval near λ304 were excited in an electron beam ion trap. Results. We observe some two hundred lines about half of which are not yet identified and included in spectral models. Conclusions. Lastly, our data provide a check on the atomic data bases underlying the spectral models that are used to interpret solar corona data. However, a multitude of mostly weak additional lines taken together represent a flux that is comparable to that of various primary lines.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2019

Assessing the Impurities in Neptunium Metal

A mass scan was previously run in the Q-ICP-MS for the Np metal dissolved samples. Peaks for Ti, Nb, I, Au, Tl, Na, Al, Sn, Ni, Zn, Y, Cr, Pb, K, Fe, Ca, Ta, Th, V, U and Pu were observed. Those elements were isotopically quantified using the HR-ICP-MS. For these determinations, the elements were split in 7 different sets, accordingly to its concentration in the sample and/or possible interferences. For each set of elements, Internal Standard (IS) was added in a similar concentration for the sample. Samples and standards were run both in 300 resolution power (RP) and 4000 RP, in order to evaluate possible interferences (except for actinides, that were run only in 300 RP). It was possible to identify the main sources of impurities, with great precision and accuracy. As the HR-ICP-MS is considered a recent instrument, specially when applied for the analysis of nuclear materials, the method development can be published in peer-reviewed journals. Np metal produced did not achieve the desired purity (which was 99.758% for the elements analyzed). The main impurities were Ta, Fe, Th, U and Pu. This information will help modify the casting process to achieve a higher purity. The maximum achievable purity with our current method is 99.990%. This is being constantly improved, now that we have enough information about the interferences on these elements analyses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Ultrasonic estimation of the contact angle of a sessile droplet

Radiation of energy by large amplitude leaky Rayleigh waves is regarded as one of the key physical mechanisms regulating the actuation and manipulation of droplets in surface acoustic wave (SAW) microfluidic devices. The interaction between a SAW and a droplet is highly complex and is presently the subject of extensive research. This paper investigates the existence of an additional interaction mechanism based on the propagation of quasi-Stoneley waves inside sessile droplets deposited on a solid substrate. In contrast with the leaky Rayleigh wave, the energy of the Stoneley wave is confined within a thin fluid layer in contact with the substrate. The hypothesis is confirmed by three-dimensional finite element simulations and ultrasonic scattering experiments measuring the reflection of Rayleigh waves from droplets of different diameters. Moreover, real-time monitoring of the droplet evaporation process reveals a clear correlation between the droplet contact angle and the spectral information of the reflected Rayleigh signal, thus paving the way for ultrasonic measurements of surface tension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Spectral Peak Enhancement by Combining Trusted Response Elements via Machine Learning (SPECTRE-ML) v0.8.0

SPECTRE-ML (Spectral Peak Enhancement by Combining Trusted Response Elements via Machine Learning) is a machine learning based program for finding optimal clusters of radiation detector segments (i.e., pixels or voxels) in order to improve spectral performance. It provides facilities for pre-processing and analyzing training datasets, running ML algorithms, and evaluating and visualizing outputs. SPECTRE-ML outperforms simpler ad-hoc segmentation methods such as uniform depth clusters, learning detector performance trends such as dead layers, edge effects, and gain shifts. Although extensible to arbitrary highly-segmented spectroscopic radiation detectors, SPECTRE-ML currently focuses on improving spectral performance in highly-segmented CdZnTe (CZT) detectors for International Atomic Energy Agency (IAEA) non-destructive assay (NDA) safeguards tasks.

Vavrek, Jayson↗

DOE OSTI2021

Spectral signatures of transverse optical modes in semiconductor nanowires

In order to build nanophotonic devices, it is important to understand and ultimately control the optical mode structure within potential components such as nanoscale waveguides. However, experimental characterization of such modes in the optical regime is difficult due to the nanoscale dimensions of such components and the perturbations that would be induced by a near-field probe. Here, we demonstrate a probe-free, far-field method to characterize the optical modes within GaN nanowires (NWs) based on a novel off-axis scanning confocal microscope system. Using this microscope, we observe spectral signatures resulting from lateral leakage of waveguide modes when they exceed their respective cutoff limits. We identify these modes within hyperspectral images using an analytical model coupled with finite element simulations. The model can also be used to predict the spectral signatures for given geometrical and dielectric parameters, which enabled us to deduce the transverse dimension of the NW from hyperspectral images with an accuracy of ∼ <#comment/> 30 nm .

Wang, Maoji↗

DOE OSTI2017

Cryo-STEM Tomography with Inpainting

Traditionally, microscopists have worked with the Nyquist-Shannon theory of sampling, which states that to be able to reconstruct the image fully it needs to be sampled at a rate of at least twice the highest spatial frequency in the image. This sampling rate assumes that the image is sampled at regular intervals and that every pixel contains information that is crucial for the image (it even assumes that noise is important). Images in general, and especially low dose S/TEM images, contain significantly less information than can be encoded by a grid of pixels (which is why image compression works). Mathematically speaking, the image data has a low dimensional or sparse representation. Through the application of compressive sensing methods [1,2,3] this representation can be found using pre-designed measurements that are usually random for implementation simplicity. These measurements and the compressive sensing reconstruction algorithms have the added benefit of reducing noise. This reconstruction approach can be extended into higher dimensions, whereby the random sampling in each 2-D image can be extended into: a sequence of tomographic projections (i.e. tilt images); a sequence of video frames (i.e. incorporating temporal resolution and dynamics); spectral resolution (i.e. energy filtering an image to see the distribution of elements); and ptychography (i.e. sampling a full diffraction image at each location in a 2-D grid across the sample). This approach has been employed experimentally for materials science samples requiring low-dose imaging [2], and can be readily applied to biological samples. Figure 1 shows the resolution possible in a complex biological system, mouse pancreatic islet beta cells [4], when tomogram slices are reconstructed using subsampling. Reducing the number of pixels (1/6 pix and 1/3*1/3) shows minimal degradation compared to the reconstructions using all pixels (all data and 1/3 tilt). Although subsampling 1/6 of the tilts (1/6 of overall dose) degrades the reconstruction to the point that the cellular structures cannot be identified. Using 1/3 of both the pixels and the tilts provides a high quality image at 1/9 the overall dose even for this most basic and rapid demonstration of the CS methods. Figure 2 demonstrates the theoretical tomogram reconstruction quality (vertical axis) as undersampling (horizontal axis) is increased; we examined subsampling pixels and tilt-angles individually and a combined approach in which both pixels and tilts are sub-sampled. Note that subsampling pixels maintains high quality reconstructions (solid lines). Using the inpainting algorithm to obtain tomograms can automatically reduce the dose applied to the system by an order of magnitude. Perhaps the best way to understand the impact is to consider that by using inpainting (and with minimal hardware changes), a sample that can normally withstand a dose of ~10 e/Å2 can potentially be imaged with an “equivalent quality” to a dose level of 103 e/Å2. To put this in perspective, this is approaching the dose level used for the most advanced images, in terms of spatial resolution, for inorganic systems. While there are issues for biological specimens beyond dose (structural complexity being the most important one), this sampling approach allows the methods that are traditionally used for materials science to be applied to biological systems [5]. References: [1] A Stevens, H Yang, L Carin et al. Microscopy 63(1), (2014), pp. 41. [2] L Kovarik, A Stevens, A Liyu et al. Appl. Phys. Lett. 109, 164102 (2016) [3] A Stevens, L Kovarik, P Abellan et al. Adv. Structural and Chemical Imaging 1(10), (2015), pp. 1. [4] MD Guay, W Czaja, MA Aronova et al. Scientific Reports 6, 27614 (2016) [5] Supported by the Chemical Imaging, Signature Discovery, and Analytics in Motion Initiatives at PNNL. PNNL is operated by Battelle Memorial Inst. for the US DOE; contract DE-AC05-76RL01830.

Stevens, Andrew↗

DOE OSTI2020

High-temperature Majorana fermions in magnet-superconductor hybrid systems

Magnet-superconductor hybrid (MSH) structures represent one of the most promising platforms to realize, control, and manipulate Majorana modes using scanning tunneling methods. By depositing either chains or islands of magnetic atoms on the surface of a conventional, elemental superconductor such as Pb or Re, topological superconducting phases can emerge. They feature either localized Majorana bound states at the chain ends or dispersing chiral Majorana modes at the island's boundary. However, some of these experiments have not reached the spectral resolution to clearly distinguish between topological Majorana and trivial Shiba states due to very small superconducting gap sizes and experiments performed at sub-Kelvin temperatures. Here we consider superconducting substrates with unconventional spin-singlet pairing, including high-temperature d -wave and extended s -wave superconductors. We derive topological phase diagrams and compute edge states for cylinder and island geometries and discuss their properties. Several time-reversal invariant topological superconducting phases of the Zhang-Kane-Mele type are found and discussed. Addiotnally, we review one-dimensional MSH structures and show that parameters to realize topologically nontrivial magnetic chains embedded into a larger, two-dimensional substrate differ from the purely one-dimensional case. Quite generally we find that unconventional superconducting substrates work as well as the conventional s -wave substrates to realize topological phases. In particular, iron-based pnictide and chalcogenide superconductors are the most promising class of substrates for future high-temperature MSH systems.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Inter-digital capacitive sensor for evaluating cable jacket and insulation aging

An inter-digital capacitive (TDC) sensor has previously been used to measure dielectric properties of cable insulation polymer material when placed in direct contact with the insulation. Often cable insulation is covered by a polymer jacket. The dielectric properties of many cable jacket and insulation polymers are known to change due to thermal and radiation exposure-related damage. These dielectric properties frequently track with other measures of cable aging, such as tensile elongation-at-break and indenter modulus that have been broadly established as cable insulation polymer assessment methods. The external jacket of a cable is likely to have a different permittivity from the underlying insulation, and frequently the jacket material exhibits more severe damage than the insulation material due to environmental exposure. Because the jacket serves primarily to guard the cable during installation, as long as the underlying insulation condition is acceptable, the jacket condition is relatively unimportant in service. As part of a continuing program to develop and evaluate nondestructive examination methods that may be applied to cable condition assessment, a set of tools has been developed including (1) a parallel-plate sensor to directly measure the permittivity spectrum of flat sheet material and (2) an TDC and fixture to measure the effect of cable polymer dielectric property change on the sensor response. The TDC consists of two fork-like electrodes facing each other with the fork tines interspersed and separated by a small gap. The electrodes are printed on one side of a flexible substrate that can be conformed to the surface of a cylindrical cable, with tines parallel to the cable axis. The electrodes are connected to a broad-frequency-spectral impedance meter that senses the capacitance between the narrowly gapped electrode tines. This capacitance is known to vary as a function of the permittivity of any material in close proximity to the electrodes. By finite element modeling (FEM) and experimentation, this study investigates the effect of tine spacing and other design parameters associated with the TDC on the voltage (potential) distribution and electric field depth of penetration. The TDC measurement of an unshielded ethylene-propylene rubber (EPR)-insulated cable is shown to track with the degree of aging and quantities obtained by established methods. For jacketed cable systems, the TDC response is dominated by the jacket but, by analyzing measurements from TDC sensors with different depths-of-field penetration into the cable under test, the influence of the chlorinated polyethylene (CPE) cable jacket material degradation can be separated from an assessment of the cable insulation thereby enabling assessment of the insulation beneath/through the jacket.

Cable NDE, Interdigital Capacitance, Cable test↗

DOE OSTI2025

Demonstration of an AI-driven workflow for dynamic x-ray spectroscopy

X-ray absorption near edge structure (XANES) spectroscopy is a powerful technique for characterizing the chemical state and symmetry of individual elements within materials, but requires collecting data at many energy points which can be time-consuming. While adaptive sampling methods exist for efficiently collecting spectroscopic data, they often lack domain-specific knowledge about the structure of XANES spectra. Here we demonstrate a knowledge-injected Bayesian optimization approach for adaptive XANES data collection that incorporates understanding of spectral features like absorption edges and pre-edge peaks. We show this method accurately reconstructs the absorption edge of XANES spectra using only 15–20% of the measurement points typically needed for conventional sampling, while maintaining the ability to determine the x-ray energy of the sharp peak after the absorption edge with errors less than 0.03 eV, the absorption edge with errors less than 0.1 eV; and overall root-mean-square errors less than 0.005 compared to traditionally sampled spectra. Our experiments on battery materials and catalysts demonstrate the method’s effectiveness for both static and dynamic XANES measurements, improving data collection efficiency and enabling better time resolution for tracking chemical changes. This approach advances the degree of automation in XANES experiments, reducing the common errors of under- or over-sampling points near the absorption edge and enabling dynamic experiments that require high temporal resolution or limited measurement time.

Bayesian optimization↗

DOE OSTI2015

Determination of a brass alloy concentration composition using calibration-free laser-induced breakdown spectroscopy

Laser-induced breakdown spectroscopy (LIBS) is a technique that can provide qualitative and quantitative measurements of the characteristics of irradiated metals. In the present work, we have calculated the parameters of the plasma produced from a brass alloy sample under the action of a pulsed Nd: YAG laser operating at 1064 nm. The emission lines of copper atoms (Cu I), zinc atoms (Zn I), and lead atoms (Pb I), which are elements of a brass alloy composition, were used to investigate the parameters of the brass plasma. The spectral profiles of Cu, Zn, and Pb lines have been used to extract the electron temperature and density of the brass alloy plasma. The characteristics of Cu, Zn, and Pb were determined quantatively by the calibration-free LIBS (CF-LIBS) method considering for accurate analysis that the laser-induced ablated plasma is optically thin in local thermodynamic equilibrium conditions and the plasma ablation is stoichiometric. The Boltzmann plot method was used to evaluate the plasma temperature, and the Stark broadened profiles were used to determine the electron density. An algorithm based on the experimentally measured values of the intensity of spectral lines and the basic laws of plasma physics was developed for the determination of Cu, Zn, and Pb concentrations in the brass sample. The concentrations C{sub CF-LIBS} calculated by CF-LIBS and the certified concentrations C{sub certified} were very close.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Determining Stellar Elemental Abundances from DESI Spectra with the Data-driven Payne

Abstract Stellar abundances for a large number of stars provide key information for the study of Galactic formation history. Large spectroscopic surveys such as the Dark Energy Spectroscopic Instrument (DESI) and LAMOST take median-to-low-resolution (R≲ 5000) spectra in the full optical wavelength range for millions of stars. However, the line-blending effect in these spectra causes great challenges for elemental abundance determination. Here we employDD-Payne, a data-driven method regularized by differential spectra from stellar physical models, to the DESI early data release spectra for stellar abundance determination. Our implementation delivers 15 labels, including effective temperatureT eff , surface gravity log g , microturbulence velocityv mic , and the abundances for 12 individual elements, namely C, N, O, Mg, Al, Si, Ca, Ti, Cr, Mn, Fe, and Ni. Given a spectral signal-to-noise ratio of 100 per pixel, the internal precisions of the label estimates are about 20 K forT eff , 0.05 dex for log g , and 0.05 dex for most elemental abundances. These results agree with the theoretical limits from the Crámer–Rao bound calculation within a factor of 2. The majority of the accreted halo stars contributed by the Gaia–Enceladus–Sausage are discernible from the disk and in situ halo populations in the resultant [Mg/Fe]–[Fe/H] and [Al/Fe]–[Fe/H] abundance spaces. We also provide distance and orbital parameters for the sample stars, which spread over a distance out to ∼100 kpc. The DESI sample has a significantly higher fraction of distant (or metal-poor) stars than the other existing spectroscopic surveys, making it a powerful data set for studying the Galactic outskirts. The catalog is publicly available.

Astronomy & Astrophysics↗

NASA NTRS1985

Coupled atmosphere/canopy model for remote sensing of plant reflectance features

Solar radiative transfer through a coupled system of atmosphere and plant canopy is modeled as a multiple-scattering problem through a layered medium of random scatterers. The radiative transfer equation is solved by the discrete-ordinates finite-element method. Analytic expressions are derived that allow the calculation of scattering and absorption cross sections for any plant canopy layer form measurable biophysical parameters such as the leaf area index, leaf angle distribution, and individual leaf reflectance and transmittance data. An expression for a canopy scattering phase function is also given. Computational results are in good agreement with spectral reflectance measurements directly above a soybean canopy, and the concept of greenness- and brightness-transforms of Landsat MSS data is reconfirmed with the computed results. A sensitivity analysis with the coupled atmosphere/canopy model quantifies how satellite-sensed spectral radiances are affected by increased atmospheric aerosols, by varying leaf area index, by anisotropic leaf scattering, and by non-Lambertian soil boundary conditions. Possible extensions to a 2-D model are also discussed.

Gerstl, S. A.↗