Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “least-squares analysis”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A function space approach to state and model error estimation for elliptic systems

An approach is advanced for the concurrent estimation of the state and of the model errors of a system described by elliptic equations. The estimates are obtained by a deterministic least-squares approach that seeks to minimize a quadratic functional of the model errors, or equivalently, to find the vector of smallest norm subject to linear constraints in a suitably defined function space. The minimum norm solution can be obtained by solving either a Fredholm integral equation of the second kind for the case with continuously distributed data or a related matrix equation for the problem with discretely located measurements. Solution of either one of these equations is obtained in a batch-processing mode in which all of the data is processed simultaneously or, in certain restricted geometries, in a spatially scanning mode in which the data is processed recursively. After the methods for computation of the optimal esimates are developed, an analysis of the second-order statistics of the estimates and of the corresponding estimation error is conducted. Based on this analysis, explicit expressions for the mean-square estimation error associated with both the state and model error estimates are then developed. While this paper focuses on theoretical developments, applications arising in the area of large structure static shape determination are contained in a closely related paper (Rodriguez and Scheid, 1982).

Rodriguez, G.↗
NASA NTRS1990

DFW microburst model based on AA-539 data

Analysis of the August 2, 1985 crash for an L-1011 jumbo jet (DL-191) on approach to the Dallas-Ft. Worth International Airport (DFW) in a thunderstorm indicates that the severe windshear microburst that caused the crash was composed not only of a strong downflow and outflow but also included several large-scale vortex rings entrained in the flowfield. This paper presents a detailed two-dimensional model of the DFW microburst based on data from the MD-80 (AA-539) that followed behind DL-191 and flew through the microburst about two minutes after the crash of DL-191. The model was developed using wind-vector and flight-path data reconstructed by NASA Ames Research Center and a combination of interactive graphics and least-squares error best fit between the modeled and measured wind vectors along the AA-539 flight path. The model indicates that the flowfield contains some significant elements and vortices not previously reported. The alternating direction of rotation of the vortices in the model suggests a microburst structure based on a von Karman vortex street rather than on a Kelvin-Helmholtz instability. The model also indicates that the reconstructed wind-vector data contain a time lag of at least one second in the horizontal winds.

Grantham, Walter J.↗
NASA NTRS1995

Calibration of the advanced microwave sounding unit-A for NOAA-K

The thermal-vacuum chamber calibration data from the Advanced Microwave Sounding Unit-A (AMSU-A) for NOAA-K, which will be launched in 1996, were analyzed to evaluate the instrument performance, including calibration accuracy, nonlinearity, and temperature sensitivity. The AMSU-A on NOAA-K consists of AMSU-A2 Protoflight Model and AMSU-A1 Flight Model 1. The results show that both models meet the instrument specifications, except the AMSU-A1 antenna beamwidths, which exceed the requirement of 3.3 +/- 10%. We also studied the instrument's radiometric characterizations which will be incorporated into the operational calibration algorithm for processing the in-orbit AMSU-A data from space. Particularly, the nonlinearity parameters which will be used for correcting the nonlinear contributions from an imperfect square-law detector were determined from this data analysis. It was found that the calibration accuracies (differences between the measured scene radiances and those calculated from a linear two-point calibration formula) are polarization-dependent. Channels with vertical polarizations show little cold biases at the lowest scene target temperature 84K, while those with horizontal polarizations all have appreciable cold biases, which can be up to 0.6K. It is unknown where these polarization-dependent cold biases originate, but it is suspected that some chamber contamination of hot radiances leaked into the cold scene target area. Further investigation in this matter is required. The existence and magnitude of nonlinearity in each channel were established and a quadratic formula for modeling these nonlinear contributions was developed. The model was characterized by a single parameter u, values of which were obtained for each channel via least-squares fit to the data. Using the best-fit u values, we performed a series of simulations of the quadratic corrections which would be expected from the space data after the launch of AMSU-A on NOAA-K. In these simulations, the cosmic background radiance corresponding to a cold space temperature 2.73K was adopted as one of the two reference points of calibration. The largest simulated nonlinear correction is about 0.3K, which occurs at channel 15 when the instrument temperature is at 38.09 deg C. Others are less than 0.2K in the remaining channels. Possible improvement for future instrument calibration is also discussed.

Mo, Tsan↗
NASA NTRS2001

