Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical validation”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Coastal Atmosphere and Sea Time Series (CoASTS)

In this document, the first three years of a time series of bio-optical marine and atmospheric measurements are presented and analyzed. These measurements were performed from an oceanographic tower in the northern Adriatic Sea within the framework of the Coastal Atmosphere and Sea Time Series (CoASTS) project, an ocean color calibration and validation activity. The data set collected includes spectral measurements of the in-water apparent (diffuse attenuation coefficient, reflectance, Q-factor, etc.) and inherent (absorption and scattering coefficients) optical properties, as well as the concentrations of the main optical components (pigment and suspended matter concentrations). Clear seasonal patterns are exhibited by the marine quantities on which an appreciable short-term variability (on the order of a half day to one day) is superimposed. This short-term variability is well correlated with the changes in salinity at the surface resulting from the southward transport of freshwater coming from the northern rivers. Concentrations of chlorophyll alpha and total suspended matter span more than two orders of magnitude. The bio-optical characteristics of the measurement site pertain to both Case-I (about 64%) and Case-II (about 36%) waters, based on a relationship between the beam attenuation coefficient at 660nm and the chlorophyll alpha concentration. Empirical algorithms relating in-water remote sensing reflectance ratios and optical components or properties of interest (chlorophyll alpha, total suspended matter, and the diffuse attenuation coefficient) are presented.

Hooker, Stanford B.↗
NASA NTRS2017

Updated MISR Dark Water Research Aerosol Retrieval Algorithm - Part 1: Coupled 1.1 km Ocean Surface Chlorophyll a Retrievals with Empirical Calibration Corrections

As aerosol amount and type are key factors in the 'atmospheric correction' required for remote-sensing chlorophyll alpha concentration (Chl) retrievals, the Multi-angle Imaging SpectroRadiometer (MISR) can contribute to ocean color analysis despite a lack of spectral channels optimized for this application. Conversely, an improved ocean surface constraint should also improve MISR aerosol-type products, especially spectral single-scattering albedo (SSA) retrievals. We introduce a coupled, self-consistent retrieval of Chl together with aerosol over dark water. There are time-varying MISR radiometric calibration errors that significantly affect key spectral reflectance ratios used in the retrievals. Therefore, we also develop and apply new calibration corrections to the MISR top-of-atmosphere (TOA) reflectance data, based on comparisons with coincident MODIS (Moderate Resolution Imaging Spectroradiometer) observations and trend analysis of the MISR TOA bidirectional reflectance factors (BRFs) over three pseudo-invariant desert sites. We run the MISR research retrieval algorithm (RA) with the corrected MISR reflectances to generate MISR-retrieved Chl and compare the MISR Chl values to a set of 49 coincident SeaBASS (SeaWiFS Bio-optical Archive and Storage System) in situ observations. Where Chl(sub in situ) less than 1.5 mg m(exp -3), the results from our Chl model are expected to be of highest quality, due to algorithmic assumption validity. Comparing MISR RA Chl to the 49 coincident SeaBASS observations, we report a correlation coefficient (r) of 0.86, a root-mean-square error (RMSE) of 0.25, and a median absolute error (MAE) of 0.10. Statistically, a two-sample Kolmogorov- Smirnov test indicates that it is not possible to distinguish between MISR Chl and available SeaBASS in situ Chl values (p greater than 0.1). We also compare MODIS-Terra and MISR RA Chl statistically, over much broader regions. With about 1.5 million MISR-MODIS collocations having MODIS Chl less than 1.5 mg m(exp -3), MISR and MODIS show very good agreement: r = 0.96, MAE = 0.09, and RMSE = 0.15. The new dark water aerosol/Chl RA can retrieve Chl in low-Chl, case I waters, independent of other imagers such as MODIS, via a largely physical algorithm, compared to the commonly applied statistical ones. At a minimum, MISR's multi-angle data should help reduce uncertainties in the MODIS-Terra ocean color retrieval where coincident measurements are made, while also allowing for a more robust retrieval of particle properties such as spectral single-scattering albedo.

Limbacher, James A.↗
NASA NTRS2011

Noise of Embedded High Aspect Ratio Nozzles

