Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural models”

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 613 records · Page 34

Low global warming potential (GWP) refrigerant supermarket refrigeration system modeling and its application

As an environmentally friendly low global warming potential (GWP) refrigerant, Carbon dioxide (CO 2 ) has continuously gained popularity and research attention as alternative refrigerant for supermarket refrigeration system. In this paper, to fulfill the increasing need of accurate Low-GWP supermarket refrigeration models for development of supervisory level control and optimization strategies, a high fidelity model is developed for CO 2 transcritical supermarket refrigeration system which includes compressor rack of low temperature (LT) compressors and medium temperature (MT) compressors, air-cooled gas cooler, evaporator, expansion valves and other auxiliary equipment. A resistance-capacity model structure is proposed to simulate the display cases. Semi-thermodynamic models are proposed to estimate reciprocating compressors volumetric efficiency and power consumption. Furthermore the zone modeling approach is used for evaporator simulation, and air-cooled gas cooler is modeled with distributed modeling method. The expansion valve simulation is based on orifice flow model. To calibrate these models, both manufacture performance data and experimental data are used. The experiments are conducted with a full instrumental CO 2 supermarket refrigeration system installed in Oak Ridge National Lab Flexible Research Platform (FRP). The simulation model can predict the system performance, including power consumption, cooling capacity, mass flow rate, temperature, and pressure, with high accuracy (within ±4%) compared to experimental data. In addition, this developed model has been used to create the system optimum high side pressure for high side expansion valve control, and to generate wide operating range simulation data for developing the batter-equivalent commercial refrigeration system model which can be used in grid interactive control development.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Atomic-Level Structure of Mesoporous Hexagonal Boron Nitride Determined by High-Resolution Solid-State Multinuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations

Mesoporous hexagonal boron nitride (p-BN) has received significant attention over the last decade as a promising candidate for water cleaning/pollutant removal and hydrogen storage applications. In this work, high-resolution solid-state NMR spectroscopy and plane-wave density-functional theory (DFT) calculations are used to obtain an atomic-level description of p-BN. 1 H– 15 N or 1 H– 14 N heteronuclear (HETCOR) correlation experiments recorded with either conventional NMR at room temperature or dynamic nuclear polarization surface-enhanced spectroscopy (DNP-SENS) at ca. 100 K reveal NB 2 H, NBH 2 , NBH 3 + species residing on the edges of BN sheets. Ultra-high field 35.2 T 11 B NMR spectroscopy was used to resolve 11 B NMR signals from BN 3 , BN 2 O x (OH) 1–x (x = 0–1), BNO x (OH) 2–x (x = 0–2), BO x (OH) 3–x (x = 0–3), and BO x (OH) 4–x – (x = 0–4). Importantly, 2D 11 B dipolar double-quantum–single-quantum homonuclear correlation spectra reveal that many pore/defect sites are composed of boron oxide/hydroxide clusters connected to the BN framework through BN 2 O units. 1D and 2D 11 B{ 15 N} HETCOR NMR experiments, in addition to plane-wave DFT calculations of nine different structural models, further confirm the assignment of all NMR signals. The detailed structure determination of the pore and edge/defect sites within p-BN should further enable the rational design and development of next-generation p-BN-based materials. In addition, the techniques outlined here should be applicable to determine structure within other porous and/or boron-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

BaCu 4/3 Si 2/3 P 2 and BaCu 2–( x + y ) Zn x Si y P 2 : Expanding the Semiconducting Landscape in the ThCr 2 Si 2 -Type Family

