Search NASA⌕ Search

SEARCH · Search NASA

Results for “bivariate”

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.

160 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Optimal Initial Perturbations for Ensemble Prediction of the Madden-Julian Oscillation during Boreal Winter

An initialization strategy, tailored to the prediction of the Madden-Julian oscillation (MJO), is evaluated using the Goddard Earth Observing System Model, version 5 (GEOS-5), coupled general circulation model (CGCM). The approach is based on the empirical singular vectors (ESVs) of a reduced-space statistically determined linear approximation of the full nonlinear CGCM. The initial ESV, extracted using 10 years (1990-99) of boreal winter hindcast data, has zonal wind anomalies over the western Indian Ocean, while the final ESV (at a forecast lead time of 10 days) reflects a propagation of the zonal wind anomalies to the east over the Maritime Continent an evolution that is characteristic of the MJO. A new set of ensemble hindcasts are produced for the boreal winter season from 1990 to 1999 in which the leading ESV provides the initial perturbations. The results are compared with those from a set of control hindcasts generated using random perturbations. It is shown that the ESV-based predictions have a systematically higher bivariate correlation skill in predicting the MJO compared to those using the random perturbations. Furthermore, the improvement in the skill depends on the phase of the MJO. The ESV is particularly effective in increasing the forecast skill during those phases of the MJO in which the control has low skill (with correlations increasing by as much as 0.2 at 20 25-day lead times), as well as during those times in which the MJO is weak.

Ham, Yoo-Geun↗
NASA NTRS2013

The Global Land-Ocean Temperature Index in Relation to Sunspot Number, the Atlantic Multidecadal Oscillation Index, the Mauna Loa Atmospheric Concentration of CO2, and Anthropogenic Carbon Emissions

Global warming/climate change has been a subject of scientific interest since the early 19th century. In particular, increases in the atmospheric concentration of carbon dioxide (CO2) have long been thought to account for Earth's increased warming, although the lack of a dependable set of observational data was apparent as late as the mid 1950s. However, beginning in the late 1950s, being associated with the International Geophysical Year, the opportunity arose to begin accurate continuous monitoring of the Earth's atmospheric concentration of CO2. Consequently, it is now well established that the atmospheric concentration of CO2, while varying seasonally within any particular year, has steadily increased over time. Associated with this rising trend in the atmospheric concentration of CO2 is a rising trend in the surface-air and sea-surface temperatures (SSTs). This Technical Publication (TP) examines the statistical relationships between 10-year moving averages (10-yma) of the Global Land-Ocean Temperature Index (GLOTI), sunspot number (SSN), the Atlantic Multidecadal Oscillation (AMO) index, and the Mauna Loa CO2 (MLCO2) index for the common interval 1964-2006, where the 10-yma values are used to indicate trends in the data. Scatter plots using the 10-yma values between GLOTI and each of the other parameters are determined, both as single-variate and multivariate fits. Scatter plots are also determined for MLCO2 using single-variate and bivariate (BV) fits, based on the GLOTI alone and the GLOTI in combination with the AMO index. On the basis of the inferred preferential fits for MLCO2, estimates for MLCO2 are determined for the interval 1885-1964, thereby yielding an estimate of the preindustrial level of atmospheric concentration of CO2. Lastly, 10-yma values of MLCO2 are compared against 10-yma estimates of the total carbon emissions (TCE) to determine the likelihood that manmade sources of carbon emissions are indeed responsible for the recent warming now being experienced. (Parametric values used in this TP are those available prior to the end of 2012.)

Wilson, Robert M.↗
NASA NTRS2012

Tree Density and Species Decline in the African Sahel Attributable to Climate

