Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Marginal Distribution”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Structure-Informed Graph Learning of Networked Dependencies for Online Prediction of Power System Transient Dynamics

Online transient analysis plays an increasingly important role in dynamic power grids as the renewable generation continues growing. Traditional numerical methods for transient analysis not only are computationally intensive but also require precise contingency information as input, and therefore, are not suitable for online applications. Existing online transient assessment studies focus on the determination of post-contingency system stability or stability margin. Here, this paper develops a novel graph-learning framework, Deep-learning Neural Representation or DNR, for online prediction, of the time-series trajectories of the system states using initial system responses that can be measured by phasor measurement units (PMUs). The proposed DNR framework consists of two sequential modules: a Network Constructor that captures network dependencies among generators, and a Dynamics Predictor that predicts the system trajectories. The key to improved prediction performance is the introduction of the spatio-temporal message-passing operations into graph neural networks with structural knowledge. Its effectiveness and scalability are validated through comparative studies, demonstrating the prediction performance under different contingency scenarios for systems of different sizes. This framework provides a solution to online predicting post-fault system dynamics based on real-time PMU measurements. Additionally, it can also be applied to facilitate the offline transient simulation without simulating the entire trajectories.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Real-Time Regional PV Spinning Reserve Estimator with AGC Look-Ahead Windows

Curtailed PV generation is a zero-marginal cost spinning reserve that can be used for a number of active power control services. However, unlike the traditional spinning reserve providers, i.e., fossil-fueled generators, who have well-defined operating characteristics, e.g., available headroom or potential high limit (PHL), PV plants have by nature variable and uncertain operating characteristics. To ensure the effective coordination between PV plants and the system operator during an active power control event, accurate forecasts of the PV PHL are essential. A novel reference-control grouping based scaling method has been proposed by NREL to estimate the PV PHL in real-time. This work further enhances the methodology by: 1) improving the model accuracy through machine learning; 2) considering look-ahead windows introduced by the computation and communication latencies; 3) applying the method to regional spinning reserve estimation. A significant performance improvement, over 99% of estimation error reduction, has been observed based on real-world data collected by CAISO and PV plant operators.

potential high limit↗
DOE OSTI2023

Real-Time Regional PV Spinning Reserve Estimator with AGC Look-Ahead Windows: Preprint

Curtailed PV generation is a zero-marginal cost spinning reserve that can be used for a number of active power control services. However, unlike the traditional spinning reserve providers, i.e., fossil-fueled generators, who have well-defined operating characteristics, e.g., available headroom or potential high limit (PHL), PV plants have by nature variable and uncertain operating characteristics. To ensure the effective coordination between PV plants and the system operator during an active power control event, accurate forecasts of the PV PHL are essential. A novel reference-control grouping based scaling method has been proposed by NREL to estimate the PV PHL in real-time. This work further enhances the methodology by: 1) improving the model accuracy through machine learning; 2) considering look-ahead windows introduced by the computation and communication latencies; 3) applying the method to regional spinning reserve estimation. A significant performance improvement, over 1.6% of estimation error reduction, has been observed based on real-world data collected by CAISO and PV plant operators.

potential high limit↗
DOE OSTI2024

Long-run Marginal Emission Rates for Electricity - Workbooks for 2023 Cambium Data

These workbooks contain modeled estimates of long-run marginal emission rates (LRMER) for the contiguous United States' electric sector. A LRMER is an estimate of the rate of emissions that would be either induced or avoided by a change in electric demand, taking into account how the change could influence both the operation as well as the structure of the grid (i.e., the building and retiring of capital assets, such as generators and transmission lines). These workbooks provide data for 18 GEA regions covering the contiguous United States. Mappings of these regions to ZIP codes and counties is given in this workbook in the corresponding tabs. For more data underlying these emissions factors, see the Cambium 2023 project at https://scenarioviewer.nrel.gov/. For more details on input assumptions and methodology see the associated report (Cambium 2023 Scenario Descriptions and Documentation, https://www.nrel.gov/docs/fy24osti/88507.pdf). Users are advised to review section 4 of the report, which discusses limitations and caveats of the data. This data is planned to be updated annually. Information on the latest versions can be found at https://www.nrel.gov/analysis/cambium.html.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2014

Service Life Extension of the Propulsion System of Long-Term Manned Orbital Stations