A family of high aspect ratio nozzles were designed to provide a parametric database of canonical embedded propulsion concepts. Nozzle throat geometries with aspect ratios of 2:1, 4:1, and 8:1 were chosen, all with convergent nozzle areas. The transition from the typical round duct to the rectangular nozzle was designed very carefully to produce a flow at the nozzle exit that was uniform and free from swirl. Once the basic rectangular nozzles were designed, external features common to embedded propulsion systems were added: extended lower lip (a.k.a. bevel, aft deck), differing sidewalls, and chevrons. For the latter detailed Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations were made to predict the thrust performance and to optimize parameters such as bevel length, and chevron penetration and azimuthal curvature. Seventeen of these nozzles were fabricated at a scale providing a 2.13 inch diameter equivalent area throat." ! The seventeen nozzles were tested for far-field noise and a few data were presented here on the effect of aspect ratio, bevel length, and chevron count and penetration. The sound field of the 2:1 aspect ratio rectangular jet was very nearly axisymmetric, but the 4:1 and 8:1 were not, the noise on their minor axes being louder than the major axes. Adding bevel length increased the noise of these nozzles, especially on their minor axes, both toward the long and short sides of the beveled nozzle. Chevrons were only added to the 2:1 rectangular jet. Adding 4 chevrons per wide side produced some decrease at aft angles, but increased the high frequency noise at right angles to the jet flow. This trend increased with increasing chevron penetration. Doubling the number of chevrons while maintaining their penetration decreased these effects. Empirical models of the parametric effect of these nozzles were constructed and quantify the trends stated above." Because it is the objective of the Supersonics Project that future design work be done more by physics-based computations and less by experiments, several codes under development were evaluated against these test cases. Preliminary results show that the RANS-based code JeNo predicts the spectral directivity of the low aspect ratio jets well, but has no capability to predict the non-axisymmetry. An effort to address this limitations, used in the RANS-based code of Leib and Goldstein, overpredicted the impact of aspect ratio. The broadband shock noise code RISN, also limited to axisymmetric assumptions, did a good job of predicting the spectral directivity of underexpanded 2:1 cold jet case but was not as successful on high aspect ratio jets, particularly when they are hot. All results are preliminary because the underlying CFD has not been validated yet. An effort using a Large Eddy Simulation code by Stanford University predicted noise that agreed with experiments to within a few dB.

Bridges, James E.↗
NASA NTRS2012

Determining Damping Trends from a Range of Cable Harness Assemblies on a Launch Vehicle Panel from Test Measurements

The team of authors at Marshall Space Flight Center (MSFC) has been investigating estimating techniques for the vibration response of launch vehicle panels excited by acoustics and/or aero-fluctuating pressures. Validation of the approaches used to estimate these environments based on ground tests of flight like hardware is of major importance to new vehicle programs. The team at MSFC has recently expanded upon the first series of ground test cases completed in December 2010. The follow on tests recently completed are intended to illustrate differences in damping that might be expected when cable harnesses are added to the configurations under test. This validation study examines the effect on vibroacoustic response resulting from the installation of cable bundles on a curved orthogrid panel. Of interest is the level of damping provided by the installation of the cable bundles and whether this damping could be potentially leveraged in launch vehicle design. The results of this test are compared with baseline acoustic response tests without cables. Damping estimates from the measured response data are made using a new software tool that employs a finite element model (FEM) of the panel in conjunction with advanced optimization techniques. This paper will report on the \damping trend differences. observed from response measurements for several different configurations of cable harnesses. The data should assist vibroacoustics engineers to make more informed damping assumptions when calculating vibration response estimates when using model based analysis approach. Achieving conservative estimates that have more flight like accuracy is desired. The paper may also assist analysts in determining how ground test data may relate to expected flight response levels. Empirical response estimates may also need to be adjusted if the measured response used as an input to the study came from a test article without flight like cable harnesses.

Smith, Andrew↗
NASA NTRS1993

Enhanced methods for determining operational capabilities and support costs of proposed space systems

This report documents the work accomplished during the first two years of research to provide support to NASA in predicting operational and support parameters and costs of proposed space systems. The first year's research developed a methodology for deriving reliability and maintainability (R & M) parameters based upon the use of regression analysis to establish empirical relationships between performance and design specifications and corresponding mean times of failure and repair. The second year focused on enhancements to the methodology, increased scope of the model, and software improvements. This follow-on effort expands the prediction of R & M parameters and their effect on the operations and support of space transportation vehicles to include other system components such as booster rockets and external fuel tanks. It also increases the scope of the methodology and the capabilities of the model as implemented by the software. The focus is on the failure and repair of major subsystems and their impact on vehicle reliability, turn times, maintenance manpower, and repairable spares requirements. The report documents the data utilized in this study, outlines the general methodology for estimating and relating R&M parameters, presents the analyses and results of application to the initial data base, and describes the implementation of the methodology through the use of a computer model. The report concludes with a discussion on validation and a summary of the research findings and results.

