Search NASA⌕ Search

SEARCH · Search NASA

Results for “sampling 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 1,639 records · Page 91

The Prompt-afterglow Connection in Gamma-ray Bursts: a Comprehensive Statistical Analysis of Swift X-ray Light-curves

We present a comprehensive statistical analysis of Swift X-ray light-curves of Gamma- Ray Bursts (GRBs) collecting data from more than 650 GRBs discovered by Swift and other facilities. The unprecedented sample size allows us to constrain the rest-frame X-ray properties of GRBs from a statistical perspective, with particular reference to intrinsic time scales and the energetics of the different light-curve phases in a common rest-frame 0.3-30 keV energy band. Temporal variability episodes are also studied and their properties constrained. Two fundamental questions drive this effort: i) Does the X-ray emission retain any kind of "memory" of the prompt γ-ray phase? ii) Where is the dividing line between long and short GRB X-ray properties? We show that short GRBs decay faster, are less luminous and less energetic than long GRBs in the X-rays, but are interestingly characterized by similar intrinsic absorption. We furthermore reveal the existence of a number of statistically significant relations that link the X-ray to prompt γ-ray parameters in long GRBs; short GRBs are outliers of the majority of these 2-parameter relations. However and more importantly, we report on the existence of a universal 3-parameter scaling that links the X-ray and the γ-ray energy to the prompt spectral peak energy of both long and short GRBs: E(sub X,iso)∝ E(sup 1.00+/-0.06)(sub γ,iso) /E(sup 0.60+/-0.10)(sub pk).

Gamma-Ray Bursts↗

First Measurement of 87 Rb( α , xn ) Cross Sections at Weak r -process Energies in Supernova ν -driven Ejecta to Investigate Elemental Abundances in Low-metallicity Stars

Observed abundances of Z ∼ 40 elements in metal-poor stars vary from star to star, indicating that the rapid and slow neutron capture processes may not contribute alone to the synthesis of elements beyond iron. The weak r-process was proposed to produce Z ∼ 40 elements in a subset of old stars. Thought to occur in the ν-driven ejecta of a core-collapse supernova, ( α, xn ) reactions would drive the nuclear flow toward heavier masses at T = 2−5 GK. However, current comparisons between modeled and observed yields do not bring satisfactory insights into the stellar environment, mainly due to the uncertainties of the nuclear physics inputs where the dispersion in a given reaction rate often exceeds 1 order of magnitude. Involved rates are calculated with the statistical model where the choice of an α -optical-model potential ( α OMP) leads to such a poor precision. The first experiment on 87 Rb( α, xn ) reactions at weak r -process energies is reported here. Total inclusive cross sections were assessed at E c.m. = 8.1−13 MeV (3.7−7.6 GK) with the active target MUlti-Sampling Ionization Chamber. With an N = 50 seed nucleus, the measured values agree with statistical model estimates using the α OMP Atomki-V2. A reevaluated reaction rate was incorporated into new nucleosynthesis calculations, focusing on ν-driven ejecta conditions known to be sensitive to this specific rate. These conditions were found to fail to reproduce the lighter heavy element abundances in metal-poor stars.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-I-33: Establishment of CT Diagnostic Reference Levels in Province Nova Scotia

Purpose: To evaluate patient radiation dose from the most frequently performed CT examinations and to establish provincial diagnostic reference levels (DRLs) as a tool for protocol optimization. Methods: The study investigated the following CT examinations: head, chest, abdomen/pelvis, and chest/abdomen/pelvis (CAP). Dose data, volume CT dose index (CTDIvol) and dose-length product (DLP), were collected from 15 CT scanners installed during 2004–2014 in 11 hospital sites of Nova Scotia. All scanners had dose modulation options and multislice capability (16–128 detector rows). The sample for each protocol included 15 average size patients (70±20 kg). Provincial DRLs were calculated as the 75th percentile of patient dose distributions. The differences in dose between hospitals were evaluated with a single factor ANOVA statistical test. Generalized linear modeling was used to determine the factors associated with higher radiation dose. A sample of 36 abdominal studies performed on three different scanners was blinded and randomized for an assessment by an experienced radiologist who graded the imaging quality of anatomic structures. Results: Data for 900 patients were collected. The DRLs were proposed using CTDIvol (mGy) and DLP (mGy*cm) values for CT head (67 and 1049, respectively), chest (12 and 393), abdomen/pelvis (16 and 717), and CAP (14 and 1034). These DRLs were lower than the published national data except for the head CTDIvol. The differences between the means of the dose distributions from each scanner were statistically significant (p<0.05) for all examinations. A very weak correlation was found between the dose and the scanner age or the number of slices with Pearson’s correlation coefficients of 0.011–0.315. The blinded analysis of image quality demonstrated no clinically significant difference except for the noise category. Conclusion: Provincial DRLs were established for typical CT examinations. The variations in dose between the hospitals suggested a large potential for optimization of examinations. Radiology Research Foundation grant.

