Search NASA⌕ Search

SEARCH · Search NASA

Results for “Yield function”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Measurement of the production of (anti)nuclei in p–Pb collisions at s NN = 8.16 TeV

Measurements of (anti)proton, (anti)deuteron, and (anti)3He production in the rapidity range -1< y < 0 as a function of the transverse momentum and event multiplicity in p–Pb collisions at a center-of-mass energy per nucleon–nucleon pair $\sqrt{S_{NN}}$ = 8.16 TeV are presented. The coalescence parameters B 2 and B 3 , measured as a function of the transverse momentum per nucleon and of the mean charged-particle multiplicity density, confirm a smooth evolution from low to high multiplicity across different collision systems and energies. The ratios between (anti)deuteron and (anti) 3 He yields and those of (anti)protons are also reported as a function of the mean charged-particle multiplicity density. A comparison with the predictions of the statistical hadronization and coalescence models for different collision systems and center-of-mass energies favors the coalescence description for the deuteron-to-proton yield ratio with respect to the canonical statistical model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Multiplicity dependence of f 0 (980) production in pp collisions at $\mathbf{\sqrt{s}}$ = 13 TeV

The dependence of f 0 (980) production on the final-state charged-particle multiplicity is reported for proton–proton (pp) collisions at the centre-of-mass energy, $\mathbf{\sqrt{s}}$ = 13 TeV. The production of (980) is measured with the ALICE detector via the decay channel in a midrapidity region of |y| < 0.5. The evolution of the integrated yields and mean transverse momentum of f 0 (980) as a function of charged-particle multiplicity measured in pp at $\mathbf{\sqrt{s}}$ = 13 TeV follows the trends observed in pp at $\mathbf{\sqrt{s}}$ = 5.02 TeV and in proton–lead (p–Pb) collisions at$\mathbf{\sqrt{s_{NN}}}$ = 5.02 TeV. Particle yield ratios of f 0 (980) to π ± and K* (892) are found to decrease with increasing charged-particle multiplicity. These particle ratios are compared with calculations from the canonical statistical thermal model as a function of charged-particle multiplicity. The thermal model calculations provide a better description of the decreasing trend of particle ratios when no strange or antistrange quark composition for f 0 (980) is assumed, which suggests that the data do not support significant hidden strangeness in the f 0 (980).

