Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained variational analysis”

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 325 records · Page 18

Deformable segmentation of 3D MR prostate images via distributed discriminative dictionary and ensemble learning

Purpose: Automatic prostate segmentation from MR images is an important task in various clinical applications such as prostate cancer staging and MR-guided radiotherapy planning. However, the large appearance and shape variations of the prostate in MR images make the segmentation problem difficult to solve. Traditional Active Shape/Appearance Model (ASM/AAM) has limited accuracy on this problem, since its basic assumption, i.e., both shape and appearance of the targeted organ follow Gaussian distributions, is invalid in prostate MR images. To this end, the authors propose a sparse dictionary learning method to model the image appearance in a nonparametric fashion and further integrate the appearance model into a deformable segmentation framework for prostate MR segmentation. Methods: To drive the deformable model for prostate segmentation, the authors propose nonparametric appearance and shape models. The nonparametric appearance model is based on a novel dictionary learning method, namely distributed discriminative dictionary (DDD) learning, which is able to capture fine distinctions in image appearance. To increase the differential power of traditional dictionary-based classification methods, the authors' DDD learning approach takes three strategies. First, two dictionaries for prostate and nonprostate tissues are built, respectively, using the discriminative features obtained from minimum redundancy maximum relevance feature selection. Second, linear discriminant analysis is employed as a linear classifier to boost the optimal separation between prostate and nonprostate tissues, based on the representation residuals from sparse representation. Third, to enhance the robustness of the authors' classification method, multiple local dictionaries are learned for local regions along the prostate boundary (each with small appearance variations), instead of learning one global classifier for the entire prostate. These discriminative dictionaries are located on different patches of the prostate surface and trained to adaptively capture the appearance in different prostate zones, thus achieving better local tissue differentiation. For each local region, multiple classifiers are trained based on the randomly selected samples and finally assembled by a specific fusion method. In addition to this nonparametric appearance model, a prostate shape model is learned from the shape statistics using a novel approach, sparse shape composition, which can model nonGaussian distributions of shape variation and regularize the 3D mesh deformation by constraining it within the observed shape subspace. Results: The proposed method has been evaluated on two datasets consisting of T2-weighted MR prostate images. For the first (internal) dataset, the classification effectiveness of the authors' improved dictionary learning has been validated by comparing it with three other variants of traditional dictionary learning methods. The experimental results show that the authors' method yields a Dice Ratio of 89.1% compared to the manual segmentation, which is more accurate than the three state-of-the-art MR prostate segmentation methods under comparison. For the second dataset, the MICCAI 2012 challenge dataset, the authors' proposed method yields a Dice Ratio of 87.4%, which also achieves better segmentation accuracy than other methods under comparison. Conclusions: A new magnetic resonance image prostate segmentation method is proposed based on the combination of deformable model and dictionary learning methods, which achieves more accurate segmentation performance on prostate T2 MR images.

60 APPLIED LIFE SCIENCES↗

The reverberation lag in the low-mass X-ray binary H1743-322

The evolution of the inner accretion flow of a black hole X-ray binary during an outburst is still a matter of active research. X-ray reverberation lags are powerful tools for constraining disk–corona geometry. We present a study of X-ray lags in the black hole transient H1743-322. We compared the results obtained from analysis of all the publicly available XMM-Newton observations. These observations were carried out during two different outbursts that occurred in 2008 and 2014. During all the observations the source was caught in the hard state and at similar luminosities (L{sub 3−10keV}/L{sub Edd}∼0.004). We detected a soft X-ray lag of ∼60 ms, most likely due to thermal reverberation. We did not detect any significant change of the lag amplitude among the different observations, indicating a similar disk–corona geometry at the same luminosity in the hard state. On the other hand, we observe significant differences between the reverberation lag detected in H1743-322 and in GX 339-4 (at similar luminosities in the hard state), which might indicate variations of the geometry from source to source.

79 ASTRONOMY AND ASTROPHYSICS↗