Increased aridity and human population have reduced tree cover in parts of the African Sahel and degraded resources for local people. Yet, tree cover trends and the relative importance of climate and population remain unresolved. From field measurements, aerial photos, and Ikonos satellite images, we detected significant 1954-2002 tree density declines in the western Sahel of 18 +/- 14% (P = 0.014, n = 204) and 17 +/- 13% (P = 0.0009, n = 187). From field observations, we detected a significant 1960-2000 species richness decline of 21 +/- 11% (P = 0.0028, n = 14) across the Sahel and a southward shift of the Sahel, Sudan, and Guinea zones. Multivariate analyses of climate, soil, and population showed that temperature most significantly (P < 0.001) explained tree cover changes. Multivariate and bivariate tests and field observations indicated the dominance of temperature and precipitation, supporting attribution of tree cover changes to climate variability. Climate change forcing of Sahel climate variability, particularly the significant (P < 0.05) 1901-2002 temperature increases and precipitation decreases in the research areas, connects Sahel tree cover changes to global climate change. This suggests roles for global action and local adaptation to address ecological change in the Sahel.

tree density declines↗
NASA NTRS2014

Recent Developments and Challenges Implementing New and Improved Stress Intensity Factor (K) Solutions in NASGRO for Damage Tolerance Analyses

Fatigue crack growth analysis software has been available to damage tolerance analysts for many years in either commercial products or via proprietary in-house codes. The NASGRO software has been publicly available since the mid-80s (known as NASA/FLAGRO up to 1999) and since 2000 has been sustained and further developed by a collaborative effort between Southwest Research Institute® (SwRI®), the NASA Johnson Space Center (JSC), and the members of the NASGRO Industrial Consortium. Since the stress intensity factor (K) is the foundation of fracture mechanics and damage tolerance analysis of aircraft structures, a significant focus of development efforts in the past fifteen years has been geared towards enhancing legacy K solutions and developing new and efficient numerical K solutions that can handle the complicated stress gradients computed by today’s analysts using detailed finite element models of fatigue critical locations. This paper provides an overview of K solutions that have been recently implemented or improved for the analysis of geometries such as two unequal through cracks at a hole and two unequal corner cracks at a hole, as well as state-of-the-art weight function models capable of computing K in the presence of univariant and/or bivariant stress gradients and complicated residual stress distributions. Some historical background is provided to review how common K solutions have evolved over the years, including selective examples from the literature and from new research. Challenges and progress in rectifying discrepancies between older legacy solutions and newer models are reviewed as well as approaches and challenges for verification and validation of K solutions. Finally, a summary of current challenges and future research and development needs is presented. A key theme throughout the presentation of this paper will be how members of the aerospace industry have collaborated with software developers to develop a practical analysis tool that is used world-wide to support new design as well as the ongoing sustainment and airworthiness of commercial and military aircraft.

Cardinal, Joseph W.↗
NASA NTRS2015

Predicting the Redshift 2 H-Alpha Luminosity Function Using [OIII] Emission Line Galaxies

Upcoming space-based surveys such as Euclid and WFIRST-AFTA plan to measure Baryonic Acoustic Oscillations (BAOs) in order to study dark energy. These surveys will use IR slitless grism spectroscopy to measure redshifts of a large number of galaxies over a significant redshift range. In this paper, we use the WFC3 Infrared Spectroscopic Parallel Survey (WISP) to estimate the expected number of H-alpha emitters observable by these future surveys. WISP is an ongoing Hubble Space Telescope slitless spectroscopic survey, covering the 0.8 - 1.65 micrometers wavelength range and allowing the detection of H-alpha emitters up to z approximately equal to 1.5 and [OIII] emitters to z approximately equal to 2.3. We derive the H-alpha-[OIII] bivariate line luminosity function for WISP galaxies at z approximately equal to 1 using a maximum likelihood estimator that properly accounts for uncertainties in line luminosity measurement, and demonstrate how it can be used to derive the H-alpha luminosity function from exclusively fitting [OIII] data. Using the z approximately equal to 2 [OIII] line luminosity function, and assuming that the relation between H-alpha and [OIII] luminosity does not change significantly over the redshift range, we predict the H-alpha number counts at z approximately equal to 2 - the upper end of the redshift range of interest for the future surveys. For the redshift range 0.7 less than z less than 2, we expect approximately 3000 galaxies per sq deg for a flux limit of 3 x 10(exp −16) ergs per sec per sq cm (the proposed depth of Euclid galaxy redshift survey) and approximately 20,000 galaxies per sq deg for a flux limit of approximately 10(exp −16) ergs per sec per sq cm (the baseline depth of WFIRST galaxy redshift survey).