One of the critical non-replaceable systems of a long-term manned orbital station is the propulsion system. Since the propulsion system operates beginning with the launch of station elements into orbit, its service life determines the service life of the station overall. Weighing almost a million pounds, the International Space Station (ISS) is about four times as large as the Russian space station Mir and about five times as large as the U.S. Skylab. Constructed over a span of more than a decade with the help of over 100 space flights, elements and modules of the ISS provide more research space than any spacecraft ever built. Originally envisaged for a service life of fifteen years, this Earth orbiting laboratory has been in orbit since 1998. Some elements that have been launched later in the assembly sequence were not yet built when the first elements were placed in orbit. Hence, some of the early modules that were launched at the inception of the program were already nearing the end of their design life when the ISS was finally ready and operational. To maximize the return on global investments on ISS, it is essential for the valuable research on ISS to continue as long as the station can be sustained safely in orbit. This paper describes the work performed to extend the service life of the ISS propulsion system. A system comprises of many components with varying failure rates. Reliability of a system is the probability that it will perform its intended function under encountered operating conditions, for a specified period of time. As we are interested in finding out how reliable a system would be in the future, reliability expressed as a function of time provides valuable insight. In a hypothetical bathtub shaped failure rate curve, the failure rate, defined as the number of failures per unit time that a currently healthy component will suffer in a given future time interval, decreases during infant-mortality period, stays nearly constant during the service life and increases at the end when the design service life ends and wear-out phase begins. However, the component failure rates do not remain constant over the entire cycle life. The failure rate depends on various factors such as design complexity, current age of the component, operating conditions, severity of environmental stress factors, etc. Development, qualification and acceptance test processes provide rigorous screening of components to weed out imperfections that might otherwise cause infant mortality failures. If sufficient samples are tested to failure, the failure time versus failure quantity can be analyzed statistically to develop a failure probability distribution function (PDF), a statistical model of the probability of failure versus time. Driven by cost and schedule constraints however, spacecraft components are generally not tested in large numbers. Uncertainties in failure rate and remaining life estimates increase when fewer units are tested. To account for this, spacecraft operators prefer to limit useful operations to a period shorter than the maximum demonstrated service life of the weakest component. Running each component to its failure to determine the maximum possible service life of a system can become overly expensive and impractical. Spacecraft operators therefore, specify the required service life and an acceptable factor of safety (FOS). The designers use these requirements to limit the life test duration. Midway through the design life, when benefits justify additional investments, supplementary life test may be performed to demonstrate the capability to safely extend the service life of the system. An innovative approach is required to evaluate the entire system, without having to go through an elaborate test program of propulsion system elements. Evaluating every component through a brute force test program would be a cost prohibitive and time consuming endeavor. ISS propulsion system components were designed and built decades ago. There are no representative ground test articles for some of the components. A 'test everything' approach would require manufacturing new test articles. The paper outlines some of the techniques used for selective testing, by way of cherry picking candidate components based on failure mode effects analysis, system level impacts, hazard analysis, etc. The type of testing required for extending the service life depends on the design and criticality of the component, failure modes and failure mechanisms, life cycle margin provided by the original certification, operational and environmental stresses encountered, etc. When specific failure mechanism being considered and the underlying relationship of that mode to the stresses provided in the test can be correlated by supporting analysis, time and effort required for conducting life extension testing can be significantly reduced. Exposure to corrosive propellants over long periods of time, for instance, lead to specific failure mechanisms in several components used in the propulsion system. Using Arrhenius model, which is tied to chemically dependent failure mechanisms such as corrosion or chemical reactions, it is possible to subject carefully selected test articles to accelerated life test. Arrhenius model reflects the proportional relationship between time to failure of a component and the exponential of the inverse of absolute temperature acting on the component. The acceleration factor is used to perform tests at higher stresses that allow direct correlation between the times to failure at a high test temperature to the temperatures to be expected in actual use. As long as the temperatures are such that new failure mechanisms are not introduced, this becomes a very useful method for testing to failure a relatively small sample of items for a much shorter amount of time. In this article, based on the example of the propulsion system of the first ISS module Zarya, theoretical approaches and practical activities of extending the service life of the propulsion system are reviewed with the goal of determining the maximum duration of its safe operation.

Kamath, Ulhas↗
NASA NTRS2000

International Space Station Evolution Data Book: Baseline Design - Volume 1