60 APPLIED LIFE SCIENCES↗

The white dwarfs within 25 pc of the Sun: Kinematics and spectroscopic subtypes

We present the fractional distribution of spectroscopic subtypes, range and distribution of surface temperatures, and kinematical properties of the white dwarfs (WDs) within 25 pc of the Sun. There is no convincing evidence of halo WDs in the total 25 pc sample of 224 WDs. There is also little to suggest the presence of genuine thick disk subcomponent members within 25 pc. It appears that the entire 25 pc sample likely belongs to the thin disk. We also find no significant kinematic differences with respect to spectroscopic subtypes. The total DA to non-DA ratio of the 25 pc sample is 1.8, a manifestation of deepening envelope convection, which transforms DA stars with sufficiently thin H surface layers into non-DAs. We compare this ratio with the results of other studies. We find that at least 11% of the WDs within 25 pc of the Sun (the DAZ and DZ stars) have photospheric metals that likely originate from accretion of circumstellar material (debris disks) around them. If this interpretation is correct, then it suggests the possibility that a similar percentage have planets, asteroid-like bodies, or debris disks orbiting them. Our volume-limited sample reveals a pileup of DC WDs at the well-known cutoff in DQ WDs at T {sub eff} ∼ 6000 K. Mindful of small number statistics, we speculate on its possible evolutionary significance. We find that the incidence of magnetic WDs in the 25 pc sample is at least 8% in our volume-limited sample, dominated by cool WDs. We derive approximate formation rates of DB and DQ degenerates and present a preliminary test of the evolutionary scenario that all cooling DB stars become DQ WDs via helium convective dredge-up with the diffusion tail of carbon extending upward from their cores.

79 ASTRONOMY AND ASTROPHYSICS↗

The NASA MSFC Earth Global Reference Atmospheric Model-2007 Version

Reference or standard atmospheric models have long been used for design and mission planning of various aerospace systems. The NASA/Marshall Space Flight Center (MSFC) Global Reference Atmospheric Model (GRAM) was developed in response to the need for a design reference atmosphere that provides complete global geographical variability, and complete altitude coverage (surface to orbital altitudes) as well as complete seasonal and monthly variability of the thermodynamic variables and wind components. A unique feature of GRAM is that, addition to providing the geographical, height, and monthly variation of the mean atmospheric state, it includes the ability to simulate spatial and temporal perturbations in these atmospheric parameters (e.g. fluctuations due to turbulence and other atmospheric perturbation phenomena). A summary comparing GRAM features to characteristics and features of other reference or standard atmospheric models, can be found Guide to Reference and Standard Atmosphere Models. The original GRAM has undergone a series of improvements over the years with recent additions and changes. The software program is called Earth-GRAM2007 to distinguish it from similar programs for other bodies (e.g. Mars, Venus, Neptune, and Titan). However, in order to make this Technical Memorandum (TM) more readable, the software will be referred to simply as GRAM07 or GRAM unless additional clarity is needed. Section 1 provides an overview of the basic features of GRAM07 including the newly added features. Section 2 provides a more detailed description of GRAM07 and how the model output generated. Section 3 presents sample results. Appendices A and B describe the Global Upper Air Climatic Atlas (GUACA) data and the Global Gridded Air Statistics (GGUAS) database. Appendix C provides instructions for compiling and running GRAM07. Appendix D gives a description of the required NAMELIST format input. Appendix E gives sample output. Appendix F provides a list of available parameters to enable the user to generate special output. Appendix G gives an example and guidance on incorporating GRAM07 as a subroutine in other programs such as trajectory codes or orbital propagation routines.

