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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Statistical framework to assess long-term spatio-temporal climate changes: East River mountainous watershed case study

Abstract Evaluation of long-term temporal and spatial climatic change in mountainous regions is a critical challenge because of the interactive effects of multiple land and climatic factors and processes. Here we present the application of the statistical framework to the assessment of changes of climatic conditions, using data from 17 meteorological stations across the East River watershed near Crested Butte, Colorado, USA, and spanning the period from 1966 to 2021. The framework is developed based on (1) a time-series analysis of daily, monthly, and yearly averaged meteorological parameters (temperature, relative humidity, precipitation, wind speed, etc.), (2) evaluation and time series analysis of potential evapotranspiration (ET o ), actual evapotranspiration (ET), aridity index (AI), standard precipitation index (SPI) and standard precipitation-evapotranspiration index (SPEI), and (3) a temporal-spatial climatic zonation of the studied area based on the hierarchical clustering and PCA analysis of the SPEI, because the SPEI can be considered an integrative characteristic of the changes of climatic conditions. The Budyko model, with the application of the Penman–Monteith equation for the estimation of ET o , was used to determine the ET. The time series analysis of the AI is used to identify the periods with energy limited and water limited conditions. Hierarchical clustering of site locations for the three temporal segments of the SPEI showed a significant temporal-spatial shifts, indicating that dynamic climatic processes drive zonation patterns. Therefore, the watershed climatic zonation requires periodic re-evaluation based on the structural time series analysis of meteorological and water balance data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Precambrian Crystalline Basement Properties From Pressure History Matching and Implications for Induced Seismicity in the US Midcontinent

Injection-induced seismicity across the US midcontinent has almost exclusively occurred in the crystalline basement that underlies the Arbuckle Group aquifer and its equivalents, the primary wastewater disposal zone in this region. However, the properties of the basement are not well known. Newly compiled data, from Class I wells in Kansas, provide a unique record of pressures in the Arbuckle and an opportunity to constrain the reservoir-scale properties of the basement such as permeability, diffusivity, and specific storage. Constraints on these parameters are critical for modeling fluid flow and pressures across the entire Arbuckle-basement system, and are necessary for accurate evaluation and prediction of injection-induced earthquakes. Here, we present a detailed, three-dimensional geological and pressure history-matched numerical model for the Arbuckle and basement, based on data from >400 wells covering a large region in south-central Kansas, where injection-induced seismicity has been concentrated since 2014. Simulations of dynamic data from 319 wells indicate that Arbuckle pressures have increased by 1.1 MPa in high injection rate areas and an overpressure of <0.1 MPa may be the cause of seismicity in the basement. Pressure-history matching also yields the likely range in porosity (0.3%–7%), permeability (0.1–0.7 mD), and diffusivity (0.004–0.07 m2 /s) for the basement. The resulting estimates suggest reservoir-scale properties of the basement are enhanced by faults and fractures. Importantly, the diffusivities determined in this study are lower than estimates derived from Kansas earthquake triggering fronts, and suggest that such seismicity-based techniques may have limitations, particularly where spacetime patterns between injection and seismicity are complex.

58 GEOSCIENCES↗
DOE OSTI2024

Discoveries in Blast-Driven Turbulence of Astrophysical Relevance

