SEARCH · Engineering Papers

Results for “spatial pattern”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Attend and Decode: 4D fMRI Task State Decoding Using Attention Models

Source code for Brain Attend and Decode paper. Functional magnetic resonance imaging (fMRI) is a neuroimaging modality that captures the blood oxygen level in a subject's brain while the subject either rests or performs a variety of functional tasks under different conditions. Given fMRI data, the problem of inferring the task, known as task state decoding, is challenging due to the high dimensionality (hundreds of million sampling points per datum) and complex spatio-temporal blood flow patterns inherent in the data. In this work, we propose to tackle the fMRI task state decoding problem by casting it as a 4D spatiotemporal classification problem. We present a novel architecture called Brain Attend and Decode (BAnD), that uses residual convolutional neural networks for spatial feature extraction and self-attention mechanisms for temporal modeling. We achieve significant performance gain compared to previous works on a 7-task benchmark from the large-scale Human Connectome Project-Young Adult (HCP-YA) dataset. We also investigate the transferability of BAnD's extracted features on unseen HCP tasks, either by freezing the spatial feature extraction layers and retraining the temporal model, or finetuning the entire model. The pre-trained features from BAnD are useful on similar tasks while finetuning them yields competitive results on unseen tasks/conditions.

Ng, BrendaM.↗

DOE OSTI2021

Evolution of microstructures in radiation fields using a coupled binary-collision Monte Carlo phase field approach

The simulation of radiation effects in materials broadly falls into two categories. At short time and length scales lies the modeling of primary radiation damage, such as point defect creation, energy deposition, and ballistic mixing. This is followed by the modeling at longer time scales of thermally activated microstructure evolution and defect reactions, such as recombination, clustering, and coarsening. The binary collision Monte Carlo method is an established, numerically efficient method for the computation of primary radiation damage. Conversely, the phase field method is a state of the art method for microstructure evolution on longer time and length scales. Here we present a concurrent coupling of these two methods, overcoming the difference between the discrete object Monte Carlo paradigm for primary radiation damage and the continuum field variable approach for microstructure evolution. The coupling is bidirectional, in which the microstructure evolution in the MOOSE finite element frame- work provides the spatial scattering data set for the charged particle trans- port and receives point defect, mass transport, and heat source terms from the simulated collision cascades that contribute to the field variable evolution. The concurrent coupling scheme is implemented in the code Magpie and demonstrated by investigating patterning for a model irradiated immiscible binary alloy. The results from the coupled binary collision Monte Carlo/phase field simulations reproduce the results of analytical models for phase separation, phase mixing, and patterning, supporting the approach and indicating its utility for modeling real materials systems.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Effects of exogenous citric acid on the concentration and spatial distribution of Ni, Zn, Co, Cr, Mn and Fe in leaves of Noccaea caerulescens grown on a serpentine soil

The aim of this study was to show the potential of citric acid in increasing the concentration of Ni, Zn, Co, Cr, Mn and Fe in leaves of the hyperaccumulator Noccaea caerulescens. Synchrotron x-ray fluorescence (μ-XRF) images were collected to assess the distribution of metals in leaves. Applying citric acid (20 mmol kg -1 ) to soil increased in 14-, 10-, 7-, 2- and 1.4- fold the concentration of Mn, Fe, Co, Ni, and Cr, respectively, compared to the control. The μ-XRF imaging revealed that Ni and Zn were not spatially correlated across the leaf. We observed a clear partitioning of Zn between veins and surrounding leaf cells while Ni was more evenly distributed between veins and leaf blade. The accumulation of metals in citric acid treated plants did not change the Ni and Zn distribution pattern in leaves but altered the Mn distribution. It seems that Mn reached toxic concentrations in leaves and we hypothesize that a mechanism driven by transpiration through the xylem was used to excrete the metal. Our results show that citric acid can enhance metal accumulation by N. caerulescens and have impact for soil remediation by either decreasing the time for clean up or increasing the access to non-labile pools of metals in soil.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Sub-10 nm Probing of Ferroelectricity in Heterogeneous Materials by Machine Learning Enabled Contact Kelvin Probe Force Microscopy