Abualrob, I. J. [University of Houston] (ORCID:000↗
DOE OSTI2023

High-resolution vibrational predissociation spectroscopy of I - ·H 2 O by single-mode CW infrared excitation in a 3D cryogenic ion trap

We describe the integration of a tunable, single-mode, continuous wave infrared laser into a cryogenic ion spectroscopy experiment to measure the rovibrational spectrum of the I - ∙H 2 O complex in the OH stretching region. These upper levels lie about 300 cm -1 above the dissociation threshold. The measurements are carried out by loading the ions in a radiofrequency ion trap at 10 Hz and cooling them to 5 K with pulsed He buffer gas. IR photodissociation (PD) of the I - ∙H 2 O complex is monitored by recording the I - product yield by time-of-flight mass spectrometry as a function of laser wavelength. Very narrow (Δν~75 MHz) rotational lines are observed throughout the spectrum, indicating long (ca. 2 ns) lifetimes for the excited metastable rovibrational levels. Rotational analysis of the band arising from the K”=1 to K’=2 transition of the free OH stretching fundamental yields the structure of the complex for the first time. Over 50% of the trapped ion ensemble in the trap can be photodissociated upon excitation of a single rotational line. This enables very high signal-to-noise in the PD spectrum, and is traced to a mechanism in which the ground state rotational levels are rapidly equilibrated by collisions with the buffer gas.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Modeling the one-dimensional imager of neutrons (ODIN) for neutron response functions at the Sandia Z facility

The one-dimensional imager of neutrons (ODIN) at the Sandia Z facility consists of a 10-cm block of tungsten with rolled edges, creating a slit imager with slit widths of either 250, 500, or 750 μm. Designed with a 1-m neutron imaging line of sight, we achieve about 4:1 magnification and 500-μm axial spatial resolution. The baseline inertial confinement fusion concept at Sandia is magnetized liner inertial fusion, which nominally creates a 1-cm line source of neutrons. ODIN was designed to determine the size, shape, and location of the neutron producing region, furthering the understanding of compression quality along the cylindrical axis of magnetized liner implosions. Challenges include discriminating neutrons from hard x-rays and gammas with adequate signal-to-noise in the 2 × 10 12 deuterium-deuterium (DD) neutron yield range, as well as understanding the point spread function of the imager to various imaging detectors (namely, CR-39). Modeling efforts were conducted with MCNP6.1 to determine neutron response functions for varying configurations in a clean DD neutron environment (without x-rays or gammas). Configuration alterations that will be shown include rolled-edge slit orientation and slit width, affecting the resolution and response function. The experiment to determine CR-39 neutron sensitivity, with and without a high density polyethylene (n, p) converter, an edge spread function, and resolution will be discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Synthesis of imidazolium-mediated Poly(benzoxazole) Ionene and composites with ionic liquids as advanced gas separation membranes

We report thermally rearranged (TR) polymers and ionic polymers are two material classes which have been employed in leading gas separation membranes. This work introduces a novel approach of combining the benzoxazole functionality associated with TR polymers with tailorable cationic groups, yielding a new type of imidazolium-mediated poly(benzoxazole) ionene polymer, “Im-PBO-Ionene” with the aim of enhanced CO 2 separation performance. The structural changes exhibited from the Coulombic interactions between the ionene backbone and the “free” ionic liquid (IL) resulted in enhanced gas separation performance, shown in fundamental characterizations, observed through increased diffusivities and more notably, retained high selectivities and 3x or 5x respective increases in CO 2 permeability upon the addition of 1 or 2 equivalents of IL per polymer repeat unit. These new high-performance ionenes demonstrate the versatility of ionene design and potential of the ionene + IL material platform for gas separation membranes with versatile incorporation of sophisticated functional and structural features.

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

Eta decay and muonic puzzles

New physics motivated by muonic puzzles (proton radius and muon g - 2 discrepancies) is studied. Using a light scalar boson Φ, assuming Yukawa interactions, accounts for these muonic puzzles simultaneously. Our previous work limits the existence of such a scalar boson's mass m Φ from about 160 keV to 60 MeV. We improve this result by including the influence of all of the possible particles that couple to the Φ in computing the decay rate. Doing this involves including the strong interaction physics, involving quarks, necessary to compute the ηπΦ vertex function. The Nambu-Jona-Lasinio model, which accounts for the spontaneous symmetry breaking that yields the constituent mass is employed to represent the relevant strong-interaction physics. We use the ηπΦ vertex function to reanalyze the electron beam dump experiments. The result is that the allowed range of m Φ lies between about 160 keV and 3.5 MeV. This narrow range represents an inviting target for ruling out or discovering this scalar boson. A possible UV completion of our phenomenological model is discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Capturing multireference excited states by constrained-density-functional theory

The computation of excited electronic states with commonly employed (approximate) methods is challenging, typically yielding states of lower quality than the corresponding ground state for a higher computational cost. In this work, we present a mean-field method that extends the previously proposed excited constrained-density-functional theory (XCDFT) from single Slater determinants to ensemble one-body reduced density matrices for computing low-lying excited states. The method still retains an associated computational complexity comparable to a semilocal density functional theory (DFT) calculation while at the same time it is capable of approaching states with multireference character. Further, we benchmark the quality of this method on well-established test sets, finding good descriptions of the electronic structure of multireference states and maintaining an overall accuracy for the predicted excitation energies comparable to semilocal time-dependent DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Report on Aging and Irradiation Theoretical Model Development for Alloy 690 in Grizzly

This report describes the current efforts to implement a model to predict the engineering scale behavior of irradiated and aged Ni 2 Cr model alloy and Alloy 690 into the Grizzly software. Grizzly is a multiphysics simulation code under development at Idaho National Laboratory to simulate aging mechanisms and their effects on the integrity of critical LWR components. Engineering scale behavior of materials are strongly impacted by the underlying microstructure: the distribution and interaction of dislocations and defects within the individual grains dictates how the material will respond to loading. The evolution of these dislocations and defects determines how the engineering scale material will perform under different loading and environmental conditions. Mesoscale models can function as a bridge to connect microstructure evolution models to engineering models of the ductile-brittle transition temperature (DBTT) curve and fracture. The toughness of a material, which governs the engineering scale fracture behavior, is a function of both elastic deformation and plastic deformation. While the yield stress measurements capture the elastic contribution to a material’s toughness, information about the dislocation movement is required to understand the plastic component of toughness. Crystal plasticity models predict dislocation density evolution and thereby connect microstructure evolution to engineering scale behavior predictions. These models calculate the plastic strain within a metal as a function of dislocations and other crystal defects, including precipitates and interstitial loops, among others, by tracking dislocation movement and the interaction of these dislocations with the crystal defects. Within NiCr alloys, long-range ordered (LRO) precipitates form under aging conditions. These precipitates act as barriers to dislocation motion and are either sheared by dislocations or pin dislocations. Additional defects, such as stacking fault tetrahedra (SFTs) from irradiation damage, also impede dislocation motion in these alloys. A crystal plasticity mode for Ni 2 Cr model alloy and Alloy 690 is being currently implemented in the Grizzly code at INL. This report includes a description of the theoretical models incorporated into the Grizzly crystal plasticity code. These theoretical models include mechanisms for dislocation glade and twinning as well as dislocation-defect interaction relationship for both LRO precipitates and SFTs. We have selected these dislocation models based on the deformation mechanisms observed in experimental micropillar compression tests conducted on the Ni 2 Cr model alloy. The models are then extended to Alloy 690 with ongoing parameter fitting.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Exchange-correlation energies of atoms from efficient density functionals: influence of the electron density

The primary goal of Kohn–Sham density functional theory is to evaluate the exchange-correlation contribution to electronic properties. However, the accuracy of a density functional can be affected by the electron density. Here we apply the nonempirical Tao–Mo (TM) semilocal functional to study the influence of the electron density on the exchange and correlation energies of atoms and ions, and compare the results with the commonly used nonempirical semilocal functionals local spin-density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE), Tao–Perdew–Staroverov–Scuseria (TPSS), and hybrid functional PBE0. We find that the spin-restricted Hartree–Fock density yields the exchange and correlation energies in good agreement with the Optimized Effective Potential method, particularly for spherical atoms and ions. However, the errors of these semilocal and hybrid functionals become larger for self-consistent densities. We further find that the quality of the electron density have greater effect on the exchange-correlation energies of kinetic energy density-dependent meta-GGA functionals TPSS and TM than on those of the LSDA and GGA, and therefore, should have greater influence on the performance of meta-GGA functionals. Lastly, we show that the influence of the density quality on PBE0 is slightly reduced, compared to that of PBE, due to the exact mixing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Exciting H 2 Molecules for Graphene Functionalization