ThCr 2 Si 2 -type layered materials are a large family of compounds with applications ranging from thermoelectricity to magnetism, with the vast majority of the members exhibiting metallic behavior. Here, in this study, we synthesized a new group of materials with Cu-Si and Cu-Zn-Si square nets with the general formula BaCu 1.33 Si 0.67 P 2 and BaCu 2–(x+y) Zn x Si y P 2 (0 ≤ x ≤ 0.9; 0.3 ≤ y ≤ 0.7). Several synthesized compounds are charge-balanced semiconductors, which are rare in the ThCr 2 Si 2 family. All the reported compounds crystallize in the ThCr 2 Si 2 -type tetragonal I4/mmm space group, with Cu/Zn/Si jointly occupying the same 4d crystallographic site. In the Zn-free composition, BaCu 1.33 Si 0.67 P 2 , Ba, and P each occupy a single crystallographic site. The introduction of Zn results in the expansion of the unit cell and splitting the Ba atomic sites along the [001] direction. Such structural displacement of the Ba atoms was confirmed by the heat capacity measurements. Band structure and density-of-states calculations on ordered hypothetical structural models reveal either a small bandgap (∼0.2 eV) or semimetallic band structures. The compounds reported here exhibit high Seebeck coefficients and ultralow thermal conductivity, making them promising candidates for the development of thermoelectric materials.

crystal structure↗

Modeling and characterizing anisotropic inclusion orientation in heterogeneous material via directional cluster functions and stochastic microstructure reconstruction

We present a framework to model and characterize the microstructure of heterogeneous materials with anisotropic inclusions of secondary phases based on the directional correlation functions of the inclusions. Specifically, we have devised an efficient method to incorporate both directional two-point correlation functions S{sub 2} and directional two-point cluster functions C{sub 2} that contain non-trivial topological connectedness information into the simulated annealing microstructure reconstruction procedure. Our framework is applied to model an anisotropic aluminum alloy and the accuracy of the reconstructed structural models is assessed by quantitative comparison with the actual microstructure obtained via x-ray tomography. We show that incorporation of directional clustering information via C{sub 2} significantly improves the accuracy of the reconstruction. In addition, a set of analytical “basis” correlation functions are introduced to approximate the actual S{sub 2} and C{sub 2} of the material. With the proper choice of basis functions, the anisotropic microstructure can be represented by a handful of parameters including the effective linear sizes of the iron-rich and silicon-rich inclusions along three orthogonal directions. This provides a general and efficient means for heterogeneous material modeling that enables one to significantly reduce the data set required to characterize the anisotropic microstructure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A Process-based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Culex pipiens/restuans Mosquito Density

While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration, and validation is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a discrete, semi-stochastic, mechanistic process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens (Diptera, Culicidae) and Culex restuans (Diptera, Culicidae) mosquito populations. This model combines known models for development and survival into a fully connected age-structured model that can reproduce mosquito population dynamics. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like West Nile virus.

59 BASIC BIOLOGICAL SCIENCES↗

The GEWEX LandFlux project: Evaluation of model evaporation using tower-based and globally gridded forcing data

Determining the spatial distribution and temporal development of evaporation at regional and global scales is required to improve our understanding of the coupled water and energy cycles and to better monitor any changes in observed trends and variability of linked hydrological processes. With recent international efforts guiding the development of long-term and globally distributed flux estimates, continued product assessments are required to inform upon the selection of suitable model structures and also to establish the appropriateness of these multi-model simulations for global application. In support of the objectives of the Global Energy and Water Cycle Exchanges (GEWEX) LandFlux project, four commonly used evaporation models are evaluated against data from tower-based eddy-covariance observations, distributed across a range of biomes and climate zones. The selected schemes include the Surface Energy Balance System (SEBS) approach, the Priestley–Taylor Jet Propulsion Laboratory (PT-JPL) model, the Penman–Monteith-based Mu model (PM-Mu) and the Global Land Evaporation Amsterdam Model (GLEAM). Here we seek to examine the fidelity of global evaporation simulations by examining the multi-model response to varying sources of forcing data. To do this, we perform parallel and collocated model simulations using tower-based data together with a global-scale grid-based forcing product. Through quantifying the multi-model response to high-quality tower data, a better understanding of the subsequent model response to the coarse-scale globally gridded data that underlies the LandFlux product can be obtained, while also providing a relative evaluation and assessment of model performance. Using surface flux observations from 45 globally distributed eddy-covariance stations as independent metrics of performance, the tower-based analysis indicated that PT-JPL provided the highest overall statistical performance (0.72; 61 W m –2 ; 0.65), followed closely by GLEAM (0.68; 64 W m –2 ; 0.62), with values in parentheses representing the R 2 , RMSD and Nash–Sutcliffe efficiency (NSE), respectively. PM-Mu (0.51; 78 W m –2 ; 0.45) tended to underestimate fluxes, while SEBS (0.72; 101 W m –2 ; 0.24) overestimated values relative to observations. A focused analysis across specific biome types and climate zones showed considerable variability in the performance of all models, with no single model consistently able to outperform any other. Results also indicated that the global gridded data tended to reduce the performance for all of the studied models when compared to the tower data, likely a response to scale mismatch and issues related to forcing quality. Rather than relying on any single model simulation, the spatial and temporal variability at both the tower- and grid-scale highlighted the potential benefits of developing an ensemble or blended evaporation product for global-scale LandFlux applications. Hence, challenges related to the robust assessment of the LandFlux product are also discussed.