Ebeling, Charles↗
NASA NTRS1990

Multiband radar characterization of forest biomes

The utility of airborne and orbital SAR in classification, assessment, and monitoring of forest biomes is investigated through analysis of orbital synthetic aperature radar (SAR) and multifrequency and multipolarized airborne SAR imagery relying on image tone and texture. Preliminary airborne SAR experiments and truck-mounted scatterometer observations demonstrated that the three dimensional structural complexity of a forest, and the various scales of temporal dynamics in the microwave dielectric properties of both trees and the underlying substrate would severely limit empirical or semi-empirical approaches. As a consequence, it became necessary to develop a more profound understanding of the electromagnetic properties of a forest scene and their temporal dynamics through controlled experimentation coupled with theoretical development and verification. The concatenation of various models into a physically-based composite model treating the entire forest scene became the major objective of the study as this is the key to development of a series of robust retrieval algorithms for forest biophysical properties. In order to verify the performance of the component elements of the composite model, a series of controlled laboratory and field experiments were undertaken to: (1) develop techniques to measure the microwave dielectric properties of vegetation; (2) relate the microwave dielectric properties of vegetation to more readily measured characteristics such as density and moisture content; (3) calculate the radar cross-section of leaves, and cylinders; (4) improve backscatter models for rough surfaces; and (5) relate attenuation and phase delays during propagation through canopies to canopy properties. These modeling efforts, as validated by the measurements, were incorporated within a larger model known as the Michigan Microwave Canopy Scattering (MIMICS) Model.

Dobson, M. Craig↗
NASA NTRS2011

Hybrid CFD/CAA Modeling for Liftoff Acoustic Predictions

This paper presents development efforts at the NASA Marshall Space flight Center to establish a hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) simulation system for launch vehicle liftoff acoustics environment analysis. Acoustic prediction engineering tools based on empirical jet acoustic strength and directivity models or scaled historical measurements are of limited value in efforts to proactively design and optimize launch vehicles and launch facility configurations for liftoff acoustics. CFD based modeling approaches are now able to capture the important details of vehicle specific plume flow environment, identifY the noise generation sources, and allow assessment of the influence of launch pad geometric details and sound mitigation measures such as water injection. However, CFD methodologies are numerically too dissipative to accurately capture the propagation of the acoustic waves in the large CFD models. The hybrid CFD/CAA approach combines the high-fidelity CFD analysis capable of identifYing the acoustic sources with a fast and efficient Boundary Element Method (BEM) that accurately propagates the acoustic field from the source locations. The BEM approach was chosen for its ability to properly account for reflections and scattering of acoustic waves from launch pad structures. The paper will present an overview of the technology components of the CFD/CAA framework and discuss plans for demonstration and validation against test data.

Strutzenberg, Louise L.↗
NASA NTRS2020

Uncertainties in CERES Top-of-Atmosphere Fluxes Caused by Changes in Accompanying Imager

The Clouds and the Earth’s Radiant Energy System (CERES) project provides observations of Earth’s radiation budget using measurements from CERES instruments on board the Terra, Aqua, Suomi National Polar-orbiting Partnership (S-NPP), and NOAA-20 satellites. The CERES top-of-atmosphere (TOA) fluxes are produced by converting radiance measurements using empirical angular distribution models, which are functions of cloud properties that are retrieved from imagers flying with the CERES instruments. As the objective is to create a long-term climate data record, not only calibration consistency of the six CERES instruments needs to be maintained for the entire time period, it is also important to maintain the consistency of other input data sets used to produce this climate data record. In this paper, we address aspects that could potentially affect the CERES TOA flux data quality. Discontinuities in imager calibration can affect cloud retrieval which can lead to erroneous flux trends. When imposing an artificial 0.6 per decade decreasing trend to cloud optical depth, which is similar to the trend difference between CERES Edition 2 and Edition 4 cloud retrievals, the decadal SW flux trend changed from −0.3 5 ± 0.18 Wm −2 to 0.61 ± 0.18 Wm −2 . This indicates that a 13% change in cloud optical depth results in about 1% change in the SW flux. Furthermore, different CERES instruments provide valid fluxes at different viewing zenith angle ranges, and including fluxes derived at the most oblique angels unique to S-NPP (>66 ∘ ) can lead to differences of 0.8 Wm −2 and 0.3 Wm −2 in global monthly mean instantaneous SW flux and LW flux. To ensure continuity, the viewing zenith angle ranges common to all CERES instruments (<66 ∘ ) are used to produce the long-term Earth’s radiation budget climate data record. The consistency of cloud properties retrieved from different imagers also needs to be maintained to ensure the TOA flux consistency.