Predicting the mineral composition of dust aerosols: Insights from elemental composition measured at the Izaña Observatory: PREDICTING DUST MINERAL COMPOSITION

Regional variations of dust mineral composition are fundamental to climate impacts but generally neglected in climate models. A challenge for models is that atlases of soil composition are derived from measurements following wet sieving, which destroys the aggregates potentially emitted from the soil. Aggregates are crucial to simulating the observed size distribution of emitted soil particles. We use an extension of brittle fragmentation theory in a global dust model to account for these aggregates. Our method reproduces the size-resolved dust concentration along with the approximately size-invariant fractional abundance of elements like Fe and Al in the decade-long aerosol record from the Izaña Observatory, off the coast of West Africa. By distinguishing between Fe in structural and free forms, we can attribute improved model behavior to aggregation of Fe and Al-rich clay particles. We also demonstrate the importance of size-resolved measurements along with elemental composition analysis to constrain models.

58 GEOSCIENCES↗

The clustering of the SDSS-IV extended Baryon Oscillation Spectroscopic Survey quasar sample: testing observational systematics on the Baryon Acoustic Oscillation measurement

ABSTRACT Baryon Acoustic Oscillations (BAO) are considered to be a very robust standard ruler against various systematics. This premise has been tested against observational systematics, but not to the level required for the next generation of galaxy surveys such as the Dark Energy Spectroscopic Instrument (DESI) and Euclid. In this paper, we investigate the effect of observational systematics on the BAO measurement of the final sample of quasars from the extended Baryon Oscillation Spectroscopic Survey Data Release 16 in order to prepare and hone a similar analysis for upcoming surveys. We employ catalogues with various treatments of imaging systematic effects using linear and neural network-based non-linear approaches and consider how the BAO measurement changes. We also test how the variations to the BAO fitting model respond to the observational systematics. As expected, we confirm that the BAO measurements obtained from the DR16 quasar sample are robust against imaging systematics well within the statistical error, while reporting slightly modified constraints that shift the line-of-sight BAO signal by less than 1.1 per cent. We use realistic simulations with similar redshift and angular distributions as the DR16 sample to conduct statistical tests for validating the pipeline, quantifying the significance of differences, and estimating the expected bias on the BAO scale in future high-precision data sets. Although we find a marginal impact for the eBOSS QSO data, the work presented here is of vital importance for constraining the nature of dark energy with the BAO feature in the new era of big data cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗

Periodic Solar Wind Structures Observed in Measurements of Elemental and Ionic Composition in situ at L1

Mesoscale periodic structures observed in solar wind plasma serve as an important diagnostic tool for constraining the processes that govern the formation of the solar wind. These structures have been observed in situ and in remote data as fluctuations in proton and electron density. However, only two events of this type have been reported regarding the elemental and ionic composition. Composition measurements are especially important in gaining an understanding of the origin of the solar wind as the composition is frozen into the plasma at the Sun and does not evolve as it advects through the heliosphere. Here, we present the analysis of four events containing mesoscale periodic solar wind structure during which the Iron and Magnesium number density data, measured by the Solar Wind Ion Composition Spectrometer (SWICS)on board the Advanced Composition Explorer spacecraft, are validated at statistically significant count levels. We use a spectral analysis method specifically designed to extract periodic signals from astrophysical time series and apply it to the SWICS 12-minute native resolution dataset. We find variations in the relative abundance of elements with low first ionization potential, mass dependencies, and charge state during time intervals in which mesoscale periodic structures are observed. These variations are linked to temporal or spatial variations in solar source regions and put constraints on the solar wind formation mechanisms that produce them. Techniques presented here are relevant for future, higher-resolution studies of data from new instruments such as Solar Orbiter’s Heavy Ion Sensor.

Irena Gershkovich↗

Azimuthal Size Scales of Solar Wind Periodic Density Structures

