Search NASA⌕ Search

SEARCH · Search NASA

Results for “R&T Report 2016”

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.

63 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Constraining Remote Oxidation Capacity with ATom Observations

The global oxidation capacity, defined as the tropospheric mean concentration of the hydroxyl radical (OH), controls the lifetime of reactive trace gases in the atmosphere such as methane and carbon monoxide (CO). Models tend to underestimate the methane lifetime and CO concentrations throughout the troposphere, which is consistent with excessive OH. Approximately half the oxidation of methane and non-methane volatile organic compounds (VOCs) is thought to occur over the oceans where oxidant chemistry has received little validation due to a lack of observational constraints. We use observations from the first two deployments of the NASA ATom aircraft campaign during July–August 2016 and January–February 2017 to evaluate the oxidation capacity over the remote oceans and its representation in the GEOS-Chem chemical transport model. The model successfully simulates the magnitude and vertical profile of remote OH within the measurement uncertainties. Comparisons against the drivers of OH production (water vapor, ozone, and NOy concentrations, ozone photolysis frequencies) also show minimal bias with the exception of wintertime NOy, for which a model overestimate may indicate insufficient wet scavenging and/or missing loss on seasalt aerosol but large uncertainties remain that require further studies of NOy partitioning and removal in the troposphere. During the ATom-1 deployment, OH reactivity (OHR) below 3 km is significantly enhanced, and this is not captured by the sum of its measured components (cOHRobs) or by the model (cOHRmod). This enhancement could suggest missing reactive VOCs but cannot be explained by new estimates of ocean VOC sources and additional modeled reactivity in this region would be difficult to reconcile with the full suite of ATom measurement constraints. The model generally reproduces the magnitude and seasonality of cOHRobs but underestimates the contribution of oxygenated VOC, mainly acetaldehyde, which is severely underestimated throughout the troposphere despite its calculated lifetime of less than a day. Missing model acetaldehyde in previous studies was attributed to measurement uncertainties that have been largely resolved. Observations of peroxyacetic acid (PAA) provide new support for remote levels of acetaldehyde. The underestimate in modeled acetaldehyde and PAA is present throughout the year in both hemispheres and peaks during Northern Hemisphere summer. The addition of ocean VOC sources in the model increases annual surface cOHRmod by 10 % and improves model-measurement agreement for acetaldehyde particularly in winter but cannot resolve the model summertime bias. Doing so would require a 100 Tg yr−1 source of a long-lived unknown precursor throughout the year with significant additional emissions in the Northern Hemisphere summer. Improving the model bias for remote acetaldehyde and PAA is unlikely to fully resolve previously reported model global biases in OH and methane lifetime, suggesting that future work should examine the sources and sinks of OH over land.

Travis, Katherine R.↗
NASA NTRS2020

The Global Methane Budget 2000–2017

Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. Atmospheric emissions and concentrations of CH4 continue to increase, making CH4 the second most important human-influenced greenhouse gas in terms of climate forcing, after carbon dioxide (CO2). The relative importance of CH4 compared to CO2 depends on its shorter atmospheric lifetime, stronger warming potential, and variations in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in reducing uncertainties in the atmospheric growth rate arise from the variety of geographically overlapping CH4 sources and from the destruction of CH4 by short-lived hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to synthesize and stimulate new research aimed at improving and regularly updating the global methane budget. Following Saunois et al. (2016), we present here the second version of the living review paper dedicated to the decadal methane budget, integrating results of top-down studies (atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up estimates (including process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). For the 2008–2017 decade, global methane emissions are estimated by atmospheric inversions (a top-down approach) to be 576 Tg CH4/yr (range 550–594, corresponding to the minimum and maximum estimates of the model ensemble). Of this total, 359 Tg CH4/yr or ∼ 60 % is attributed to anthropogenic sources, that is emissions caused by direct human activity (i.e. anthropogenic emissions; range 336–376 Tg CH4/yr or 50 %–65 %). The mean annual total emission for the new decade (2008–2017) is 29 Tg CH4/yr larger than our estimate for the previous decade (2000–2009), and 24 Tg CH4/yr larger than the one reported in the previous budget for 2003–2012 (Saunois et al., 2016). Since 2012, global CH4 emissions have been tracking the warmest scenarios assessed by the Intergovernmental Panel on Climate Change. Bottom-up methods suggest almost 30 % larger global emissions (737 Tg CH4/yr, range 594–881) than top-down inversion methods. Indeed, bottom-up estimates for natural sources such as natural wetlands, other inland water systems, and geological sources are higher than top-down estimates. The atmospheric constraints on the top-down budget suggest that at least some of these bottom-up emissions are overestimated. The latitudinal distribution of atmospheric observation-based emissions indicates a predominance of tropical emissions (∼ 65 % of the global budget, < 30° N) compared to mid-latitudes (∼ 30 %, 30–60° N) and high northern latitudes (∼ 4 %, 60–90° N). The most important source of uncertainty in the methane budget is attributable to natural emissions, especially those from wetlands and other inland waters. Some of our global source estimates are smaller than those in previously published budgets (Saunois et al., 2016; Kirschke et al., 2013). In particular wetland emissions are about 35 Tg CH4/yr lower due to improved partition wetlands and other inland waters. Emissions from geological sources and wild animals are also found to be smaller by 7 Tg CH4/yr by 8 Tg CH4/yr, respectively. However, the overall discrepancy between bottom-up and top-down estimates has been reduced by only 5 % compared to Saunois et al. (2016), due to a higher estimate of emissions from inland waters, highlighting the need for more detailed research on emissions factors. Priorities for improving the methane budget include (i) a global, high-resolution map of water-saturated soils and inundated areas emitting methane based on a robust classification of different types of emitting habitats; (ii) further development of process-based models for inland-water emissions; (iii) intensification of methane observations at local scales (e.g., FLUXNET-CH4 measurements) and urban-scale monitoring to constrain bottom-up land surface models, and at regional scales (surface networks and satellites) to constrain atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) development of a 3D variational inversion system using isotopic and/or co-emitted species such as ethane to improve source partitioning.

