Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical change detection”

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 181 records · Page 10

Analysis of Long-Term Quality Control Data for a 137 Cs Dosimetry Calibration Source

Strict quality assurance programs are required for many radiological applications, but these seldom exist for verifying dosimetry calibration sources. After initial characterization of a dosimetry calibration facility, quality control procedures are recommended to ensure the early detection of any changes or malfunctions. These also result in refined knowledge about average dose rate and experimental variations in dose delivery. This paper describes the implementation of a phase I quality control protocol for a 137 Cs dosimetry calibration source and includes an analysis of the resulting data collected over a 24-mo period. During this time, substantial data was collected to establish trial control limits. Air kerma rate measurements were obtained using an ion chamber and were adjusted for decay, corrected for ambient temperature, pressure and humidity, and then analyzed using quality control charts. Three variations of rational subgrouping methods were used in order to find assignable causes of error, and Nelson's Rules were followed to detect any non-random statistical variations. Measurements were subgrouped according to same-day measurements in order to detect positional errors as well as atmospheric correction errors. Additionally, measurements were subgrouped according to analogous experimental setups in order to detect failure in equipment or incorrect settings. Both were analyzed using the X-bar and R chart method. Similarly, individuals and moving ranges charts were used to carefully examine each position in order to observe any situational errors that may occur which include timing, positional, or interference errors. Each method was successful in identifying unique out-of-control data points that occurred during the phase I application of forming control limits. Furthermore, over the 24-mo period, enough data points were deemed in-control to establish reliable trial limits. Future experiments will include the phase II application of gaining more reliable measurements in order to fine-tune the limits, as well as performing a designed experiment, where variables are purposefully changed in order to test the variation of the data.

61 RADIATION PROTECTION AND DOSIMETRY↗

Cross-Layered Distributed Data-Driven Framework for Enhanced Smart Grid Cyber-Physical Security

Smart Grid (SG) research and development has drawn much attention from academia, industry and government due to the great impact it will have on society, economics and the environment. Securing the SG is a considerably significant challenge due the increased dependency on communication networks to assist in physical process control, exposing them to various cyber-threats. In addition to attacks that change measurement values using False Data Injection (FDI) techniques, attacks on the communication network may disrupt the power system's real-time operation by intercepting messages, or by flooding the communication channels with unnecessary data. Addressing these attacks requires a cross-layer approach. In this paper a cross-layered strategy is presented, called Cross-Layer Ensemble CorrDet with Adaptive Statistics(CECD-AS), which integrates the detection of faulty SG measurement data as well as inconsistent network inter-arrival times and transmission delays for more reliable and accurate anomaly detection and attack interpretation. Numerical results show that CECD-AS can detect multiple False Data Injections, Denial of Service (DoS) and Man In The Middle (MITM) attacks with a high F1-score compared to current approaches that only use SG measurement data for detection such as the traditional physics-based State Estimation, Ensemble CorrDet with Adaptive Statistics strategy and other machine learning classification-based detection schemes.

cyber-physical security↗

A study of the feasibility of statistical analysis of airport performance simulation

The feasibility of conducting a statistical analysis of simulation experiments to study airport capacity is investigated. First, the form of the distribution of airport capacity is studied. Since the distribution is non-Gaussian, it is important to determine the effect of this distribution on standard analysis of variance techniques and power calculations. Next, power computations are made in order to determine how economic simulation experiments would be if they are designed to detect capacity changes from condition to condition. Many of the conclusions drawn are results of Monte-Carlo techniques.

Myers, R. H.↗

Volcanic eruptions and solar activity

The historical record of large volcanic eruptions from 1500 to 1980 is subjected to detailed time series analysis. In two weak but probably statistically significant periodicities of about 11 and 80 yr, the frequency of volcanic eruptions increases (decreases) slightly around the times of solar minimum (maximum). Time series analysis of the volcanogenic acidities in a deep ice core from Greenland reveals several very long periods ranging from about 80 to about 350 yr which are similar to the very slow solar cycles previously detected in auroral and C-14 records. Solar flares may cause changes in atmospheric circulation patterns that abruptly alter the earth's spin. The resulting jolt probably triggers small earthquakes which affect volcanism.

