Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coupled inversion”

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 271 records · Page 15

The double universal joint wrist on a manipulator: Solution of inverse position kinematics and singularity analysis

This paper presents three methods to solve the inverse position kinematics position problem of the double universal joint attached to a manipulator: (1) an analytical solution for two specific cases; (2) an approximate closed form solution based on ignoring the wrist offset; and (3) an iterative method which repeats closed form position and orientation calculations until the solution is achieved. Several manipulators are used to demonstrate the solution methods: cartesian, cylindrical, spherical, and an anthropomorphic articulated arm, based on the Flight Telerobotic Servicer (FTS) arm. A singularity analysis is presented for the double universal joint wrist attached to the above manipulator arms. While the double universal joint wrist standing alone is singularity-free in orientation, the singularity analysis indicates the presence of coupled position/orientation singularities of the spherical and articulated manipulators with the wrist. The cartesian and cylindrical manipulators with the double universal joint wrist were found to be singularity-free. The methods of this paper can be implemented in a real-time controller for manipulators with the double universal joint wrist. Such mechanically dextrous systems could be used in telerobotic and industrial applications, but further work is required to avoid the singularities.

Williams, Robert L., III↗

Kinetics model for initiation and promotion for describing tumor prevalence from HZE radiation

A kinetics model for cellular repair and misrepair for multiple radiation-induced lesions (mutation-inactivation) is coupled to a two-mutation model of initiation and promotion in tissue to provide a parametric description of tumor prevalence in the Harderian gland in a mouse. Dose-response curves are described for gamma-rays and relativistic ions. The effects of nuclear fragmentation are also considered for high-energy proton and alpha particle exposures The model described provides a parametric description of age-dependent cancer induction for a wide range of radiation fields. We also consider the two hypotheses that radiation acts either solely as an initiator or as both initiator and promoter and make model calculations for fractionation exposures from gamma-rays and relativistic Fe ions. For fractionated Fe exposures, an inverse dose-rate effect is provided by a promotion hypothesis using a mutation rate for promotion typical of single-gene mutations.

Cucinotta, Francis A.↗

Benchmark solutions for the galactic ion transport equations with spatial and energy coupling

In order to anticipate future space shielding requirements, NASA has initiated an effort to formulate computational methods to simulate radiation effects in space. As part of the program, numerical transport algorithms have been developed for the deterministic Boltzman equation describing galactic cosmic ray (GCR) interactions with matter. It thus becomes necessary to assess the accuracy of proposed deterministic algorithms. For this reason, analytical benchmark solutions to mathematically tractable galactic cosmic ray equations have recently been obtained. Even though these problems involve simplifying assumptions of the associated physics, they still contain the essential features of the basic transport processes. The solutions obtained are features of the basic transport processes. The solutions obtained are compared to results from numerical algorithms in order to ensure proper coding and to provide a measure of the accuracy of the numerical methods used in the algorithm. For the first time, mathematical methods have been applied to the galactic ion transport (GIT) equations in the straight ahead approximation with constant nuclear properties. The approach utilizes a Laplace transforms inversion yielding a closed form benchmark solution which is also computationally efficient.

Ganapol, Barry D.↗

Causes of slowing‐down seasonal CO2 amplitude at Mauna Loa

Changing amplitude of the seasonal cycle of atmospheric CO2 (SCA) in the northern hemisphere is an emerging carbon cycle property. Mauna Loa (MLO) station (20°N, 156°W), which has the longest continuous northern hemisphere CO2 record, shows an increasing SCA before the 1980s (p < .01), followed by no significant change thereafter. We analyzed the potential driving factors of SCA slowing‐down, with an ensemble of dynamic global vegetation models (DGVMs) coupled with an atmospheric transport model. We found that slowing‐down of SCA at MLO is primarily explained by response of net biome productivity (NBP) to climate change, and by changes in atmospheric circulations. Through NBP, climate change increases SCA at MLO before the 1980s and decreases it afterwards. The effect of climate change on the slowing‐down of SCA at MLO is mainly exerted by intensified drought stress acting to offset the acceleration driven by CO2 fertilization. This challenges the view that CO2 fertilization is the dominant cause of emergent SCA trends at northern sites south of 40°N. The contribution of agricultural intensification on the deceleration of SCA at MLO was elusive according to land–atmosphere CO2 flux estimated by DGVMs and atmospheric inversions. Our results also show the necessity to adequately account for changing circulation patterns in understanding carbon cycle dynamics observed from atmospheric observations and in using these observations to benchmark DGVMs.