methane budget↗
NASA NTRS2020

Understanding and Improving Model Representation of Aerosol Optical Properties for a Chinese Haze Event Measured During KORUS-AQ

KORUS-AQ was an international cooperative air quality field study in South Korea that measured local and remote sources of air pollution affecting the Korean Peninsula during May–June 2016. Some of the largest aerosol mass concentrations were measured during a Chinese haze transport event (24 May). Air quality forecasts using the WRF-Chem model with aerosol optical depth (AOD) data assimilation captured AOD during this pollution episode but overpredicted surface particulate matter concentrations in South Korea, especially PM2.5, often by a factor of 2 or larger. Analysis revealed multiple sources of model deficiency related to the calculation of optical properties from aerosol mass that explain these discrepancies. Using in situ observations of aerosol size and composition as inputs to the optical properties calculations showed that using a low-resolution size bin representation (four bins) underestimates the efficiency with which aerosols scatter and absorb light (mass extinction efficiency). Besides using finer-resolution size bins (8–16 bins), it was also necessary to increase the refractive indices and hygroscopicity of select aerosol species within the range of values reported in the literature to achieve better consistency with measured values of the mass extinction efficiency (6.7 m2 g−1 observed average) and light-scattering enhancement factor (f(RH)) due to aerosol hygroscopic growth (2.2 observed average). Furthermore, an evaluation of the optical properties obtained using modeled aerosol properties revealed the inability of sectional and modal aerosol representations in WRF-Chem to properly reproduce the observed size distribution, with the models displaying a much wider accumulation mode. Other model deficiencies included an underestimate of organic aerosol density (1.0 g cm−3 in the model vs. observed average of 1.5 g cm−3) and an overprediction of the fractional contribution of submicron inorganic aerosols other than sulfate, ammonium, nitrate, chloride, and sodium corresponding to mostly dust (17 %–28 % modeled vs. 12 % estimated from observations). These results illustrate the complexity of achieving an accurate model representation of optical properties and provide potential solutions that are relevant to multiple disciplines and applications such as air quality forecasts, health impact assessments, climate projections, solar power forecasts, and aerosol data assimilation.

Pablo E. Saide↗
NASA NTRS2022

A Coupled Evaluation of Operational MODIS and Model Aerosol Products for Maritime Environments Using Sun Photometry: Evaluation of the Fine and Coarse Mode

