Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “delta method”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 847 records · Page 47

Bias correcting regional scale Earth system model projections: novel approach using empirical mode decomposition

Bias correction is a crucial step in using Earth system model outputs for assessments, as it adjusts systematic errors by comparing the model to observations. However, standard methods – ranging from mean-based linear scaling to distribution-based quantile mapping typically treat bias correction as a single-scale process, overlooking the fact that biases can manifest differently across daily, seasonal, and annual timescales. In this study, we propose a novel, timescale-aware bias-correction approach built on Empirical Mode Decomposition. By decomposing the meteorological signal into multiple oscillatory components and aggregating them to represent distinct timescales, we apply targeted corrections to each component, thereby preserving both short- and long-term structure in the data. Experimental illustrations show that the timescale-aware EMDBC framework matches the performance of conventional quantile-delta mapping (QDM) at the native daily scale and achieves progressively larger bias reductions at bi-weekly, seasonal, and annual scales. As a result, the proposed approach offers a more robust path to accurate and reliable Earth system projections, strengthening their utility for resilience and adaptation planning.

Ganguli, Arkaprabha [Argonne National Laboratory (↗

Ground-Based Localization of Mars Rovers

The document discusses a procedure for localizing the Mars rovers in site frame, a locally defined reference frame on the Martian surface. MER onboard position within a site frame is estimated onboard and is based on wheel odometry. Odometry estimation of rover position is only reliable over relatively short distances assuming no wheel slip, sinkage, etc. As the rover traverses, its onboard estimate of position in the current site frame accumulates errors and will need to be corrected on occasions via relocalization on the ground (mission operations). The procedure provides a systematic process for ground operators to localize the rover. The method focuses on analysis of acquired images used to declare a site frame and images acquired post-drive. Target selection is performed using two main steps. In the first step, the user identifies features of interest from the images used to declare the current site. Each of the selected target s position in site frame is recorded. In the second step, post-traverse measurements of the selected features positions are recorded again, this time in rover frame, using images acquired post-traverse. In the third step, we transform the post-traverse target s positions to local level frame. In the fourth step, we compute the delta differences in the pre- and post-traverse target s position. In the fifth step, we analyze the delta differences with techniques that compute their statistics to determine the rover s position in the site frame.

Trebi-Ollennu, Ashitey↗

An analysis of energy deposition in a tissue equivalent proportional counter onboard the space shuttle

An improved prediction for space radiations in the lower earth orbits measured by the shuttle TEPC is obtained when energy loss straggling and chord length distribution of the detector are considered. A generalized analytic model is used to describe the energy deposition of direct ion interaction events in a micron-size detector. The transport calculation accounting for the shuttle configuration is accomplished by using a new version of HZETRN that has been extensively verified with laboratory and flight data. The agreement of predicted and measured lineal energy spectra is within 70% for the region above 2 keV/micrometer but within a factor of 2.3 underpredicted for the region below this value. The inclusion of indirect delta ray events in the model is needed before possible causes for the underprediction below 2 keV/micrometer can be assessed.

Flight Experiment↗

Where do the fish go in winter? A year of observing seasonal changes in Sequim Bay’s nearshore fish community

The use of environmental DNA (eDNA) sampling has been proposed as a complementary method to monitor fish species in marine environments. eDNA offers a non-invasive, cost-effective, and scalable way to detect aquatic species. It is appealing in environments where traditional methods are limited by access or visibility, especially in complex or sensitive habitats such as tidal channels and other marine energy sites. Before eDNA can be fully relied upon, we must verify its accuracy against established methods, like underwater photography. In this study, we collected eDNA samples and concurrently deployed a 360-degree camera near the floating dock of PNNL-Sequim in the tidal channel of Sequim Bay once a month for twelve consecutive months. During deployments, the camera remained on the seafloor for several hours overlapping slack tide and captured time-lapse photographs at ten second intervals. Counts and identifications of fish species observed in the images were used to calculate monthly Shannon diversity and Pielou evenness indices. These values were compared across the months using a Kruskal-Wallis test paired with a Conover-Iman post-hoc test. Cliff’s Delta was also calculated to quantify the effect size of the monthly differences. No fish were observed from December through April, likely due to seasonal behavior changes within the local fish community. Fish returned to the shoreline in May, with the greatest diversity and evenness recorded in September. These findings reveal substantial seasonal variation in nearshore fish communities. Many monthly comparisons were found to contain statistically significant differences within the diversity and evenness, and even more were found to have large effect sizes, signifying large ecological changes throughout the seasons. The absence of fish observed during the winter months is a key outcome of this survey that will hopefully be reflected in the eDNA results still to come, which would help validate the eDNA approach.

59 BASIC BIOLOGICAL SCIENCES↗

Aging attenuates the vestibulosympathetic reflex in humans

BACKGROUND: The vestibular system contributes to sympathetic activation by engagement of the otolith organs. However, there is a significant loss of vestibular function with aging. Therefore, the purpose of the present study was to determine if young and older individuals differ in their cardiovascular and sympathetic responses to otolithic stimulation (ie, head-down rotation, HDR). We hypothesized that responses to otolithic stimulation would be attenuated in older adults because of morphological and physiological alterations that occur in the vestibular system with aging. METHODS AND RESULTS: Arterial blood pressure, heart rate, muscle sympathetic nerve activity (MSNA), and head rotation were measured during HDR in 11 young (26 +/- 1 years) and 11 older (64 +/- 1 years) subjects in the prone posture. Five older subjects performed head rotation (chin to chest) in the lateral decubitus position, which simulates HDR but does not alter afferent inputs from the vestibular system. MSNA responses to HDR were significantly attenuated in older as compared with young subjects (P<0.01). MSNA increased in the older subjects by only 12 +/- 5% as compared with 85 +/- 16% in the young. Furthermore, HDR elicited significant reductions in mean arterial blood pressure in older (Delta-6 +/- 1 mm Hg; P<0.01) but not young subjects (Delta1 +/- 1 mm Hg). In contrast to HDR, head rotation performed in the lateral decubitus position did not elicit hypotension. MSNA responses to baroreceptor unloading and the cold pressor test were not different between the age groups. CONCLUSIONS: These data indicate that aging attenuates the vestibulosympathetic reflex in humans and may contribute to the increased prevalence of orthostatic hypotension with age.

NASA Discipline Neuroscience↗

Methodology and Results of the Near-Earth Object (NEO) Human Space Flight (HSF) Accessible Targets Study (NHATS)

Near-Earth Asteroids (NEAs) have been identified by the current administration as potential destinations for human explorers during the mid-2020s. While the close proximity of these objects' orbits to Earth's orbit creates a risk of highly damaging or catastrophic impacts, it also makes some of these objects particularly accessible to spacecraft departing Earth, and this presents unique opportunities for solar system science and humanity's first ventures beyond cislunar space. Planning such ambitious missions first requires the selection of potentially accessible targets from the growing population of nearly 7,800 NEAs. To accomplish this, NASA is conducting the Near-Earth Object (NEO) Human Space Flight (HSF) Accessible Targets Study (NHATS). Phase I of the NHATS was executed during September of 2010, and Phase II was completed by early March of 2011. The study is ongoing because previously undetected NEAs are being discovered constantly, which has motivated an effort to automate the analysis algorithms in order to provide continuous monitoring of NEA accessibility. The NHATS analysis process consists of a trajectory filter and a minimum maximum estimated size criterion. The trajectory filter employs the method of embedded trajectory grids to compute all possible ballistic round-trip mission trajectories to every NEA in the Jet Propulsion Laboratory (JPL) Small-Body Database (SBDB) and stores all solutions that satisfy the trajectory filter criteria. An NEA must offer at least one qualifying trajectory solution to pass the trajectory filter. The Phase II NHATS filter criteria were purposely chosen to be highly inclusive, requiring Earth departure date between January 1st, 2015 and December 31st, 2040, total round-trip flight time <= 450 days, stay time at the NEA >= 8 days, Earth departure C(sub 3) energy <= 60 km(exp 2)/s(exp 2), total mission delta-v <= 12 km/s (including an Earth departure maneuver from a 400 km altitude circular parking orbit), and a maximum atmospheric re-entry speed of 12 km/s. After determining which NEAs offer at least one trajectory solution meeting the criteria, the estimated size constraint is then imposed whereby those NEAs may only be considered NHATS- qualifying NEAs if their maximum estimated size is >= 30 m. This corresponds to an absolute magnitude H <= 26.5 with an assumed albedo p = 0:05. The following is a brief high-level summary of the Phase II study results. Of the 7,665 NEAs in the SBDB as of February 3rd, 2011, 765 NEAs passed the trajectory filter and yielded a total of 79,157,604 trajectory solutions. The trajectory solutions for each NEA are post-processed into Pork Chop Contour (PCC) plots which show total mission delta-v as a function of Earth departure date and total mission duration. Although the PCC plots necessarily compress a very multi-dimensional design space into a two-dimensional plot, they permit rapid assessment of the breadth and quality of an NEA's available Earth departure season and clearly indicate the regions of the trajectory design space which warrant further analysis and optimization. The PCC plot for the NEA with the greatest number of NHATS-qualifying trajectory solutions, 2000 SG-344, is shown. Of the 765 NEAs which passed the Phase II trajectory filter, a total of 590 NEAs also satisfied the further constraint of maximum estimated size >= 30 m. The distributions of osculating heliocentric orbital semi-major axis (a), eccentricity (e), and inclination (i), for those 590 NEAs are shown. Note that the semi-latus rectum used is equal to alpha (1-e(exp 2)). To further our understanding of round-trip trajectory accessibility dynamics, it is instructive to examine the distribution of the NHATS-Qualifying NEAs according to orbit classification. NEAs are grouped into four orbit families: Atiras (aphelion < 0.983 AU), Atens (aphelion > 0.983 AU, alpha < 1.0 AU), Apollos (perihelion < 1.017 AU, alpha > 1.0 AU), and Amors (1.017 < perihelion < 1.3 AU). Of the 765 NEAhich satisfied the NHATS trajectory criteria, none are Atiras, 193 are Atens (31% of known Atens), 456 are Apollos (11% of known Apollos), and 116 are Amors (4% of known Amors). While Apollos comprise 60% of the NEAs which pass the NHATS trajectory filter and Atens comprise only 25%, the percentages according to orbit family are perhaps more relevant. Note that only 11% of known Apollos passed the trajectory filter while 31% of known Atens passed. These simple statistics alone strongly suggest that Aten orbits possess features which tend to enhance their round-trip trajectory accessibility as compared to Apollos or Amors. This is significant because Atens' orbits cause them to spend considerable time in Earth's daytime sky, making them difficult to discover and track using ground-based observing assets. In this paper we will detail the NHATS analysis algorithms, present and analyze all NHATS results to date, and discuss aspects of HSF mission architecture design for future NEA missions.

Barbee, Brent↗

Modelling of SAR polarisation phase difference from trees

The data for polarization phase difference Delta Phi between the HH- and VV-polarized backscattered waves from tree-covered fields were obtained with an airborne synthetic aperture radar at 1.225 GHz. The mean values over tree-covered fields were derived from the images of the phase difference and were examined as a function of incident beam angle from 15 to 55 deg. A theoretical model for simulating these data, based on the electromagnetic wave scatterings from the tree trunk and its branches, both of which are assumed as very long dielectric cylinders was developed. The radius and direction of a tree branch are taken as random variables and are chosen by a Monte Carlo method to encounter the incident waves in producing the scattering events. The Monte Carlo simulated results are in good agreement with the observations within experimental uncertainty.

Mo, Tsan↗

Lee-side flow over delta wings at supersonic speeds

An experimental investigation of the lee-side flow on sharp leading-edge delta wings at supersonic speeds has been conducted. Pressure data were obtained at Mach numbers from 1.5 to 2.8, and three types of flow-visualization data (oil-flow, tuft, and vapor-screen) were obtained at Mach numbers from 1.7 to 2.8 for wing leading-edge sweep angles from 52.5 deg to 75 deg. From the flow-visualization data, the lee-side flows were classified into seven distinct types and a chart was developed that defines the flow mechanism as a function of the conditions normal to the wing leading edge, specifically, angle of attack and Mach number. Pressure data obtained experimentally and by a semiempirical prediction method were employed to investigate the effects of angle of attack, leading-edge sweep, and Mach number on vortex strength and vortex position. In general, the predicted and measured values of vortex-induced normal force and vortex position obtained from experimental data have the same trends with angle of attack, Mach number, and leading-edge sweep; however, the vortex-induced normal force is underpredicted by 15 to 30 percent, and the vortex spanwise location is overpredicted by approximately 15 percent.

Miller, D. S.↗

Optimization of active control systems to suppress flutter and minimize turbulence response

A method for optimizing an active control system which will suppress flutter and minimize the response of a lifting surface to atmospheric turbulence is presented. A mathematical search method is developed which will find a control law which will cause an active control system to flutter at a specified freestream velocity, air density, and Mach number. With the flutter velocity of the system held constant, the control law is then modified in such a way that the peak output power spectral density function of the angular response of a lifting surface (as a result of atmospheric turbulence) is minimized for a specified flight velocity which is less than the flutter velocity. The von Karman generalized power spectrum for the transverse components of turbulence is used in an example problem to increase the flutter velocity and minimize the turbulence response of a simplified delta-wing model which has leading and trailing edge control surfaces.

Rudisill, C. S.↗

Analytical solution for the drift behavior of Dynamics Explorer-B

The Dynamics Explorer (DE) mission was designed to explore the earth's magnetosphere, ionosphere, and plasmasphere. The DE mission employs two spacecraft, including DE-A and DE-B. The spacecraft were launched from the Western Test Range on Aug. 3, 1981, onboard the same Delta launch vehicle. The two spacecraft were placed in a coplanar polar orbit, with DE-A in a high altitude orbit (perigee = 570 km and apogee = 23174 km) and DE-B in a low altitude orbit (perigee = 309 km and apogee = 1013 km). In the present investigation, an attempt has been made to validate the flight data for DE-B with analytic and simulation results. The external torques acting on the spacecraft are represented in terms of tractable mathematical functions. A piecewise linear model of the superrotation of the upper atmosphere is assumed. The effect of individual torques on the long-term pitch axis motion is investigated using analytic and simulation methods. The results are found to be in very good agreement with the available flight data.

Sellappan, R. G.↗

Noninvasive estimation of transmitral pressure drop across the normal mitral valve in humans: importance of convective and inertial forces during left ventricular filling

OBJECTIVES: We hypothesized that color M-mode (CMM) images could be used to solve the Euler equation, yielding regional pressure gradients along the scanline, which could then be integrated to yield the unsteady Bernoulli equation and estimate noninvasively both the convective and inertial components of the transmitral pressure difference. BACKGROUND: Pulsed and continuous wave Doppler velocity measurements are routinely used clinically to assess severity of stenotic and regurgitant valves. However, only the convective component of the pressure gradient is measured, thereby neglecting the contribution of inertial forces, which may be significant, particularly for nonstenotic valves. Color M-mode provides a spatiotemporal representation of flow across the mitral valve. METHODS: In eight patients undergoing coronary artery bypass grafting, high-fidelity left atrial and ventricular pressure measurements were obtained synchronously with transmitral CMM digital recordings. The instantaneous diastolic transmitral pressure difference was computed from the M-mode spatiotemporal velocity distribution using the unsteady flow form of the Bernoulli equation and was compared to the catheter measurements. RESULTS: From 56 beats in 16 hemodynamic stages, inclusion of the inertial term ([deltapI]max = 1.78+/-1.30 mm Hg) in the noninvasive pressure difference calculation significantly increased the temporal correlation with catheter-based measurement (r = 0.35+/-0.24 vs. 0.81+/-0.15, p< 0.0001). It also allowed an accurate approximation of the peak pressure difference ([deltapc+I]max = 0.95 [delta(p)cathh]max + 0.24, r = 0.96, p<0.001, error = 0.08+/-0.54 mm Hg). CONCLUSIONS: Inertial forces are significant components of the maximal pressure drop across the normal mitral valve. These can be accurately estimated noninvasively using CMM recordings of transmitral flow, which should improve the understanding of diastolic filling and function of the heart.

Non-NASA Center↗

Spatially Resolved, In Situ Carbon Isotope Analysis of Archean Organic Matter

Spatiotemporal variability in the carbon isotope composition of sedimentary organic matter (OM) preserves information about the evolution of the biosphere and of the exogenic carbon cycle as a whole. Primary compositions, and imprints of the post-depositional processes that obscure them, exist at the scale of individual sedimentary grains (mm to micron). Secondary ion mass spectrometry (SIMS) (1) enables analysis at these scales and in petrographic context, (2) permits morphological and compositional characterization of the analyte and associated minerals prior to isotopic analysis, and (3) reveals patterns of variability homogenized by bulk techniques. Here we present new methods for in situ organic carbon isotope analysis with sub-permil precision and spatial resolution to 1 micron using SIMS, as well as new data acquired from a suite of Archean rocks. Three analytical protocols were developed for the CAMECA ims1280 at WiscSIMS to analyze domains of varying size and carbon concentration. Average reproducibility (at 2SD) using a 6 micron spot size with two Faraday cup detectors was 0.4 %, and 0.8 % for analyses using 1 micron and 3 micron spot sizes with a Faraday cup (for C-12) and an electron multiplier (for C-13). Eight coals, two ambers, a shungite, and a graphite were evaluated for micron-scale isotopic heterogeneity, and LCNN anthracite (delta C-13 = -23.56 +/- 0.1 %, 2SD) was chosen as the working standard. Correlation between instrumental bias and H/C was observed and calibrated for each analytical session using organic materials with H/C between 0.1 and 1.5 (atomic), allowing a correction based upon a C-13H/C-13 measurement included in every analysis. Matrix effects of variable C/SiO2 were evaluated by measuring mm to sub-micron graphite domains in quartzite from Bogala mine, Sri Lanka. Apparent instrumental bias and C-12 count rate are correlated in this case, but this may be related to a crystal orientation effect in graphite. Analyses of amorphous Archean OM suggest that instrumental bias is consistent for 12C count rates as low as 10% relative to anthracite. Samples from the ABDP-9 (n=3; Mount McRae Shale, approximately 2.5 Ga), RHDH2a (n=2; Carrawine Dolomite and Jeerinah Fm, approximately 2.6 Ga), WRL1 (n=3; Wittenoom Fm, Marra Mamba Iron Formation, and Jeerinah Fm, approximately 2.6 Ga), and SV1 (n=1; Tumbiana Fm, approximately 2.7 Ga) drill cores, each previously analyzed for bulk organic carbon isotope composition, yielded 100 new, in situ data from Neoarchean sedimentary OM. In these samples, delta C-13 varies between -53.1 and -28.3 % and offsets between in situ and bulk compositions range from -8.3 to 18.8%. In some cases, isotopic composition and mode of occurrence (e.g. morphology and mineral associations) are statistically correlated, enabling the identification of distinct reservoirs of OM. Our results support previous evidence for gradients of oxidation with depth in Neoarchean environments driven by photosynthesis and methane metabolism. The relevance of these findings to questions of bio- and syngenicity as well as the alteration history of previously reported Archean OM will be discussed.

Williford, Kenneth H.↗

Theoretical studies on flapped delta wings

The effects of leading edge flaps on the aerodynamic characteristics of a low aspect-ratio delta wing are studied theoretically. As an extension of the classical crossflow plane analysis and in order to include separated shear layers, an analogy between three dimensional steady conical and two dimensional unsteady self-similar flows is explored. This analogy provides a simple steady-unsteady relationship. The criteria for the validity of the steady-unsteady analogy are also examined. Two different theoretical techniques are used to represent the separated shear layers based on the steady-unsteady analogy, neglecting the trailing edge effect. In the first approach, each vortex system is represented by a pair of concentrated vortices connected to the separation points by straight feeding sheets. In the second approach, the vortex cloud method is adopted for simulating the flow field in the crossflow plane. The separated shear layers are replaced with a cloud of discrete vortices and the boundary element method is employed to represent the wing trace by a vorticity distribution. A simple merging scheme is used to model the core region of the vortical flow as a single vortex by imposing a restriction on the shear layer rotation angle. The results are compared with experiments and with results from 3-D panel calculations.

Oh, S.↗

Examination of a Practical Aerobraking Guidance Algorithm

A practical real time guidance algorithm has been developed for aerobraking vehicles that minimizes the post-aeropass Delta V requirements for orbit insertion while nearly minimizing the maximum heating rate and the maximum structural loads. The algorithm is general in the sense that a minimum of assumptions is made, thus greatly reducing the number of parameters that must be determined prior to a given mission. An interesting feature is that in-plane guidance performance is tuned by adjusting one mission-dependent parameter, the bank margin; similarly, the out-of-plane guidance performance is tuned by adjusting a plane controller time constant. Other features of the algorithm are simplicity, efficiency, and ease of use. The algorithm is designed for, but not restricted to, a trimmed vehicle with bank angle modulation as the method of trajectory control. Performance of this guidance algorithm during flight in Earth's atmosphere is examined by its use in an aerobraking testbed program. The performance inquiry extends to a wide range of entry speeds covering a number of potential mission applications. Favorable results have been obtained with a minimum of development effort, and directions for improvement of performance are indicated.

Evans, Steven W.↗

Revised (Mixed-Effects) Estimation for Forest Burning Emissions of Gases and Smoke, Fire/Emission Factor Typology, and Potential Remote Sensing Classification of Types for Ozone and Black-Carbon Simulation

We summarize recent progress (a) in correcting biomass burning emissions factors deduced from airborne sampling of forest fire plumes, (b) in understanding the variability in reactivity of the fresh plumes sampled in ARCTAS (2008), DC3 (2012), and SEAC4RS (2013) airborne missions, and (c) in a consequent search for remotely sensed quantities that help classify forest-fire plumes. Particle properties, chemical speciation, and smoke radiative properties are related and mutually informative, as pictures below suggest (slopes of lines of same color are similar). (a) Mixed-effects (random-effects) statistical modeling provides estimates of both emission factors and a reasonable description of carbon-burned simultaneously. Different fire plumes will have very different contributions to volatile organic carbon reactivity; this may help explain differences of free NOx(both gas- and particle-phase), and also of ozone production, that have been noted for forest-fire plumes in California. Our evaluations check or correct emission factors based on sequential measurements (e.g., the Normalized Ratio Enhancement and similar methods). We stress the dangers of methods relying on emission-ratios to CO. (b) This work confirms and extends many reports of great situational variability in emissions factors. VOCs vary in OH reactivity and NOx-binding. Reasons for variability are not only fuel composition, fuel condition, etc., but are confused somewhat by rapid transformation and mixing of emissions. We use "unmixing" (distinct from mixed-effects) statistics and compare briefly to approaches like neural nets. We focus on one particularly intense fire the notorious Yosemite Rim Fire of 2013. In some samples, NOx activity was not so suppressed by binding into nitrates as in other fires. While our fire-typing is evolving and subject to debate, the carbon-burned delta(CO2+CO) estimates that arise from mixed effects models, free of confusion by background-CO2 variation, should provide a solid base for discussion. (c) We report progress using promising links we find between emissions-related "fire types" and promising features deducible from remote observations of plumes, e.g., single scatter albedo, Angstrom exponent of scattering, Angstrom exponent of absorption, (CO column density)/(aerosol optical depth).

Remote sensing↗

Comparison between the magnetic and transport critical current densities in high critical current density melt-textured yttrium barium copper-oxide

Using a recently developed pulsed critical current density (Jc) measuring system, the Jc of the high-Jc melt-textured YBa2Cu3O(7-delta) (Y123) bulk samples has been determined. I-V curves with a voltage resolution of 0.5 microV were obtained, and transport Jc's along the a-b plane as high as 7.2 x 10 to the 4th A/sq cm were extracted. These results are comparable to the values obtained magnetically. On the other hand, transport Jc along the c axis were found to be two orders of magnitude smaller, even though the magnetic Jc along the c axis is only about five times smaller than Jc along the a-b plane. It is suggested that for the high-temperature superconducting materials which are highly anisotropic, caution should be taken when using the nontransport magnetic methods to determine Jc.

Gao, L.↗

Aerodynamic Characteristics in Pitch of a Series of Cruciform-Wing Missiles with Canard Controls at a Mach Number of 2.01

An investigation has been conducted in the Langley 4- by 4-foot supersonic pressure tunnel at a Mach number of 2.01 to determine the static longitudinal stability and control characteristics of a series of missile configurations with canard controls at angles of attack up to about 28 deg. The missiles had cruciform wings and canard surfaces of delta plan form with 70 deg.swept leading edges. Five bodies having fineness ratios of 19.1, 17.7, 16.7, 15.7, and 14.8 were investigated. The results of the investigation indicated a large nonlinear variation of pitching moment with angle of attack for the body of largest fineness ratio that was progressively reduced by decreasing the fineness ratio until it was essentially eliminated for a body of fineness ratio 14.8. The increased linearity of the moment curve would make it possible to reduce the margin of stability so that, for a given canard size and deflection, a higher trim angle of attack might be obtained for the shortest missile than for the longest missile. The pitching-moment results indicated that methods of prediction which assumed linear variations with angle of attack for the wing-alone and wing-plus-interference characteristics were adequate for angles of attack up to about 12 deg. At higher angles of attack it was evident that the characteristics of these components were nonlinear and that more refined methods would be required for adequate prediction.

Spearman, M. Leroy↗

The Impact of Soil Reflectance on the Quantification of the Green Vegetation Fraction from NDVI

The green vegetation fraction (Fg) is an important climate and hydrologic model parameter. A common method to calculate Fg is to create a simple linear mixing rnodeP between two NDVI endmembers: bare soil NDVI (NDVI(sub o)) and full vegetation NDVI (NDVI(sub infinity)). Usually it is assumed that NDVI(sub o), is close to zero (NDVI(sub o) approx.-0.05) and is generally chosen from the lowest observed NDVI values. However, the mean soil NDVI computed from 2906 samples is much larger (NDVI=0.2) and is highly variable (standard deviation=O. 1). We show that the underestimation of NDVI(sub o) yields overestimations of Fg. The largest errors occur in grassland and shrubland areas. Using parameters for NDVI(sub o) and NDVI(sub infinity) derived from global scenes yields overestimations of Fg ((Delta) Fg*) that are larger than 0.2 for the majority of U.S. land cover types when pixel NDVI values are 0.2<NDVI(sub pixel)<0.4. Figure 1 shows how the Fg overestimation varies for the most common land cover types in the conterminous U.S. for typical seasonal NDVI values. When using conterminous U.S. scenes to derive NDV(sub o) and NDVI(sub infinity), the overestimation is less (0.10-0.17 for 0.2<NDVI(sub pixel)<0.4). As a result, parts of the conterminous U.S. are affected at different times of the year depending on the local seasonal NDVI cycle. We propose using global databases of NDVI(sub o) along with information on historical NDVI(sub pixel) values to compute a statistically most-likely estimate of Fg (Fg*). Using in situ measurements made at the Sevilleta LTER, we show that this approach yields better estimates of Fg than using global invariant NDVI(sub o) values estimated from whole scenes (Figure 2). At the two studied sites, the Fg estimate was adjusted by 52% at the grassland and 86% at the shrubland. More significant advances will require information on spatial distribution of soil reflectance.

Montandon, L. M.↗