Leslie, F.W.↗

Partial Least Squares, Experimental Design, and Near-Infrared Spectrophotometry for the Remote Quantification of Nitric Acid Concentration and Temperature

Near-infrared spectrophotometry and partial least squares regression (PLSR) were evaluated to create a pleasantly simple yet effective approach for measuring HNO3 concentration with varying temperature levels. A training set, which covered HNO3 concentrations (0.1–8 M) and temperature (10–40 °C), was selected using a D-optimal design to minimize the number of samples required in the calibration set for PLSR analysis. The top D-optimal-selected PLSR models had root mean squared error of prediction values of 1.4% for HNO 3 and 4.0% for temperature. The PLSR models built from spectra collected on static samples were validated against flow tests including HNO 3 concentration and temperature gradients to test abnormal conditions (e.g., bubbles) and the model performance between sample points in the factor space. Based on cross-validation and prediction modeling statistics, the designed near-infrared absorption approach can provide remote, quantitative analysis of HNO 3 concentration and temperature for production-oriented applications in facilities where laser safety challenges would inhibit the implementation of other optical techniques (e.g., Raman spectroscopy) and in which space, time, and/or resources are constrained. The experimental design approach effectively minimized the number of samples in the training set and maintained or improved PLSR model performance, which makes the described chemometric approach more amenable to nuclear field applications.

07 ISOTOPE AND RADIATION SOURCES↗

The Atacama Cosmology Telescope: Sunyaev-Zel'dovich selected galaxy clusters at 148 GHz from three seasons of data

We present a catalog of 68 galaxy clusters, of which 19 are new discoveries, detected via the Sunyaev-Zel'dovich effect (SZ) at 148 GHz in the Atacama Cosmology Telescope (ACT) survey on the celestial equator. With this addition, the ACT collaboration has reported a total of 91 optically confirmed, SZ detected clusters. The 504 square degree survey region includes 270 square degrees of overlap with SDSS Stripe 82, permitting the confirmation of SZ cluster candidates in deep archival optical data. The subsample of 48 clusters within Stripe 82 is estimated to be 90% complete for M{sub 500c} > 4.5 × 10{sup 14}M{sub s}un and redshifts 0.15 < z < 0.8. While a full suite of matched filters is used to detect the clusters, the sample is studied further through a ''Profile Based Amplitude Analysis'' using a statistic derived from a single filter at a fixed θ{sub 500} = 5.'9 angular scale. This new approach incorporates the cluster redshift along with prior information on the cluster pressure profile to fix the relationship between the cluster characteristic size (R{sub 500}) and the integrated Compton parameter (Y{sub 500}). We adopt a one-parameter family of ''Universal Pressure Profiles'' (UPP) with associated scaling laws, derived from X-ray measurements of nearby clusters, as a baseline model. Three additional models of cluster physics are used to investigate a range of scaling relations beyond the UPP prescription. Assuming a concordance cosmology, the UPP scalings are found to be nearly identical to an adiabatic model, while a model incorporating non-thermal pressure better matches dynamical mass measurements and masses from the South Pole Telescope. A high signal to noise ratio subsample of 15 ACT clusters with complete optical follow-up is used to obtain cosmological constraints. We demonstrate, using fixed scaling relations, how the constraints depend on the assumed gas model if only SZ measurements are used, and show that constraints from SZ data are limited by uncertainty in the scaling relation parameters rather than sample size or measurement uncertainty. We next add in seven clusters from the ACT Southern survey, including their dynamical mass measurements, which are based on galaxy velocity dispersions and thus are independent of the gas physics. In combination with WMAP7 these data simultaneously constrain the scaling relation and cosmological parameters, yielding 68% confidence ranges described by σ{sub 8} = 0.829 ± 0.024 and Ω{sub m} = 0.292 ± 0.025.. We consider these results in the context of constraints from CMB and other cluster studies. The constraints arise mainly due to the inclusion of the dynamical mass information and do not require strong priors on the SZ scaling relation parameters. The results include marginalization over a 15% bias in dynamical masses relative to the true halo mass. In an extension to ΛCDM that incorporates non-zero neutrino mass density, we combine our data with WMAP7, Baryon Acoustic Oscillation data, and Hubble constant measurements to constrain the sum of the neutrino mass species to be Σ{sub ν}m{sub ν} < 0.29 eV (95% confidence limit)