Although satellite retrievals and data assimilation have progressed to where there is a good skill for monitoring maritime Aerosol Optical Depth (AOD), there remains uncertainty in achieving further degrees of freedom, such as distinguishing fine and coarse mode dominated species in maritime environments (e.g., coarse mode sea salt and dust versus fine mode terrestrial anthropogenic emissions, biomass burning, and maritime secondary production). For the years 2016 through 2019, we performed an analysis of 550 nm total AOD550, fine mode AOD (FAOD550; also known as FM AOD in the literature), coarse mode AOD (CAOD550), and fine mode fraction (η550) between Moderate Resolution Spectral Imaging Radiometer (MODIS) V6.1 MOD/MYD04 dark target aerosol retrievals and the International Cooperative for Aerosol Prediction (ICAP) core four multi-model consensus (C4C) of analyses/short term forecasts that assimilate total MODIS AOD550. Differences were adjudicated by the global shipboard Maritime Aerosol Network (MAN) and selected island AERONET sun photometer observations with the application of the spectral deconvolution algorithm (SDA). Through a series of conditional and regional analyses, we found divergence included regions of terrestrial influence and latitudinal dependencies in the remote oceans. Notably, MODIS and the C4C and its members, while having good correlations overall, have a persistent +0.04 to +0.02 biases relative to MAN and AERONET for typical AOD550 values (84th% < 0.28), with the C4C underestimating significant events thereafter. Second, high biases in AOD550 are largely associated with the attribution of the fine mode in satellites and models alike. Thus, both MODIS and C4C members are systematically overestimating AOD550 and FAOD550 but perform better in characterizing the CAOD550. Third, for MODIS, findings are consistent with previous reports of a high bias in the retrieved Ångström Exponent, and we diagnosed both the optical model and cloud masking as likely causal factors for the AOD550 and FAOD550 high bias, whereas for the C4C, it is likely from secondary overproduction and perhaps numerical diffusion. Fourth, while there is no wind-speed-dependent bias for surface winds <12 m s−1, the C4C and MODIS AOD550s also overestimate CAOD550 and FAOD550, respectively, for wind speeds above 12 m/s. Finally, sampling bias inherent in MAN, as well as other circumstantial evidence, suggests biases in MODIS are likely MODIS and the C4C products have their own strengths and challenges for a given climate application and discuss needed research. even larger than what was diagnosed here. We conclude with a discussion on how

Aerosols↗
NASA NTRS2024

ISS Technology Demonstrations for Future Spaceflight Medical Systems

Throughout the history of human spaceflight, crewmembers have experienced various in-flight medical conditions including illness and injury. Planned missions to the Moon and Mars will require capabilities to maintain the health of future space travelers. Mass, power, and volume available in the vehicles and habitats for these missions will be severely constrained; resupply of resources will be limited or non-existent, as will opportunities for evacuation to Earth. Furthermore, ground-based support will be hampered by communication latencies and blackouts. These vehicle and mission constraints will necessitate a medical system that has been efficiently planned, providing on-board procedural guidance in addition to a variety of medical devices and consumable resources. Medical capabilities required for the diagnosis and treatment of potential medical conditions during future spaceflight missions may include real-time health monitoring, medical imaging, and biomarker analyses ( e.g., blood or urine). Terrestrial medicine shares these needs, thus many of these medical capabilities could likely be satisfied by Commercial-Off-The-Shelf (COTS) devices and methodologies; however, in some cases the unique space environment and increased mission duration will drive the need to modify technologies and the way care is provided. NASA’s Human Research Program (HRP) Exploration Medical Capability (ExMC) Element and Mars Campaign Office’s Exploration Medical Integrated Product Team (XMIPT) are working together to decrease medical risk during exploration missions. Flight-tested medical diagnostic and treatment technologies are necessary to effectively manage medical conditions relevant to exploration missions while meeting vehicle constraints, integrating with medical decision-support tools, and enabling increasingly Earth-independent operations. Several projects have leveraged the ISS as a testbed for exploration, including 1) i n- situ blood analysis, 2) medical inventory, 3) intravenous fluid generation, and 4) autonomous medical procedure guidance. Management of several in-flight medical conditions, such as bacterial and viral infections and acute radiation syndrome, is dramatically improved with ability to assess blood cell populations, electrolytes, and metabolites. I n December 2020 and January 2021 ExMC performed an ISS technology demonstration (Tech Demo) of the HemoCue® WBC DIFF analyzer (HemoCue, Brea, CA), a COTS device that was modified to enable functionality in a spaceflight environment. This Tech Demo marked the first time that hematology measurements were successfully performed real-time in microgravity. Also modified and demonstrated was the reusable Handheld Electrolytes and Lab Technology for Humans (rHEALTH) ONE analyzer (rHEALTH, Bedford, MA), which uses flow cytometry and sheath-based hydrodynamic focusing methodologies. The rHEALTH ONE ISS Tech Demo in May 2022 demonstrated test results obtained in-flight matched those on the ground. NASA currently relies on crew self-reporting to manage and maintain medical inventory on ISS.The ability to maintain an accurate inventory becomes more critical during long duration missions since the crew will need to be increasingly autonomous in finding and utilizing medical items, including those scenarios when alternative treatments need to be considered due to limited or no resupply. HRP’s Medical Consumables Tracking (MCT) project was developed by ZIN Technologies, Inc. (Cleveland, OH), and demonstrated real-time tracking of medical supplies aboard the ISS between December 2016 and July 2018. The MCT system design utilized Radio Frequency Identification Device (RFID) technology to perform automated inventory and was installed in the Crew Health Care System (CHeCS) Resupply Stowage Rack (RSR). The challenge of limited shelf life, exacerbated by the lack of resupply opportunities, affects a plethora of medical system components including consumables, pharmaceuticals, and intravenous (IV) fluid. In 2010, ExMC funded ZIN Technologies, Inc. (Cleveland, OH), to develop the Intravenous Fluid Generation (IVGEN) system. IV fluids were successfully generated with IVGEN using the potable water supply on ISS during ISS Expedition 23. The XMIPT is in the process of developing a miniaturized version of the original IVGEN hardware for a future Tech Demo aboard the ISS. Current ISS medical operations rely heavily on preflight training and real-time remote guidance, both of which become impractical or impossible for exploration missions. The primary goal of the Autonomous Medical Officer Support (AMOS) Software Tech Demos on ISS was to confirm telemedical proof-of-concept for autonomous medical imaging in an operational setting. This novel software tool shifts emphasis from preflight training and real-time remote guidance to in-flight just-in-time instruction, a new and necessary paradigm for crew medical autonomy. AMOS introduces a novel, streamlined skill management concept for exploration missions featuring comprehensive training and guidance modules for ultrasound examinations using the ISS Ultrasound 2 (a modified GE Vivid-q™; General Electric HealthCare, Chicago, IL). With no prior crew training or remote guidance, two Tech Demos on the ISS (April 2020 and June 2022) resulted in high quality, clinically useful image sets. We will provide a review of historical, current, and planned medical devices and technologies considered for inclusion within future spaceflight medical systems and summarize hardware development activities and medical device tech demos conducted on the ISS.