58 GEOSCIENCES↗

Development of the FHR advanced natural circulation analysis code and application to FHR safety analysis

The University of California, Berkeley (UCB) is performing thermal hydraulics safety analysis to develop the technical basis for design and licensing of fluoride-salt-cooled, high-temperature reactors (FHRs). FHR designs investigated by UCB use natural circulation for emergency, passive decay heat removal when normal decay heat removal systems fail. The FHR advanced natural circulation analysis (FANCY) code has been developed for assessment of passive decay heat removal capability and safety analysis of these innovative system designs. The FANCY code uses a one-dimensional, semi-implicit scheme to solve for pressure-linked mass, momentum and energy conservation equations. Graph theory is used to automatically generate a staggered mesh for complicated pipe network systems. Heat structure models have been implemented for three types of boundary conditions (Dirichlet, Neumann and Robin boundary conditions). Heat structures can be composed of several layers of different materials, and are used for simulation of heat structure temperature distribution and heat transfer rate. Control models are used to simulate sequences of events or trips of safety systems. A proportional-integral controller is also used to automatically make thermal hydraulic systems reach desired steady state conditions. A point kinetics model is used to model reactor kinetics behavior with temperature reactivity feedback. The underlying large sparse linear systems in these models are efficiently solved by using direct and iterative solvers provided by the SuperLU code on high performance machines. Input interfaces are designed to increase the flexibility of simulation for complicated thermal hydraulic systems. In conclusion, this paper mainly focuses on the methodology used to develop the FANCY code, and safety analysis of the Mark 1 pebble-bed FHR under development at UCB is performed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Massively Parallel Capability in Sierra/SD for Simulation Vibration with Piezoelectrics

Sierra/SD is an engineering structural dynamics code that provides Sandia and other customers a tool to model structural and acoustic physics on large complex physical systems using massively parallel processing. This report provides a detailed overview on Sierra/SD’s most recent physics package: coupled electro-mechanical physics. This capability uses the finite element method to model coupled electro-mechanical physics exhibited by piezoelectric materials. This report provides an applications overview, theory overview, and verification examples demonstrating the electro-mechanical physics modeling capabilities of Sierra/SD.

97 MATHEMATICS AND COMPUTING↗

Molecular interaction between K-Ras and H-REV107 in the Ras signaling pathway

Ras proteins are small GTPases that serve as master moderators of a large number of signaling pathways involved in various cellular processes. Activating mutations in Ras are found in about one-third of cancers. H-REV107, a K-Ras binding protein, plays an important role in determining K-Ras function. H-REV107 is a member of the HREV107 family of class II tumor suppressor genes and a growth inhibitory Ras target gene that suppresses cellular growth, differentiation, and apoptosis. Expression of H-REV107 was strongly reduced in about 50% of human carcinoma cell lines. However, the specific molecular mechanism by which H-REV107 inhibits Ras is still unknown. In the present study, we suggest that H-REV107 forms a strong complex with activating oncogenic mutation Q61H K-Ras from various biochemical binding assays and modeled structures. In addition, the interaction sites between K-Ras and H-REV107 were predicted based on homology modeling. Here, we found that some structure-based mutants of the K-Ras disrupted the complex formation with H-REV107. Finally, a novel molecular mechanism describing K-Ras and H-REV107 binding is suggested and insights into new K-Ras effector target drugs are provided. - Highlights: • H-REV107 interacts strongly with active GNP-bound K-Ras mutation. • H-REV107 rotated dramatically to interact with K-Ras Q61H (GNP) from the binding position with K-Ras WT (GDP). • In the K-Ras and H-REV107 interaction, E31 of K-Ras is an important binding residue. • H-REV107 can be a stable drug target for oncogenic K-Ras mutation.

