Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic error”

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 361 records · Page 20

Quantitative capabilities of STXM to measure spatially resolved organic volume fractions of mixed organic/inorganic particles

Scanning transmission X-ray microscopy coupled with near-edge X-ray absorption and fine structure (STXM-NEXAFS) spectroscopy can be used to characterize the morphology and composition of aerosol particles. Here, two inorganic/organic systems are used to validate the calculation of organic volume fraction (OVF) and determine the level of associated error by using carbon K-edge STXM data at 278, 285.4, 288.6, and 320 eV. Using the mixture of sodium chloride and sucrose as one system and ammonium sulfate and sucrose as another, three solutions were made with 10:1, 1:1, and 1:10 mass ratios (inorganic to organic). The OVFs of the organic-rich aerosols of both systems deviated from the bulk OVF by less than 1%, while the inorganic-rich aerosols deviated by approximately 1 %. Aerosols from the equal mass mixture deviated more (about 4 %) due to thick inorganic regions exceeding the linear range of Beer's law. These calculations were performed after checking the data for poor image alignment, defocusing issues, and particles too thick to be analyzed. The potential for systematic error in the OVF calculation was also tested by assuming the incorrect composition. There is a small (about 0.5 %) OVF difference if the organic is erroneously assumed to be adipic acid rather than the known organic, sucrose. A much larger difference (up to 25 %) is seen if sodium chloride is assumed instead of ammonium sulfate. These results show that the OVF calculations are fairly insensitive to the organic while being much more sensitive to the choice of inorganic.

47 OTHER INSTRUMENTATION↗

Determining the bubble nucleation efficiency of low-energy nuclear recoils in superheated C$_3$F$_8$ dark matter detectors

The bubble nucleation efficiency of low-energy nuclear recoils in superheated liquids plays a crucial role in interpreting results from direct searches for weakly interacting massive particle (WIMP) dark matter. The PICO Collaboration presents the results of the efficiencies for bubble nucleation from carbon and fluorine recoils in superheated C$_3$F$_8$ from calibration data taken with 5 distinct neutron spectra at various thermodynamic thresholds ranging from 2.1 keV to 3.9 keV. Instead of assuming any particular functional forms for the nuclear recoil efficiency, a generalized piecewise linear model is proposed with systematic errors included as nuisance parameters to minimize model-introduced uncertainties. A Markov-Chain Monte-Carlo (MCMC) routine is applied to sample the nuclear recoil efficiency for fluorine and carbon at 2.45 keV and 3.29 keV thermodynamic thresholds simultaneously. The nucleation efficiency for fluorine was found to be $\geq 50\, \%$ for nuclear recoils of 3.3 keV (3.7 keV) at a thermodynamic Seitz threshold of 2.45 keV (3.29 keV), and for carbon the efficiency was found to be $\geq 50\, \%$ for recoils of 10.6 keV (11.1 keV) at a threshold of 2.45 keV (3.29 keV). Simulated data sets are used to calculate a p-value for the fit, confirming that the model used is compatible with the data. The fit paradigm is also assessed for potential systematic biases, which although small, are corrected for. Additional steps are performed to calculate the expected interaction rates of WIMPs in the PICO-60 detector, a requirement for calculating WIMP exclusion limits.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

THE DISK WIND IN THE RAPIDLY SPINNING STELLAR-MASS BLACK HOLE 4U 1630–472 OBSERVED WITH NuSTAR

We present an analysis of a short NuSTAR observation of the stellar-mass black hole and low-mass X-ray binary 4U 1630–472. Reflection from the inner accretion disk is clearly detected for the first time in this source, owing to the sensitivity of NuSTAR. With fits to the reflection spectrum, we find evidence for a rapidly spinning black hole, a{sub ∗}=0.985{sub −0.014}{sup +0.005} (1σ statistical errors). However, archival data show that the source has relatively low radio luminosity. Recently claimed relationships between jet power and black hole spin would predict either a lower spin or a higher peak radio luminosity. We also report the clear detection of an absorption feature at 7.03 ± 0.03 keV, likely signaling a disk wind. If this line arises in dense, moderately ionized gas (log ξ=3.6{sub −0.3}{sup +0.2}) and is dominated by He-like Fe XXV, the wind has a velocity of v/c=0.043{sub −0.007}{sup +0.002} (12900{sub −2100}{sup +600} km s{sup –1}). If the line is instead associated with a more highly ionized gas (log ξ=6.1{sub −0.6}{sup +0.7}), and is dominated by Fe XXVI, evidence of a blueshift is only marginal, after taking systematic errors into account. Our analysis suggests the ionized wind may be launched within 200-1100 Rg, and may be magnetically driven.

