Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 289 records · Page 16

A description of the solar wind-magnetosphere coupling based on nonlinear filters

A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Vassiliadis, D.↗

U.S. Coast Guard Boat Recovery Simulation at the NASA Ames Vertical Motion Simulator

The U.S. Coast Guard routinely uses the Over-The-Horizon (OTH-IV), a cutter deployed, rigid-hulled inflatable boat for rescue and law enforcement operations. The term "cutter" refers to a Coast Guard vessel 65 feet in length or greater with accommodations for crew to live aboard and the ability to deploy smaller boats including the OTH-IV. A deployment method employed by some cutter class vessels is a davit crane system and specialized hook mechanism to connect the smaller boat. The launch and recovery operation requires several crewmembers, where one is tasked with releasing and connecting the hook to the small boat. Manipulating the hook can pose a risk to the crewmember depending on sea conditions combined with fatigue level or task expertise. Equipment and crewmembers are tested extensively due to the inherent risk of the task, but the sea conditions cannot be controlled. To address the sea condition variable, the U.S. Coast Guard and NASA leveraged the Vertical Motion Simulator (VMS) as a platform to test new equipment or train crewmembers in varying sea conditions. With 6 degrees of freedom and a vertical displacement of +/- 22ft, the VMS is capable of simulating conditions up to sea-state 5. The proof of concept experiment took place in September 2018 and demonstrated that the VMS can accurately simulate varying sea-states, collect performance data, and design a reliable and safe system for participants. The U.S. Coast Guard provided boat displacement data from the Large Amplitude Motion Program (LAMP). The program generates multi-directional waves with a cosine squared spreading function which produced a time history response of displacement data (surge, sway, heave, roll, pitch, yaw) of the OTH-IV's center of gravity (C.G). The displacement data was differentiated into acceleration evaluated at the crewmember's position at the front (bow) of the boat. Those accelerations were transformed to the VMS coordinate system and used as command inputs to the motion system's washout filters. The VMS uses Interchangeable Cabs (ICab) to perform flight simulations on a variety of aircraft or spacecraft configurations. An ICab typically contains seats, restraints, pilot controls, instrument panel, and out-the-window projectors all enclosed in a familiar cockpit configuration. To replicate the launch and recovery of a small boat, all the conventional features of a flight simulator and canopy were removed to be replaced with a to-scale OTH-IV bow mockup. The bow section was constructed from aluminum, and a connection point known as the davit ring. Additionally, a hook replica was suspended by a cable from the ceiling as if it were hanging from a davit crane aboard a cutter vessel. The last step was providing an auxiliary control to close the distance between suspended hook and the davit ring on the bow mockup. Without control the boat would simulate the wave motion, but with no guarantee that the crewmember will be in reach of the hook. Since the wave motion was repeatable, an additional acceleration command sent to the washout filters to move the simulator in the surge, sway, and heave directions, independently. The commands were rate limited such that the additional movement was not detectable with respect to the wave motion. The results showed that for each wave profile, a unique set of acceleration commands provided opportunities for crewmembers to gain control of the hook and connect it to the davit ring.

Riccobono, Nicholas S.↗

On the passage of radiation through inhomogeneous, moving media. XII - Polarized waves in a plane, sheared medium

Using the complete set of Maxwell equations, the ray paths and Stokes parameters of a polarized wave are followed which passes through a slab of differentially shearing material of variable refractive index. The polarization state of the signal leaving the top of the slab is related to the polarization state of the initial wave. In general, the two are not the same. These calculations were made to illustrate the point that the polarization properties of observed pulsar signals may be very different from those of the emitted signals due to propagation effects in the differentially shearing magnetospheres that most pulsars are believed to possess. While it is recognized that the precise amount of polarization variations is model dependent, the calculations show that even for simple situations a noticeable polarization variation occurs. Accordingly, the calculations reported here are an educative device: they show that it is, perhaps, unwise to assume that the observed polarization properties of pulsar signals are the same as the emitted signals.

