Search NASA⌕ Search

SEARCH · Search NASA

Results for “temporal reasoning”

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 325 records · Page 18

Global Retrieval of Aerosol Properties over Desert and Semi-Desert Regions from SeaWiFS and MODIS

Mineral aerosols (dust) play an important role in both climate forcing and oceanic productivity throughout the entire year. Due to the relatively short lifetime (a few hours to about a week), the distributions of these airborne dust particles vary extensively in both space and time. Consequently, satellite observations are needed over both source and sink regions for continuous temporal and spatial sampling of dust properties. However, despite their importance, the high spatial resolution satellite measurements of dust near its source have been lacking. In this paper, we will demonstrate the capability of a new satellite algorithm to retrieve aerosol optical thickness and single scattering albedo over bright-reflecting surfaces such as urban areas and deserts. Such retrievals have been difficult to perform using previously available algorithms that use wavelengths from the mid-visible to the near IR because they have trouble separating the aerosol signal from the contribution due to the bright surface reflectance. The new algorithm, called Deep Blue, utilizes blue-wavelength measurements from instruments such as SeaWiFS and MODIS to infer the properties of aerosols, sinre the stirfare reflectance nver land in the blue part of the spectrum is much lower than for longer wavelength channels. We have validated the satellite retrieved aerosol optical thickness with data from AERONET sunphotometers over desert and semi-desert regions. The comparisons show reasonable agreements between these two. These new satellite products will allow scientists to determine quantitatively the aerosol properties near sources using high spatial resolution measurements from SeaWiFS and MODIS-like instruments.

Hsu, Nai-Yung↗

A Comparison Between Gravity Wave Momentum Fluxes in Observations and Climate Models

For the first time, a formal comparison is made between gravity wave momentum fluxes in models and those derived from observations. Although gravity waves occur over a wide range of spatial and temporal scales, the focus of this paper is on scales that are being parameterized in present climate models, sub-1000-km scales. Only observational methods that permit derivation of gravity wave momentum fluxes over large geographical areas are discussed, and these are from satellite temperature measurements, constant-density long-duration balloons, and high-vertical-resolution radiosonde data. The models discussed include two high-resolution models in which gravity waves are explicitly modeled, Kanto and the Community Atmosphere Model, version 5 (CAM5), and three climate models containing gravity wave parameterizations,MAECHAM5, Hadley Centre Global Environmental Model 3 (HadGEM3), and the Goddard Institute for Space Studies (GISS) model. Measurements generally show similar flux magnitudes as in models, except that the fluxes derived from satellite measurements fall off more rapidly with height. This is likely due to limitations on the observable range of wavelengths, although other factors may contribute. When one accounts for this more rapid fall off, the geographical distribution of the fluxes from observations and models compare reasonably well, except for certain features that depend on the specification of the nonorographic gravity wave source functions in the climate models. For instance, both the observed fluxes and those in the high-resolution models are very small at summer high latitudes, but this is not the case for some of the climate models. This comparison between gravity wave fluxes from climate models, high-resolution models, and fluxes derived from observations indicates that such efforts offer a promising path toward improving specifications of gravity wave sources in climate models.

mathematical models↗

Global Retrieval of Aerosol Properties from Sources to Sinks By MODIS

Mineral dust and smoke aerosols play an important role in both climate forcing and oceanic productivity throughout the entire year. Due to the relatively short lifetime (a few hours to about a week), the distributions of these airborne particles vary extensively in both space and time. Consequently, satellite observations are needed over both source and sink regions for continuous temporal and spatial sampling of dust and smoke properties. However, despite their importance, the high spatial resolution satellite measurements of these aerosols near their sources have been lacking, In this paper, we will demonstrate the capability of a new satellite algorithm to retrieve aerosol optical thickness and single scattering albedo over bright-reflecting surfaces such as urban areas and deserts. Such retrievals have been difficult to perform using previously available algorithms that use wavelengths from the mid-visible to the near IR because they have trouble separating the aerosol signal from the contribution due to the bright surface reflectance. The new algorithm, called Deep Blue, utilizes blue-wavelength measurements from instruments such as MODIS and SeaWiFS to infer the properties of aerosols, since the surface reflectance over land in the blue part of the spectrum is much lower than for longer wavelength channels. We have validated the satellite retrieved aerosol optical thickness with data from AERONET sunphotometers over land, including desert and semi-desert regions. The comparisons show reasonable agreements between these two. Our results show that the dust plumes lifted from the deserts near India/Pakistan border, and over Afghanistan, and the Arabian Peninsula are often observed by MODIS to be transported along the Indo-Gangetic Basin and mixed with the fine mode pollution particles generated by anthropogenic activities in this region, particularly during the pre-monsoon season (April-May). These new satellite products will allow scientists to determine quantitatively the aerosol properties near sources using high spatial resolution measurements from SeaWiFS and MODIS-like instruments.