Periodic Density Structures (PDSs) are quasi-periodic variations of solar wind density ranging from a few minutes to a few hours. PDSs advect with the solar wind and have radial length scales (Lx) of tens to several thousand Mm, thus belonging to the class of “mesoscale structures”. Current interplanetary multi-spacecraft observations are not at spatial separations capable of directly measuring the 3D size scale of PDSs or other mesoscale structures. Instead, previous investigations estimated characteristic spatial scales in solar wind parameters using cross-correlation and/or coherence analysis applied to multi-spacecraft observations. For solar wind density and interplanetary magnetic field (IMF) intensity, the reported size scales perpendicular to the Sun-Earth line (Ly) ranged between ≈30 and ≈200 Earth Radii (RE). Here, we implemented a similar approach on the same parameters but focused on high density slow solar wind intervals with PDSs observed by the Wind and ARTEMIS-P1 spacecraft. Additionally, this is the first statistical study of the IMF intensity periodicities in relation to PDSs. We identified intervals in which the two spacecraft observed the same periodicity obtaining two PDS groups based on their radial length scale: Lx1 ≈86RE and Lx2 ≈35RE. Then, we classified the events based on the periodic variations’ coherence level. Reproducing the results with simulations of the PDSs’ transit, we inferred the Ly order of magnitudes for the two PDS groups: Ly1 ≈340RE and Ly2 ≈187RE. Knowing the PDSs’ size scales is fundamental to constrain models aimed at reproducing these structures and is critical to better understand the PDS-magnetosphere coupling.

Interplanetary physics↗

AGE AND MASS FOR 920 LARGE MAGELLANIC CLOUD CLUSTERS DERIVED FROM 100 MILLION MONTE CARLO SIMULATIONS

We present new age and mass estimates for 920 stellar clusters in the Large Magellanic Cloud (LMC) based on previously published broadband photometry and the stellar cluster analysis package, MASSCLEANage. Expressed in the generic fitting formula, d{sup 2} N/dMdt{proportional_to}M{sup {alpha}} t{sup {beta}}, the distribution of observed clusters is described by {alpha} = -1.5 to -1.6 and {beta} = -2.1 to -2.2. For 288 of these clusters, ages have recently been determined based on stellar photometric color-magnitude diagrams, allowing us to gauge the confidence of our ages. The results look very promising, opening up the possibility that this sample of 920 clusters, with reliable and consistent age, mass, and photometric measures, might be used to constrain important characteristics about the stellar cluster population in the LMC. We also investigate a traditional age determination method that uses a {chi}{sup 2} minimization routine to fit observed cluster colors to standard infinite-mass limit simple stellar population models. This reveals serious defects in the derived cluster age distribution using this method. The traditional {chi}{sup 2} minimization method, due to the variation of U, B, V, R colors, will always produce an overdensity of younger and older clusters, with an underdensity of clusters in the log (age/yr) = [7.0, 7.5] range. Finally, we present a unique simulation aimed at illustrating and constraining the fading limit in observed cluster distributions that includes the complex effects of stochastic variations in the observed properties of stellar clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

Slip and Dilation Tendency Analysis of the Tuscarora Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the Tuscarora geothermal field was calculated based on the faults mapped Tuscarora area (Dering, 2013). The Tuscarora area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the Tuscarora area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Tuscarora is defined by a left-step in a major north- to-north northeast striking, west-dipping range-bounding normal fault system. Faults within the broad step define an anticlinal accommodation ...

15 GEOTHERMAL ENERGY↗

Interpreting mixing relationships in energetic melts to estimate vapor contribution and composition