Astronaut health and performance↗
NASA NTRS2016

First Flight of the Gamma-Ray Imager Polarimeter for Solar Flares (GRIPS) Instrument

The Gamma-Ray Imager/Polarimeter for Solar ares (GRIPS) instrument is a balloon-borne telescope designed to study solar-flare particle acceleration and transport. We describe GRIPS's first Antarctic long-duration flight in January 2016 and report preliminary calibration and science results. Electron and ion dynamics, particle abundances and the ambient plasma conditions in solar flares can be understood by examining hard X-ray (HXR) and gamma-ray emission (20 keV to 10 MeV). Enhanced imaging, spectroscopy and polarimetry of flare emissions in this energy range are needed to study particle acceleration and transport questions. The GRIPS instrument is specifically designed to answer questions including: What causes the spatial separation between energetic electrons producing hard X-rays and energetic ions producing gamma-ray lines? How anisotropic are the relativistic electrons, and why can they dominate in the corona? How do the compositions of accelerated and ambient material vary with space and time, and why? GRIPS's key technological improvements over the current solar state of the art at HXR/gamma-ray energies, the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), include 3D position-sensitive germanium detectors (3D-GeDs) and a single-grid modulation collimator, the multi-pitch rotating modulator (MPRM). The 3D-GeDs have spectral FWHM resolution of a few hundred keV and spatial resolution less than 1cu mm. For photons that Compton scatter, usually greater or equal to 150 keV, the energy deposition sites can be tracked, providing polarization measurements as well as enhanced background reduction through Compton imaging. Each of GRIPS's detectors has 298 electrode strips read out with ASIC/FPGA electronics. In GRIPS's energy range, indirect imaging methods provide higher resolution than focusing optics or Compton imaging techniques. The MPRM grid-imaging system has a single-grid design which provides twice the throughput of a bi-grid imaging system like RHESSI. The grid is composed of 2.5 cm deep tungsten-copper slats, and quasi-continuous FWHM angular coverage from 12.5-162 arcsecs are achieved by varying the slit pitch between 1-13 mm. This angular resolution is capable of imaging the separate magnetic loop footpoint emissions in a variety of are sizes. In comparison, RHESSI's 35-arcsec resolution at similar energies makes the footpoints resolvable in only the largest ares.

Duncan, Nicole↗
NASA NTRS2019

