Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Implicit”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Polarization and properties of low-frequency waves in warm magnetized two-fluid plasma

This paper presents the derivation of a general wave dispersion relation for warm magnetized plasma under the two-fluid formalism. The discussion is quite general except for the assumption of low frequency and slow phase speed, for which the displacement current is negligible, under the implicit assumption that the plasma is sufficiently dense to satisfy the condition ωpe>ωce, where ωpe and ωce denote the plasma oscillation frequency and electron gyro frequency, respectively. The present discussion does not invoke charge neutrality associated with the fluctuations although it is implicitly satisfied. The resulting dispersion relation that includes the fluid thermal effects shows that there are three eigen modes, which include those corresponding to ideal MHD, namely, fast, slow, and kinetic Alfvén waves, as well as higher-frequency modes including the ion and electron cyclotron and lower-hybrid resonances. The fluid effects in the ideal MHD wave branches are influenced by the finite Larmor radius scales, and when the wave number in the cross field direction is comparable to these values, the fluid effects become significant. It is found that the Larmor radius should be interpreted in the sense as ion-acoustic gyro-radius instead of ion thermal gyro radius only. That is, it is found that the electrons also contribute to the non-ideal effect associated with the kinetic Alfvén wave. A comprehensive explanation of the polarization of each mode is also presented. The present findings indicate that the polarity may change its sign only for the kinetic Alfvén mode branch and that such a transition is based on the propagation angle. When such a change does take place, it is found that the kinetic Alfvén wave transits to an ion-acoustic mode. For each branch, it is also found that the electric field along the ambient magnetic field is purely transverse.

Physics↗
DOE OSTI2022

Driving force and pathway in polyelectrolyte complex coacervation

There is notable discrepancy between experiments and coarse-grained model studies regarding the thermodynamic driving force in polyelectrolyte complex coacervation: experiments find the free energy change to be dominated by entropy, while simulations using coarse-grained models with implicit solvent usually report a large, even dominant energetic contribution in systems with weak to intermediate electrostatic strength. Here, using coarse-grained, implicit-solvent molecular dynamics simulation combined with thermodynamic analysis, we study the potential of mean force (PMF) in the two key stages on the coacervation pathway for symmetric polyelectrolyte mixtures: polycation–polyanion complexation and polyion pair–pair condensation. We show that the temperature dependence in the dielectric constant of water gives rise to a substantial entropic contribution in the electrostatic interaction. By accounting for this electrostatic entropy, which is due to solvent reorganization, we find that under common conditions (monovalent ions, room temperature) for aqueous systems, both stages are strongly entropy-driven with negligible or even unfavorable energetic contributions, consistent with experimental results. Furthermore, for weak to intermediate electrostatic strengths, this electrostatic entropy, rather than the counterion-release entropy, is the primary entropy contribution. From the calculated PMF, we find that the supernatant phase consists predominantly of polyion pairs with vanishingly small concentration of bare polyelectrolytes, and we provide an estimate of the spinodal of the supernatant phase. Finally, we show that prior to contact, two neutral polyion pairs weakly attract each other by mutually induced polarization, providing the initial driving force for the fusion of the pairs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Validation of a Hybrid Domain Overlapping Coupling Between SAM and CFD Against the TALL-3D Transients

