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.

SourceResearch connectionFrom yearThrough year

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

At least 1,207 records · Page 67

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Nonparametric statistical modeling of binary star separations

We develop a comprehensive statistical model for the distribution of observed separations in binary star systems, in terms of distributions of orbital elements, projection effects, and distances to systems. We use this model to derive several diagnostics for estimating the completeness of imaging searches for stellar companions, and the underlying stellar multiplicities. In application to recent imaging searches for low-luminosity companions to nearby M dwarf stars, and for companions to young stars in nearby star-forming regions, our analyses reveal substantial uncertainty in estimates of stellar multiplicity. For binary stars with late-type dwarf companions, semimajor axes appear to be distributed approximately as a(exp -1) for values ranging from about one to several thousand astronomical units. About one-quarter of the companions to field F and G dwarf stars have semimajor axes less than 1 AU, and about 15% lie beyond 1000 AU. The geometric efficiency (fraction of companions imaged onto the detector) of imaging searches is nearly independent of distances to program stars and orbital eccentricities, and varies only slowly with detector spatial limitations.

Heacox, William D.↗

NASA NTRS1993

A study of the influence of satellite data on GLA analyses over the Atlantic Ocean during a period of blocking anticyclone development

This paper presents a statistical diagnosis of the differences between two Goddard Laboratory for Atmospheres analysis representations of temperature and wind speed, one containing satellite data influences (SAT analyses) and one without (NOSAT analyses), for a case of oceanic blocking anticyclone development. Results, obtained using area means and standard deviations of the two fields and correlation coefficients and root-mean-square differences between the two fields, indicate that the inclusion of satellite data can have significant impact on ocean-domain analyses. This is especially evident if one examines higher-order fields that represent gradients in the basic variable fields or contain covariances between two different variables (e.g., advections). More specifically, the inclusion of satellite data resulted in a cold (warm) temperature bias at low (high) levels and weaker temperature gradients, stronger winds, weaker vertical wind shears, and a warm-air advection bias at most levels. The SAT analyses also exhibited larger standard deviations than the NOSAT for wind speed, relative vorticity, vorticity advection, temperature advection, and 500-mb height tendencies, suggesting that for this case the spatial variability of the circulation features were enhanced by the inclusion of satellite data. Even so, since height tendencies forced by vorticity and temperature advection were less sensitive to the addition of satellite data than were the advection quantities themselves, the dynamics of the blocking system was apparently less influenced by the presence of satellite data than was the structure of the system.

Lamberty, Gregory L.↗

NASA NTRS2009

Hybrid NN/SVM Computational System for Optimizing Designs

A computational method and system based on a hybrid of an artificial neural network (NN) and a support vector machine (SVM) (see figure) has been conceived as a means of maximizing or minimizing an objective function, optionally subject to one or more constraints. Such maximization or minimization could be performed, for example, to optimize solve a data-regression or data-classification problem or to optimize a design associated with a response function. A response function can be considered as a subset of a response surface, which is a surface in a vector space of design and performance parameters. A typical example of a design problem that the method and system can be used to solve is that of an airfoil, for which a response function could be the spatial distribution of pressure over the airfoil. In this example, the response surface would describe the pressure distribution as a function of the operating conditions and the geometric parameters of the airfoil. The use of NNs to analyze physical objects in order to optimize their responses under specified physical conditions is well known. NN analysis is suitable for multidimensional interpolation of data that lack structure and enables the representation and optimization of a succession of numerical solutions of increasing complexity or increasing fidelity to the real world. NN analysis is especially useful in helping to satisfy multiple design objectives. Feedforward NNs can be used to make estimates based on nonlinear mathematical models. One difficulty associated with use of a feedforward NN arises from the need for nonlinear optimization to determine connection weights among input, intermediate, and output variables. It can be very expensive to train an NN in cases in which it is necessary to model large amounts of information. Less widely known (in comparison with NNs) are support vector machines (SVMs), which were originally applied in statistical learning theory. In terms that are necessarily oversimplified to fit the scope of this article, an SVM can be characterized as an algorithm that (1) effects a nonlinear mapping of input vectors into a higher-dimensional feature space and (2) involves a dual formulation of governing equations and constraints. One advantageous feature of the SVM approach is that an objective function (which one seeks to minimize to obtain coefficients that define an SVM mathematical model) is convex, so that unlike in the cases of many NN models, any local minimum of an SVM model is also a global minimum.