Stothers, Richard B.↗

Small-Magnitude Seismic Swarms in Central Utah (US): Interactions of Regional Tectonics, Local Structures and Hydrothermal Systems

Swarms in Central Utah are situated in the complex transition between the Basin and Range (BR) province and the Colorado Plateau. Transecting transverse structures, volcanic deposits, and hydrothermal systems complicate the extensional BR horst and graben structures and provide a multitude of plausible triggering mechanisms. Revisiting the catalog of the University of Utah Seismograph Stations (1981–2022), we analyze spatio-temporal patterns and characteristic features of seismic sequences. Swarms with alternating seismicity rates, bursts, and longer swarms with persistent moment release exhibit a remarkable diversity in temporal evolution. Swarm durations do not scale with cumulative seismic moment: swarms lasting less than 1 day can have similar cumulative seismic moments as month-long swarms. Here, we observe stationary swarms re-occurring for years (e.g., Mineral Mountains), as well as singular swarms in low-seismicity areas (e.g., activating a local structure). The swarms show a pronounced heterogeneity in triggering and driving mechanisms, observed in the detailed analysis of exemplary sequences (detections, relocations, moment tensors, waveform-based clustering, and repeater analysis). The 2022 Sevier Valley sequence activated a BR-related normal fault, the first resolved fault plane in the valley since 1983. The 2011 Circleville sequence is interpreted as a swarm triggered by mainshock-aftershock activity characterized by increasing magnitudes, changing rupture mechanisms, and a concentration of highly similar events in the second part of the sequence. By jointly discussing exemplary sequences and catalog statistics, we draw a comprehensive picture of swarm activity and its relation to geothermal and tectonic activity.

58 GEOSCIENCES↗

Planck 2018 results

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E -mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and optimal modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following final results: f NL local = -0.9 ± 5.1; f NL equil = -26 ± 47; and f NL ortho = -38 ± 24 (68% CL, statistical). These results include low-multipole (4 ≤ ℓ < 40) polarization data that are not included in our previous analysis. The results also pass an extensive battery of tests (with additional tests regarding foreground residuals compared to 2015), and they are stable with respect to our 2015 measurements (with small fluctuations, at the level of a fraction of a standard deviation, which is consistent with changes in data processing). Polarization-only bispectra display a significant improvement in robustness; they can now be used independently to set primordial NG constraints with a sensitivity comparable to WMAP temperature-based results and they give excellent agreement. In addition to the analysis of the standard local, equilateral, and orthogonal bispectrum shapes, we consider a large number of additional cases, such as scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is, however, detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 σ . Beyond estimates of individual shape amplitudes, we also present model-independent reconstructions and analyses of the Planck CMB bispectrum. Our final constraint on the local primordial trispectrum shape is g NL local = (-5.8 ± 6.5) × 10 4 (68% CL, statistical), while constraints for other trispectrum shapes are also determined. Exploiting the tight limits on various bispectrum and trispectrum shapes, we constrain the parameter space of different early-Universe scenarios that generate primordial NG, including general single-field models of inflation, multi-field models (e.g. curvaton models), models of inflation with axion fields producing parity-violation bispectra in the tensor sector, and inflationary models involving vector-like fields with directionally-dependent bispectra. Our results provide a high-precision test for structure-formation scenarios, showing complete agreement with the basic picture of the ΛCDM cosmology regarding the statistics of the initial conditions, with cosmic structures arising from adiabatic, passive, Gaussian, and primordial seed perturbations.

79 ASTRONOMY AND ASTROPHYSICS↗

Resolution dependence of the turbulent atmospheric boundary layer in global storm-resolving climate simulations