Hydrogen functionalization of graphene by exposure to vibrationally excited H 2 molecules is investigated by combined scanning tunneling microscopy, high-resolution electron energy loss spectroscopy, X-ray photoelectron spectroscopy measurements, and density functional theory calculations. The measurements reveal that vibrationally excited H 2 molecules dissociatively adsorb on graphene on Ir(111) resulting in nanopatterned hydrogen functionalization structures. Calculations demonstrate that the presence of the Ir surface below the graphene lowers the H 2 dissociative adsorption barrier and allows for the adsorption reaction at energies well below the dissociation threshold of the H–H bond. The first reacting H 2 molecule must contain considerable vibrational energy to overcome the dissociative adsorption barrier. However, this initial adsorption further activates the surface resulting in reduced barriers for dissociative adsorption of subsequent H 2 molecules. This enables functionalization by H 2 molecules with lower vibrational energy, yielding an avalanche effect for the hydrogenation reaction. These results provide an example of a catalytically active graphene-coated surface and additionally set the stage for a re-interpretation of previous experimental work involving elevated H 2 background gas pressures in the presence of hot filaments.

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

Explosive Yield Estimation Using Regional Seismic Moment Tensors

Here, we use the Pasyanos and Chiang (2022) data set to calculate the seismic moment M 0 for each explosion and use the measured explosive yield W to validate the W~M 0 relationship in Denny and Johnson (1991; hereafter, DJ91). The M 0 is corrected by transforming to a potency tensor and applying more appropriate near-source geophysical parameter values in the moment estimate. The mean residual between observed and predicted yield is near zero; however, the standard deviation of the residuals results in an F-value (a 95% confidence factor) of about 5. We re-estimate the coefficients in the DJ91 model and find similar values and only a slight improvement in the F-value. Next, we embark on a similar model selection process as DJ91, allowing for non-cube-root yield scaling and other plausible near-source elastic moduli. As was found by DJ91, the yield dependence is not significantly different from unity, and a cube root assumption is valid. Therefore, we yield scale the seismic moment and test the significance of all plausible explanatory variables. Isotropic moment performs better in the response variable than total moment. The preference for isotropic moment could be due to its relationship to volume change, which would be more directly affected by explosive yield. Surprisingly, we find that the overburden pressure, which is a function of depth, is not a significant parameter in the model. We hypothesize that this is due to the competing depth effects on source asymmetry and the incorporation of depth in the Green’s functions used to calculate the seismic moment tensors. Importantly, this emphasizes that only seismic moment tensor-derived moments should be used in these models. After removing insignificant model parameters, we are left with a simple model to predict explosive yield $\widehat{W}$ in kt from isotropic moment M I in N·m, $\widehat{W}$=κ –1.4132 10 0.035626GP M I , in which κ and GP are the near-source bulk modulus and gas porosity in Pa and %, respectively. The F-value for this model is approximately 3.

