Search NASA⌕ Search

SEARCH · Search NASA

Results for “constrained variational analysis”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

The Nanograv Nine-Year Data Set: Measurement and Analysis of Variations in Dispersion Measures

We analyze dispersion measure(DM) variations of 37 millisecond pulsars in the nine-year North American Nanohertz Observatory for Gravitational Waves (NANOGrav) data release and constrain the sources of these variations. DM variations can result from a changing distance between Earth and the pulsar, inhomogeneities in the interstellar medium, and solar effects. Variations are significant for nearly all pulsars, with characteristic timescales comparable to or even shorter than the average spacing between observations. Five pulsars have periodic annual variations, 14 pulsars have monotonically increasing or decreasing trends, and 14 pulsars show both effects. Of the four pulsars with linear trends that have line-of-sight velocity measurements, three are consistent with a changing distance and require an overdensity of free electrons local to the pulsar. Several pulsars show correlations between DM excesses and lines of sight that pass close to the Sun. Mapping of the DM variations as a function of the pulsar trajectory can identify localized interstellar medium features and, in one case, an upper limit to the size of the dispersing region of 4 au. Four pulsars show roughly Kolmogorov structure functions (SFs), and another four show SFs less steep than Kolmogorov. One pulsar has too large an uncertainty to allow comparisons. We discuss explanations for apparent departures from a Kolmogorov-like spectrum, and we show that the presence of other trends and localized features or gradients in the interstellar medium is the most likely cause.

pulsars: general↗
NASA NTRSDate not supplied

Contribution to Solar Brightness of Small-Size Magnetic Elements

Irradiance variability is mostly driven by surface magnetism, each magnetic feature contributing in a different manner. The contribution of small-size magnetic elements observed ubiquitously at high-spatial resolution on the solar photosphere, is still debated, as such features are mostly unresolved on full-disk images employed to model irradiance variability. Understanding the contribution to solar brightness of small-size magnetic elements, especially in quiet regions, will help understanding irradiance variability on the decadal and longer temporal scales, which, in turn, are fundamental to understand the role of the Sun on the Earth’s climate. We present a preliminary study of the brightness of small-size magnetic elements using high-spatial resolution spectro-polarimetric observations acquired with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) during Commissioning Phase 1. We focus on the Intensity contrast vs magnetic field relation, which is a fundamental observable in irradiance studies. Previous studies mostly focused on the FeI 630.1/630.2 nm spectral range and so to contextualize our results we focus here on this spectral region as well. By comparing DKIST results with MURaM magneto-hydrodynamic simulations, and with previous results published in the literature, we conclude that our results are in agreement with studies conducted on HINODE observations. We plan to extend the analysis to UV and IR ViSP observations and VBI blue continuum in order to better constrain the quiet sun contribution to spectral irradiance variations.

SMD↗
DOE OSTI2020

Enriched East Asian oxygen isotope of precipitation indicates reduced summer seasonality in regional climate and westerlies

Speleothem oxygen isotope records over East Asia reveal apparently large and rapid paleoclimate changes over the last several hundred thousand years. However, what the isotopic variation actually represent in terms of the regional climate and circulation is debated. We present an answer that emerges from an analysis of the interannual variation in amount-weighted annual δ 18 O of precipitation over East Asia as simulated by an isotope-enabled model constrained by large-scale atmospheric reanalysis fields. 18 O-enriched years have reduced summer seasonality both in terms of precipitation isotopes and in the large-scale circulation. Changes occur between June and October, where the δ 18 O of precipitation (δ 18 O p ) transitions from the isotopically heavier winter to the lighter summer regime. For 18 O-enriched years, this transition is less pronounced. Variations in precipitation amount alone are insufficient to explain the amount-weighted annual δ 18 O p between 18 O-enriched and 18 O-depleted years. Reduced summer seasonality is also expressed in the low-level monsoonal southerlies and upper-level westerlies; for the latter, the northward migration across the Tibetan Plateau in the summer is less pronounced. Our result thus implicates the westerlies across the plateau as the proximate cause of East Asian paleomonsoon changes, manifested as a modulation of its summer peak.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Seasonal variations of snow depth on Mars

