Search NASA⌕ Search

SEARCH · Search NASA

Results for “fault estimation”

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 703 records · Page 39

An analysis of ridge axis gravity roughness and spreading rate

Fast and slow spreading ridges have radically different morphologic and gravimetric characteristics. In this study, altimeter measurements from the Geosat Exact Repeat Mission are used to investigate spreading rate dependence of the ridge axis gravity field. Gravity roughness provides an estimate of the amplitude of the gravity anomaly and is robust to small errors in the location of the ridge axis. Gravity roughness as a weighted root mean square of the vertical deflection at 438 ridge crossings on the mid-ocean ridge system is computed. Ridge axis gravity anomalies show a decrease in amplitude with increasing spreading rate up to an intermediate rate of about 60-80 mm/yr and almost no change at higher rates; overall the roughness decreases by a factor of 10 between the lowest and highest rates. In addition to the amplitude decrease, the range of roughness values observed at a given spreading rate shows a similar order of magnitude decrease with transition between 60 and 80 mm/yr. The transition of ridge axis gravity is most apparent at three relatively unexplored locations on the Southeast Indian Ridge and the Pacific-Antarctic Rise; on these intermediate rate ridges the transition occurs abruptly across transform faults.

Small, Christopher↗

Instrumentation System Diagnoses a Thermocouple

An improved self-validating thermocouple (SVT) instrumentation system not only acquires readings from a thermocouple but is also capable of detecting deterioration and a variety of discrete faults in the thermocouple and its lead wires. Prime examples of detectable discrete faults and deterioration include open- and short-circuit conditions and debonding of the thermocouple junction from the object, the temperature of which one seeks to measure. Debonding is the most common cause of errors in thermocouple measurements, but most prior SVT instrumentation systems have not been capable of detecting debonding. The improved SVT instrumentation system includes power circuitry, a cold-junction compensator, signal-conditioning circuitry, pulse-width-modulation (PWM) thermocouple-excitation circuitry, an analog-to-digital converter (ADC), a digital data processor, and a universal serial bus (USB) interface. The system can operate in any of the following three modes: temperature measurement, thermocouple validation, and bonding/debonding detection. The software running in the processor includes components that implement statistical algorithms to evaluate the state of the thermocouple and the instrumentation system. When the power is first turned on, the user can elect to start a diagnosis/ monitoring sequence, in which the PWM is used to estimate the characteristic times corresponding to the correct configuration. The user also has the option of using previous diagnostic values, which are stored in an electrically erasable, programmable read-only memory so that they are available every time the power is turned on.

Perotti, Jose↗

A dual-processor multi-frequency implementation of the FINDS algorithm

This report presents a parallel processing implementation of the FINDS (Fault Inferring Nonlinear Detection System) algorithm on a dual processor configured target flight computer. First, a filter initialization scheme is presented which allows the no-fail filter (NFF) states to be initialized using the first iteration of the flight data. A modified failure isolation strategy, compatible with the new failure detection strategy reported earlier, is discussed and the performance of the new FDI algorithm is analyzed using flight recorded data from the NASA ATOPS B-737 aircraft in a Microwave Landing System (MLS) environment. The results show that low level MLS, IMU, and IAS sensor failures are detected and isolated instantaneously, while accelerometer and rate gyro failures continue to take comparatively longer to detect and isolate. The parallel implementation is accomplished by partitioning the FINDS algorithm into two parts: one based on the translational dynamics and the other based on the rotational kinematics. Finally, a multi-rate implementation of the algorithm is presented yielding significantly low execution times with acceptable estimation and FDI performance.

Godiwala, Pankaj M.↗

Flexural rigidity of the Basin and Range-Colorado Plateau-Rocky Mountain transition from coherence analysis of gravity and topography