Rai, Man Mohan↗

NASA NTRS2020

The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Aliaksandr Siniuk↗

NASA NTRS2002

A quantitative comparison of simultaneous BOLD fMRI and NIRS recordings during functional brain activation

Near-infrared spectroscopy (NIRS) has been used to noninvasively monitor adult human brain function in a wide variety of tasks. While rough spatial correspondences with maps generated from functional magnetic resonance imaging (fMRI) have been found in such experiments, the amplitude correspondences between the two recording modalities have not been fully characterized. To do so, we simultaneously acquired NIRS and blood-oxygenation level-dependent (BOLD) fMRI data and compared Delta(1/BOLD) (approximately R(2)(*)) to changes in oxyhemoglobin, deoxyhemoglobin, and total hemoglobin concentrations derived from the NIRS data from subjects performing a simple motor task. We expected the correlation with deoxyhemoglobin to be strongest, due to the causal relation between changes in deoxyhemoglobin concentrations and BOLD signal. Instead we found highly variable correlations, suggesting the need to account for individual subject differences in our NIRS calculations. We argue that the variability resulted from systematic errors associated with each of the signals, including: (1) partial volume errors due to focal concentration changes, (2) wavelength dependence of this partial volume effect, (3) tissue model errors, and (4) possible spatial incongruence between oxy- and deoxyhemoglobin concentration changes. After such effects were accounted for, strong correlations were found between fMRI changes and all optical measures, with oxyhemoglobin providing the strongest correlation. Importantly, this finding held even when including scalp, skull, and inactive brain tissue in the average BOLD signal. This may reflect, at least in part, the superior contrast-to-noise ratio for oxyhemoglobin relative to deoxyhemoglobin (from optical measurements), rather than physiology related to BOLD signal interpretation.

Non-NASA Center↗

NASA NTRS1993

Mapping of the Ronda peridotite massif (Spain) from AVIRIS spectro-imaging survey: A first attempt

In both AVIRIS and ISM data, through the use of mixing models, geological boundaries of the Ronda massif are identified with respect to the surrounding rocks. We can also yield first-order vegetation maps. ISM and AVIRIS instruments give consistent results. On the basis of endmember fraction images, it is then possible to discard areas highly vegetated or not belonging to the peridotite massif. Within the remaining part of the mosaic, spectro-mixing analysis reveals spectral variations in the peridotite massif between the well-exposed areas. Spatially organized units are depicted, related to differences in the relative depth of the absorption band at 1 micron, and it may be due to a different pyroxene content. At this stage, it is worth noting that, although mineralogical variations observed in the rocks are at a sub-pixel scale for the airborne analysis, we see an emerging spatial pattern in the distribution of spectral variations across the massif which might be prevailingly related to mineralogy. Although it is known from fieldwork that the Ronda peridotite massif exhibits mineralogical variations at local scale in the content of pyroxene, and at regional scale in different mineral facies, ranging from garnet-, to spinel- to plagioclase-lherzolites, no attempt has been done yet to produce a synoptic map relating the two scales of analysis. The present work is a first attempt to reach this objective, though a lot more work is still required. In particular, for the purpose of mineralogical interpretation, it is critical to relate the airborne observation to field work and laboratory spectra of Ronda rocks already obtained, with the use of image endmembers and associated reference endmembers. Also, the pretty rough linear mixing model used here is taken as a 'black-box' process which does not necessarily apply correctly to the physical situation at the sub-pixel level. One may think of using the ground-truth observations bearing on the sub-pixel statistical characteristics (texture, structural pattern, surface distribution and vegetation contribution (grass,..)) to produce a more advanced mixing model, physically appropriate to the geologic and environmental contexts.

