Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial statistics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71

Seasonal Variability of Middle Latitude Ozone in the Lowermost Stratosphere Derived from Probability Distribution Functions

We present a study of the distribution of ozone in the lowermost stratosphere with the goal of understanding the relative contribution to the observations of air of either distinctly tropospheric or stratospheric origin. The air in the lowermost stratosphere is divided into two population groups based on Ertel's potential vorticity at 300 hPa. High [low] potential vorticity at 300 hPa suggests that the tropopause is low [high], and the identification of the two groups helps to account for dynamic variability. Conditional probability distribution functions are used to define the statistics of the mix from both observations and model simulations. Two data sources are chosen. First, several years of ozonesonde observations are used to exploit the high vertical resolution. Second, observations made by the Halogen Occultation Experiment [HALOE) on the Upper Atmosphere Research Satellite [UARS] are used to understand the impact on the results of the spatial limitations of the ozonesonde network. The conditional probability distribution functions are calculated at a series of potential temperature surfaces spanning the domain from the midlatitude tropopause to surfaces higher than the mean tropical tropopause [approximately 380K]. Despite the differences in spatial and temporal sampling, the probability distribution functions are similar for the two data sources. Comparisons with the model demonstrate that the model maintains a mix of air in the lowermost stratosphere similar to the observations. The model also simulates a realistic annual cycle. By using the model, possible mechanisms for the maintenance of mix of air in the lowermost stratosphere are revealed. The relevance of the results to the assessment of the environmental impact of aircraft effluence is discussed.

Cerniglia, M. C.↗

Visualizing Distributions from Multi-Return Lidar Data to Understand Forest Structure

Spatially distributed probability density functions (pdfs) are becoming relevant to the Earth scientists and ecologists because of stochastic models and new sensors that provide numerous realizations or data points per unit area. One source of these data is from multi-return airborne lidar, a type of laser that records multiple returns for each pulse of light sent towards the ground. Data from multi-return lidar is a vital tool in helping us understand the structure of forest canopies over large extents. This paper presents several new visualization tools that allow scientists to rapidly explore, interpret and discover characteristic distributions within the entire spatial field. The major contribution from-this work is a paradigm shift which allows ecologists to think of and analyze their data in terms of the distribution. This provides a way to reveal information on the modality and shape of the distribution previously not possible. The tools allow the scientists to depart from traditional parametric statistical analyses and to associate multimodal distribution characteristics to forest structures. Examples are given using data from High Island, southeast Alaska.

Kao, David L.↗

The Galex Time Domain Survey. I. Selection And Classification of Over a Thousand Ultraviolet Variable Sources