79 ASTRONOMY AND ASTROPHYSICS↗

Missing data in multi-omics integration: Recent advances through artificial intelligence

Biological systems function through complex interactions between various ‘omics (biomolecules), and a more complete understanding of these systems is only possible through an integrated, multi-omic perspective. This has presented the need for the development of integration approaches that are able to capture the complex, often non-linear, interactions that define these biological systems and are adapted to the challenges of combining the heterogenous data across ‘omic views. A principal challenge to multi-omic integration is missing data because all biomolecules are not measured in all samples. Due to either cost, instrument sensitivity, or other experimental factors, data for a biological sample may be missing for one or more ‘omic techologies. Recent methodological developments in artificial intelligence and statistical learning have greatly facilitated the analyses of multi-omics data, however many of these techniques assume access to completely observed data. A subset of these methods incorporate mechanisms for handling partially observed samples, and these methods are the focus of this review. We describe recently developed approaches, noting their primary use cases and highlighting each method's approach to handling missing data. We additionally provide an overview of the more traditional missing data workflows and their limitations; and we discuss potential avenues for further developments as well as how the missing data issue and its current solutions may generalize beyond the multi-omics context.

97 MATHEMATICS AND COMPUTING↗

On-line estimation of error covariance parameters for atmospheric data assimilation

A simple scheme is presented for on-line estimation of covariance parameters in statistical data assimilation systems. The scheme is based on a maximum-likelihood approach in which estimates are produced on the basis of a single batch of simultaneous observations. Simple-sample covariance estimation is reasonable as long as the number of available observations exceeds the number of tunable parameters by two or three orders of magnitude. Not much is known at present about model error associated with actual forecast systems. Our scheme can be used to estimate some important statistical model error parameters such as regionally averaged variances or characteristic correlation length scales. The advantage of the single-sample approach is that it does not rely on any assumptions about the temporal behavior of the covariance parameters: time-dependent parameter estimates can be continuously adjusted on the basis of current observations. This is of practical importance since it is likely to be the case that both model error and observation error strongly depend on the actual state of the atmosphere. The single-sample estimation scheme can be incorporated into any four-dimensional statistical data assimilation system that involves explicit calculation of forecast error covariances, including optimal interpolation (OI) and the simplified Kalman filter (SKF). The computational cost of the scheme is high but not prohibitive; on-line estimation of one or two covariance parameters in each analysis box of an operational bozed-OI system is currently feasible. A number of numerical experiments performed with an adaptive SKF and an adaptive version of OI, using a linear two-dimensional shallow-water model and artificially generated model error are described. The performance of the nonadaptive versions of these methods turns out to depend rather strongly on correct specification of model error parameters. These parameters are estimated under a variety of conditions, including uniformly distributed model error and time-dependent model error statistics.

Dee, Dick P.↗

A Measurement of Gravitational Lensing of the Cosmic Microwave Background by Galaxy Clusters Using Data from the South Pole Telescope

Clusters of galaxies are expected to gravitationally lens the cosmic microwave background (CMB) and thereby generate a distinct signal in the CMB on arcminute scales. Measurements of this effect can be used to constrain the masses of galaxy clusters with CMB data alone. Here we present a measurement of lensing of the CMB by galaxy clusters using data from the South Pole Telescope (SPT). We also develop a maximum likelihood approach to extract the CMB cluster lensing signal and validate the method on mock data. We quantify the effects on our analysis of several potential sources of systematic error and find that they generally act to reduce the best-fit cluster mass. It is estimated that this bias to lower cluster mass is roughly 0.85σ in units of the statistical error bar, although this estimate should be viewed as an upper limit. Furthermore, we apply our maximum likelihood technique to 513 clusters selected via their Sunyaev–Zeldovich (SZ) signatures in SPT data, and rule out the null hypothesis of no lensing at 3.1σ. The lensing-derived mass estimate for the full cluster sample is consistent with that inferred from the SZ flux: M 200,lens = 0.83 +0.38 -0.37 M 200,SZ (68% C.L., statistical error only).

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization and CST Batch Contact Equilibrium Testing of Modified Tank 9H Process Supernate Samples in Support of TCCR

