Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model fitting”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS2006

Optimal Design of Passive Flow Control for a Boundary-Layer-Ingesting Offset Inlet Using Design-of-Experiments

This research will investigate the use of Design-of-Experiments (DOE) in the development of an optimal passive flow control vane design for a boundary-layer-ingesting (BLI) offset inlet in transonic flow. This inlet flow control is designed to minimize the engine fan face distortion levels and first five Fourier harmonic half amplitudes while maximizing the inlet pressure recovery. Numerical simulations of the BLI inlet are computed using the Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, developed at NASA. These simulations are used to generate the numerical experiments for the DOE response surface model. In this investigation, two DOE optimizations were performed using a D-Optimal Response Surface model. The first DOE optimization was performed using four design factors which were vane height and angles-of-attack for two groups of vanes. One group of vanes was placed at the bottom of the inlet and a second group symmetrically on the sides. The DOE design was performed for a BLI inlet with a free-stream Mach number of 0.85 and a Reynolds number of 2 million, based on the length of the fan face diameter, matching an experimental wind tunnel BLI inlet test. The first DOE optimization required a fifth order model having 173 numerical simulation experiments and was able to reduce the DC60 baseline distortion from 64% down to 4.4%, while holding the pressure recovery constant. A second DOE optimization was performed holding the vanes heights at a constant value from the first DOE optimization with the two vane angles-of-attack as design factors. This DOE only required a second order model fit with 15 numerical simulation experiments and reduced DC60 to 3.5% with small decreases in the fourth and fifth harmonic amplitudes. The second optimal vane design was tested at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel in a BLI inlet experiment. The experimental results showed a 80% reduction of DPCPavg, the circumferential distortion level at the engine fan face.

Allan, Brian G.↗
NASA NTRS2007

Extinction-to-Backscatter Ratios of Lofted Aerosol Layers Observed During the First Three Months of CALIPSO Measurements

Case studies from the first three months of the Cloud and Aerosol Lidar and Infrared Pathfinder Spaceborne Observations (CALIPSO) measurements of lofted aerosol layers are analyzed using transmittance [Young, 1995] and two-wavelength algorithms [Vaughan et al., 2004] to determine the aerosol extinction-to-backscatter ratios at 532 and 1064 nm. The transmittance method requires clear air below the layer so that the transmittance through the layer can be determined. Suitable scenes are selected from the browse images and clear air below features is identified by low 532 nm backscatter signal and confirmed by low depolarization and color ratios. The transmittance and two-wavelength techniques are applied to a number of lofted layers and the extinction-to-backscatter ratios are compared with values obtained from the CALIPSO aerosol models [Omar et al., 2004]. The results obtained from these studies are used to adjust the aerosol models and develop observations based extinction-to-backscatter ratio look-up tables and phase functions. Values obtained by these techniques are compared to Sa determinations using other independent methods with a goal of developing probability distribution functions of aerosol type-specific extinction to backscatter ratios. In particular, the results are compared to values determined directly by the High Spectral Resolution Lidar (HSRL) during the CALIPSO CloudSat Validation Experiments (CCVEX) and Sa determined by the application of the two-wavelength lidar Constrained Ratio Aerosol Model-fit (CRAM) retrieval approach [Cattrall et al., 2005; Reagan et al., 2004] to the HSRL data. The results are also compared to values derived using the empirical relationship between the multiple-scattering fraction and the linear depolarization ratio by using Monte Carlo simulations of water clouds [Hu et al., 2006].

Omar, Ali H.↗
NASA NTRS2010

Assessing the Hydrologic Performance of the EPA's Nonpoint Source Water Quality Assessment Decision Support Tool Using North American Land Data Assimilation System (Products)