We present the selection and classification of over a thousand ultraviolet (UV) variable sources discovered in approximately 40 deg(exp 2) of GALEX Time Domain Survey (TDS) NUV images observed with a cadence of 2 days and a baseline of observations of approximately 3 years. The GALEX TDS fields were designed to be in spatial and temporal coordination with the Pan-STARRS1 Medium Deep Survey, which provides deep optical imaging and simultaneous optical transient detections via image differencing. We characterize the GALEX photometric errors empirically as a function of mean magnitude, and select sources that vary at the 5 sigma level in at least one epoch. We measure the statistical properties of the UV variability, including the structure function on timescales of days and years. We report classifications for the GALEX TDS sample using a combination of optical host colors and morphology, UV light curve characteristics, and matches to archival X-ray, and spectroscopy catalogs. We classify 62% of the sources as active galaxies (358 quasars and 305 active galactic nuclei), and 10% as variable stars (including 37 RR Lyrae, 53 M dwarf flare stars, and 2 cataclysmic variables). We detect a large-amplitude tail in the UV variability distribution for M-dwarf flare stars and RR Lyrae, reaching up to absolute value(Δm) = 4.6 mag and 2.9 mag, respectively. The mean amplitude of the structure function for quasars on year timescales is five times larger than observed at optical wavelengths. The remaining unclassified sources include UV-bright extragalactic transients, two of which have been spectroscopically confirmed to be a young core-collapse supernova and a flare from the tidal disruption of a star by dormant supermassive black hole. We calculate a surface density for variable sources in the UV with NUV less than 23 mag and absolute value(Δm) greater than 0.2 mag of approximately 8.0, 7.7, and 1.8 deg(exp −2) for quasars, active galactic nuclei, and RR Lyrae stars, respectively. We also calculate a surface density rate in the UV for transient sources, using the effective survey time at the cadence appropriate to each class, of approximately 15 and 52 deg(exp −2 yr−1 for M dwarfs and extragalactic transients, respectively.

Gezari, S.↗

Visual Analytics of Crosstalk in Quantum Hardware

Crosstalk remains a major obstacle to building scalable and fault-tolerant quantum computers. Conventional diagnostic techniques-often based on numerical simulation or statistical modeling-struggle to scale with hardware complexity and offer limited interpretability. In this work, we present a visual analytics framework for diagnosing qubit crosstalk using lightweight, circuit-based models integrated with an interactive user interface. Our approach quantifies correlations between active and idle qubits under parameterized single- and twoqubit operations, enabling detection of both spatial and gateinduced crosstalk. The system incorporates qubit topology and gate performance data to support sector-based exploration and correlation mapping. This tool assists users in identifying correlated error sources, informing qubit placement strategies, and guiding noise-aware circuit design.

Chae, Junghoon [ORNL] (ORCID:0000000206016746)↗

Dynamically Downscaled Hourly Future Weather Data with 12-km Resolution Covering Most of North America

This is an hourly future weather dataset for energy modeling applications. The dataset is primarily based on the output of a regional climate model (RCM), i.e., the Weather Research and Forecasting (WRF) model version 3.3.1. The WRF simulations are driven by the output of a general circulation model (GCM), i.e., the Community Climate System Model version 4 (CCSM4). This dataset is in the EPW format, which can be read or translated by more than 25 building energy modeling programs (e.g., EnergyPlus, ESP-r, and IESVE), energy system modeling programs (e.g., System Advisor Model (SAM)), indoor air quality analysis programs (e.g., CONTAM), and hygrothermal analysis programs (e.g., WUFI). It contains 13 weather variables, which are the Dry-Bulb Temperature, Dew Point Temperature, Relative Humidity, Atmospheric Pressure, Horizontal Infrared Radiation Intensity from Sky, Global Horizontal Irradiation, Direct Normal Irradiation, Diffuse Horizontal Irradiation, Wind Speed, Wind Direction, Sky Cover, Albedo, and Liquid Precipitation Depth. The weather data is created for two emissions scenarios: RCP4.5 and RCP8.5 and spans two 10-year time slices in the future: 2045 - 2054 and 2085 - 2094. It offers a spatial resolution of 12 km by 12 km with extensive coverage across most of North America. Due to the enormous size of the entire dataset, in the first stage of its distribution, we provide 20 years of future weather data for the centroid of each Public Use Microdata Area (PUMA), excluding Hawaii. PUMAs are non-overlapping, statistical geographic areas that partition each state or equivalent entity into geographic areas containing no fewer than 100,000 people each. The 2,378 PUMAs as a whole cover the entirety of the U.S. The weather data can be utilized alongside the large-scale energy analysis tools, ResStock and ComStock, developed by National Renewable Energy Laboratory, whose smallest resolution is at the PUMA scale. The data for RCP4.5 is still being processed and will be published soon.

Array↗

Engineering studies related to the GEOS-C radar altimeter

Experiment requirements, technical characteristics, and GEOS-C radar altimeter related analyses are discussed along with results of a study on engineering test data requirements. Statistical analyses related to determination of wave height resolution achievable as a function of system characteristics and averaging period are described, in addition to a discussion on the desirability of using computer procedures to compensate for altitude tracker time-jitter. Data processing considerations for the GEOS-C system are examined. An extensive analysis of the spatial filter effect is given and results of a computation of geoidal power spectral density, based on Skylab altimeter data, is displayed and interpreted in terms of projected GEOS-C random errors. This information is then used in deriving minimum-mean-square filter procedures for both geoid undulation and slope data. The characteristics of mean received waveforms as a function of off-nadir angle are used to obtain tracker bias as a function of sea state and pointing angle. The angle estimation process proposed by the GEOS-C hardware contractor is also investigated from a standpoint of achievable angular resolution.

Miller, L. S.↗

Computer analysis of X-band radar data

The effectiveness of using currently available computer techniques for interpretation of MSS data to interpret SAR imagery for forest monitoring was assessed. Data were gathered with NASA's airborne APQ-102 dual-polarized, X-band SAR in a flight at 60,000 ft. Microdensitometry was employed to digitize the HH- and HV-polarized imagery. A ground spatial resolution of 15 m was obtained, control points were identified, a second order biquadratic transformation was applied to compensate for orientation, and rms errors were calculated. A second data set was taken with 30 m resolution in order to simulate thematic mapper operation. Classification was performed with pixel-by-pixel and textural classification algorithms. A statistical analysis was also carried out to find any significant differences between classifiers in a data set for a given classifier. Each polarization featured an independent distortion which required appropriate preprocessing to correct. Further studies are recommended with multiple frequencies viewing and multiple polarizations and look angles to define the actual forest classifications that can be made with the SAR imagery.

Knowlton, D. J.↗

Stratified Magnetically Driven Accretion-Disk Winds and Their Relations To Jets

We explore the poloidal structure of two-dimensional magnetohydrodynamic (MHD) winds in relation to their potential association with the X-ray warm absorbers (WAs) and the highly ionized ultra-fast outflows (UFOs) in active galactic nuclei (AGNs), in a single unifying approach. We present the density n(r, theta), ionization parameter xi(r, theta), and velocity structure v(r, theta) of such ionized winds for typical values of their fluid-to-magnetic flux ratio, F, and specific angular momentum, H, for which wind solutions become super-Alfvenic. We explore the geometrical shape of winds for different values of these parameters and delineate the values that produce the widest and narrowest opening angles of these winds, quantities necessary in the determination of the statistics of AGN obscuration. We find that winds with smaller H show a poloidal geometry of narrower opening angles with their Alfv´en surface at lower inclination angles and therefore they produce the highest line of sight (LoS) velocities for observers at higher latitudes with the respect to the disk plane. We further note a physical and spatial correlation between the X-ray WAs and UFOs that form along the same LoS to the observer but at different radii, r, and distinct values of n, xi, and v consistent with the latest spectroscopic data of radio-quiet Seyfert galaxies. We also show that, at least in the case of 3C 111, the winds' pressure is sufficient to contain the relativistic plasma responsible for its radio emission. Stratified MHD disk winds could therefore serve as a unique means to understand and unify the diverse AGN outflows.

accretion↗

Observations of Turbulence in a Kelvin-Helmholtz Event on 8 September 2015 by the Magnetospheric Multiscale Mission

Spatial and high-time-resolution properties of the velocities, magnetic field, and 3-D electric field within plasma turbulence are examined observationally using data from the Magnetospheric Multiscale mission. Observations from a Kelvin-Helmholtz instability (KHI) on the Earth's magnetopause are examined, which both provides a series of repeatable intervals to analyze, giving better statistics, and provides a first look at the properties of turbulence in the KHI. For the first time direct observations of both the high-frequency ion and electron velocity spectra are examined, showing differing ion and electron behavior at kinetic scales. Temporal spectra exhibit power law behavior with changes in slope near the ion gyrofrequency and lower hybrid frequency. The work provides the first observational evidence for turbulent intermittency and anisotropy consistent with quasi two-dimensional turbulence in association with the KHI. The behavior of kinetic-scale intermittency is found to have differences from previous studies of solar wind turbulence, leading to novel insights on the turbulent dynamics in the KHI.

Stawarz, J. E.↗

A Bayesian Analysis of SDSS J0914+0853, a Low-mass Dual AGN Candidate

We present the first results from Bayesian AnalYsis of Multiple AGN in X-rays (BAYMAX), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual active galactic nuclei (AGNs) is to use X-ray observations, only sources that are widely separated relative to the instrumentʼs point-spread function are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandraʼs angular resolution (<1″). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ∼ 0"3 from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find a very low probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0"3. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, as well as studying the dual AGN population where past spatial resolution limits have prevented systematic analyses.

Active galaxies↗

Automated Tracking of Shallow Maritime Clouds on Geostationary Imagery to Extract Lifecycle Characteristics

Shallow moist convection is ubiquitous throughout the tropics and represents a key player in boundary layer processes. Satellites have provided many statistics on shallow clouds, such as size, structure, and geographical coverage, from static views of recurring cloud fields. But determining why certain cloud features appear and persist for different periods requires a time-evolving view of their behaviors. Geostationary satellites provide a unique opportunity to follow the time evolution of individual convective features, given their enhanced spatial and temporal sampling. A cloud-tracking tool was developed to identify properties of cloud lifecycle from the NASA Cloud, Aerosol, and Monsoon Processes Philippines Experiment (CAMP2EX) field campaign of 2019. The mission conducted intensive sampling of shallow cumulus in the West Pacific Ocean, in tandem with Rapid Scan imagery from the Advanced Himawari Imager (AHI) on the Japan Meteorological Agency’s (JMA) Himawari-8 satellite. Shallow cumulus was segmented according to thresholds in 0.5-km visible reflectance and with blurring techniques. Despite being limited to daytime hours, the segmentations yielded the best resolution possible for capturing cloud initiation and decay. The tracking procedure is based on a computer vision package that includes Kalman filters for motion prediction, object overlap search, and the Hungarian (or Kuhn-Munkres) matching algorithm for track designation. AHI radiances available within the tracked cloud boundaries are assembled to form individual spectral histories. The resulting catalog provides thousands of cloud histories for domains measuring only a few degrees in latitude and longitude. We present an overview of the cloud-tracking tool, strategies to identify development stages from cloud tracks, and preliminary results that document cumulus lifecycle properties from satellite. The application of AHI 0.5-km reflectance has both strengths and limitations when attempting to track lifecycles of the smallest resolvable clouds. We show that by aggregating cloud tracks from a few case studies of CAMP2EX, we can discern differences in cloud lifetime and development according to ensembles selected from areas of interest. The results demonstrate an ability to quantify lifetimes and assess rates of change in cloud characteristics that are likely controlled by the surrounding environment and meteorology.

Cloud Tracking↗

Comparing Satellite Measurements of Volcanic SO2 Mass from OMI, OMPS and TROPOMI

Sulfur dioxide (SO2) is a major air pollutant that contributes to acid rain and aerosol formation (e.g., sulfates), adversely affects the environment and human health, and explosive volcanic SO2 emissions can impact climate. The majority of SO2 emissions are related to anthropogenic processes (e.g., fossil fuel burning, metal ore smelting operations), although natural processes such as volcanic eruptions and degassing also play an important role as anthropogenic SO2 emissions continue to decline. Generally, the most interest in volcanoes occurs during major eruptions. We will focus on comparing volcanic SO2 outgassing that occurs on an almost daily basis from lesser known volcanoes using satellite data. At NASA’s Global Sulfur Dioxide Monitoring Home page (https://so2.gsfc.nasa.gov/), we have been posting daily SO2 maps from 40 volcanic and industrial regions around the world using measurements from three satellite instruments; the Ozone Monitoring Instrument (OMI) onboard NASA’s Earth Observing System Aura satellite, the Ozone Monitoring and Profiler Suite (OMPS) onboard the NASA-NOAA Suomi National Polar-orbiting Partnership (NPP) satellite, and the TROPOspheric Monitoring Instrument (TROPOMI) onboard the ESA/Copernicus Sentinel-5 Precursor satellite. These instruments in low Earth sun-synchronous polar orbits with 1:30-2:00 pm equator crossing local time provide daily SO2 maps at different spatial resolutions: 13 x 24 km2, 50 x 50 km2 and 5.5 x 3.5 km2 for OMI, OMPS and TROPOMI respectively. Data from OMI are available since October 2004 (partial coverage since 2008), from OMPS since 2012 and from TROPOMI since 2018. We will present comparative SO2 mass time-series (see Hunga-Tonga plot) and statistical analyses of recent eruptions that have data from all the instruments.

SO2↗

Human-caused long-term changes in global aridity

Widespread aridification of the land surface causes substantial environmental challenges and is generally well documented. However, the mechanisms underlying increased aridity remain relatively underexplored. Here, we investigated the anthropogenic and natural factors affecting long-term global aridity changes using multisource observation-based aridity index, factorial simulations from the Coupled Model Intercomparison Project phase 6 (CMIP6), and rigorous detection and attribution (D&A) methods. Our study found that anthropogenic forcings, mainly rising greenhouse gas emissions (GHGE) and aerosols, caused the increased aridification of the globe and each hemisphere with high statistical confidence for 1965–2014; the GHGE contributed to drying trends, whereas the aerosol emissions led to wetting tendencies; moreover, the bias-corrected CMIP6 future aridity index based on the scaling factors from optimal D&A demonstrated greater aridification than the original simulations. These findings highlight the dominant role of human effects on increasing aridification at broad spatial scales, implying future reductions in aridity will rely primarily on the GHGE mitigation.

54 ENVIRONMENTAL SCIENCES↗

Boundary-layer receptivity due to distributed surface imperfections of a deterministic or random nature

Acoustic receptivity of a Blasius boundary layer in the presence of distributed surface irregularities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface irregularities, only a small band of Fourier components can lead to an efficient conversion of the acoustic input at any given frequency to an unstable eigenmode of the boundary layer flow. The location, and width, of this most receptive band of wavenumbers corresponds to a relative detuning of O(R sub l.b.(exp -3/8)) with respect to the lower-neutral instability wavenumber at the frequency under consideration, R sub l.b. being the Reynolds number based on a typical boundary-layer thickness at the lower branch of the neutral stability curve. Surface imperfections in the form of discrete mode waviness in this range of wavenumbers lead to initial instability amplitudes which are O(R sub l.b.(exp 3/8)) larger than those caused by a single, isolated roughness element. In contrast, irregularities with a continuous spatial spectrum produce much smaller instability amplitudes, even compared to the isolated case, since the increase due to the resonant nature of the response is more than that compensated for by the asymptotically small band-width of the receptivity process. Analytical expressions for the maximum possible instability amplitudes, as well as their expectation for an ensemble of statistically irregular surfaces with random phase distributions, are also presented.

Choudhari, Meelan↗

An Extended Objective Evaluation of the 29-km Eta Model for Weather Support to the United States Space Program

This report describes the Applied Meteorology Unit's objective verification of the National Centers for Environmental Prediction 29-km eta model during separate warm and cool season periods from May 1996 through January 1998. The verification of surface and upper-air point forecasts was performed at three selected stations important for 45th Weather Squadron, Spaceflight Meteorology Group, and National Weather Service, Melbourne operational weather concerns. The statistical evaluation identified model biases that may result from inadequate parameterization of physical processes. Since model biases are relatively small compared to the random error component, most of the total model error results from day-to-day variability in the forecasts and/or observations. To some extent, these nonsystematic errors reflect the variability in point observations that sample spatial and temporal scales of atmospheric phenomena that cannot be resolved by the model. On average, Meso-Eta point forecasts provide useful guidance for predicting the evolution of the larger scale environment. A more substantial challenge facing model users in real time is the discrimination of nonsystematic errors that tend to inflate the total forecast error. It is important that model users maintain awareness of ongoing model changes. Such changes are likely to modify the basic error characteristics, particularly near the surface.

Nutter, Paul↗

The EPIC-MOS Particle-Induced Background Spectrum

We have developed a method for constructing a spectrum of the particle-induced instrumental background of the XMM-Newton EPIC MOS detectors that can be used for observations of the diffuse background and extended sources that fill a significant fraction of the instrument field of view. The strength and spectrum of the particle-induced background, that is, the background due to the interaction of particles with the detector and the detector surroundings, is temporally variable as well as spatially variable over individual chips. Our method uses a combination of the filter-wheel-closed data and a database of unexposed-region data to construct a spectrum of the "quiescent" background. We show that, using this method of background subtraction, the differences between independent observations of the same region of "blank sky" are consistent with the statistical uncertainties except when there is clear evidence of solar wind charge exchange emission. We use the blank sky observations to show that contamination by SWCX emission is a strong function of the solar wind proton flux, and that observations through the flanks of the magnetosheath appear to be contaminated only at much higher solar wind fluxes. We have also developed a spectral model of the residual soft proton flares, which allows their effects to be removed to a substantial degree during spectral fitting.

Kuntz, K. D.↗

Supporting Information for ”The 2019-2020 Australian Drought and Bushfires Altered the Partitioning of Hydrological Fluxes

The model configuration employs the modified International Geosphere Biosphere Programme (IGBP) MODIS 20 category landcover data, soil parameters derived from the International Soil Reference and Information Centre, and the Shuttle Radar Topography Mission based elevation, slope, and aspect data. Statistical downscaling approaches are used to transform the coarse resolution MERRA2 meteorological inputs to 1km. The input meteorological fields of air temperature, humidity, surface pressure, wind, downward shortwave radiation, and downward longwave radiation are downscaled to 1km by adjusting for terrain differences in elevation, slope, and aspect. The high resolution monthly precipitation climatology from WorldClim is used to spatially disaggregate input MERRA2 precipitation to 1km. The initial conditions for the model simulations are generated from a long spinup of NoahMP starting from year 2000. All model integrations and evaluations are conducted using the NASA Land Information System and the Land surface Verification Toolkit.

Sujay V Kumar↗

Evaluating Approaches Relating Ecosystem Productivity with Desis Spectral Information

Data from the DLR Earth Sensing Imaging Spectrometer (DESIS), mounted on the International Space Station (ISS), were used to develop and test algorithms for remotely retrieving ecosystem productivity. Twenty DESIS images were used from three widely separated forested study sites representing deciduous and conifer forests. Gross primary production (GPP) values from eddy covariance flux towers at the sites were matched with DESIS spectral reflectances collected on the same days. Multiple algorithms were successful relating spectral reflectance with GPP, including: spectral vegetation indices (SVI) sensitive to chlorophyll content, SVI used in a photosynthetic light-use efficiency model framework, spectral shape characteristics through spectral derivatives and absorption feature analysis, and statistical models leading to multiband hyperspectral indices from partial least squares regression. Successful algorithms were able to achieve R2 better than 0.7 using a diverse set of observations combining data from different sites from multiple years and at multiple times during the year. The demonstrated robustness of the algorithms provides some confidence in using DESIS imagery to map spatial patterns of GPP.

K F Huemmrich↗