Heterogeneous, glassy debris can be produced by energetic events such as near surface nuclear explosions and asteroid impacts. In these events, the resultant debris comprises a rapidly-cooled mixture of the projectile or bomb material with proximate materials such as soil and anthropogenic structures. Glasses formed in these events preserve inter and/or intra sample compositional variation reflecting complex mixtures of the precursor materials or their chemically fractionated products. Understanding formation mechanisms and chemical processes relating to energetic glass formation based on study of the end products requires an objective method for calculating the precursor components that mixed together to produce the final product. Here, we compare a set of multivariate linear-least squares approaches to address the unmixing problem of heterogeneous glasses including classical least squares (CLS), principal component analysis (PCA) and two multivariate curve resolution – alternating least squares (MCR-ALS) approaches. We generate a synthetic data set representing spatially resolved, compositional data in a set of heterogeneously mixed glasses to understand model sensitivity to precursor abundance and composition. We explore the impact of incomplete knowledge of the system prior to the formation of the melts on the resulting confidence and accuracy of the solution (i.e., the number, composition, and abundances of the precursor materials). We conclude that a closure and equality constrained MCR-ALS approach succeeds in identifying unrecognized precursor components and is a suitable alternative to CLS approaches when the investigator has incomplete knowledge of the starting precursor chemical compositions. Finally, we demonstrate application of this MCR-ALS approach to measurements of Trinitite.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Ground-based detectability of terrestrial and Jovian extrasolar planets: observations of CM Draconis at Lick Observatory

The detection of terrestrial-sized extrasolar planets from the ground has been thought to be virtually impossible due to atmospheric scintillation limits. However, we show that this is not the case especially selected (but nevertheless main sequence) stars, namely small eclipsing binaries. For the smallest of these systems, CM Draconis, several months to a few years of photometric observations with 1-m-class telescopes will be sufficient to detect the transits of any short-period planets of sizes > or = 1.5 Earth radii (RE), using cross-correlation analysis with moderately good photometry. Somewhat larger telescopes will be needed to extend this detectability to terrestrial planets in larger eclipsing binary systems. (We arbitrarily define "terrestrial planets" herein as those whose disc areas are closer to that of Earth's than Neptune's i.e., less than about 2.78 RE.) As a "spin-off" of such observations, we will also be able to detect the presence of Jovian-mass planets without transits using the timing of the eclipse minima. Eclipse minima will drift in time as the binary system is offset by a sufficiently massive planet (i.e., one Jupiter mass) about the binary/giant-planet barycenter, causing a periodic variation in the light travel time to the observer. We present here an outline of present observations taking place at the University of California Lick Observatory using the Crossley 0.9-m telescope in collaboration with other observatories (in South Korea, Crete, France, Canary Islands, and New York) to detect or constrain the existence of terrestrial planets around main sequence eclipsing binary star systems, starting with CM Draconis. We demonstrate the applicability of photometric data to the general detection of gas giant planets via eclipse minima timings in many other small-mass eclipsing binary systems as well.

NASA Program Exobiology↗

Slip and Dilation Tendency Analysis of the San Emidio Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the San Emidio geothermal field was calculated based on the faults mapped Tuscarora area (Rhodes, 2011). The San Emidio area lies in the Basin and Range Province, as such we applied a normal faulting stress regime to the San Emidio area faults, with a minimum horizontal stress direction oriented 115, based on inspection of local and regional stress determinations, as explained above. This is consistent with the shmin determined through inversion of fault data by Rhodes (2011). Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Interesting, the San Emidio geothermal field lies in an area of primarily north striking faults, w...

15 GEOTHERMAL ENERGY↗

Slip and Dilation Tendency Analysis of the Salt Wells Geothermal Area