The System Thermal Hydraulics (STH) code SAM has been coupled to the Computational Fluid Dynamics (CFD) code Simcenter STAR-CCM+ utilizing a hybrid domain overlapping method with an explicit coupling in time. The coupling aims to extend the STH code’s applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional (3D) mixing. The coupling method’s numerical stability was verified in the past against two closed-loop configurations, and it was validated against a double T-junction experiment with 3D scalar mixing. In the present work, the coupling method is validated against the TALL-3D STH/CFD coupling benchmark facility. TALL-3D is a three-legged, liquid-metal facility with a large, pool-type enclosure (test section) that exhibits 3D flow effects to be modeled by a CFD code. The rest of the system exhibits approximately 1D behavior well-predicted by an STH code. First, the present STAR-CCM+ CFD model of the 3D test section is validated against experimental data. Then, the SAM-STARCCM+ coupled model is validated against six different TALL-3D steady states, including SAM standalone model results for comparison. Lastly, the SAM-STARCCM+ coupled model is validated against two TALL-3D transients, one exhibiting flow reversal in the test section and one exhibiting nonlinear, Limit Cycle Oscillations (LCO). For the first transient, the SAM-STARCCM+ coupled model properly predicts an increase in the test section’s inlet temperature during flow reversal, and this is not predicted by the SAM standalone model. Following flow reversal, the SAM-STARCCM+ coupled model better-predicts the initial flow recovery and following oscillations as the system approaches a final natural circulation state. For the second transient, no true final steady state is observed due to LCO. Neither the SAM-STARCCM+ coupled model nor the SAM standalone model can perfectly capture the experiment’s changing oscillation frequency during the transient. However, the SAM-STARCCM+ coupled model does reproduce the oscillatory feedback observed in the system. This is a significant achievement as the SAM-STARCCM+ coupled model only uses an explicit coupling in time, as opposed to a semi-implicit coupling. In comparison, previous STH/CFD coupling efforts of the TALL-3D facility required semi-implicit coupling to obtain similar results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Potential improvements in brain dose estimates for internal emitters

Background: Element-specific biokinetic models are used to reconstruct doses to systemic tissues from internal emitters. Typically, a systemic model for a radionuclide explicitly depicts only its dominant repositories. Remaining tissues and fluids are aggregated into a pool called Other tissue in which the radionuclide is assumed to be uniformly distributed. In the systemic biokinetic models used in radiation protection, the brain usually is addressed as an implicit mass fraction of Other tissue rather than an explicitly depicted repository. Additionally, due to increasing interest in radiation effects on the brain, efforts are underway to improve brain dosimetry for internal radiation sources. Methods: We assessed potential improvements in brain dosimetry for internal emitters by explicitly modelling brain kinetics rather than treating the brain as a mass fraction of Other tissue. We selected ten elements for which brain kinetics can be modeled using published biokinetic data. Injection dose coefficients were calculated for a relatively long-lived radioisotope of each element using each of two versions of the ICRP’s latest systemic biokinetic model for each element, the original version and a modified version differing only in the treatment of brain. If the ICRP model contained an explicit brain pool, the modified version depicted brain instead as a mass fraction of Other tissue. If the ICRP model included brain in Other tissue, the modified version included an explicit brain pool with kinetics based on best available brain-specific data.Results: The result for a given radionuclide is expressed as a ratio A:B, where A and B are the dose coefficients based on the versions of the model with and without an explicit brain pool, respectively. The following ratios A:B were obtained for the 10 radionuclides addressed here: 241 Am, 0.13; 207 Bi, 0.57; 234 U, 0.81; 239 Pu, 0.96; 203 Hg (vapor), 1.4; 134 Cs, 1.5; 54 Mn, 1.7; 210 Po, 1.7; 226 Ra, 1.9; 210 Pb, 3.3. These ratios indicate that a dose estimate for brain based on a biokinetic model with brain implicitly contained in Other tissue may substantially underestimate or substantially overestimate a dose estimate that reflects best available brain-specific biokinetic data. Of course, the reliability of the latter estimate depends on the quality of the underlying biokinetic data.Conclusions: Finally, where feasible, the brain should be depicted explicitly in biokinetic models used in epidemiological studies addressing adverse effects of ionizing radiation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

A spectral deferred correction strategy for low Mach number reacting flows subject to electric fields