The current generation of state-of-the-art global climate models are being run at increasingly higher horizontal resolutions, with the goal of resolving organised deep convection explicitly. How a kilometre-scale resolution impacts the representation of the atmospheric boundary layer is, however, not well known. Using statistical analysis on global fields as well as high-frequency data at selected locations, produced with the Integrated Forecasting System (IFS) model for the Next Generation Earth-system Models (nextGEMS) project, we investigate the horizontal resolution dependence of some boundary-layer processes. We find that a change in resolution from 9 to 2.8 km causes no substantial changes to boundary-layer properties and processes at most of the locations studied, although some global changes are detected that indicate circulation changes. Small changes to the boundary-layer depth and structure are found in the Tropics. The short simulation length and lack of data for optimal boundary-layer analysis limits the conclusions, especially in relation to the connection between the boundary layer and the atmosphere general circulation.

54 ENVIRONMENTAL SCIENCES↗

Can ERA5 Be Used to Study Mesoscale Convective System Climatological Characteristics?

Mesoscale convective systems (MCSs) produce more than half of tropical rainfall and are central to the global hydrologic cycle. As the climate warms, environments favorable for MCSs may become more common; however, limited observational records hamper understanding of how MCSs respond to variations and changes in their environments. Here, we evaluate how well MCSs are represented in ERA5, a widely used global high‐resolution reanalysis product. Using PyFLEXTRKR, which jointly tracks top‐of‐atmosphere infrared brightness temperature and surface precipitation, we identified MCSs in ERA5 and compared them with those identified in satellite observations using the same detection algorithm. This comparison analysis spans 2007–2020 using hourly data at 0.25° horizontal resolution focusing over the tropics. ERA5 reproduces observed brightness‐temperature statistics and captures the geographic distribution and seasonal and diurnal cycles of MCS cold cloud shields. However, ERA5 precipitation exhibits an intensity bias—too much light rain and too little heavy rain—which shifts the rain‐rate distribution and reduces the frequency of MCSs relative to observations. Within MCSs, ERA5 precipitation exhibits the same pattern of bias, yielding a systematic underestimation of MCS precipitation intensity. Consistent with these biases, ERA5 underestimates the contribution of MCS to tropical rainfall by 25%–34% in key regions. Overall, ERA5 is suitable for studying MCS cold cloud‐shield climatology and evolution, but precipitation‐based MCS characteristics (including event‐level precipitation features and the geospatial distribution of MCS precipitation) should be interpreted with caution. These findings clarify which aspects of MCS behavior are robustly represented in ERA5 for climatological applications.

mesoscale convection↗

First observations of solar halo gamma rays over a full solar cycle

We analyze 15 years of Fermi-LAT data and produce a detailed model of the Sun’s inverse-Compton scattering emission (solar halo), which is powered by interactions between ambient cosmic-ray electrons and positrons with sunlight. By developing a novel analysis method to analyze moving sources, we robustly detect the solar halo at energies between 31.6 MeV and 100 GeV, and angular extensions up to 45° from the Sun, providing new insight into spatial regions where there are no direct measurements of the Galactic cosmic-ray flux. The large statistical significance of our signal allows us to subdivide the data and provide the first 𝛾-ray probes into the time variation and azimuthal asymmetry of the solar modulation potential, finding time-dependent changes in solar modulation both parallel and perpendicular to the ecliptic plane. Our results are consistent with (but with independent uncertainties from) local cosmic-ray measurements, unlocking new probes into astrophysical processes near the solar surface.

79 ASTRONOMY AND ASTROPHYSICS↗

Longitudinal analysis of host protein serum signatures of treatment and recovery in pulmonary tuberculosis