60 APPLIED LIFE SCIENCES↗

Divergence in land surface modeling: linking spread to structure

Divergence in land carbon cycle simulation is persistent and widespread. Regardless of model intercomparison project, results from individual models diverge significantly from each other and, in consequence, from reference datasets. Here we link model spread to structure using a 15-member ensemble of land surface models from the Multi-scale synthesis and Terrestrial Model Intercomparison Project (MsTMIP) as a test case. Our analysis uses functional benchmarks and model structure as predicted by model skill in a machine learning framework to isolate discrete aspects of model structure associated with divergence. We also quantify how initial conditions prejudice present-day model outcomes after centennial-scale transient simulations. Overall, the functional benchmark and machine learning exercises emphasize the importance of ecosystem structure in correctly simulating carbon and water cycling, highlight uncertainties in the structure of carbon pools, and advise against hard parametric limits on ecosystem function. We also find that initial conditions explain 90% of variation in global satellite-era values—initial conditions largely predetermine transient endpoints, historical environmental change notwithstanding. As MsTMIP prescribes forcing data and spin-up protocol, the range in initial conditions and high levels of predetermination are also structural. Our results suggest that methodological tools linking divergence to discrete aspects of model structure would complement current community best practices in model development.

58 GEOSCIENCES↗

Fundamental problems in astroparticle physics and cosmology (Final Technical Report)

In the standard cosmological model, structure forms when gas cools and condenses to the center of dark matter halos. The onset of structure is governed by rapid star formation and subsequent evolution of the first massive stars. One of the endpoints of stellar evolution is the formation of a black hole. We showed that all black hole merger events must occur after the universe is about 85 million years old. This prediction is important because it means that if a future experiment makes a detection of a black hole merger event at an earlier time then our understanding of cosmology and the onset of structure formation is not fully developed. In addition, we explored dark matter models that have either self-interactions or experience decays. In the first case, self-interactions may lead to the formation of a thick disk that can be co-planar to the Galactic disc due to correlations in angular momentum between dark matter and baryons that are predicted by tidal torque theory. We showed that for heavy dark matter particles, either the self interaction may not be strong enough to solve the small scale structure motivation or a dark disc cannot be present in the Milky Way. In the second-case, we showed that if dark matter experiences decays into two particles, one of which behaves like radiation, then it is possible to relieve the tension that exists between measurements of the local expansion rate. Finally, we developed methods on how to distinguish gamma-ray excesses from promising dark matter targets as compared to conventional astrophysical sources. This is an important toolset that is applicable to any suspected dark matter source.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

X-ray Spectroscopy Fingerprints of Pristine and Functionalized Graphene

In this work, we demonstrate how to identify and characterize the atomic structure of pristine and functionalized graphene materials from a combination of computational simulation of X-ray spectra, on the one hand, and computer-aided interpretation of experimental spectra, on the other. Despite the enormous scientific and industrial interest, the precise structure of these 2D materials remains under debate. As we show in this study, a wide range of model structures from pristine to heavily oxidized graphene can be studied and understood with the same approach. We move systematically from pristine to highly oxidized and defective computational models, and we compare the simulation results with experimental data. Comparison with experiments is valuable also the other way around; this method allows us to verify that the simulated models are close to the real samples, which in turn makes simulated structures amenable to several computational experiments. Our results provide ab initio semiquantitative information and a new platform for extended insight into the structure and chemical composition of graphene-based materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mol* Viewer: modern web app for 3D visualization and analysis of large biomolecular structures