Stochastic inversion for flexural loads and flexural rigidity of the continental elastic layer can be accomplished most effectively by using the coherence of gravity and topography. However, the spatial resolution of coherence analysis has been limited by use of two-dimensional periodogram spectra from very large (greater than 10(exp 5)sq km) windows that generally include multiple tectonic features. Using a two-dimensional spectral estimator based on the maximum entropy method, the spatial resolution of flexural proerties can be enhanced by a factor of 4 or more, enabling more detailed analysis at the scale of individual tectonic features. This new approach is used to map the spatial variation of flexural rigidity along the Basin and Range transition to the Colorado Plateau and Middle Rocky Mountains physiographic provinces. Large variations in flexural isostatic responses are found, with rigidities ranging from as low as 8.7 x 10(exp 20) N m (elastic thickness (T(sub e) = 4.6 km) in the Basin and Range to as high as 4.1 x 10(exp 24) N m T(sub e) = 77 km) in the Middle Rocky Mountains. These results compare favorably woith independent determinations of flexural rigidity in the region. Areas of low flexural rigidity correlate strongly with areas of high surface heat flow, as is expected from the contingence of flexural rigidity on a temperature-dependent flow law. Also, late Cenozoic normal faults with large displacements are found primarily in area of low flexural rigidity region. The highest flexural rigidity is found within the Archean Wyoming craton, where evidence suggests that deeply rooted cratonic lithosphere may play a role in determining the distribution of tectonism at the surface.

Lowery, Anthony R.↗

Research Supporting Satellite Communications Technology

This report describes the second year of research effort under the grant Research Supporting Satellite Communications Technology. The research program consists of two major projects: Fault Tolerant Link Establishment and the design of an Auto-Configurable Receiver. The Fault Tolerant Link Establishment protocol is being developed to assist the designers of satellite clusters to manage the inter-satellite communications. During this second year, the basic protocol design was validated with an extensive testing program. After this testing was completed, a channel error model was added to the protocol to permit the effects of channel errors to be measured. This error generation was used to test the effects of channel errors on Heartbeat and Token message passing. The C-language source code for the protocol modules was delivered to Goddard Space Flight Center for integration with the GSFC testbed. The need for a receiver autoconfiguration capability arises when a satellite-to-ground transmission is interrupted due to an unexpected event, the satellite transponder may reset to an unknown state and begin transmitting in a new mode. During Year 2, we completed testing of these algorithms when noise-induced bit errors were introduced. We also developed and tested an algorithm for estimating the data rate, assuming an NRZ-formatted signal corrupted with additive white Gaussian noise, and we took initial steps in integrating both algorithms into the SDR test bed at GSFC.

Horan Stephen↗

Physics informed neural network can retrieve rate and state friction parameters from acoustic monitoring of laboratory stick-slip experiments

Various machine learning (ML) and deep learning (DL) techniques have been recently applied to the forecasting of laboratory earthquakes from friction experiments. The magnitude and timing of shear failures in stick-slip cycles are predicted using features extracted from the recorded ultrasonic or acoustic emission (AE) signals. In addition, the Rate and State Friction (RSF) constitutive laws are extensively used to model the frictional behavior of faults. In this work, we use data from shear experiments coupled with passive acoustic (variance, kurtosis, and AE rate) interleaved with active source ultrasonic monitoring (transmitted wave amplitude) to develop physics-informed neural network (PINN) models incorporating the RSF law and AE rate generation equation with wave amplitude serving as a proxy for friction state variable. This PINN framework allows learning RSF parameters from stick-slip experiments rather than measuring them through a series of velocity step experiments. We observe that when the stick-slip cycles are irregular, the PINN models outperform the data-driven DL models. Transfer learning (TL) PINN models are also developed by pre-training on data collected at one normal stress level followed by forecasting shear failures and retrieving RSF parameters at other stress levels (i.e., with different recurrence intervals) after retraining on a limited amount of new data. Our findings suggest that TL models perform better compared to standalone models. Both standalone and TL PINN-estimated RSF parameters and their ground truth values show excellent agreements thus demonstrating that RSF parameters can be retrieved from laboratory stick-slip experiments using the corresponding acoustic data and that the transmitted wave amplitude provides a good representation of the evolving frictional state during stick-slips.

58 GEOSCIENCES↗

Efficient Measurement-Driven Eigenenergy Estimation with Classical Shadows

Quantum algorithms exploiting real-time evolution under a target Hamiltonian have demonstrated remarkable efficiency in extracting key spectral information. However, the broader potential of these methods, particularly beyond ground-state calculations, is underexplored. In this work, we introduce the framework of multiobservable dynamic mode decomposition (MODMD), which combines the observable dynamic mode decomposition (DMD), a measurement-driven eigensolver tailored for near-term implementation, with classical shadow tomography. MODMD leverages random scrambling in the classical shadow technique to construct, with exponentially reduced resource requirements, a signal subspace that encodes rich spectral information. Notably, we replace typical Hadamard-test circuits with a protocol designed to predict low-rank observables, thereby broadening the use of classical shadow tomography for predicting many low-rank observables. We establish theoretical guarantees on the spectral approximation from MODMD, taking into account distinct sources of error. In the ideal case, we prove that the spectral error scales as exp (−Δ⁢𝐸⁢𝑡 max ), where Δ⁢𝐸 is the Hamiltonian spectral gap and 𝑡 max is the maximal simulation time. This analysis provides a rigorous justification of the rapid convergence observed across simulations. To demonstrate the utility of our framework, we consider its application to fundamental tasks, such as determining the low-lying, i.e., ground or excited, energies of representative many-body systems. Our work paves the path for efficient designs of measurement-driven algorithms on near-term and early fault-tolerant quantum devices.