Using topography collected over one martian year from the Mars Orbiter Laser Altimeter on the Mars Global Surveyor (MGS) spacecraft, we have measured temporal changes in the elevation of the martian surface that correlate with the seasonal cycle of carbon dioxide exchange between the surface and atmosphere. The greatest elevation change (1.5 to 2 meters) occurs at high latitudes ( above 80 degrees ), whereas the bulk of the mass exchange occurs at lower latitudes (below 75 degrees N and below 73 degrees S). An unexpected period of sublimation was observed during northern hemisphere autumn, coincident with dust storms in the southern hemisphere. Analysis of MGS Doppler tracking residuals revealed temporal variations in the flattening of Mars that correlate with elevation changes. The combined changes in gravity and elevation constrain the average density of seasonally deposited carbon dioxide to be 910 +/- 230 kilograms per cubic meter, which is considerably denser than terrestrial snow.

Carbon Dioxide↗
DOE OSTI2021

A Constrained Data Assimilation Algorithm Based on GSI Hybrid 3D-EnVar and Its Application

Abstract Data assimilation (DA) at mesoscales is important for severe weather forecasts, yet the techniques of data assimilation at this scale remain a challenge. This study introduces dynamical constraints in the Gridpoint Statistical Interpolation (GSI) three-dimensional ensemble variational (3D-EnVar) data assimilation algorithm to enable the use of high-resolution surface observations of precipitation to improve atmospheric analysis at mesoscales. The constraints use the conservations of mass and moisture. Mass constraint suppresses the unphysical high-frequency oscillation, while moisture conservation constrains the atmospheric states to conform with the observed high-resolution precipitation. We show that the constrained data assimilation (CDA) algorithm significantly reduced the spurious residuals of the mass and moisture budgets compared to the original data assimilation (ODA). A case study is presented for a squall line over the Southern Great Plains on 20 May 2011 during Midlatitude Continental Convective Clouds Experiment (MC3E) of the Atmospheric Radiation Measurement (ARM) program by using ODA or CDA analysis as initial condition of forecasts. The state variables, and the location and intensity of the squall line are better simulated in the CDA experiment. Results show how surface observation of precipitation can be used to improve atmospheric analysis through data assimilation by using the dynamical constraints of mass and moisture conservations.

Wang, Jia↗
NASA NTRS2007

Evidence for Interannual to Decadal Variations in Hadley and Walker Circulations and Links to Water and Energy Fluxes

Mass and energy transports associated with the Hadley and Walker circulations are important components of the earth s climate system and are strongly linked to hydrologic processes. Interannual to decadal variation in these flows likely signify a combination of natural climate noise as well as a response to anthropgenic forcing. There remains considerable uncertainty in quantifying variations in these flows. Evidence in the surface pressure record supports a weakening of the Walker circulation over the Pacific in recent decades. Conversely the NCEP / NCAR and ERA 40 reanalyses indicate that the Hadley circulation has increased in strength over the last two decades, though these analyses depict significantly different mass circulation changes. Interestingly, the NCEP - II / DOE reanalysis contains essentially no Hadley circulation changes. Most climate model integrations anticipate a weakening of both tropical circulations associated with stronger static stability. Clearly there is much uncertainty not only with the mass transports, but also how they are linked to water and energy balance of the planet through variations in turbulent heat and radiative fluxes and horizontal exports / imports of energy. Here we examine heat and water budget variations from a number of reanalysis products and focus on the linear and nonlinear response of ENSO warm and cold events as opportunities to study budget variations over the past 15-20 years. Our analysis addresses such questions as To what extent do Hadley and Walker Cell variations compensate each other on mass and energy transport? Do static stability adjustments appear to constrain fractional precipitation response vs. fractional water vapor response? We appeal to constraints offered by GPCP precipitation, SSWI ocean evaporation estimates, and ISCCP-FD radiative fluxes, and other satellite data sets to interpret and confirm reanalysis-based diagnostics. Using our findings we also attempt to place in context the recent findings that tropical ocean evaporation increased by order 5% or more during the 1990s, reconciling this with GPCP precipitation variations.