Mauna Loa↗

Handling Qualities Evaluations of Low Complexity Model Reference Adaptive Controllers for Reduced Pitch and Roll Damping Scenarios

National Aeronautics and Space Administration (NASA) researchers have conducted a series of flight experiments designed to study the effects of varying levels of adaptive controller complexity on the performance and handling qualities of an aircraft under various simulated failure or damage conditions. A baseline, nonlinear dynamic inversion controller was augmented with three variations of a model reference adaptive control design. The simplest design consisted of a single adaptive parameter in each of the pitch and roll axes computed using a basic gradient-based update law. A second design was built upon the first by increasing the complexity of the update law. The third and most complex design added an additional adaptive parameter to each axis. Flight tests were conducted using NASA s Full-scale Advanced Systems Testbed, a highly modified F-18 aircraft that contains a research flight control system capable of housing advanced flight controls experiments. Each controller was evaluated against a suite of simulated failures and damage ranging from destabilization of the pitch and roll axes to significant coupling between the axes. Two pilots evaluated the three adaptive controllers as well as the non-adaptive baseline controller in a variety of dynamic maneuvers and precision flying tasks designed to uncover potential deficiencies in the handling qualities of the aircraft, and adverse interactions between the pilot and the adaptive controllers. The work was completed as part of the Integrated Resilient Aircraft Control Project under NASA s Aviation Safety Program.

Hanson, Curt↗

Investigation of Cloud Interior Structure Using a Tomographic Technique Based on Radon Transform. Application to the Measurements Made by the Research Scanning Polarimeter During CAMP2Ex

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP's uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid-water cloud droplet sizes using the rainbow structure observed in the polarized reflectance measurements. The RSP observations also provide for derivation of 2D field of extinction coefficient inside the cloud using a tomographic technique. This technique utilizes the family of cloud shapes derived using "cutout" technique and corresponding to a number of thresholds in total reflectance (separating bright cloud from its darker background). These shapes can be interpreted as level curves of an abstract "reflectance density", which is then used for derivation of the directional cloud optical thickness (dCOT) tomogram. The latter is a collection of dCOTs parameterized by the angles and offsets of the corresponding view rays (chords) relative to the cloud center. After this, the inverse Radon Transform (the mathematical basis of the X-ray com-puter tomography) is applied to the dCOT tomogram yielding 2D spatial distribution of the extinction coefficient. The later can be converted into droplet number concentration using the droplet size profiles derived from the RSP's polarized reflectance measurements. After tests on synthetic data this technique was applied to real RSP measurements from NASA's Cloud, Aerosol and Monsoon Processes Philippines Experiment (CAMP2Ex) conducted in the vicinity of the Philippines during the Southwest Monsoon (August–September 2019). We have investigated the interiors of a number of clouds observed during CAMP2Ex focusing on Cu and CuCg (Tcu) cases. Our retrievals were routinely compared with the correlative measurements by lidar (HSRL-2) and cloud radar (APR-3) that were deployed on the same airborne platform (NASA P-3B) during this field experiment.

Research Scanning Polarimeter↗

Modeling Atmosphere-Ocean Radiative Transfer: A PACE Mission Perspective

The research frontiers of radiative transfer (RT) in coupled atmosphere-ocean systems are explored to enable new science and specifically to support the upcoming Plankton, Aerosol, Cloud ocean Ecosystem (PACE) satellite mission. Given (i) the multitude of atmospheric and oceanic constituents at any given moment that each exhibits a large variety of physical and chemical properties and (ii) the diversity of light-matter interactions (scattering, absorption, and emission), tackling all outstanding RT aspects related to interpreting and/or simulating light reflected by atmosphere-ocean systems becomes impossible. Instead, we focus on both theoretical and experimental studies of RT topics important to the science threshold and goal questions of the PACE mission and the measurement capabilities of its instruments. We differentiate between (a) forward (FWD) RT studies that focus mainly on sensitivity to influencing variables and/or simulating data sets, and (b) inverse (INV) RT studies that also involve the retrieval of atmosphere and ocean parameters. Our topics cover (1) the ocean (i.e., water body): absorption and elastic/inelastic scattering by pure water (FWD RT) and models for scattering and absorption by particulates (FWD RT and INV RT); (2) the air-water interface: variations in ocean surface refractive index (INV RT) and in whitecap reflectance (INV RT); (3) the atmosphere: polarimetric and/or hyperspectral remote sensing of aerosols (INV RT) and of gases (FWD RT); and (4) atmosphere-ocean systems: benchmark comparisons, impact of the Earth's sphericity and adjacency effects on space-borne observations, and scattering in the ultraviolet regime (FWD RT). We provide for each topic a summary of past relevant (heritage) work, followed by a discussion (for unresolved questions) and RT updates.