The Tank Closure Cesium Removal (TCCR) system uses ion exchange columns filled with crystalline silicotitanate (CST) media to process radioactive waste solutions for the removal of ¹³⁷Cs. TCCR currently focuses on dissolving Savannah River Site (SRS) radioactive tank waste (primarily sodium saltcake solids) within Tank 10H followed by at-tank ion exchange column treatment. Two supernate batches from Tank 10H have been processed through the TCCR unit and processing of a third batch is expected soon. After processing of this supernate batch, plans are to replace the CST columns and process dissolved salt solution from Tank 9H through Tank 10H and then through the new CST columns installed in the TCCR unit. The new columns are expected to contain either a media similar to an archived CST batch (IE-911) or the R9120-B CST media used in the current TCCR columns (the two materials are fundamentally the same; just different specifications, product names, and preconditioning steps). Samples of Tank 9H dissolved saltcake were received at the Savannah River National Laboratory (SRNL) and characterized. The Tank 9H supernate contained a high sodium concentration (~9.6 M Na⁺) and will require dilution to near 6 M [Na⁺] prior to processing through the TCCR unit. The cesium concentration in Tank 9H is currently significantly higher than was observed with Tank 10H. Three dilutions of the Tank 9H supernate were conducted to mimic possible dilutions that could be conducted in the tank farm prior to TCCR processing using inhibited water, sodium hydroxide, and sodium nitrate solutions. Dilution #1 was prepared by diluting the Tank 9H supernate by a factor of 1.6 with inhibited water, Dilution #2 was prepared by diluting the Tank 9H supernate by a factor of 3.7 with a mixed sodium hydroxide/sodium nitrate diluent, and Dilution #3 was prepared by also diluting by a factor of 3.7 but with sodium hydroxide only. All three dilutions have similar Na⁺ concentrations (~6 M), but Dilutions #2 and #3 contain significantly less cesium. Major and key minor components of the diluted Tank 9H supernate samples are provided in Table ES-1. Batch contact equilibrium tests were conducted with the Tank 9H dilutions and the two different CST media batches being considered for use in the new TCCR columns. Results are summarized in Table ES-2 and compared to ZAM model predictions. The highest cesium distribution coefficient and percent removal were observed with Dilution #3. IE-911 CST (an archived CST media batch) was more effective at removing cesium than the more recently prepared R9120-B CST media, though both media samples removed >88% of the cesium and the differences may not be statistically different considering the overall uncertainty. Based on the results, maximum cesium loadings were calculated for each CST media type and Tank 9H dilution. In general, maximum cesium loadings from dilutions of this supernate batch are quite high (approaching 0.1 mmol total Cs⁺/g CST for Dilution #1). The highest calculated maximum ¹³⁷Cs loading for the Tank 9H dilutions using the ZAM model with input of the tank compositions and batch contact results was 207 Ci/kg CST (Dilution #1 with IE-911 CST). In all cases, higher maximum cesium loading values were predicted for IE-911 CST versus R9120-B, although the differences varied considerably between the dilutions. The maximum loading value for IE-911 CST with Tank 9H Dilution #1 was only 7% higher than the maximum loading for R9120-B. The maximum loading value for IE-911 CST with Dilution #3 was 24% higher than the maximum loading for R9120-B. Dilution #2 was intermediate between these values. Note that these maximum loading values are the theoretical calculated values, and actual operating conditions will cause differences. A CST binder dilution (correction) factor near 0.7 (relative to pure powder CST) was required for each batch contact test with IE-911 engineered CST using the three Tank 9H dilutions. Correction factors calculated for R9120-B CST ranged from 0.56 to 0.66 for the three dilutions. Following the batch contact tests, the CST samples were isolated from the Tank 9H solution, washed, dried, and digested in acid following established procedures. The analysis results are provided in Table ES-3. Total cesium loading values determined by CST digestion were similar to the calculated loading values based on solution analysis (Table ES-2) for all samples. In addition, the CST was observed to load calcium (R9120-B sample only), iron, strontium, lead, uranium, and plutonium after contact with the waste supernate, as has been observed previously.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Adventures in Parallel Processing: Entry, Descent and Landing Simulation for the Genesis and Stardust Missions