quantum algorithms & computation↗

Profiling the Quito basin (Ecuador) using seismic ambient noise

SUMMARY Quito, the capital of Ecuador, with more than 2.5 M inhabitants, is exposed to a high seismic hazard due to its proximity to the Pacific subduction zone and active crustal faults, both capable of generating significant earthquakes. Furthermore, the city is located in an intermontane piggy-back basin prone to seismic wave amplification. To understand the basin’s seismic response and characterize its geological structure, 20 broad and medium frequency band seismic stations were deployed in Quito’s urban area between May 2016 and July 2018 that continuously recorded ambient seismic noise. We first compute horizontal-to-vertical spectral ratios to determine the resonant frequency distribution in the entire basin. Secondly, we cross-correlate seismic stations operating simultaneously to retrieve interstations surface-wave Green’s functions in the frequency range of 0.1–2 Hz. We find that Love waves travelling in the basin’s longitudinal direction (NNE–SSW) show much clearer correlograms than those from Rayleigh waves. We then compute Love wave phase-velocity dispersion curves and invert them in conjunction with the HVSR curves to obtain shear-wave velocity profiles throughout the city. The inversions highlight a clear difference in the basin’s structure between its northern and southern parts. In the centre and northern areas, the estimated basin depth and mean shear-wave velocity are about 200 m and 1800 ms−1, respectively, showing resonance frequency values between 0.6 and 0.7 Hz. On the contrary, the basement’s depth and shear-wave velocity in the southern part are about 900 m and 2500 ms−1, having a low resonance frequency value of around 0.3 Hz. This difference in structure between the centre-north and the south of the basin explains the spatial distribution of low-frequency seismic amplifications observed during the Mw 7.8 Pedernales earthquake in April 2016 in Quito.

Pacheco, D. (ORCID:0000000282355096)↗

Detecting Short Circuits: Post Accident Electric Vehicle Battery Safety Check

Fast and accurate detection of soft short circuits (SCs) in the battery packs of damaged electric vehicles is needed by first responders and mechanics to mitigate the potential risk from battery fires that may occur hours, days, or weeks after an accident. Here, this paper presents an SC-detection algorithm for potentially damaged lithium-ion batteries that works quickly and without a priori knowledge of the battery-pack chemistry, capacity, state of charge, or state of health. The proposed universal SC-detection algorithm is designed to be implemented on an inexpensive handheld device that can connect to and monitor the voltages of all cells in a pack. Transient filtering and linear-quadratic state observation provide estimates of normalized SC current for every cell in the pack. Cells with SC-current estimates outside a sigma-based threshold are detected. Simulations, experiments, and electric vehicle (EV) crash data are used to verify the speed, sensitivity, and accuracy of the method, demonstrating 96% accurate detection of 0.0027 C SCs in under 1 h for 5S cell groups in the lab and no false positives for crashed Volkswagen, Chevrolet, and Tesla vehicles without SCs.

25 - ENERGY STORAGE↗

The effects of history and predictive information on the ability of the transport aircraft pilot to predict an alert

The early detection of a developing aircraft-subsystem fault has the potential to lessen its ultimate severity. The lack of capability for such early detection is becoming critical in the aviation community. In the commercial sector, for example, twin-engine aircraft are being used for extended transport operations over water. One method to decrease the severity of a developing problem is to predict its behavior and to take appropriate corrective action. In order to investigate pilots' ability to predict events, an experiment was conducted where eighteen airline pilots predicted the time to an alert using three different displays of dials and three different time profile complexities. The three displays of dials were as follows: standard, resembling current aircraft dial presentations; history, indicating the value five seconds in the past; and predictive, indicating the value five seconds into the future. The time profiles describing the behavior of the parameter consisted of constant velocity profiles, decelerating profiles, and accelerating then decelerating profiles. Although pilots indicated that they preferred the predictive dial, the objective data did not support its use. The objective data did show that the time profiles had the most significant effect on performance in estimating the time to an alert.