Reducing the dimensions of ferroelectric materials down to the nanoscale has strong implications on the ferroelectric polarization pattern and on the ability to switch the polarization. As the size of ferroelectric domains shrinks to the nanometer scale, the heterogeneity of the polarization pattern becomes increasingly pronounced, enabling a large variety of possible polar textures in nanocrystalline and nanocomposite materials. Critical to the understanding of fundamental physics of such materials and hence their applications in electronic nanodevices is the ability to investigate their ferroelectric polarization at the nanoscale in a nondestructive way. We show that contact Kelvin probe force microscopy (cKPFM) combined with a k-means response clustering algorithm enables to measure the ferroelectric response at a mapping resolution of 8 nm. In a BaTiO 3 thin film on silicon composed of tetragonal and hexagonal nanocrystals, we determine a nanoscale lateral distribution of discrete ferroelectric response clusters, fully consistent with the nanostructure determined by transmission electron microscopy. Moreover, we apply this data clustering method to the cKPFM responses measured at different temperatures, which allows us to follow the corresponding change in the polarization pattern as the Curie temperature is approached and across the phase transition. This work opens up perspectives for mapping complex ferroelectric polarization textures such as curled/swirled polar textures that can be stabilized in epitaxial heterostructures and more generally for mapping the polar domain distribution of any spatially highly heterogeneous ferroelectric materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Autonomous Multistate Nanoencoding Using Combinatorial Ferroelectric Closure Domains in BiFeO 3

Recent advances in ferroic materials have identified topological defects as promising candidates for enabling additional functionalities in future electronic systems. The generation of stable and customizable polar topologies is needed to achieve multistates that enable beyond-binary device architectures. Here, in this study, we show how to autonomously pattern on-demand highly tunable striped closure domains in pristine rhombohedral-phase BiFeO 3 thin films through precise scanning of a biased atomic force microscopy tip along carefully designed paths. By employing this strategy, we generate and manipulate closed-loop structures with high spatial resolution in an automated manner, allowing the creation of highly tunable and intricate topological domain structures that exhibit distinct polarization configurations without the need for electrode deposition or complex heterostructure growth. As a proof-of-concept for ferroelectric beyond-binary memory devices, we use such topological domains as multistates, engineering an alphabet and automating the symbolic writing/reading process using autonomous microscopy. The resulting information density is compared with that of current commercially available memory devices, demonstrating the potential of ferroelectric topological domains for multistate information storage applications.

BiFeO3↗

DOE OSTI2021

Cellular fluidics

The natural world provides many examples of multiphase transport and reaction processes that have been optimized by evolution. These phenomena take place at multiple length and time scales and typically include gas–liquid–solid interfaces and capillary phenomena in porous media. Many biological and living systems have evolved to optimize fluidic transport. However, living things are exceptionally complex and very difficult to replicate, and human-made microfluidic devices (which are typically planar and enclosed) are highly limited for multiphase process engineering. In this paper, we introduce the concept of cellular fluidics: a platform of unit-cell-based, three-dimensional structures—enabled by emerging 3D printing methods—for the deterministic control of multiphase flow, transport and reaction processes. We show that flow in these structures can be ‘programmed’ through architected design of cell type, size and relative density. We demonstrate gas–liquid transport processes such as transpiration and absorption, using evaporative cooling and CO 2 capture as examples. We design and demonstrate preferential liquid and gas transport pathways in three-dimensional cellular fluidic devices with capillary-driven and actively pumped liquid flow, and present examples of selective metallization of pre-programmed patterns. Our results show that the design and fabrication of architected cellular materials, coupled with analytical and numerical predictions of steady-state and dynamic behaviour of multiphase interfaces, provide deterministic control of fluidic transport in three dimensions. Cellular fluidics may transform the design space for spatial and temporal control of multiphase transport and reaction processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Simplex‐based model for nanoparticle grain identification in four‐dimensional scanning transmission electron microscopy data

