Search NASA⌕ Search

SEARCH · Search NASA

Results for “Localization Error”

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 1,009 records · Page 56

Estimation of Crop Gross Primary Production (GPP): Do Scaled (MODIS) Vegetation Indices Improve Performance? - 2

Satellite remote sensing estimates of Gross Primary Production (GPP) have routinely been made using spectral Vegetation Indices (VIs) over the past two decades. The Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), the green band Wide Dynamic Range Vegetation Index (WDRVIgreen), and the green band Chlorophyll Index (CIgreen) have been employed to estimate GPP under the assumption that GPP is proportional to the product of VI and photosynthetically active radiation (PAR) (where VI is one of four VIs: NDVI, EVI, WDRVIgreen, or CIgreen). However, the empirical regressions between VI*PAR and GPP measured locally at flux towers do not pass through the origin (i.e., the zero X-Y value for regressions). Therefore they are somewhat difficult to interpret and apply. This study investigates (1) what are the scaling factors and offsets (i.e., regression slopes and intercepts) between the fraction of PAR absorbed by chlorophyll of a canopy (fAPARchl) and the VIs, and (2) whether the scaled VIs developed in (1) can eliminate the deficiency and improve the accuracy of GPP estimates. Three AmeriFlux maize and soybean fields were selected for this study, two of which are irrigated and one is rainfed. The four VIs and fAPARchl of the fields were computed with the MODerate resolution Imaging Spectroradiometer (MODIS) satellite images. The GPP estimation performance for the scaled VIs was compared to results obtained with the original VIs and evaluated with standard statistics: the coefficient of determination (R2), the root mean square error (RMSE), and the coefficient of variation (CV). Overall, the scaled EVI obtained the best performance. The performance of the scaled NDVI, EVI and WDRVIgreen was improved across sites, crop types and soil/background wetness conditions. The scaled CIgreen did not improve results, compared to the original CIgreen. The scaled green band indices (WDRVIgreen, CIgreen) did not exhibit superior performance to either the scaled EVI or NDVI in estimating crop daily GPP at these agricultural fields. The scaled VIs are more physiologically meaningful than original un-scaled VIs, but scaling factors and offsets may vary across crop types and surface conditions.

Vegetation Index↗

Odyssey Mars Orbiter - Thirteen Years of On-Orbit Navigation

The Odyssey spacecraft has been in Mars orbit since October 24, 2001 and has nearly completed 61,490 orbits. Navigation operational objectives include the following: Control the local mean solar time for science observations; for most of the mission, this varied from 3:45 pm to 5:20 pm. Currently, an orbit trim maneuver planned for November 10, 2015 will place Odyssey at 6:45 pm/6:45 am at equator crossings in order to observe early morning ground frost, fog and clouds. Initially, Odyssey was late by 42 minutes for an over-flight of the critical seven minutes of Phoenix's entry, descent and landing (EDL). Odyssey was successfully positioned for this over-flight using the Delta V from angular momentum desaturations (AMD). Similar results for the Mars Science Laboratory's EDL and Comet Siding Spring's minimum risk location will be presented. Odyssey has and continues to relay significant quantities of rover data. Navigation successfully models frequent AMD Delta Vs in order to generate accurate sixty-day trajectory predictions; a typical timing error is 25 seconds after 60 days. However, unexpected events, such as safe-mode entries with their larger and more frequent thrusting, severely impact that trajectory accuracy. Impacted trajectories can have timing errors ranging from a few minutes to ten-to-fifteen minutes after sixty-days. Other analyses (briefly stated) include: a) the offset of the orbital ground track pattern after an initial cycle of 30 days or 362 orbits and b) an operations environment of continuous thrusting if/when one of the three remaining reaction wheels fails.

trajectory accuracy↗

A feasibility study of hand kinematics for EVA analysis using magnetic resonance imaging

