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 415 records · Page 23

The effect of the dynamic wet troposphere on VLBI measurements

Calculations using a statistical model of water vapor fluctuations yield the effect of the dynamic wet troposphere on Very Long Baseline Interferometry (VLBI) measurements. The statistical model arises from two primary assumptions: (1) the spatial structure of refractivity fluctuations can be closely approximated by elementary (Kolmogorov) turbulence theory, and (2) temporal fluctuations are caused by spatial patterns which are moved over a site by the wind. The consequences of these assumptions are outlined for the VLBI delay and delay rate observables. For example, wet troposphere induced rms delays for Deep Space Network (DSN) VLBI at 20-deg elevation are about 3 cm of delay per observation, which is smaller, on the average, than other known error sources in the current DSN VLBI data set. At 20-deg elevation for 200-s time intervals, water vapor induces approximately 1.5 x 10 to the minus 13th power s/s in the Allan standard deviation of interferometric delay, which is a measure of the delay rate observable error. In contrast to the delay error, the delay rate measurement error is dominated by water vapor fluctuations. Water vapor induced VLBI parameter errors and correlations are calculated. For the DSN, baseline length parameter errors due to water vapor fluctuations are in the range of 3 to 5 cm. The above physical assumptions also lead to a method for including the water vapor fluctuations in the parameter estimation procedure, which is used to extract baseline and source information from the VLBI observables.

Treuhaft, R. N.↗

Spectral characteristics and the extent of paleosols of the Palouse formation

Thematic mapping data was analyzed and verified by comparison to previously gathered transect samples and to aerial photographs. A bare-soil field with exposed paleosols characterized by slight enrichment of iron was investigated. Spectral relationships were first investigated statistically by creating a data set with DN values spatially matched as nearly as possible to field sample points. Chemical data for each point included organic carbon, free iron oxide, and amorphous iron content. The chemical data, DN values, and various band ratios were examined with the program package Statistix in order to find the combinations of reflectance data most likely to show a relationship which would dependably separate the exposed paleosols from the other soils. Cluster analysis and Fastclas classification procedures were applied to the most promising of the band ratio combinations.

Frazier, B. E.↗

Contextual classification on the massively parallel processor

Classifiers are often used to produce land cover maps from multispectral Earth observation imagery. Conventionally, these classifiers have been designed to exploit the spectral information contained in the imagery. Very few classifiers exploit the spatial information content of the imagery, and the few that do rarely exploit spatial information content in conjunction with spectral and/or temporal information. A contextual classifier that exploits spatial and spectral information in combination through a general statistical approach was studied. Early test results obtained from an implementation of the classifier on a VAX-11/780 minicomputer were encouraging, but they are of limited meaning because they were produced from small data sets. An implementation of the contextual classifier is presented on the Massively Parallel Processor (MPP) at Goddard that for the first time makes feasible the testing of the classifier on large data sets.

Tilton, James C.↗

Time-cumulated visible and infrared histograms used as descriptor of cloud cover

To study the statistical behavior of clouds for different climate regimes, the spatial and temporal stability of VIS-IR bidimensional histograms is tested. Also, the effect of data sampling and averaging on the histogram shapes is considered; in particular the sampling strategy used by the International Satellite Cloud Climatology Project is tested.

Seze, G.↗

Asymptotic modal analysis and statistical energy analysis

The sound field of a structural-acoustic enclosure was subject to experimental analysis and theoretical description in order to develop an efficient and accurate method for predicting sound pressure levels in enclosures such as aircraft fuselages. Asymptotic Modal Analysis (AMA) is the method under investigation. AMA is derived from classical modal analysis (CMA) by considering the asymptotic limit of the sound pressure level as the number of acoustic and/or structural modes approaches infinity. Using AMA, results identical to those of Statistical Energy Analysis (SEA) were obtained for the spatially-averaged sound pressure levels in the interior. AMA is systematically derived from CMA and therefore the degree of generality of the end result can be adjusted through the choice of appropriate simplifying assumptions. For example, AMA can be used to obtain local sound pressure levels at particular points inside the enclosure, or to include the effects of varying the size and/or location of the sound source. AMA theoretical results were compared with CMA theory and also with experiment for the case where the structural-acoustic enclosure is a rectangular cavity with part of one wall flexible and vibrating, while the rest of the cavity is rigid.

Dowell, Earl H.↗

Pervasive variability in the quiet solar transition region

Extreme UV spectroheliograms from an experimental time series are employed to investigate the nature of the quiet solar transition region with respect to temporal variability. A statistical treatment is developed to analyze the fraction of spatial elements that yield intensity variations significantly higher than the signal noise. The EUV intensity in every spatial resolution varies on a time scale of minutes by about 10-30 percent, and the fractional amplitude of the temporal variations is found to be nearly independent of mean intensity. The paper concludes that the quiet solar transition region is probably generated and modulated by small-scale magnetic activity, since the temporal variability in the transition region is spatially uniform and rather pervasive.