Abstract Large biomolecular structures are being determined experimentally on a daily basis using established techniques such as crystallography and electron microscopy. In addition, emerging integrative or hybrid methods (I/HM) are producing structural models of huge macromolecular machines and assemblies, sometimes containing 100s of millions of non-hydrogen atoms. The performance requirements for visualization and analysis tools delivering these data are increasing rapidly. Significant progress in developing online, web-native three-dimensional (3D) visualization tools was previously accomplished with the introduction of the LiteMol suite and NGL Viewers. Thereafter, Mol* development was jointly initiated by PDBe and RCSB PDB to combine and build on the strengths of LiteMol (developed by PDBe) and NGL (developed by RCSB PDB). The web-native Mol* Viewer enables 3D visualization and streaming of macromolecular coordinate and experimental data, together with capabilities for displaying structure quality, functional, or biological context annotations. High-performance graphics and data management allows users to simultaneously visualise up to hundreds of (superimposed) protein structures, stream molecular dynamics simulation trajectories, render cell-level models, or display huge I/HM structures. It is the primary 3D structure viewer used by PDBe and RCSB PDB. It can be easily integrated into third-party services. Mol* Viewer is open source and freely available at https://molstar.org/.

Sehnal, David↗

Call for Participation: Collaborative Benchmarking of Functional-Structural Root Architecture Models. The Case of Root Water Uptake

Three-dimensional models of root growth, architecture and function are becoming important tools that aid the design of agricultural management schemes and the selection of beneficial root traits. However, while benchmarking is common in many disciplines that use numerical models, such as natural and engineering sciences, functional-structural root architecture models have never been systematically compared. The following reasons might induce disagreement between the simulation results of different models: different representation of root growth, sink term of root water and solute uptake and representation of the rhizosphere. Presently, the extent of discrepancies is unknown, and a framework for quantitatively comparing functional-structural root architecture models is required. We propose, in a first step, to define benchmarking scenarios that test individual components of complex models: root architecture, water flow in soil and water flow in roots. While the latter two will focus mainly on comparing numerical aspects, the root architectural models have to be compared at a conceptual level as they generally differ in process representation. Therefore, defining common inputs that allow recreating reference root systems in all models will be a key challenge. In a second step, benchmarking scenarios for the coupled problems are defined. We expect that the results of step 1 will enable us to better interpret differences found in step 2. This benchmarking will result in a better understanding of the different models and contribute toward improving them. Improved models will allow us to simulate various scenarios with greater confidence and avoid bugs, numerical errors or conceptual misunderstandings. This work will set a standard for future model development.

Plant Sciences↗

Cryo-based Structural Characterization and Growth Model of Salt Film on Metal

Metal salt films play a critical role in high rate dissolution processes such as within an active pit. However, structural and compositional characterization of salt films is challenging owing to the dynamic nature of the dissolving interface. Here, the salt film formed on SS304 one-dimensional artificial pit surface was directly characterized using flash freezing and cryo-based focused ion beam/scanning electron microscopy. The salt film exhibits a porous structure. Ni is co-located with Fe in the salt film, indicating the formation of co-precipitated salt containing Fe2+ and Ni2+. Additionally, a salt film growth model is presented based on ion transport.

Salt film, Pitting corrosion, Cryogenic, Cryo-FIB/↗

Temperature induced phase transition of CaMn{sub 0.5}Zr{sub 1.5}(PO{sub 4}){sub 3} phosphate

In this work we investigated the structural behaviour of a CaMn{sub 0.5}Zr{sub 1.5}(PO{sub 4}){sub 3}. Due to the presence of divalent Mn{sup 2+} cations this compound can possess interesting luminescence properties. It was recently understood that this phosphate undergoes a temperature induced irreversible phase transition in the range of 800–875 °C. It has also been shown that the 3d–3d luminescence of Mn{sup 2+} increases 10 fold for the high temperature polymorph. To determine the Mn environment structural investigations of both phases have been performed by the X-ray powder diffraction and Raman spectroscopy methods. The low temperature modification adopts the trigonal NZP structure type with a slightly lower symmetry (space group R32, a=8.7850(2) Å, c=22.6496(7) Å, V=1514.8(1) Å{sup 3}). The high temperature form in turn has orthorhombic symmetry (space group Pnma, a=6.2350(3) Å, b=6.6281(3) Å, c=14.4731(6) Å, V=598.13(5) Å{sup 3}). Both structures were solved ab-initio from powder data and structural analysis was performed. In-situ and RT Raman spectra are consistent with the XRD derived structural model. Mn{sup 2+} cations occupy different types of positions in these structures and a change in Mn coordination number (6 for LT phase, 7 for HT phase) results in different Mn–O bond lengths. These differences may explain the change in the optical properties between the polymorphs. - Graphical abstract: The compound CaMn{sub 0.5}Zr{sub 1.5}(PO{sub 4}){sub 3} was synthesized in order to create a material with enhanced luminescent properties. The goal of present studies is to define Mn{sup 2+} environment and its changes due to the structural transformations of the phosphate along phase transition at the T range of 800–875 °C. It was found that LT modification adopts the trigonal NZP structure type, sp.gr. R32, the HT form in turn exhibits orthorhombic symmetry sp.gr. Pnma. Mn2+ cations occupy different types of positions in those structures and a change in coordination number of Mn (6 for LT phase, 7 for HT phase) results in a change in Mn–O bond lengths.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Pseudo-Polymorphism in Layered FeS Intercalates: A Competition between Charged and Neutral Guest Species