A new method of analyzing the kinematics of joint motion is developed. Magnetic Resonance Imaging (MRI) offers several distinct advantages. Past methods of studying anatomic joint motion have usually centered on four approaches. These methods are x-ray projection, goniometric linkage analysis, sonic digitization, and landmark measurement of photogrammetry. Of these four, only x-ray is applicable for in vivo studies. The remaining three methods utilize other types of projections of inter-joint measurements, which can cause various types of error. MRI offers accuracy in measurement due to its tomographic nature (as opposed to projection) without the problems associated with x-ray dosage. Once the data acquisition of MR images was complete, the images were processed using a 3D volume rendering workstation. The metacarpalphalangeal (MCP) joint of the left index finger was selected and reconstructed into a three-dimensional graphic display. From the reconstructed volumetric images, measurements of the angles of movement of the applicable bones were obtained and processed by analyzing the screw motion of the MCP joint. Landmark positions were chosen at distinctive locations of the joint at fixed image threshold intensity levels to ensure repeatability. The primarily two dimensional planar motion of this joint was then studied using a method of constructing coordinate systems using three (or more) points. A transformation matrix based on a world coordinate system described the location and orientation of a local target coordinate system. Future research involving volume rendering of MRI data focusing on the internal kinematics of the hand's individual ligaments, cartilage, tendons, etc. will follow. Its findings will show the applicability of MRI to joint kinematics for gaining further knowledge of the hand-glove (power assisted) design for extravehicular activity (EVA).

Dickenson, Rueben D.↗

An Algorithm for Efficient Maximum Likelihood Estimation and Confidence Interval Determination in Nonlinear Estimation Problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The algorithm was developed for airplane parameter estimation problems but is well suited for most nonlinear, multivariable, dynamic systems. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort. MNRES determines the sensitivities with less computational effort than using either a finite-difference method or integrating the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, thus eliminating algorithm reformulation with each new model and providing flexibility to use model equations in any format that is convenient. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. It is observed that the degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. The CR bounds were found to be close to the bounds determined by the search when the degree of nonlinearity was small. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels for the parameter confidence limits. The primary utility of the measure, however, was found to be in predicting the degree of agreement between Cramer-Rao bounds and search estimates.

Murphy, Patrick Charles↗

Comparison of ISCCP and other cloud amounts

A new 8-year global cloud climatology has been produced by the International Satellite Cloud Climatology Project (ISCCP) that provides information every 3 h at 280-km spatial resolution covering the period from July 1983 through June 1991. If cloud detection errors and differences in area sampling are neglected, individual ISCCP cloud amounts agree with individual surface observations to within 15% rms with biases of only a few percent. When measurements of small-scale, broken clouds are isolated in the comparison, the rms differences between satellite and surface cloud amounts are about 25%, similar to the rms difference between ISCCP and Landsat determinations of cloud amount. For broken clouds, the average ISCCP cloud amounts are about 5% smaller than estimated by surface observers (difference between earth cover and sky cover), but about 5% larger than estimated from very high spatial resolution satellite observations (overestimate due to low spatial resolution offset by underestimate due to finite radiance thresholds). Detection errors caused by errors in the ISCCP clear-sky radiances or incorrect radiance threshold magnitudes, are the dominant source of error in monthly average cloud amounts. The ISCCP cloud amounts appear to be too low over land by about 10%, somewhat less in summer and somewhat more in winter, and about right (maybe slightly low) over oceans. In polar regions, ISCCP cloud amounts are probably too low by about 15%-25% in summer and 5%-10% in winter. Comparison of the ISCCP climatology to three other cloud climatologies shows excellent agreement in the geographic distribution and seasonal variation of cloud amounts; there is little agreement about day/night contrasts in cloud amount. Notable results from ISCCP are that the gobal annual mean cloud amount is about 63%, being about 23% higher over oceans than over land, that it varies by less than 1% rms from month to month, and that it has varied by about 4% on a time scale approximately equal to 2-4 years. The magnitude of interannual variations of local (280-km scale) monthly mean cloud amounts is about 7%-9%. Longitudinal contrasts in cloud amount are just as large as latitudinal contrasts. The largest seasonal variation of cloud amount occurs in the tropics, being larger in summer than in winter; the seasonal variation in the middle latitudes has the opposite phase. Polar regions may have little seasonable variability in cloud amount. The ISCCP results show slightly more nightime than daytime cloud amount over oceans and more daytime than nightime cloud amount over land.

Rossow, William B.↗