Rabin, Douglas↗

Investigating the Performance of the Jason-2/OSTM Radar Altimeter Over Lakes and Reservoirs

Many inland water investigations utilize archival and near-real time radar altimetry data to enable observation of the variation in surface water level. A multi-altimeter approach allows a more global outlook with improved spatial resolution, and combined long-term observations improve statistical analyses. Central to all programs is a performance assessment of each instrument. Here, we focus on data quantity and quality pertaining to the Poseidon-3 radar altimeter onboard the Jason-2/OSTM satellite.Utilizing an interim data set (IGDR), studies show that the new on-board DIODE/median and DIODE/DEM tracking modes are performing well, acquiring and maintaining the majority of lake and reservoir surfaces in varying terrains. The 20-Hz along-track resolution of the data, and particularly the availability of the range output from the ice-retracker algorithm, also improves the number of valid height measurements. Based on test-case lakes and reservoirs, output from the ice-retracker algorithm is also seen to have a clear advantage over the ocean-retracker having better height stability across calm and icy surfaces, a greater ability to gain coastline waters, and less sensitivity to loss of water surface when there is island contamination in the radar echo. Such on-board tracking and postprocessing retracking enables the lake waters to be quickly gained after coastline crossing. Values can range from <0.1 s to 2.5 s, but the majority of measurements are obtained in less than 0.4 s or <2.3 km from the coast. Validation exercises reveal that targets of 150 km2 surface area and 0.8 km width are able to be monitored offering greater potential to acquire lakes in the 100 C300 km2 size-category. Time series of height variations are also found to be accurate to 3 to 33 cm rms depending on target size and the presence of winter ice. These findings are an improvement over the IGDR/GDR results from the predecessor Jason-1 and TOPEX/Poseidon missions and can satisfy the accuracy requirements of both the science-related and operational lake study programs.

Birkett, C. M.↗

The role of high resolution observations in determining energy release and transport processes

With present observations from the Solar Maximum Mission, we are seeing a very selective sample of spatially large flares in the soft part (10 to 50 keV) of the hard X-ray spectrum. The spatial resolution is at best 5600 km with a corresponding time resolution of 4.5 s for adequate count statistics. This resolution gives rise to the following problems: We cannot resolve the minor radii of the loops involved or tell where and how the energy release occurs. The manner in which loops interact and the relationship between the soft and hard or approx 00 keV) parts of hard X-rays remains elusive. We cannot see how energy propagates in most cases. Thus it is desirable to determine the minimum increase in spatial and temporal resolution required to solve these problems. Spatially we need to resolve the minor radius of a small loop which is about 800 km or 1 arc s. Upper limits to observed speeds of conduction fronts and shocks are approx km/s with theoretical limits running about a factor of 2 higher. Thus, a compatible minimum time resolution is in the range of 0.2 to 0.4 s. With these spatial and temporal resolutions, sufficient count statistics are required to go up to approx 120 keV with a sufficient number of energy bands to obtain spectra.

Smith, D. F.↗

Intraseasonal oscillations in the global atmosphere. II - Southern Hemisphere

A systematic examination of oscillatory modes in the global atmosphere is completed by studying 12 years of 500 mb geopotential heights in the Southern Hemisphere. The data were band-pass filtered to focus on intraseasonal (IS) phenomena, and spatial empirical orthogonal functions (EOFs) were obtained. The leading principal components were subjected to singular spectrum analysis in order to identify nonlinear IS oscillations with high statistical confidence. In the Southern Hemisphere, the dominant mode has a period of 23 days, with spatial patterns carried by the second and third winter EOF of the IS band. It has a zonal wavenumber-four structure. The 40-day mode is second, and is dominated by wavenumbers three and four, while a 16-day mode is too weak to separate its spatial behavior from the previous two. The IS dynamics in the Southern Hemisphere is more complex and dominated by shorter wavenumbers than the Northern Hemisphere. No statistically significant correlations between the Southern Hemisphere and the tropics or the Northern Hemisphere are apparent in the IS band.

Ghil, Michael↗

Fine-scale motion in the central Orion nebula

The inner region of M42 has been mapped in the O III N1 forbidden line at an unprecedented combination of velocity and spatial resolution and photometric accuracy. Three separate velocity systems along the line of sight are identified. Differential mass flow is present in two of these systems. All three systems show random fluctuations in velocity well above their measurement errors. The point-by-point velocities are used to determine the diagnostic statistical parameter, the structure function. The random motions show correlation with spatial separation, as expected for turbulence; but comparison with the predictions of classical theory for nebular turbulence shows poor agreement. If Kolmogorov (1941, 1942) theory applies, then the energy must be input at many scales.