The International Space Station (ISS) will provide an Earth-orbiting facility that will accommodate engineering experiments as well as research in a microgravity environment for life and natural sciences. The ISS will distribute resource utilities and support permanent human habitation for conducting this research and experimentation in a safe and habitable environment. The objectives of the ISS program are to develop a world-class, international orbiting laboratory for conducting high-value scientific research for the benefit of humans on Earth; to provide access to the microgravity environment; to develop the ability to live and work in space for extended periods; and to provide a research test bed for developing advanced technology for human and robotic exploration of space. The current design and development of the ISS has been achieved through the outstanding efforts of many talented engineers, designers, technicians, and support personnel who have dedicated their time and hard work to producing a state-of-the-art Space Station. Despite these efforts, the current design of the ISS has limitations that have resulted from cost and technology issues. Regardless, the ISS must evolve during its operational lifetime to respond to changing user needs and long-term national and international goals. As technologies develop and user needs change, the ISS will be modified to meet these demands. The design and development of these modifications should begin now to prevent a significant lapse in time between the baseline design and the realization of future opportunities. For this effort to begin, an understanding of the baseline systems and current available opportunities for utilization needs to be achieved. Volume I of this document provides the consolidated overview of the ISS baseline systems. It also provides information on the current facilities available for pressurized and unpressurized payloads. Information on current plans for crew availability and utilization; resource timelines and margin summaries including power, thermal, and storage volumes; and an overview of the ISS cargo traffic and the vehicle traffic model is also included.

Jorgensen, Catherine A.↗
DOE OSTI2022

A computationally efficient algorithm for computing convex hull prices

Electricity markets worldwide allow participants to bid non-convex production offers. While non-convex offers can more accurately reflect a resource's capabilities, they create challenges for market clearing processes. For example, system operators may be required to execute side payments to participants whose costs are not covered through energy sales as determined via traditional locational marginal pricing schemes. Convex hull pricing minimizes this and other types of side payments while providing uniform (i.e., locationally and temporally consistent) prices. Computing convex hull prices involves solving either a large-scale linear program or the Lagrangian dual of the corresponding non-convex scheduling problem. Further, the former approach requires explicit descriptions of market participants' convex hulls. While linear programs for computing convex hull prices are large, their structure is naturally decomposable by generators. Here, in this work, we propose and empirically analyze a Benders decomposition approach to computing convex hull prices that leverages recent advances in convex hull formulations for thermal generating units. We demonstrate across a large set of test instances that our decomposition approach only requires modest computational effort, obtaining solutions at least an order of magnitude faster than the equivalent large-scale linear programming approach. Overall, we provide a computationally feasible method for computing convex hull prices for industrial scale market clearing problems, enabling the possibility of practical adoption of this advanced pricing mechanism.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

2022 Summer Resource Adequacy in the ERCOT Region

This report is a 2022 update to the National Energy Technology Laboratory (NETL) report that examined the potential system performance of the Electric Reliability Council of Texas (ERCOT) system during the summer of 2021. Anticipated summer reserve margins in ERCOT have risen year-over-year from the 15.7 percent reported in ERCOT’s Final Seasonal Assessment of Resource Adequacy (SARA) for Summer 2021 to 22.8 percent in the Final SARA for Summer 2022, driven mostly by increases in wind and solar generation. Using data from ERCOT’s 2022 Summer SARA, Capacity, Demand and Reserves (CDR) reports, and historical performance data obtained from ERCOT and Hitachi Energy Velocity Suite, six economic dispatch scenarios using PROMOD were developed to assess the risk to ERCOT of a load shedding event for Summer 2022. The results suggests that ERCOT may be facing a very tight resource adequacy situation, with the potential for a serious shortfall during the summer peak, usually the end of July through mid-August, if demand and generator outage rates significantly exceed peak levels. Under these conditions, ERCOT will need to call on its operating tools, including Energy Emergency Alerts (EEA) and calls to conserve electricity, to maintain reliability and continue serving load. Those measures alone may not be sufficient to avert a shedding of load.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Impacts of hybridization and forecast errors on the probabilistic capacity credit of batteries

Battery storage is increasingly identified as being among the least-cost mix of technologies in the evolving U.S. electricity mix. This study explores the marginal capacity credit of batteries using a probabilistic, reliability-based, effective firm capacity method, which we apply for multiple battery power ratings, durations, coupling types, deployment locations, and dispatch profiles within a test system that is based on the Texas Interconnection in the year 2024. We find that the capacity credits for all battery durations depend on their ability to predict the timing of reliability events. Even 1-2 h forecast errors - resulting in early or delayed battery discharging relative to the onset of a reliability event - lead pronounced capacity credit reductions, especially for 4-h duration batteries. Coupling batteries with solar mitigates the uncertainty associated with a shorter-duration battery's availability during reliability events, primarily due to the relatively high solar capacity credit in our test system. Coupled (or hybrid) system designs with oversized solar arrays, the ability to charge the coupled battery with grid energy, and larger batteries lead to the greatest capacity credit benefits of hybridization. We do not see evidence that the hybrid capacity credit exceeds the sum of the separate battery and solar capacity credits.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