A better understanding of treatment progression and recovery in pulmonary tuberculosis (TB) infectious disease is crucial. This study analyzed longitudinal serum samples from pulmonary TB patients undergoing interventional treatment to identify surrogate markers for TB-related outcomes. Serum that was collected at baseline and 8, 17, 26, and 52 weeks from 30 TB patients experiencing durable cure were evaluated and compared using a sensitive LC-MS/MS proteomic platform for the detection and quantification of differential host protein signatures relative to timepoint. The global proteome signature was analyzed for statistical differences across the time course and between disease severity and treatment groups. A total of 676 proteins showed differential expression in the serum over these timepoints relative to baseline. Comparisons to understand serum protein dynamics at 8 weeks, treatment endpoints at 17 and 26 weeks, and post-treatment at 52 weeks were performed. The largest protein abundance changes were observed at 8 weeks as the initial effects of antibiotic treatment strongly impacted inflammatory and immune modulated responses. However, the largest number of proteome changes was observed at the end of treatment time points 17 and 26 weeks respectively. Post-treatment 52-week results showed an abatement of differential proteome signatures from end of treatment, though interestingly those proteins uniquely significant at post-treatment were almost exclusively downregulated. Patients were additionally stratified based upon disease severity and compared across all timepoints, identifying 461 discriminating proteome signatures. These proteome signatures collapsed into discrete expression profiles with distinct pathways across immune activation and signaling, hemostasis, and metabolism annotations. Insulin-like growth factor (IGF) and Integrin signaling maintained a severity signature through 52 weeks, implying an intrinsic disease severity signature well into the post-treatment timeframe. Previous proteome studies have primarily focused on the 8-week timepoint in relation to culture conversion status. While this study confirms previous observations, it also highlights some differences. The inclusion of additional end of treatment and post-treatment time points offers a more comprehensive assessment of treatment progression within the serum proteome. Examining the expression dynamics at these later time periods will help in the investigation of relapse patients and has provided indicative markers of response and recovery.

59 BASIC BIOLOGICAL SCIENCES↗

Searching for hard X-ray directivity during the rise, peak, and decay phases of solar flares

We have identified 72 large solar flares (peak counting rates more than 1000 counts/s) observed by Hard X-ray Burst Spectroscopy (HXRBS) on-board the Solar Maximum Mission (SMM). Using a database of these flares, we have studied hard X-ray (50-850 keV) spectral center-to-limb variation and its evolution with time. The major results are the following: (1) During the rise phase, the center-to-limb spectral variation is small, with a hardness of delta delta = 0.02 +/- 0.25, and a statistical significance of 0.1 sigma. (2) During the peak phase, the center-to-limb variation is delta delta = 0.13 +/- 0.13, with a statistical significance of 1 sigma. (3) During the decay phase, the center-to-limb variation changes to softening. The softness is relatively large with delta delta = -0.25 +/- 0.21, and a statistical significance of 1.2 sigma. (4) The linear least-squares fits to the spectral center-to-limb variations do not have slopes significantly different from zero during all those three phases. (5) The center events and limb events spectral distributions are shown to be not different by using Kolmogorov-Smirnov two-samples test. (6) The fraction of events detected near the limb is marginally consistent with that expected from isotropically emitting flares. (7) On average, flares evolve as soft-hard-soft. These results suggest that there is no statistically significant evidence for hard X-ray directivity during the rise, peak, and decay phases of solar flares. The hard X-ray radiation pattern at those energies is almost isotropic during all those phases. This lack of directivity (or anisotropy) found in this study is not in agreement with the results discovered by Vestrand et al. (1987) in which they found energetic photon source is anisotropic, using SMM Gamma-Ray Spectrometer (GRS) data at a much higher energy band of 0.3-1 MeV. If we want to interpret the results of Vestrand et al. (1987) and our present results in a self-consistent way, we must conclude that at energies of greater than 50 keV, the Compton backscattering suppresses the hard X-ray anisotropy sufficiently so that no significant directivity can be measured during the rise, peak, and decay phases of solar flares. during the rise, peak, and decay phases of solar flares.

Li, Peng↗

Optimizing Science Return from Plato's Observations of the Kepler Field