The fluid mixing caused by variable-density instabilities is important in a wide variety of scenarios from ocean mixing and astrophysical phenomena to nuclear fusion techniques and atomic weapons. This thesis explores the mixing resulting from a specific instability known as the Blast Driven Instability (BDI). This work investigates the variable density mixing in an explosively driven environment due to the fluid instabilities at the material interfaces. Specifically, diverging Richtmyer-Meshkov (impulsive-acceleration environment) and Rayleigh-Taylor (variable-acceleration environment) instabilities (present in supernova and inertial confinement fusion) are studied using advanced high-speed diagnostics in carefully designed laboratory experiments. The BDI morphology is presented through a time development of Mie scattering images, and steps through the parameter space (varying density ratio and driver speed), highlighting the development of the structures that form during mixing. A scaling criterion is used to relate the two systems of vastly different spatiotemporal scales. Velocity fields in the BDI have been captured for the first time using the high temporal resolution PIV technique. Subsequent analysis of the dynamics of the instability from the velocity fields illustrates the distribution of kinetic energy, the transition to turbulence, and the characteristic growth of the instability are discussed. This study furthers understanding of how blast-driven instability pertains to supernova and inertial confinement fusion science. The morphology of the BDI has been characterized for the first time. This work steps through the parameter space covered in Mie scattering experiments, and how the different parameters contribute to development of structures and mixing. It also examines a scaling of the Atwood number for expanding predictive capabilities to other experimental conditions and simulations. The first collection of velocity fields acquired for the BDI are recorded, and subsequent analysis evaluating the distribution of kinetic energy throughout space and time for two density ratios from the overall parameter space, as well as the transition to turbulence, estimated from a Reynolds number calculated based on momentum mixing are all presented. This information is useful in advancing the development of models to predict physics of high energy density applications where experiments are not always readily available. This research has successfully demonstrated understanding for the time criteria defining regimes where the shock driven (Richtmyer-Meshkov instability) and the buoyancy driven (Rayleigh-Taylor instability) dominates through a parametric study of density variation (Atwood number) and driver speed (Mach number). All this furthers understanding of how the BDI pertains to SN and ICF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evaluating alternative ebullition models for predicting peatland methane emission and its pathways via data–model fusion

Abstract. Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Design stress evaluation based on strain-rate sensitivity analysis for nickel alloys used in the very-high temperature nuclear system

Both Alloy 617 and Alloy 230 have been considered the most promising structural materials for the Very High Temperature Reactor (VHTR). In this study, mechanical properties of both alloys were examined by performing tensile tests at three different strain rates and at temperatures up to 1000 deg.C. This range covers time-dependent (plasticity) to time-independent (creep) deformations. Strain-rate sensitivity analysis for each alloy was conducted in order to approximate the long-term flow stresses. The strain-rate sensitivities for the 0.2% flow stress were found to be temperature independent (m {approx_equal} 0) at temperatures ranging from room temperature to 700 deg.C due to dynamic strain aging. At elevated temperatures (800-1000 deg.C), the strain-rate sensitivity significantly increased (m > 0.1). Compared to Alloy 617, Alloy 230 displayed higher strain-rate sensitivities at high temperatures. This leads to a lower estimated long-term flow stresses. Results of this analysis were used to evaluate current American Society of Mechanical Engineers (ASME) allowable design limits. According to the comparison with the estimated flow stresses, the allowable design stresses in ASME B and PV Code for either alloy did not provide adequate degradation estimation for the possible long-term service life in VHTR. However, rupture stresses for Alloy 617, developed in ASME code case N-47-28, can generally satisfy the safety margin estimated in the study following the strain-rate sensitivity analysis. Nevertheless, additional material development studies might be required, since the design parameters for rupture stresses are constrained such that current VHTR conceptual designs cannot satisfy the limits. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Enhancing T c in a composite superconductor/metal bilayer system: A dynamical cluster approximation study

It has been proposed that the superconducting transition temperature T c of an unconventional superconductor with a large pairing scale but strong phase fluctuations can be enhanced by coupling it to a metal. However, the general efficacy of this approach across different parameter regimes remains an open question. Using the dynamical cluster approximation, we study this question in a system composed of an attractive Hubbard layer in the intermediate coupling regime, where the magnitude of the attractive Coulomb interaction |U| is slightly larger than the bandwidth W, hybridized with a noninteracting metallic layer. We find that while the superconducting transition becomes more mean-field-like with increasing interlayer hopping, the superconducting transition temperature T c exhibits a nonmonotonic dependence on the strength of the hybridization t ⊥ . Furthermore, this behavior arises from a reduction of the effective pairing interaction in the correlated layer that outcompetes the growth in the intrinsic pair-field susceptibility induced by the coupling to the metallic layer. We find that the largest T c inferred here for the composite system is comparable to the maximum value currently estimated for the isolated negative-U Hubbard model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Static and dynamic properties of the frustrated spin-$\frac{1}{2}$ depleted-kagome antiferromagnet Cu 7 (TeO 3 ) 2 (SO 4 ) 2 (OH) 6