79 ASTRONOMY AND ASTROPHYSICS↗

Systematic problems with using dark matter simulations to model stellar halos

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be 'painted' onto dark matter (DM) particles in progenitor satellites; (2) pure DM simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted DM particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure DM simulation that does not form a baryonic disk. We find that the 'painting' simplification reduces the halo concentration and internal structure, predominantly because painted DM particles have different kinematics from star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure DM simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on DM-only simulations when such differences are less than an order of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

Accurate determination of the electron spin polarization in magnetized iron and nickel foils for Møller polarimetry

The Møller polarimeter in Hall A at Jefferson Lab in Newport News, VA, has provided reliable measurements of electron beam polarization for the past two decades. Past experiments have typically required polarimetry at the 1% level of absolute uncertainty which the Møller polarimeter has delivered. However, the upcoming proposed experimental program including MOLLER and SoLID have stringent requirements on beam polarimetry precision at the level of 0.4% (The MOLLER Collaboration, 2014; The SoLID collaboration, 2019), requiring a systematic re-examination of all the contributing uncertainties. Møller polarimetry uses the double polarized scattering asymmetry of a polarized electron beam on a target with polarized atomic electrons. The target is a ferromagnetic material magnetized to align the spins in a given direction. In Hall A, the target is a pure iron foil aligned perpendicular to the beam and magnetized out of plane parallel or antiparallel to the beam direction. The acceptance of the detector is engineered to collect scattered electrons close to 90° in the center of mass frame where the analyzing power is a maximum (-7/9). One of the leading systematic errors comes from determination of the target foil polarization. Polarization of a magnetically saturated target foil requires knowledge of both the saturation magnetization and g', the electron g-factor which includes components from both spin and orbital angular momentum from which the spin fraction of magnetization is determined. Target foil polarization has been previously addressed in a 1997 publication “A precise target for Møller polarimetry” by de Bever et al. (1997) at a level of precision sufficient for experiments up to this point. Several shortcomings with the previous published value require revisiting the result prior to MOLLER. Here, we utilize the existing world data to provide a best estimate for target polarization for both nickel and iron foils including uncertainties in magnetization, high-field and temperature dependence, and fractional contribution to magnetization from orbital effects. We determine the foil electron spin polarization at 294 K to be 0.08020 ± 0.00018 (@4 T applied field) for iron and 0.018845+/-0.000053 (@2 T applied field) for nickel. Finally, we conclude with a brief discussion of additional systematic uncertainties to Møller polarimetry using this technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

A Hard Look at the Neutron Stars and Accretion Disks in 4U 1636-53, GX 17+2, and 4U 1705-44 with NuStar

We present NuSTAR observations of neutron star (NS) low-mass X-ray binaries: 4U 1636-53, GX 17+2, and 4U 1705-44. We observed 4U 1636-53 in the hard state, with an Eddington fraction, F{sub Edd}, of 0.01; GX 17+2 and 4U 1705-44 were in the soft state with fractions of 0.57 and 0.10, respectively. Each spectrum shows evidence for a relativistically broadened Fe K{sub α} line. Through accretion disk reflection modeling, we constrain the radius of the inner disk in 4U 1636-53 to be R{sub in}=1.03±0.03 ISCO (innermost stable circular orbit), assuming a dimensionless spin parameter a{sub ∗}=cJ/GM{sup 2}=0.0, and R{sub in}=1.08±0.06 ISCO for a{sub ∗}=0.3 (errors quoted at 1σ). This value proves to be model independent. For a{sub ∗}=0.3 and M=1.4 M{sub ⊙}, for example, 1.08 ± 0.06 ISCO translates to a physical radius of R=10.8±0.6 km, and the NS would have to be smaller than this radius (other outcomes are possible for allowed spin parameters and masses). For GX 17+2, R{sub in}=1.00--1.04 ISCO for a{sub ∗}=0.0 and R{sub in}=1.03--1.30 ISCO for a{sub ∗}=0.3. For a{sub ∗}=0.3 and M=1.4 M{sub ⊙}, R{sub in}=1.03--1.30 ISCO translates to R=10.3--13.0 km. The inner accretion disk in 4U 1705-44 may be truncated just above the stellar surface, perhaps by a boundary layer or magnetosphere; reflection models give a radius of 1.46–1.64 ISCO for a{sub ∗}=0.0 and 1.69–1.93 ISCO for a{sub ∗}=0.3. We discuss the implications our results may have on the equation of state of ultradense, cold matter and our understanding of the innermost accretion flow onto NSs with low surface magnetic fields, and systematic errors related to the reflection models and spacetime metric around less idealized NSs.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-E-T-642: Safety Procedures for Error Elimination in Cyberknife Stereotactic Radiosurgery (SRS) and Stereotactic Body Radiation Therapy (SBRT)