The Plato mission could substantially enhance its scientific yield by focusing a significant fraction of its observations with all 24 cameras viewing the Kepler Field of View (FOV). Kepler, while exceptionally productive, did not accomplish its goal of providing a direct measurement of the abundance of Earth-analog planet around sunlike stars. On this issue, Plato observations centered on the Kepler FOV would enable planet searches using the combined data sets to achieve more precise estimates of planetary occurrence rates than obtainable using data from either mission alone. The ecliptic latitude of the Kepler FOV is too low for Plato to observe with 24 cameras year-round, so current plan’s for Plato’s Northern Stare places the Kepler FOV in the periphery, reducing the number of photons collected from Kepler targets by factors of 2 or 4. We propose instead four ~6 - 8 month Plato stares with 24 cameras observing the Kepler FOV during the Plato primary mission, instead of a single 2 year stare with the Kepler FOV located to the periphery. This change would greatly increase the amount of 'Plato glass’ dedicated to the continued study of the Kepler FOV, while still enabling Plato to observe other parts of the sky during the majority of its prime mission. In addition to discovering planets below Kepler's detection threshold, this change would yield transit follow-up of thousands of known exoplanets with an effective throughput of a 0.6 m space telescope. It would also extend the observing baseline of the Kepler FOV to timescales of decades, providing key data for stellar cycles, the analysis of transit timing variations, and other decade-timescale variations. Such a campaign enhances the likelihood of the detection of true Earth-analog transiting planet candidates around sunlike stars, provide the best statistical study of exoplanet systems with transit timing variations and duration variations, and open up the study of stellar variability to new time-domains with precision on large numbers of faint targets that only 24 Plato cameras can provide.

Jack Lissauer↗

Preventing Failures By Dataset Shift Detection in Safety-Critical Graph Applications

Dataset shift refers to the problem where the input data distribution may change over time (e.g., between training and test stages). Since this can be a critical bottleneck in several safety-critical applications such as healthcare, drug-discovery, etc., dataset shift detection has become an important research issue in machine learning. Though several existing efforts have focused on image/video data, applications with graph-structured data have not received sufficient attention. Therefore, in this paper, we investigate the problem of detecting shifts in graph structured data through the lens of statistical hypothesis testing. Specifically, we propose a practical two-sample test based approach for shift detection in large-scale graph structured data. Our approach is very flexible in that it is suitable for both undirected and directed graphs, and eliminates the need for equal sample sizes. Using empirical studies, we demonstrate the effectiveness of the proposed test in detecting dataset shifts. We also corroborate these findings using real-world datasets, characterized by directed graphs and a large number of nodes.

97 MATHEMATICS AND COMPUTING↗

Automated thematic mapping and change detection of ERTS-A images

The author has identified the following significant results. For the recognition of terrain types, spatial signatures are developed from the diffraction patterns of small areas of ERTS-1 images. This knowledge is exploited for the measurements of a small number of meaningful spatial features from the digital Fourier transforms of ERTS-1 image cells containing 32 x 32 picture elements. Using these spatial features and a heuristic algorithm, the terrain types in the vicinity of Phoenix, Arizona were recognized by the computer with a high accuracy. Then, the spatial features were combined with spectral features and using the maximum likelihood criterion the recognition accuracy of terrain types increased substantially. It was determined that the recognition accuracy with the maximum likelihood criterion depends on the statistics of the feature vectors. Nonlinear transformations of the feature vectors are required so that the terrain class statistics become approximately Gaussian. It was also determined that for a given geographic area the statistics of the classes remain invariable for a period of a month but vary substantially between seasons.

Gramenopoulos, N.↗

Star Tracker Accuracy Improvement and Optimization for Attitude Measurement in Three-Axis

High precision attitude measurement systems obviate the need for the beacon from the receiver making it possible for the spacecraft to beam a laser communications signal to a ground station without the ground station advertising its location. The research presented targets new detection and estimation methods to improve the accuracy in locating stars on a focal plane detector, and an understanding of the effects of changes in the optics design parameters and aberration, including defocus, on the navigation solution itself. This understanding can lead to an optimization of the attitude solution with respect to those optics realm parameter changes. The methodology discussed includes the development of a model of a current star tracker system. Using this model, multiple algorithms are implemented, including a multi-hypothesis method (MHT), to detect and estimate the position of the stars on the focal plane detector. It will be shown that using the MHT for detection and estimation, a greater accuracy can be found for each star estimation from more traditional detection and estimation algorithms. The approach then uses the model to develop statistics of the star tracker and the attitude estimation outputs to understand the accuracy, or variance, of the system's attitude solution. This solution is repeated for a range of defocus aberration, and a lower limit to the variance of the attitude solution is shown. A Cramer Rao lower bound solution is derived for the star tracker system and the results are compared to the Monte Carlo analysis from the model and shown to correlate very well. The approach uses a star image not as a Gaussian spot on the focal plane as done in previous work, and use of an image that includes the effects of aberrations of the optic system, and the effects of under-sampling and noise from the focal plane detector as well. Analysis includes exploring a star tracker's accuracy improvement through the combination of focus error and under-sampling effects alone, possibly contradicting conventional wisdom and approaches.