Critically stressed fault segments have a relatively high likelihood of acting as fluid flow conduits (Sibson, 1994). As such, the tendency of a fault segment to slip (slip tendency; Ts; Morris et al., 1996) or to dilate (dilation tendency; Td; Ferrill et al., 1999) provides an indication of which faults or fault segments within a geothermal system are critically stressed and therefore likely to transmit geothermal fluids. The slip tendency of a surface is defined by the ratio of shear stress to normal stress on that surface: Ts = T / on (Morris et al., 1996). Dilation tendency is defined by the stress acting normal to a given surface: Td = (o1-on) / (o1-o3) (Ferrill et al., 1999). Slip and dilation were calculated using 3DStress (Southwest Research Institute). Slip and dilation tendency are both unitless ratios of the resolved stresses applied to the fault plane by ambient stress conditions. Values range from a maximum of 1, a fault plane ideally oriented to slip or dilate under ambient stress conditions to zero, a fault plane with no potential to slip or dilate. Slip and dilation tendency values were calculated for each fault in the focus study areas at, McGinness Hills, Neal Hot Springs, Patua, Salt Wells, San Emidio, and Tuscarora on fault traces. As dip is not well constrained or unknown for many faults mapped in within these we made these calculations using the dip for each fault that would yield the maximum slip tendency or dilation tendency. As such, these results should be viewed as maximum tendency of each fault to slip or dilate. The resulting along-fault and fault-to-fault variation in slip or dilation potential is a proxy for along fault and fault-to-fault variation in fluid flow conduit potential. Stress Magnitudes and directions Stress field variation within each focus area was approximated based on regional published data and the world stress database (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2010; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012; Moeck et al., 2010; Moos and Ronne, 2010 and Reinecker et al., 2005) as well as local stress information if applicable. For faults within these focus systems we applied either a normal faulting stress regime where the vertical stress (sv) is larger than the maximum horizontal stress (shmax) which is larger than the minimum horizontal stress (sv>shmax>shmin) or strike-slip faulting stress regime where the maximum horizontal stress (shmax) is larger than the vertical stress (sv) which is larger than the minimum horizontal stress (shmax >sv>shmin) depending on the general tectonic province of the system. Based on visual inspection of the limited stress magnitude data in the Great Basin we used magnitudes such that shmin/shmax = .527 and shmin/sv= .46, which are consistent with complete and partial stress field determinations from Desert Peak, Coso, the Fallon area and Dixie valley (Hickman et al., 2000; Hickman et al., 1998 Robertson-Tait et al., 2004; Hickman and Davatzes, 2011; Davatzes and Hickman, 2006; Blake and Davatzes 2011; Blake and Davatzes, 2012). Slip and dilation tendency for the Salt Wells geothermal field was calculated based on the faults mapped in the Bunejug Mountains quadrangle (Hinz et al., 2011). The Salt Wells area lies in the Basin and Range Province (N. Hinz personal comm.) As such we applied a normal faulting stress regime to the Salt Wells area faults, with a minimum horizontal stress direction oriented 105, based on inspection of local and regional stress determinations. Under these stress conditions north-northeast striking, steeply dipping fault segments have the highest dilation tendency, while north-northeast striking 60 degrees dipping fault segments have the highest tendency to slip. Several such faults intersect in high density in the core of the accommodation zone in the Bunejug Mountains and local to the Salt Wells geothermal .

15 GEOTHERMAL ENERGY↗

Neutron Absorption Measurements Constrain Eucrite-Diogenite Mixing in Vesta's Regolith

The NASA Dawn Mission s Gamma Ray and Neutron Detector (GRaND) [1] acquired mapping data during 5 months in a polar, low altitude mapping orbit (LAMO) with approx.460-km mean radius around main-belt asteroid Vesta (264-km mean radius) [2]. Neutrons and gamma rays are produced by galactic cosmic ray interactions and by the decay of natural radioelements (K, Th, U), providing information about the elemental composition of Vesta s regolith to depths of a few decimeters beneath the surface. From the data acquired in LAMO, maps of vestan neutron and gamma ray signatures were determined with a spatial resolution of approx.300 km full-width-at-half-maximum (FWHM), comparable in scale to the Rheasilvia impact basin (approx.500 km diameter). The data from Vesta encounter are available from the NASA Planetary Data System. Based on an analysis of gamma-ray spectra, Vesta s global-average regolith composition was found to be consistent with the Howardite, Eucrite, and Diogenite (HED) meteorites, reinforcing the HED-Vesta connection [2-7]. Further, an analysis of epithermal neutrons revealed variations in the abundance of hydrogen on Vesta s surface, reaching values up to 400 micro-g/g [2]. The association of high concentrations of hydrogen with equatorial, low-albedo surface regions indicated exogenic delivery of hydrogen by the infall of carbonaceous chondrite (CC) materials. This finding was buttressed by the presence of minimally-altered CC clasts in howardites, with inferred bulk hydrogen abundances similar to that found by GRaND, and by studies using data from Dawn s Framing Camera (FC) and VIR instruments [8-10]. In addition, from an analysis of neutron absorption, spatial-variations in the abundance of elements other than hydrogen were detected [2].