Robertson, Franklin↗
NASA NTRS2011

Design and Evaluation of the Terminal Area Precision Scheduling and Spacing System

This paper describes the design, development and results from a high fidelity human-in-the-loop simulation of an integrated set of trajectory-based automation tools providing precision scheduling, sequencing and controller merging and spacing functions. These integrated functions are combined into a system called the Terminal Area Precision Scheduling and Spacing (TAPSS) system. It is a strategic and tactical planning tool that provides Traffic Management Coordinators, En Route and Terminal Radar Approach Control air traffic controllers the ability to efficiently optimize the arrival capacity of a demand-impacted airport while simultaneously enabling fuel-efficient descent procedures. The TAPSS system consists of four-dimensional trajectory prediction, arrival runway balancing, aircraft separation constraint-based scheduling, traffic flow visualization and trajectory-based advisories to assist controllers in efficient metering, sequencing and spacing. The TAPSS system was evaluated and compared to today's ATC operation through extensive series of human-in-the-loop simulations for arrival flows into the Los Angeles International Airport. The test conditions included the variation of aircraft demand from a baseline of today's capacity constrained periods through 5%, 10% and 20% increases. Performance data were collected for engineering and human factor analysis and compared with similar operations both with and without the TAPSS system. The engineering data indicate operations with the TAPSS show up to a 10% increase in airport throughput during capacity constrained periods while maintaining fuel-efficient aircraft descent profiles from cruise to landing.

Swenson, Harry N.↗
NASA NTRS2001

Characteristic Ultraviolet/Optical Timescales in Active Galactic Nuclei

The UV/optical light curves of active galactic nuclei (AGNs) may be roughly characterized by a periodic fractional flux variations of approx. 10% on timescales of 1 month. The physical mechanism(s) responsible remain ill defined. We present a structure function analysis, i.e., measure the power distribution over a range of timescales tau, of 13 AGNs to constrain the origin of UV/optical emission. On timescales tau approx. 5- 60 days, the mean UV and optical power density spectra (PDS) are equivalent. This may suggest that the underlying energy generating mechanism is identical. The combined W/optical PDS is P(f) proportional to f(sup -alpha)- with alpha = 2.13(sup +0.22, sub -0.06). For sources with measured X-ray PDS indices, we find they are indistinguishable from their UV/optical counterparts. This supports scenarios whereby X-rays are generated via Compton upscattering of UV photons, to later radiatively drive optical variations. At the same time, we present evidence for characteristic variability timescales tau(sub char) of approx. 5-100 days in 10 sources. These variability timescales combined with reverberation based masses M suggest a M - tau(sub char) relationship; higher mass systems have larger characteristic timescales. The UV tau(sub char) may possibly reflect dynamical or accretion disk thermal timescales. We find suggestive evidence for a dichotomy, at tau approx. 30 days and M approx. l0(exp 7) solar mass, between short- and long-time scale optical variations. These optical variations may be attributable to dynamical and accretion disk thermal or starburst activity timescales, respectively.

Collier, Stefan↗
DOE OSTI2017

Assimilation of repeated woody biomass observations constrains decadal ecosystem carbon cycle uncertainty in aggrading forests