We propose an algorithm for low Mach number reacting flows subjected to electric field that includes the chemical production and transport of charged species. This work is an extension of a multi-implicit spectral deferred correction (MISDC) algorithm designed to advance the conservation equations in time at scales associated with advective transport. The fast and nontrivial interactions of electrons with the electric field are treated implicitly using a Jacobian-Free Newton Krylov approach for which a preconditioning strategy is developed. Within the MISDC framework, this enables a close and stable coupling of diffusion, reactions and dielectric relaxation terms with advective transport and is shown to exhibit second-order convergence in space and time. The algorithm is then applied to a series of steady and unsteady problems to demonstrate its capability and stability. Although developed in a one-dimensional case, the algorithmic ingredients are carefully designed to be amenable to multidimensional applications.

42 ENGINEERING↗
DOE OSTI2013

A NUMERICAL TREATMENT OF ANISOTROPIC RADIATION FIELDS COUPLED WITH RELATIVISTIC RESISTIVE MAGNETOFLUIDS

We develop a numerical scheme for solving fully special relativistic, resistive radiation magnetohydrodynamics. Our code guarantees conservation of total mass, momentum, and energy. The radiation energy density and the radiation flux are consistently updated using the M-1 closure method, which can resolve an anisotropic radiation field, in contrast to the Eddington approximation, as well as the flux-limited diffusion approximation. For the resistive part, we adopt a simple form of Ohm's law. The advection terms are explicitly solved with an approximate Riemann solver, mainly the Harten-Lax-van Leer scheme; the HLLC and HLLD schemes are also solved for some tests. The source terms, which describe the gas-radiation interaction and the magnetic energy dissipation, are implicitly integrated, relaxing the Courant-Friedrichs-Lewy condition even in an optically thick regime or a large magnetic Reynolds number regime. Although we need to invert 4 Multiplication-Sign 4 matrices (for the gas-radiation interaction) and 3 Multiplication-Sign 3 matrices (for the magnetic energy dissipation) at each grid point for implicit integration, they are obtained analytically without preventing massive parallel computing. We show that our code gives reasonable outcomes in numerical tests for ideal magnetohydrodynamics, propagating radiation, and radiation hydrodynamics. We also applied our resistive code to the relativistic Petschek-type magnetic reconnection, revealing the reduction of the reconnection rate via radiation drag.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

The classification of ambiguity in polarimetric reconstruction of coronal mass ejection

The Thomson scattering theory indicates that there exist explicit and implicit ambiguities in polarimetric analyses of coronal mass ejection (CME) observations. We suggest a classification for these ambiguities in CME reconstruction. Three samples, including double explicit, mixed, and double implicit ambiguity, are shown with the polarimetric analyses of STEREO CME observations. These samples demonstrate that this classification is helpful for improving polarimetric reconstruction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

AN ALGORITHM FOR RADIATION MAGNETOHYDRODYNAMICS BASED ON SOLVING THE TIME-DEPENDENT TRANSFER EQUATION

We describe a new algorithm for solving the coupled frequency-integrated transfer equation and the equations of magnetohydrodynamics in the regime that light-crossing time is only marginally shorter than dynamical timescales. The transfer equation is solved in the mixed frame, including velocity-dependent source terms accurate to O(v/c). An operator split approach is used to compute the specific intensity along discrete rays, with upwind monotonic interpolation used along each ray to update the transport terms, and implicit methods used to compute the scattering and absorption source terms. Conservative differencing is used for the transport terms, which ensures the specific intensity (as well as energy and momentum) are conserved along each ray to round-off error. The use of implicit methods for the source terms ensures the method is stable even if the source terms are very stiff. To couple the solution of the transfer equation to the MHD algorithms in the ATHENA code, we perform direct quadrature of the specific intensity over angles to compute the energy and momentum source terms. We present the results of a variety of tests of the method, such as calculating the structure of a non-LTE atmosphere, an advective diffusion test, linear wave convergence tests, and the well-known shadow test. We use new semi-analytic solutions for radiation modified shocks to demonstrate the ability of our algorithm to capture the effects of an anisotropic radiation field accurately. Since the method uses explicit differencing of the spatial operators, it shows excellent weak scaling on parallel computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Fully coupled simulation of cosmic reionization. I. numerical methods and tests