Global Carbon Budget 2019

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (E(FF)) are based on energy statistics and cement production data, while emissions from land use change (E(LUC)), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (G(ATM)) is computed from the annual changes in concentration. The ocean CO2 sink (S(OCEAN)) and terrestrial CO2 sink (S(LAND)) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (B(IM)), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), E(FF) was 9.5±0.5 GtC/yr, E(LUC) 1.5±0.7 GtC/yr, G(ATM) 4.9±0.02 GtC/yr (2.3±0.01 ppm/yr), S(OCEAN) 2.5±0.6 GtC/yr, and S(LAND) 3.2±0.6 GtC/yr, with a budget imbalance B(IM) of 0.4 GtC/yr indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in E(FF) was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC/yr, reaching 10 GtC/yr for the first time in history, E(LUC) was 1.5±0.7 GtC/yr, for total anthropogenic CO2 emissions of 11.5±0.9 GtC/yr (42.5±3.3 GtCO2). Also for 2018, G(ATM) was 5.1±0.2 GtC/yr (2.4±0.1 ppm/yr), S(OCEAN) was 2.6±0.6 GtC/yr, and S(LAND) was 3.5±0.7 GtC/yr, with a B(IM) of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in E(FF) of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC/yr persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013).

global carbon budget 2019↗
NASA NTRSDate not supplied

A Combined Al-Mg/Pb-Pb Age of the Solar System