This paper will describe the Entry, Descent and Landing simulation tradeoffs and techniques that were used to provide the Monte Carlo data required to approve entry during a critical period just before entry of the Genesis Sample Return Capsule. The same techniques will be used again when Stardust returns on January 15, 2006. Only one hour was available for the simulation which propagated 2000 dispersed entry states to the ground. Creative simulation tradeoffs combined with parallel processing were needed to provide the landing footprint statistics that were an essential part of the Go/NoGo decision that authorized release of the Sample Return Capsule a few hours before entry.

Lyons, Daniel T.↗

Porting Adaptive Ensemble Molecular Dynamics Workflows to the Summit Supercomputer

Molecular dynamics (MD) simulations must take very small (femtosecond) integration steps in simulation-time to avoid numerical errors. Efficient use of parallel programming models and accelerators in state-of-the art MD programs now is pushing Moore’s limit for time-per-MD step. As a result, directly simulating timescales beyond milliseconds will not be attainable directly, even at exascale. However, concepts from statistical physics can be used to combine many parallel simulations to provide information about longer timescales and to adequately sample the simulation space, while preserving details about the dynamics of the system. Implementing such an approach requires a workflow program that allows adaptable steering of task assignments based on extensive statistical analysis of intermediate results. Here we report the implementation of such an adaptable workflow program to drive simulations on the Summit IBM Power System AC922, a pre-exascale supercomputer at the Oak Ridge Leadership Computing Facility (OLCF). We compare to experiences on Titan, Summit’s predecessor, report the performance of the workflow and its components, and describe the porting process. We find that using a workflow program managed by a Mongo database can provide the fault tolerance, scalable performance, task dispatch rate, and reconfigurability required for robust and portable implementation of ensemble simulations such as are used in enhanced-sampling molecular dynamics. This type of workflow generator can also be used to provide adaptive steering of ensemble simulations for other applications in addition to MD.

Ossyra, John R.↗

Statistical Characterization of Hot Jupiter Atmospheres Using Spitzer's Secondary Eclipses

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6 and 4.5 μm using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for James Webb Space Telescope such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit (e cos w = 0.00614 ± 0.00064), and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan & Agol (2011). Planets having maximum day-side temperatures exceeding ∼2200 K are consistent with having zero albedo and a distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 ± 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800–2500 K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profiles of real planetary atmospheres as compared to models.

Exoplanet astronomy↗

In Situ and Lidar Observations of Tropopause Subvisible Cirrus Clouds During TC4

During the Tropical Composition, Clouds, and Climate Coupling (TC4) experiment in July-August 2007, the NASA WB-57F and ER-2 aircraft made coordinated flights through a tropopause subvisible cirrus (SVC) layer off the Pacific Coast of Central America. The ER-2 aircraft was equipped with a remote sensing payload that included the cloud physics lidar (CPL). The WB-57F payload included cloud microphysical and trace gas measurements, and the aircraft made four vertical profiles through the SVC layer shortly after the ER-2 flew over. The in situ and remotely sensed data are used to quantify the meteorological and microphysical properties of the SVC layer, and these data are compared to the limited set of SVC measurements that have previously been made. It is found that the layer encountered was particularly tenuous, with optical depths (tau) between about 10(exp -4) and 10(exp -3). From the in situ and other meteorological data, radiative heating rate perturbations of approx.0.05-0.1 K/day are calculated. These heating rates are smaller than previous estimates for tropopause SVC, consistent with the smaller tau in the present study. Coverage statistics based on CPL data from other TC4 flights indicate that this cloud was not an outlier among the sampled population. SVC with properties similar to the one presented here are below the detection limit of space \based lidars such as CALIPSO, and a comparison with the TC4 statistics suggests that a majority (>50%) of tropopause SVC (with tau < 0.01) could be unaccounted for in studies using CALIPSO data.

