Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic parameter estimation”

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 721 records · Page 40

Self-interacting asymmetric dark matter coupled to a light massive dark photon

Dark matter (DM) with sizeable self-interactions mediated by a light species offers a compelling explanation of the observed galactic substructure; furthermore, the direct coupling between DM and a light particle contributes to the DM annihilation in the early universe. If the DM abundance is due to a dark particle-antiparticle asymmetry, the DM annihilation cross-section can be arbitrarily large, and the coupling of DM to the light species can be significant. We consider the case of asymmetric DM interacting via a light (but not necessarily massless) Abelian gauge vector boson, a dark photon. In the massless dark photon limit, gauge invariance mandates that DM be multicomponent, consisting of positive and negative dark ions of different species which partially bind in neutral dark atoms. We argue that a similar conclusion holds for light dark photons; in particular, we establish that the multi-component and atomic character of DM persists in much of the parameter space where the dark photon is sufficiently light to mediate sizeable DM self-interactions. We discuss the cosmological sequence of events in this scenario, including the dark asymmetry generation, the freeze-out of annihilations, the dark recombination and the phase transition which gives mass to the dark photon. We estimate the effect of self-interactions in DM haloes, taking into account this cosmological history. We place constraints based on the observed ellipticity of large haloes, and identify the regimes where DM self-scattering can affect the dynamics of smaller haloes, bringing theory in better agreement with observations. Moreover, we estimate the cosmological abundance of dark photons in various regimes, and derive pertinent bounds.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-year incubation experiments boost confidence in model projections of long-term soil carbon dynamics

Abstract Global soil organic carbon (SOC) stocks may decline with a warmer climate. However, model projections of changes in SOC due to climate warming depend on microbially-driven processes that are usually parameterized based on laboratory incubations. To assess how lab-scale incubation datasets inform model projections over decades, we optimized five microbially-relevant parameters in the Microbial-ENzyme Decomposition (MEND) model using 16 short-term glucose (6-day), 16 short-term cellulose (30-day) and 16 long-term cellulose (729-day) incubation datasets with soils from forests and grasslands across contrasting soil types. Our analysis identified consistently higher parameter estimates given the short-term versus long-term datasets. Implementing the short-term and long-term parameters, respectively, resulted in SOC loss (–8.2 ± 5.1% or –3.9 ± 2.8%), and minor SOC gain (1.8 ± 1.0%) in response to 5 °C warming, while only the latter is consistent with a meta-analysis of 149 field warming observations (1.6 ± 4.0%). Comparing multiple subsets of cellulose incubations (i.e., 6, 30, 90, 180, 360, 480 and 729-day) revealed comparable projections to the observed long-term SOC changes under warming only on 480- and 729-day. Integrating multi-year datasets of soil incubations (e.g., > 1.5 years) with microbial models can thus achieve more reasonable parameterization of key microbial processes and subsequently boost the accuracy and confidence of long-term SOC projections.

54 ENVIRONMENTAL SCIENCES↗

Multimodal sensor fusion for real-time standoff estimation in directed energy deposition

In Laser Powder-based Direct Energy Deposition (LP-DED) systems, achieving consistency, precision and quality of produced parts requires tight control over printing parameters. One of the critical parameters is the standoff distance. Maintaining an optimal standoff height is crucial for achieving correct laser power density and powder catchment efficiency, as both laser and powder streams are focused at this distance. Here, this study introduces a novel approach using multimodal sensor fusion to predict standoff height in real-time. The proposed system integrates two low-profile, cost-effective sensors: an RGB coaxial camera and a high frequency and high dynamic range microphone. By utilizing a simple fully connected neural network, trained on a limited dataset, data fusion of these sensors allowed for the real-time prediction of the standoff height. The results demonstrate high resolution and accuracy of the predictions across multiple geometries and a wide range of standoff heights. This approach offers a simple, and cost-effective solution for real-time standoff height monitoring and lays the groundwork for future integration into commercial LP-DED systems.

42 ENGINEERING↗

A finite scale model for shock structure

In this paper we explore the use of a finite scale model of fluid dynamics to predict the finite structure of a shock wave in a perfect gas. We begin by documenting the history and issues that have arisen when Navier–Stokes theory is applied to the shock structure problem, and continue by motivating the improvement that finite scale theory might provide from its representation of inviscid (anomalous) dissipation. Furthermore, our primary results include the formulation of a traveling wave equation from finite scale theory, an analysis of the solutions of that equation as regards shock width and monotonicity properties, and an estimation of the critical parameter of the theory from experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Discovery of a faint outer halo milky way star cluster in the southern sky