Castaneda, Hector O.↗

Statistical Quality Control of Moisture Data in GEOS DAS

A new statistical quality control algorithm was recently implemented in the Goddard Earth Observing System Data Assimilation System (GEOS DAS). The final step in the algorithm consists of an adaptive buddy check that either accepts or rejects outlier observations based on a local statistical analysis of nearby data. A basic assumption in any such test is that the observed field is spatially coherent, in the sense that nearby data can be expected to confirm each other. However, the buddy check resulted in excessive rejection of moisture data, especially during the Northern Hemisphere summer. The analysis moisture variable in GEOS DAS is water vapor mixing ratio. Observational evidence shows that the distribution of mixing ratio errors is far from normal. Furthermore, spatial correlations among mixing ratio errors are highly anisotropic and difficult to identify. Both factors contribute to the poor performance of the statistical quality control algorithm. To alleviate the problem, we applied the buddy check to relative humidity data instead. This variable explicitly depends on temperature and therefore exhibits a much greater spatial coherence. As a result, reject rates of moisture data are much more reasonable and homogeneous in time and space.

Dee, D. P.↗

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗

Predicting weather impacts on corn production in a data-limited region using a transfer learning approach

The stability of food supply and prices may depend more on annual changes in yields from year-to-year variability in weather than on longer-term average changes from changing climatic conditions. However, the absence of high-quality data on crop yields at fine spatial resolutions in many regions of the world makes it challenging to statistically model their response to interannual variability in weather patterns. Therefore, there is a need for empirical methods that can project annual crop yield changes even in limited data regions. Here, we propose a transfer learning algorithm that uses high spatial resolution data from one region to project yields in another region with more limited data. The goal of our work is to understand what data types can be beneficial for transferring learning from a source region to a very different target region with more limited data. We utilize Long Short-Term Memory to develop a transfer learning model that is trained on historical county-level corn yield in the United States and predicts district-level corn yield variations in India. Even using smaller amounts of data in India, simulating a data-scarce region, we achieve an average root mean square error of 0.48 bu acre−1 in predicting interannual yield variations. Using Shapley values to interpret results, we explore the contribution of the different weather parameters to interannual yield variability and find a larger influence of precipitation-related variables. Our study demonstrates the usefulness of this method for transferring models of weather impacts on crop yields trained on a data-rich country to one with more limited data. It suggests the potential of applying the transfer learning model to mitigate the need for extensive raw data globally.