Two-Component Fitting of Coronal-Hole and Quiet-Sun He I 1083 Spectra

We present reduction techniques and first results for detailed fitting of solar spectra obtained with the NASA/National Solar Observatory Spectromagnetograph (NASA/NSO SPM over a 2 nm bandpass centered on the He 1 1083 nm line. The observation for this analysis was a spectra-spectroheliogram obtained at the NSO/Kitt Peak Vacuum Telescope (KPVT) on 00 Apr 17 at 21:46 UT spanning an area of 512 x 900 arc-seconds; the field of view included a coronal hole near disk center as well as surrounding quiet sun. Since the He I line is very weak and blended with nearby solar and telluric lines, accurate determination of the continuum intensity as a function of wavelength is crucial. We have modified the technique of Malanushenko {\it et al.) (1992; {\it AA) (\bf 259), 567) to tie regions of continuua and the wings of spectral lines which show little variation over the image to standard reference spectra such as the NSO Fourier Transform Spectrometer atlas (Wallace {\it et al). 1993; NSO Tech Report \#93-001). We performed detailed least-squares fits of spectra from selected areas, accounting for all the known telluric and solar absorbers in the spectral bandpass. The best physically consistent fits to the Helium lines were obtained with Gaussian profiles from two components (one ''cool'', characteristic of the upper chromosphere; one ''hot'', representing the cool transition region at 2-3 x 10$^{4)$ K). In the coronal hole, the transition-region component, shifted by 6-7 km/s to the blue, is mildly dominant, consistent with mass outflow as suggested by Dupree {\it et all. (1996; {\it Ap. J.}-{\bf 467), 121). In quiet-sun spectra there is less evidence of outward flow, and the chromospheric component is more important. All our fitted spectra show a very weak unidentified absorption feature at 1082.880 nm in the red wing of the nearby Si I line.

Jones, Harrison P.↗
NASA NTRS1989

The relationship between IR, optical, and UV extinction

An analysis is presented for the variability of absolute IR, optical, and UV extinction, A(sub lambda), derived through the ratio of total-to-selective extinction, R, for 31 lines of sight for which reliable UV extinction parameters were derived. These data sample a wide range of environments and are characterized by 2.5 is less than or equal to R is less than or equal to 6.0. It was found that there is a strong linear dependence between extinction expressed as A(sub lambda)/A(sub V) and 1/R for 1.25 micron is less than or equal to lambda is less than or equal to 0.12 micron. Differences in the general shape of extinction curves are largely due to variations in shape of optical/near-UV extinction corresponding to changes in R, with A(sub lambda)/A(sub V) decreasing for increasing R. From a least-squares fit of the observed R-dependence as a function of wavelength for 0.8/micron is less than or greater than 1/lambda is less than or equal to 8.3/micron, an analytic expression was generated from which IR, optical, and UV extinction curves of the form A(sub lambda)/A(sub V) can be reproduced with reasonable accuracy from a knowledge of R. It was also found that the absolute bump strength normalized to A(sub V) shows a general decrease with increasing R, suggesting that some fraction of bump grains may be selectively incorporated into coagulated grains. Finally, it was found that absolute extinction normalized by suitably chosen color indices results in a minimization of the R-dependence of portions of the UV curve, allowing A(sub lambda) to be estimated for these wavelengths independent of R.

Cardelli, Jason A.↗
NASA NTRS1992

Extracting galactic structure parameters from multivariated density estimation

Multivariate statistical analysis, including includes cluster analysis (unsupervised classification), discriminant analysis (supervised classification) and principle component analysis (dimensionlity reduction method), and nonparameter density estimation have been successfully used to search for meaningful associations in the 5-dimensional space of observables between observed points and the sets of simulated points generated from a synthetic approach of galaxy modelling. These methodologies can be applied as the new tools to obtain information about hidden structure otherwise unrecognizable, and place important constraints on the space distribution of various stellar populations in the Milky Way. In this paper, we concentrate on illustrating how to use nonparameter density estimation to substitute for the true densities in both of the simulating sample and real sample in the five-dimensional space. In order to fit model predicted densities to reality, we derive a set of equations which include n lines (where n is the total number of observed points) and m (where m: the numbers of predefined groups) unknown parameters. A least-square estimation will allow us to determine the density law of different groups and components in the Galaxy. The output from our software, which can be used in many research fields, will also give out the systematic error between the model and the observation by a Bayes rule.

Chen, B.↗
NASA NTRS1989

An airborne study of microwave surface sensing and boundary layer heat and moisture fluxes for Project FIFE

This is the third semiannual status report describing the University of Wyoming participation in project FIFE, under NASA Grant NAG 5-913. The reader is referred to the first and second semiannual status reports for this grant for descriptions of University of Wyoming participation in the field operations of 1987. This report is concerned with continuing analysis of field data. Most of the instrument calibration effort has concentrated on interpretation and use of the high-rate water vapor measurements with the University of Wyoming Lyman-alpha device. In addition, fluxes of water vapor have been corrected for sensible heat flux, and all the high-rate (10 Hz for the Wyoming King Air) data have been filtered before the flux calculations. For each flight the Lyman-alpha data (10 Hz) are calibrated against the Cambridge dew-point hygrometer (1 Hz), by using data from the take-off sounding. The reasoning here is that the takeoff sounding (held to 500 fpm rate-of-climb specifically for the calibration) covers the widest range of humidity conditions of any sequence of the flight. For each takeoff sounding, then, the Lyman-alpha data, averaged to 1 Hz, are matched to the dewpoint data with a second-order least-squares fit. The resulting quadratic equation is then used in all subsequent analyses to convert the Lyman-alpha voltages to vapor desities for that specific flight. An alternative method is to use the range of fluctuation for the same two devices from horizontal passes to establish the calibration function. However, since the Lyman-alpha data are used only in the flux calculations, i.e., only the fluctuation and not the absolute values are used, this method is thought unnecessary.

Kelly, Robert D.↗
NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗
NASA NTRS1999

A Post-Pathfinder Evaluation of Areocentric Solar Coordinates with Improved Timing Recipes for Mars Seasonal/Diurnal Climate Studies

The accurate determination of the Mars pole vector derived from Pathfinder and Viking Lander radio data, together with the VSOP87 representation of planetary orbits, have been applied to a new evaluation of the right ascension of the "fictitious mean sun" (FMS) at Mars. With DELTA t (sub J2000) the elapsed time in days from the J2000 epoch (J.D.2451545.0 (sup TT), alpha FMS = 270 degrees.3863 + 0.52403840(degrees/d) (raised dot) DELTA T (sub j2000) - 4 x 10 (exp -13) (degrees/d (sup 2)) (raised dot) DELTA t (sup 2) (sub J2000) represents a best least-squares quadratic fit of the FMS, including aberration, to each instance of the four equinox and solstice passages for each of 134 Mars orbits spanning the calendar years 1874-2127. The implied tropical orbit period for Mars, 686.9726 (sup d), closely agrees with the recent evaluations. Together with the Pathfinder radio determination of the Mars sidereal rotation, the derived FMS rate corresponds to a mean solar day (or "sol") of 1.027491251 (sup d). The new FMS determination would serve to define the Mean Solar Time at Mars to the nearest tenth-second, according to historical conventions originally established for terrestrial time keeping, once the Mars prime meridian defined by the crater Airy-O is navigated to the same accuracy. For convenient reference to current epochs, 2000 Jan 06 00:00 UTC (= MJD 51549.000 (sup UTC)) corresponds to a coincidence of (alpha (sub FMS)) and the rotation angle of the crater Airy-O measured with respect to the Mars equinox (i.e. "mean solar midnight" on the planet's prime meridian), to within the current uncertainty of several seconds in the locational definition of the planet's cartographic grid. As a further result of the analysis, the consistently derived Mars obliquity of date is epsilon = 25 degrees.192 + 3.45 x l0 (exp -7)(degrees/d)(raised dot) DELTA t (sub J2000). An improved analytic recipe for the calculation of the solar areocentric longitude (L (sub s)) of Mars to an accuracy of 0 degrees.01 is also provided, accounting for the primary perturbations of Earth, Jupiter, and Venus, which may in turn be applied to an efficient evaluation of Mars local true solar time (LTST) to within the uncertainty of the inertial position of the Mars prime meridian. For specific applications to the data archives for landed Mars spacecraft, simple conversion formulae are given for the determination of the Viking "Local Lander Time" and the Pathfinder "Local True Solar Time" in terms of the terrestrial calendar date and UTC.

Allison, Michael↗
DOE OSTI2025

Rapid monitoring of fermentations: a feasibility study on biological 2,3-butanediol production

2,3-butanediol (2,3-BDO) is an economically important platform chemical that can be produced by the fermentation of sugars using an engineered strain of Zymomonas mobilis . These fermentations require continuous monitoring and modification of fermentation conditions to maximize 2,3-BDO yields and minimize the production of the undesired coproducts glycerol and acetoin. Because of the time required for sampling and off-line chromatographic measurement of fermentation samples, the ability of fermentation scientists to modify fermentation conditions in a timely manner is limited. The goal of this study was to test if near-infrared spectroscopy (NIRS) along with multivariate statistics could reduce the time needed for this analysis and enable real-time monitoring and control of the fermentation. In this work we developed partial least squares (PLS) calibration models to predict the concentrations of glucose, xylose, 2,3-BDO, acetoin, and glycerol in fermentations via NIRS using two different spectrometers and two different spectroscopy modalities. We first evaluated the feasibility of rapid NIRS monitoring through experiments where we measured the signals from each analyte of interest and built NIRS-based PLS models using spectra from synthetic samples containing uncorrelated concentrations of these analytes. All analytes showed unique spectral signatures, and this initial modeling showed that all analytes could be detected simultaneously. We then began work with samples from laboratory fermentation experiments and tested the feasibility of regression model development across two spectral collection modalities (at-line and on-line) and two instruments: a laboratory-grade instrument and a low-cost instrument with a more limited spectral range. All modalities showed promise in the ability to monitor Z. mobilis fermentations of glucose and xylose to 2,3-BDO. The low-cost instrument displayed a lower signal-to-noise ratio than the laboratory-grade instrument, which led to comparatively lower performance overall, but still provided sufficient accuracy to monitor fermentation trends. While the ease of use of on-line monitoring systems was favored as compared to at-line systems due to the lack of sampling required and potential for automated process control, we observed some decrease in performance due to the additional complexity of the sample matrix. We have demonstrated that NIRS combined with multivariate analysis can be used for at-line and on-line monitoring of the concentrations of glucose, xylose, 2,3-BDO, acetoin, and glycerol during Z. mobilis fermentations. The decrease in signal-to-noise ratio when using a low-cost spectrometer led to greater prediction error than the laboratory-grade spectrometer for at-line monitoring. The on-line monitoring modality showed great promise for real time process control via NIRS.

09 BIOMASS FUELS↗
NASA NTRS1993

A finite-element-based perturbation model for the rotordynamic analysis of shrouded pump impellers: Part 2: User's guide

This report describes the computational steps involved in executing a finite-element-based perturbation model for computing the rotor dynamic coefficients of a shrouded pump impeller or a simple seal. These arise from the fluid/rotor interaction in the clearance gap. In addition to the sample cases, the computational procedure also applies to a separate category of problems referred to as the 'seal-like' category. The problem, in this case, concerns a shrouded impeller, with the exception that the secondary, or leakage, passage is totally isolated from the primary-flow passage. The difference between this and the pump problem is that the former is analytically of the simple 'seal-like' configuration, with two (inlet and exit) flow-permeable stations, while the latter constitutes a double-entry / double-discharge flow problem. In all cases, the problem is that of a rotor clearance gap. The problem here is that of a rotor excitation in the form of a cylindrical whirl around the housing centerline for a smooth annular seal. In its centered operation mode, the rotor is assumed to give rise to an axisymmetric flow field in the clearance gap. As a result, problems involving longitudinal or helical grooves, in the rotor or housing surfaces, go beyond the code capabilities. Discarding, for the moment, the pre- and post-processing phases, the bulk of the computational procedure consists of two main steps. The first is aimed at producing the axisymmetric 'zeroth-order' flow solution in the given flow domain. Detailed description of this problem, including the flow-governing equations, turbulence closure, boundary conditions, and the finite-element formulation, was covered by Baskharone and Hensel. The second main step is where the perturbation model is implemented, with the input being the centered-rotor 'zeroth-order' flow solution and a prescribed whirl frequency ratio (whirl frequency divided by the impeller speed). The computational domain, in the latter case, is treated as three dimensional, with the number of computational planes in the circumferential direction being specified a priori. The reader is reminded that the deformations in the finite elements are all infinitesimally small because the rotor eccentricity itself is a virtual displacement. This explains why we have generically termed the perturbation model the 'virtually' deformable finite-element category. The primary outcome of implementing the perturbation model is the tangential and radial components, F(sub theta)(sup *) and F(sub r)(sup *) of the fluid-exerted force on the rotor surface due to the whirling motion. Repetitive execution of the perturbation model subprogram over a sufficient range of whirl frequency ratios, and subsequent interpolation of these fluid forces, using the least-square method, finally enable the user to compute the impeller rotor dynamic coefficients of the fluid/rotor interaction. These are the direct and cross-coupled stiffness, damping, and inertia effects of the fluid/rotor interaction.

Baskharone, Erian A.↗
DOE OSTI2025

Tracking seasonal variability in plant traits from spaceborne PRISMA and NEON AOP across forest types and ecoregions

Plant traits serve as critical indicators of how plants adapt to environmental changes and influence ecosystem functions. While airborne hyperspectral remote sensing effectively maps plant traits through detailed reflectance properties, it is limited by cost and scale, making large-scale and temporal studies challenging. The recently launched spaceborne hyperspectral imager, PRecursore IperSpettrale della Missione Applicativa (PRISMA), offers frequent, large scale and high-fidelity observations on a spatial resolution of 30 m and a revisit time of around 29 days, making it suitable for large-scale seasonal trait mapping. However, their potential remains largely unexplored. This study developed a multi-stage framework by leveraging the PRISMA spaceborne hyperspectral data and National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) hyperspectral data to investigate the seasonal dynamics of four key plant traits — chlorophyll content, carotenoid content, equivalent water thickness, and nitrogen content — across eleven NEON sites representing diverse forest types and ecoregions in the contiguous U.S. Our results demonstrated that PRISMA hyperspectral data can reliably track seasonal variability in plant traits, achieving overall R 2 values ranging from 0.78 to 0.88 and normalized root mean square error (NRMSE) values ranging from 5.4% to 8.4% for the four traits. Seasonal patterns revealed bell-shaped trajectories for chlorophyll and carotenoids, while equivalent water thickness decreased steadily across most sites, driven by structural changes during leaf maturation and senescence. Nitrogen content exhibited less pronounced seasonal variation but followed expected nutrient resorption patterns. Analysis of environmental drivers showed that seasonal variability is primarily controlled by solar radiation and day length in northern sites, vapor pressure in semi-arid regions, and temperature in mid-southeastern sites. Spatial variability, meanwhile, was primarily driven by soil properties, particularly during the peak growing season. However, the influence of soil variables slightly declines toward the end of the season at several sites, as climatic factors become more prominent. This study highlights the capability of PRISMA, and potentially other similar spaceborne hyperspectral data for large-scale, time-series plant trait mapping and provides valuable insights into the interactions between plant traits and environmental factors. In conclusion, these findings contribute to advancing our understanding of plant functional ecology and improving predictions of ecosystem responses to environmental changes.

Environmental drivers↗
NASA NTRS1982

Least-squares sequential parameter and state estimation for large space structures

This paper presents the formulation of simultaneous state and parameter estimation problems for flexible structures in terms of least-squares minimization problems. The approach combines an on-line order determination algorithm, with least-squares algorithms for finding estimates of modal approximation functions, modal amplitudes, and modal parameters. The approach combines previous results on separable nonlinear least squares estimation with a regression analysis formulation of the state estimation problem. The technique makes use of sequential Householder transformations. This allows for sequential accumulation of matrices required during the identification process. The technique is used to identify the modal prameters of a flexible beam.

Thau, F. E.↗
NASA NTRS1989

Optimal focal-plane restoration

Image restoration can be implemented efficiently by calculating the convolution of the digital image and a small kernel during image acquisition. Processing the image in the focal-plane in this way requires less computation than traditional Fourier-transform-based techniques such as the Wiener filter and constrained least-squares filter. Here, the values of the convolution kernel that yield the restoration with minimum expected mean-square error are determined using a frequency analysis of the end-to-end imaging system. This development accounts for constraints on the size and shape of the spatial kernel and all the components of the imaging system. Simulation results indicate the technique is effective and efficient.

Reichenbach, Stephen E.↗
NASA NTRS1991

Optimal least-squares finite element method for elliptic problems

An optimal least squares finite element method is proposed for two dimensional and three dimensional elliptic problems and its advantages are discussed over the mixed Galerkin method and the usual least squares finite element method. In the usual least squares finite element method, the second order equation (-Delta x (Delta u) + u = f) is recast as a first order system (-Delta x p + u = f, Delta u - p = 0). The error analysis and numerical experiment show that, in this usual least squares finite element method, the rate of convergence for flux p is one order lower than optimal. In order to get an optimal least squares method, the irrotationality Delta x p = 0 should be included in the first order system.

Jiang, Bo-Nan↗
NASA NTRS1993

Active control of wake/blade-row interaction noise through the use of blade surface actuators

A combined analytical/computational approach for controlling of the noise generated by wake/blade-row interaction through the use of anti-sound actuators on the blade surfaces is described. A representative two-dimensional section of a fan stage, composed of an upstream fan rotor and a downstream fan exit guide vane (FEGV), is examined. An existing model for the wakes generated by the rotor is analyzed to provide realistic magnitudes for the vortical excitations imposed at the inlet to the FEGV. The acoustic response of the FEGV is determined at multiples of the blade passing frequency (BPF) by using the linearized unsteady flow analysis, LINFLO. Acoustic field contours are presented at each multiple of BPF illustrating the generated acoustic response disturbances. Anti-sound is then provided by placing oscillating control surfaces, whose lengths and locations are specified arbitrarily, on the blades. An analysis is then conducted to determine the complex amplitudes required for the control surface motions to best reduce the noise. It is demonstrated that if the number of acoustic response modes to be controlled is equal to the number of available independent control surfaces, complete noise cancellation can be achieved. A weighted least squares minimization procedure for the control equations is given for cases in which the number of acoustic modes exceeds the number of available control surfaces. The effectiveness of the control is measured by the magnitude of a propagating acoustic response vector, which is related to the circumferentially averaged sound pressure level (SPL), and is minimized by a standard least-squares minimization procedure.

Kousen, Kenneth A.↗
NASA NTRS2003

MOBY, A Radiometric Buoy for Performance Monitoring and Vicarious Calibration of Satellite Ocean Color Sensors: Measurement and Data Analysis Protocols

The Marine Optical Buoy (MOBY) is the centerpiece of the primary ocean measurement site for calibration of satellite ocean color sensors based on independent in situ measurements. Since late 1996, the time series of normalized water-leaving radiances L(sub WN)(lambda) determined from the array of radiometric sensors attached to MOBY are the primary basis for the on-orbit calibrations of the USA Sea-viewing Wide Field-of-view Sensor (SeaWiFS), the Japanese Ocean Color and Temperature Sensor (OCTS), the French Polarization Detection Environmental Radiometer (POLDER), the German Modular Optoelectronic Scanner on the Indian Research Satellite (IRS1-MOS), and the USA Moderate Resolution Imaging Spectrometer (MODIS). The MOBY vicarious calibration L(sub WN)(lambda) reference is an essential element in the international effort to develop a global, multi-year time series of consistently calibrated ocean color products using data from a wide variety of independent satellite sensors. A longstanding goal of the SeaWiFS and MODIS (Ocean) Science Teams is to determine satellite-derived L(sub WN)(labda) with a relative combined standard uncertainty of 5 %. Other satellite ocean color projects and the Sensor Intercomparison for Marine Biology and Interdisciplinary Oceanic Studies (SIMBIOS) project have also adopted this goal, at least implicitly. Because water-leaving radiance contributes at most 10 % of the total radiance measured by a satellite sensor above the atmosphere, a 5 % uncertainty in L(sub WN)(lambda) implies a 0.5 % uncertainty in the above-atmosphere radiance measurements. This level of uncertainty can only be approached using vicarious-calibration approaches as described below. In practice, this means that the satellite radiance responsivity is adjusted to achieve the best agreement, in a least-squares sense, for the L(sub WN)(lambda) results determined using the satellite and the independent optical sensors (e.g. MOBY). The end result of this approach is to implicitly absorb unquantified, but systematic, errors in the atmospheric correction, incident solar flux, and satellite sensor calibration into a single correction factor to produce consistency with the in situ data.

Clark, Dennis K.↗