58 GEOSCIENCES↗
DOE OSTI2017

Characteristics of Solar Energetic Ions as a Function of Longitude

Since the 2006 launch of STEREO , multi-spacecraft studies have yielded several surprising results regarding the spread of solar energetic particles (SEPs) within the inner heliosphere. We have investigated the role of energy and ridigity, using ACE and STEREO 10 MeV n{sup −1} oxygen data to identify 41 large SEP events observed by two or three spacecraft. We calculated fluence spectra from ∼0.1 to >10 MeV n{sup −1} for H, He, O, and Fe for each event at the observing spacecraft (including SOHO and GOES ). The particle fluences at 0.3, 1, and 10 MeV n{sup −1} were examined as a function of the distance between the associated solar flare longitude and the spacecraft magnetic footpoints at the Sun to determine the longitudinal spread of particles and study how the distribution centers and widths depend on energy and charge-to-mass (Q/M) for the first time. On average, the three-spacecraft event distributions were centered at 22 ± 4° west of the flare site and were 43 ± 1° wide, though there was substantial variability, while the fit to the aggregate of the two-spacecraft event fluences yielded significantly wider distributions at 0.3 and 1 MeV n{sup −1}. The widths derived from both the three- and two-spacecraft events show an energy dependence with distributions narrowing with increasing energy, consistent with lower energy ions experiencing more field line co-rotation, or being accelerated over a larger portion of the CME-driven shock or for longer times as the shock expands. Surprisingly, no clear evidence was found for a Q/M dependence to the widths or centers suggesting that rigidity-related processes are not the dominant means of spreading particles in longitude.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

CYP76BK1 orthologs catalyze furan and lactone ring formation in clerodane diterpenoids across the mint family

The Lamiaceae (mint family) is the largest known source of furanoclerodanes, a subset of clerodane diterpenoids with broad bioactivities including insect antifeedant properties. The Ajugoideae subfamily, in particular, accumulates significant numbers of structurally related furanoclerodanes. The biosynthetic capacity for formation of these diterpenoids is retained across most Lamiaceae subfamilies, including the early-diverging Callicarpoideae which forms a sister clade to the rest of Lamiaceae. VacCYP76BK1, a cytochrome P450 monooxygenase from Vitex agnus-castus, was previously found to catalyze the formation of the proposed precursor to furan and lactone-containing labdane diterpenoids. Through transcriptome-guided pathway exploration, we identified orthologs of VacCYP76BK1 in Ajuga reptans and Callicarpa americana. Functional characterization demonstrated that both could catalyze the oxidative cyclization of clerodane backbones to yield a furan ring. Subsequent investigation revealed a total of 10 CYP76BK1 orthologs across six Lamiaceae subfamilies. Through analysis of available chromosome-scale genomes, we identified four CYP76BK1 members as syntelogs within a conserved syntenic block across divergent subfamilies. This suggests an evolutionary lineage that predates the speciation of the Lamiaceae. Functional characterization of the CYP76BK1 orthologs affirmed conservation of function, as all catalyzed furan ring formation. Additionally, some orthologs yielded two novel lactone ring moieties. The presence of the CYP76BK1 orthologs across Lamiaceae subfamilies closely overlaps with the distribution of reported furanoclerodanes. Together, the activities and distribution of the CYP76BK1 orthologs identified here support their central role in furanoclerodane biosynthesis within the Lamiaceae family. Our findings lay the groundwork for biotechnological applications to harness the economic potential of this promising class of compounds.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Maize yield under a changing climate: The hidden role of vapor pressure deficit