Davis, Sean↗

Orbital Trade Study for the PREFIRE Mission

The Polar Radiant Energy in the Far-Infrared Experiment (PREFIRE) is a low-cost cubesat-based Earth Science mission selected for rapid implementation through NASAs Earth Ventures-Instrument program. Two cubesats carrying thermal infrared spectrometers (TIRS) in polar orbits will test long-held assumptions regarding radiative balance at long wavelengths where the cold-dry Arctic and Antarctic surfaces expel significant amounts of the earth's radiant energy. Thermal radiometry with better than one micron sampling at wavelengths beyond 15 microns will provide critical observational bounds for inputs into glacier melt and climate change models through improvements in surface emissivities, longwave greenhouse effect parameterizations, and cloud presence/types. Additional information on potential far-infrared surface/atmosphere feedbacks during rapid melt processes may be captured with sub-diurnal revisits of targeted areas. A baseline mission with two cubesats in different orbits provides a means for optimization of sub-diurnal sampling for science purposes. In this presentation we provide details of orbit studies that inform expected sampling metrics. Candidate orbits are modeled in SGP and then TIRS sampling is projected into surface footprints with appropriate geo-location. Matching algorithms then compile statistics for time-delayed revisits between and within the two orbits. Binning of results by latitude allows for abstraction of precession issues as well as for correlations with surface types as defined by radiometric classification schemes. Multiple acceptable pairs of science-driven orbits were found that optimize sampling in areas known to experience rapid melt events. The final choice of PREFIRE orbital parameters includes the systems trade-space, particularly the power budget, as well as the availability of launch opportunities.

"L'Ecuyer, Tristan"↗

Modeling the Atmospheric Phase Effects of a Digital Antenna Array Communications System

In an antenna array system such as that used in the Deep Space Network (DSN) for satellite communication, it is often necessary to account for the effects due to the atmosphere. Typically, the atmosphere induces amplitude and phase fluctuations on the transmitted downlink signal that invalidate the assumed stationarity of the signal model. The degree to which these perturbations affect the stationarity of the model depends both on parameters of the atmosphere, including wind speed and turbulence strength, and on parameters of the communication system, such as the sampling rate used. In this article, we focus on modeling the atmospheric phase fluctuations in a digital antenna array communications system. Based on a continuous-time statistical model for the atmospheric phase effects, we show how to obtain a related discrete-time model based on sampling the continuous-time process. The effects of the nonstationarity of the resulting signal model are investigated using the sample matrix inversion (SMI) algorithm for minimum mean-squared error (MMSE) equalization of the received signal

Tkacenko, A.↗

Some Mathematical and Geophysical Considerations in Radioisotope Dating Applications

Some quality considerations for use in isotopic dating are presented to identify and correct heretofore unidentified overestimate scenarios. These include to a lesser degree the statistical interpretation issues with linear-least-squares fitting results but more importantly the isotope effect in the individual components of the isochron coefficient ratios. By taking into consideration the isotope effect (differential mass diffusion rates) when measuring isotopic ratios from very old samples, the distribution dependency in the coefficient ratios will cause a bias if isotopic diffusion rates are not identical throughout a sample. The isotope effect is that isotopes having a smaller atomic mass will diffuse faster throughout a medium than will their heavier counterparts causing concentration gradients of their ratios even when there are no contributions from radioactive decay. The application to Rb/Sr dating is evaluated and shown to result in expected age overestimates when isotopic ratios are employed to linearize the isochron. A suggested method to test for this effect is argued to require rigorous statistical analysis. An associated optimal sampling technique would involve using single-grain etching. It is also shown that the only method to fully eliminate the isotope effect is to not use isotopic ratios at all in radioisotopic dating as the physics do not require the use of isotopic ratios for geochronological dating. However, without the ratios, the data are inherently noisy.

74 ATOMIC AND MOLECULAR PHYSICS↗