Chowdhary, Jacek↗

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 4: Computer user's manual for UAAP turboprop aeroacoustic code

The Unified AeroAcoustic Program (UAAP) code calculates the airloads on a single rotation prop-fan, or propeller, and couples these airloads with an acoustic radiation theory, to provide estimates of near-field or far-field noise levels. The steady airloads can also be used to calculate the nonuniform velocity components in the propeller wake. The airloads are calculated using a three dimensional compressible panel method which considers the effects of thin, cambered, multiple blades which may be highly swept. These airloads may be either steady or unsteady. The acoustic model uses the blade thickness distribution and the steady or unsteady aerodynamic loads to calculate the acoustic radiation. The users manual for the UAAP code is divided into five sections: general code description; input description; output description; system description; and error codes. The user must have access to IMSL10 libraries (MATH and SFUN) for numerous calls made for Bessel functions and matrix inversion. For plotted output users must modify the dummy calls to plotting routines included in the code to system-specific calls appropriate to the user's installation.

Menthe, R. W.↗

On Boundary-Value Problems for RANS Equations and Two-Equation Turbulence Models

Currently, in engineering computations for high Reynolds number turbulent flows, turbulence modeling continues to be the most frequently used approach to represent the effects of turbulence. Such models generally rely on solving either one or two transport equations along with the Reynolds-Averaged Navier–Stokes (RANS) equations. The solution of the boundary-value problem of any system of partial differential equations requires the complete delineation of the equations and the boundary conditions, including any special restrictions and conditions. In the literature, such a description is often incomplete, neglecting important details related to the boundary conditions and possible restrictive conditions, such as how to ensure satisfying prescribed values of the dependent variables of the transport equations in the far field of a finite domain. In this article, we discuss the possible influence of boundary values, as well as near-field and far-field behavior, on the solution of the RANS equations coupled with transport equations for turbulence modeling. In so doing, we defne the concept of a welldefined boundary-value problem. Additionally, a three-dimensional, rather than a simpler one-dimensional analysis is performed to analyze the near-wall and far-field behavior of the turbulence model variables. This allows an assessment of the decay rate of these variables required to realize the boundary conditions in the far field. This paper also addresses the impact of various transformations of two-equation models (e.g., the model of Wilcox) to remove the singular behavior of the dissipation rate (ω) at the surface boundary. Finally, the issue of well-posedness regarding the governing equations is considered. A compelling argument (although not a proof) for ill-posedness is made for both direct and inverse problems.

Navier–Stokes equations↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Robert S. Schrom↗

Organic chemistry on Titan

Observations of nonequilibrium phenomena on the Saturn satellite Titan indicate the occurrence of organic chemical evolution. Greenhouse and thermal inversion models of Titan's atmosphere provide environmental constraints within which various pathways for organic chemical synthesis are assessed. Experimental results and theoretical modeling studies suggest that the organic chemistry of the satellite may be dominated by two atmospheric processes: energetic-particle bombardment and photochemistry. Reactions initiated in various levels of the atmosphere by cosmic ray, Saturn wind, and solar wind particle bombardment of a CH4 - N2 atmospheric mixture can account for the C2-hydrocarbons, the UV-visible-absorbing stratospheric haze, and the reddish color of the satellite. Photochemical reactions of CH4 can also account for the presence of C2-hydrocarbons. In the lower Titan atmosphere, photochemical processes will be important if surface temperatures are sufficiently high for gaseous NH3 to exist. Hot H-atom reactions initiated by photo-dissociation of NH3 can couple the chemical reactions of NH3 and CH4 and produce organic matter.

Chang, S.↗

Circadian Entrainment, Sleep-Wake Regulation and Neurobehavioral Performance During Extended Duration Space Flight