Hsu, N. Christina↗

Summary of Quantitative Interpretation of Image Far Ultraviolet Auroral Data

Direct imaging of the magnetosphere by instruments on the IMAGE spacecraft is supplemented by simultaneous observations of the global aurora in three far ultraviolet (FUV) wavelength bands. The purpose of the multi-wavelength imaging is to study the global auroral particle and energy input from thc magnetosphere into the atmosphere. This paper describes provides the method for quantitative interpretation of FUV measurements. The Wide-Band Imaging Camera (WIC) provides broad band ultraviolet images of the aurora with maximum spatial and temporal resolution by imaging the nitrogen lines and bands between 140 and 180 nm wavelength. The Spectrographic Imager (SI), a dual wavelength monochromatic instrument, images both Doppler-shifted Lyman alpha emissions produced by precipitating protons, in the SI-12 channel and OI 135.6 nm emissions in the SI-13 channel. From the SI-12 Doppler shifted Lyman alpha images it is possible to obtain the precipitating proton flux provided assumptions are made regarding the mean energy of the protons. Knowledge of the proton (flux and energy) component allows the calculation of the contribution produced by protons in the WIC and SI-13 instruments. Comparison of the corrected WIC and SI-13 signals provides a measure of the electron mean energy, which can then be used to determine the electron energy fluxun-. To accomplish this reliable modeling emission modeling and instrument calibrations are required. In-flight calibration using early-type stars was used to validate the pre-flight laboratory calibrations and determine long-term trends in sensitivity. In general, very reasonable agreement is found between in-situ measurements and remote quantitative determinations.

Frey, H. U.↗

Overall Traveling-Wave-Tube Efficiency Improved By Optimized Multistage Depressed Collector Design