Pinet, P. C.↗

NASA NTRS2014

Using a Support Vector Machine and a Land Surface Model to Estimate Large-Scale Passive Microwave Temperatures over Snow-Covered Land in North America

A support vector machine (SVM), a machine learning technique developed from statistical learning theory, is employed for the purpose of estimating passive microwave (PMW) brightness temperatures over snow-covered land in North America as observed by the Advanced Microwave Scanning Radiometer (AMSR-E) satellite sensor. The capability of the trained SVM is compared relative to the artificial neural network (ANN) estimates originally presented in [14]. The results suggest the SVM outperforms the ANN at 10.65 GHz, 18.7 GHz, and 36.5 GHz for both vertically and horizontally-polarized PMW radiation. When compared against daily AMSR-E measurements not used during the training procedure and subsequently averaged across the North American domain over the 9-year study period, the root mean squared error in the SVM output is 8 K or less while the anomaly correlation coefficient is 0.7 or greater. When compared relative to the results from the ANN at any of the six frequency and polarization combinations tested, the root mean squared error was reduced by more than 18 percent while the anomaly correlation coefficient was increased by more than 52 percent. Further, the temporal and spatial variability in the modeled brightness temperatures via the SVM more closely agrees with that found in the original AMSR-E measurements. These findings suggest the SVM is a superior alternative to the ANN for eventual use as a measurement operator within a data assimilation framework.

Modeling↗

NASA NTRS2001

Error Estimation of An Ensemble Statistical Seasonal Precipitation Prediction Model

This NASA Technical Memorandum describes an optimal ensemble canonical correlation forecasting model for seasonal precipitation. Each individual forecast is based on the canonical correlation analysis (CCA) in the spectral spaces whose bases are empirical orthogonal functions (EOF). The optimal weights in the ensemble forecasting crucially depend on the mean square error of each individual forecast. An estimate of the mean square error of a CCA prediction is made also using the spectral method. The error is decomposed onto EOFs of the predictand and decreases linearly according to the correlation between the predictor and predictand. Since new CCA scheme is derived for continuous fields of predictor and predictand, an area-factor is automatically included. Thus our model is an improvement of the spectral CCA scheme of Barnett and Preisendorfer. The improvements include (1) the use of area-factor, (2) the estimation of prediction error, and (3) the optimal ensemble of multiple forecasts. The new CCA model is applied to the seasonal forecasting of the United States (US) precipitation field. The predictor is the sea surface temperature (SST). The US Climate Prediction Center's reconstructed SST is used as the predictor's historical data. The US National Center for Environmental Prediction's optimally interpolated precipitation (1951-2000) is used as the predictand's historical data. Our forecast experiments show that the new ensemble canonical correlation scheme renders a reasonable forecasting skill. For example, when using September-October-November SST to predict the next season December-January-February precipitation, the spatial pattern correlation between the observed and predicted are positive in 46 years among the 50 years of experiments. The positive correlations are close to or greater than 0.4 in 29 years, which indicates excellent performance of the forecasting model. The forecasting skill can be further enhanced when several predictors are used.

Shen, Samuel S. P.↗

NASA NTRS2013

GASPS--A Herschel Survey of Gas and Dust in Protoplanetary Disks: Summary and Initial Statistics