Comparison of Five System Identification Algorithms for Rotorcraft Higher Harmonic Control

This report presents an analysis and performance comparison of five system identification algorithms. The methods are presented in the context of identifying a frequency-domain transfer matrix for the higher harmonic control (HHC) of helicopter vibration. The five system identification algorithms include three previously proposed methods: (1) the weighted-least- squares-error approach (in moving-block format), (2) the Kalman filter method, and (3) the least-mean-squares (LMS) filter method. In addition there are two new ones: (4) a generalized Kalman filter method and (5) a generalized LMS filter method. The generalized Kalman filter method and the generalized LMS filter method were derived as extensions of the classic methods to permit identification by using more than one measurement per identification cycle. Simulation results are presented for conditions ranging from the ideal case of a stationary transfer matrix and no measurement noise to the more complex cases involving both measurement noise and transfer-matrix variation. Both open-loop identification and closed- loop identification were simulated. Closed-loop mode identification was more challenging than open-loop identification because of the decreasing signal-to-noise ratio as the vibration became reduced. The closed-loop simulation considered both local-model identification, with measured vibration feedback and global-model identification with feedback of the identified uncontrolled vibration. The algorithms were evaluated in terms of their accuracy, stability, convergence properties, computation speeds, and relative ease of implementation.

Jacklin, Stephen A.↗

Phase Retrieval Implementation for the WFIRST Coronagraph Development Testbed

The Wide-Field Infrared Survey Telescope (WFIRST) is a 2.4m diameter space telescope NASA program. The payload will include two instruments, including a coronagraph instrument (CGI). The CGI will be the first active, high-contrast stellar coronagraph in space intended for imaging, discovery, and spectral characterization of exoplanets and debris discs. The CGI designs under development use deformable mirrors (DM) to create a point spread function (PSF) with a dark region around the obscured star object. Each DM consists of an array of actuators, and each actuator causes a local displacement to a thin, mirror-coated, glass face sheet. The actuator displacements are designed to create the PSF. Typically, phase retrieval is used to set the DM actuator displacements to their designed pattern, and then the actuators are further modified by an iterative control loop to correct for system non-idealities. To observe an exoplanet near the star object, the contrast in the annular region must be better than 10^8 (10^9 after post processing) relative to the unobstructed image of the star [1,2]. The coronagraph testbed at Jet Propulsion Laboratory is used to test coronagraph design and control algorithms. The optical system (Fig. 1) is maintained in a vacuum tank so that the PSF is stable for long integration times. Phase retrieval is an important tool for the testbed and serves several purposes, in addition to setting initial DM actuator displacements. Other uses include evaluation of systematic wavefront error introduced by system components, and alignment of the DMs and pupil masks. For example, the location of specific DM actuators relative to the Lyot mask can be determined by displacing the actuator, then measuring the location of the phase disturbance in the pupil plane with phase retrieval. In the next section we further describe the optical train of the testbed and the collection of source and pupil images available as inputs for phase retrieval, as well as one phase retrieval algorithm employed. Section 3 shows example phase retrieval results.

Marx, David↗

Rotation Rate of Saturn's Magnetosphere using CAPS Plasma Measurements

We present the present status of an investigation of the rotation rate of Saturn 's magnetosphere using a 3D velocity moment technique being developed at Goddard which is similar to the 2D version used by Sittler et al. (2005) [1] for SOI and similar to that used by Thomsen et al. (2010). This technique allows one to nearly cover the full energy range of the CAPS IMS from 1 V less than or equal to E/Q less than 50 kV. Since our technique maps the observations into a local inertial frame, it does work during roll manoeuvres. We have made comparisons with Wilson et al. (2008) [2] (2005-358 and 2005-284) who performs a bi-Maxwellian fit to the ion singles data and our results are nearly identical. We will also make comparisons with results by Thomsen et al. (2010) [3]. Our analysis uses ion composition data to weight the non-compositional data, referred to as singles data, to separate H+, H2+ and water group ions (W+) from each other. The ion data set is especially valuable for measuring flow velocities for protons, which are more difficult to derive using singles data within the inner magnetosphere, where the signal is dominated by heavy ions (i.e., proton peak merges with W+ peak as low energy shoulder). Our technique uses a flux function, which is zero in the proper plasma flow frame, to estimate fluid parameter uncertainties. The comparisons investigate the experimental errors and potential for systematic errors in the analyses, including ours. The rolls provide the best data set when it comes to getting 4PI coverage of the plasma but are more susceptible to time aliasing effects. Since our analysis is a velocity moments technique it will work within the inner magnetosphere where pickup ions are important and velocity distributions are non-Maxwellian. So, we will present results inside Enceladus' L shell and determine if mass loading is important. In the future we plan to make comparisons with magnetic field observations, use Saturn ionosphere conductivities as presently known and the field aligned currents necessary for the planet to enforce corotation of the rotating plasma.