Long-duration manned space flight requires crew members to maintain a high level of cognitive performance and vigilance while operating and monitoring sophisticated instrumentation. However, the reduction in the strength of environmental synchronizers in the space environment leads to misalignment of circadian phase among crew members, coupled with restricted time available to sleep, results in sleep deprivation and consequent deterioration of neurobehavioral function. Crew members are provided, and presently use, long-acting benzodiazepine hypnotics on board the current, relatively brief space shuttle missions to counteract such sleep disruption, a situation that is only likely to worsen during extended duration missions. Given the known carry-over effects of such compounds on daytime performance, together with the reduction in emergency readiness associated with their use at night, NASA has recognized the need to develop effective but safe countermeasures to allow crew members to obtain an adequate amount of sleep. Over the past eight years, we have successfully implemented a new technology for shuttle crew members involving bright light exposure during the pre-launch period to facilitate adaptation of the circadian timing system to the inversions of the sleep-wake schedule often required during dual shift missions. However for long duration space station missions it will be necessary to develop effective and attainable countermeasures that can be used chronically to optimize circadian entrainment. Our current research effort is to study the effects of light-dark cycles with reduced zeitgeber strength, such as are anticipated during long-duration space flight, on the entrainment of the endogenous circadian timing system and to study the effects of a countermeasure that consists of scheduled brief exposures to bright light on the human circadian timing system. The proposed studies are designed to address the following Specific Aims: (1) test the hypothesis that synchronization of the human circadian pacemaker will be disturbed in men and women by the reduction in LD cycle strength. (2) test the hypothesis that this disturbed circadian synchronization will result in the secretion of the sleep-promoting hormone melatonin during the waking day, disturbed sleep, reduced growth hormone secretion, and impaired performance and daytime alertness; (3) as a countermeasure, test the hypothesis that brief daily exposures to bright light (10,000 lux) will reestablish normal entrained circadian phase, resulting in improved sleep consolidation, normalized sleep structure and endogenous growth hormone secretion and enhanced daytime performance. To date, we have carried out twelve experiments to address Hypotheses I and 2 and data analyses are in progress. The results of the current research may have important implications for the treatment of circadian rhythm sleep disorders, such as delayed sleep phase syndrome and shift-work dyssomnia, which are anticipated to have a high incidence and prevalence during extended duration space flight such as planned for the International Space Station and manned missions to Mars.

Czeisler, Charles A.↗

Improved Ionospheric Electrodynamic Models and Application to Calculating Joule Heating Rates

Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Weimer, D. R.↗

Analysis of the Impedance Resonance of Piezoelectric Multi-Fiber Composite Stacks