Forest carbon sink strengths are governed by plant growth, mineralization of dead organic matter, and disturbance. Across landscapes, remote sensing can provide information about aboveground states of forests and this information can be linked to models to estimate carbon cycling in forests close to steady state. For aggrading forests this approach is more challenging and has not been demonstrated. Here in this paper we apply a Bayesian approach, linking a simple model to a range of data, to evaluate their information content, for two aggrading forests. We compare high information content analyses using local observations with retrievals using progressively sparser remotely sensed information (repeated, single, and no woody biomass observations). The net biome productivity of both forests is constrained to be a net sink with <2 Mg C ha -1 yr -1 variation across the range of inputs. However, the sequestration of particular carbon pool(s) varies with assimilated biomass information. Assimilation of repeated biomass observations reduces uncertainty and/or bias in all ecosystem C pools not just wood, compared to analyses using single or no stock information. As verification, our repeated biomass analysis explains 78–86% of variation in litter dynamics at one forest, while at the second forest total dead organic matter estimates are within observational uncertainty. The uncertainty of retrieved ecosystem traits in the repeated biomass analysis is reduced by up to 50% compared to analyses with less biomass information. This study quantifies the importance of repeated woody observations in constraining the dynamics of both wood and dead organic matter, highlighting the benefit of proposed remote sensing missions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Data and Code for: Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits

This repository contains the simulation outputs and processing scripts associated with the study of winter wheat traits across the United States, utilizing the Ecosys agroecosystem model. The dataset includes model results for both rainfed and irrigated winter wheat systems, supporting the findings presented in the manuscript titled "Observation-constrained agroecosystem model inversion reveals continental-scale variation of winter wheat traits." Data includes the original Ecosys simulation outputs (archived in .db format within the compressed .zip files) and extracted analysis data (stored in .pkl files for efficient processing). Python code for data processing and figure generation is provided in a Jupyter notebook. External Observational Datasets should refer to the following official repositories for the input and validation data used in this study. The eddy covariance data from the AmeriFlux network (https://ameriflux.lbl.gov/). Climate-forcing data of NLDAS-2 from NASA LDAS (https://ldas.gsfc.nasa.gov/nldas/nldas-2-forcing-data). Soil data from the Gridded Soil Survey Geographic Database (gSSURGO), available at (https://www.nrcs.usda.gov/resources/data-and-reports/gridded-soil-survey-geographic-gssurgo-database). Crop yields, planting and harvest dates from the USDA public databases (https://quickstats.nass.usda.gov/; https://webapp.rma.usda.gov/apps/actuarialinformationbrowser/CropCriteria.aspx). Satellite-derived SLOPE GPP data from ORNL DAAC (https://daac.ornl.gov/cgi-bin/dsviewer.pl?ds_id=1786). Land use and crop progress information from the USDA Crop Data Layer and Crop Progress and Condition Gridded Layers (https://www.nass.usda.gov/Research_and_Science/). The Ecosys model code is available online at https://github.com/jinyun1tang/ECOSYS.

Wheat↗
DOE OSTI2021

Variational Quantum Algorithms for Semidefinite Programming

A semidefinite program (SDP) is a particular kind of convex optimization problem with applications in operations research, combinatorial optimization, quantum information science, and beyond. In this work, we propose variational quantum algorithms for approximately solving SDPs. For one class of SDPs, we provide a rigorous analysis of their convergence to approximate locally optimal solutions, under the assumption that they are weakly constrained (i.e., N " M, where N is the dimension of the input matrices and M is the number of constraints). We also provide algorithms for a more general class of SDPs that requires fewer assumptions. Finally, we numerically simulate our quantum algorithms for applications such as MaxCut, and the results of these simulations provide evidence that convergence still occurs in noisy settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Variational Quantum Algorithms for Semidefinite Programming

A semidefinite program (SDP) is a particular kind of convex optimization problem with applications in operations research, combinatorial optimization, quantum information science, and beyond. In this work, we propose variational quantum algorithms for approximately solving SDPs. For one class of SDPs, we provide a rigorous analysis of their convergence to approximate locally optimal solutions, under the assumption that they are weakly constrained (i.e., N$\gg$M, where N is the dimension of the input matrices and M is the number of constraints). We also provide algorithms for a more general class of SDPs that requires fewer assumptions. Finally, we numerically simulate our quantum algorithms for applications such as MaxCut, and the results of these simulations provide evidence that convergence still occurs in noisy settings.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1998

Magnetospheric Image Unfolding

The Grant was a three year grant funded under the Space Physics Supporting Research and Technology and Suborbital Program. Our objective was to develop automated techniques needed to unfold or "invert" global images of the magnetospheric ion populations obtained by the new magnetospheric imaging techniques (ENA, EUV) in anticipation of future missions such as the Magnetospheric Imager and, now, IMAGE. Our focus on the present three year grant is to determine the degree to which such images can quantitatively constrain the global electromagnetic properties of the magnetosphere. In a previous three year grant period we successfully automated a forward modeling inversion algorithm, demonstrated that these inversions are robust in the face of realistic instrumental considerations such as counting statistics and backgrounds, applied error analysis techniques to the extracted parameters using variational procedures, implemented very realistic magnetospheric test images to test the inversion algorithms using the Rice University Magnetospheric Specification Model, and began the process of generating parametric models with the flexibility to handle the realistic magnetospheric images (e.g. Roelof et al, 1992; 1993). Our plan for the present 3 year grant period was to complete the development of the inversion tools needed to handle realistic magnetospheric images, assess the degree to which global electrodynamics is quantitatively constrained by ENA images of the magnetosphere, and bring the inversion of EUV images up to the maturity that we will have achieved for the ENA imaging. Below the accomplishments of our three year effort are present followed by a list of our presentations and publications. The accomplishments of all three years are presented here, and thus some of these items appeared on interim progress reports.

Source record↗
DOE OSTI2017

REFINED SYSTEM PARAMETERS AND TTV STUDY OF TRANSITING EXOPLANETARY SYSTEM HAT-P-20

We report new photometric observations of the transiting exoplanetary system HAT-P-20, obtained using CCD cameras at Yunnan Observatories and Ho Koon Nature Education cum Astronomical Centre, China, from 2010 to 2013, and Observatori Ca l’Ou, Sant Marti Sesgueioles, Spain, from 2013 to 2015. The observed data are corrected for systematic errors according to the coarse de-correlation and SYSREM algorithms, so as to enhance the signal of the transit events. In order to consistently model the star spots and transits of this exoplanetary system, we develop a highly efficient tool STMT based on the analytic models of Mandel and Agol and Montalto et al. The physical parameters of HAT-P-20 are refined by homogeneously analyzing our new data, the radial velocity data, and the earlier photometric data in the literature with the Markov chain Monte Carlo technique. New radii and masses of both host star and planet are larger than those in the discovery paper due to the discrepancy of the radius among K-dwarfs between predicted values by standard stellar models and empirical calibration from observations. Through the analysis of all available mid-transit times calculated with the normal model and spotted model, we conclude that the periodic transit timing variations in these transit events revealed by employing the normal model are probably induced by spot crossing events. From the analysis of the distribution of occulted spots by HAT-P-20b, we constrain the misaligned architecture between the planetary orbit and the spin of the host star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Spectral characteristics of background error covariance and multiscale data assimilation

The spatial resolutions of numerical atmospheric and oceanic circulation models have steadily increased over the past decades. Horizontal grid spacing down to the order of 1 km is now often used to resolve cloud systems in the atmosphere and sub-mesoscale circulation systems in the ocean. These fine resolution models encompass a wide range of temporal and spatial scales, across which dynamical and statistical properties vary. In particular, dynamic flow systems for small scales can become spatially localized and temporarily intermittent. An analysis shows that the background correlation length scale is larger than 75 km for streamfunctions, even for a 2-km resolution model, and larger than 25 km for water vapor mixing ratios. The theoretical analyses suggest that such correlation length scales prevent the currently used data assimilation schemes from constraining spatial scales smaller than 150 km for streamfunctions and 50 km for water vapor mixing ratios. These results highlight the necessity of fundamentally modifying the currently used data assimilation algorithm for assimilating high-resolution observations into the aforementioned fine resolution models. A multiscale methodology based on scale decomposition is suggested, and challenges are discussed.

Li, Zhijin↗
NASA NTRS2008

Modeling the Spin Motor Current of the International Space Station's Control Moment Gyroscopes

The International Space Station (ISS) attitude control is provided by two means: The Russian Segment uses thrusters and the U.S. Segment uses double-gimbaled control moment gyroscopes (CMG). CMGs are used as momentum exchange devices, providing non propulsive attitude control for the vehicle. The CMGs are very important for the ISS program because, first, they save propellant - which needs to be transferred to the Station in special cargo vehicles - and, second, they provide the microgravity environment on the Station - which is necessary for scientific experiments planned for the ISS mission. Since 2002, when one of the CMG on the ISS failed, all CMGs are closely monitored. High gimbal rates, vibration spikes, unusual variations of spin motor current and bearing temperatures are of great concern, since these parameters are the CMG health indicators. The telemetry analysis of these and some other CMG parameters is used to determine constrains and make changes to the CMGs operation on board. These CMG limitations, in turn, may limit the ISS attitude control capabilities and may be critical to ISS operation. Therefore, it is important to know whether the CMG parameter is nominal or out of family, and why. The goal of this project is to analyze an important CMG parameter - spin motor current. Some operational decisions are made now based on the spin motor current signatures. The spin motor current depends on gimbal rates, ISS rates, and spin bearing friction. The spin bearing friction in turn depends on the bearing temperatures, wheel rates, normal load - which is a function of gimbal and wheel rates - lubrication, etc. The first task of this project is to create a spin motor current mathematical model based on CMG dynamics model and the current knowledge on bearing friction in microgravity.

Pereira, Miguel A.↗
DOE OSTI2022

Misinterpreting modified gravity as dark energy: a quantitative study

Standard cosmological data analyses typically constrain simple phenomenological dark-energy parameters, for example the present-day value of the equation of state parameter, w 0 , and its variation with scale factor, w a . However, results from such an analysis cannot easily indicate the presence of modified gravity. Even if general relativity does not hold, experimental data could still be fit sufficiently well by a phenomenological w 0 w a CDM, unmodified-gravity model. Hence, it would be very useful to know if there are generic signatures of modified gravity in standard analyses. Here we present, for the first time to our knowledge, a quantitative mapping showing how modified gravity models look when (mis)interpreted within the standard unmodified-gravity analysis. Scanning through a broad space of modified-gravity (Horndeski) models, and assuming a near-future survey consisting of CMB, BAO, and SNIa observations, we report values of the best-fit set of cosmological parameters including (w 0 , w a ) that would be inferred if modified gravity were at work. We find that modified gravity models that can masquerade as standard gravity lead to very specific biases in standard-parameter spaces. Here, we also comment on implications for measurements of the amplitude of mass fluctuations described by the parameter S 8 .

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Variational four-dimensional analysis using quasi-geostrophic constraints

A variational four-dimensional analysis technique using quasi-geostrophic models as constraints is examined using gridded fields as data. The analysis method uses a standard iterative nonlinear minimization technique to find the solution to the constraining forecast model which best fits the data as measured by a predefined functional. The minimization algorithm uses the derivative of the functional with respect to each of the initial condition values. This derivative vector is found by inserting the weighted differences between the model solution and the inserted data into a backwards integrating adjoint model. The four-dimensional analysis system was examined by applying it to fields created from a primitive equations model forecast and to fields created from satellite retrievals. The results show that the technique has several interesting characteristics not found in more traditional four-dimensional assimilation techniques. These features include a close fit of the model solution to the observations throughout the analysis interval and an insensitivity to the frequency of data insertion or the amount of data. The four-dimensional analysis technique is very versatile and can be extended to more complex problems with little theoretical difficulty.

Derber, John C.↗