Here, we describe an extension of the Enzo code to enable fully coupled radiation hydrodynamical simulation of inhomogeneous reionization in large similar to(100 Mpc)(3) cosmological volumes with thousands to millions of point sources. We solve all dynamical, radiative transfer, thermal, and ionization processes self-consistently on the same mesh, as opposed to a postprocessing approach which coarse-grains the radiative transfer. But, we employ a simple subgrid model for star formation which we calibrate to observations. The numerical method presented is a modification of an earlier method presented in Reynolds et al. differing principally in the operator splitting algorithm we use to advance the system of equations. Radiation transport is done in the gray flux-limited diffusion (FLD) approximation, which is solved by implicit time integration split off from the gas energy and ionization equations, which are solved separately. This results in a faster and more robust scheme for cosmological applications compared to the earlier method. The FLD equation is solved using the hypre optimally scalable geometric multigrid solver from LLNL. By treating the ionizing radiation as a grid field as opposed to rays, our method is scalable with respect to the number of ionizing sources, limited only by the parallel scaling properties of the radiation solver. We test the speed and accuracy of our approach on a number of standard verification and validation tests. We show by direct comparison with Enzo's adaptive ray tracing method Moray that the well-known inability of FLD to cast a shadow behind opaque clouds has a minor effect on the evolution of ionized volume and mass fractions in a reionization simulation validation test. Finally, we illustrate an application of our method to the problem of inhomogeneous reionization in a 80 Mpc comoving box resolved with 3200(3) Eulerian grid cells and dark matter particles.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Shape-driven, emergent behavior in active particle mixtures

Abstract Active particle systems can vary greatly from one-component systems of spheres to mixtures of particle shapes at different composition ratios. We investigate computationally the combined effect of anisotropy and stoichiometry on the collective behavior of two-dimensional active colloidal mixtures of polygons. We uncover three emergent phenomena not yet reported in active Brownian particle systems. First, we find that mixtures containing hexagons exhibit micro-phase separation with large grains of hexagonal symmetry. We quantify a measurable, implicit ‘steric attraction’ between the active particles as a result of shape anisotropy and activity. This calculation provides further evidence that implicit interactions in active systems, even without explicit attraction, can lead to an effective preferential attraction between particles. Next, we report stable fluid clusters in mixtures containing one triangle or square component. We attribute the fluidization of the dense cluster to the interplay of cluster destabilizing particles, which introduce grain boundaries and slip planes into the system, causing solid-like clusters to break up into fluid clusters. Third, we show that fluid clusters can coexist with solid clusters within a sparse gas of particles in a steady state of three coexisting phases. Our results highlight the potential for a wide variety of behavior to be accessible to active matter systems and establish a route to control active colloidal systems through simple parameter designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

The different varieties of the Suyama-Yamaguchi consistency relation and its violation as a signal of statistical inhomogeneity