The structural and magnetic properties of the two-dimensional spin-$\frac{1}{2}$ depleted-kagome compound Cu 7 (TeO 3 ) 2 (SO 4 ) 2 (OH) 6 are investigated using x-ray diffraction, magnetization, heat capacity, and 1 H Nuclear Magnetic Resonance (NMR) measurements. From the analysis of magnetic susceptibility, we found a large Curie-Weiss temperature [𝜃 CW = −50⁢(2)⁢K] and the coexistence of antiferromagnetic and ferromagnetic interactions. The value of 𝜃 CW gives an estimate of the average nearest-neighbor antiferromagnetic interaction of 𝐽/𝑘 B ≃ 66K. The NMR relaxation rates (1/𝑇 1 and 1/𝑇 2 ) exhibit a peak, providing evidence for a magnetic long-range order at 𝑇* ≃ 4K which appears to be canted antiferromagnetic type. Heat capacity also features a broad maximum at 𝑇* that moves towards higher temperatures with increasing magnetic field, reflecting defect induced Schottky anomaly. The frustration parameter 𝑓 𝑟 = |𝜃 CW |/𝑇* ≃ 12.5 renders the compound a highly frustrated low-dimensional magnet.

Akshay, K. U. [Indian Institute of Science Educati↗
DOE OSTI2016

Hierarchical mixture of experts and diagnostic modeling approach to reduce hydrologic model structural uncertainty

In most water resources applications, a single model structure might be inadequate to capture the dynamic multi-scale interactions among different hydrological processes. Calibrating single models for dynamic catchments, where multiple dominant processes exist, can result in displacement of errors from structure to parameters, which in turn leads to over-correction and biased predictions. An alternative to a single model structure is to develop local expert structures that are effective in representing the dominant components of the hydrologic process and adaptively integrate them based on an indicator variable. In this study, the Hierarchical Mixture of Experts (HME) framework is applied to integrate expert model structures representing the different components of the hydrologic process. Various signature diagnostic analyses are used to assess the presence of multiple dominant processes and the adequacy of a single model, as well as to identify the structures of the expert models. The approaches are applied for two distinct catchments, the Guadalupe River (Texas) and the French Broad River (North Carolina) from the Model Parameter Estimation Experiment (MOPEX), using different structures of the HBV model. Furthermore, the results show that the HME approach has a better performance over the single model for the Guadalupe catchment, where multiple dominant processes are witnessed through diagnostic measures. Whereas, the diagnostics and aggregated performance measures prove that French Broad has a homogeneous catchment response, making the single model adequate to capture the response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Three-body capture, ejection, and the demographics of bound objects in binary systems

ABSTRACT We study the capture of light objects of arbitrary velocity by binary systems. Extending results for the capture of comets in the Solar system, we develop a simple geometric characterization of the capture cross-section, leading directly to the distribution of orbital parameters of captured objects. We use the same framework to study the lifetimes of these bound orbits prior to ejection, and find that a simplified version of the Öpik–Arnold approach readily yields a closed-form estimate for the ejection rate that agrees well with numerical experiments. Without any detailed-balance assumptions, our results make manifest the characteristics of close encounters leading to capture and ejection. As an application of our results, we demonstrate the estimation of the equilibrium population of captured dark matter particles in a binary system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

iTOUGH2 v7.1

iTOUGH2 (inverse TOUGH2) provides inverse modeling capabilities for TOUGH2, a simulator for multi-dimensional, multi-phase, multi-component, non-isothermal flow and transport in fractured porous media. iTOUGH2 performs sensitivity analyses, data-worth analyses, parameter estimation, and uncertainty propagation analyses in geosciences and reservoir engineering and other application areas. iTOUGH2 supports a number of different combinations of fluids and components (equation-of-state (EOS) modules). In addition, the optimization routines implemented in iTOUGH2 can also be used for sensitivity analysis, automatic model calibration, and uncertainty quantification of any external code that uses text-based input and output files using the PEST protocol. iTOUGH2 solves the inverse problem by minimizing a non-linear objective function of the weighted differences between model output and the corresponding observations. Multiple minimization algorithms (derivative-free, gradient-based, and second-order; local and global) are available. iTOUGH2 also performs Latin Hypercube Monte Carlo simulations for uncertainty propagation analyses. A detailed residual and error analysis is provided. This upgrade includes (a) global sensitivity analysis methods, (b) dynamic memory allocation (c) additional input features and output analyses, (d) increased forward simulation capabilities, (e) parallel execution on multicore PCs and Linux clusters, and (f) bug fixes.

Finsterle, Stefan↗
DOE OSTI2019

Identification of aggregate building thermal dynamic model and unmeasured internal heat load from data

This paper is on the problem of simultaneously identifying the parameters of an aggregate thermal dynamic model of a multi-zone building and unknown disturbances from input-output data. An aggregate model is a single-zone equivalent of a multi-zone building, and is useful for many purposes, including model based control of large heating, ventilation and air conditioning (HVAC) equipment that delivers thermal energy to the entire building. A key challenge in identification is the presence of unknown disturbance since it is not measurable but non-negligible.We first present a principled method to aggregate a multizone building model into a single zone model. We then provide a method to identify thermal parameters and the unknown disturbance for this aggregate (single-zone) model. Finally, we test our proposed identification algorithm to data generated from a virtual building. A key insight provided by the aggregation method allows us to recognize under what conditions the estimation of the disturbance signal will be necessarily poor and uncertain.

Guo, Zhong↗
DOE OSTI2025

Determination of Constraint-Independent Crack Tip Opening Angle for Stable Crack Growth in High-Strength Ductile Steels

The crack tip opening angle (CTOA) is one of fracture toughness parameters that has been used for decades in describing large stable crack growth in thin-walled aerospace structures under the low-constraint conditions. Recently, the pipeline industry has developed a growing interest in using the CTOA parameter to serve as the minimum required fracture toughness to arrest dynamic crack propagation in modern gas transmission pipelines made of high-strength ductile steel. To meet this industrial need, the CTOA test standard ASTM E3039 was therefore developed for measuring the constant critical CTOA. ASTM E3039 recommends a drop weight tearing test (DWTT) specimen with a shallow crack for standard CTOA testing, but its CTOA may depend on the low constraint condition of the DWTT specimen at the crack tip. Verifying the constraint independence of the DWTT-measured CTOA thus becomes indispensable for applying CTOA toughness to the running fracture control in the pipeline design. For this purpose, the present paper evaluates critical CTOA values in a set of fracture toughness tests on single-edge notched bend (SENB) specimens with shallow and deep cracks, based on four CTOA estimation models. Among these, the Ln(P)-LLD linear fit model is similar to that recommended by ASTM E3039 for CTOA calculation. Fracture test data for X80 pipeline steel and HY80 structural steel were considered in the CTOA evaluation. The results showed that the four CTOA models were able to determine a constraint independent CTOA value for stable crack growth in the SENB specimens. As a result, a single, reliable, constant CTOA value could be determined regardless of the specimen geometry or the crack-tip constraint conditions. Therefore, the CTOA measured using ASTM E3039 is constraint-independent and transferable to use in cases of actual cracks propagating in gas transmission pipelines.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The Innograph irradiations; HTR graphite material properties from 0 to 25 DPA - 18541

In the fifth, sixth, and seventh Framework Program of the European Commission, the Nuclear Research and Consultancy Group (NRG) have conducted a large graphite irradiation program for the development of High Temperature Reactor technology. The graphite irradiation program comprises five irradiations executed in the High Flux Reactor (HFR) in Petten, resulting in design curves at irradiation temperatures between 650 deg. C and 950 deg. C. The neutron doses of these irradiations, called INNOGRAPH, range from 1 dpa to 23 dpa. The irradiation facilities contained present day available graphite grades: NBG-10, NBG-17, NBG-18, NBG-20, NBG-25, PCEA, PPEA, PCIB, LPEB, IG-110 and IG-430. The post-irradiation examination of the INNOGRAPH samples is focused on dimensional changes, dynamic Young s modulus, coefficient of thermal expansion (CTE) and coefficient of thermal conductivity. Irradiation behaviour of HTR candidate graphite grades are compared with respect to manufacturing process and raw coke material. Accuracy and reliability of irradiation parameters and post irradiation examination methods are essential for safety and performance relevant material property determination, in particular for probabilistic material behaviour such as graphite. Considerations on quality requirements with regards to graphite material property determination are therefore included in the paper. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Pore-scale and Continuum Simulations of Solute Transport Micromodel Benchmark Experiments

Four sets of micromodel nonreactive solute transport experiments were conducted with flow velocity, grain diameter, pore-aspect ratio, and flow focusing heterogeneity as the variables. The data sets were offered to pore-scale modeling groups to test their simulators. Each set consisted of two learning experiments, for which all results was made available, and a challenge experiment, for which only the experimental description and base input parameters were provided. The experimental results showed a nonlinear dependence of the dispersion coefficient on the Peclet number, a negligible effect of the pore-aspect ratio on transverse mixing, and considerably enhanced mixing due to flow focusing. Five pore-scale models and one continuum-scale model were used to simulate the experiments. Of the pore-scale models, two used a pore-network (PN) method, two others are based on a lattice-Boltzmann (LB) approach, and one employed a computational fluid dynamics (CFD) technique. The learning experiments were used by the PN models to modify the standard perfect mixing approach in pore bodies into approaches to simulate the observed incomplete mixing. The LB and CFD models used these experiments to appropriately discretize the grid representations. The continuum model use published non-linear relations between transverse dispersion coefficients and Peclet numbers to compute the required dispersivity input values. Comparisons between experimental and numerical results for the four challenge experiments show that all pore-scale models were all able to satisfactorily simulate the experiments. The continuum model underestimated the required dispersivity values and, resulting in less dispersion. The PN models were able to complete the simulations in a few minutes, whereas the direct models needed up to several days on supercomputers to resolve the more complex problems.

benchmark↗
DOE OSTI2014

The density and mass of unshocked ejecta in Cassiopeia a through low frequency radio absorption

Characterizing the ejecta in young supernova remnants is a requisite step toward a better understanding of stellar evolution. In Cassiopeia A the density and total mass remaining in the unshocked ejecta are important parameters for modeling its explosion and subsequent evolution. Low frequency (<100 MHz) radio observations of sufficient angular resolution offer a unique probe of unshocked ejecta revealed via free-free absorption against the synchrotron emitting shell. We have used the Very Large Array plus Pie Town Link extension to probe this cool, ionized absorber at 9'' and 18.''5 resolution at 74 MHz. Together with higher frequency data we estimate an electron density of 4.2 cm{sup –3} and a total mass of 0.39 M {sub ☉} with uncertainties of a factor of ∼2. This is a significant improvement over the 100 cm{sup –3} upper limit offered by infrared [S III] line ratios from the Spitzer Space Telescope. Our estimates are sensitive to a number of factors including temperature and geometry. However using reasonable values for each, our unshocked mass estimate agrees with predictions from dynamical models. We also consider the presence, or absence, of cold iron- and carbon-rich ejecta and how these affect our calculations. Finally we reconcile the intrinsic absorption from unshocked ejecta with the turnover in Cas A's integrated spectrum documented decades ago at much lower frequencies. These and other recent observations below 100 MHz confirm that spatially resolved thermal absorption, when extended to lower frequencies and higher resolution, will offer a powerful new tool for low frequency astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Compartment modeling of dynamic brain PET—The impact of scatter corrections on parameter errors

Purpose: The aim of this study was to investigate the effect of scatter and its correction on kinetic parameters in dynamic brain positron emission tomography (PET) tumor imaging. The 2-tissue compartment model was used, and two different reconstruction methods and two scatter correction (SC) schemes were investigated. Methods: The GATE Monte Carlo (MC) software was used to perform 2 × 15 full PET scan simulations of a voxelized head phantom with inserted tumor regions. The two sets of kinetic parameters of all tissues were chosen to represent the 2-tissue compartment model for the tracer 3′-deoxy-3′-({sup 18}F)fluorothymidine (FLT), and were denoted FLT{sub 1} and FLT{sub 2}. PET data were reconstructed with both 3D filtered back-projection with reprojection (3DRP) and 3D ordered-subset expectation maximization (OSEM). Images including true coincidences with attenuation correction (AC) and true+scattered coincidences with AC and with and without one of two applied SC schemes were reconstructed. Kinetic parameters were estimated by weighted nonlinear least squares fitting of image derived time–activity curves. Calculated parameters were compared to the true input to the MC simulations. Results: The relative parameter biases for scatter-eliminated data were 15%, 16%, 4%, 30%, 9%, and 7% (FLT{sub 1}) and 13%, 6%, 1%, 46%, 12%, and 8% (FLT{sub 2}) for K{sub 1}, k{sub 2}, k{sub 3}, k{sub 4}, V{sub a}, and K{sub i}, respectively. As expected, SC was essential for most parameters since omitting it increased biases by 10 percentage points on average. SC was not found necessary for the estimation of K{sub i} and k{sub 3}, however. There was no significant difference in parameter biases between the two investigated SC schemes or from parameter biases from scatter-eliminated PET data. Furthermore, neither 3DRP nor OSEM yielded the smallest parameter biases consistently although there was a slight favor for 3DRP which produced less biased k{sub 3} and K{sub i} estimates while OSEM resulted in a less biased V{sub a}. The uncertainty in OSEM parameters was about 26% (FLT{sub 1}) and 12% (FLT{sub 2}) larger than for 3DRP although identical postfilters were applied. Conclusions: SC was important for good parameter estimations. Both investigated SC schemes performed equally well on average and properly corrected for the scattered radiation, without introducing further bias. Furthermore, 3DRP was slightly favorable over OSEM in terms of kinetic parameter biases and SDs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Assimilation of multiple datasets results in large differences in regional- to global-scale NEE and GPP budgets simulated by a terrestrial biosphere model

In spite of the importance of land ecosystems in offsetting carbon dioxide emissions released by anthropogenic activities into the atmosphere, the spatiotemporal dynamics of terrestrial carbon fluxes remain largely uncertain at regional to global scales. Over the past decade, data assimilation (DA) techniques have grown in importance for improving these fluxes simulated by terrestrial biosphere models (TBMs), by optimizing model parameter values while also pinpointing possible parameterization deficiencies. Although the joint assimilation of multiple data streams is expected to constrain a wider range of model processes, their actual benefits in terms of reduction in model uncertainty are still under-researched, also given the technical challenges. In this study, we investigated with a consistent DA framework and the ORCHIDEE-LMDz TBM–atmosphere model how the assimilation of different combinations of data streams may result in different regional to global carbon budgets. To do so, we performed comprehensive DA experiments where three datasets (in situ measurements of net carbon exchange and latent heat fluxes, spaceborne estimates of the normalized difference vegetation index, and atmospheric CO 2 concentration data measured at stations) were assimilated alone or simultaneously. We thus evaluated their complementarity and usefulness to constrain net and gross C land fluxes. We found that a major challenge in improving the spatial distribution of the land C sinks and sources with atmospheric CO 2 data relates to the correction of the soil carbon imbalance.

54 ENVIRONMENTAL SCIENCES↗