Sittler, E.↗

Bayesian Estimation of Precipitation from Satellite Passive Microwave Observations Using Combined Radar-Radiometer Retrievals

Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Grecu, Mircea↗

Galileo EM-2 contributions to the lunar control network

A local control network is being developed using Galileo images that cover the region north of the Apollo area and lie between 10 and 100 degrees east longitude. This network is tied to the Apollo control network and will have a positional accuracy of approximately 500-1500 m. This region has been photographed by Earth-based telescopes and the Mariner 10 and Lunar Orbiter spacecraft, but the Gallileo images are preferred for control because of their superior viewing angles, resolution, and Galileo's geometrically stable sensor. Based on Davies et al., 1987, the potential accuracy of the near-side Apollo network is estimated to be between 50 and 300 m. The laser ranging retroreflector locations at the Apollo 11, 14, and 15 sites and the Lunikhod 2 site have been determined with an accuracy of about 10 m. The locations of the ALSEP transmitters relative to the retroreflectors have also been measured. These known coordinates are used to estimate errors in the Apollo network. The telescopic control is approximately bounded by 75 degrees north and south latitude and 75 degrees east and west longitudes. The accuracy of the telescopic network is thought to be about one to two km. Thus, the Galileo network might improve coordinates in the telescopic area as well as north and east of it.

Davies, M. E.↗

Version 3 of the SMAP Level 4 Soil Moisture Product

The NASA Soil Moisture Active Passive (SMAP) Level 4 Soil Moisture (L4_SM) product provides 3-hourly, 9-km resolution, global estimates of surface (0-5 cm) and root zone (0-100 cm) soil moisture as well as related land surface states and fluxes from 31 March 2015 to present with a latency of 2.5 days. The ensemble-based L4_SM algorithm is a variant of the Goddard Earth Observing System version 5 (GEOS-5) land data assimilation system and ingests SMAP L-band (1.4 GHz) Level 1 brightness temperature observations into the Catchment land surface model. The soil moisture analysis is non-local (spatially distributed), performs downscaling from the 36-km resolution of the observations to that of the model, and respects the relative uncertainties of the modeled and observed brightness temperatures. Prior to assimilation, a climatological rescaling is applied to the assimilated brightness temperatures using a 6 year record of SMOS observations. A new feature in Version 3 of the L4_SM data product is the use of 2 years of SMAP observations for rescaling where SMOS observations are not available because of radio frequency interference, which expands the impact of SMAP observations on the L4_SM estimates into large regions of northern Africa and Asia. This presentation investigates the performance and data assimilation diagnostics of the Version 3 L4_SM data product. The L4_SM soil moisture estimates meet the 0.04 m3m3 (unbiased) RMSE requirement. We further demonstrate that there is little bias in the soil moisture analysis. Finally, we illustrate where the assimilation system overestimates or underestimates the actual errors in the system.

GEOS-↗

GRAIL Refinements to Lunar Seismic Structure