Grain identification in polycrystalline nanoparticles, for example, determining which crystal phases are present at each spatial location, is fundamental to materials characterisation. This is particularly challenging when grains overlap extensively, as commonly occurs in four-dimensional scanning transmission electron microscopy (4D-STEM) datasets. We propose a simplex-based model (SBM) in which each simplex vertex represents the diffraction pattern (DP) of a pure grain, and the simplex edges and interior represent overlapping grains. Our SBM grain identification algorithm operates on the Bragg disk (BD) data matrix distilled from the 4D-STEM data to identify the grain membership at each scan position, together with a BD feature matrix whose columns represent the DPs for each constituent grain, which is important for identifying the crystal structure of each grain. We solve the model using a two-stage algorithm. In Stage 1, we adapt a linear mixing algorithm to estimate an initial BD feature matrix whose columns represent DPs of potentially overlapping grains. Our Stage 2 algorithm incorporates sparsity considerations to transform the initial BD feature matrix so that its columns represent DPs of pure grains. Using simulated datasets with various grain configurations, we demonstrate that SBM recovers both the BD feature matrix and membership maps more accurately than existing methods, even when a grain lacks any pure region and completely overlaps with other grains.

4D-STEM segmentation↗

DOE OSTI2021

GOTHIC Aerosol Source Depletion Studies

Pacific Northwest National Laboratory has continued work to develop an aerosol-laden flow modeling capability with the Generation of Thermal Hydraulic Information in Containment (GOTHIC TM ) computer code to perform simulations for thermal hydraulic conditions and aerosol transport and deposition in spent fuel casks. This report describes our recent work to expand our model that was originally developed in 2019 to allow for thermal characterization, carrier gas flow determination, and tracking of particulate behavior throughout the entire canister volume. This was achieved through conversion of decay heat source in all fuel tubes within the model from heaters to thermal conductors in GOTHIC as well as remeshing the internal volume of the canister. These model improvements allowed for simulations of base case scenarios for comparison of similar efforts by our collaborators with alternative tools at Sandia and Oak Ridge. Several characteristics of the GOTHIC code were elucidated by this effort and were documented for consideration in how the code is used for this effort moving forward. For instance, treatment of particle size distributions, initial spatial distributions, mesh gradients, and minimum volume concentrations were identified as important issues for consideration. Ultimately, the current version of the code is capable of tracking temperatures, flow rates, and particle behavior throughout the canister internal volume. This report presents predicted depletion times and preferred deposition patterns within the canister. Additionally, this report documents code performance and insights for further model development and use for comparisons to other codes and for predicting system behavior in experiments being performed and planned by our collaborators.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Longitudinal, Lateral, Vertical, and Temporal Thermal Heterogeneity in a Large Impounded River: Implications for Cold-Water Refuges

Dam operations can affect mixing of the water column, thereby influencing thermal heterogeneity spatially and temporally. This occurs by restricting or eliminating connectivity in longitudinal, lateral, vertical, and temporal dimensions. We examined thermal heterogeneity across space and time and identified potential cold-water refuges for salmonids in a large impounded river in inland northwestern USA. To describe these patterns, we used thermal infrared (TIR) imagery, in situ thermographs, and high-resolution, 3-D hydraulic mapping. We explained the median water temperature and probability of occurrence of cool-water areas using generalized additive models (GAMs) at reach and subcatchment scales, and we evaluated potential cold-water refuge occurrence in relation to these patterns. We demonstrated that (1) lateral contributions from tributaries dominated thermal heterogeneity, (2) thermal variability at confluences was approximately an order of magnitude greater than of the main stem, (3) potential cold-water refuges were mostly found at confluences, and (4) the probability of occurrence of cool areas and median water temperature were associated with channel geomorphology and distance from dam. These findings highlight the importance of using multiple approaches to describe thermal heterogeneity in large, impounded rivers and the need to incorporate these types of rivers in the understanding of thermal riverscapes because of their limited representation in the literature.

Mejia, Francine H.↗

DOE OSTI2020

Intra-Annual and Interannual Dynamics of Evaporation Over Western Lake Erie