We present the different consistency relations that can be seen as variations of the well known Suyama-Yamaguchi (SY) consistency relation τ{sub NL}≥((6/5)f{sub NL}){sup 2}, the latter involving the levels of non-gaussianity f{sub NL} and τ{sub NL} in the primordial curvature perturbation ζ. It has been (implicitly) claimed that the following variation: τ{sub NL}(k{sub 1},k{sub 3})≥((6/5)){sup 2}f{sub NL}(k{sub 1})f{sub NL}(k{sub 3}), which we call ''the fourth variety'', in the collapsed (for τ{sub NL}) and squeezed (for f{sub NL}) limits is always satisfied independently of any physics; however, the proof depends sensitively on the assumption of scale-invariance (expressing this way the fourth variety of the SY consistency relation as τ{sub NL}≥((6/5)f{sub NL}){sup 2}) which only applies for cosmological models involving Lorentz-invariant scalar fields (at least at tree level), leaving room for a strong violation of this variety of the consistency relation when non-trivial degrees of freedom, for instance vector fields, are in charge of the generation of the primordial curvature perturbation. With this in mind as a motivation, we explicitly state, in the first part of this work, under which conditions the SY consistency relation has been claimed to hold in its different varieties (implicitly) presented in the literature since its inception back in 2008; as a result, we show for the first time that the variety τ{sub NL}(k{sub 1},k{sub 1})≥((6/5)f{sub NL}(k{sub 1})){sup 2}, which we call ''the fifth variety'', is always satisfied even when there is strong scale-dependence and high levels of statistical anisotropy as long as statistical homogeneity holds: thus, an observed violation of this specific variety would prevent the comparison between theory and observation, shaking this way the foundations of cosmology as a science. In the second part, we concern about the existence of non-trivial degrees of freedom, concretely vector fields for which the levels of non-gaussianity have been calculated for very few models; among them, and by making use of the δN formalism at tree level, we study a class of models that includes the vector curvaton scenario, vector inflation, and the hybrid inflation with coupled vector and scalar ''waterfall field'' where ζ is generated at the end of inflation, finding that the fourth variety of the SY consistency relation is indeed strongly violated for some specific wavevector configurations while the fifth variety continues to be well satisfied. Finally, as a byproduct of our investigation, we draw attention to a quite recently demonstrated variety of the SY consistency relation: τ{sup iso}{sub NL}≥((6/5)f{sup iso}{sub NL}){sup 2}, in scenarios where scalar and vector fields contribute to the generation of the primordial curvature perturbation; this variety of the SY consistency relation is satisfied although the isotropic pieces of the non-gaussianity parameters receive contributions from the vector fields. We discuss further implications for observational cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Locally adaptive activation functions with slope recovery for deep and physics-informed neural networks

Here we propose two approaches of locally adaptive activation functions namely, layer-wise and neuron-wise locally adaptive activation functions, which improve the performance of deep and physics-informed neural networks. The local adaptation of activation function is achieved by introducing a scalable parameter in each layer (layer-wise) and for every neuron (neuron-wise) separately, and then optimizing it using a variant of stochastic gradient descent algorithm. In order to further increase the training speed, an activation slope-based slope recovery term is added in the loss function, which further accelerates convergence, thereby reducing the training cost. On the theoretical side, we prove that in the proposed method, the gradient descent algorithms are not attracted to sub-optimal critical points or local minima under practical conditions on the initialization and learning rate, and that the gradient dynamics of the proposed method is not achievable by base methods with any (adaptive) learning rates. We further show that the adaptive activation methods accelerate the convergence by implicitly multiplying conditioning matrices to the gradient of the base method without any explicit computation of the conditioning matrix and the matrix–vector product. The different adaptive activation functions are shown to induce different implicit conditioning matrices. Furthermore, the proposed methods with the slope recovery are shown to accelerate the training process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Electromagnetic and neutral-weak response functions of 4 He and 12 C

A major goal of nuclear theory is to understand the strong interaction in nuclei as it manifests itself in terms of two- and many-body forces among the nuclear constituents, the protons and neutrons, and the interactions of these constituents with external electroweak probes via one- and many-body currents. The objective of the present work is to calculate the quasielastic electroweak response functions in light nuclei within the realistic dynamical framework outlined above. These response functions determine the inclusive cross section as function of the lepton momentum and energy transfers. Their ab initio calculation is a very challenging quantum many-body problem, since it requires summation over the entire excitation spectrum of the nucleus and inclusion in the electroweak currents of one- and many-body terms. Green's functions Monte Carlo methods allow one to circumvent both difficulties by computing the response in imaginary time (the so-called Euclidean response) and hence summing implicitly over the bound and continuum states of the nucleus, and by implementing specific algorithms designed to deal with the complicated spin-isospin structure of nuclear many-body operators. Theoretical predictions for 4 He and 12 C, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasielastic to the dip region and beyond. In conclusion, these results challenge the conventional picture of quasielastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Fast explicit solutions for neutrino-electron scattering: Explicit asymptotic methods