Lichter, Michael J.↗

Soft Decision Analyzer

The Soft Decision Analyzer (SDA) is an instrument that combines hardware, firmware, and software to perform realtime closed-loop end-to-end statistical analysis of single- or dual- channel serial digital RF communications systems operating in very low signal-to-noise conditions. As an innovation, the unique SDA capabilities allow it to perform analysis of situations where the receiving communication system slips bits due to low signal-to-noise conditions or experiences constellation rotations resulting in channel polarity in versions or channel assignment swaps. SDA s closed-loop detection allows it to instrument a live system and correlate observations with frame, codeword, and packet losses, as well as Quality of Service (QoS) and Quality of Experience (QoE) events. The SDA s abilities are not confined to performing analysis in low signal-to-noise conditions. Its analysis provides in-depth insight of a communication system s receiver performance in a variety of operating conditions. The SDA incorporates two techniques for identifying slips. The first is an examination of content of the received data stream s relation to the transmitted data content and the second is a direct examination of the receiver s recovered clock signals relative to a reference. Both techniques provide benefits in different ways and allow the communication engineer evaluating test results increased confidence and understanding of receiver performance. Direct examination of data contents is performed by two different data techniques, power correlation or a modified Massey correlation, and can be applied to soft decision data widths 1 to 12 bits wide over a correlation depth ranging from 16 to 512 samples. The SDA detects receiver bit slips within a 4 bits window and can handle systems with up to four quadrants (QPSK, SQPSK, and BPSK systems). The SDA continuously monitors correlation results to characterize slips and quadrant change and is capable of performing analysis even when the receiver under test is subjected to conditions where its performance degrades to high error rates (30 percent or beyond). The design incorporates a number of features, such as watchdog triggers that permit the SDA system to recover from large receiver upsets automatically and continue accumulating performance analysis unaided by operator intervention. This accommodates tests that can last in the order of days in order to gain statistical confidence in results and is also useful for capturing snapshots of rare events.

Steele, Glen↗

Identifying environmental features for land management decisions

The benefits of changes in management organization and facilities for the Center for Remote Sensing and Cartography in Utah are reported as well as interactions with and outreach to state and local agencies. Completed projects are described which studied (1) Unita Basin wetland/land use; (2) Davis County foothill development; (3) Farmington Bay shoreline fluctuation; (4) irrigation detection; and (5) satellite investigation of snow cover/mule deer relationships. Techniques developed for composite computer mapping, contrast enhancement, U-2 CIR/LANDSAT digital interface; factor analysis, and multivariate statistical analysis are described.

Source record↗

Exploring variability in seasonal average and extreme precipitation using unsupervised machine learning.

Focal Area(s): We will use unsupervised machine learning methods to identify and quantify the influence of large scale natural modes of climate variability to gain insight into the observed and simulated seasonal average and extreme precipitation changes. Science Challenge: A recent paper, led by co-PI Mark Risser, finds that although much of the variability in seasonal average and extreme precipitation over CONUS is unforced, the effect of large-scale modes of circulation variability (such as ENSO, AMO, PNA, etc.) can be detected and attributed. However, it is unclear whether or not unsupervised learning methods can (a) replicate this finding or (b) yield insight into possible nonlinear behavior that was not captured in the initial statistical analysis. Further work would entail extending this framework to other global land areas.

54 ENVIRONMENTAL SCIENCES↗