We describe a large-scale far-infrared line and continuum survey of protoplanetary disk through to young debris disk systems carried out using the ACS instrument on the Herschel Space Observatory. This Open Time Key program, known as GASPS (Gas Survey of Protoplanetary Systems), targeted approx. 250 young stars in narrow wavelength regions covering the [OI] fine structure line at 63 micron the brightest far-infrared line in such objects. A subset of the brightest targets were also surveyed in [OI]145 micron, [CII] at 157 μm, as well as several transitions of H2O and high-excitation CO lines at selected wavelengths between 78 and 180 micron. Additionally, GASPS included continuum photometry at 70, 100 and 160 micron, around the peak of the dust emission. The targets were SED Class II– III T Tauri stars and debris disks from seven nearby young associations, along with a comparable sample of isolated Herbig AeBe stars. The aim was to study the global gas and dust content in a wide sample of circumstellar disks, combining the results with models in a systematic way. In this overview paper we review the scientific aims, target selection and observing strategy of the program. We summarize some of the initial results, showing line identifications, listing the detections, and giving a first statistical study of line detectability. The [OI] line at 63 micron was the brightest line seen in almost all objects, by a factor of ∼10. Overall [OI] 63 micron detection rates were 49%, with 100% of HAeBe stars and 43% of T Tauri stars detected. A comparison with published disk dust masses (derived mainly from sub-mm continuum, assuming standard values of the mm mass opacity) shows a dust mass threshold for [OI] 63 μm detection of approx.10(exp -5) Solar M.. Normalizing to a distance of 140 pc, 84% of objects with dust masses ≥10 (exp -5) Solar M can be detected in this line in the present survey; 32% of those of mass 10(exp -6) – 10 (exp -5) Solar M, and only a very small number of unusual objects with lower masses can be detected. This is consistent with models with a moderate UV excess and disk flaring. For a given disk mass, [OI] detectability is lower for M stars compared with earlier spectral types. Both the continuum and line emission was, in most systems, spatially and spectrally unresolved and centered on the star, suggesting that emission in most cases was from the disk. Approximately 10 objects showed resolved emission, most likely from outflows. In the GASPS sample, [OI] detection rates in T Tauri associations in the 0.3–4 Myr age range were approx. 50%. For each association in the 5–20 Myr age range, approx. 2 stars remain detectable in [OI] 63 micron, and no systems were detected in associations with age >20 Myr. Comparing with the total number of young stars in each association, and assuming a ISM-like gas/dust ratio, this indicates that approx. 18% of stars retain a gas-rich disk of total mass approx. Jupiter- M for 1–4 Myr, 1–7% keep such disks for 5–10 Myr, but none are detected beyond 10–20 Myr. The brightest [OI] objects from GASPS were also observed in [OI]145 micron, [CII]157 micron and CO J = 18- 17, with detection rates of 20–40%. Detection of the [CII] line was not correlated with disk mass, suggesting it arises more commonly from a compact remnant envelope.

Herschel Space Observatory↗

NASA NTRS2020

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗

NASA NTRS2012

An Automated Algorithm to Screen Massive Training Samples for a Global Impervious Surface Classification

An algorithm is developed to automatically screen the outliers from massive training samples for Global Land Survey - Imperviousness Mapping Project (GLS-IMP). GLS-IMP is to produce a global 30 m spatial resolution impervious cover data set for years 2000 and 2010 based on the Landsat Global Land Survey (GLS) data set. This unprecedented high resolution impervious cover data set is not only significant to the urbanization studies but also desired by the global carbon, hydrology, and energy balance researches. A supervised classification method, regression tree, is applied in this project. A set of accurate training samples is the key to the supervised classifications. Here we developed the global scale training samples from 1 m or so resolution fine resolution satellite data (Quickbird and Worldview2), and then aggregate the fine resolution impervious cover map to 30 m resolution. In order to improve the classification accuracy, the training samples should be screened before used to train the regression tree. It is impossible to manually screen 30 m resolution training samples collected globally. For example, in Europe only, there are 174 training sites. The size of the sites ranges from 4.5 km by 4.5 km to 8.1 km by 3.6 km. The amount training samples are over six millions. Therefore, we develop this automated statistic based algorithm to screen the training samples in two levels: site and scene level. At the site level, all the training samples are divided to 10 groups according to the percentage of the impervious surface within a sample pixel. The samples following in each 10% forms one group. For each group, both univariate and multivariate outliers are detected and removed. Then the screen process escalates to the scene level. A similar screen process but with a looser threshold is applied on the scene level considering the possible variance due to the site difference. We do not perform the screen process across the scenes because the scenes might vary due to the phenology, solar-view geometry, and atmospheric condition etc. factors but not actual landcover difference. Finally, we will compare the classification results from screened and unscreened training samples to assess the improvement achieved by cleaning up the training samples. Keywords:

Tan, Bin↗

NASA NTRS2021

Unsteady Pressure-Sensitive-Paint Shot Noise Reduction

In the study of unsteady Pressure-Sensitive Paint (uPSP) image data sets, pixel intensity values collected by Phantom cameras from multiple perspectives are used to calculate the surface pressure of an experimental model. This paper focuses on mitigation of shot noise which is introduced into the image data set by the cameras. Shot noise impacts the quality of data collected, increasing the system error in its pressure value calculations. This paper characterises statistical methods to reduce shot noise within experimental data by taking an average or weighted average across an area of pixels of their respective counts. Each combined-pixel has a higher total effective intensity value and therefore exhibits reduced noise throughout its time history. This is demonstrated by calculating the standard deviations of the combined pixel values through time, which decrease as the combined pixel area increases. Finally, these methods are compared across experimental wind tunnel conditions to examine their effectiveness at different Mach numbers. Ultimately, the results are used to recommend a blur filter type and size which will be applied to the raw image data prior to processing, as well as a minimum camera resolution for future experiments. Downstream of the filter application, the software to convert the 2D image data sets into surface pressure readings on a 3D grid use a ratio of decimation that assigns multiple pixels to each 3D grid node. This decimation ratio will also be aligned with the size of the blur filter, resulting in a more manageable data set size and compromised spatial resolution. Combined with the effect of the blur filter, the resulting pixel intensity readings will exhibit significantly reduced shot noise, reducing the magnitude of uncertainty and error in the final calculated pressure. For future tests, the camera spatial resolution must be sufficient to capture small flow and model features even after the recommended blur filter is applied. The image data set is taken from a uPSP experiment of the Space Launch System Block 1 cargo vehicle model in September 2019.

Lucy Tang↗

NASA NTRS2021

Unsteady Pressure-Sensitive Paint (UPSP) Shot Noise Reduction

Lucy Tang↗

NASA NTRS2016

Evaluation of Different MODIS AOD Retrieval Algorithms for PM (sub 2.5) Estimation in the Western, Midwestern and Southeastern United States with Implications for Public Health

This study was part of the research activities of the Center for Applied Atmospheric Research and Education (CAARE) funded by the NASA MUREP (Minority University Research and Education Project) Institutional Research Opportunity (MIRO) Program. Satellite measurements of Aerosol Optical Depth (AOD) have been shown to be correlated with ground measurements of fine particulate matter less than 2.5 microns PM (sub 2.5), which in turn has been linked to respiratory and heart diseases. The strength of the correlation between AOD and PM (sub 2.5) varies for different AOD retrieval algorithms and geographic regions. We evaluated several Moderate Resolution Imaging Spectrometer (MODIS) AOD products from different satellites (Aqua vs. Terra), retrieval algorithms (Dark Target versus Deep Blue), Collections (5.1 versus 6) and spatial resolutions (10-kilometers versus 3-kilometers) for cities in the Western, Midwestern and Southeastern U.S. We developed and validated PM (sub 2.5) prediction models using remotely-sensed AOD data, which were improved by incorporating meteorological variables (temperature, relative humidity, precipitation, wind speed, and wind direction) from the North American Land Data Assimilation System Phase 2 (NLDAS-2). Adding these meteorological data significantly improved the predictive power of all the PM (sub 2.5) models, especially in the Western U.S. Temperature, relative humidity and wind speed were the most significant meteorological variables throughout the year in the Western U.S. Wind speed was the most significant meteorological variable for the cold season while temperature was the most significant variable for the warm season in the Midwestern and Southeastern U.S. Our study re-establishes the connection between PM (sub 2.5) and public health concerns including respiratory and cardiovascular diseases (asthma, high blood pressure, coronary heart disease, heart attack, and stroke). Using PM (sub 2.5) data and health data from the Centers for Disease Control and Prevention (CDC)'s Behavioral Risk Factor Surveillance System (BRFSS), our statistical analysis showed that heart attack and stroke occurrences had the strongest correlations with PM (sub 2.5).

Environmental Public Health↗