On the Rate Constant for NH 2 +HO 2 and Third-Body Collision Efficiencies for NH 2 +H(+M) and NH 2 +NH 2 (+M)

In low-temperature flash photolysis of NH 3 /O 2 /N 2 mixtures, the NH 2 consumption rate and the product distribution is controlled by the reactions NH 2 + HO 2 → products (R1), NH 2 + H (+M) → NH 3 (+M) (R2), and NH 2 + NH 2 (+M) → N 2 H 4 (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH 2 + H (+N 2 ) and NH 2 + NH 2 (+N 2 ) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N 2 for R2 and R3 were calculated for O 2 and NH 3 as well as for other selected molecules. We report results were in good agreement with the limited experimental data. Based on reported NH 2 decay rates in flash photolysis of NH 3 /O 2 /N 2 , a rate constant for NH 2 + HO 2 → NH 3 + O 2 (R1a) of $k_{1\text{a}}$ = 1.5(±0.5) × 10 14 cm 3 mol –1 s –1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N 2 O yields indicates that NH 2 + HO 2 → H 2 NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH4 indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH 2 + HO 2 → HNO + H 2 O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate $k_{1\text{c}}$ = 2.5 × 10 13 cm 3 mol –1 s –1 and $k_{1\text{b}}$ = 2.5 × 10 12 cm 3 mol –1 s –1 at 295 K, with significant uncertainty margins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

To What Extent Will Decarbonization Deepen the Conversation Between Industry and the Grid?

Decarbonization - the transition away from un-mitigated fossil fuel combustion throughout the economy - requires big changes from both power and process systems. On the power system side, those changes are expected to include large increases in variable generation, e.g., from wind and solar, which has near-zero marginal costs and at large shares can produce infrequent but consequential energy droughts. On the process systems side, industries are investigating their options for direct and indirect electrification, the latter exemplified by replacing fossil fuel inputs with zero-carbon, energy-carrying chemicals like hydrogen and ammonia produced via electrochemical processes. The economic features of these changes within the larger context of power and process systems suggest that their realization could be accompanied by a paradigm shift in how industrial facilities interact with the grid. For example, the dominant type of demand participation in power markets could change from today's focus on load reductions at peak times to a new focus on shifting electricity use, enabled in part by large-scale product storage, to take advantage of renewable energy that would otherwise be curtailed and to avoid consumption during high-price energy droughts. This talk will describe these and other possible design and operational approaches from grid and industrial economic perspectives, culminating in an enumeration of open problems that lie at the interface of today and tomorrow's power and process systems.

co-design↗
DOE OSTI2021

Identification and Analysis of Backwater Nursery Habitats in the Middle Green River, Utah, during 2014, 2017, and 2018 Using High-Resolution Optical Remotely Sensed Imagery

Low-velocity channel-margin habitats, known as backwaters, serve as important nursery habitats for the endangered Colorado pikeminnow (Ptychocheilus lucius) in the middle Green River between Jensen and Ouray, Utah. The backwater synthesis report for the 1987–2013 period (Grippo et al. 2017) reported a decreased number of backwaters per river mile (RM) and an increased mean backwater size across the middle Green River. Information about backwaters is critical for understanding habitat characteristics that are important for Colorado pikeminnow recruitment. The goal of this study was to determine backwater number and size in the middle Green River using high-resolution imagery for 2018 to add to our understanding of long-term trends of backwater habitat availability across the reach. For comparative analysis of backwater habitats between 2004 and 2018, information for 2004, 2006, 2013, 2014 and 2017 from our previous studies was also utilized (Grippo et al. 2017; Hamada et al. 2017, 2021, and in review). Across the 2004, 2006, 2013, 2014, 2017, and 2018 study years, mean daily flow at the Jensen gage during image collection ranged from 1,220 cfs (2004) to 2,940 cfs (2014), which was approximately 140% greater in 2014 than in 2004. The number of delineated backwaters ranged from 88 (2013) to 245 (2017), and total backwater area ranged from 118,938 m 2 (2006) to 209,611 m 2 (2014). Mean backwater size ranged from 822 m 2 (2017) to 1,508 m 2 (2018), but there was considerable variability in backwater size for each of the study years. Backwater areas differed significantly among the study years based on a Kruskal-Wallis Rank Sum test ( P < 0.001). Pairwise Wilcoxon tests indicated that the median backwater areas in 2017 and 2014 were significantly lower than in the other sample years ( P ≤ 0.05), even though these years had the highest total backwater area. 2018 had the greatest median backwater area, although the distribution of backwater area values was not statistically different from backwater areas measured in 2013 ( P > 0.05). These two years also had the lowest maximum backwater area values among all sample years. The distribution of backwater area values was not significantly different ( P > 0.05) among the 2006, 2004, and 2013 sample years.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Swarm Flyby Gravimetry

This study describes a new technology for discerning the gravity fields and mass distribution of a solar system small body, without requiring dedicated orbiters or landers. Instead of a lander, a spacecraft releases a collection of small, simple probes during a flyby past an asteroid or comet. By tracking those probes from the host spacecraft, one can estimate the asteroid's gravity field and infer its underlying composition and structure. This approach offers a diverse measurement set,equivalent to planning and executing many independent and unique flyby encounters of a single spacecraft. This report assesses a feasible hardware implementation, derives the underlying models,and analyzes the performance of this concept via simulation.In terms of hardware, a small, low mass, low cost implementation is presented, which consists of a dispenser and probes. The dispenser contains roughly 12 probes in a tube and has a total size commensurate with a 6U P-Pod. The probes are housed in disc shaped sabots. When commanded,the dispenser ejects the top-most probe using a linear motor. The ejected probe separates from its sabots and unfolds using internal springs. There are two types of probes, each designed for a particular tracking modality. The reflective probe type, tracked by a telescope, unfolds to forma diffusely reflective sphere. The retroreflector probe type, tracked by a lidar, unfolds to form a corner-cube retroreflector assembly. Both types are designed to spherical so that their attitude doesn't affect the spacecraft's tracking performance.This analysis indicates that the point-mass term of small bodies larger than roughly 500 m in diameter can be observed from a host spacecraft that tracks locally deployed probes throughout a flyby to an uncertainty of better than 5%. The conditions by which this measurement is possible depends on the characteristics of the asteroid (size, type), the flyby velocity, and the type of tracking available (angles-only or angles+ranging). For most encounters, a few (1-3) well placed probes can be very effective, with marginal improvement for additional probes. Given realistic deployment errors, an encounter may require roughly 10-12 probes to ensure that 1-3 achieve their target. Long duration tracking of probes flying by large asteroids (>5 km diameter) can sometimes provide observability of the gravity field's first spherical harmonic, J( sub 2). In summary, this method offers a feasible, affordable approach to enabling or augmenting flyby science.

Atchison, Justin↗
DOE OSTI2023

ML-Based Power System Stability Assessment Considering Network Topology Changes: WECC 20,000+ Bus System Case Study

Modern power grids are fast-changing and thus require real-time monitoring and online stability assessment. With the rapid development of machine learning (ML) techniques, using data-driven models to provide fast and accurate estimations of power system stability marginal information, such as frequency nadir for frequency stability and critical clearing time (CCT) for transient stability, have become possible. However, despite the numerous research on ML-based methods for frequency nadir and CCT prediction, there is limited work on the impact of different network topology changes. Furthermore, most previous studies only focused on small or synthetic systems, and there is a lack of research on actual large power system models. In this paper, the above issues are addressed by studying the actual U.S. Western Electricity Coordinating Council (WECC) system model with more than 20,000 buses. Massive simulations are conducted in PowerWorld Simulator to study the impact of various topology change scenarios on both frequency stability and transient stability. System operating information is extracted from the success dispatch cases of various network topologies to generate a comprehensive dataset for ML-based models. Two ML methods, random forest (RF) and multilayer perceptron (MLP) neural network, are trained and tested for both frequency nadir prediction and CCT prediction. Test results have proven the models are capable of online stability assessment for large power networks such as the WECC system with sufficient accuracy.

critical clearing time↗
NASA NTRS2011

The Impact of Rise of the Andes and Amazon Landscape Evolution on Diversification of Lowland terra-firme Forest Birds

Since the 19th Century, the unmatched biological diversity of Amazonia has stimulated a diverse set of hypotheses accounting for patterns of species diversity and distribution in mega-diverse tropical environments. Unfortunately, the evidence supporting particular hypotheses to date is at best described as ambiguous, and no generalizations have emerged yet, mostly due to the lack of comprehensive comparative phylogeographic studies with thorough trans-Amazonian sampling of lineages. Here we report on spatial and temporal patterns of diversification estimated from mitochondrial gene trees for 31 lineages of birds associated with upland terra-firme forest, the dominant habitat in modern lowland Amazonia. The results confirm the pervasive role of Amazonian rivers as primary barriers separating sister lineages of birds, and a protracted spatio-temporal pattern of diversification, with a gradual reduction of earlier (1st and 2nd) and older (> 2 mya) splits associated with each lineage in an eastward direction. (The easternmost tributaries of the Amazon, the Xingu and Tocantins Rivers, are not associated with any splits older than > 2 mya). For the suboscine passerines, maximum-likelihood estimates of rates of diversification point to an overall constant rate over the past 5 my (up to a significant downturn at 300,000 y ago). This "younging-eastward" pattern may have an abiotic explanation related to landscape evolution. Triggered by a new pulse of Andean uplift, it has been proposed that modern Amazon basin landscapes may have evolved successively eastward, away from the mountain chain, starting approximately 10 mya. This process was likely based on the deposition of vast fluvial sediment masses, known as megafans, that may have extended progressively and in series eastward from Andean sources. This process plausibly explains the progressive extinction of original Pebas wetland of western-central Amazonia by the present fluvial landsurfaces of a more terra-firme type. The youngest landsurfaces thus lie furthest from the mountains. In this scenario major drainages were also reoriented in wholesale fashion away from a northerly orientation generally towards the east and an Atlantic Ocean outlet. The advance of megafans is best seen by the location of axial rivers such as the Orinoco and Mamore which lie against the cratonic margins furthest from the Andes, at the distal ends of major megafan ramparts. More importantly, other major river courses in western-central Amazonia will have been established at progressively younger dates with distance eastward. If this landscape-sequence scenario is accurate, it parallels the progressive younging of the passerine lineages. The bird DNA data appears to confirm strongly the pervasive role of Amazonian rivers--as primary barriers separating sister lineages of birds, and thus probably as facilitaters of bird speciation. We show for the first time that a general spatio-temporal pattern of diversification for terra-firme lineages in the Amazon is associated with rivers ("younging-eastward"), and furthermore parallels a specific scenario of regional drainage evolution.

Aleixo, Alexandre↗
NASA NTRS2016

Soil Moisture Active Passive Mission L4_SM Data Product Assessment (Version 2 Validated Release)

During the post-launch SMAP calibration and validation (Cal/Val) phase there are two objectives for each science data product team: 1) calibrate, verify, and improve the performance of the science algorithm, and 2) validate the accuracy of the science data product as specified in the science requirements and according to the Cal/Val schedule. This report provides an assessment of the SMAP Level 4 Surface and Root Zone Soil Moisture Passive (L4_SM) product specifically for the product's public Version 2 validated release scheduled for 29 April 2016. The assessment of the Version 2 L4_SM data product includes comparisons of SMAP L4_SM soil moisture estimates with in situ soil moisture observations from core validation sites and sparse networks. The assessment further includes a global evaluation of the internal diagnostics from the ensemble-based data assimilation system that is used to generate the L4_SM product. This evaluation focuses on the statistics of the observation-minus-forecast (O-F) residuals and the analysis increments. Together, the core validation site comparisons and the statistics of the assimilation diagnostics are considered primary validation methodologies for the L4_SM product. Comparisons against in situ measurements from regional-scale sparse networks are considered a secondary validation methodology because such in situ measurements are subject to up-scaling errors from the point-scale to the grid cell scale of the data product. Based on the limited set of core validation sites, the wide geographic range of the sparse network sites, and the global assessment of the assimilation diagnostics, the assessment presented here meets the criteria established by the Committee on Earth Observing Satellites for Stage 2 validation and supports the validated release of the data. An analysis of the time average surface and root zone soil moisture shows that the global pattern of arid and humid regions are captured by the L4_SM estimates. Results from the core validation site comparisons indicate that "Version 2" of the L4_SM data product meets the self-imposed L4_SM accuracy requirement, which is formulated in terms of the ubRMSE: the RMSE (Root Mean Square Error) after removal of the long-term mean difference. The overall ubRMSE of the 3-hourly L4_SM surface soil moisture at the 9 km scale is 0.035 cubic meters per cubic meter requirement. The corresponding ubRMSE for L4_SM root zone soil moisture is 0.024 cubic meters per cubic meter requirement. Both of these metrics are comfortably below the 0.04 cubic meters per cubic meter requirement. The L4_SM estimates are an improvement over estimates from a model-only SMAP Nature Run version 4 (NRv4), which demonstrates the beneficial impact of the SMAP brightness temperature data. L4_SM surface soil moisture estimates are consistently more skillful than NRv4 estimates, although not by a statistically significant margin. The lack of statistical significance is not surprising given the limited data record available to date. Root zone soil moisture estimates from L4_SM and NRv4 have similar skill. Results from comparisons of the L4_SM product to in situ measurements from nearly 400 sparse network sites corroborate the core validation site results. The instantaneous soil moisture and soil temperature analysis increments are within a reasonable range and result in spatially smooth soil moisture analyses. The O-F residuals exhibit only small biases on the order of 1-3 degrees Kelvin between the (re-scaled) SMAP brightness temperature observations and the L4_SM model forecast, which indicates that the assimilation system is largely unbiased. The spatially averaged time series standard deviation of the O-F residuals is 5.9 degrees Kelvin, which reduces to 4.0 degrees Kelvin for the observation-minus-analysis (O-A) residuals, reflecting the impact of the SMAP observations on the L4_SM system. Averaged globally, the time series standard deviation of the normalized O-F residuals is close to unity, which would suggest that the magnitude of the modeled errors approximately reflects that of the actual errors. The assessment report also notes several limitations of the "Version 2" L4_SM data product and science algorithm calibration that will be addressed in future releases. Regionally, the time series standard deviation of the normalized O-F residuals deviates considerably from unity, which indicates that the L4_SM assimilation algorithm either over- or under-estimates the actual errors that are present in the system. Planned improvements include revised land model parameters, revised error parameters for the land model and the assimilated SMAP observations, and revised surface meteorological forcing data for the operational period and underlying climatological data. Moreover, a refined analysis of the impact of SMAP observations will be facilitated by the construction of additional variants of the model-only reference data. Nevertheless, the “Version 2” validated release of the L4_SM product is sufficiently mature and of adequate quality for distribution to and use by the larger science and application communities.