The present ]day internal structure of the Moon provides insight not only into its own formation and evolution, but also that of all rocky planetary bodies. The most direct way to probe a planet fs interior structure is through seismology. As part of the Apollo lunar missions, four seismometers were deployed on the nearside surface of the Moon between the years 1969 and 1972. These instruments operated continuously until 1977, forming the only substantial extraterrestrial seismic data set in existence. These data have been used to constrain various aspects of the seismic velocity and density structure of the Moon. Typical 1-D models recognize a 30-60 km thick crust overlying a nearly constant ]velocity mantle, and extend to a depth of approximately 1000 km, below which the lack of penetrating moonquake ray ]paths precludes the seismic determination of deeper structure. Previously, the lack of observed moonquakes from the far side of the Moon has been used to infer the presence of a highly attenuating (possibly molten) core. Indirect geophysical measurements such as moment of inertia, magnetic induction, lunar laser ranging, and elemental abundances of mare basalts also place varying constraints on core size and state. In combination with seismic studies, these indirect measurements have been used to arrive at a commonly accepted model of the Moon's deepest interior that includes a solid inner and fluid outer core, overlain by a partial melt boundary layer. We recently applied modern array seismology techniques to the Apollo data and revealed detailed core structure, including the first direct confirmation of the presence of a solid inner core. Our study focused on the identification of core ]reflected phases in deep moonquake seismograms. The resulting model of the Moon fs innermost structure was found to be consistent with the commonly accepted model. However, the modeled layer radii may vary by tens of kilometers, as is expected when accounting for uncertainties such as moonquake location, timing errors, and potential seismic heterogeneities. In addition, the modeled velocities may vary with a 1-to-1 trade ]off with the modeled reflector depth. The GRAIL (Gravity Recovery and Interior Laboratory) mission, launched in Sept. 2011, placed two nearly identical spacecraft in lunar orbit. The two satellites make extremely high-resolution measurements of the lunar gravity field, which can be used to constrain the interior structure of the Moon using a "crust to core" approach. GRAIL fs constraints on crustal thickness, mantle structure, core radius and stratification, and core state (solid vs. molten) will complement seismic investigations in several ways. Here we present a progress report on our efforts to advance our knowledge of the Moon fs internal structure using joint gravity and seismic analyses. We will focus on methodology, including 1) refinements to the seismic core constraint accomplished through array processing of Apollo seismic data, made by applying a set of travel time corrections based on GRAIL structure estimates local to each Apollo seismic station; 2) modeling deep lunar structure through synthetic seismograms, to test whether the seismic core model can reproduce the core reflections observed in the Apollo seismograms; and 3) a joint seismic and gravity inversion in which we attempt to fit a family of seismic structure models with the gravity constraints from GRAIL, resulting in maps of seismic velocities and densities that vary from a nominal model both laterally and with depth.

Weber, Renee↗

Accessing Martian Fluvial and Lacustrine Sediments by Landing in Holden Crater, Margaritifer Sinus

Rover missions to the surface of Mars after MER 2003, are likely to be centered around focused geologic field mapping. One objective with high priority in selecting landing sites for these missions will be to characterize the nature, spatial distribution, internal structure, composition, and depositional history of exposed sedimentary layered deposits by visiting a number of distributed outcrops identified previously (and with a high degree of certainty) from orbit. These deposits may contain prebiotic material, even fossil organisms, but their primary value will be to enable an assessment of the planet's climate at the time they were emplaced. High resolution imaging from a mobile rover will enable the detailed study of these deposits over a wide area, their internal structure and mineralogy at distributed localities, and could resolve biologically-derived structures (such as stromatolite-like textures) if they are present. With the addition of a spectrometer, it should be possible to ascertain the presence of carbonates, sulfates, organics, water (liquid, frost, and bound water), as well as a variety of silicate minerals in the context of the collected imagery. Such a mission approach is directly relevant to future exploration of Mars, because it provides the geologic context comparable to what a field geologist visiting a site for the first time would acquire. Rover missions after MER will likely have much better targeting and hazard avoidance landing systems, enabling access to planimetrically-challenged sites of high scientific interest. These vehicles will also likely have greater mobility than MER, capable of driving greater distances in a shorter amount of time. Many scientists and mission planners have realized the need to design a rover whose mobility can be comparable to the dimensions of its 3-sigma landing error ellipse.

Parker, T. J.↗

Aircraft electric field measurements: Calibration and ambient field retrieval