Temperatures over the next century are expected to rise to levels detrimental to crop growth and yield. As the atmosphere warms without additional water vapor input, vapor pressure deficit (VPD) increases as well. Increased temperatures and accompanied elevated VPD levels can both lead to negative impacts on crop yield. The independent importance of VPD, however, is often neglected or conflated with that from temperature due to a tight correlation between the two climate factors. We used a coupled process-based crop (MAIZSIM) and soil (2DSOIL) model to gain a mechanistic understanding of the independent roles temperature and VPD play in crop yield projections, as well as their interactions with rising CO 2 levels and changing precipitation patterns. We found that by separating out the VPD effect from rising temperatures, VPD increases had a greater negative impact on yield (12.9 ± 1.8%, increase in VPD associated with 2 °C warming) compared to that from warming (8.5 ± 1.4%, the direct effect of 2 °C warming). The negative impact of these two factors varied with precipitation levels and influenced yield through separate mechanisms. Warmer temperatures caused yield loss mainly through shortening the growing season, while elevated VPD increased water loss and triggered several water stress responses such as reduced photosynthetic rates, lowered leaf area development, and shortened growing season length. Elevated CO 2 concentrations partially alleviated yield loss under warming or increased VPD conditions through water savings, but the impact level varied with precipitation levels and was most pronounced under drier conditions. These results demonstrate the key role VPD plays in crop growth and yield, displaying a magnitude of impact comparative to temperature and CO 2 . A mechanistic understanding of the function of VPD and its relation with other climate factors and management practices is critical to improving crop yield projections under a changing climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

High-shade dryland agrivoltaic conditions enhanced carbon uptake and water-use efficiency in zucchini ( Cucurbita pepo )

Introduction: The increasing global demand for food and energy is intensifying land-use competition. Agrivoltaic systems are a multifunctional land-use approach that vertically integrates the production of agricultural crops and solar power on the same land area. Most food crops are adapted to full-sun conditions, and the physiological responses of these crops to the novel microclimate under solar panels remain poorly understood. We hypothesized that the microclimate beneath the high-density photovoltaic system would influence carbon uptake, water use, and yield outcomes of zucchini summer squash.Methods: We conducted a field experiment in a hot, semi-arid climate on zucchini (Cucurbita pepo). Plants were grown under an agrivoltaic system with a 75% ground cover ratio (GCR) and in a full-sun control plot, each with two irrigation regimes (100 and 50%). We measured leaf-level photosynthesis, microclimate variables, and fruit yield at plant maturity and throughout the growing season.Results: The agrivoltaic array reduced photosynthetically active radiation (PAR) by ~79%, resulting in a cooler (−1.1 °C), more humid environment with higher soil moisture. These microclimatic conditions enhanced midday photosynthesis and daily cumulative carbon uptake. However, fruit yield was consistently lower under the panels, indicating a shift in carbon allocation toward vegetative growth. Photosynthesis was primarily driven by PAR across treatments, while soil moisture significantly influenced photosynthesis only in the control plots, suggesting water limitation was alleviated under the panels.Discussion: These findings highlight a trade-off between improved physiological performance and reduced yield under high-density agrivoltaics. While the system buffered heat and drought stress and improved overall plant function, excessive shade reduced reproductive output. Optimizing panel density or selecting crops cultivated for non-fruit yields will be essential for balancing food production and energy generation in dryland agrivoltaic settings.

14 SOLAR ENERGY↗
DOE OSTI2017

Uranyl nitrate hexahydrate solubility in nitric acid and its crystallization selectivity in the presence of nitrate salts