Trujillo, Anna C.↗

Mineralogy and composition of Archean Crust, Greenland: A pilot study

The Portable Instant Display and Analysis Spectrometer (PIDAS) was taken to southwestern Greenland to investigate in situ the potential application of AVIRIS to estimate the mineralogy and composition of rocks exposed in Archean terranes. The goal was to determine the feasibility of using a high spectral resolution scanner to find and study pristine rocks, those that have not been altered by subsequent deformation and metamorphism. The application of AVIRIS data to the problems in Greenland is logical. However, before a costly deployment of the U-2 aircraft to Greenland is proposed, this study was undertaken to acquire the spectral data necessary to verify that mineralogical mapping in the environmental conditions found there is possible. Although field conditions were far from favorable, all of the major objectives of the study were addressed. One of the major concerns was that lichens would obscure the rock surfaces. It was found that the spectral signature of the lichens was distinct from the underlying rocks. Thus, a spectrum of a rock outcrop, with its partial cover of lichens, can be un-mixed into rock and lichen components. The data acquired during the course of this study supports the conclusion that areas of pristine Archean crust can be differentiated from that which has experienced low grade alteration associated with Proterizoic faulting.

Goetz, Alexander F. H.↗

Classification and regression models of audio and vibration signals for machine state monitoring in precision machining systems

Here we present a data-driven method for monitoring machine status in manufacturing processes. Audio and vibration data from precision machining are used for inference in two operating scenarios: (a) variable machine health states (anomaly detection); and (b) settings of machine operation (state estimation). Audio and vibration signals are first processed through Fast Fourier Transform and Principal Component Analysis to extract transformed and informative features. These features are then used in the training of classification and regression models for machine state monitoring. Specifically, three classifiers (K-nearest neighbors, convolutional neural networks and support vector machines) and two regressors (support vector regression and neural network regression) were explored, in terms of their accuracy in machine state prediction. It is shown that the audio and vibration signals are sufficiently rich in information about the machine that 100% state classification accuracy could be accomplished. Data fusion was also explored, showing overall superior accuracy of data-driven regression models.