The accuracy of streamflow predictions in the EPA's BASINS (Better Assessment Science Integrating Point and Nonpoint Sources) decision support tool is affected by the sparse meteorological data contained in BASINS. The North American Land Data Assimilation System (NLDAS) data with high spatial and temporal resolutions provide an alternative to the NOAA National Climatic Data Center (NCDC)'s station data. This study assessed the improvement of streamflow prediction of the Hydrological Simulation Program-FORTRAN (HSPF) model contained within BASINS using the NLDAS 118 degree hourly precipitation and evapotranspiration estimates in seven watersheds of the Chesapeake Bay region. Our results demonstrated consistent improvements of daily streamflow predictions in five of the seven watersheds when NLDAS precipitation and evapotranspiration data was incorporated into BASINS. The improvement of using the NLDAS data is significant when watershed's meteorological station is either far away or not in a similar climatic region. When the station is nearby, using the NLDAS data produces similar results. The correlation coefficients of the analyses using the NLDAS data were greater than 0.8, the Nash-Sutcliffe (NS) model fit efficiency greater than 0.6, and the error in the water balance was less than 5%. Our analyses also showed that the streamflow improvements were mainly contributed by the NLDAS's precipitation data and that the improvement from using NLDAS's evapotranspiration data was not significant; partially due to the constraints of current BASINS-HSPF settings. However, NLDAS's evapotranspiration data did improve the baseflow prediction. This study demonstrates the NLDAS data has the potential to improve stream flow predictions, thus aid the water quality assessment in the EPA nonpoint water quality assessment decision tool.

Lee, S.↗
NASA NTRS2011

The RHESSI Microflare Height Distribution

We present the first in-depth statistical survey of flare source heights observed by RHESSI. Flares were found using a flare-finding algorithm designed to search the 6-10 keV count-rate when RHESSI's full sensitivity was available in order to find the smallest events (Christe et al., 2008). Between March 2002 and March 2007, a total of 25,006 events were found. Source locations were determined in the 4-10 keV, 10-15 keV, and 15-30 keV energy ranges for each event. In order to extract the height distribution from the observed projected source positions, a forward-fit model was developed with an assumed source height distribution where height is measured from the photosphere. We find that the best flare height distribution is given by g (h) oc exp(-h/lambda) where lambda = 6.1 plus or minus 0.3 Mm is the scale height. A power law height distribution with a negative power law index, gamma = 3.1 plus or minus 0.3 is also consistent with the data. Interpreted as thermal loop top sources, these heights are compared to loops generated by a potential field model (PFSS). The measured flare heights distribution are found to be much steeper than the potential field loop height distribution which may be a signature of the flare energization process.

Christe, P.↗
NASA NTRS2011

Gamma-Ray Pulsar Light Curves in Offset Polar Cap Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres, used to model high-energy light curves have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profile~. We find that. corn pared to the profile:-; derived from :-;ymmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines. increases significantly relative to the off-peak emission. formed along leading field lines. The enhanced contrast produces greatly improved slot gap model fits to Fermi pulsar light curves like Vela, which show very little off-peak emIssIon.

Harding, Alice K.↗
NASA NTRS2010

Selection Effects in Identifying Magnetic Clouds and the Importance of the Closest Approach Parameter