An aircraft locally distorts the ambient thundercloud electric field. In order to determine the field in the absence of the aircraft, an aircraft calibration is required. In this work a matrix inversion method is introduced for calibrating an aircraft equipped with four or more electric field sensors and a high-voltage corona point that is capable of charging the aircraft. An analytic, closed form solution for the estimate of a (3 x 3) aircraft calibration matrix is derived, and an absolute calibration experiment is used to improve the relative magnitudes of the elements of this matrix. To demonstrate the calibration procedure, we analyze actual calibration date derived from a Lear jet 28/29 that was equipped with five shutter-type field mill sensors (each with sensitivities of better than 1 V/m) located on the top, bottom, port, starboard, and aft positions. As a test of the calibration method, we analyze computer-simulated calibration data (derived from known aircraft and ambient fields) and explicitly determine the errors involved in deriving the variety of calibration matrices. We extend our formalism to arrive at an analytic solution for the ambient field, and again carry all errors explicitly.

Koshak, William J.↗

Wind Prediction Accuracy for Air Traffic Management Decision Support Tools

The performance of Air Traffic Management and flight deck decision support tools depends in large part on the accuracy of the supporting 4D trajectory predictions. This is particularly relevant to conflict prediction and active advisories for the resolution of conflicts and the conformance with of traffic-flow management flow-rate constraints (e.g., arrival metering / required time of arrival). Flight test results have indicated that wind prediction errors may represent the largest source of trajectory prediction error. The tests also discovered relatively large errors (e.g., greater than 20 knots), existing in pockets of space and time critical to ATM DST performance (one or more sectors, greater than 20 minutes), are inadequately represented by the classic RMS aggregate prediction-accuracy studies of the past. To facilitate the identification and reduction of DST-critical wind-prediction errors, NASA has lead a collaborative research and development activity with MIT Lincoln Laboratories and the Forecast Systems Lab of the National Oceanographic and Atmospheric Administration (NOAA). This activity, begun in 1996, has focussed on the development of key metrics for ATM DST performance, assessment of wind-prediction skill for state of the art systems, and development/validation of system enhancements to improve skill. A 13 month study was conducted for the Denver Center airspace in 1997. Two complementary wind-prediction systems were analyzed and compared to the forecast performance of the then standard 60 km Rapid Update Cycle - version 1 (RUC-1). One system, developed by NOAA, was the prototype 40-km RUC-2 that became operational at NCEP in 1999. RUC-2 introduced a faster cycle (1 hr vs. 3 hr) and improved mesoscale physics. The second system, Augmented Winds (AW), is a prototype en route wind application developed by MITLL based on the Integrated Terminal Wind System (ITWS). AW is run at a local facility (Center) level, and updates RUC predictions based on an optimal interpolation of the latest ACARS reports since the RUC run. This paper presents an overview of the study's results including the identification and use of new large mor wind-prediction accuracy metrics that are key to ATM DST performance.

Cole, Rod↗

A Wind-Driven Snow Redistribution Module for Alpine3d V3.3.0: Adaptations Designed for Downscaling Ice Sheet Surface Mass Balance

Ice sheets gain mass via snow accumulation at the ice sheet surface, which is the primary component of surface mass balance. On the Antarctic ice sheet, winds redistribute snow resulting in surface mass balance that is variable in both space and time. Representing wind-driven snow redistribution processes in models is critical for local assessments of surface mass balance, repeat altimetry studies, and interpretation of ice core accumulation records. To this end, we have adapted Alpine3D, an existing distributed snow modeling framework, to downscale Antarctic surface mass balance to horizontal resolutions up to 1 km. In particular, we have introduced a new two-dimensional advection-based wind-driven snow redistribution module that is driven by an offline coupling between WindNinja, a wind downscaling model, and Alpine3D. We then show that large accumulation variability can be at least partially explained by terrain-induced wind speed variations which subsequently redistribute snow around rolling topography. By comparing Alpine3D to airborne-derived snow accumulation measurements within a testing domain over Pine Island Glacier in West Antarctica, we demonstrate that our Alpine3D downscaling approach improves surface mass balance estimates when compared to MERRA-2, a global atmospheric reanalysis which we use as atmospheric forcing. In particular, when compared to MERRA-2, Alpine3D reduces simulated surface mass balance root mean squared error by 23.4 mm w.e.yr−1 (13%) and increases variance explained by 24%. Despite these improvements, Alpine3D still underestimates observed accumulation variability, thus providing an opportunity for future model improvement.