We report the discovery of a new, low-luminosity star cluster in the outer halo of the Milky Way. High-quality gr photometry is presented, from which a color–magnitude diagram is constructed, and estimates of age, [Fe/H], [α/Fe], and distance are derived. The star cluster, which we designate as Kim 2, lies at a heliocentric distance of ∼105 kpc. With a half-light radius of ∼12.8 pc and ellipticity of ϵ ∼ 0.12, it shares the properties of outer halo globular clusters, except for at higher metallicity ([Fe/H] ∼ −1.0) and lower luminosity (M{sub V}∼−1.5). These parameters are similar to those for the globular cluster AM 4, which is considered to be associated with the Sagittarius dwarf spheroidal galaxy. We find evidence of dynamical mass segregation and the presence of extra-tidal stars that suggests that Kim 2 is most likely a star cluster. Spectroscopic observations for radial-velocity membership and chemical abundance measurements are needed to further understand the nature of the object.

79 ASTRONOMY AND ASTROPHYSICS↗

DynaShape: Dynamic Shaping of Grid Response with Inverters

We consider the problem of controlling the frequency of low-inertia power systems via inverter-based resources (IBRs) that are weakly connected to the grid. We propose a novel grid-forming control strategy, the so-called frequency shaping control, that aims to shape the frequency response of synchronous generators (SGs) to load perturbations so as to efficiently arrest sudden frequency drops. Our solution relaxes several existing assumptions in the literature and is able to navigate trade-offs between peak power requirements and maximum frequency deviations. Finally, we analyze the robustness to imperfect knowledge of network parameters, while particularly highlighting the importance of accurate estimation of these parameters.

frequency-shaping control↗

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Diffusion and perfusion MR parameters to assess preoperative short-course radiotherapy response in locally advanced rectal cancer: a comparative explorative study among Standardized Index of Shape by DCE-MRI, intravoxel incoherent motion- and diffusion kurtosis imaging-derived parameters

Purpose: To assess preoperative short-course radiotherapy (SCR) tumor response in locally advanced rectal cancer (LARC) by means of Standardized Index of Shape (SIS) by dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI), apparent diffusion coefficient (ADC), intravoxel incoherent motion (IVIM) and diffusion kurtosis imaging (DKI) parameters derived from diffusion-weighted MRI (DW-MRI). Materials and methods: Thirty-four patients with LARC who underwent MRI scans before and after SCR followed by delayed surgery, retrospectively, were enrolled. SIS, ADC, IVIM parameters [tissue diffusion (D{sub t}), pseudo-diffusion (D{sub p}), perfusion fraction (f{sub p})] and DKI parameters [mean diffusivity (MD), mean of diffusional kurtosis (MK)] were calculated for each patient. IVIM parameters were estimated using two methods, namely conventional biexponential fitting (CBFM) and variable projection (VARPRO). After surgery, the pathological TNM and tumor regression grade (TRG) were estimated. For each parameter, percentage changes between before and after SCR were evaluated. Furthermore, an artificial neural network was trained for outcome prediction. Nonparametric sample tests and receiver operating characteristic curve (ROC) analysis were performed. Results: Fifteen patients were classified as responders (TRG ≤ 2) and 19 as not responders (TRG > 3). Seven patients had TRG 1 (pathological complete response, pCR). Mean and standard deviation values of pre-treatment CBFM D{sub p} and mean value of VARPRO D{sub p} pre-treatment showed statistically significant differences to predict pCR. (p value at Mann–Whitney test was 0.05, 0.03 and 0.008, respectively.) Exclusively SIS percentage change showed significant differences between responder and non-responder patients after SCR (p value << 0.001) and to assess pCR after SCR (p value << 0.001). The best results to predict pCR were obtained by VARPRO Fp mean value pre-treatment with area under ROC of 0.84, a sensitivity of 96.4%, a specificity of 71.4%, a positive predictive value (PPV) of 92.9%, a negative predictive value (NPV) of 83.3% and an accuracy of 91.2%. The best results to assess after treatment complete pathological response were obtained by SIS with an area under ROC of 0.89, a sensitivity of 85.7%, a specificity of 92.6%, a PPV of 75.0%, a NPV of 96.1% and an accuracy of 91.2%. Moreover, the best results to differentiate after treatment responders vs. non-responders were obtained by SIS with an area under ROC of 0.94, a sensitivity of 93.3%, a specificity of 84.2%, a PPV of 82.4%, a NPV of 94.1% and an accuracy of 88.2%. Promising initial results were obtained using a decision tree tested with all ADC, IVIM and DKI extracted parameter: we reached high accuracy to assess pathological complete response after SCR in LARC (an accuracy of 85.3% to assess pathological complete response after SCR using VARPRO D{sub p} mean value post-treatment, ADC standard deviation value pre-treatment, MD standard deviation value post-treatment). Conclusion: SIS is a hopeful DCE-MRI angiogenic biomarker to assess preoperative treatment response after SCR with delayed surgery. Furthermore, an important prognostic role was obtained by VARPRO F{sub p} mean value pre-treatment and by a decision tree composed by diffusion parameters derived by DWI and DKI to assess pathological complete response.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Development and validation of a second-order thermal network model for residential buildings