Prettyman, T. H.↗

Influence of physical and chemical hydrology on bioremediation of a U-contaminated aquifer informed by reactive transport modeling incorporating 238 U/ 235 U ratios

Microbially-catalyzed reductive immobilization of aqueous uranium (U) as a solid phase has been proposed as a U remediation technique. Both laboratory and field experiments have demonstrated that this reduction reaction alters the 238 U/ 235 U ratio, producing a 238 U-enriched U(IV) solid. In contrast, other major U reactive transport processes fractionate these isotopes much less. This suggests the potential to quantify the extent of bioreduction occurring in groundwater containing U using the 238 U/ 235 U ratio, which would substantially improve upon current practices largely relying on U concentration measurements alone. Current reactive transport models for uranium dynamics include only concentration measurements, which are strongly influenced by highly coupled reactive transformation and aqueous transport processes. The complex physical and chemical behavior of U potentially compromises such quantitative analysis of U storage and release. Here we report the first numerical reactive transport model which explicitly incorporates variations in the 238 U/ 235 U ratio of U and demonstrates improved interpretation of the principal chemical reactions and groundwater transport processes affecting the subsurface mobility and distribution of this widespread contaminant. Recent U bioreduction studies performed in a contaminated aquifer in Rifle, Colorado, USA applied Rayleigh distillation models to interpret U stable isotope fractionation observed as a result of acetate amendment. These simplified models were unable to resolve the spatiotemporal pattern of U isotope fractionation recorded in the aqueous solutes. Here, we employ the multi-component, isotope-enabled CrunchTope reactive transport software to interpret these measured U isotope ratios, and demonstrate accurate reproduction of observed trends in both geochemistry and 238 U/ 235 U ratios for two consecutive years of field experiments. Overall, our results indicate that accurately modeling both the U concentration and isotope ratio distributions greatly constrains the parameter space of the model. We find that the transport properties of U in the Rifle aquifer are governed by the presence of low-permeability regions, which the isotopes are uniquely sensitive to. When U reduction is spatially constrained by these low-permeability regions, the shift in the 238 U/ 235 U ratio becomes more muted. Accurate modeling of observed U isotope ratios thus provides a powerful means to better understand bioremediation, and the current study serves to advance the application of this novel method.

238U/235U↗

Large reductions in Permian Basin methane intensity shown in multi-year comparison of aerially-visible methane emissions