Depressed Collector Design The microwave traveling wave tube (TWT) is used widely for space communications and high-power airborne transmitting sources. One of the most important features in designing a TWT is overall efficiency. Yet, overall TWT efficiency is strongly dependent on the efficiency of the electron beam collector, particularly for high values of collector efficiency. For these reasons, the NASA Glenn Research Center developed an optimization algorithm based on simulated annealing to quickly design highly efficient multistage depressed collectors (MDC's). Simulated annealing is a strategy for solving highly nonlinear combinatorial optimization problems. Its major advantage over other methods is its ability to avoid becoming trapped in local minima. Simulated annealing is based on an analogy to statistical thermodynamics, specifically the physical process of annealing: heating a material to a temperature that permits many atomic rearrangements and then cooling it carefully and slowly, until it freezes into a strong, minimum-energy crystalline structure. This minimum energy crystal corresponds to the optimal solution of a mathematical optimization problem. The TWT used as a baseline for optimization was the 32-GHz, 10-W, helical TWT developed for the Cassini mission to Saturn. The method of collector analysis and design used was a 2-1/2-dimensional computational procedure that employs two types of codes, a large signal analysis code and an electron trajectory code. The large signal analysis code produces the spatial, energetic, and temporal distributions of the spent beam entering the MDC. An electron trajectory code uses the resultant data to perform the actual collector analysis. The MDC was optimized for maximum MDC efficiency and minimum final kinetic energy of all collected electrons (to reduce heat transfer). The preceding figure shows the geometric and electrical configuration of an optimized collector with an efficiency of 93.8 percent. The results show the improvement in collector efficiency from 89.7 to 93.8 percent, resulting in an increase of three overall efficiency points. In addition, the time to design a highly efficient MDC was reduced from a month to a few days. All work was done in-house at Glenn for the High Rate Data Delivery Program. Future plans include optimizing the MDC and TWT interaction circuit in tandem to further improve overall TWT efficiency.

Vaden, Karl R.↗

X-Ray Fluctuation Power Spectral Densities of Seyfert 1 Galaxies

By combining complementary monitoring observations spanning long, medium and short time scales, we have constructed power spectral densities (PSDs) of six Seyfert 1 galaxies. These PSDs span approx. greater than 4 orders of magnitude in temporal frequency, sampling variations on time scales ranging from tens of minutes to over a year. In at least four cases, the PSD shows a "break," a significant departure from a power law, typically on time scales of order a few days. This is similar to the behavior of Galactic X-ray binaries (XRBs), lower mass compact systems with breaks on time scales of seconds. NGC 3783 shows tentative evidence for a doubly-broken power law, a feature that until now has only been seen in the (much better-defined) PSDs of low-state XRBs. It is also interesting that (when one previously-observed object is added to make a small sample of seven), an apparently significant correlation is seen between the break time scale T and the putative black hole mass M(sub BH), while none is seen between break time scale and luminosity. The data are consistent with the linear relation T = M(sub BH) /10(exp 6.5) solar mass; extrapolation over 6-7 orders of magnitude is in reasonable agreement with XRBs. All of this strengthens the case for a physical similarity between Seyfert 1s and XRBs.

Markowitz, A.↗

Discrete Aurora on the Nightside of Mars: Occurrence Location and Probability

This paper represents the first attempt to predict the occurrence location and probability of discrete electron aurora on the nightside of Mars. We run a 3-D time-dependent magnetohydrodynamic model to characterize the spatial and temporal dynamics of magnetic field and plasma distributions over the course of one planetary rotation. We perform eight simulation cases under solar minimum quiet-solar-wind conditions (four equinox/solstice seasons, each with two interplanetary magnetic field polarities) and in an actual interplanetary coronal mass ejection (ICME) case to assess quiet and space weather situations, respectively. The occurrence of detectable discrete aurora is subject to the combination of the probabilities that (a) the ionosphere is magnetically connected with high altitudes through open field lines and (b) precipitating energy fluxes of >30 eV electrons exceed 0.1 erg/cm2/s. Our results show that during quiet solar activity, discrete aurora occurs likely on small-scale patches embedded inside strong crustal magnetic field regions (with a magnitude greater than 50 nT at 150 km), and the overall chance across the globe is ∼0.77%. The higher probability over strong crustal field regions is attributed to the stronger magnetic field convergence. Modeling shows the occurrence probability dramatically increases during the ICME event, particularly by more than an order of magnitude in weak crustal field regions. Our model results reasonably agree with NASA Mars Atmosphere and Volatile EvolutioN and Mars Express observations. Our study suggests that nightside discrete electron aurora is not caused by the direct entry of magnetosheath plasma in a cusp-like process but due to the recycling of nightside magnetospheric electrons.

Xiaohua Fang↗

A Synthesis of Light Absorption Properties of the Arctic Ocean: Application to Semi-analytical Estimates of Dissolved Organic Carbon Concentrations from Space

The light absorption coefficients of particulate and dissolved materials are the main factors determining the light propagation of the visible part of the spectrum and are, thus, important for developing ocean color algorithms. While these absorption properties have recently been documented by a few studies for the Arctic Ocean [e.g., Matsuoka et al., 2007, 2011; Ben Mustapha et al., 2012], the datasets used in the literature were sparse and individually insufficient to draw a general view of the basin-wide spatial and temporal variations in absorption. To achieve such a task, we built a large absorption database at the pan-Arctic scale by pooling the majority of published datasets and merging new datasets. Our results showed that the total non-water absorption coefficients measured in the Eastern Arctic Ocean (EAO; Siberian side) are significantly higher 74 than in the Western Arctic Ocean (WAO; North American side). This higher absorption is explained 75 by higher concentration of colored dissolved organic matter (CDOM) in watersheds on the Siberian 76 side, which contains a large amount of dissolved organic carbon (DOC) compared to waters off 77 North America. In contrast, the relationship between the phytoplankton absorption (a()) and chlorophyll a (chl a) concentration in the EAO was not significantly different from that in the WAO. Because our semi-analytical CDOM absorption algorithm is based on chl a-specific a() values [Matsuoka et al., 2013], this result indirectly suggests that CDOM absorption can be appropriately erived not only for the WAO but also for the EAO using ocean color data. Derived CDOM absorption values were reasonable compared to in situ measurements. By combining this algorithm with empirical DOC versus CDOM relationships, a semi-analytical algorithm for estimating DOC concentrations for coastal waters at the Pan-Arctic scale is presented and applied to satellite ocean color data.

Pan-Arctic↗

Advanced Energy Scale Correction Techniques for the X-ray Transition Edge Sensors of the Athena mission

The X-ray Integral Field Unit (X-IFU) onboard the future European X-ray telescope Athena will be the first space instrument carrying an array of more than a thousand transition edge sensors. One of the key challenges of the X-IFU is the measurement of narrow X-ray atomic lines to determine velocity shifts at an unprecedented level of accuracy. For this reason, the energy scale of the instrument needs to be known with extreme accuracy, of 0.4 eV (1σ) up to 7 keV. The energy scale will be measured on the ground through a dedicated calibration campaign using fiducial X-ray sources. Though calibrated, the energy scale is extremely sensitive to the environmental conditions around the TES array, and drifts in the readout chain electronics. Uncorrected, the energy scale can naturally drift up to hundreds of eVs. Changes of the TES gain will be monitored via onboard X-ray calibration sources, and the energy scale will be corrected either per pixel, or within a small groups of pixels. Although simulations show that a 0.4 eV level can be achieved, the very high accuracy required by the X-IFU calls for experimental validation. A dedicated measurement campaign has been performed by NASA Goddard Space Flight Center to characterize the energy scale of a prototype kilo-pixel array of X-IFU-representative TESs. The analysis of the data demonstrated the ability to correct for various drifts using two fiducial lines to track the temporal gain variation. In this paper, we propose to extend this study on the same data set by investigating multi-parameter correction techniques based on both the pulse-height of the fiducial line and the prepulse baseline level, using the knowledge of the TES energy scale at reference temperature/magnetic field set points acquired on the ground. Investigations on the co-adding of pixels to perform a joint correction over pools of pixels is also explored.

79 ASTRONOMY AND ASTROPHYSICS↗

Interference Lattice-based Loop Nest Tilings for Stencil Computations

A common method for improving performance of stencil operations on structured multi-dimensional discretization grids is loop tiling. Tile shapes and sizes are usually determined heuristically, based on the size of the primary data cache. We provide a lower bound on the numbers of cache misses that must be incurred by any tiling, and a close achievable bound using a particular tiling based on the grid interference lattice. The latter tiling is used to derive highly efficient loop orderings. The total number of cache misses of a code is the sum of (necessary) cold misses and misses caused by elements being dropped from the cache between successive loads (replacement misses). Maximizing temporal locality is equivalent to minimizing replacement misses. Temporal locality of loop nests implementing stencil operations is optimized by tilings that avoid data conflicts. We divide the loop nest iteration space into conflict-free tiles, derived from the cache miss equation. The tiling involves the definition of the grid interference lattice an equivalence class of grid points whose images in main memory map to the same location in the cache-and the construction of a special basis for the lattice. Conflicts only occur on the boundaries of the tiles, unless the tiles are too thin. We show that the surface area of the tiles is bounded for grids of any dimensionality, and for caches of any associativity, provided the eccentricity of the fundamental parallelepiped (the tile spanned by the basis) of the lattice is bounded. Eccentricity is determined by two factors, aspect ratio and skewness. The aspect ratio of the parallelepiped can be bounded by appropriate array padding. The skewness can be bounded by the choice of a proper basis. Combining these two strategies ensures that pathologically thin tiles are avoided. They do not, however, minimize replacement misses per se. The reason is that tile visitation order influences the number of data conflicts on the tile boundaries. If two adjacent tiles are visited successively, there will be no replacement misses on the shared boundary. The iteration space may be covered with pencils larger than the size of the cache while avoiding data conflicts if the pencils are traversed by a scanning-face method. Replacement misses are incurred only on the boundaries of the pencils, and the number of misses is minimized by maximizing the volume of the scanning face, not the volume of the tile. We present an algorithm for constructing the most efficient scanning face for a given grid and stencil operator. In two dimensions it is based on a continued fraction algorithm. In three dimensions it follows Voronoi's successive minima algorithm. We show experimental results of using the scanning face, and compare with canonical loop orderings.

VanderWijngaart, Rob F.↗

Evaluation of Global Surface Water Temperature Data Sets for use in Passive Remote Sensing of Soil Moisture

Inland open water bodies often pose a systematic error source in the passive remote sensing retrievals of soil moisture. Water temperature is a necessary variable used to compute water emissions that is required to be subtracted from satellite observation to yield actual emissions from the land portion, which in turn generates accurate soil moisture retrievals. Therefore, overestimation of soil moisture can often be corrected using concurrent water temperature data in the overall mitigation procedure. In recent years, several data sets of lake water temperature have become available, but their specifications and accuracy have rarely been investigated in the context of passive soil moisture remote sensing on a global scale. For this reason, three lake temperature products were evaluated against in-situ measurements from 2007 to 2011. The data sets include the lake surface water temperature (LSWT) from Global Observatory of Lake Responses to Environmental Change (GloboLakes), the Copernicus Global Land Operations Cryosphere and Water (C-GLOPS), as well as the lake mix-layer temperature (LMLT) from the European Centers for Medium-Range Weather Forecast (ECMWF) ERA5 Land Reanalysis. GloboLakes, C-GLOPS, and ERA5 Land have overall comparable performance with Pearson correlations (R) of 0.87, 0.92 and 0.88 in comparison with in-situ measurements. LSWT products exhibit negative median biases of −0.27 K (GloboLakes) and −0.31 K (C-GLOPS), whereas the median bias of LMLT is 1.56 K. When mapped from their respective native resolutions to a common 9 km Equal-Area Scalable Earth (EASE) Grid 2.0 projection, similar relative performance was observed. LMLT and LSWT data are closer in performance over the 9 km grid cells that exhibit a small range of lake cover fractions (0.05–0.5). Despite comparable relative performance, ERA5 Land shows great advantages in spatial coverage and temporal resolution. In summary, an integrated evaluation on data accuracy, long-term availability, global coverage, temporal resolution, and regular forward processing with modest data latency led us to conclude that LMLT from the ERA5 Land Reanalysis product represents the most optimal path for use in the development of a long-term soil moisture product.

Runze Zhang↗

X-HAB 2020: AR Field Treks Summary and Conclusions

As part of the FY20 X-Hab Challenge, BLiSS sought to create an Augmented Reality (AR) toolkit to help with analog field trek operations under the supervision of the Solar System Exploration Research Virtual Institute (SSERVI). These treks are operational and technical demonstrations at space-like destinations on Earth to test current extra-vehicular activity (EVA) techniques. While BLiSS as an organization has experience studying operational tasks such as this, it has never developed AR software at this scale. For that reason, another team at the University was brought on to work in parallel. The Collaborative Lab for Advancing Work in Space (CLAWS) is a veteran group of the NASA Spacesuit User Interface Technologies for Students (SUITS) challenge in which Hololens displays for astronauts are created within a year. The operational and technological pairing was ideally suited for tackling this problem. The team divided its responsibilities so that BLiSS would handle the research required to shape the project. As this deliverable had an end user, it was decided that interviewing these field geologists and operations specialists would provide the best insight. These interviews paired with literature review would reveal niche applications for AR that remained within feasible bounds. These science-driven EVAs in unknown terrain require more flexible tools than the current generation of EVA assistants. Rather than focus on sequential instructions, there instead needs to be a broad toolkit that's only called upon in specific instances. This AR Toolkit for Lunar Astronauts and Scientists (ATLAS) became the development goal of the project: create a non-intrusive assembly of tools that could be accessed in AR on the field. The current ATLAS design makes use of a geospatially and temporally annotated eld note system called GeoNotes. This allows for data to be collected and coordinated in a way that's synchronized across time, space, and different users. A Mission Control Center (MCC) and Mobile Support Equipment (MSE) were all needed to transport the AR headset into the field with the user. A network infrastructure was designed and set up within the University to enable this functionality. The software is based on a Protocol-Module structure that allows for modular development of each capability. A Protocol Manager coordinates different protocols that make use of modules. Each module tackles a different individual task while the protocol puts each one to use. The protocol manager coordinates when these are called to be used. This software is hosted on a head-mounted display (HMD) with the MCC acting as support from afar. While the software would be unit-tested at each level and each hardware component verified, a final demonstration would serve to prove the system's capabilities: an analog field trek. The team would prepare to support a user in a remote location from the MCC back at the University. A local area near campus would be tested before going out to do sample field geology further away. This unfortunately became impossible with the arrival of COVID-19. Access to all of the facilities to complete the project as planned were shut down. Our team was scattered across the globe and forced to complete the rest virtually. Adjustments were made to produce a small virtual concept in Adobe XD in the meantime. Even digital surveys were created based on the NASA task-load index (TLX) originally intended for testing actual users. The goal shifted towards completing software and getting feedback on the user interfaces (UI) and user experiences (UX). This team has reformed in response to COVID and its focus has shifted to what can be done remotely. There is still an intention to finish the original deliverable described in this report. The work has been expanded beyond the original X-Hab challenge and has instead become its own research e ort to be continued afterwards. This report collects the processes and knowledge gained from a year of studying and working at this problem with two teams. It should preserve it for the time until the world returns to normal and work can resume. CLAWS will be taking over full responsibility from that point forward, eventually surpassing the original needs of the project. While this document captures the work done towards an eventual end, the CLAWS team has written their own proposal alongside it. It outlines a new future for ATLAS beyond X-Hab, BLiSS, and hopefully beyond COVID-19. This project began as a vague goal hoping to place a new technology into the unique setting of exploration science. The project has since comfortably taken root and will hopefully bloom over the next year.

Alex Sena↗

The use of aircraft and satellite remote sensing of phytoplankton chlorophyll concentrations in case 2 estuarine waters of the Chesapeake Bay

Two projects using remote sensing of phytoplankton chlorophyll concentrations in the Chesapeake Bay estuary were proposed. The first project used aircraft remote sensing with a compact radiometer system developed at NASA's Goddard Space Flight Center (GSFC), the Ocean Data Acquisition System (ODAS). ODAS includes three radiometers at 460, 490, and 520 nm, an infrared temperature sensor (PRT-5), Loran-C for navigation, and a data acquisition system using a PC and mass storage device. This instrument package can be flown in light aircraft at relatively low expense, permitting regular and frequent flights. Sixteen flights with ODAS were completed using the Virginia Institute of Marine Science's De Havilland 'Beaver'. The goal was to increase spatial and temporal resolution in assaying phytoplankton pigment concentrations in the Chesapeake. At present, analysis is underway of flight data collected between March and July 1989. The second project focused on satellite data gathered with the Nimbus-7 Coastal Zone Color Scanner (CZSC) between late 1978 and mid 1986. The problem in using CZSC data for the Chesapeake Bay is that the optical characteristics of this (and many) coastal and estuarine waters are distinct from those of the open ocean for which algorithms for computing pigment concentrations were developed. The successful use of CZCS data for the estuary requires development of site-specific algorithms and analytical approaches. Of principal importance in developing site-specific procedures is the availability of in-situ data on pigment concentrations. A significant data set was acquired from EPA's Chesapeake Bay Program in Annapolis, Maryland, and clear satellite scenes are being analyzed for which same-day sea truth measurements of pigment were obtained. Both the University of Miami and GSFC Seapak systems are being used in this effort. The main finding to date is an expected one, i.e., the algorithms developed for oceanic waters are inadequate to compute pigment concentrations for the Case 2 waters of the Chesapeake Bay. One reason is the overestimation of aerosol radiances by assuming that water-leaving radiance in Band 4 of CZCS (670 nm) is zero, an assumption that is invalid for the Bay. This prompted any attempts to iterative procedures for estimating the proportion of the Band 4 radiance that is actually attributable to aerosol by estimating the water-leaving component using optical data. A cruise on the Chesapeake the week of 7 August 1989 was conducted to collect additional optical data necessary to this task.

Harding, Lawrence W., Jr.↗

Aerial Vehicle Routing and Scheduling for UAS Traffic Management: A Monte Carlo Tree Search Approach

Numerous unmanned aircraft systems operating at low altitudes to deliver goods and services may one day become ubiquitous in our cities. In the Unmanned Aircraft Systems (UAS) Traffic Management (UTM) framework, such a concept is envisioned, where aerial vehicles operate beyond visual line of sight (BVLOS) within specifically reserved and time stamped “corridors” in the airspace. For example, these corridors or operational intent volumes can connect an aerial vehicle’s origin site to its destination site for package delivery operations. There may also be more than one corridor available for an aerial vehicle to choose from and often different corridors may intersect with one another. Thus, it is imperative to ensure flight trajectories belonging to different aerial vehicles are not in conflict. Per the UTM CONOPs, we assume that a vehicle almost always stays inside its corridor or operational volume. This work provides a framework for strategic deconfliction of UTM or package delivery drones, where we schedule the departure time of all vehicles subject to various temporal constraints (including the corridor deconfliction at the intersections). We present the “multi-route weighted package delivery problem” which serves as an exemplifying model for strategic deconfliction in UTM. In the multi-route weighted package delivery problem, a graph network is given which consists of a set of depots (source) and drop-off (destination) nodes, with multiple routes (defined as a sequence of waypoints) connecting the depots to drop-off nodes. In addition, routes are weighted by the associated ground risk and total travel distance for package delivery. The goal is for a known set of aerial vehicles to depart from the depots, choose a route and take off time, while avoiding conflicts with other aerial vehicles, and minimizing both risk and distance traveled. We provide a mixed integer linear programming (MILP) formulation of the problem, as well as a heuristic solution based on Monte Carlo Tree Search (MCTS) – a method used in game theory and artificial intelligence – to overcome limitations inherent to optimal solvers. Computational results show the advantages of using MCTS over the MILP formulation; the former can provide a sub-optimal solution quickly, and may sometimes even reach an optimal solution, whereas the latter may not even produce a solution in reasonable time. Furthermore, results from both the MILP formulation and MCTS methods were validated using a preliminary agent-based simulator implementing the UTM concept of operations. Thus, the MCTS method can be seen as a scalable solution to the complex multi-route weighted package delivery problem and may possibly be extended to similar complex optimization problems.

Kenny Chour↗

NPP Clouds and the Earth's Radiant Energy System (CERES) Predicted Sensor Performance Calibration and Preliminary Data Product Performance

Continuation of the Earth Radiation Budget (ERB) Climate Data Record (CDR) has been identified as critical in the 2007 NRC Decadal Survey, the Global Climate Observing System WCRP report, and in an assessment titled Impacts of NPOESS Nunn-McCurdy Certification on Joint NASA-NOAA Climate Goals. In response, NASA, NOAA and NPOESS agreed in early 2008 to fly the final existing CERES Flight Model (FM-5) on the NPP spacecraft for launch in 2010. Future opportunities for ERB CDR continuity consist of procuring an additional CERES Sensor with modest performance upgrades for flight on the NPOESS C1 spacecraft in 2013, followed by a new CERES follow-on sensor for flight in 2018 on the NPOESS C3 spacecraft. While science goals remain unchanged for the long-term ERB Climate Data Record, it is now understood that the task of achieving these goals is more difficult for two reasons. The first is an increased understanding of the dynamics of the Earth/atmosphere system which demonstrates that rigorous separation of natural variability from anthropogenic change on decadal time scales requires higher accuracy and stability than originally envisioned. Secondly, future implementation scenarios involve less redundancy in flight hardware (1 vs. 2 orbits and operational sensors) resulting in higher risk of loss of continuity and reduced number of independent observations to characterize performance of individual sensors. Although EOS CERES CDR's realize a factor of 2 to 4 improvement in accuracy and stability over previous ERBE CDR's, future sensors will require an additional factor of 2 improvement to answer rigorously the science questions moving forward. Modest investments, defined through the CERES Science Team s 30-year operational history of the EOS CERES sensors, in onboard calibration hardware and pre-flight calibration and test program will ensure meeting these goals while reducing costs in re-processing scientific datasets. The CERES FM-5 pre-flight radiometric characterization program benefited from the 30-year operational experience of the CERES EOS sensors, as well as a stronger emphasis of radiometric characterization in the Statement of Work with the sensor provider. Improvements to the pre-flight program included increased spectral, spatial, and temporal sampling under vacuum conditions as well as additional tests to characterize the primary and transfer standards in the calibration facility. Future work will include collaboration with NIST to further enhance the understanding of the radiometric performance of this equipment prior to flight. The current effort summarizes these improvements to the CERES FM-5 pre-flight sensor characterization program, as well as modifications to inflight calibration procedures and operational tasking. In addition, an estimate of the impacts to the system level accuracy and traceability is presented.

Priestly, Kory↗

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗

Deep Ocean Warming Assessed from Altimeters, GRACE, 3 In-situ Measurements, and a Non-Boussinesq OGCM

Observational surveys have shown significant oceanic bottom water warming, but they are too spatially and temporally sporadic to quantify the deep ocean contribution to the present-day sea level rise (SLR). In this study, altimetry sea surface height (SSH), Gravity Recovery and Climate Experiment (GRACE) ocean mass, and in situ upper ocean (0-700 m) steric height have been assessed for their seasonal variability and trend maps. It is shown that neither the global mean nor the regional trends of altimetry SLR can be explained by the upper ocean steric height plus the GRACE ocean mass. A non-Boussinesq ocean general circulation model (OGCM), allowing the sea level to rise as a direct response to the heat added into the ocean, is then used to diagnose the deep ocean steric height. Constrained by sea surface temperature data and the top of atmosphere (TOA) radiation measurements, the model reproduces the observed upper ocean heat content well. Combining the modeled deep ocean steric height with observational upper ocean data gives the full depth steric height. Adding a GRACE-estimated mass trend, the data-model combination explains not only the altimetry global mean SLR but also its regional trends fairly well. The deep ocean warming is mostly prevalent in the Atlantic and Indian oceans, and along the Antarctic Circumpolar Current, suggesting a strong relation to the oceanic circulation and dynamics. Its comparison with available bottom water measurements shows reasonably good agreement, indicating that deep ocean warming below 700 m might have contributed 1.1 mm/yr to the global mean SLR or one-third of the altimeter-observed rate of 3.11 +/- 0.6 mm/yr over 1993-2008.

sea level↗

VIC-Global Parameter Dataset Sensitivity with the Variable Infiltration Capacity Model: Evaluating the importance of dynamic land surface parameters when using the VIC-Global parameter dataset

Accurate prediction of runoff is essential to water resources management, flood risk assessment, and ecosystem protection. However, many hydrological models still have relatively substantial limitations when representing the influence of land use and land cover (LULC) on runoff generation and routing. Changes in LULC, such as deforestation, urban expansion, agricultural intensification, and wetland loss, have been shown to alter the water balance at the land surface through fundamental hydrologic processes (e.g., interception, infiltration, evapotranspiration, and soil storage). However, it remains an open question what the exact magnitude and timing of these impacts are for the spatial and temporal scales commonly used in engineering applications. In this analysis we focus on one aspect of recent LULC change for assessing human impacts, which is urbanization. Specifically we seek to determine the impacts of urbanization on the magnitude and timing of surface runoff and baseflow in HUC-12 basins in Clark County, Nevada which has experienced rapid urbanization. We use the Variable Infiltration Capacity (VIC) hydrology model with a widely used off-the-shelf dataset of land surface parameters, VIC-Global, both of which have been commonly used in the past for water and energy balance modeling for large scale hydrologic studies. We examine two scenarios where the first scenario removes all urbanized land cover and parameterizes those areas of the basins as barren or open shrubland. The second scenario tests the opposite case where all areas of the basins are classified as urban regardless of their present classification. The results from the VIC model show there is a low sensitivity for daily surface runoff between scenarios. The daily baseflow values indicate similar low sensitivity to the classification change during specific periods, but then have substantial differences during other period when large precipitation events are occurring. This is likely due to the assumed parameter values for the urban land cover classification made by the VIC-Global dataset. Using a static land cover parameterization is reasonable for large domain hydrology models that are being used for near-term planning horizons (<30 years). However, longer planning horizons where feedbacks between the atmosphere and land surface are important, especially in transient climate situations, considerations for how to update land surface parameters should be incorporated.

42 ENGINEERING↗