Systematic synthesis studies of the formation of tetrahedral FeS-ethylenediamine intercalates resulted in the synthesis of a new compound, [Fe 9.4(2) S 10 ][Fe(en) 3 ] 0.6(1) ·en 0.9 ( 3 ). The composition and complex crystal structure were determined based on a synergistic combination of elemental composition, decomposition behavior, high-resolution synchrotron X-ray diffraction and total scattering, 57 Fe Mössbauer spectroscopy, and electron diffraction. The structural model was derived based on a systematic comparison to the previously reported structures [Fe 8 S 10 ][Fe(en) 3 ] 1 ·en 0.5 and tetragonal FeS. The new compound has flat Fe 9.4 S 10 layers, analogous to those in superconducting binary FeS. In the crystal structure of [Fe 9.4 S 10 ][Fe(en) 3 ] 0.6 ·en 0.9 , the interlayer space is occupied by [Fe(en) 3 ] 2+ complexes and neutral ethylenediamine molecules in a ~2:3 ratio. Interlayer species are not randomly oriented but ordered as evidenced by superstructural diffraction peaks in both high-resolution X-ray diffraction and electron diffraction patterns. Magnetic studies reveal no superconducting transition down to 2 K, indicating that the presence of minute amounts (~6%) of iron vacancies at the Fe-S layer in [Fe 9.4 S 10 ][Fe(en) 3 ] 0.6 ·en 0.9 is still sufficient to shift the position of the Fermi level resulting in an adjustment of the properties. Here, our work shows the importance of detailed characterization of the crystal structure of intercalated compounds to understand the origin of the observed properties and develop proper structure–property relationships.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fusion of Multiple Models for Improving Gross Primary Production Estimation With Eddy Covariance Data Based on Machine Learning

Abstract Terrestrial gross primary production (GPP) represents the magnitude of CO 2 uptake through vegetation photosynthesis, and is a key variable for carbon cycles between the biosphere and atmosphere. Light use efficiency (LUE) models have been widely used to estimate GPP for its physiological mechanisms and availability of data acquisition and implementation, yet each individual GPP model has exhibited large uncertainties due to input errors and model structure, and further studies of systematic validation, comparison, and fusion of those models with eddy covariance (EC) site data across diverse ecosystem types are still needed in order to further improve GPP estimation. We here compared and fused five GPP models (VPM, EC‐LUE, GOL‐PEM, CHJ, and C‐Fix) across eight ecosystems based on FLUXNET2015 data set using the ensemble methods of Bayesian Model Averaging (BMA), Support Vector Machine (SVM), and Random Forest (RF) separately. Our results showed that for individual models, EC‐LUE gave a better performance to capture interannual variability of GPP than other models, followed by VPM and GLO‐PEM, while CHJ and C‐Fix were more limited in their estimation performance. We found RF and SVM were superior to BMA on merging individual models at various plant functional types (PFTs) and at the scale of individual sites. On the basis of individual models, the fusion methods of BMA, SVM, and RF were examined by a five‐fold cross validation for each ecosystem type, and each method successfully improved the average accuracy of estimation by 8%, 18%, and 19%, respectively.

Environmental Sciences & Ecology↗