NASA NTRS2015

Extension and Statistical Analysis of the GACP Aerosol Optical Thickness Record.

The primary product of the Global Aerosol Climatology Project (GACP) is a continuous record of the aerosol optical thickness (AOT) over the oceans. It is based on channel-1 and -2 radiance data from the Advanced Very High Resolution Radiometer (AVHRR) instruments flown on successive National Oceanic and Atmospheric Administration (NOAA) platforms. We extend the previous GACP dataset by four years through the end of 2009 using NOAA-17 and -18 AVHRR radiances recalibrated against MODerate resolution Imaging Spectroradiometer (MODIS) radiance data, thereby making the GACP record almost three decades long. The temporal overlap of over three years of the new NOAA-17 and the previous NOAA-16 record reveals an excellent agreement of the corresponding global monthly mean AOT values, thereby confirming the robustness of the vicarious radiance calibration used in the original GACP product. The temporal overlap of the NOAA-17 and -18 instruments is used to introduce a small additive adjustment to the channel-2 calibration of the latter resulting in a consistent record with increased data density. The Principal Component Analysis (PCA) of the newly extended GACP record shows that most of the volcanic AOT variability can be isolated into one mode responsible for ~12% of the total variance. This conclusion is confirmed by a combined PCA analysis of the GACP, MODIS, andMulti-angle Imaging SpectroRadiometer (MISR) AOTs during the volcano-free period fromFebruary 2000 to December 2009.We show that the modes responsible for the tropospheric AOT variability in the three datasets agree well in terms of correlation and spatial patterns. A previously identified negative AOT trend which started in the late 1980s and continued into the early 2000s is confirmed. Its magnitude and duration indicate that it was caused by changes in tropospheric aerosols. The latest multi-satellite segment of the GACP record shows that this trend tapered off, with no noticeable AOT change after 2002. This result is consistent with the MODIS andMISR AOT records as well as with the recent gradual reversal frombrightening to dimming revealed by surface flux measurements in many aerosol producing regions. Thus the robustness of the GACP record is confirmed, increasing our confidence in the validity of the negative trend. Although the nominal negative GACP AOT trend could partially be an artifact of increasing aerosol absorption, we argue that the time dependence of the GACP record, including the latest flat period, is more consistent with the actual decrease in the tropospheric AOT.

aerosols↗

NASA NTRS2011

Phase-Sensitive Coherence and the Classical-Quantum Boundary in Ghost Imaging

The theory of partial coherence has a long and storied history in classical statistical optics. the vast majority of this work addresses fields that are statistically stationary in time, hence their complex envelopes only have phase-insensitive correlations. The quantum optics of squeezed-state generation, however, depends on nonlinear interactions producing baseband field operators with phase-insensitive and phase-sensitive correlations. Utilizing quantum light to enhance imaging has been a topic of considerable current interest, much of it involving biphotons, i.e., streams of entangled-photon pairs. Biphotons have been employed for quantum versions of optical coherence tomography, ghost imaging, holography, and lithography. However, their seemingly quantum features have been mimicked with classical-sate light, questioning wherein lies the classical-quantum boundary. We have shown, for the case of Gaussian-state light, that this boundary is intimately connected to the theory of phase-sensitive partial coherence. Here we present that theory, contrasting it with the familiar case of phase-insensitive partial coherence, and use it to elucidate the classical-quantum boundary of ghost imaging. We show, both theoretically and experimentally, that classical phase-sensitive light produces ghost imaging most closely mimicking those obtained in biphotons, and we derived the spatial resolution, image contrast, and signal-to-noise ratio of a standoff-sensing ghost imager, taking into account target-induced speckle.

coherence theory↗

NASA NTRS2007

Utility of AIRS Retrievals for Climate Studies

