On the a priori information in sequential estimation problems.
Error effect in a priori information on sequential estimate variance of linear system states, noting optimal filter synthesis
SEARCH · Search NASA
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.
Error effect in a priori information on sequential estimate variance of linear system states, noting optimal filter synthesis
A beta-plane model of the stratosphere is used to investigate the planetary-wave amplitude vacillations first reported by Holton and Mass (1976). This model differs from theirs in allowing more horizontal modes. For low surface wave amplitudes, a new class of solutions is found which exhibits a stationary, partially reflecting critical line at steady state. The critical line equilibrates at lower altitudes as the wave forcing is increased. Vacillating solutions occur when the steady state critical line occurs near the lower boundary. The maximum wave amplitude and the maximum steady-state wave amplitude found in the model are in the ratio of 2:1, in good agreement with theoretical predictions. The maximum wave amplitude never exceeds 2200 gpm which is quite close to the saturation limit predicted by Schoeberl (1982). An analysis of the statistics of slowly and rapidly vacillating flows shows that both the wave and zonal mean variances are important in determining the time mean, zonal mean dynamics of the upper stratosphere.
The potential of short wave infrared (SWIR) measurements in vegetation discrimination is further substantiated through a discussion of field studies and an examination of the physical bases which cause SWIR measurements to vary with the vegetation type observed. The research reported herein supported the AGRISTARS program objective to incorporate TM measurements in the analysis of agricultural activity. Field measurements on corn and soybeans in Iowa were conducted, and the mean and variance of canopy reflectance were computed for each observation date. The Suits canopy reflectance model was used to evaluate possible explanations of the observed corn/soybeans reflectance patterns /39/. The SWIR measurements were shown to effectively discriminate corn and soybeans on the basis of leaf absorption properties.
A general procedure for error detection in complex systems, called the data block capture and analysis monitoring process, is described and analyzed. It is assumed that, in addition to being exposed to potential external fault sources, a complex system will in general always contain embedded hardware and software fault mechanisms which can cause the system to perform incorrect computations and/or produce incorrect output. Thus, in operation, the system continuously moves back and forth between error and no-error states. These external fault sources or internal fault mechanisms are extremely difficult to detect. The data block capture and analysis monitoring process is concerned with detecting deviations from the normal performance of the system, known as errors, which are symptomatic of fault conditions. The process consists of repeatedly recording a fixed amount of data from a set of predetermined observation lines of the system being monitored (i.e., capturing a block of data) and then analyzing the captured block in an attempt to determine whether the system is functioning correctly. The performances of linear, quadratic, and logarithmic data analysis algorithms are rigorously characterized in terms of the probability of correctly detecting an error, the expectation and variance of the number of false alarms per error, and the expectation and variance of the latency in detection of errors. Insight into the nature of the general problem of error detection is obtained.
We use Monte Carlo analysis to estimate the uncertainty in predictions of total O3 trends between 1979 and 1995 made by the Goddard Space Flight Center (GSFC) two-dimensional (2D) model of stratospheric photochemistry and dynamics. The uncertainty is caused by gas-phase chemical reaction rates, photolysis coefficients, and heterogeneous reaction parameters which are model inputs. The uncertainty represents a lower bound to the total model uncertainty assuming the input parameter uncertainties are characterized correctly. Each of the Monte Carlo runs was initialized in 1970 and integrated for 26 model years through the end of 1995. This was repeated 419 times using input parameter sets generated by Latin Hypercube Sampling. The standard deviation (a) of the Monte Carlo ensemble of total 03 trend predictions is used to quantify the model uncertainty. The 34% difference between the model trend in globally and annually averaged total O3 using nominal inputs and atmospheric trends calculated from Nimbus 7 and Meteor 3 total ozone mapping spectrometer (TOMS) version 7 data is less than the 46% calculated 1 (sigma), model uncertainty, so there is no significant difference between the modeled and observed trends. In the northern hemisphere midlatitude spring the modeled and observed total 03 trends differ by more than 1(sigma) but less than 2(sigma), which we refer to as marginal significance. We perform a multiple linear regression analysis of the runs which suggests that only a few of the model reactions contribute significantly to the variance in the model predictions. The lack of significance in these comparisons suggests that they are of questionable use as guides for continuing model development. Large model/measurement differences which are many multiples of the input parameter uncertainty are seen in the meridional gradients of the trend and the peak-to-peak variations in the trends over an annual cycle. These discrepancies unambiguously indicate model formulation problems and provide a measure of model performance which can be used in attempts to improve such models.
In the Project for Intercomparison of Land-Surface Parameterization Schemes phase 2a experiment, meteorological data for the year 1987 from Cabauw, the Netherlands, were used as inputs to 23 land-surface flux schemes designed for use in climate and weather models. Schemes were evaluated by comparing their outputs with long-term measurements of surface sensible heat fluxes into the atmosphere and the ground, and of upward longwave radiation and total net radiative fluxes, and also comparing them with latent heat fluxes derived from a surface energy balance. Tuning of schemes by use of the observed flux data was not permitted. On an annual basis, the predicted surface radiative temperature exhibits a range of 2 K across schemes, consistent with the range of about 10 W/m in predicted surface net radiation. Most modeled values of monthly net radiation differ from the observations by less than the estimated maximum monthly observational error (+/- 10 W/sq m). However, modeled radiative surface temperature appears to have a systematic positive bias in most schemes; this might be explained by an error in assumed emissivity and by models' neglect of canopy thermal heterogeneity. Annual means of sensible and latent heat fluxes, into which net radiation is partitioned, have ranges across schemes of 30 W/sq m and 25 W/sq m, respectively. Annual totals of evapotranspiration and runoff, into which the precipitation is partitioned, both have ranges of 315 mm. These ranges in annual heat and water fluxes were approximately halved upon exclusion of the three schemes that have no stomatal resistance under non-water-stressed conditions. Many schemes tend to underestimate latent heat flux and overestimate sensible heat flux in summer, with a reverse tendency in winter. For six schemes, root-mean-square deviations of predictions from monthly observations are less than the estimated upper bounds on observation errors (5 W/m for sensible beat flux and 10 W/m for latent heat flux). Actual runoff at the site is believed to be dominated by vertical drainage to ground water, but several schemes produced significant amounts of runoff as overland flow or interflow. There is a range across schemes of 184 mm (40% of total pore volume) in the simulated annual mean root-zone soil moisture. Unfortunately, no measurements of soil moisture were available for model evaluation. A theoretical analysis suggested that differences in boundary conditions used in various schemes are not sufficient to explain the large variance in soil moisture. However, many of the extreme values of soil moisture could be explained in terms of the particulars of experimental setup or excessive evapotranspiration.
The traditional method of applying a digital optimal filter to measure X-ray pulses from transition-edge sensor (TES) devices does not achieve the best energy resolution when the signals have a highly non-linear response to energy, or the noise is non-stationary during the pulse. We present an implementation of a method to analyze X-ray data from TESs, which is based upon principal component analysis (PCA). Our method separates the X-ray signal pulse into orthogonal components that have the largest variance. We typically recover pulse height, arrival time, differences in pulse shape, and the variation of pulse height with detector temperature. These components can then be combined to form a representation of pulse energy. An added value of this method is that by reporting information on more descriptive parameters (as opposed to a single number representing energy), we generate a much more complete picture of the pulse received. Here we report on progress in developing this technique for future implementation on X-ray telescopes. We used an 55Fe source to characterize Mo/Au TESs. On the same dataset, the PCA method recovers a spectral resolution that is better by a factor of two than achievable with digital optimal filters.
Observability of satellite state vectors, using only SST tracking data was investigated by covariance analysis under a variety of satellite and station configurations. These results indicate very precarious observability in most short arc cases. The consequences of this are large variances on many state components, such as the downrange component of the relay satellite position. To illustrate the impact of observability problems, an example is given of two distinct satellite orbit pairs generating essentially the same data arc. The physical bases for unobservability are outlined and related to proposed TDRSS configurations. Results are relevant to any mission depending upon TDRSS to determine satellite state. The required mathematical analysis and the software used is described.
This report explores the magnitude of the potential buildout of new solar development on Bureau of Land Management (BLM) land in Arizona based on market conditions in 2023. The BLM currently has 18 solar energy zones (SEZs), including three in Arizona, where land has undergone preliminary reviews and is suitable for expedited environmental review and permitting for solar energy projects. The BLM also administers variance areas that are outside SEZs and may consider right-of-way (ROW) applications for utility-scale solar development on such variance lands. Within variance areas in Arizona, the BLM has further identified Renewable Energy Development Area (REDA) lands that are likely well suited to solar or wind development. The purpose of the analysis reported here is to explore the anticipated extent of new solar development on BLM land in the state by 2035.
The frequency distribution of overshoots in a stationary Gaussian stochastic process is analyzed. The primary processes involved in this analysis are computer simulation and statistical estimation. Computer simulation is used to simulate stationary Gaussian stochastic processes that have selected autocorrelation functions. An analysis of the simulation results reveals a frequency distribution for overshoots with a functional dependence on the mean and variance of the process. Statistical estimation is then used to estimate the mean and variance of a process. It is shown that for an autocorrelation function, the mean and the variance for the number of overshoots, a frequency distribution for overshoots can be estimated.
The carrier-phase recovery performance of a Costas loop for a nonlinear satellite channel with uplink and downlink noise is analyzed, and the extension of these results to the case where uplink radio-frequency interference (RFI) is present is considered. The signal format here is binary phase-shift-keyed (BPSK) or unbalanced quaternary phase-shift-keyed (UQPSK), the signals being either the non-return-to-zero (NRZ) or biphase (Manchester) type. With UQPSK, at least one data stream must be PN-spread. The channel comprises additive white Gaussian noise (WGN), a wideband filter, a bandpass nonlinearity and more additive WGN. The RFI being considered is pulsed continuous-wave, pulsed WGN, or a combination of the two. It is noted that the Tracking and Data Relay Satellite System (TDRSS) will be subjected to RFI from certain locations, especially in the S-band of frequencies. The characteristics of this link are generalized to form the assumptions underlying this analysis. The analytical results given include the S-curve of the equivalent loop and the phase-error variance of the linearized loop.
Backscattered ultraviolet (BUV) data from the Nimbus-4 spacecraft for the period 1970-1977 have been processed to a refined level. The seasonal and interannual variations of total ozone are examined on a global scale, using daily zonal means of 10 deg latitude bands and a time-latitude cross section. A harmonic analysis was performed on the daily zonal means and the amplitude, days of peak ozone values, and percentage of variance were computed for the annual, semiannual and higer harmonics for several years and each year. A clear quasi-biennial oscillation (QBO) was revealed from the tropics to midlatitudes after removing the mean annual wave. Asymmetries in the annual wave in the two hemispheres found earlier for the period 1970-1972 persist through the entire observation period. Interannual variations appear to be the result of the QBO of ozone from low to midlatitudes. Asymmetries in the QBO amplitude phase, and period were also detected in the two hemispheres.
Recently, Feynman and Ruzmaikin (1994) showed that IMF parameters for the 1973 to 1990 period are not log-normally distributed as previously suggested by Burlaga and King (1979) for the data obtained over a shorter time period (1963-75). They studied the first four moments, namely: mean, variance, skewness, and kurtosis. For a Gaussian distribution, moments higher than the variance should vanish. In particular, Feynman and Ruzmaikin obtained very high values of kurtosis during some periods of their analysis. We note that the coverage for IMF parameters is very uneven for the period analyzed by them, ranging from less than 40% to greater than 80%. So a question arises as to whether the amount of coverage is a factor in their analysis. We decided to test this for the B(sub z) component of IMF, since it is an effective geoactive parameter for short term disturbances. Like them, we used 1-hour averaged data available on the Omnitape. We studied the scatter plots of the annual mean values of B(sub z)(nT) and its kurtosis versus the percent coverage for the year. We obtain a correlation coefficient of 0.48 and 0.42 respectively for the 1973-90 period. The probability for a chance occurrence of these correlation coefficients for 18 pair of points is less than 8%. As a rough measure of skewness, we determined the percent asymmetry between the areas of the histograms representing the distributions of the positive and the negative values of B(sub z) and studied its correlation with the coverage for the year. This analysis yields a correlation coefficient of 0.41 When we extended the analysis for the whole period for which IMF data are available (1963-93) the corresponding correlation coefficients are 0.59, 0.14, and 0.42. Our findings will be presented and discussed
The annual cycles of upward and downward longwave fluxes at the Earth s surface are investigated by use of the NASA/GEWEX Surface Radiation Budget Data Set. Because of the immense difference between the heat capacity of land and ocean, the surface of Earth is partitioned into these two categories. Principal component analysis is used to quantify the annual cycles. Over land, the first principal component describes over 95% of the variance of the annual cycle of the upward and downward longwave fluxes. Over ocean the first term describes more than 87% of these annual cycles. Empirical orthogonal functions show the corresponding geographical distributions of these cycles. Phase plane diagrams of the annual cycles of upward longwave fluxes as a function of net shortwave flux show the thermal inertia of land and ocean.
A near-global chart of surface elevations associated with the stationary M2 internal tide is empirically constructed from multi-mission satellite altimeter data. An advantage of a strictly empirical mapping approach is that results are independent of assumptions about ocean wave dynamics and, in fact, can be used to test such assumptions. A disadvantage is that present-day altimeter coverage is only marginally adequate to support mapping such short-wavelength features. Moreover, predominantly north-south ground-track orientations and contamination from nontidal oceanographic variability can lead to deficiencies in mapped tides. Independent data from Cryosphere Satellite-2 (CryoSat-2) and other altimeters are used to test the solutions and show positive reduction in variance except in regions of large mesoscale variability. The tidal fields are subjected to two-dimensional wavenumber spectral analysis, which allows for the construction of an empirical map of modal wavelengths. Mode-1 wavelengths show good agreement with theoretical wavelengths calculated from the ocean's mean stratification, with a few localized exceptions (e.g., Tasman Sea). Mode-2 waves are detectable in much of the ocean, with wavelengths in reasonable agreement with theoretical expectations, but their spectral signatures grow too weak to map in some regions.
Compilation describes characteristics of truncated and censored samples, and presents six illustrations of practical use of tables in computing mean and variance estimates for normal distribution using selected samples.
The author has identified the following significant results. A quantitative measure of the registration processor accuracy in terms of the variance of the registration error was derived. With the appropriate assumptions, the variance was shown to be inversely proportional to the square of the effective bandwidth times the signal to noise ratio. The final expressions were presented to emphasize both the form and simplicity of their representation. In the situation where relative spatial distortions exist between images to be registered, expressions were derived for estimating the loss in output signal to noise ratio due to these spatial distortions. These results are in terms of a reduction factor.
Rainfall data collected by radar in the vicinity of Darwin, Australia, have been analyzed in terms of their mean, variance, autocorrelation of area-averaged rain rate, and diurnal variation. It is found that, when compared with the well-studied GATE (Global Atmospheric Research Program Atlantic Tropical Experiment) data, Darwin rainfall has larger coefficient of variation (CV), faster reduction of CV with increasing area size, weaker temporal correlation, and a strong diurnal cycle and intermittence. The coefficient of variation for Darwin rainfall has larger magnitude and exhibits larger spatial variability over the sea portion than over the land portion within the area of radar coverage. Stationary, and nonstationary models have been used to study the sampling errors associated with space-based rainfall measurement. The nonstationary model shows that the sampling error is sensitive to the starting sampling time for some sampling frequencies, due to the diurnal cycle of rain, but not for others. Sampling experiments using data also show such sensitivity. When the errors are averaged over starting time, the results of the experiments and the stationary and nonstationary models match each other very closely. In the small areas for which data are available for I>oth Darwin and GATE, the sampling error is expected to be larger for Darwin due to its larger CV.