Evaporation (E) is a critical component of the water and energy budget in lake systems yet is challenging to quantify directly and continuously. We examined the magnitude and changes of E and its drivers over Lake Erie - the shallowest and most southern lake of the Laurentian Great Lakes. We deployed two eddy-covariance tower sites in the western Lake Erie Basin - one located nearshore (CB) and one offshore (LI) - from September 2011 through May 2016. Monthly E varied from 5 to 120 mm, with maximum E occurring in August–October. The annual E was 635 ± 42 (±SD) mm at CB and 604 ± 32 mm at LI. Mean winter (October–March) E was 189 ± 61 mm at CB and 178 ± 25 mm at LI, accounting for 29.8% and 29.4% of annual E. Mean daily E was 1.8 mm during the coldest month (January) and 7.4 mm in the warmest month (July). Monthly E exhibited a strong positive linear relationship to the product of wind speed and vapor pressure deficit. Pronounced seasonal patterns in surface energy fluxes were observed with a 2-month lag in E from R n , due to the lake's heat storage. This lag was shorter than reports regarding other Great Lakes. Difference in E between the offshore and nearshore sites reflected within-lake spatial heterogeneity, likely attributable to climatic and bathymetric differences between them. These findings suggest that predictive models need to consider lake-specific heat storage and spatial heterogeneity in order to accurately simulate lake E and its seasonal dynamics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Spatial ecology of translocated raccoons

Raccoons (Procyon lotor) are routinely translocated both legally and illegally to mitigate conflicts with humans, which has contributed to the spread of rabies virus across eastern North America. The movement behavior of translocated raccoons has important ramifications for disease transmission yet remains understudied and poorly quantified. To examine the spatial ecology of raccoons following experimental translocation, we performed reciprocal 16 km-distance translocations of 30 raccoons between habitats of high and low raccoon density (bottomland hardwood and upland pine, respectively) across the Savannah River Site (SRS) in Aiken, South Carolina, USA (2018–2019). Translocation influenced patterns of raccoon space use, with translocated animals exhibiting a 13-fold increase in 95% utilization distributions (UDs) post- compared to pre-translocation (mean 95% UD 35.8 ± 36.1 km 2 vs 1.96 ± 1.17 km 2 ). Raccoons originating from upland pine habitats consistently had greater space use and larger nightly movement distances post-translocation compared to raccoons moved from bottomland hardwood habitats, whereas these differences were generally not observed prior to translocation. Estimated home ranges of male raccoons were twice the area as estimated for female raccoons, on average, and this pattern was not affected by translocation. After a transient period lasting on average 36.5 days (SD = 30.0, range = 3.25–92.8), raccoons often resumed pre-experiment movement behavior, with 95% UD sizes not different from those prior to translocation (mean = 2.27 ± 1.63km 2 ). Most animals established new home ranges after translocation, whereas three raccoons moved > 16 km from their release point back to the original capture location. Four animals crossed a 100-m wide river within the SRS post-translocation, but this behavior was not documented among collared raccoons prior to translocation. Large increases in space use combined with the crossing of geographic barriers such as rivers may lead to elevated contact rates with conspecifics, which can heighten disease transmission risks following translocation. These results provide additional insights regarding the potential impacts of raccoon translocation towards population level risks of rabies outbreaks and underscore the need to discourage mesocarnivore translocations to prevent further spread of wildlife rabies.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Relative reproductive phenology and synchrony affect neonate survival in a nonprecocial ungulate

Abstract Degree of reproductive synchronization in prey is hypothesized as a predator defense strategy reducing prey risk via predator satiation or predator avoidance. Species with precocial young, especially those exposed to specialist predators, should be highly synchronous to satiate predators (predator satiation hypothesis), while prey with nonprecocial (i.e. altricial) young, especially those exposed to generalist predators, should become relatively asynchronous to avoid predator detection (predator avoidance hypothesis). The white‐tailed deer Odocoileus virginianus in North America is an example of a nonprecocial ungulate that uses the hider strategy early in life; its primary predator (coyote; Canis latrans ) is a generalist, making white‐tailed deer a good model species to test the predator avoidance hypothesis. We used birth dates and known fates of white‐tailed deer neonates ( n = 1,032) across nine study sites varying in relative synchrony and predator assemblages to test the predator avoidance hypothesis. We predicted that relative birthing asynchrony of the population would increase relative survival at the population level; therefore, at the individual scale, neonate birth date nearer to mean birthing date in a respective population would not influence individual survival. Coyotes were responsible for the majority of predation events, and survival of those neonates increased the closer the individual was born to peak birthing season in each respective population. Also, at the population level, reproductive asynchronization negatively affected survival. Contrary to the predator avoidance hypothesis, our data indicate patterns in neonate survival for white‐tailed deer better support the predator satiation hypothesis at the individual and population level. Additionally, coyotes may present a selective force great enough to shift reproductive synchrony such that predator satiation may become a feasible defense strategy for neonates at local spatial scales. Our results indicate that synchronizing reproduction may still be the most effective strategy to reduce individual predation risk from generalist predators, particularly when the window of heightened resource availability to the prey is narrow. A free Plain Language Summary can be found within the Supporting Information of this article.