Purpose: Cyberknife system is used for providing stereotactic radiosurgery (SRS) and stereotactic body radiation therapy (SBRT) hypofractionation scheme. The whole treatment delivery is based on live imaging of the patient. The minor error made at any stage may bring severe radiation injury to the patient or damage to the system itself. Several safety measures were taken to make the system safer. Methods: The radiation treatment provided thru a 6MV linac attached to Kuka robot (Cyberknife G4, Accuray Inc. Sunnyvale, CA, USA). Several possible errors were identified related to patient alignment, treatment planning, dose delivery and physics quality assurance. During dose delivery, manual and visual checks were introduced to confirm pre and intra-treatment imaging to reduce possible errors. One additional step was introduced to confirm that software tracking-tools had worked correctly with highest possible confidence level. Robotic head move in different orientations over and around the patient body, the rigidity of linac-head cover and other accessories was checked periodically. The vender was alerted when a tiny or bigger piece of equipment needed additional interlocked support. Results: As of our experience treating 525 patients on Cyberknife during the last four years, we saw on and off technical issues. During image acquisition, it was made essential to follow the site-specific imaging protocols. Adequate anatomy was contoured to document the respective doses. Followed by auto-segmentation, manual tweaking was performed on every structure. The calculation box was enclosing the whole image during the final calculation. Every plan was evaluated on slice-by slice basis. To review the whole process, a check list was maintained during the physics 2nd-check. Conclusion: The implementation of manual and visual additional checks introduced along with automated checks for confirmation was found promising in terms of reduction in systematic errors and making the system safer than before.

60 APPLIED LIFE SCIENCES↗

Resolve survey photometry and volume-limited calibration of the photometric gas fractions technique

We present custom-processed ultraviolet, optical, and near-infrared photometry for the REsolved Spectroscopy of a Local VolumE (RESOLVE) survey, a volume-limited census of stellar, gas, and dynamical mass within two subvolumes of the nearby universe (RESOLVE-A and RESOLVE-B). RESOLVE is complete down to baryonic mass ∼10{sup 9.1−9.3} M{sub ⊙}, probing the upper end of the dwarf galaxy regime. In contrast to standard pipeline photometry (e.g., SDSS), our photometry uses optimal background subtraction, avoids suppressing color gradients, and employs multiple flux extrapolation routines to estimate systematic errors. With these improvements, we measure brighter magnitudes, larger radii, bluer colors, and a real increase in scatter around the red sequence. Combining stellar mass estimates based on our optimized photometry with the nearly complete H i mass census for RESOLVE-A, we create new z = 0 volume-limited calibrations of the photometric gas fractions (PGF) technique, which predicts gas-to-stellar mass ratios (G/S) from galaxy colors and optional additional parameters. We analyze G/S-color residuals versus potential third parameters, finding that axial ratio is the best independent and physically meaningful third parameter. We define a “modified color” from planar fits to G/S as a function of both color and axial ratio. In the complete galaxy population, upper limits on G/S bias linear and planar fits. We therefore model the entire PGF probability density field, enabling iterative statistical modeling of upper limits and prediction of full G/S probability distributions for individual galaxies. These distributions have two-component structure in the red color regime. Finally, we use the RESOLVE-B 21 cm census to test several PGF calibrations, finding that most systematically under- or overestimate gas masses, but the full probability density method performs well.

79 ASTRONOMY AND ASTROPHYSICS↗

Constraints on Neutrino Oscillation Parameters from Neutrinos and Antineutrinos with Machine Learning