42 ENGINEERING↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optormechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by, the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Wrinkle ridges on Mercury and the Moon within and outside of mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER)spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 min relief with a mean of ~350m (median=~340m, n=150) on Mercury and from ~47 to 678 min relief with a mean of ~198 m (median=~168m, n=150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Lisa S. Schleicher↗

The Single Event Effect Characteristics of the 486-DX4 Microprocessor

This research describes the development of an experimental radiation testing environment to investigate the single event effect (SEE) susceptibility of the 486-DX4 microprocessor. SEE effects are caused by radiation particles that disrupt the logic state of an operating semiconductor, and include single event upsets (SEU) and single event latchup (SEL). The relevance of this work can be applied directly to digital devices that are used in spaceflight computer systems. The 486-DX4 is a powerful commercial microprocessor that is currently under consideration for use in several spaceflight systems. As part of its selection process, it must be rigorously tested to determine its overall reliability in the space environment, including its radiation susceptibility. The goal of this research is to experimentally test and characterize the single event effects of the 486-DX4 microprocessor using a cyclotron facility as the fault-injection source. The test philosophy is to focus on the "operational susceptibility," by executing real software and monitoring for errors while the device is under irradiation. This research encompasses both experimental and analytical techniques, and yields a characterization of the 486-DX4's behavior for different operating modes. Additionally, the test methodology can accommodate a wide range of digital devices, such as microprocessors, microcontrollers, ASICS, and memory modules, for future testing. The goals were achieved by testing with three heavy-ion species to provide different linear energy transfer rates, and a total of six microprocessor parts were tested from two different vendors. A consistent set of error modes were identified that indicate the manner in which the errors were detected in the processor. The upset cross-section curves were calculated for each error mode, and the SEU threshold and saturation levels were identified for each processor. Results show a distinct difference in the upset rate for different configurations of the on-chip cache, as well as proving that one vendor is superior to the other in terms of latchup susceptibility. Results from this testing were also used to provide a mean-time-between-failure estimate of the 486-DX4 operating in the radiation environment for the International Space Station.

Kouba, Coy↗

Wrinkle Ridges on Mercury and the Moon Within and Outside of Mascons

Found on all terrestrial planets, wrinkle ridges are anticlines formed by thrust faulting and folding resulting from crustal shortening. The MErcury Surface, Space Environment, Geochemistry, and Ranging (MESSENGER) spacecraft's orbital phase returned high resolution images and topographic data of the previously unimaged northern high latitudes of Mercury where there are large expanses of smooth plains deformed by wrinkle ridges. Concurrently, the Lunar Reconnaissance Orbiter (LRO) is obtaining high resolution images and topographic data covering lunar mare wrinkle ridges. These data allow quantitative comparison of the scale of wrinkle ridges in smooth plains volcanic units on Mercury with mare wrinkle ridges. We evaluate the topographic relief of 300 wrinkle ridges within and outside of mascon basins on the Moon and Mercury. Measured wrinkle ridges range from ~112 to 776 m in relief with a mean of ~350 m (median = ~340 m, n = 150) on Mercury and from ~47 to 678 m in relief with a mean of ~198 m (median = ~168 m, n = 150) on the Moon. Wrinkle ridges on Mercury thus are approximately twice as large in mean relief compared to their counterparts on the Moon. The larger scale of Mercury's wrinkle ridges suggests that their formation can be attributed, in part, to global contraction. As global contraction on the Moon is estimated to be an order of magnitude smaller than on Mercury, the smaller scale of lunar wrinkle ridges suggests they most likely form primarily by load induced subsidence of the mare basalt. Wrinkle ridges located in lunar mascon basins and in the Caloris mascon on Mercury are not statistically significantly different in relief than ridges in non-mascon regions, suggesting comparable levels of contractional strain. The fact that mascon basins do not host wrinkle ridges with greater structural relief relative to non-mascon units may indicate the critical role lithospheric thickness plays in controlling subsidence and contraction of thick volcanic sequences on the Moon and Mercury.

Global contraction↗

Understanding software faults and their role in software reliability modeling

This study is a direct result of an on-going project to model the reliability of a large real-time control avionics system. In previous modeling efforts with this system, hardware reliability models were applied in modeling the reliability behavior of this system. In an attempt to enhance the performance of the adapted reliability models, certain software attributes were introduced in these models to control for differences between programs and also sequential executions of the same program. As the basic nature of the software attributes that affect software reliability become better understood in the modeling process, this information begins to have important implications on the software development process. A significant problem arises when raw attribute measures are to be used in statistical models as predictors, for example, of measures of software quality. This is because many of the metrics are highly correlated. Consider the two attributes: lines of code, LOC, and number of program statements, Stmts. In this case, it is quite obvious that a program with a high value of LOC probably will also have a relatively high value of Stmts. In the case of low level languages, such as assembly language programs, there might be a one-to-one relationship between the statement count and the lines of code. When there is a complete absence of linear relationship among the metrics, they are said to be orthogonal or uncorrelated. Usually the lack of orthogonality is not serious enough to affect a statistical analysis. However, for the purposes of some statistical analysis such as multiple regression, the software metrics are so strongly interrelated that the regression results may be ambiguous and possibly even misleading. Typically, it is difficult to estimate the unique effects of individual software metrics in the regression equation. The estimated values of the coefficients are very sensitive to slight changes in the data and to the addition or deletion of variables in the regression equation. Since most of the existing metrics have common elements and are linear combinations of these common elements, it seems reasonable to investigate the structure of the underlying common factors or components that make up the raw metrics. The technique we have chosen to use to explore this structure is a procedure called principal components analysis. Principal components analysis is a decomposition technique that may be used to detect and analyze collinearity in software metrics. When confronted with a large number of metrics measuring a single construct, it may be desirable to represent the set by some smaller number of variables that convey all, or most, of the information in the original set. Principal components are linear transformations of a set of random variables that summarize the information contained in the variables. The transformations are chosen so that the first component accounts for the maximal amount of variation of the measures of any possible linear transform; the second component accounts for the maximal amount of residual variation; and so on. The principal components are constructed so that they represent transformed scores on dimensions that are orthogonal. Through the use of principal components analysis, it is possible to have a set of highly related software attributes mapped into a small number of uncorrelated attribute domains. This definitively solves the problem of multi-collinearity in subsequent regression analysis. There are many software metrics in the literature, but principal component analysis reveals that there are few distinct sources of variation, i.e. dimensions, in this set of metrics. It would appear perfectly reasonable to characterize the measurable attributes of a program with a simple function of a small number of orthogonal metrics each of which represents a distinct software attribute domain.

Munson, John C.↗