Michel, Eric S.↗

DOE OSTI2024

Polarization phenomenon in heavy-ion collisions

The strongly interacting system created in ultrarelativistic nuclear collisions behaves almost as an ideal fluid with rich patterns of the velocity field exhibiting strong vortical structure. Vorticity of the fluid, via spin-orbit coupling, leads to particle spin polarization. Due to the finite orbital momentum of the system, the polarization on average is not zero. It depends on the particle momenta reflecting the spatial variation of the local vorticity. In the last few years, this field experienced a rapid growth due to experimental discoveries of the global and local polarizations. Recent measurements triggered further development of the theoretical description of the spin dynamics and suggestions of several new mechanisms for particle polarization. In this review, we focus mostly on the experimental results. We compare the measurements with the existing theoretical calculations but try to keep the discussion of possible underlying physics at the qualitative level. Future measurements and how they can help to answer open theoretical questions are also discussed. We pay a special attention to the employed experimental methods, as well as to the detector effects and associated corrections to the measurements.

Physics↗

DOE OSTI2025

OpenCRUMS USA: An Open Machine Learning Framework for Characterizing Variability in Aerosol Reanalysis Data

Advances in artificial intelligence (AI) have called for exploring how these techniques can be used for exploring patterns in large climate datasets. To that regard, the U.S. Department of Energy AI for Earth System Predictability (AI4ESP) supported a pilot initiative called the Open Classification of Regimes in the Southeast USA (OpenCRUMS USA) project to explore how AI can be used to characterize modes of spatial variability in large climate datasets. For this study, we focus on comparing two methods for characterizing the modes of spatial variability of surface aerosol concentration over the Houston region: empirical orthogonal functions (EOFs) and layerwise relevance propagation (LRP) applied to a convolutional neural network (CNN) classifier. We show that EOF analysis typically attributes spatial variability modes that span all of southeast Texas, prohibiting the attribution of spatial variability to localized regions. However, using LRP on the CNN classifier resolves the explanatory parameters at a finer spatial resolution than EOFs. This allows for the attribution of the spatial variability of surface aerosols to local regions of organic carbon which was not possible using EOFs. In addition, the LRP analysis also suggests that synoptic-scale transport of dust is most prevalent during anticyclonic and pretrough synoptic conditions as categorized by self-organizing maps.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Fermilab Booster loss modelling and rebalancing using Bayesian methods

To meet PIP-II upgrade requirements, Fermilab Booster losses need to be reduced by 50% compared to present levels. So far, simulations are not good enough to predict loss patterns. Thus, an extensive Booster tune up will be necessary to achieve required performance. In this paper we present an effort to build a data-driven loss model using Bayesian techniques, and subsequently to rebalance losses for higher trip margins. We first created several sets of spatially and temporally isolated orbit and optics knobs, and trained Gaussian process models for each beam loss monitor as well as beam current. Novel techniques of uncertainty constraints and approximate GP fitting were introduced to handle safety and timing requirements. We then performed single and multi-objective tuning using scalarized objectives comprised of critical beam loss locations. We achieved significant rebalancing of losses, increasing margins by 25%, as well as an overall improvement in transmission efficiency of 0.4%. Automated data collection is being developed so that more accurate surrogate models can be trained over time.

Kuklev, N. [Fermilab]↗

DOE OSTI2020

Drought Drives Spatial Variation in the Millet Root Microbiome