Wenying Su↗
NASA NTRS2011

A Source-Term Based Boundary Layer Bleed/Effusion Model for Passive Shock Control

A modeling framework for boundary layer effusion has been developed based on the use of source (or sink) terms instead of the usual practice of specifying bleed directly as a boundary condition. This framework allows the surface boundary condition (i.e. isothermal wall, adiabatic wall, slip wall, etc.) to remain unaltered in the presence of bleed. This approach also lends itself to easily permit the addition of empirical models for second order effects that are not easily accounted for by simply defining effective transpiration values. Two effusion models formulated for supersonic flows have been implemented into this framework; the Doerffer/Bohning law and the Slater formulation. These models were applied to unit problems that contain key aspects of the flow physics applicable to bleed systems designed for hypersonic air-breathing propulsion systems. The ability of each model to predict bulk bleed properties was assessed, as well as the response of the boundary layer as it passes through and downstream of a porous bleed system. The model assessment was performed with and without the presence of shock waves. Three-dimensional CFD simulations that included the geometric details of the porous plate bleed systems were also carried out to supplement the experimental data, and provide additional insights into the bleed flow physics. Overall, both bleed formulations fared well for the tests performed in this study. However, the sample of test problems considered in this effort was not large enough to permit a comprehensive validation of the models.

Baurle, Robert A.↗
NASA NTRS1999

The Prospect for Remote Sensing of Cirrus Clouds with a Submillimeter-Wave Spectrometer

Given the substantial radiative effects of cirrus clouds and the need to validate cirrus cloud mass in climate models, it is important to measure the global distribution of cirrus properties with satellite remote sensing. Existing cirrus remote sensing techniques, such as solar reflectance methods, measure cirrus ice water path (IWP) rather indirectly and with limited accuracy. Submillimeter/wave radiometry is an independent method of cirrus remote sensing based on ice particles scattering the upwelling radiance emitted by the lower atmosphere. A new aircraft instrument, the Far Infrared Sensor for Cirrus (FIRSC), is described. The FIRSC employs a Fourier Transform Spectrometer (FTS). which measures the upwelling radiance across the whole submillimeter region (0.1 1.0-mm wavelength). This wide spectral coverage gives high sensitivity to most cirrus particle sizes and allows accurate determination of the characteristic particle size. Radiative transfer modeling is performed to analyze the capabilities of the submillimeter FTS technique. A linear inversion analysis is done to show that cirrus IWP, particle size, and upper-tropospheric temperature and water vapor may be accurately measured, A nonlinear statistical algorithm is developed using a database of 20000 spectra simulated by randomly varying most relevant cirrus and atmospheric parameters. An empirical orthogonal function analysis reduces the 500-point spectrum (20 - 70/cm) to 15 "pseudo-channels" that are then input to a neural network to retrieve cirrus IWP and median particle diameter. A Monte Carlo accuracy study is performed with simulated spectra having realistic noise. The retrieval errors are low for IWP (rms less than a factor of 1.5) and for particle sizes (rins less than 30%) for IWP greater than 5 g/sq m and a wide range of median particle sizes. This detailed modeling indicates that there is good potential to accurately measure cirrus properties with a submillimeter FTS.

Evans, K. Franklin↗
NASA NTRS1993

Establishment of computerized numerical databases on thermophysical and other properties of molten as well as solid materials and data evaluation and validation for generating recommended reliable reference data