NOvA is a two detector, long baseline neutrino oscillation experiment measuring the oscillations of muon neutrinos from the \numi neutrino beam over a baseline of \SI{810}{km}. The experiment uses four oscillation channels, $\numu \rightarrow \numu$, $\numubar \rightarrow \numubar$, $\numu \rightarrow \nue$, and $\numubar \rightarrow \nuebar$, with a peak neutrino energy of \SI{1.8}{GeV}. This dissertation describes the analysis of these channels using a dataset of $13.6\times10^{20}$ protons on target neutrino beam mode and $12.5\times10^{20}$ protons on target antineutrino beam mode. The analysis makes use of improvements in the treatment of systematic uncertainties and machine learning techniques to reconstruct neutrino interactions. A technique for decorrelating systematic errors using principle component analysis was utilized to reduce and optimize neutrino cross section and beam related uncertainties. The improved machine learning algorithms make use of convolutional ne ural net works for neutrino event classification, particle classification, and instance segmentation. The selection of neutrino signal events utilizing the neutrino event classifier shows an efficiency of 63\% for the selection of electron neutrinos in neutrino beam mode and 75\% for electron antineutrinos in antineutrino beam mode. Using this algorithm, 82 appearing electron neutrino candidates and 33 appearing electron antineutrino candidates were observed with expected backgrounds of 26.8 and 14.0 respectively. In addition, 211 surviving muon neutrino candidates and 105 muon antineutrino candidates were identified with a purity of more than 96\% using the same neutrino event classifier. Fitting these data to the three flavor neutrino oscillation model, using constraints on \thetaonetwo, \thetaonethree, and \dmsqonetwo from solar and reactor neutrino experiments, the oscillation parameters are measured to be $\sintwothree = 0.57^{+0.04}_{-0.03}$, $\dmsqthreetwo = \SI[parse-numbers= false]{+ 2.41\pm0.07 \times 10^{-3}}{eV^2}$, and $\dcp=0.82^{+0.27}_{-0.87}\pi$ with a preference for the normal neutrino mass hierarchy. Leading systematic uncertainties for these measurements come from detector calibration and neutrino interaction models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cluster cosmology without cluster finding

ABSTRACT We propose that observations of supermassive galaxies contain cosmological statistical constraining power similar to conventional cluster cosmology, and we provide promising indications that the associated systematic errors are comparably easier to control. We consider a fiducial spectroscopic and stellar mass complete sample of galaxies drawn from the Dark Energy Spectroscopic Instrument (DESI) and forecast how constraints on Ωm–σ8 from this sample will compare with those from number counts of clusters based on richness λ. At fixed number density, we find that massive galaxies offer similar constraints to galaxy clusters. However, a mass-complete galaxy sample from DESI has the potential to probe lower halo masses than standard optical cluster samples (which are typically limited to λ ≳ 20 and Mhalo ≳ 1013.5 M⊙ h−1); additionally, it is straightforward to cleanly measure projected galaxy clustering wp for such a DESI sample, which we show can substantially improve the constraining power on Ωm. We also compare the constraining power of M*-limited samples to those from larger but mass-incomplete samples [e.g. the DESI Bright Galaxy Survey (BGS) sample]; relative to a lower number density M*-limited samples, we find that a BGS-like sample improves statistical constraints by 60 per cent for Ωm and 40 per cent for σ8, but this uses small-scale information that will be harder to model for BGS. Our initial assessment of the systematics associated with supermassive galaxy cosmology yields promising results. The proposed samples have a ∼10 per cent satellite fraction, but we show that cosmological constraints may be robust to the impact of satellites. These findings motivate future work to realize the potential of supermassive galaxies to probe lower halo masses than richness-based clusters and to potentially avoid persistent systematics associated with optical cluster finding.

(cosmology): large-scale structure of Universe↗

TYPE Ia SUPERNOVA COLORS AND EJECTA VELOCITIES: HIERARCHICAL BAYESIAN REGRESSION WITH NON-GAUSSIAN DISTRIBUTIONS