Redshift↗
NASA NTRS2015

Continental Spatio-Temporal Data Analysis with Linear Spectral Mixture Model Using FOSS

This work demonstrates the development and implementation of a Fully Constrained Least Squares (FCLS) unmixing model developed in C++ programming language with OpenCV package and boost C++ libraries in the NASA Earth Exchange (NEX). Visualization of the results is supported by GRASS GIS and statistical analysis is carried in R in a Linux system environment. FCLS was first tested on computer simulated data with Gaussian noise of various signal-to-noise ratio, and Landsat data of an agricultural scenario and an urban environment using a set of global end members of substrate (soils, sediments, rocks, and non-photosynthetic vegetation), vegetation that includes green photosynthetic plants and dark objects which encompasses absorptive substrate materials, clear water, deep shadows, etc. For the agricultural scenario, a spectrally diverse collection of 11 scenes of Level 1 terrain corrected, cloud free Landsat-5 TM data of Fresno, California, USA were unmixed and the results were validated with the corresponding ground data. To study an urbanized landscape, a clear sky Landsat-5 TM data were unmixed and validated with coincident World View-2 abundance maps (of 2 m spatial resolution) for an area of San Francisco, California, USA. The results were evaluated using descriptive statistics, correlation coefficient, RMSE, probability of success, boxplot and bivariate distribution function. Finally, FCLS was used for sub-pixel land cover analysis of the monthly WELD (Wen-enabled Landsat data) repository from 2008 to 2011 of North America. The abundance maps in conjunction with DMSP-OLS nighttime lights data were used to extract the urban land cover features and analyze their spatial-temporal growth.

Landsat Satellites↗
NASA NTRS2017

New Mathematical Functions for Vacuum System Analysis

A new bivariate function has been found that provides solutions of integrals having the form u (sup minus eta) e (sup u) du which arise when developing predictions for the behavior of pressure within a rigid volume under high vacuum conditions in the presence of venting as well as sources characterized by power law transient decay over the range [0,1] for eta and for u greater than or equal to 0. A few properties of the new function are explored in this work. For instance the eta equals 1/2 case reproduces the Dawson function. In addition, a slight variation of the solution technique reproduces the exponential integral for eta equals 1. The technique used to generate these functions leads to an approach for solving a more general class of nonlinear ordinary differential equations, with the potential for identifying other new functions that solve other integrals.

Dawson Function↗
NASA NTRS2017

The Impact of Anthropogenic Land Use and Land Cover Change on Regional Climate Extremes

Recent research highlights the role of land surface processes in heat waves, droughts, and other extreme events. Here we use an earth system model (ESM) from the Geophysical Fluid Dynamics Laboratory (GFDL) to investigate the regional impacts of historical anthropogenic land useland cover change (LULCC) on combined extremes of temperature and humidity. A bivariate assessment allows us to consider aridity and moist enthalpy extremes, quantities central to human experience of near-surface climate conditions. We show that according to this model, conversion of forests to cropland has contributed to much of the upper central US and central Europe experiencing extreme hot, dry summers every 2-3 years instead of every 10 years. In the tropics, historical patterns of wood harvesting, shifting cultivation and regrowth of secondary vegetation have enhanced near surface moist enthalpy, leading to extensive increases in the occurrence of humid conditions throughout the tropics year round. These critical land use processes and practices are not included in many current generation land models, yet these results identify them as critical factors in the energy and water cycles of the midlatitudes and tropics.