This study is motivated by the unusually low number of magnetic clouds (MCs) that are strictly identified within interplanetary coronal mass ejections (ICMEs), as observed at 1 AU; this is usually estimated to be around 30% or lower. But a looser definition of MCs may significantly increase this percentage. Another motivation is the unexpected shape of the occurrence distribution of the observers' "closest approach distances" (measured from a MC's axis, and called CA) which drops off somewhat rapidly as |CA| (in % of MC radius) approaches 100%, based on earlier studies. We suggest, for various geometrical and physical reasons, that the |CA|-distribution should be somewhere between a uniform one and the one actually observed, and therefore the 30% estimate should be higher. So we ask, When there is a failure to identify a MC within an ICME, is it occasionally due to a large |CA| passage, making MC identification more difficult, i.e., is it due to an event selection effect? In attempting to answer this question we examine WIND data to obtain an accurate distribution of the number of MCs vs. |CA| distance, whether the event is ICME-related or not, where initially a large number of cases (N=98) are considered. This gives a frequence distribution that is far from uniform, confirming earlier studies. This along with the fact that there are many ICME identification-parameters that do not depend on |CA| suggest that, indeed an MC event selection effect may explain at least part of the low ratio of (No. MCs)/(No. ICMEs). We also show that there is an acceptable geometrical and physical consistency in the relationships for both average "normalized" magnetic field intensity change and field direction change vs. |CA| within a MC, suggesting that our estimates of |CA|, B(sub 0) (magnetic field intensity on the axis), and choice of a proper "cloud coordinate" system (all needed in the analysis) are acceptably accurate. Therefore the MC fitting model (Lepping et al., 1990) is adequate, on average, for our analysis. However, this selection effect is not likely to completely answer our original question, on the unexpected ratio of MCs to ICMEs, so we must look for other factors, such as peculiarities of CME birth conditions. As a by-product of this analysis, we determine that the first order structural effects within a MC due to its interaction with the solar wind, plus the MC's usual expansion at 1 AU (i.e., the non-force free components of the MC's field) are, on average, weakly dependent on radial distance from the MC's axis; that is, in the outer reaches of a typical MC the non-force free effects show up, but even there they are rather weak. Finally, we show that it is not likely that a MC's size distribution statistically controls the occurrence distribution of the estimated |CA|s.

Lepping, R. P.↗
NASA NTRS2011

Dynamic Black-Level Correction and Artifact Flagging for Kepler Pixel Time Series

Methods applied to the calibration stage of Kepler pipeline data processing [1] (CAL) do not currently use all of the information available to identify and correct several instrument-induced artifacts. These include time-varying crosstalk from the fine guidance sensor (FGS) clock signals, and manifestations of drifting moire pattern as locally correlated nonstationary noise, and rolling bands in the images which find their way into the time series [2], [3]. As the Kepler Mission continues to improve the fidelity of its science data products, we are evaluating the benefits of adding pipeline steps to more completely model and dynamically correct the FGS crosstalk, then use the residuals from these model fits to detect and flag spatial regions and time intervals of strong time-varying black-level which may complicate later processing or lead to misinterpretation of instrument behavior as stellar activity.

Kolodziejczak, J. J.↗
NASA NTRS2011

Gamma-Ray Pulsar Light Curves in Vacuum and Force-Free Geometry

Recent studies have shown that gamma-ray pulsar light curves are very sensitive to the geometry of the pulsar magnetic field. Pulsar magnetic field geometries, such as the retarded vacuum dipole and force-free magnetospheres have distorted polar caps that are offset from the magnetic axis in the direction opposite to rotation. Since this effect is due to the sweepback of field lines near the light cylinder, offset polar caps are a generic property of pulsar magnetospheres and their effects should be included in gamma-ray pulsar light curve modeling. In slot gap models (having two-pole caustic geometry), the offset polar caps cause a strong azimuthal asymmetry of the particle acceleration around the magnetic axis. We have studied the effect of the offset polar caps in both retarded vacuum dipole and force-free geometry on the model high-energy pulse profiles. We find that, compared to the profiles derived from symmetric caps, the flux in the pulse peaks, which are caustics formed along the trailing magnetic field lines, increases significantly relative to the off-peak emission, formed along leading field lines. The enhanced contrast produces improved slot gap model fits to Fermi pulsar light curves like Vela, with vacuum dipole fits being more favorable.

Harding, Alice K.↗
NASA NTRS2011

The WIRED Survey: Infrared Excesses in the SDSS DR7 White Dwarf Catalog - 2