Vishwakarma, Srishti [ORNL] (ORCID:000000031674419↗

Path Profiles of Cn2 Derived from Radiometer Temperature Measurements and Geometrical Ray Tracing

Atmospheric turbulence has significant impairments on the operation of Free-Space Optical (FSO) communication systems, in particular temporal and spatial intensity fluctuations at the receiving aperture resulting in power surges and fades, changes in angle of arrival, spatial coherence degradation, etc. The refractive index structure parameter Cn2 is a statistical measure of the strength of turbulence in the atmosphere and is highly dependent upon vertical height. Therefore to understand atmospheric turbulence effects on vertical FSO communication links such as space-to-ground links, it is necessary to specify Cn2 profiles along the atmospheric propagation path. To avoid the limitations on the applicability of classical approaches, propagation simulation through geometrical ray tracing is applied. This is achieved by considering the atmosphere along the optical propagation path as a spatial distribution of spherical bubbles with varying relative refractive index deviations representing turbulent eddies. The relative deviations of the refractive index are statistically determined from altitude-dependent and time-varying temperature fluctuations, as measured by a microwave profiling radiometer. For each representative atmosphere ray paths are analyzed using geometrical optics, which is particularly advantageous in situations of strong turbulence where there is severe wavefront distortion and discontinuity. The refractive index structure parameter is then determined as a function of height and time.

Atmospheric turbulence↗

Estimated Accuracy of Three Common Trajectory Statistical Methods

Three well-known trajectory statistical methods (TSMs), namely concentration field (CF), concentration weighted trajectory (CWT), and potential source contribution function (PSCF) methods were tested using known sources and artificially generated data sets to determine the ability of TSMs to reproduce spatial distribution of the sources. In the works by other authors, the accuracy of the trajectory statistical methods was estimated for particular species and at specified receptor locations. We have obtained a more general statistical estimation of the accuracy of source reconstruction and have found optimum conditions to reconstruct source distributions of atmospheric trace substances. Only virtual pollutants of the primary type were considered. In real world experiments, TSMs are intended for application to a priori unknown sources. Therefore, the accuracy of TSMs has to be tested with all possible spatial distributions of sources. An ensemble of geographical distributions of virtual sources was generated. Spearman s rank order correlation coefficient between spatial distributions of the known virtual and the reconstructed sources was taken to be a quantitative measure of the accuracy. Statistical estimates of the mean correlation coefficient and a range of the most probable values of correlation coefficients were obtained. All the TSMs that were considered here showed similar close results. The maximum of the ratio of the mean correlation to the width of the correlation interval containing the most probable correlation values determines the optimum conditions for reconstruction. An optimal geographical domain roughly coincides with the area supplying most of the substance to the receptor. The optimal domain s size is dependent on the substance decay time. Under optimum reconstruction conditions, the mean correlation coefficients can reach 0.70 0.75. The boundaries of the interval with the most probable correlation values are 0.6 0.9 for the decay time of 240 h and 0.5 0.95 for the decay time of 12 h. The best results of source reconstruction can be expected for the trace substances with a decay time on the order of several days. Although the methods considered in this paper do not guarantee high accuracy they are computationally simple and fast. Using the TSMs in optimum conditions and taking into account the range of uncertainties, one can obtain a first hint on potential source areas.

Kabashnikov, Vitaliy P.↗

A stochastic precipitation disaggregation scheme for GCM applications

In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Gao, Xiaogang↗

Potential Multi-Component Structure of the Debris Disk Around HIP 17439 Revealed by Herschel DUNES

Context. The dust observed in debris disks is produced through collisions of larger bodies left over from the planet/planetesimal formation process. Spatially resolving these disks permits to constrain their architecture and thus that of the underlying planetary/planetesimal system. Aims. Our Herschel open time key program DUNES aims at detecting and characterizing debris disks around nearby, sun-like stars. In addition to the statistical analysis of the data, the detailed study of single objects through spatially resolving the disk and detailed modeling of the data is a main goal of the project. Methods. We obtained the first observations spatially resolving the debris disk around the sun-like star HIP 17439 (HD 23484) using the instruments PACS and SPIRE on board the Herschel Space Observatory. Simultaneous multi-wavelength modeling of these data together with ancillary data from the literature is presented. Results. A standard single component disk model fails to reproduce the major axis radial profiles at 70 μm, 100 μm, and 160 μm simultaneously. Moreover, the best-fit parameters derived from such a model suggest a very broad disk extending from few au up to few hundreds of au from the star with a nearly constant surface density which seems physically unlikely. However, the constraints from both the data and our limited theoretical investigation are not strong enough to completely rule out this model. An alternative, more plausible, and better fitting model of the system consists of two rings of dust at approx. 30 au and 90 au, respectively, while the constraints on the parameters of this model are weak due to its complexity and intrinsic degeneracies. Conclusions. The disk is probably composed of at least two components with different spatial locations (but not necessarily detached), while a single, broad disk is possible, but less likely. The two spatially well-separated rings of dust in our best-fit model suggest the presence of at least one high mass planet or several low-mass planets clearing the region between the two rings from planetesimals and dust.

sun-like stars↗

FilDReaMS: II. Application to the analysis of the relative orientations between filaments and the magnetic field in four Herschel fields

Context. Both simulations and observations of the interstellar medium show that the study of the relative orientations between filamentary structures and the magnetic field can bring new insight into the role played by magnetic fields in the formation and evolution of filaments and in the process of star formation. Aims. We provide a first application of FilDReaMS, the new method presented in the companion paper to detect and analyze filaments in a given image. The method relies on a template that has the shape of a rectangular bar with variable width. Our goal is to investigate the relative orientations between the detected filaments and the magnetic field. Methods. We apply FilDReaMS to a small sample of four Herschel fields (G210, G300, G82, G202) characterized by different Galactic environments and different evolutionary stages. First, we look for the most prevalent bar widths, and we examine the networks formed by filaments of different bar widths as well as their hierarchical organization. Second, we compare the filament orientations to the magnetic field orientation inferred from Planck polarization data and, for the first time, we study the statistics of the relative orientation angle as functions of both spatial scale and H2 column density. Results. We find preferential relative orientations in the four Herschel fields: small filaments with low column densities tend to be slightly more parallel than perpendicular to the magnetic field; in contrast, large filaments, which all have higher column densities, are oriented nearly perpendicular (or, in the case of G202, more nearly parallel) to the magnetic field. In the two nearby fields (G210 and G300), we observe a transition from mostly parallel to mostly perpendicular relative orientations at an H 2 column density ≃ 1.1 × 10 21 cm -2 and 1.4 × 10 21 cm -2 , respectively, consistent with the results of previous studies. Conclusions. Our results confirm the existence of a coupling between magnetic fields at cloud scales and filaments at smaller scale. They also illustrate the potential of combining Herschel and Planck observations, and they call for further statistical analyses with our dedicated method.

79 ASTRONOMY AND ASTROPHYSICS↗