Near-surface↗
NASA NTRS2018

Prediction and Predictability of the Madden Julian Oscillation in the NASA GEOS-5 Seasonal-to-Subseasonal System

In this study, we examine the prediction skill and predictability of the Madden Julian Oscillation (MJO) in a recent version of the NASA GEOS-5 atmosphere-ocean coupled model run at at 1/2 degree horizontal resolution. The results are based on a suite of hindcasts produced as part of the NOAA SubX project, consisting of seven ensemble members initialized every 5 days for the period 1999-2015. The atmospheric initial conditions were taken from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and the ocean and the sea ice were taken from a GMAO ocean analysis. The land states were initialized from the MERRA-2 land output, which is based on observation-corrected precipitation fields. We investigated the MJO prediction skill in terms of the bivariate correlation coefficient for the real-time multivariate MJO (RMM) indices. The correlation coefficient stays at or above 0.5 out to forecast lead times of 26-36 days, with a pronounced increase in skill for forecasts initialized from phase 3, when the MJO convective anomaly is located in the central tropical Indian Ocean. A corresponding estimate of the upper limit of the predictability is calculated by considering a single ensemble member as the truth and verifying the ensemble mean of the remaining members against that. The predictability estimates fall between 35-37 days (taken as forecast lead when the correlation reaches 0.5) and are rather insensitive to the initial MJO phase. The model shows slightly higher skill when the initial conditions contain strong MJO events compared to weak events, although the difference in skill is evident only from lead 1 to 20. Similar to other models, the RMM-index-based skill arises mostly from the circulation components of the index. The skill of the convective component of the index drops to 0.5 by day 20 as opposed to day 30 for circulation fields. The propagation of the MJO anomalies over the Maritime Continent does not appear problematic in the GEOS-5 hindcasts implying that the Maritime Continent predictability barrier may not be a major concern in this model. Finally, the MJO prediction skill in this version of GEOS-5 is superior to that of the current seasonal prediction system at the GMAO; this could be partly attributed to a slightly better representation of the MJO in the free running version of this model and partly to the improved atmospheric initialization from MERRA-2.

Achuthavarier, Deepthi↗
NASA NTRSDate not supplied

Utilizing Open-Source Earth Observations to Inform the Toa Baja Municipality’s Flood Risk Mitigation Efforts and Educate the Public

Global climate changes contribute to more intense and frequent tropical storms, subjecting places like Toa Baja, Puerto Rico to critical damage. Known as “the underwater city” due to its propensity to flood, residents of Toa Baja face constant flood risk. During extreme tropical storm events, such as Hurricane Maria in 2017, residents experienced up to 20 feet of inundation. The NASA DEVELOP National Program collaborated with the Municipio Autónomo de Toa Baja, ResilientSEE-PR, and the MIT Urban Risk Lab to supplement 2018 FEMA HEC-RAS flood maps that designate 63% of Toa Baja as a flood plain. This analysis provides a high-resolution interpretation of flood risk through two lenses; susceptibility and vulnerability. For this analysis, susceptibility consists of nine weighted layers: NDVI, landcover, slope, elevation, topographic wetness index, height above nearest drainage, saturated hydraulic conductivity, distance to water, and storm surge. These factors are consistently used to evaluate susceptibility to flood, but their weights vary by analysis. Vulnerability consists of population, informal settlements, and building density, which were given equal weight. Susceptibility and vulnerability were combined to map flood risk. This analysis used a bivariate legend to understand the different levels of risk along a spectrum from low susceptibility and low vulnerability (low risk) to high susceptibility and high vulnerability (high risk). Data processed in Google Earth Engine, which identified historical inundation on various occasions, were used to validate the flood susceptibility layers. Results showed 89% of areas designated as high susceptibility are located within the floodway designated by the FEMA HEC-RAS maps. The eastern region of Toa Baja is most at risk for flooding due to high susceptibility to flooding along with a high density of population, buildings, and informal settlements. The resulting map also reveals the presence of smaller high-risk areas all around the municipality. This analysis provides scientific evidence for flood risk mitigation in Toa Baja by highlighting areas that might be impacted by strong floods in the future. Additionally, these results are communicated in an Esri ArcGIS StoryMap, an accessible platform that can easily inform the public about the flood risk in their neighborhood.