Efforts to boost crop yield and meet global food demands while striving to reach sustainability goals are hindered by the increasingly severe impacts of abiotic stress, such as drought. One strategy for alleviating drought stress in crops is to utilize rootassociated bacteria, yet knowledge concerning the relationship between plant hosts and their microbiomes during drought remain under-studied. One broad pattern that has recently been reported in a variety of monocot and dicot species from both native and agricultural environments, is the enrichment of Actinobacteria within the drought-stressed root microbiome. In order to better understand the causes of this phenomenon, we performed a series of experiments in millet plants to explore the roles of drought severity, drought localization, and root development in provoking Actinobacteria enrichment within the root endosphere. Through 16S rRNA ampliconbased sequencing, we demonstrate that the degree of drought is correlated with levels of Actinobacterial enrichment in four species of millet. Additionally, we demonstrate that the observed drought-induced enrichment of Actinobacteria occurs along the length of the root, but the response is localized to portions of the root experiencing drought. Finally, we demonstrate that Actinobacteria are depleted in the dead root tissue of Japanese millet, suggesting saprophytic activity is not the main cause of observed shifts in drought-treated root microbiome structure. Collectively, these results help narrow the list of potential causes of drought-induced Actinobacterial enrichment in plant roots by showing that enrichment is dependent upon localized drought responses but not root developmental stage or root death.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

In Situ Transmission Electron Microscopy of High-Temperature Inconel-625 Corrosion by Molten Chloride Salts

This paper describes an approach to monitor high temperature molten chloride (MgCl 2 -NaCl-KCl) salt corrosion of Inconel-625 alloy in real time at high spatial resolution. The approach is based on a micro-environmental-cell assembly integrated into a transmission-electron-microscope goniometer to examine in situ the salt-alloy interface during corrosion, employing real time electron diffraction and imaging. It establishes procedures to minimize incorporation of H 2 O or O 2 from atmosphere in the chloride salts during sample fabrication and corrosion, which is critical to understanding the fundamental corrosion mechanisms. A clustering algorithm and a 2D Gaussian fit function are used to determine diffraction spot intensities in in situ diffraction patterns, to quantify alloy corrosion. This facilitates quantitative observation of the evolution of individual grains, in contrast to conventional macroscopic corrosion rate quantification. The isothermal corrosion rate of Inconel-625 in an anhydrous, unoxidized salt-stack is 220 ± 30 μm year -1 at 700 °C and 350 ± 20 μm year -1 at 800 °C. However, the corrosion rate at 700 °C increases five-fold to 1000 ± 170 μm year -1 when the salt stack is air-exposed, indicating the dominant effects of hydrated or oxidized impurities on corrosion acceleration. Furthermore, real time imaging of the microstructure evolution suggests that corrosion is initiated at grain boundaries.

14 SOLAR ENERGY↗

DOE OSTI2021

Long-term leaf C:N ratio change under elevated CO 2 and nitrogen deposition in China: Evidence from observations and process-based modeling

Climate change, elevating atmosphere CO 2 (eCO 2 ) and increased nitrogen deposition (iNDEP) are altering the biogeochemical interactions between plants, microbes and soils, which further modify plant leaf carbon-nitrogen (C:N) stoichiometry and their carbon assimilation capability. Many field experiments have observed large sensitivity of leaf C:N ratio to eCO 2 and iNDEP. However, the large-scale pattern of this sensitivity is still unclear, because eCO 2 and iNDEP drive leaf C:N ratio toward opposite directions, which are further compounded by the complex processes of nitrogen acquisition and plant-and-microbial nitrogen competition. Here, we attempt to map the leaf C:N ratio spatial variation in the past 5 decades in China with a combination of data-driven model and process-based modeling. These two approaches showed consistent results. Over different regions, we found that leaf C:N ratio had significant but uneven changes between 2 time periods (1960-1989 and 1990-2015): a 5% ± 8% increase for temperate grasslands in northern China, a 3% ± 6% increase for boreal grasslands in western China, and by contrast, a 7% ± 6% decrease for temperate forests in southern China, and a 3% ± 5% decrease for boreal forests in northeastern China. Additionally, the structural equation models indicated that the leaf C:N change was sensitive to ΔNDEP, ΔCO 2 and ΔMAT rather than ΔMAP and ecosystem types. In this work, process-based modeling suggested that iNDEP was the main source of soil mineral nitrogen change, dominating leaf C:N ratio change in most areas in China, while eCO 2 led to leaf C:N ratio increase in low iNDEP area. This study also indicates that the long-term leaf C:N ratio acclimation was dominated by climate constraint, especially temperature, but was constrained by soil N availability over decade scale.

54 ENVIRONMENTAL SCIENCES↗