Lee, M. A.↗

Mission Simulation of Space Lidar Measurements for Seasonal and Regional CO2 Variations

Results of mission simulation studies are presented for a laser-based atmospheric [82 sounder. The simulations are based on real-time carbon cycle process modeling and data analysis. The mission concept corresponds to the Active Sensing of [82 over Nights, Days, and Seasons (ASCENDS) recommended by the US National Academy of Sciences Decadal Survey of Earth Science and Applications from Space. One prerequisite for meaningful quantitative sensor evaluation is realistic CO2 process modeling across a wide range of scales, i.e., does the model have representative spatial and temporal gradients? Examples of model comparison with data will be shown. Another requirement is a relatively complete description of the atmospheric and surface state, which we have obtained from meteorological data assimilation and satellite measurements from MODIS and [ALIPS0. We use radiative transfer model calculations, an instrument model with representative errors ' and a simple retrieval approach to complete the cycle from "nature" run to "pseudo-data" CO2, Several mission and instrument configuration options are examined/ and the sensitivity to key design variables is shown. We use the simulation framework to demonstrate that within reasonable technological assumptions for the system performance, relatively high measurement precision can be obtained, but errors depend strongly on environmental conditions as well as instrument specifications. Examples are also shown of how the resulting pseudo - measurements might be used to address key carbon cycle science questions.

Kawa, Stephan↗

Detecting Risk and Anomalies in Airplane Dynamics Through Entropic Analysis of Time Series Data

Despite recent efforts to move away from traditional threshold exceedance detection methods for aircraft state monitoring, modern aircraft still rely on safety thresholds to communicate to pilots the identification of an anomaly in the aircraft when a threshold is surpassed. Current anomaly detection methods mainly depend on uninterpretable machine learning models to learn complex patterns and relationships contained in the time series data of aircraft. Although these methods are capable of identifying known anomalies, their deficiency in interpretability presents a challenge when translating them to different aircraft. To overcome this deficiency, entropic analysis of aircraft dynamics seeks to characterize the complexity, or lack thereof, of the aircraft dynamics prior to the development of a risk scenario. This complexity characterization provides a more straightforward summary of state changes in the dynamics of flight variables. To build a foundation for entropic analysis, we analyzed the complexity of unstable approaches, an anomalous event present in many of today’s aviation accidents. The analysis revealed a statistically significant difference in the complexity distribution of flight variables under a stable approach versus an unstable approach. These differences in complexity were especially notable minutes before an approach was identified as unstable. Moreover, the multiscale entropic analysis revealed the presence of signal complexity at multiple time scales across multiple time windows before landing. By capturing state changes and corrections in the aircraft dynamics using entropy, advanced, yet still interpretable, sensor systems based on entropic frameworks from this study can be constructed in the future using classical machine learning approaches.

Risk detection↗

Radar determination of Mars surface properties

Radar studies of Mars have provided measurements of surface texture on scales of centimeters to hundreds of meters and measurements of surface material properties. Texture (rms surface tilts and estimates of small-scale roughness) may be inferred from dispersion and/or polarization of the radar echo; material properties (reflectivity of dielectric constant) are derived from echo strength. Mars is a diverse target; depending on location, rms surface tilts have been found to vary over the range 0.25 to 10 deg while reflectivity covers at least 3 to 13 percent. Plains units are the most variable, having both the smoothest and roughest surfaces, cratered terrain can be considered predictable and 'average' by comparison. Recent data identify scattering by small structures (perhaps rocks on or near the surface) as playing a more important role than previously recognized. Scattering by the residual ice cap near Mars' south pole is particularly unusual. The present state of radar surface studies is summarized.

Simpson, Richard A.↗

A rapid implicit-explicit solution to the two-dimensional time dependent incompressible Navier-Stokes equations

A second-order time-accurate and spatially factored algorithm was used in a finite difference scheme for the numerical solution of the time-dependent, incompressible, two dimensional Navier-Stokes equations in conservation-law form using vorticity and stream function variables. The systems of equations are solved at each time step by an iterative technique. Numerical results were obtained for a circular cylinder at a Reynolds number of 15, and an NACA 0012 airfoil at zero angle of attack at Reynolds numbers of 10 to the third and 10 to the fourth powers. The results are in agreement with another numerical technique, and the computing time required to obtain the steady state solution at the Reynolds number of 10 to the 4th power was 49.7 sec on CDC 7600 computer using a 65 x 84 computational grind.

Davis, J. E.↗

Magnetic susceptibility of Inconel alloys 718, 625, and 600 at cryogenic temperatures

After a hydrogen fuel bleed valve problem on the Discovery Space Shuttle was traced to the strong magnetization of Inconel 718 in the armature of the linear variable differential transformer near liquid hydrogen temperatures, the ac magnetic susceptibility of three samples of Inconel 718 of slightly different compositions, one sample of Inconel 625, and on sample of Inconel 600 were measured as a function of temperature. Inconel 718 alloys are found to exhibit a spin glass state below 16 K. Inconel 600 exhibits three different magnetic phases, the lowest-temperature state (below 6 K) being somewhat similar to that of Inconel 718. The magnetic states of the Inconel alloys and their magnetic susceptibilities appear to be strongly dependent on the exact composition of the alloy.

Goldberg, Ira B.↗

Arctic Tropospheric Ozone: Assessment of Current Knowledge and Model Performance

As the third most important greenhouse gas (GHG) after carbon dioxide (CO2) and methane (CH4), tropospheric ozone (O3) is also an air pollutant causing damage to human health and ecosystems. This study brings together recent research on observations and modeling of tropospheric O3 in the Arctic, a rapidly warming and sensitive environment. At different locations in the Arctic, the observed surface O3 seasonal cycles are quite different. Coastal Arctic locations, for example, have a minimum in the springtime due to O3 depletion events resulting from surface bromine chemistry. In contrast, other Arctic locations have a maximum in the spring. The 12 state-of-the-art models used in this study lack the surface halogen chemistry needed to simulate coastal Arctic surface O3 depletion in the springtime; however, the multi-model median (MMM) has accurate seasonal cycles at non-coastal Arctic locations. There is a large amount of variability among models, which has been previously reported, and we show that there continues to be no convergence among models or improved accuracy in simulating tropospheric O3 and its precursor species. The MMM underestimates Arctic surface O3 by 5 % to 15 % depending on the location. The vertical distribution of tropospheric O3 is studied from recent ozonesonde measurements and the models. The models are highly variable, simulating free-tropospheric O3 within a range of ±50 % depending on the model and the altitude. The MMM performs best, within ±8 % for most locations and seasons. However, nearly all models overestimate O3 near the tropopause (∼300 hPa or ∼8 km), likely due to ongoing issues with underestimating the altitude of the tropopause and excessive downward transport of stratospheric O3 at high latitudes. For example, the MMM is biased high by about 20 % at Eureka. Observed and simulated O3 precursors (CO, NOx, and reservoir PAN) are evaluated throughout the troposphere. Models underestimate wintertime CO everywhere, likely due to a combination of underestimating CO emissions and possibly overestimating OH. Throughout the vertical profile (compared to aircraft measurements), the MMM underestimates both CO and NOx but overestimates PAN. Perhaps as a result of competing deficiencies, the MMM O3 matches the observed O3 reasonably well. Our findings suggest that despite model updates over the last decade, model results are as highly variable as ever and have not increased in accuracy for representing Arctic tropospheric O3.

Arctic↗

Influence of exit impedance on finite difference solutions of transient acoustic mode propagation in ducts

The time-dependent governing acoustic-difference equations and boundary conditions are developed and solved for sound propagation in an axisymmetric (cylindrical) hard-wall duct without flow and with spinning acoustic modes. The analysis begins with a harmonic sound source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the steady solutions of the acoustic field. The transient method did not converge to the steady-state solution for cutoff acoustic duct modes. This has implications as to its use in a variable-area duct, where modes may become cutoff in the smal-area portion of the duct. For single cutoff mode propagation the steady-state impedance boundary condition produced acoustic reflections during the initial transient that caused finite instabilities in the numerical calculations. The stability problem is resolved by reformulating the exit boundary condition. Example calculations show good agreement with exact analytical and numerical results for forcing frequencies above, below, and nearly at the cutoff frequency.

Baumeister, K. J.↗

Variability in Comet P/Swift-Tuttle

Spatial profiles of the coma of Comet P/Swift-Tuttle perpendicular to the projected Sun-comet line were obtained for the emission bands of CN, C2, and C3 as well as for two continuum bands from spectrophotometric observations taken from Oct. 5 to 8, 1992. The intensities were converted into emissivities per sq km in the cases of the continua and into column densities for the emission band profiles. Spatial and temporal variabilities have been found in all five investigated components, which are consistent with the rotational period of the nucleus determined from the observations of the comet during its last perihelion passage in 1862. The emission band profiles were fitted with the vectorial model and the production rates of CN, C2, and C3 were determined. One half of the profiles was fitted adequately under steady state conditions for the production rate in a first approximation, whereas the other half showed prominent bumps, which could only be explained by introducing a time-dependent production rate. Further investigations showed evidence for the presence of gaseous as well as dust jets in the coma, which indicated the presence of at least two active areas on the surface of the nucleus. The projected radial expansion velocities of two different features were determined to be (430 +/- 100) and (460 +/- 100) m/sec, respectively.

Schulz, Rita↗

Predictability of Ocean Heat Content From Electrical Conductance

Ocean heat content (OHC) is a key climate variable that needs to be monitored to know how Earth's energy imbalance is changing, yet observing OHC remains a challenge. The present study examines whether a depth integral of the ocean's electrical conductivity (“conductance”), which may be inferred from both in situ methods and satellite magnetometers over the global ocean, could help monitor OHC. The ocean's electrical conductivity locally depends on temperature, salinity, and pressure, but it is not as well known how the conductance depends on OHC and ocean salt content. By examining the output of an ocean state estimate shown to agree well with observations that have not been assimilated, this study evaluates the fundamental limitations of using perfectly known ocean conductance to predict OHC, rather than the challenges associated with accounting for observational error. It is found that the ocean's conductance and OHC fields are nonlinearly related but nevertheless highly correlated. A statistical framework tends to predict OHC more accurately than ocean salt content from ocean conductance in regions where conductivity is more sensitive to salinity than temperature. The annually (bidecadally) averaged OHC can be predicted from a combination of conductance and depth‐averaged conductivity ocean fields to within nearly 0.1% (1%) error globally and even more accurately in many poorly observed (e.g., ice‐covered) regions. Practical application of this statistical framework to monitor OHC requires examination of the effect of uncertainties in the observed bathymetry and ocean conductance, which vary with application.

Trossman, D. S.↗

Comprehensive Comparison of Airborne and Spaceborne SAR and LiDAR Estimates of Forest Structure in the Tallest Mangrove Forest on Earth

A recent suite of new global-scale satellite sensors and regional-scale airborne campaigns are providing a wealth of remote sensing data capable of dramatically advancing our current understanding of the spatial distribution of forest structure and carbon stocks. However, a baseline for forest stature and biomass estimates has yet to be established for the wide array of available remote sensing products. At present, it remains unclear how the estimates from these sensors compare to one another in terrestrial forests, with a clear dearth of studies in high carbon density mangrove ecosystems. In the tallest mangrove forest on Earth (Pongara National Park, Gabon), we leverage the data collected during the AfriSAR campaign to evaluate 17 state-of-the-art sensor data products across the full range of height and biomass known to exist globally in mangrove forest ecosystems, providing a much-needed baseline for sensor performance. Our major findings are: (Houghton, Hall, Goetz) height estimates are not consistent across products, with opposing trends in relative and absolute errors, highlighting the need for an adaptive approach to constraining height estimates (Panet al., 2011); radar height estimates had the lowest calibration error and bias, with further improvements using LiDAR fusion (Bonan, 2008); biomass variability and uncertainty strongly depends on forest stature, with variation across products increasing with canopy height, while relative biomass variation was highest in low-stature stands (Le Quéréet al., 2017); a remote sensing product's sensitivity to variations in canopy structure is more important than the absolute accuracy of height estimates (Mitchardet al., 2014); locally-calibrated area-wide totals are more representative than generalized global biomass models for high-precision biomass estimates. The findings presented here provide critical baseline expectations for height and biomass predictions across the full range of mangrove forest stature, which can be directly applied to current (TanDEM-X, GEDI, ICESat-2) and future (NISAR, BIOMASS) global-scale forest monitoring missions.

Atticus E.L. Stovall↗

Impacts of Interannual Climate Variability on Agricultural and Marine Ecosystems

The El Nino - Southern Oscillation (ENSO) is the dominant mode of global interannual climate variability, and seems to be the only mode for which current prediction methods are more skillful than climatology or persistence. The Zebiak and Cane intermediate coupled ocean-atmosphere model has been in use for ENSO prediction for more than a decade, with notable success. However, the sole dependence of its original initialization scheme and the improved initialization on wind fields derived from merchant ship observations proved to be a liability during 1997/1998 El Nino event: the deficiencies of wind observations prevented the oceanic component of the model from reaching the realistic state during the year prior to the event, and the forecast failed. Our work on the project was concentrated on the use of satellite data for improving various stages of ENSO prediction technology: model initialization, bias correction, and data assimilation. Close collaboration with other teams of the IDS project was maintained throughout.

Cane, M. A.↗

Increasing Accuracy in Computed Inviscid Boundary Conditions

A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.

Dyson, Roger↗

Modeling Whitecaps on Global Scale

Whitecaps play an important role in the surface-atmosphere interactions across the ocean. They are directly linked to the energy dissipation rate during wave breaking and transfer of heat, momentum, and gas/aerosol exchange from the wind to the sea surface. Although the first models of W were dependent only on wind speeds, a large number of diverse models based on wind and sea state which include wave height, wave age, friction velocity, and stability effects have been proposed since then. However, it is recognized that most of the proposed W models have strong systematic (e.g., zonal bias) and random errors when compared against observations. This is partly due to the differences in environmental conditions, measurement techniques, and geographical locations among these studies. But, some of these biases are linked to the inability of the proposed models to capture the variability in W in certain wind/wave regimes. Despite the knowledge of existing biases, W residual relationships from the models with wind and wave fields remain highly uncertain, with residual trends varying between the published studies. Here, we take advantage of the availability of relatively dense observations of W from WindSat microwave satellite retrievals in combination with the University of Miami wave model which was recently incorporated within the NASA GMAO/GEOS system (GEOS-UMWM). We use Windsat W retrievals to assess and constrain the previously published W models and understand the relationships of residuals from models in different wind/wave regimes. We link these unexplained residual variations to additional factors such as swell index, drag coefficient etc and add information to the existing whitecap models. Since Windsat retrievals cover wide range of environmental conditions, it helps to reduce the uncertainties associated with differences in measurement techniques. Regression of wind-wave fields against all Windsat data points (CTL) results in larger residuals for lower wave age and W is overestimated upto ~4% for wave age < 10 and underestimated by upto ~2% as wave age increases. We attest to this bias by considering two approaches. One is to perform regression separately for different stages of wave development such as developing sea, fully developed, and wind sea regimes thereby understanding the sensitivity of regression coefficients to sea state (EXP1). Another is to derive coefficients of W models in EXP1 as a function of additional wind/wave factors such as swell index, drag coefficient, and mean squared slope, deriving more nonlinear W models (EXP2). EXP2 provides reduction in Root Mean Squared Error (RMSE) by 0.1-0.3%. Sea surface drag has a stronger relationship with regression coefficients compared to swell index.These additional factors provide improved parameterizations in different wind and wave age regimes, with smaller unexplained/residual variations in W that has been a major concern in the W community.

Raman, Aishwarya↗

Impacts of Aerosol and Cloud Variability on Satellite-Derived All-skies Direct Aerosol Radiative Effects (DARE) over the North and Southeast Atlantic

The Top-of-Atmosphere (TOA) Shortwave (SW) Direct Aerosol Radiative Effect (DARE) in all-sky conditions (i.e., aerosols in clear skies and aerosols above and below all types of clouds) is the global change in upwelling radiative flux due to aerosols. It is one of the strongest indicators of global climate change due to aerosols. SW DARE at TOA depends on the Earth’s surface albedo, cloud fraction, cloud optical properties, and aerosol optical properties, which are all challenging to accurately characterize from space. The overarching goals of our project are to provide state-of-the-art observational all-sky TOA SW DARE, along with guidance on which aerosol and/or cloud properties are the most important to measure, and at which spatio-temporal scales, for accurate DARE observations. We compute all-sky DARE based on state-of-the-art cloud and aerosol retrieval algorithms from CALIOP (Cloud–Aerosol Lidar with Orthogonal Polarization) and MODIS (Moderate Resolution Imaging spectroradiometer) satellite sensors, as well as aerosol intensive properties from MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, version 2) simulations over three specific regions of the Atlantic Ocean from 2012 to 2016. In this symposium, we present a characterization of the cloud and aerosol optical and physical properties and the observational TOA SW DARE along three individual satellite tracks in the Southeast Atlantic region -- on September 18 and 20, 2016 and August 13, 2017. We then quantify the impact of assuming homogeneous cloud or aerosol fields in space-based TOA SW DARE calculations. And, finally, the resulting space-based all-skies aerosol vertical distribution and DARE calculations are evaluated using remote sensing observations from coincident suborbital flights from the NASA ORACLES (the ObseRvations of Aerosols above CLouds and their intEractionS, ORACLES field campaign) field campaign. The NASA Atmospheric Observing System (AOS) mission addresses the NASA Aerosol, Cloud, Convection and Precipitation (ACCP) designated observables and proposes, as one of its aerosol objectives, to reduce uncertainties in estimates of global mean all-sky SW DARE at TOA. Well characterizing clouds, aerosol vertical distributions, aerosol types and associated all-skies DARE over the Atlantic Ocean will inform the AOS community on where, when, how, and how often the satellite retrievals should be performed to estimate DARE and reduce all-skies DARE uncertainties most accurately. These comprehensive characterizations will also identify the key regions and times when the AOS (or other ACCP-related) suborbital missions should be conducted to evaluate and improve the AOS space-based observations and retrievals.

aerosol↗

Effective Interpolation of Incomplete Satellite-Derived Leaf-Area Index Time Series for the Continental United States

Many earth science modeling applications employ continuous input data fields derived from satellite data. Environmental factors, sensor limitations and algorithmic constraints lead to data products of inherently variable quality. This necessitates interpolation of one form or another in order to produce high quality input fields free of missing data. The present research tests several interpolation techniques as applied to satellite-derived leaf area index, an important quantity in many global climate and ecological models. The study evaluates and applies a variety of interpolation techniques for the Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf-Area Index Product over the time period 2001-2006 for a region containing the conterminous United States. Results indicate that the accuracy of an individual interpolation technique depends upon the underlying land cover. Spatial interpolation provides better results in forested areas, while temporal interpolation performs more effectively over non-forest cover types. Combination of spatial and temporal approaches offers superior interpolative capabilities to any single method, and in fact, generation of continuous data fields requires a hybrid approach such as this.

Jasinski, Michael F.↗