Spanning the US states of Texas and New Mexico, the Permian Basin has been a hotspot of methane emissions from oil and natural gas activity 1–4, although studies disagree over the magnitude of these emissions. The most comprehensive measurement campaigns published were conducted in 2019 1–5 .There have been large changes in the energy industry since then, including in the prices of oil and gas, both state and federal regulatory environments, investor and activist pressure over methane emissions, and the adoption of new technologies and policies by energy operators. Understanding how any or all of these might influence methane emissions is important for policy makers, oil and gas operators, and other stakeholders. We characterize the time evolution of Permian Basin methane emissions using a series of comprehensive aerial surveys conducted every year from 2020-2023 and compare them to the 2019 results cited above. To maintain comparability, all the data sets are from surveys using Insight M point source methane sensing technology. The scope of these surveys expanded over time: from 33-46% of wells, oil production, and gas production in 2020 to 60-65% in 2021, to 84% of wells and over 90% of both oil and gas production in 2023. These surveys by Insight M also include hundreds of gas processing plants and compressor stations as well as 1000s of km of gathering and transmission pipelines. Considering only the aerially detected portion of emissions (typically the majority of the total in such surveys 4), we find reductions of more than 70% in methane emissions intensity compared to the 2019 New Mexico-only Insight M survey, with variation depending on the year 3,4. Notably, although sources below 100 kg/hr contributed less than 10% of aerially measured emissions the 2019 New Mexico survey 3, these smaller sources constitute a larger proportion of total aerially measured emissions (although not the majority) in 2020-2023. Production facilities and gathering pipelines are responsible for the larges shares of total emissions, followed by compressor stations and gas processing plants. Permian methane emissions were also measured in a comprehensive 2019 Permian-wide survey by the Carbon Mapper team 2. That analysis led to a lower total emissions estimate at the time 4. These new Insight M-based emission rates are still roughly 30-70% lower than the aerially measured portion of the 2019 Carbon Mapper-based estimates 4. Further work is needed to harmonize these surveys in space and time to create the most intercomparable numbers possible 5. Additional analysis is needed to compare our findings to the more spatially constrained 2020, 2021, and 2023 Carbon Mapper surveys in the Permian 4,6. The evidence is strong from these two survey teams that emissions intensity has declined significantly since 2019. Reasons for this trend are currently unclear but point to possible success of emissions control programs. Future work investigating frequency, source, and operator-specific intensities could provide insights into the causes of this promising trend.

methane, oil and gas, data science, remote sensing↗

HOD-dependent systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we estimate the level of systematics induced in the DESI BAO analysis due the assumed Halo Occupation Distribution (HOD) model for the Emission Line Galaxy (ELG) tracer. We make use of mock galaxy catalogs constructed by fitting various HOD models to early DESI data, namely the One-Percent survey data. Our analysis includes typical HOD models for the ELG tracer used in the literature as well as extensions to the baseline models. Among the extensions, we consider various recipes for galactic conformity and assembly bias. We use 25 AbacusSummit simulations under the ΛCDM cosmology for each HOD model and perform independent analyses in Fourier space and in configuration space. To recover the BAO signal from our mocks we perform BAO reconstruction and apply the control variates technique to reduce sample variance noise. Our BAO analyses can recover the isotropic BAO parameter α iso within 0.1% and the Alcock Paczynski parameter α AP within 0.3%. Overall, we find that the systematic error due to the HOD dependence is below 0.17%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis, for the assumptions made.

79 ASTRONOMY AND ASTROPHYSICS↗

MODELING GALACTIC EXTINCTION WITH DUST AND 'REAL' POLYCYCLIC AROMATIC HYDROCARBONS

We investigate the remarkable apparent variety of galactic extinction curves by modeling extinction profiles with core-mantle grains and a collection of single polycyclic aromatic hydrocarbons. Our aim is to translate a synthetic description of dust into physically well-grounded building blocks through the analysis of a statistically relevant sample of different extinction curves. All different flavors of observed extinction curves, ranging from the average galactic extinction curve to virtually 'bumpless' profiles, can be described by the present model. We prove that a mixture of a relatively small number (54 species in 4 charge states each) of polycyclic aromatic hydrocarbons can reproduce the features of the extinction curve in the ultraviolet, dismissing an old objection to the contribution of polycyclic aromatic hydrocarbons to the interstellar extinction curve. Despite the large number of free parameters (at most the 54 Multiplication-Sign 4 column densities of each species in each ionization state included in the molecular ensemble plus the 9 parameters defining the physical properties of classical particles), we can strongly constrain some physically relevant properties such as the total number of C atoms in all species and the mean charge of the mixture. Such properties are found to be largely independent of the adopted dust model whose variation provides effects that are orthogonal to those brought about by the molecular component. Finally, the fitting procedure, together with some physical sense, suggests (but does not require) the presence of an additional component of chemically different very small carbonaceous grains.

79 ASTRONOMY AND ASTROPHYSICS↗