Astrophysical models of planet formation and protoplanetary disk evolution demand precise and accurate timing of the sequence of events in the solar nebula, relative to a time t=0, usually taken to be during the short epoch of CAI (Ca-rich, Al-rich inclusion) formation. Most CAIs formed withlive26Al (mean-life τ26= 1.034 Myr [1]), with an abundance 26Al/27Al ≈ (26Al/27Al)SS= 5.23 × 10-5[2]. We adopt this as the widespread level of 26Al in the solar nebula at t=0. Assuming spatial homogeneity of 26Al, an inclusion that had less 26Al, (26Al/27Al)0, formed a time Δt26= τ26ln[(26Al/27Al)SS/ (26Al/27Al)0] after t=0.These ages are typical precise to within ±0.1 Myr. Igneous bulk meteorites and inclusions can be relatively dated by the Al-Mg chronometer, but only ifΔt26<6 Myr. The Pb-Pb system is useful as a longer relative chronometer. It yields absolute ages tPb using 207Pb/206Pb, 206Pb/204Pb, and 238U/235U ratios measured indifferent portions of a sample, assuming certain half-lives [4]. These absolute ages are uncertain to within ±9 Myr due to uncertainties in the 235U half-life[3], but times of formation ΔtPb= tCAI–tPb relative to t=0, are more precise(±0.5Myr),iftCAI can be found. Here, tCAI means the Pb-Pb age that would be measured in CAIs using the half-lives the community typically uses, if they achieved isotopic closure at t=0. Unfortunately, direct Pb-Pb dating of CAIs has not definitively determined tCAI. Based on four CAIs with canonical (26Al/27Al)0,[5,6] found tPb= 4567.30 ± 0.16 Myr. No other CAI ages with measured 238U/235U have been reported in the refereed literature, but there are hints of other CAIs with ages tPb= 4568.0 ± 0.2 Myr [7] and tPb= 4568.3 ± 0.2 Myr [8].It is unclear whether anyof these igneous type B CAIs isotopically closed at t=0 or represents tCAI. Instead of measurements, we advocate finding tCAI by minimizing the discrepancies between the Al-Mg and Pb-Pb chronometers. Assuming Δt26=ΔtPb, we find the implied t’CAI= tPb+Δt26, then define t*CAIas the weighted mean of the t’CAI. t*CAIis the best guess for the Pb-Pb age of t=0; the assumption of homogeneity is justified if the t’CAI cluster within errors around t*CAI. This statistical approach is similar to, but improves on, that of[9]. We find t*CAI= 4568.73 ± 0.16 Myr. Below we discuss our methodology and the implications of this age for CAIs, 1.4 Myr older than the reported and typically used age 4567.30±0.16 Myr. Methods: We base our estimate of t*CAIon five achondrites for which published (26Al/27Al)0and Pb-Pb ages exist: the quenched angrites D’Orbigny, Sahara 99555 (SAH 99555), and Northwest Africa (NWA) 1670; the pseudo-eucrite Asuka 881394; and the inner disk achondrite. All are “NC” (non-carbonaceous) achondrites that likely cooled quickly enough that the Al-Mg and Pb-Pb systems achieved isotopic closure simultaneously. We also considered the “CC” (carbonaceous chondrite-like) achondrites NWA 2796 and NWA 6704, butdo not include them in our fit. Al-Mg and Pb-Pb seem not to have closed simultaneously, possibly because formation in the outer disk from volatile-rich composition led to slower cooling. Of the 8 chondrules from NWA 5697 measured by [20,21], we also consider the 4 for which 238U/235U was measured: 2-C1, 5-C2, 3-C5, 11-C1.Depending on their post-formation thermal histories, the Al-Mg and Pb-Pb systems in chondrules may or may not have closed simultaneously. Table 1: (26Al/27Al)0, Pb-Pb ages of selected samplesSample(26Al/27Al)0/ 10-6RefPb-PbRefD’Orbigny3.98±0.15104563.43±0.19♮10-12SAH 995553.64±0.18104563.88±0.2712NWA16705.92±0.59104564.39±0.24*10Asuka 88139413.1±0.5613-154564.98±0.1715NWA 73253.03±0.14164563.4±2.616NWA 27963.94±0.16174562.89±0.5917NWA 67043.03±0.14184562.76±0.26192-C17.56±1.53204567.57±0.56*215-C27.04±1.51204567.54±0.52*213-C58.85±1.83204566.20±0.63*2111-C15.55±1.84204565.84±0.72*21*regression based on one subset of data points ♮weighted mean of two datasets Pb-Pb ages are proportional to the intercept of the line formed by linear regression of 207Pb/206Pb vs. 206Pb/204Pb data from various washes, leachates and residues of acid dissolution of a sample. Because contamination by terrestrial or primordial Pb is pervasive, some fractions must be excluded from regressions to ensure a fit with acceptable mean squares weighted deviation (MSWD). Usually points are excluded based on low [Pb], or low 206Pb/204Pb ratio(low radiogenic component), with single outliers identified [11,12,15,16,17]. In the starred examples (Table 1)and the case of 3 CAI Pb-Pb ages [5], up to half the points were excluded solely because did not fit a pre-determined line. This approach is vulnerable to confirmation bias and produces fits with low MSWD and too-low Pb-Pb age uncertainty. Regressing the same data points as [10], were produce the Pb-Pb age of NWA 1670 of 4564.39±0.24 Myr. But selecting other combinations of data points, other, equally valid, isochrons yield ages from 4563.77±0.21 Myr to 4564.64±0.23 Myr. Similar arguments apply to the Pb-Pb isochrons built by [21] for chondrules 2-C1 (we find 4567.33±0.44 to 4567.85±0.46 Myr), 5-C2 (4566.84±0.53 to 4567.70±0.44 Myr), 3-C5 (4565.84±0.54to 4567.04±0.54)and 11-C1 (4565.36±0.51 to 4565.74±0.45 Myr). Our adopted ages for these and NWA 1670 are listed in Table 2.Table 2. tCAI estimated from various components, using our regressions for the chondrules & NWA 1670.SampleΔt26(Myr)tPb(Myr)t’CAI(Myr)D’Orbigny5.05±0.044563.43±0.194568.48±0.19SAH 995555.14±0.054563.88±0.274569.02±0.27NWA16704.64±0.104564.21±0.634568.85±0.67Asuka 8813943.81±0.044564.98±0.174568.79±0.17NWA 73255.33±0.054563.4±2.64568.7±2.6NWA 27965.06±0.044562.89±0.594567.95±0.59NWA 67045.29±0.134562.76±0.264568.05±0.292-C12.00±0.214567.59±0.704569.59±0.725-C22.07±0.224567.23±0.914569.30±0.933-C51.84±0.214566.44±1.124568.28±1.1411-C12.32±0.344565.52±0.664567.84±0.73achondrite4568.72±0.16chondrules4568.76±0.58combined4568.73±0.16A weighted average of the five NC achondrites(or just D’Orbigny, SAH 99555 and Asuka 881394)yields t*CAI= 4568.72 ± 0.16Myr. All are consistent with this value to within 1.8σ, and MSWD=1.5. Including the 4 U-corrected chondrules, t*CAI= 4568.73± 0.16Myrwith MSWD=1.66, which is statistically significant. All chondrules and NC achondrites are consistent with this to within 1.8σ, (Figure 1).Figure 1. Al-Mg formation times after t=0 vs. Pb-Pb ages. The five NC achondrites and four chondrules are consistent with a Pb-Pb age of t=0 of 4568.7 Myr. Discussion: The data from achondrites and chondrules are consistent with a single Pb-Pb age at t=0, justifying the assumption of 26Al homogeneity. The age, 4568.7 Myr, is ≈1.4 Myr older than the commonly accepted Pb-Pb age of CAIs that formed with canonical 26Al/27Al at t=0 [3]. Others have interpreted the discrepancy to signify 26Al heterogeneity in the CAI-forming region[5,21]. We suggest instead that CAIs were exposed to transient heating events that reset the Pb-Pb system without disturbing the Al-Mg system. Notably, chondrules typically experienced transient heating at these times in the nebula [22]. If so, direct measurements of CAIs will not yield as reliable a Pb-Pb age of t=0 as statistical approaches like this and that of [9].References:[1] Auer et al. 2009. [2] Jacobsen, B et al. 2008, EPSL 272, 353-364. [3] Tissot, Fet al. 2017, GCA 213, 593-617.[4] Villa, I et al. 2016, GCA 172, 387-392.[5] Amelin, Y et al. 2010, EPSL 300, 343-350.[6] Connelly, Jet al. 2012, Science 338, 651.[7] Bouvier, Aet al. 2011, LPICo 1639, 9054. [8] Bouvier, Aand Wadhwa, M2010, Nat Geosci 3, 637-641.[9] Nyquist, Let al. 2009, GCA 73, 5115-5136.[10] Schiller et al. 2015. [11] Wadhwa & Brennecka 2012. [12] Tissot et al. 2017. [13] Nyquist et al. 2003. [14] Wadhwa et al. 2009. [15] Wimpenny et al. 2019, GCA 244, 478-501.[16] Koefoed et al. 2016, GCA 183, 31-45. [17] Bouvier, A et al. 2011, GCA 75, 5310-5323. [18] Sanborn, Met al. 2019, GCA 245, 577-596. [19] Amelin, Yet al. 2019, GCA 245, 628-642.[20] Bollard, Jet al. 2017, Sci Adv 3 ,e1700407. [21] Bollard, Jet al. 2019, GCA 260, 62-83.[22] Villeneuve, J et al. 2009, Science 325, 985