The Center for Information and Numerical Data Analysis and Synthesis, (CINDAS), measures and maintains databases on thermophysical, thermoradiative, mechanical, optical, electronic, ablation, and physical properties of materials. Emphasis is on aerospace structural materials especially composites and on infrared detector/sensor materials. Within CINDAS, the Department of Defense sponsors at Purdue several centers: the High Temperature Material Information Analysis Center (HTMIAC), the Ceramics Information Analysis Center (CIAC) and the Metals Information Analysis Center (MIAC). The responsibilities of CINDAS are extremely broad encompassing basic and applied research, measurement of the properties of thin wires and thin foils as well as bulk materials, acquisition and search of world-wide literature, critical evaluation of data, generation of estimated values to fill data voids, investigation of constitutive, structural, processing, environmental, and rapid heating and loading effects, and dissemination of data. Liquids, gases, molten materials and solids are all considered. The responsibility of maintaining widely used databases includes data evaluation, analysis, correlation, and synthesis. Material property data recorded on the literature are often conflicting, diverging, and subject to large uncertainties. It is admittedly difficult to accurately measure materials properties. Systematic and random errors both enter. Some errors result from lack of characterization of the material itself (impurity effects). In some cases assumed boundary conditions corresponding to a theoretical model are not obtained in the experiments. Stray heat flows and losses must be accounted for. Some experimental methods are inappropriate and in other cases appropriate methods are carried out with poor technique. Conflicts in data may be resolved by curve fitting of the data to theoretical or empirical models or correlation in terms of various affecting parameters. Reasons (e.g. phase transitions) must be found for unusual dependence or any anomaly. Such critical evaluation involves knowledge of theory, experience in measurement, familiarity with metallurgy (microstructural behavior) and not inconsiderable judgment. An examination of typical data compiled and analyzed by CINDAS shows that the thermal conductivity of a material reported in the literature may vary by a factor of two of more; the range of reported values increases as temperature increases reflecting the difficulty of high temperature measurements. Often only estimates of melt behavior are available, despite the importance of melt properties in modeling, welding, or other solidification processes. There may be only a few measurements available for properties such as kinematic viscosity, even for widely used materials such as stainless steel. In the face of such a paucity of existing data and in a national environment where too few new data are being generated it is nonetheless the responsibility of CINDAS to select and disseminate recommended values of a wide variety of thermophysical properties.

Ho, C. Y.↗
NASA NTRS1999

A Critical Examination of Wind-Wave Spectral Functional Form

Traditionally, data from random ocean waves are presented in spectral functions. The spectra are the result of Fourier analysis. Fourier spectral analysis has dominated data analysis for, at least, the last hundred years. It has been the standard method for is examining the global amplitude-frequency distributions. Although Fourier transform valid under extremely general conditions, there are some crucial restrictions for the Fourier spectral analysis. The system must be linear, and the data must be stationary- otherwise, the resulting spectrum will make little physical sense. The stationarity requirement is also a common required criterion for most of other available data analysis methods. Nevertheless, few, if any, natural phenomena are linear and stationary. To compound these complications is the imperfection of our probes or numerical schemes the interactions of the imperfect probes even with a perfect linear system can make the final data nonlinear. Furthermore, all the available data are usually of finite duration. Under these conditions, Fourier analysis is of limited use, For lack of alternatives, however, Fourier analysis is still used to process such data. The loose application of Fourier analysis and the insouciant adoption of the stationary and linear assumptions may lead to misleading conclusions. Ocean waves are know to be nonlinear, and the wind system generating the wave field are seldom stationary- As a result, the traditional examination of the spectral form hardly made physical sense. A new method for analyzing nonlinear and nonstationary data has been developed. The key part is the Empirical Mode Decomposition (EMD) method with which any complicated data set can be decomposed into a finite and often small number of Intrinsic Mode Functions (IMF) that serve as the basis of the representation of the data, This decomposition method is adaptive, and, therefore, highly efficient. The IMFs admit well-behaved Hilbert transforms, and yield instantaneous energy and frequency as functions of time that give sharp identifications of imbedded structures. The final presentation of the results is an energy-frequency-time distribution, designated as the Hilbert Spectrum. Among the main conceptual innovations is the introduction of the instantaneous frequencies for complicated data sets, which eliminate the need of spurious harmonics to represent nonlinear and nonstationary signals. The spectral forms of the ocean waves are very different. This new method offers the first physical meaningful view of the wave spectrum. Data from laboratory and field will be presented to illustrate the differences.