We investigate the statistical dependence of the peak intrinsic colors of Type Ia supernovae (SNe Ia) on their expansion velocities at maximum light, measured from the Si II λ6355 spectral feature. We construct a new hierarchical Bayesian regression model, accounting for the random effects of intrinsic scatter, measurement error, and reddening by host galaxy dust, and implement a Gibbs sampler and deviance information criteria to estimate the correlation. The method is applied to the apparent colors from BVRI light curves and Si II velocity data for 79 nearby SNe Ia. The apparent color distributions of high-velocity (HV) and normal velocity (NV) supernovae exhibit significant discrepancies for B – V and B – R, but not other colors. Hence, they are likely due to intrinsic color differences originating in the B band, rather than dust reddening. The mean intrinsic B – V and B – R color differences between HV and NV groups are 0.06 ± 0.02 and 0.09 ± 0.02 mag, respectively. A linear model finds significant slopes of –0.021 ± 0.006 and –0.030 ± 0.009 mag (10{sup 3} km s{sup –1}){sup –1} for intrinsic B – V and B – R colors versus velocity, respectively. Because the ejecta velocity distribution is skewed toward high velocities, these effects imply non-Gaussian intrinsic color distributions with skewness up to +0.3. Accounting for the intrinsic-color-velocity correlation results in corrections to A{sub V} extinction estimates as large as –0.12 mag for HV SNe Ia and +0.06 mag for NV events. Velocity measurements from SN Ia spectra have the potential to diminish systematic errors from the confounding of intrinsic colors and dust reddening affecting supernova distances.

79 ASTRONOMY AND ASTROPHYSICS↗

Lattice QCD Calculation of Nucleon Structure