With the launch of the Wide-field Infrar.ed Survey Explorer (WISE), a new era of detecting planetary debris and brown dwarfs (BDs) around white dwarfs (WDs) has begun with the WISE InfraRed Excesses around Degenerates (WIRED) Survey. The WIRED Survey is sensitive to substellar objects and dusty debris around WDs out to distances exceeding 100 pc, well beyond the completeness level of local WDs. In this paper, we present a cross-correlation of the preliminary Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) WD catalog between the WISE, Two-Micron All Sky Survey (2MASS), UKIRT Infrared Deep Sky Survey (UKIDSS), and SDSS DR7 photometric catalogs. From -18,000 input targets, there are WISE detections comprising 344 "naked" WDs (detection of the WD photosphere only), 1020 candidate WD+M dwarf binaries, 42 candidate WD+BD systems, 52 candidate WD+dust disk systems, and 69 targets with indeterminate infrared excess. We classified all of the detected targets through spectral energy distribution model fitting of the merged optical, near-IR, and WISE photometry. Some of these detections could be the result of contaminating sources within the large (approx. 6") WISE point-spread function; we make a preliminary estimate for the rates of contamination for our WD+BD and WD+disk candidates and provide notes for each target of interest. Each candidate presented here should be confirmed with higher angular resolution infrared imaging or infrared spectroscopy. We also present an overview of the observational characteristics of the detected WDs in the WISE photometric bands, including the relative frequencies of candidate WD+M, WD+BD, and WD+disk systems.

Debes, John H.↗
NASA NTRS2014

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x‐acceleration peak response correlated with the human response at 8 and 10‐G 100 ms but not 10‐G 70 ms. The phase lagged the human response. Head z‐acceleration was not correlated. Chest x‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x‐acceleration was not well correlated. Head and chest z‐acceleration was in phase but had a higher peak response. Chest z‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z‐displacement but the THOR was in phase and was comparable to the mean peak response. Head xand z‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, N.↗
NASA NTRS2011

Relation Between the 3D-Geometry of the Coronal Wave and Associated CME During the 26 April 2008 Event

We study the kinematical characteristics and 3D geometry of a large-scale coronal wave that occurred in association with the 26 April 2008 flare-CME event. The wave was observed with the EUVI instruments aboard both STEREO spacecraft (STEREO-A and STEREO-B) with a mean speed of approx 240 km/s. The wave is more pronounced in the eastern propagation direction, and is thus, better observable in STEREO-B images. From STEREO-B observations we derive two separate initiation centers for the wave, and their locations fit with the coronal dimming regions. Assuming a simple geometry of the wave we reconstruct its 3D nature from combined STEREO-A and STEREO-B observations. We find that the wave structure is asymmetric with an inclination toward East. The associated CME has a deprojected speed of approx 750 +/- 50 km/s, and it shows a non-radial outward motion toward the East with respect to the underlying source region location. Applying the forward fitting model developed by Thernisien, Howard, and Vourlidas we derive the CME flux rope position on the solar surface to be close to the dimming regions. We conclude that the expanding flanks of the CME most likely drive and shape the coronal wave.

CMEs↗
NASA NTRS2013

Daily MODIS 500 m Reflectance Anisotropy Direct Broadcast (DB) Products for Monitoring Vegetation Phenology Dynamics

Land surface vegetation phenology is an efficient bio-indicator for monitoring ecosystem variation in response to changes in climatic factors. The primary objective of the current article is to examine the utility of the daily MODIS 500 m reflectance anisotropy direct broadcast (DB) product for monitoring the evolution of vegetation phenological trends over selected crop, orchard, and forest regions. Although numerous model-fitted satellite data have been widely used to assess the spatio-temporal distribution of land surface phenological patterns to understand phenological process and phenomena, current efforts to investigate the details of phenological trends, especially for natural phenological variations that occur on short time scales, are less well served by remote sensing challenges and lack of anisotropy correction in satellite data sources. The daily MODIS 500 m reflectance anisotropy product is employed to retrieve daily vegetation indices (VI) of a 1 year period for an almond orchard in California and for a winter wheat field in northeast China, as well as a 2 year period for a deciduous forest region in New Hampshire, USA. Compared with the ground records from these regions, the VI trajectories derived from the cloud-free and atmospherically corrected MODIS Nadir BRDF (bidirectional reflectance distribution function) adjusted reflectance (NBAR) capture not only the detailed footprint and principal attributes of the phenological events (such as flowering and blooming) but also the substantial inter-annual variability. This study demonstrates the utility of the daily 500 m MODIS reflectance anisotropy DB product to provide daily VI for monitoring and detecting changes of the natural vegetation phenology as exemplified by study regions comprising winter wheat, almond trees, and deciduous forest.

Vegetation↗
NASA NTRS2014

High-Resolution Submillimeter and Near-Infrared Studies of the Transition Disk Around Sz 91

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(3--2) observations with the Submillimeter Array (approximately 1" - 3" resolution), and high-resolution imaging of polarized intensity at the K(sub s) -band by using the HiCIAO instrument on the Subaru Telescope (0.25" resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 AU and 170 AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H 2 mass of 2.4×10(exp −3) M(solar mass) in the cold (T less than 30 K) outer part at 65 less than r less than 170 AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount (greater than 3×10(exp −9) M(solar mass)) of hot (T approximately 180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(3--2) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

canonical gas-to-dust mass ratio↗
NASA NTRS2014

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post-­‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi-­‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two-­‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x-­‐acceleration peak response correlated with the human response at 8 and 10-­‐G 100 ms but not 10-­‐G 70 ms. The phase lagged the human response. Head z-­‐acceleration was not correlated. Chest x-­‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x-­‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x-­‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x-­‐acceleration was not well correlated. Head and chest z-­‐acceleration was in phase but had a higher peak response. Chest z-­‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z-­‐displacement but the THOR was in phase and was comparable to the mean peak response. Head x-­‐ and z-­‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, Nate↗
NASA NTRS2014

High-Resolution Submillimeter and Near-Infrared Studies of the Transition Disk around Sz 91

To reveal the structures of a transition disk around a young stellar object in Lupus, Sz 91, we have performed aperture synthesis 345 GHz continuum and CO(32) observations with the Submillimeter Array ( 13 resolution), and high-resolution imaging of polarized intensity at the Ks-band by using the Hi-CIAO instrument on the Subaru Telescope (0.25 resolution). Our observations successfully resolved the inner and outer radii of the dust disk to be 65 and 170AU, respectively, which indicates that Sz 91 is a transition disk source with one of the largest known inner holes. The model fitting analysis of the spectral energy distribution reveals an H2 mass of 2.4 103 M in the cold (T 30 K) outer part at 65 r 170 AU by assuming a canonical gas-to-dust mass ratio of 100, although a small amount ( 3109 M) of hot (T 180 K) dust possibly remains inside the inner hole of the disk. The structure of the hot component could be interpreted as either an unresolved self-luminous companion body (not directly detected in our observations) or a narrow ring inside the inner hole. Significant CO(32) emission with a velocity gradient along the major axis of the dust disk is concentrated on the Sz 91 position, suggesting a rotating gas disk with a radius of 420 AU. The Sz 91 disk is possibly a rare disk in an evolutionary stage immediately after the formation of protoplanets because of the large inner hole and the lower disk mass than other transition disks studied thus far.

Hi-CIAO instrument↗
NASA NTRS2013

Studies in Astronomical Time Series Analysis. VI. Bayesian Block Representations

This paper addresses the problem of detecting and characterizing local variability in time series and other forms of sequential data. The goal is to identify and characterize statistically significant variations, at the same time suppressing the inevitable corrupting observational errors. We present a simple nonparametric modeling technique and an algorithm implementing it-an improved and generalized version of Bayesian Blocks [Scargle 1998]-that finds the optimal segmentation of the data in the observation interval. The structure of the algorithm allows it to be used in either a real-time trigger mode, or a retrospective mode. Maximum likelihood or marginal posterior functions to measure model fitness are presented for events, binned counts, and measurements at arbitrary times with known error distributions. Problems addressed include those connected with data gaps, variable exposure, extension to piece- wise linear and piecewise exponential representations, multivariate time series data, analysis of variance, data on the circle, other data modes, and dispersed data. Simulations provide evidence that the detection efficiency for weak signals is close to a theoretical asymptotic limit derived by [Arias-Castro, Donoho and Huo 2003]. In the spirit of Reproducible Research [Donoho et al. (2008)] all of the code and data necessary to reproduce all of the figures in this paper are included as auxiliary material.

signal detection↗
NASA NTRS2015

Mantle Water Contents Beneath the Rio Grande Rift (NM, USA): FTIR Analysis of Rio Puerco and Kilbourne Hole Peridotite Xenoliths

Peridotite xenoliths from the Rio Grande Rift (RGR) are being analyzed for H (sub 2) O contents by FTIR (Fourier Transform Infrared) as well as for major and trace element compositions. Nine samples are from the Rio Puerco Volcanic Field (RP) which overlaps the central RGR and southeastern Colorado Plateau; seventeen samples are from Kilbourne Hole (KH) in the southern RGR. Spinel Cr# (Cr/(Cr+Al)) (0.08-0.46) and olivine Mg# (Mg/(Mg plus Fe)) (0.883-0.911) of all RGR samples fall within the olivine-spinel mantle array from [1], an indicator that peridotites are residues of partial melting. Pyroxene H (sub 2) O in KH correlate with bulk rock and pyroxene Al (sub 2) O (sub 3).The KH clinopyroxene rare earth element (REE) variations fit models of 0-13 percent fractional melting of a primitive upper mantle. Most KH peridotites have bulk-rock light REE depleted patterns, but five are enriched in light REEs consistent with metasomatism. Variation in H (sub 2) O content is unrelated to REE enrichment. Metasomatism is seen in RP pyroxenite xenoliths [2] and will be examined in the peridotites studied here. Olivine H (sub 2) O contents are low (less than or equal to 15 parts per million), and decrease from core to rim within grains. This is likely due to H loss during xenolith transport by the host magma [3]. Diffusion models of H suggest that mantle H (sub 2) O contents are still preserved in cores of KH olivine, but not RP olivine. The average H (sub 2) O content of Colorado Plateau clinopyroxene (670 parts per million) [4] is approximately 300 parts per million higher than RGR clinopyroxene (350 parts per million). This upholds the hypothesis that hydration-induced lithospheric melting occurred during flat-slab subduction of the Farallon plate [5]. Numerical models indicate hydration via slab fluids is possible beneath the plateau, approximately 600 kilometers from the paleo-trench, but less likely approximately 850 kilometers away beneath the rift [6].

Schaffer, L. A.↗
NASA NTRS2015

Modeling the Swift BAT Trigger Algorithm with Machine Learning

To draw inferences about gamma-ray burst (GRB) source populations based on Swift observations, it is essential to understand the detection efficiency of the Swift burst alert telescope (BAT). This study considers the problem of modeling the Swift BAT triggering algorithm for long GRBs, a computationally expensive procedure, and models it using machine learning algorithms. A large sample of simulated GRBs from Lien et al. (2014) is used to train various models: random forests, boosted decision trees (with AdaBoost), support vector machines, and artificial neural networks. The best models have accuracies of approximately greater than 97% (approximately less than 3% error), which is a significant improvement on a cut in GRB flux which has an accuracy of 89:6% (10:4% error). These models are then used to measure the detection efficiency of Swift as a function of redshift z, which is used to perform Bayesian parameter estimation on the GRB rate distribution. We find a local GRB rate density of eta(sub 0) approximately 0.48(+0.41/-0.23) Gpc(exp -3) yr(exp -1) with power-law indices of eta(sub 1) approximately 1.7(+0.6/-0.5) and eta(sub 2) approximately -5.9(+5.7/-0.1) for GRBs above and below a break point of z(sub 1) approximately 6.8(+2.8/-3.2). This methodology is able to improve upon earlier studies by more accurately modeling Swift detection and using this for fully Bayesian model fitting. The code used in this is analysis is publicly available online.

gamma rays: general↗