SMAP L4_SM↗
NASA NTRS1995

Determination of Eros Physical Parameters for Near Earth Asteroid Rendezvous Orbit Phase Navigation

Navigation of the orbit phase of the Near Earth steroid Rendezvous (NEAR) mission will re,quire determination of certain physical parameters describing the size, shape, gravity field, attitude and inertial properties of Eros. Prior to launch, little was known about Eros except for its orbit which could be determined with high precision from ground based telescope observations. Radar bounce and light curve data provided a rough estimate of Eros shape and a fairly good estimate of the pole, prime meridian and spin rate. However, the determination of the NEAR spacecraft orbit requires a high precision model of Eros's physical parameters and the ground based data provides only marginal a priori information. Eros is the principal source of perturbations of the spacecraft's trajectory and the principal source of data for determining the orbit. The initial orbit determination strategy is therefore concerned with developing a precise model of Eros. The original plan for Eros orbital operations was to execute a series of rendezvous burns beginning on December 20,1998 and insert into a close Eros orbit in January 1999. As a result of an unplanned termination of the rendezvous burn on December 20, 1998, the NEAR spacecraft continued on its high velocity approach trajectory and passed within 3900 km of Eros on December 23, 1998. The planned rendezvous burn was delayed until January 3, 1999 which resulted in the spacecraft being placed on a trajectory that slowly returns to Eros with a subsequent delay of close Eros orbital operations until February 2001. The flyby of Eros provided a brief glimpse and allowed for a crude estimate of the pole, prime meridian and mass of Eros. More importantly for navigation, orbit determination software was executed in the landmark tracking mode to determine the spacecraft orbit and a preliminary shape and landmark data base has been obtained. The flyby also provided an opportunity to test orbit determination operational procedures that will be used in February of 2001. The initial attitude and spin rate of Eros, as well as estimates of reference landmark locations, are obtained from images of the asteroid. These initial estimates are used as a priori values for a more precise refinement of these parameters by the orbit determination software which combines optical measurements with Doppler tracking data to obtain solutions for the required parameters. As the spacecraft is maneuvered; closer to the asteroid, estimates of spacecraft state, asteroid attitude, solar pressure, landmark locations and Eros physical parameters including mass, moments of inertia and gravity harmonics are determined with increasing precision. The determination of the elements of the inertia tensor of the asteroid is critical to spacecraft orbit determination and prediction of the asteroid attitude. The moments of inertia about the principal axes are also of scientific interest since they provide some insight into the internal mass distribution. Determination of the principal axes moments of inertia will depend on observing free precession in the asteroid's attitude dynamics. Gravity harmonics are in themselves of interest to science. When compared with the asteroid shape, some insight may be obtained into Eros' internal structure. The location of the center of mass derived from the first degree harmonic coefficients give a direct indication of overall mass distribution. The second degree harmonic coefficients relate to the radial distribution of mass. Higher degree harmonics may be compared with surface features to gain additional insight into mass distribution. In this paper, estimates of Eros physical parameters obtained from the December 23,1998 flyby will be presented. This new knowledge will be applied to simplification of Eros orbital operations in February of 2001. The resulting revision to the orbit determination strategy will also be discussed.