It is emphasized in the 2015 NSAC Long Range Plan that "understanding the structure of hadrons in terms of QCD's quarks and gluons is one of the central goals of modern nuclear physics." Over the last three decades, lattice QCD has developed into a powerful tool for ab initio calculations of strong-interaction physics. Up until now, it is the only theoretical approach to solving QCD with controlled statistical and systematic errors. Since 1985, we have proposed and carried out first-principles calculations of nucleon structure and hadron spectroscopy using lattice QCD which entails both algorithmic development and large-scale computer simulation. We started out by calculating the nucleon form factors -- electromagnetic, axial-vector, πNN, and scalar form factors, the quark spin contribution to the proton spin, the strangeness magnetic moment, the quark orbital angular momentum, the quark momentum fraction, and the quark and glue decomposition of the proton momentum and angular momentum. The first round of calculations were done with Wilson fermions in the `quenched' approximation where the dynamical effects of the quarks in the sea are not taken into account in the Monte Carlo simulation to generate the background gauge configurations. Beginning in 2000, we have started implementing the overlap fermion formulation into the spectroscopy and structure calculations. This is mainly because the overlap fermion honors chiral symmetry as in the continuum. It is going to be more and more important to take the symmetry into account as the simulations move closer to the physical point where the u and d quark masses are as light as a few MeV only. We began with lattices which have quark masses in the sea corresponding to a pion mass at ~ 300 MeV and obtained the strange form factors, charm and strange quark masses, the charmonium spectrum and the D s meson decay constant f D s , the strangeness and charmness, the meson mass decomposition and the strange quark spin from the anomalous Ward identity. Recently, we have started to include multiple lattices with different lattice spacings and different volumes including large lattices at the physical pion mass point. We are getting quite close to being able to calculate the hadron structure at the physical point and to do the continuum and large volume extrapolations, which is our ultimate aim. We have now finished several projects which have included these systematic corrections. They include the leptonic decay width of the ρ, the πN sigma and strange sigma terms, and the strange quark magnetic moment. Over the years, we have also studied hadron spectroscopy with lattice calculations and in phenomenology. These include Roper resonance, pentaquark state, charmonium spectrum, glueballs, scalar mesons a 0 (1450) and σ(600) and other scalar mesons, and the 1 -+ meson. In addition, we have employed the canonical approach to explore the first-order phase transition and the critical point at finite density and finite temperature. We have also discovered a new parton degree of freedom -- the connected sea partons, from the path-integral formulation of the hadronic tensor, which explains the experimentally observed Gottfried sum rule violation. Combining experimental result on the strange parton distribution, the CT10 global fitting results of the total u and d anti-partons and the lattice result of the ratio of the momentum fraction of the strange vs that of u or d in the disconnected insertion, we have shown that the connected sea partons can be isolated. In this final technical report, we shall present a few representative highlights that have been achieved in the project.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Kiloampere, Variable-Temperature, Critical-Current Measurements of High-Field Superconductors

We review variable-temperature, transport critical-current ($I_\text{c}$) measurements made on commercial superconductors over a range of critical currents from less than 0.1 A to about 1 kA. We have developed and used a number of systems to make these measurements over the last 15 years. Two exemplary variable-temperature systems with coil sample geometries will be described: a probe that is only variable-temperature and a probe that is variable-temperature and variable-strain. The most significant challenge for these measurements is temperature stability, since large amounts of heat can be generated by the flow of high current through the resistive sample fixture. Therefore, a significant portion of this review is focused on the reduction of temperature errors to less than ±0.05 K in such measurements. A key feature of our system is a pre-regulator that converts a flow of liquid helium to gas and heats the gas to a temperature close to the target sample temperature. The pre-regulator is not in close proximity to the sample and it is controlled independently of the sample temperature. This allows us to independently control the total cooling power, and thereby fine tune the sample cooling power at any sample temperature. The same general temperature-control philosophy is used in all of our variable temperature systems, but the addition of another variable, such as strain, forces compromises in design and results in some differences in operation and protocol. These aspects are analyzed to assess the extent to which the protocols for our systems might be generalized to other systems at other laboratories. Our approach to variable-temperature measurements is also placed in the general context of measurement-system design, and the perceived advantages and disadvantages of design choices are presented. To verify the accuracy of the variable-temperature measurements, we compared critical-current values obtained on a specimen immersed in liquid helium (“liquid” or $I_{\text{c}}$ $_{\text{liq}}$) at 5 K to those measured on the same specimen in flowing helium gas (“gas” or Ic gas) at the same temperature. These comparisons indicate the temperature control is effective over the superconducting wire length between the voltage taps, and this condition is valid for all types of sample investigated, including Nb-Ti, Nb 3 Sn, and MgB 2 wires. The liquid/gas comparisons are used to study the variable-temperature measurement protocol that was necessary to obtain the “correct” critical current, which was assumed to be the $I_{\text{c}}$ $_{\text{liq}}$. We also calibrated the magnetoresistance effect of resistive thermometers for temperatures from 4 K to 35 K and magnetic fields from 0 T to 16 T. This calibration reduces systematic errors in the variable-temperature data, but it does not affect the liquid/gas comparison since the same thermometers are used in both cases.

47 OTHER INSTRUMENTATION↗

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models. Current status and future directions

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and loss from soil accounts for a large pro portion of land-atmosphere C exchange. Due to large pool size and variable residence time from years to millennia, even small changes in soil organic C(SOC) have substantial effects on the terrestrial C budget, thereby affecting atmospheric carbon dioxide (CO2)concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain and identifying major driving forces controlling soil C storage and fluxes remains a key research challenge his study has compiled century-long (1901-2010)estimates of SOC storage and heterotrophic respiration (Rh) from ten terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project (MsTMIP) and two observation based datasets. The ten-TBM ensemble shows that global SOC estimate range from 4 to 2111 Pg C (1 Pg = 10 15 g) with a median value of 1158 Pg C33 in 2010. Modeling approach estimates a broad range of Rh from 35 to 69 Pg C yr -1 with a median value of 51Pg C yr -1 during 200–2010. The largest uncertainty in SOC stocks exists in the 40–65°N latitude band while Rh differences are the largest in the tropics. All the models agreed that climate and land use changes have decreased SOC stocks while elevated CO 2 and atmospheric nitrogen deposition have increased SOC stocks though the response varied significantly among models. Model representations of temperature and moisture sensitivity,nutrient limitation and land use partially explain the divergent estimates of global SOC stocks and soil fluxes in this study. In addition, major sources of uncertainty from model estimation include exclusion of SOC storage in wetlands and peatlands as well as C storage in deep soil layers.

54 ENVIRONMENTAL SCIENCES↗

CMIP5 multi-model hindcasts for the mid-1970s shift and early 2000s hiatus and predictions for 2016-2035: Meehl and Teng: Initialized hindcasts and predictions

Compared to uninitialized climate change projections, a multi-model ensemble from the CMIP5 10 year decadal prediction experiments produces more warming during the mid-1970s climate shift and less warming in the early 2000s hiatus in both the tropical Indo-Pacific region and globally averaged surface air temperature (TAS) in closer agreement with observations. Assuming bias in TAS has stabilized in the 10 year predictions, after bias adjustment, TAS anomalies for the 2016–2035 period in the 30 year predictions initialized in 2006 are about 16% less than the uninitialized projections. One contributing factor for the improved climate simulation is the bias adjustment, which corrects the models’ systematic errors and higher-than-observed decadal warming trend. Another important factor is the initialization with observations which constrains the ocean such that the starting points of the initialized simulations are close to the observed initial states.

54 ENVIRONMENTAL SCIENCES↗

Trends and sensitivities of low streamflow extremes to discharge timing and magnitude in Pacific Northwest mountain streams

Abstract Path analyses of historical streamflow data from the Pacific Northwest indicate that the precipitation amount has been the dominant control on the magnitude of low streamflow extremes compared to the air temperature‐affected timing of snowmelt runoff. The relative sensitivities of low streamflow to precipitation and temperature changes have important implications for adaptation planning because global circulation models produce relatively robust estimates of air temperature changes but have large uncertainties in projected precipitation amounts in the Pacific Northwest U.S. Quantile regression analyses indicate that low streamflow extremes from the majority of catchments in this study have declined from 1948 to 2013, which may significantly affect terrestrial and aquatic ecosystems, and water resource management. Trends in the 25th percentile of mean annual streamflow have declined and the center of timing has occurred earlier. We quantify the relative influences of total precipitation and air temperature on the annual low streamflow extremes from 42 stream gauges using mean annual streamflow as a proxy for precipitation amount effects and streamflow center of timing as a proxy for temperature effects on low flow metrics, including 7q10 summer (the minimum 7 day flow during summer with a 10 year return period), mean August , mean September , mean summer , 7q10 winter , and mean winter flow metrics. These methods have the benefit of using only readily available streamflow data, which makes our results robust against systematic errors in high elevation distributed precipitation data. Winter low flow metrics are weakly tied to both mean annual streamflow and center of timing .

Kormos, Patrick R.↗

Evidence of Free-Bound Transitions in Warm Dense Matter and Their Impact on Equation-of-State Measurements

Warm dense matter (WDM) is now routinely created and probed in laboratories around the world, providing unprecedented insights into conditions achieved in stellar atmospheres, planetary interiors, and inertial confinement fusion experiments. However, the interpretation of these experiments is often filtered through models with systematic errors that are difficult to quantify. Due to the simultaneous presence of quantum degeneracy and thermal excitation, transitions in which free electrons are de-excited into thermally unoccupied bound states transferring momentum and energy to a scattered x-ray photon become viable. Here we show that such free-bound transitions are a particular feature of WDM and vanish in the limits of cold and hot temperatures. The inclusion of these transitions into the analysis of recent X-ray Thomson Scattering experiments on WDM at the National Ignition Facility and the Linac Coherent Light Source is required to obtain a physically consistent temperature from the Chihara decomposition. This interpretation is corroborated by agreement with a recently developed model-free thermometry technique and presents an important step for precisely characterizing and understanding the complex WDM state of matter.

Boehme, M [HZDR Dresden]↗

Evaluation of energy balance closure adjustment methods by independent evapotranspiration estimates from lysimeters and hydrological simulations

Here, we report non-closure of the surface energy balance is a frequently observed phenomenon of hydrometeorological field measurements, when using the eddy-covariance method, which can be ascribed to an underestimation of the turbulent fluxes. Several approaches have been proposed in order to adjust the measured fluxes for this apparent systematic error. However, there are uncertainties about partitioning of the energy balance residual between the sensible and latent heat flux and whether such a correction should be applied on 30-minute data or longer time scales. The data for this study originate from two grassland sites in southern Germany, where measurements from weighable lysimeters are available as reference. The adjusted evapotranspiration rates are also compared with joint energy and water balance simulations using a physically-based distributed hydrological model. We evaluate two adjustment methods: the first one preserves the Bowen ratio and the correction factor is determined on a daily basis. The second one attributes a smaller portion of the residual energy to the latent heat flux than to the sensible heat flux for closing the energy balance for every 30-minute flux integration interval. Both methods lead to an improved agreement of the eddy-covariance based fluxes with the independent lysimeter estimates and the physically-based model simulations. The first method results in a better comparability of evapotranspiration rates, and the second method leads to a smaller overall bias. These results are similar between both sites despite considerable differences in terrain complexity and grassland management. Moreover, we found that a daily adjustment factor leads to less scatter than a complete partitioning of the residual for every half-hour time interval. Lastly, the vertical temperature gradient in the surface layer and friction velocity were identified as important predictors for a potential future parameterization of the energy balance residual.

58 GEOSCIENCES↗