S. J. Desch↗
NASA NTRS2018

On the 2018 Outburst of the Accreting Millisecond X-Ray Pulsar Swift J1756.9-2508 As Seen with NICER

We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar SwiftJ1756.92508during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combiningour NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we alsostudied the long-term spin and orbital evolution of this source. We find that the binary system is well describedby a constant orbital period model, with an upper limit on the orbital period derivative of Pb < 7.4 ´ 10-13 ss1.Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times,finding R.A.=17h56m57 18±0 08 and decl.=25°0627 8±3 5, while simultaneously measuring thelong-term spin frequency derivative as n = -7.3 ´ 10-16 Hzs1. We briefly discuss the implications of thesemeasurements in the context of the wider population of accreting millisecond pulsars. We reported on the coherent timing analysis of the 2018 outburst of Swift J1756 as observed with NICER. Consistent with analyses of the previous outbursts (Krimm et al. 2007b; Patruno et al. 2010), we find that the X-ray pulsations have energy dependent amplitudes; the fractional amplitude of the fundamental increases with energy, whereas the fractional amplitude of the harmonic shows a slight decline with energy. This energy dependent behavior is not unusual in AMXPs (Patruno & Watts 2012) and can be interpreted in terms of the thermal emission from the stellar hotspot and reprocessing in the accretion column (e.g., Gierliński et al. 2002; Ibragimov & Poutanen 2009). The pulse arrival times of the 2018 outburst are well described by a timing model consisting of a circular orbit with a constant spin frequency. The pulse phases with respect to this model do not show spurious residuals with time or orbital phase, and no evidence is found that the pulse arrival times exhibit an additional delay associated with passing through the gravitational well of the companion star (Shapiro delay). We note, however, that the expected Shapiro delay is given as (Shapiro et al. 1971) Equation (5) where Φ is the orbital phase, G is the gravitational constant, c is the speed of light, and i is the inclination. Even for the maximum allowed companion mass, ${M}_{C}=0.030\,{M}_{\odot }$ (Krimm et al. 2007b, but see Section 4.2 for more details) and an inclination of 90°, the largest delay we can expect is only 4 μs. As this time-delay is smaller than the uncertainty on our phase residuals by nearly two orders of magnitude (see Figure 1), we are not sensitive to Shapiro delays in Swift J1756. Comparing our measurements for the 2018 outburst with those of the 2007 and 2009 outbursts as observed with RXTE, we analyzed the long-term evolution of this source. We found that the binary system is consistent with having a constant orbital period and that the pulsar shows a spin frequency derivative of $\dot{\nu }=-7.3\times {10}^{-16}\,\mathrm{Hz}\,{{\rm{s}}}^{-1}$. 4.1. Spin-down Evolution The long-term spin frequency derivative measured in Swift J1756 is of the same order as the spin frequency derivatives measured in other AMXPs (Hartman et al. 2008; Patruno 2010; Riggio et al. 2011). This frequency change is most likely driven by the neutron star's loss of rotational energy. If so, then the spin-down luminosity is given as Equation (6) where I represents the neutron star moment of inertia. The long-term spin-down of a neutron star is usually assumed to be dominated by the braking torque associated with a spinning magnetic field. Assuming this mechanism is responsible for the observed spin-down in Swift J1756, we can compute the magnetic dipole moment as (Spitkovsky 2006) Equation (7) where α is the misalignment angle between the rotational and magnetic poles. Considering α = 0°–90°, we then find a magnetic field strength of $B\simeq (4\mbox{--}6)\times {10}^{8}$ G at the stellar magnetic poles. This magnetic field strength estimate is in line with those obtained for other accreting millisecond pulsars (see Mukherjee et al. 2015 and references therein). 4.2. Orbit Evolution The observed long-term binary evolution of Swift J1756 is consistent with this source having a constant orbital period and a lower limit on the evolutionary timescale of Equation (8) Binary evolution theory predicts that systems of this type evolve due to angular momentum loss through gravitational radiation (Kraft et al. 1962; Rappaport et al. 1982; Verbunt 1993). For conservative mass transfer, the binary period derivative is given by di Salvo et al. (2008), Equation (9) where MNS is the neutron star mass, $q={M}_{C}/{M}_{\mathrm{NS}}$ is the binary mass ratio, and −1/3 < n < 1 is the mass–radius index of the companion star. Depending on the source inclination, Krimm et al. (2007b) derived a companion mass of ${M}_{C}\,=0.007\mbox{--}0.022\,{M}_{\odot }$ for a neutron star mass of 1.4 ${M}_{\odot }$. For a neutron star mass of 2.2 ${M}_{\odot }$, the allowed range increased to ${M}_{C}=0.009\mbox{--}0.030\,{M}_{\odot }$. In both cases, they assumed an upper limit on the inclination of i < 85°, motivated by the fact that Swift J1756 does not show eclipses in its light curve. Accounting for the extreme cases of stellar masses and n, the binary may either be contracting or expanding. In either case, however, the rate of change is limited to $| {\dot{P}}_{b}| \lesssim 7\times {10}^{-14}$ s s−1, which is well below the upper limit obtained in this work. Although the binary evolution timescale we obtain for Swift J1756 is consistent with theory, it is worth noting that this is not generally true for low-mass X-ray binaries (see Patruno et al. 2017, for a comprehensive discussion). The AMXP SAX J1808.4–3658, in particular, has been found to evolve on a much shorter timescale, with a first derivative on the orbital period of $3.5\times {10}^{-12}$ s s−1 (Hartman et al. 2008; Patruno et al. 2012; Sanna et al. 2017a). Two models have been proposed to explain this discrepancy: highly nonconservative mass transfer due to irradiation of the companion star by the pulsar (di Salvo et al. 2008; Burderi et al. 2009), and spin–orbit coupling in the companion star (Hartman et al. 2008, 2009). While the latter depends on the companion star, and may vary from source to source, the former should operate in all AMXPs (see also Patruno 2017; Sanna et al. 2017c), including Swift J1756. The spin-down luminosity impinging on the companion star can be estimated as Equation (10) where ${\dot{E}}_{\mathrm{abl}}$ is the ablation luminosity, RL2 is the Roche lobe radius of the companion (Eggleton 1983), and a the binary separation. The irradiation fraction is $f={\dot{E}}_{\mathrm{abl}}/{\dot{E}}_{\mathrm{sd}}$, which, accounting for the range of allowed neutron star and companion masses, evaluates to f = 0.15%–0.35%. The associated mass loss for the companion is given by Equation (11) such that, assuming an efficiency of η = 100%, ${\dot{M}}_{C}\,\sim -3\times {10}^{-10}\,{M}_{\odot }$ yr−1. The effect of this mass loss on the orbital period follows through the relation (Frank et al. 2002) Equation (12) giving a period derivative due to mass loss of ${\dot{P}}_{b,\mathrm{ML}}\,=5\times {10}^{-12}$ s s−1. This value is well above our limit on the period derivative. Hence, in order for this mechanism to be consistent with our observations of Swift J1756, the efficiency at which the companion star converts the incident luminosity into mass loss must be η < 15%. This value is very different from the 40% required in SAX J1808.4–3658 (Patruno et al. 2016) and is instead in line with the <5% efficiency determined for IGR J00291+5934 (Patruno 2017). This work was supported by NASA through the NICER mission and the Astrophysics Explorers Program, and made use of data and software provided by the High Energy Astrophysics Science Archive Research Center (HEASARC). P.B. was supported by an NPP fellowship at NASA Goddard Space Flight Center. D.A. acknowledges support from the Royal Society.

Bult, Peter↗