The solubility of uranyl nitrate hexahydrate was determined as a function of nitric acid concentration and temperature, and the crystallization yield was calculated. Results showed an increase in crystal formation at lower initial nitric acid concentrations upon cooling a saturated solution of U(VI) from 50 °C to 2 °C, with over 70% recovery of U(VI) mass at all nitric acid concentrations and nearly quantitative recovery starting at 4 M HNO 3 . A direct correlation between the change in mother liquor volume percent and U mass removal percent was observed. By reducing the cooling rate from roughly 4.0 °C/min to 0.22 °C/min, the separation factor was increased from 3.88 to 15.7 to greater than 81 for the separation of U(VI) from Sr, Cs, and Nd. At the slower cooling rate, the separation factors were measured as a function of acidity for 2.0–4.3 M HNO 3 , showing a decrease in selectivity with a decrease in the acidity. There was also no indication that tetravalent metal double-salt precipitation occurred with either Zr 4+ or Ce 4+ . Here, these results indicate that a high-yield, high purity hexavalent actinide crystallization scheme may offer attractive benefits for nuclear-fuel recycle in that only a single very simple and well-understood technology is employed, and the use of organic compounds and solvents is avoided.

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

Milestone 1.2.9: Radiolytic Gas Generation Measurements from Helium Backfilled Samples of AA1100 and AA6061 Coupons

Nearly 18 metric tons of aluminum-clad spent nuclear fuel (ASNF) is safely managed by the U.S. Department of Energy (DOE). These assemblies are currently in interim storage, with the intention of extended storage (>50 years) until final disposal. Strategies for the continued safe storage of this material are under evaluation, of which a key criterion is the extent of molecular hydrogen gas (H2) formation from the radiolysis of hydrated (oxy)hydroxide aluminum corrosion layers arising from in-reactor and wet storage conditions. Radiation-induced H2 formation has the potential to compromise cladding and storage canister integrity, in addition to promoting the formation of unfavorable gaseous environments. Consequently, understanding this radiation-induced phenomenon is essential for the development of predictive modeling capabilities to support technical considerations and the identification of radiation related challenges for the extended storage of ASNF. Here, we report radiolytic H2 yields (G-values, G(H2)) from the gamma irradiation of ‘pristine’ and pre-corroded aluminum coupons in helium (He) environments as a function of alloy composition (AA1100 and AA6061), relative humidity, and absorbed gamma dose. Measured yields were lower than corresponding values reported for argon environments, a positive result for proposed extended dry storage strategies that would employ helium as a backfill gas. Interestingly, the presented G(H2)He values are comparable to those previously measured in nitrogen environments, suggesting a He mediated H2 inhibition process, attributed here to Penning ionization.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

Artificial reforestation produces less diverse soil nitrogen–cycling genes than natural restoration

Reforestation is effective in restoring ecosystem functions and enhancing ecosystem services of degraded land. The three most commonly employed reforestation methods of natural reforestation, artificial reforestation with native Masson pine, and introduced slash pine plantations were equally successful in biomass yield in southern China. However, it is not known whether soil ecosystem functions, such as nitrogen (N) cycling, are also successfully restored. Here, we employed a functional microarray to illustrate soil N-cycling. The composition of N-cycling genes in soils varied significantly with reforestation method and varied with constructive species identity and plant diversity. Artificial reforestation had less superior organization of N-cycling genes than natural reforestation, as indicated by the less diverse and less stable pathways to perform the biogeochemical N-cycling processes. Besides, artificial reforestation had lower functional potential (abundance of ammonification, denitrification, assimilatory, and dissimilatory nitrate reduction to ammonium genes) in soils than natural method. Evaluations of the abundance and interactions of N-cycling genes in soils showed that plantations, especially artificial reforestation with slash pine plantations, possessed a smaller range of ecosystem functions that provide a less diverse array of N-related substrates and nutrients to microbial communities compared with natural restoration. This might lead to a lower independence of N-cycling, which indicated a higher risk of N release in plantations. The unfavorable N-cycling conditions in plantations were corroborated by the lower contents of available N, ammonium N, and nitrate N. These findings demonstrate that reforestation methods could have broad regional and possibly global implications for N-cycling.

59 BASIC BIOLOGICAL SCIENCES↗