Here, we present results of explicit asymptotic approximations applied to neutrino-electron scattering in a representative model of neutrino population evolution under conditions characteristic of core-collapse supernova explosions or binary neutron star mergers. It is shown that this approach provides stable solutions of these stiff systems of equations, with accuracy and time stepping comparable to that for standard implicit treatments such as backward Euler, fixed point iteration, and Anderson-accelerated fixed point iteration. Because each time step can be computed more rapidly with the explicit asymptotic approximation than with implicit methods, this suggests that algebraically stabilized explicit integration methods could be used to compute neutrino evolution coupled to hydrodynamics more efficiently in stellar explosions and mergers than the methods currently in use.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A Big Data Analytics Framework for HPC Log Data: Three Case Studies Using the Titan Supercomputer Log

Reliability, availability and serviceability (RAS) logs of high performance computing (HPC) resources, when closely investigated in spatial and temporal dimensions, can provide invaluable information regarding system status, performance, and resource utilization. These data are often generated from multiple logging systems and sensors that cover many components of the system. The analysis of these data for finding persistent temporal and spatial insights faces two main difficulties: the volume of RAS logs makes manual inspection difficult and the unstructured nature and unique properties of log data produced by each subsystem adds another dimension of difficulty in identifying implicit correlation among recorded events. To address these issues, we recently developed a multi-user Big Data analytics framework for HPC log data at Oak Ridge National Laboratory (ORNL). This paper introduces three in-progress data analytics projects that leverage this framework to assess system status, mine event patterns, and study correlations between user applications and system events. We describe the motivation of each project and detail their workflows using three years of log data collected from ORNL's Titan supercomputer.

Park, Byung↗
DOE OSTI2015

Grandmaster: Interactive text-based analytics of social media

People use social media resources like Twitter, Facebook, forums etc. to share and discuss various activities or topics. By aggregating topic trends across many individuals using these services, we seek to construct a richer profile of a person’s activities and interests as well as provide a broader context of those activities. This profile may then be used in a variety of ways to understand groups as a collection of interests and affinities and an individual’s participation in those groups. Our approach considers that much of these data will be unstructured, free-form text. By analyzing free-form text directly, we may be able to gain an implicit grouping of individuals with shared interests based on shared conversation, and not on explicit social software linking them. In this paper, we discuss a proof-of-concept application called Grandmaster built to pull short sections of text, a person’s comments or Twitter posts, together by analysis and visualization to allow a gestalt understanding of the full collection of all individuals: how groups are similar and how they differ, based on their text inputs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A substantial role of soil erosion in the land carbon sink and its future changes

Realistic representation of land carbon sink in climate models is vital for predicting carbon climate feedbacks in a changing world. Although soil erosion that removes land organic carbon has increased substantially since the onset of agriculture, it is rarely included in the current generation of climate models. Using an Earth system model with soil erosion represented, we estimated that on average soil erosion displaces 5% of newly fixed land organic carbon downslope annually in the continental US. In the lower Mississippi river basin and the Cascades, the fraction can be as large as 40%. About 12% of the eroded organic carbon is eventually exported to inland waters, which is equal to 14% of the simulated net carbon gain by terrestrial ecosystems. By comparing the eroded organic carbon export to rivers with the particulate organic carbon export to oceans, we demonstrated that a large fraction of the carbon export to rivers could have been mineralized in inland waters. Importantly, with a direct comparison of eroded and exported soil organic carbon and land net carbon uptake, we found that Earth system models that ignore soil erosion likely offset the erosional carbon loss by increasing heterotrophic respiration implicitly. But as soil erosion and heterotrophic respiration respond differently to a warming climate, this unrealistic compensation would lead to biased predictions of future land carbon sink.

Tan, Zeli↗