Huang, Norden E.↗
NASA NTRS2014

Control and Non-Payload Communications Generation 1 Prototype Radio Flight Test Report

Unmanned aircraft (UA) represent a new capability that will provide a variety of services in the Government (public) and commercial (civil) aviation sectors. The growth of this potential industry has not yet been realized because of the lack of a common understanding of what is required to safely operate Unmanned Aircraft Systems in the National Airspace System (UAS in the NAS). The desire and ability to fly UA is of increasing urgency. The application of UA to perform national security, defense, scientific, and emergency management are driving the critical need for less restrictive access by UA to the NAS. Existing Federal Aviation Regulations, procedures, and technologies do not allow routine UA access to the NAS. Access to the NAS is hampered by challenges such as the lack of an onboard pilot to see and avoid other aircraft; the ability of a single pilot or operator to control multiple UA; the reliance on command and control (C2) links; the altitudes, speeds, and duration at which the aircraft fly; and the wide variation in UA size and performance. NASA is working with other Government agencies to provide solutions that reduce technical barriers and make access to the NAS routine. This goal will be accomplished through system-level integration of key concepts, technologies, or procedures and through demonstrations of these integrated capabilities in an operationally relevant environment. This project provides an opportunity to transition the acquired empirical data and knowledge to the Federal Aviation Administration and other stakeholders to help them define the requirements for routine UA access to the NAS.Radio communications channels for UA are currently managed through exceptions and use either Department of Defense frequencies for line-of-sight (LOS) and satellite-based communications links, low-power LOS links in amateur bands, or unlicensed Industrial/Scientific/Medical (ISM) frequencies. None of these frequency bands are designated for safety and regularity of flight. Only recently has radiofrequency (RF) spectrum been allocated by the International Telecommunications Union specifically for commercial UA C2, LOS communication (L-Band: 960 to 1164 MHz, and C-Band: 5030 to 5091 MHz). The safe and efficient integration of UA into the NAS requires the use of protected RF spectrum allocations and a new data communications system that is both secure and scalable to accommodate the potential growth of these new aircraft. Data communications for UA-referred to as control and non-payload communications (CNPC)-will be used to exchange information between a UA and a ground station (GS) to ensure safe, reliable, and effective UA flight operation. The focus of this effort is on validating and allocating new RF spectrum and data link communications to enable civil UA integration into the NAS. Through a cost-sharing cooperative agreement with Rockwell Collins, Inc., the NASA Glenn Research Center is exploring and performing the necessary development steps to realize a prototype UA CNPC system. These activities include investigating signal waveforms and access techniques, developing representative CNPC radio hardware, and executing relevant testing and validation activities. There is no intent to manufacture the CNPC end product, rather the goals are to study, demonstrate, and validate a typical CNPC system that will allow safe and efficient communications within the L-Band and C-Band spectrum allocations. The system is addressing initial "seed" requirements from RTCA, Inc., Special Committee 203 (SC-203) and is on a path to Federal Aviation Administration certification. This report provides results from the flight testing campaign of the Rockwell Collins Generation 1 prototype radio, referred hereafter as the "radio." The radio sets operate within the 960- to 977-MHz frequency band with both air and ground radios using identical hardware. Flight tests involved one aircraft and one GS. Results include discussion of aircraft flight paths and associated radio performance.

Shalkhauser, Kurt A.↗
NASA NTRS2004

Time Changes of the European Gravity Field from GRACE: A Comparison with Ground Measurements from Superconducting Gravimeters and with Hydrology Model Predictions

We investigate the time-variable gravity changes in Europe retrieved from the initial GRACE monthly solutions spanning a 18 month duration from April 2002 to October 2003. Gravity anomaly maps are retrieved in Central Europe from the monthly satellite solutions we compare the fields according to various truncation levels (typically between degree 10 and 20) of the initial fields (expressed in spherical harmonics to degree 120). For these different degrees, an empirical orthogonal function (EOF) decomposition of the time-variable gravity field leads us to its main spatial and temporal characteristics. We show that the dominant signal is found to be annual with an amplitude and a phase both in agreement with predictions in Europe modeled using snow and soil-moisture variations from recent hydrology models. We compare these GRACE gravity field changes to surface gravity observations from 6 superconducting gravimeters of the GGP (Global Geodynamics Project) European sub-network, with a special attention to loading corrections. Initial results suggest that all 3 data sets (GRACE, hydrology and GGP) are responding to annual changes in near-surface water in Europe of a few microGal (at length scales of approx.1000 km) that show a high value in winter and a summer minimum. We also point out that the GRACE gravity field evolution seems to indicate that there is a trend in gravity between summer 2002 and summer 2003 which can be related to the 2003 heatwave in Europe and its hydrological consequences (drought). Despite the limited time span of our analysis and the uncertainties in retrieving a regional solution from the network of gravimeters, the calibration and validation aspects of the GRACE data processing based on the annual hydrology cycle in Europe are in progress.