Heating, Ventilation, and Air Conditioning (HVAC) systems can maintain the space air temperature of residential buildings, either directly by heating/cooling the air, or indirectly via heat transfer to and from the building structure that acts as a thermal mass. Hence, HVAC systems can help achieve load shifting, peak load reduction, and/or energy cost saving, thus enabling grid-interactive HVAC operation. A home thermal model that can accurately reflect the dynamics of the space air and interior wall surface temperatures, is therefore valuable. This paper develops such a model using the standard RC (resistance-capacitance) approach. The model contains a virtual envelope node and an internal space node and is thus second-order. A hybrid parameter identification scheme, made up of the least-squares and optimal search methods, is also developed. The proposed model and scheme were validated using data collected from a test home. It was found that a modest amount of training data was sufficient to yield reliable parameter estimates and accurate prediction. It was also found that when making 24-hour-ahead prediction of the space air temperature, both methods had comparable performances when the training data began in a transition season. However, when they began in an HVAC season, the optimal search method performed better. Furthermore, the least-squares method is recommended during a transition season due to its lower computational burden, while the optimal search method is recommended during an HVAC season due to its better estimation performance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Scale-up of microwave-assisted, continuous flow, liquid phase reactors: Application to 5-Hydroxymethylfurfural production

Microwave (MW) technology can be powerful for the electrification and intensification of chemical manufacturing. However, very few examples of scaled MW processes exist. In this work, we build a scaled, MW-assisted, continuous flow reactor for the processing of liquid phase chemistries with specific application to 5-hydroxymethylfurfural (HMF) production. Here, we demonstrate a corresponding computational fluid dynamics model to simulate the reactor’s temperature profile and performance, both of which are in good agreement with experiments. We construct a surrogate model to relate operating parameters to performance for guiding experiments, without demanding simulations, via active learning. Heat recirculation is demonstrated in this scaled reactor, further extending its scale and reducing the energy demand. An HMF yield of ~55 % and a productivity of 0.1 kg/hr, 8x higher than any other reactor, at a flowrate >20x than prior work, are demonstrated while maintaining energy efficiency of >98 %. A basic economic analysis estimates a $1.85/kg cost of HMF from fructose and >60 % reduction in CO 2 emissions than a conventional system without considering the impact of green electricity. We demonstrate that MW technology is well-applied to systems larger than the laboratory scale, and strategies to further reduction in cost and CO 2 are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THE EVIL-MC MODEL FOR ELLIPSOIDAL VARIATIONS OF PLANET-HOSTING STARS AND APPLICATIONS TO THE HAT-P-7 SYSTEM

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a {approx} Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 {+-} 0.06) Multiplication-Sign 10{sup -3}, and we have revised the planet's dayside brightness temperature to 2680{sup +10}{sub -20} K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates that this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

79 ASTRONOMY AND ASTROPHYSICS↗

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, ( x ) u - d , helicity moment, ( x ) Δ u - Δ d , and the transversity moment, ( x ) δu - δd , of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass M π . The final results are obtained using a simultaneous fit in the lattice spacing a , pion mass M π and the finite volume parameter M π L keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the $ \overline{\mathrm{MS}} $ scheme at 2 GeV, are ( x ) u - d = 0 . 155(17)(20), ( x ) Δ u - Δ d = 0 . 183(14)(20) and ( x ) δu - δd = 0 . 220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Evidence for pressure induced unconventional quantum criticality in the coupled spin ladder antiferromagnet C 9 H 18 N 2 CuBr 4

Quantum phase transitions in quantum matter occur at zero temperature between distinct ground states by tuning a nonthermal control parameter. Often, they can be accurately described within the Landau theory of phase transitions, similarly to conventional thermal phase transitions. However, this picture can break down under certain circumstances. Here, we present a comprehensive study of the effect of hydrostatic pressure on the magnetic structure and spin dynamics of the spin-1/2 ladder compound C 9 H 18 N 2 CuBr 4 . Single-crystal heat capacity and neutron diffraction measurements reveal that the Néel-ordered phase breaks down beyond a critical pressure of Pc~1.0 GPa through a continuous quantum phase transition. Estimates of the critical exponents suggest that this transition may fall outside the traditional Landau paradigm. The inelastic neutron scattering spectra at 1.3 GPa are characterized by two well-separated gapped modes, including one continuum-like and another resolution-limited excitation in distinct scattering channels, which further indicates an exotic quantum-disordered phase above Pc.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗

Entropy and Gravitation—From Black Hole Computers to Dark Energy and Dark Matter

We show that the concept of entropy and the dynamics of gravitation provide the linchpin in a unified scheme to understand the physics of black hole computers, spacetime foam, dark energy, dark matter and the phenomenon of turbulence. We use three different methods to estimate the foaminess of spacetime, which, in turn, provides a back-door way to derive the Bekenstein-Hawking formula for black hole entropy and the holographic principle. Generalizing the discussion for a static spacetime region to the cosmos, we find a component of dark energy (resembling an effective positive cosmological constant of the correct magnitude) in the current epoch of the universe. The conjunction of entropy and gravitation is shown to give rise to a phenomenological model of dark matter, revealing the natural emergence, in galactic and cluster dynamics, of a critical acceleration parameter related to the cosmological constant; the resulting mass profiles are consistent with observations. Unlike ordinary matter, the quanta of the dark sector are shown to obey infinite statistics. This property of dark matter may lead to some non-particle phenomenology and may explain why dark matter particles have not been detected in dark matter search experiments. We also show that there are deep similarities between the problem of “quantum gravity” (more specifically, the holographic spacetime foam) and turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

MASS RATIO DETERMINATION FROM H{sub α} LINES IN BLACK HOLE X-RAY TRANSIENTS

We find that the mass ratio q in quiescent black hole (BH) X-ray transients is tightly correlated with the ratio of the double-peak separation (DP) to the full width half maximum (FWHM) of the H{sub α} emission line, logq=−6.88−23.2 log(DP/FWHM). This correlation is explained through the efficient truncation of the outer disk radius by the 3:1 resonance with the companion star. This is the dominant tidal interaction for extreme mass ratios q=M{sub 2}/M{sub 1}≲0.25, the realm of BH (and some neutron star) X-ray transients. Mass ratios can thus be estimated with a typical uncertainty of ≈32%, provided that the H{sub α} profile used to measure DP/FWHM is an orbital phase average. We apply the DP/FWHM–q relation to the three faint BH transients XTE J1650–500, XTE J1859+226, and Swift J1357–0933 and predict q=0.026{sub −0.007}{sup +0.038}, 0.049{sub −0.012}{sup +0.023} and 0.040{sub −0.005}{sup +0.003}, respectively. This new relation, together with the FWHM--K{sub 2} correlation presented in Paper I, allows the extraction of fundamental parameters from very faint targets and, therefore, the extension of dynamical BH studies to much deeper limits than was previously possible. As an example, we combine our mass ratio determination for Swift J1357–0933 with previous reported values to yield a BH mass of 12.4 ± 3.6 M {sub ⊙}. This confirms Swift J1357–0933 as one of the most massive BH low-mass X-ray binaries in the Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Linking big models to big data: efficient ecosystem model calibration through Bayesian model emulation

Data-model integration plays a critical role in assessing and improving our capacity to predict ecosystem dynamics. Similarly, the ability to attach quantitative statements of uncertainty around model forecasts is crucial for model assessment and interpretation and for setting field research priorities. Bayesian methods provide a rigorous data assimilation framework for these applications, especially for problems with multiple data constraints. However, the Markov chain Monte Carlo (MCMC) techniques underlying most Bayesian calibration can be prohibitive for computationally demanding models and large datasets. We employ an alternative method, Bayesian model emulation of sufficient statistics, that can approximate the full joint posterior density, is more amenable to parallelization, and provides an estimate of parameter sensitivity. Analysis involved informative priors constructed from a meta-analysis of the primary literature and specification of both model and data uncertainties, and it introduced novel approaches to autocorrelation corrections on multiple data streams and emulating the sufficient statistics surface. We report the integration of this method within an ecological workflow management software, Predictive Ecosystem Analyzer (PEcAn), and its application and validation with two process-based terrestrial ecosystem models: SIPNET and ED2. In a test against a synthetic dataset, the emulator was able to retrieve the true parameter values. A comparison of the emulator approach to standard brute-force MCMC involving multiple data constraints showed that the emulator method was able to constrain the faster and simpler SIPNET model's parameters with comparable performance to the brute-force approach but reduced computation time by more than 2 orders of magnitude. The emulator was then applied to calibration of the ED2 model, whose complexity precludes standard (brute-force) Bayesian data assimilation techniques. Both models are constrained after assimilation of the observational data with the emulator method, reducing the uncertainty around their predictions. Performance metrics showed increased agreement between model predictions and data. Our study furthers efforts toward reducing model uncertainties, showing that the emulator method makes it possible to efficiently calibrate complex models.

58 GEOSCIENCES↗