Adriana Le Compte↗
NASA NTRS2022

Understanding the Microphysical Control and Spatial-Temporal Variability of Warm Rain Probability Using CloudSat and MODIS Observations

By combining measurements from MODIS and the CloudSat radar, we develop a parameterization scheme to quantify the combined microphysical controls by liquid water path (LWP) and cloud droplet number concentration (CDNC) of the probability of precipitation (PoP) in marine low cloud over tropical oceans. We demonstrate that the spatial-temporal variation of grid-mean in-cloud can be largely explained by the variation of the joint probability density function of LWP and CDNC in the phase space specified by the bivariate PoP (LWP and CDNC) function. Through a series of sensitivity tests guided by this understanding, we find that in the Southeastern Pacific and Atlantic the stratocumulus to cumulus transition of the is mainly due to the variation of CDNC while the annual cycle is mainly due to the variation of LWP. The results of this study provide a viable way to diagnose the root cause of warm rain problems in global climate models.

Zhibo Zhang↗
NASA NTRS2023

Great Salt Lake Health and Air Quality: Monitoring Lakebed Exposure and its Impact on Air Quality and Environmental Hazards in the Great Salt Lake Watershed

Water flow into the Great Salt Lake has declined rapidly over the last forty years due to human withdrawals and climate change. As a result of declining lake levels, over 50% of the lakebed is now exposed. Dust storms may grow in frequency and intensity across Northern Utah as lakebed dust becomes airborne under specific meteorological conditions. In our research project, we utilized satellite imagery from Terra and Aqua, Sentinel-5P, CALIPSO, Landsat 5 TM, Landsat 7 ETM+, Landsat 8 OLI-2, Suomi NPP, ground sensor environmental data, and demographic data to understand the relationship between lake desiccation and dust, and the impact of pollution upon the communities surrounding the Great Salt Lake. By plotting changes in Lake Surface Area against Aerosol Optical Depth (AOD) over our study period (2010-2022), we found an inverse relationship (R2=0.3423) between lake surface area and dust levels within our study area. We conducted a Vertical Feature Mask (VFM) and Extinction Coefficient Plot, from which we identified that during dust events, the aerosol type is mainly polluted dust and the aerosol height is 200 meters from the surface. Lastly, we created bivariate choropleth maps, which demonstrate which census tracts within our study area are most vulnerable to AOD (a proxy for PM2.5 from dust), NO2 and HCHO (precursors to ozone). In summary, our findings revealed that declining lake levels are associated with an increase in intensity of dust events, and these dust events will particularly impact residents of Tooele County and the west side of Salt Lake City. Project resources support partner needs by informing targeted air monitoring efforts, lakebed management practices, and advocacy efforts for GSL stewardship.

Terminal Saline Lake↗
NASA NTRS2024

Oculometric Biomarkers of Visuomotor Deficits in Clinically Asymptomatic Patients With Systemic Lupus Erythematosus Undergoing Long-Term Hydroxychloroquine Treatment