Miller, J. K.↗
NASA NTRS2003

Platinum-group elements (PGE) and Rhenium in Marine Sediments across the Cretaceous-Tertiary Boundary: Constraints on Re-PGE Transport in the Marine Environment

The nature of Re-platinum-group element (PGE; Pt, Pd, Ir, Os, Ru) transport in the marine environment was investigated by means of marine sediments at and across the Cretaceous-Tertiary boundary (KTB) at two hemipelagic sites in Europe and two pelagic sites in the North and South Pacific. A traverse across the KTB in the South Pacific pelagic clay core found elevated levels of Re, Pt, Ir, Os, and Ru, each of which is approximately symmetrically distributed over a distance of approx. 1.8 m across the KTB. The Re-PGE abundance patterns are fractionated from chondritic relative abundances: Ru, Pt, Pd, and Re contents are slightly subchondritic relative to Ir, and Os is depleted by approx. 95% relative to chondritic Ir proportions. A similar depletion in Os (approx. 90%) was found in a sample of the pelagic KTB in the North Pacific, but it is enriched in Ru, Pt, Pd, and Re relative to Ir. The two hemipelagic KTB clays have near-chondritic abundance patterns. The approx. 1.8-m-wide Re-PGE peak in the pelagic South Pacific section cannot be reconciled with the fallout of a single impactor, indicating that postdepositional redistribution has occurred. The elemental profiles appear to fit diffusion profiles, although bioturbation could have also played a role. If diffusion had occurred over approx. 65 Ma, the effective diffusivities are approx. 10(exp -13)sq cm/s, much smaller than that of soluble cations in pore waters (approx. 10(exp -5) sq cm/s). The coupling of Re and the PGEs during redistribution indicates that postdepositional processes did not significantly fractionate their relative abundances. If redistribution was caused by diffusion, then the effective diffusivities are the same. Fractionation of Os from Ir during the KTB interval must therefore have occurred during aqueous transport in the marine environment. Distinctly subchondritic Os/Ir ratios throughout the Cenozoic in the South Pacific core further suggest that fractionation of Os from Ir in the marine environment is a general process throughout geologic time because most of the inputs of Os and Ir into the ocean have OsAr ratios greater than or = 1. Mass balance calculations show that Os and Re burial fluxes in pelagic sediments account for only a small fraction of the riverine Os (less than 10%) and Re (less than 0.1%) inputs into the oceans. In contrast, burial of Ir in pelagic sediments is similar to the riverine Ir input, indicating that pelagic sediments are a much larger repository for Ir than for Os and Re. If all of the missing Os and Re is assumed to reside in anoxic sediments in oceanic margins, the calculated burial fluxes in anoxic sediments are similar to observed burial fluxes. However, putting all of the missing Os and Re into estuarine sediments would require high concentrations to balance the riverine input and would also fail to explain the depletion of Os at pelagic KTB sites, where at most approx. 25% of the K-T impactor's Os could have passed through estuaries. If Os is preferentially sequestered in anoxic marine environments, it follows that the OsAr ratio of pelagic sediments should be sensitive to changes in the rates of anoxic sediment deposition. There is thus a clear fractionation of Os and Re from Ir in precipitation out of sea water in pelagic sections. Accordingly, it is inferred here that Re and Os are removed from sea water in anoxic marine depositional regimes.

Lee, Cin-Ty Aeolus↗