Hinderer, J.↗
NASA NTRSDate not supplied

On the use of DSCOVR EPIC to monitor the visible calibration stability of polar orbiting imagers to improve the next edition of the CERES climate data record.

The NASA CERES project has provided the climate quality observed TOA and computed surface fluxes to the scientific community. CERES instruments are onboard the Terra, Aqua, NPP, and NOAA-20 spacecraft. CERES uses MODIS and VIIRS imagers to retrieve cloud properties needed to convert CERES footprint radiance observations into fluxes using empirically derived angular directional models obtained during early CERES record. CERES utilizes geostationary sensors to infer the regional diurnal fluxes in between Terra and Aqua CERES observations. The imagers are also used to radiometrically scale the geostationary sensors (GEO) radiances to the imager calibration reference to ensure that the GEO derived cloud properties and broadband TOA fluxes are consistent in both space and time. Both the imager and GEO retrieved cloud properties are used to compute the surface fluxes. The Aqua-MODIS, SNPP-VIIRS, and NOAA-20 VIIRS afternoon imagers will also need to be radiometrically scaled to the same common calibration reference. Although all three imagers employ onboard solar diffusers, the calibration is not consistent over time due to instrument anomalies and ageing of the onboard calibrators. The CERES imager and GEO calibration group (IGCG) has been tasked to remove the imager channel calibration drifts for the next CERES reprocessing effort. The team will rely primarily on deep convective clouds, desert, and polar ice invariant targets to monitor the imager channel stability. The team uses ray-matched radiance pairs to radiometrically scale the SNPP and NOAA-20 VIIRS sensors to Aqua-MODIS. The radiometric scaling is further validated using geostationary imagers as transfer radiometers. The DSCOVR satellite was launched on February 25, 2015 and orbits around the L1 Lagrange point directly between the Earth and the sun. The EPIC instrument onboard DSCOVR employs a CCD array to image the Earth approximately every 2-hours. The EPIC sensor contains no onboard calibration systems. However, multiple inter-calibration studies have shown that the EPIC imager is very stable in time. This is likely due to the DSCOVR satellite being located about a 1.5M km from the Earth, where very little Earth reflected solar radiation degrades the optics. The excellent radiometric stability of EPIC allows the CERES IGCG to utilize the EPIC observations as a stable reference for monitoring the calibration stability of the three afternoon imagers, as well as to validate the radiometric scaling factors between them. Examples of the imager relative calibration using EPIC before and after radiometric scaling will be shown along with the results from the use of invariant targets to remove imager calibration drifts.

DSCOVR-EPIC↗
NASA NTRS2018

Estimation of Lower-Body Kinetics from Loading Profile and Kinematics Alone, Without Measured Ground Reaction Forces

Biomechanical models of human motion can estimate kinetic outcomes, such as joint moments, joint forces and muscle forces. Typically, one performs an inverse dynamics (ID) analysis to compute joint moments from joint angles and measured external forces. Sometimes it is impractical to measure ground reaction forces and moments (GRF&M). We devised an empirical method for performing ID analysis of resistance exercises without measured GRF&M. The method solves the multibody dynamics equations of motion with four key assumptions about the GRF&M that reduce the number of unknowns. The assumptions are 1) negligible ground reaction moments, 2) fixed lateral/medial location of the center of pressure (COP), 3) equal fore/aft location of the COP between the feet, and 4) constant angle of the GRF vector relative to the vertical axis in the frontal plane. We used evaluation trials from a spaceflight countermeasure resistance training device to test this approach. Four participants performed squat and deadlift exercises at various loads. We compared results from traditional ID analysis to results without measured GRF&M using our method. We found that joint moment trajectories in the sagittal plane were qualitatively similar in shape between the two methods, and the amount of root mean squared error (RMSE), measured by difference in joint moment impulse, was typically under 10 percent. Non-sagittal joint moment trajectories, which are much lower in overall magnitude, were not qualitatively similar in shape between the two methods. Non-sagittal moments displayed much higher RMSE, with typical values well over 50 percent. These findings were further supported by validation metrics (Sprague and Geers' P and M metrics, Pearson's r correlation coefficient). Based on these findings, we concluded that useful kinetic results are obtained from ID analysis of squat and deadlift exercises, even when GRF&M are not measured, as long as the outcomes of interest lie in the sagittal plane.