Introduction: This study examines a set of oculomotor measurements, or “oculometric” biomarkers, as potential early indicators of visual and visuomotor deficits due to retinal toxicity in asymptomatic Systemic Lupus Erythematosus (SLE) patients on long-term hydroxychloroquine (HCQ) treatment. The aim is to identify subclinical functional impairments that are otherwise undetectable by standard clinical tests and to link them to structural retinal changes. Methods: We measured oculomotor responses in a cohort of SLE patients on chronic HCQ therapy using a previously established behavioral task and analysis technique. We also examined the relationship between oculometrics, OCT measures of retinal thickness, and standard clinical perimetry measures of visual function in our patient group using Bivariate Pearson Correlation and a Linear Mixed-Effects Model (LMM). Results: Significant visual and visuomotor deficits were found in 12 asymptomatic SLE patients on long-term HCQ therapy compared to a cohort of 17 age-matched healthy controls. Notably, six oculometrics were significantly different. The median initial pursuit acceleration was 22%, steady-state pursuit gain 16%, proportion smooth 7%, and target speed responsiveness 31% lower, while catch-up saccade amplitude was 46% and fixation error 46% larger. Excluding the two patients with diagnosed mild toxicity, four oculometrics, all but fixation error and proportion smooth, remained significantly impaired compared to controls. Across our population of 12 patients (24 retinae), we found that pursuit latency, initial acceleration, steady-state gain, and fixation error were linearly related to retinal thickness even when age was accounted for, while standard measures of clinical function (Mean Deviation and Pattern Standard Deviation) were not. Discussion: Our data show that specific oculometrics are sensitive early biomarkers of functional deficits in SLE patients on HCQ that could be harnessed to assist in the early detection of HCQ-induced retinal toxicity and other visual pathologies, potentially providing early diagnostic value beyond standard visual field and OCT evaluations.

eye movements↗
DOE OSTI2024

Data analytics for intermodal freight transportation applications

With the growth of intermodal freight transportation, it is important that transportation planners and decision-makers are knowledgeable about freight flow data to make informed decisions. This is particularly true with Intelligent Transportation Systems (ITS) offering new capabilities for intermodal freight transportation. Specifically, ITS enables access to multiple different data sources, but they have different formats, resolutions, and time scales. Thus, knowledge of data science is essential to be successful in future ITS-enabled intermodal freight transportation systems. This chapter discusses the commonly used descriptive and predictive data analytic techniques in intermodal freight transportation applications. These techniques cover the entire spectrum of univariate, bivariate, and multivariate analyses. In addition to illustrating how to apply these techniques manually, this chapter will also show how to apply them using the statistical software R. Additional exercises are provided for those who wish to apply the described techniques to more complex problems.

Huynh, Nathan↗
DOE OSTI2026

Efficient Routing of Quantum LDPC Codes on Programmable 2D Toric Architectures

Quantum low-density parity-check codes are promising candidates towards scalable fault-tolerant quantum computation. Among these, bivariate bicycle (BB) codes offer superior encoding rates and large code distance compared to surface codes. However, their requirement on long-range stabilizer measurements poses significant challenges for implementation on realistic hardware with limited connectivity, such as superconducting circuit platforms. In this work, we introduce a novel hardware-software co-design that leverages a programmable communication network architecture to address these limitations. Our approach utilizes a 2D toric network of oscillators as a flexible communication fabric linking qubits at each site. Such architecture significantly reduces the number of long-range couplers required from O ( n ) to O (√ n ). Dual-rail qubits, along with native gates including Swap-Wait-Swap gates and beamsplitter SWAPs, ensure that long-range two-qubit gates can be executed with high fidelity and low latency. To further enhance performance, our qubit layout and routing algorithm utilize symmetries of the codes and enable maximum parallelism for long-range two-qubit gates, maintaining a low syndrome extraction cycle duration and scalability over the code length. We perform circuit-level simulation with realistic noise modeling based on experimental hardware parameters, observing an logical error rate per logical qubit per cycle of 3.06% for [[18,4,4]] BB code, 2.6× less than the existing experimental result. These findings provide a practical roadmap and identify key technological advancements needed to achieve low-overhead fault-tolerant quantum computing at scale.

Liu, Kun [Yale Univ., New Haven, CT (United States↗
DOE OSTI2025

Multitaper Magnitude‐Squared Coherence for Time Series With Missing Data: Understanding Oscillatory Processes Traced by Multiple Observables

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic data sets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman α flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

Astro-statistics techniques (1886)↗