Eric Keenan↗

Follow That Satellite: EO-1 Maneuvers into Closed Formation With Landsat-7

As the Landsat-7 spacecraft continued NASA's historic program of earth imaging, begun over three decades ago, NASA launched the Earth Observing-1 (EO-1) spacecraft carrying examples of the next generation of Landsat-7 instruments. The validation method for these instruments was to have EO-1 fly in a close formation behind Landsat-7 on the same World Reference System path. From that formation hundreds of near coincident images would be taken by each spacecraft and compared to evaluate improvements in the EO-1 instruments. This paper will address the mission analysis required to launch and maneuver EO-1 into the formation with Landsat-7 where instrument validation was to occur plus a summary of completing the formation acquisition. EO-1 is required to operate one minute +/- 6 seconds behind Landsat-7 during the period of co-fly imaging with a cross track separation of within + 3 kilometers. This separation time can also be stated as a one minute +/- 6 seconds time difference in the Mean Local Time (MLT) at the descending nodes. Achieving the required MLT is heavily dependent on the time of launch. The EO-1 launch window, which had to accommodate the dual payloads of EO-1 and SAC-C, was very limited ranging from 0 to 22 seconds over the 16 day Landsat-7 WRS repeat cycle during which EO-1 was launched. Each EO-1 launch opportunity that occurred on a different day of a Landsat-7 16 day repeat cycle required a separate and distinct maneuver profile. These profiles varied significantly in duration and amount of onboard propellant required to achieve them. EO-1 launched on a day judged to have "medium" resource requirements for achieving the formation with Landsat-7. To phase EO-1 one minute behind Landsat-7 in the along track direction, a series of altitude adjusts separated by specific drift intervals were executed. Additional maneuvers slightly changed the EO-1 inclination to maintain the MLT requirements. Orbit maneuvers were planned and executed within errors of less than 1.5 percent and propellant usage was near nominal consuming 3.4 kilograms out of a launch and early orbit budget of 11 kilograms. The pre-launch 3 -sigma propellant budget allowed for 1+ years of EO-1 mission life. The success of the EO-1 launch and early orbit operations provided sufficient propellant for nearly 4 years of on orbit operations. Special action taken during the EO-1 maneuver period involved some maneuver re-planning to reduce concerns about a potential close approach between EO-1 and Landsat-7. The result of this re-planning was a safer close approach and improvement in future formation acquisition planning.

DeFazio, Robert↗

MODIS 3 Km Aerosol Product: Applications over Land in an Urban/suburban Region

MODerate resolution Imaging Spectroradiometer (MODIS) instruments aboard the Terra and Aqua satellites have provided a rich dataset of aerosol information at a 10 km spatial scale. Although originally intended for climate applications, the air quality community quickly became interested in using the MODIS aerosol data. However, 10 km resolution is not sufficient to resolve local scale aerosol features. With this in mind, MODIS Collection 6 is including a global aerosol product with a 3 km resolution. Here, we evaluate the 3 km product over the Baltimore/Washington D.C., USA, corridor during the summer of 2011, by comparing with spatially dense data collected as part of the DISCOVER-AQ campaign these data were measured by the NASA Langley Research Center airborne High Spectral Resolution Lidar (HSRL) and a network of 44 sun photometers (SP) spaced approximately 10 km apart. The HSRL instrument shows that AOD can vary by up to 0.2 within a single 10 km MODIS pixel, meaning that higher resolution satellite retrievals may help to characterize aerosol spatial distributions in this region. Different techniques for validating a high-resolution aerosol product against SP measurements are considered. Although the 10 km product is more statistically reliable than the 3 km product, the 3 km product still performs acceptably, with more than two-thirds of MODIS/SP collocations falling within the expected error envelope with high correlation (R > 0.90). The 3 km product can better resolve aerosol gradients and retrieve closer to clouds and shorelines than the 10 km product, but tends to show more significant noise especially in urban areas. This urban degradation is quantified using ancillary land cover data. Overall, we show that the MODIS 3 km product adds new information to the existing set of satellite derived aerosol products and validates well over the region, but due to noise and problems in urban areas, should be treated with some degree of caution.

MODIS↗