Multi-Fiber CompositesTM (MFC's) produced by Smart Materials Corp behave essentially like thin planar stacks where each piezoelectric layer is composed of a multitude of fibers. We investigate the suitability of using previously published inversion techniques for the impedance resonances of monolithic co-fired piezoelectric stacks to the MFCTM to determine the complex material constants from the impedance data. The impedance equations examined in this paper are those based on the derivation. The utility of resonance techniques to invert the impedance data to determine the small signal complex material constants are presented for a series of MFC's. The technique was applied to actuators with different geometries and the real coefficients were determined to be similar within changes of the boundary conditions due to change of geometry. The scatter in the imaginary coefficient was found to be larger. The technique was also applied to the same actuator type but manufactured in different batches with some design changes in the non active portion of the actuator and differences in the dielectric and the electromechanical coupling between the two batches were easily measureable. It is interesting to note that strain predicted by small signal impedance analysis is much lower than high field stains. Since the model is based on material properties rather than circuit constants, it could be used for the direct evaluation of specific aging or degradation mechanisms in the actuator as well as batch sorting and adjustment of manufacturing processes.

Multi-fiber composites↗

Investigations of Clouds and Aerosols on Mars, Venus and Titan

For Titan our interest during the past few years was to explain the observed asymmetry in the albedo. We suggested earlier, from one-dimensional modeling studies, that vertical transport rates were comparable to particle fall speeds. Since heating of the upper atmosphere, which drives dynamical motions, is largely due to the aerosols, a non- linear interaction between dynamics, radiative heating and particle microphysics is possible. We pursued this interaction in a two-dimensional model. We showed that the observed variations in the albedo between the two hemispheres and over an orbital cycle, could be due to dynamical motions suspending particles so that particle sizes and optical depths vary across the planet. In the Cassini time frame, future studies of this interaction between dynamics, radiation and microphysics may be worthwhile using the strong modeling base that others, and we have developed. In our recently approved proposal, however, we plan to extend our modeling to hydrocarbon clouds that lie at lower levels. We know very little about such clouds, and numerical models for their properties are non-existent. These clouds may be observed by the Huygen's Probe, and by the Cassini orbiter, so predictions of their properties should help in the analysis of Cassini data. We have also developed a sophisticated model for the lower, condensational, clouds on Venus. In this model we explored the water vapor budget on Venus, and the properties of the clouds such as particle size distribution. Most researchers have investigated the upper clouds on Venus, which are essentially a photochemical aerosol with a long lifetime. The lower clouds, however, are similar to stratus clouds on Earth. These clouds have short lifetimes, and are tightly coupled to the dynamics at the base of the Venus cloud deck. We believe that the holes in these clouds seen at near infrared wavelengths by Galileo are related to some interaction between dynamics and cloud physics. One goal of our recently approved proposed work is to better understand this interaction, and to better understand how these clouds may vary if the climate of Venus were to vary. During the past few years we have applied our model of the water ice clouds on Mars to new data sets from Pathfinder, and Mars Global Surveyor. We have compared our predictions of cloud properties with those seen by Pathfinder, and found reasonable agreement. More recently we have improved our model by including a radiative transfer algorithm and a boundary layer transport scheme . The goal of this work is to understand the thermal inversions seen in Mars Global Surveyor data. We find that the inversions are created by clouds radiatively cooling the atmospheric layers in which the clouds form. We have also initiated laboratory work on the physics of carbon dioxide particle formation. In our recently approved work we propose to use these laboratory data, and to extend our modeling to carbon dioxide clouds in the Martian atmosphere.

Toon, Owen B.↗

Mass Spectrometric and Langmuir Probe Measurements in Inductively Coupled Plasmas in Ar, CHF3/Ar and CHF3/Ar/O2 Mixtures

Absolute fluxes and energy distributions of ions in inductively coupled plasmas of Ar, CHF3/Ar, and CHF3/Ar/O2 have been measured. These plasmas were generated in a Gaseous Electronics Conference (GEC) cell modified for inductive coupling at pressures 10-50 mTorr and 100-300 W of 13.56 MHz radio frequency (RF) power in various feedgas mixtures. In pure Ar plasmas, the Ar(+) flux increases linearly with pressure as well as RF-power. Total ion flux in CHF3 mixtures decreases with increase in pressure and also CHF3 concentration. Relative ion fluxes observed in the present studies are analyzed with the help of available cross sections for electron impact ionization and charge-exchange ion-molecule reactions. Measurements of plasma potential, electron and ion number densities, electron energy distribution function, and mean electron energy have also been made in the center of the plasma with a RF compensated Langmuir probe. Plasma potential values are compared with the mean ion energies determined from the measured ion energy distributions and are consistent. Electron temperature, plasma potential, and mean ion energy vary inversely with pressure, but increase with CHF3 content in the mixture.

Kim, J. S.↗

OSOAA: A Vector Radiative Transfer Model of Coupled Atmosphere-Ocean System for a Rough Sea Surface Application to the Estimates of the Directional Variations of the Water Leaving Reflectance to Better Process Multi-angular Satellite Sensors Data Over the Ocean

In this study, we present a radiative transfer model, so-called OSOAA, that is able to predict the radiance and degree of polarization within the coupled atmosphere-ocean system in the presence of a rough sea surface. The OSOAA model solves the radiative transfer equation using the successive orders of scattering method. Comparisons with another operational radiative transfer model showed a satisfactory agreement within 0.8%. The OSOAA model has been designed with a graphical user interface to make it user friendly for the community. The radiance and degree of polarization are provided at any level, from the top of atmosphere to the ocean bottom. An application of the OSOAA model is carried out to quantify the directional variations of the water leaving reflectance and degree of polarization for phytoplankton and mineral-like dominated waters. The difference between the water leaving reflectance at a given geometry and that obtained for the nadir direction could reach 40%, thus questioning the Lambertian assumption of the sea surface that is used by inverse satellite algorithms dedicated to multi-angular sensors. It is shown as well that the directional features of the water leaving reflectance are weakly dependent on wind speed. The quantification of the directional variations of the water leaving reflectance obtained in this study should help to correctly exploit the satellite data that will be acquired by the current or forthcoming multi-angular satellite sensors.

Ocean surface↗