Satellites provide an ideal platform to study the Earth-atmosphere system on practically all spatial and temporal scales. Thus, one may expect that their rapidly growing datasets could provide crucial insights not only for short-term weather processes/predictions but into ongoing and future climate change processes as well. Though Earth-observing satellites have been around for decades, extracting climatically reliable information from their widely varying datasets faces rather formidable challenges. AIRS/AMSU is a state of the art infrared/microwave sounding system that was launched on the EOS Aqua platform on May 4, 2002, and has been providing operational quality measurements since September 2002. In addition to temperature and atmospheric constituent profiles, outgoing longwave radiation and basic cloud parameters are also derived from the AIRS/AMSU observations. However, so far the AIRS products have not been rigorously evaluated and/or validated on a large scale. Here we present preliminary assessments of monthly and 8-day mean AIRS "Version 4.0" retrieved products (available to the public through the DAAC at NASA/GSFC) to assess their utility for climate studies. First we present "consistency checks" by evaluating the time series of means, and "anomalies" (relative to the first 4 full years' worth of AIRS "climate statistics") of several climatically important retrieved parameters. Finally, we also present preliminary results regarding interrelationships of some of these geophysical variables, to assess to what extent they are consistent with the known physics of climate variability/change. In particular, we find at least one observed relationship which contradicts current general circulation climate (GCM) model results: the global water vapor climate feedback which is expected to be strongly positive is deduced to be slightly negative (shades of the "Lindzen effect"?). Though the current AIRS climatology covers only -4.5 years, it will hopefully extend much further into the future.

Molnar, Guyla I.↗

NASA NTRS2000

Optical Performance Modeling of FUSE Telescope Mirror

We describe the Metrology Data Processor (METDAT), the Optical Surface Analysis Code (OSAC), and their application to the image evaluation of the Far Ultraviolet Spectroscopic Explorer (FUSE) mirrors. The FUSE instrument - designed and developed by the Johns Hopkins University and launched in June 1999 is an astrophysics satellite which provides high resolution spectra (lambda/Delta(lambda) = 20,000 - 25,000) in the wavelength region from 90.5 to 118.7 nm The FUSE instrument is comprised of four co-aligned, normal incidence, off-axis parabolic mirrors, four Rowland circle spectrograph channels with holographic gratings, and delay line microchannel plate detectors. The OSAC code provides a comprehensive analysis of optical system performance, including the effects of optical surface misalignments, low spatial frequency deformations described by discrete polynomial terms, mid- and high-spatial frequency deformations (surface roughness), and diffraction due to the finite size of the aperture. Both normal incidence (traditionally infrared, visible, and near ultraviolet mirror systems) and grazing incidence (x-ray mirror systems) systems can be analyzed. The code also properly accounts for reflectance losses on the mirror surfaces. Low frequency surface errors are described in OSAC by using Zernike polynomials for normal incidence mirrors and Legendre-Fourier polynomials for grazing incidence mirrors. The scatter analysis of the mirror is based on scalar scatter theory. The program accepts simple autocovariance (ACV) function models or power spectral density (PSD) models derived from mirror surface metrology data as input to the scatter calculation. The end product of the program is a user-defined pixel array containing the system Point Spread Function (PSF). The METDAT routine is used in conjunction with the OSAC program. This code reads in laboratory metrology data in a normalized format. The code then fits the data using Zernike polynomials for normal incidence systems or Legendre-Fourier polynomials for grazing incidence systems. It removes low order terms from the metrology data, calculates statistical ACV or PSD functions, and fits these data to OSAC models for the scatter analysis. In this paper we briefly describe the laboratory image testing of FUSE spare mirror performed in the near and vacuum ultraviolet at John Hopkins University and OSAC modeling of the test setup performed at NASA/GSFC. The test setup is a double-pass configuration consisting of a Hg discharge source, the FUSE off-axis parabolic mirror under test, an autocollimating flat mirror, and a tomographic imaging detector. Two additional, small fold flats are used in the optical train to accommodate the light source and the detector. The modeling is based on Zernike fitting and PSD analysis of surface metrology data measured by both the mirror vendor (Tinsley) and JHU. The results of our models agree well with the laboratory imaging data, thus validating our theoretical model. Finally, we predict the imaging performance of FUSE mirrors in their flight configuration at far-ultraviolet wavelengths.

Saha, Timo T.↗