computational modeling↗
NASA NTRS2021

Spectral attenuation coefficients from measurements of light transmission in bare ice on the Greenland Ice Sheet

Light transmission into bare glacial ice affects surface energy balance, biophotochemistry, and light detection and ranging (lidar) laser elevation measurements but has not previously been reported for the Greenland Ice Sheet. We present measurements of spectral transmittance at 350–900 nm in bare glacial ice collected at a field site in the western Greenland ablation zone (67.15∘ N, 50.02∘ W). Empirical irradiance attenuation coefficients at 350–750 nm are ∼ 0.9–8.0 m−1 for ice at 12–124 cm depth. The absorption minimum is at ∼ 390–397 nm, in agreement with snow transmission measurements in Antarctica and optical mapping of deep ice at the South Pole. From 350–530 nm, our empirical attenuation coefficients are nearly 1 order of magnitude larger than theoretical values for optically pure ice. The estimated absorption coefficient at 400 nm suggests the ice volume contained a light-absorbing particle concentration equivalent to ∼ 1–2 parts per billion (ppb) of black carbon, which is similar to pre-industrial values found in remote polar snow. The equivalent mineral dust concentration is ∼ 300–600 ppb, which is similar to values for Northern Hemisphere warm periods with low aeolian activity inferred from ice cores. For a layer of quasi-granular white ice (weathering crust) extending from the surface to ∼ 10 cm depth, attenuation coefficients are 1.5 to 4 times larger than for deeper bubbly ice. Owing to higher attenuation in this layer of near-surface granular ice, optical penetration depth at 532 nm is 14 cm (20 %) lower than asymptotic attenuation lengths for optically pure bubbly ice. In addition to the traditional concept of light scattering on air bubbles, our results imply that the granular near-surface ice microstructure of weathering crust is an important control on radiative transfer in bare ice on the Greenland Ice Sheet ablation zone, and we provide new values of flux attenuation, absorption, and scattering coefficients to support model development and validation.

IceSat-2↗
NASA NTRS2001

GCM Simulations of the Aerosol Indirect Effect: Sensitivity to Cloud Parameterization and Aerosol Burden

We describe the coupling of the Goddard Institute for Space Studies (GISS) general circulation model (GCM) to an online sulfur chemistry model and source models for organic matter and sea-salt that is used to estimate the aerosol indirect effect. The cloud droplet number concentration is diagnosed empirically from field experiment datasets over land and ocean that observe droplet number and all three aerosol types simultaneously; corrections are made for implied variations in cloud turbulence levels. The resulting cloud droplet number is used to calculate variations in droplet effective radius, which in turn allows us to predict aerosol effects on cloud optical thickness and microphysical process rates. We calculate the aerosol indirect effect by differencing the top-of-the-atmosphere net cloud radiative forcing for simulations with present-day vs. pre-industrial emissions. Both the first (radiative) and second (microphysical) indirect effects are explored. We test the sensitivity of our results to cloud parameterization assumptions that control the vertical distribution of cloud occurrence, the autoconversion rate, and the aerosol scavenging rate, each of which feeds back significantly on the model aerosol burden. The global mean aerosol indirect effect for all three aerosol types ranges from -1.55 to -4.36 W/sq m in our simulations. The results are quite sensitive to the pre-industrial background aerosol burden, with low pre-industrial burdens giving strong indirect effects, and to a lesser extent to the anthropogenic aerosol burden, with large burdens giving somewhat larger indirect effects. Because of this dependence on the background aerosol, model diagnostics such as albedo-particle size correlations and column cloud susceptibility, for which satellite validation products are available, are not good predictors of the resulting indirect effect.

Menon, Surabi↗