Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Isotope impact on Alfvén eigenmodes and fast ion transport in DIII-D

Abstract Measurements of beam driven Alfvén Eigenmode (AE) activity in matched deuterium (D) and hydrogen (H) DIII-D plasmas show a dramatic difference in unstable mode activity and fast ion transport for a given injected beam power. The dependence of the unstable AE spectrum in reversed magnetic shear plasmas on beam and thermal species is investigated in the current ramp by varying beam power in a sequence of discharges for fixed thermal and beam species at fixed density. In general, a spectrum of Reversed Shear Alfvén Eigenmodes (RSAEs) and Toroidal Alfvén Eigenmodes (TAEs) are driven unstable with sub-Alfvénic D beam injection while primarily only RSAEs are driven unstable for the H beam cases investigated. Further, for a given beam power, the driven AE amplitude is always reduced with H beams relative to D and for H thermal plasma relative to pure D or mixed D/H plasmas. Estimates of the fast ion stored energy combined with modeling using the hybrid kinetic-MHD code MEGA indicate that the dominant mechanism contributing to the difference between H and D beam drive is the faster classical slowing down of H beam ions relative to D and the resultant lower beam ion pressure. Calculations of the AE induced stored energy deficits using the reduced critical gradient model TGLFEP show quantitative agreement with the observed dependencies on injected power, isotope and minimum safety factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Ionization and temperature measurements in warm dense copper using x-ray absorption spectroscopy

We detail experimental results inferring ionization and temperature for warm dense copper plasmas at several times solid density (15–25 g/cm 3 ) and temperatures of 10–21 eV. Experiments performed at the OMEGA Laser Facility generate uniform warm dense matter conditions via symmetric shock compression of a buried copper layer. The plasma is probed with a laser-generated x-ray source to collect the K-shell x-ray absorption spectrum. Fitting bound-bound absorption contributions from constituent charge states of copper provides an estimated $\overline{Z}$ of approximately 4–7 for these warm dense copper plasmas. We find that these partially ionized plasmas have K-edge shifts of 12–30 eV and bound-bound resonance 1s → 3p absorption shifts of 4–26 eV with respect to the cold K edge. This study provides necessary experimental data to improve ionization and opacity models in the warm dense matter regime.

42 ENGINEERING↗
DOE OSTI2023

Charge order in the kagome lattice Holstein model: a hybrid Monte Carlo study

The Holstein model is a paradigmatic description of the electron-phonon interaction, in which electrons couple to local dispersionless phonon modes, independent of momentum. The model has been shown to host a variety of ordered ground states such as charge density wave (CDW) order and superconductivity on several geometries, including the square, honeycomb, and Lieb lattices. In this work, we study CDW formation in the Holstein model on the kagome lattice, using a recently developed hybrid Monte Carlo simulation method. We present evidence for √3 × √3 CDW order at an average electron filling of < n > = 2/3 per site, with an ordering wavevector at the K-points of the Brillouin zone. We estimate a phase transition occurring at T c ≈ t/18, where t is the nearest-neighbor hopping parameter. Our simulations find no signature of CDW order at other electron fillings or ordering momenta for temperatures T ≥ t/20.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electrostatic gate-controlled quantum interference in a high-mobility two-dimensional electron gas at the (La 0.3 ⁢Sr 0.7 )(Al 0.65 ⁢Ta 0.35 )O 3 /SrTiO 3 interface

Here, we report quantum oscillations in magnetoresistance that are periodic in magnetic field (𝐵), observed at the interface between (La 0.3 ⁢Sr 0.7 )(Al 0.65 ⁢Ta 0.35 )O 3 and SrTiO 3 . Unlike Shubnikov–de Haas oscillations, which appear at magnetic fields > 7 T and diminish quickly as the temperature rises, these 𝐵-periodic oscillations emerge at low fields and persist up to 10 K. Their amplitude decays exponentially with both temperature and field, specifying dephasing of quantum interference. Increasing the carrier density through electrostatic gating results in a systematic reduction in both the amplitude and frequency of the oscillations, with complete suppression beyond a certain gate voltage. We attribute these oscillations to the Altshuler-Aronov-Spivak effect, likely arising from naturally formed closed-loop paths due to the interconnected quasi-one-dimensional conduction channels along SrTiO 3 domain walls. The relatively long phase coherence length (≃ 1.8 µ⁢m at 0.1 K), estimated from the oscillation amplitude, highlights the potential of complex oxide interfaces as a promising platform for exploring quantum interference effects and advancing device concepts in quantum technologies, such as mesoscopic interferometers and quantum sensors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The Anatomy of an Unusual Edge-on Protoplanetary Disk. I. Dust Settling in a Cold Disk

As the earliest stage of planet formation, massive, optically thick, and gas-rich protoplanetary disks provide key insights into the physics of star and planet formation. When viewed edge-on, high-resolution images offer a unique opportunity to study both the radial and vertical structures of these disks and relate this to vertical settling, radial drift, grain growth, and changes in the midplane temperatures. In this work, we present multi-epoch Hubble Space Telescope and Keck scattered light images, and an Atacama Large Millimeter/submillimeter Array 1.3 mm continuum map for the remarkably flat edge-on protoplanetary disk SSTC2DJ163131.2–242627, a young solar-type star in ρ Ophiuchus. We model the 0.8 μm and 1.3 mm images in separate Markov Chain Monte Carlo (MCMC) runs to investigate the geometry and dust properties of the disk using the MCFOST radiative transfer code. In scattered light, we are sensitive to the smaller dust grains in the surface layers of the disk, while the submillimeter dust continuum observations probe larger grains closer to the disk midplane. An MCMC run combining both data sets using a covariance-based log-likelihood estimation was marginally successful, implying insufficient complexity in our disk model. The disk is well characterized by a flared disk model with an exponentially tapered outer edge viewed nearly edge-on, though some degree of dust settling is required to reproduce the vertically thin profile and lack of apparent flaring. A colder than expected disk midplane, evidence for dust settling, and residual radial substructures all point to a more complex radial density profile to be probed with future, higher-resolution observations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Combining Planck and SPT Cluster Catalogs: Cosmological Analysis and Impact on the Planck Scaling Relation Calibration

We provide the first combined cosmological analysis of the South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the νΛ CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find ${\alpha }_{\mathrm{SZ}}={1.49}_{-0.10}^{+0.07}$ and ${\left(1-b\right)}_{\mathrm{SZ}}={0.69}_{-0.14}^{+0.07}$, respectively. The results for the mass slope show a ~4 σ departure from the self-similar evolution, α SZ ~ 1.8. This shift is mainly driven by the matter density value preferred by SPT data, Ωm = 0.30 ± 0.03, lower than the one obtained by Planck data alone, ${{\rm{\Omega }}}_{m}={0.37}_{-0.06}^{+0.02}$. The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, (1 – b) ~ 0.8, and with results required by the Planck cosmic microwave background cosmology, (1 – b) ~ 0.6. From this analysis, we obtain a new catalog of Planck cluster masses M 500 . We estimate the ratio between the published Planck M SZ masses and our derived masses M 500 , as a "measured mass bias," ${\left(1-b\right)}_{M}$. We analyze the mass, redshift, and detection noise dependence of ${\left(1-b\right)}_{M}$, finding an increasing trend toward high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass, high-mass, low-z, and high-z regimes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Elucidation of Structure and Physical Properties of Pyrolytic Sugar Oligomers Derived from Cellulose Depolymerization/Dehydration Reactions: A Density Functional Theory Study

Fast pyrolysis of lignocellulosic materials is a promising research area to produce renewable fuels and chemicals. Dehydration is known to be among the most important reaction families during cellulose pyrolysis; water is the most important product. Together with water, dehydration reactions also form a range of poorly known oligomer species of varying molecular sizes, often collected as part of bio-oil water-soluble (WS) fraction. In this work, we used electronic structure calculations to evaluate the relative thermodynamic stabilities of several oligomer species resulting from up to three consecutive dehydration events from cellulose depolymerization intermediates. A library of the thermodynamically favored candidate molecular structures was compiled. Results revealed that most of the water molecules are eliminated from the non-reducing end, forming thermodynamically more stable conjugated compounds. This is consistent with results reported by other researchers in literature where dehydration reactions occur preferably at the non-reducing ends of oligomers. The physical-chemical properties of the proposed structures were estimated using quantitative structure-property relationships (QSPRs) and quantitative property-property relationships (QPPRs). The anhydro-sugars derived from cellulose are often blamed for coke formation during bio-oil hydrotreatment. Understanding their chemical structure could help to develop rational strategies to mitigate coke formation. Furthermore, the thermo-physical properties reported (boiling point, melting point, Gibb’s free energy of formation, enthalpy of formation, and solubility parameters among others) are also fundamental to conducting first principle engineering calculations to design and analyze new pyrolysis reactors and bio-oil up-grading units.

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

Prospects for Neutron Reactions on Excited States in High-Density Plasmas

With the reactions of high flux neutrons, such as in a DT plasma, there is a prospect of seeing new kinds of neutron-nucleus reactions for the first time. If neutrons excite a heavy nucleus, for example, there is possibility of a second neutron reacting on excited states of the residual nucleus before that nucleus has de-excited to its ground state. The possibility of such reactions on excited states has rarely been considered. The cross section for neutron induced fission on the isomeric state of 235 U has been measured (D’Eer et al., Phys. Rev. C, 1988, 38: 1270–1276) and calculated (Younes et al., 2003, Maslov, 2007), and reactions on rotationally excited nuclear states has been calculated (Kawano et al., Phys. Rev. C, 2009, 80: 024611). In high flux plasmas, however, a much wider range of reactions is possible. We therefore need to consider excited states at much higher-energies than previously modeled, and then estimate whether second neutrons are likely to rescatter on those excited states. To determine the likelihood of such rescattering events, we first need to know the probable time series of nuclear decays of those excited states. The lifetimes of many low-lying states have been measured experimentally, but now we need to know the lifetimes of the many higher excited states that could be produced from incident 14 MeV neutrons. These are too numerous to be measured and also too numerous to be calculated individually, so statistical Hauser-Feshbach decay models are used. I show some lifetime calculations for 89 Y, 169 Tm, and 197 Au targets, and predictions for the number of rescattering events in plausible plasma scenarios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Study of $\langle {p}_{\text{T}}\rangle$ and its higher moments, and extraction of the speed of sound in Pb-Pb collisions with ALICE

Ultrarelativistic heavy-ion collisions produce a state of hot and dense strongly interacting QCD matter called quark-gluon plasma (QGP). On an event-by-event basis, the volume of the QGP in ultracentral collisions is mostly constant, while its total entropy can vary significantly due to quantum fluctuations, leading to variations in the temperature of the system. Exploiting this unique feature of ultracentral collisions allows for the interpretation of the correlation of the mean transverse momentum ($\langle$p T $\rangle$) of produced charged hadrons and the number of charged hadrons as a measure for the speed of sound, c s . This speed is related to the rate at which compression waves travel in the QGP and is determined by fitting the relative increase in $\langle$p T $\rangle$ with respect to the relative change in the average charged-particle density ($\langle$dN ch /dη$\rangle$) measured at mid-rapidity. This study reports the event-average $\langle$p T $\rangle$ of charged particles as well as the variance, skewness, and kurtosis of the event-by-event transverse momentum per charged particle ([p T ]) distribution in ultracentral Pb-Pb collisions at a center-of-mass energy of 5.02 TeV per nucleon pair using the ALICE detector. Different centrality estimators based on charged-particle multiplicity or the transverse energy of the event are used to select ultracentral collisions. By ensuring a pseudorapidity gap between the region used to define the centrality and the region used to perform the measurement, the influence of biases and their potential effects on the rise of the mean transverse momentum is tested. The measured c$^{2}_{s}$ is found to strongly depend on the exploited centrality estimator and ranges between 0.1146±0.0028 (stat.)±0.0065 (syst.) and 0.4374±0.0006 (stat.)±0.0184 (syst.) in natural units. The self-normalized variance shows a steep decrease towards ultracentral collisions, while the self-normalized skewness variables show a maximum, followed by a fast decrease. These non-Gaussian features are understood in terms of the vanishing of the impact-parameter fluctuations contributing to the event-to-event [p T ] distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Voltammetric Measurement of Rates and Energetics for Surface Methoxylation of Si(100) in Methanol with Dissolved Electron Acceptors Using Si Ultramicroelectrodes

The steady-state voltammetric responses of n-type Si(100) semiconductor ultramicroelectrodes (SUMEs) immersed in air- and water-free methanolic electrolytes have been measured. The response characteristics of these SUMEs in the absence of illumination were modeled and understood through a framework that describes the distribution of the applied potential across the semiconductor/electrolyte contact using four discrete regions: the semiconductor space charge, surface, Helmholtz, and diffuse layers. The latter region was described by the full Gouy–Chapman model. This framework afforded insight on how relevant parameters such as the semiconductor band edge potentials, the reorganization energies for charge transfer, the standard potential of redox species in solution, the density and energy of surface state populations, and the presence of an insulating (tunneling) layer individually and collectively dictate the observable current–potential responses. With this information, the methoxylation of Si surfaces was evaluated by analysis of the change in voltammetric responses during the course of prolonged immersion in methanol. Here, the electrochemical data were consistent with a surface methoxylation mechanism that depended on the standard potential of redox species dissolved in solution. Estimates of the enthalpies of adsorption as well as the potential-dependent rate constant for surface methoxylation were obtained. Collectively, these measurements supported the contention that the rates of Si surface reactions can be systematically tuned by exposure to dissolved outer-sphere electron acceptors. Moreover, the data represent the quantitative utility of voltammetry with SUMEs for the measurement of semiconductor/liquid contacts.

14 SOLAR ENERGY↗
DOE OSTI2026

Towards a theoretical understanding of excitonic properties of phthalocyanine thin films. I. Low-temperature exciton absorption spectra

Phthalocyanine (Pc)-based molecular thin films have emerged in recent years as a promising class of organic semiconductor materials for optoelectronic applications owing to their long exciton coherence length and fast exciton diffusion. However, the dependence of their exciton properties on the dimensionality and thermodynamic conditions, presence of metal ions, and effects of chemical modifications to the PC systems is not yet fully understood. As a first step towards a more comprehensive theoretical understanding of the excitonic properties of Pc thin films, we model their low-temperature exciton absorption spectra by employing the Frenkel Hamiltonian. The latter is derived from quantum-chemical estimates of site energies and exciton–exciton couplings. The predicted exciton absorption spectra of octabutoxy phthalocyanine (H 2 OBPc) is found to be strongly dependent on the dimensionality of the model as well as the distance cutoff for including the monomer–monomer exciton coupling. We also caution that the widely used dipole–dipole approximation could substantially overestimate the excitonic coupling between different monomers compared to a more accurate evaluation using the respective transition densities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique

Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.

fluctuations correlations↗
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗
DOE OSTI2025

Data-Driven Closures and Assimilation for Stiff Multiscale Random Dynamics

Here, we introduce a data-driven and physics-informed framework for propagating uncertainty in stiff, multiscale random ordinary differential equations (RODEs) driven by correlated (colored) noise. Unlike systems subjected to Gaussian white noise, a deterministic equation for the joint probability density function (PDF) of RODE state variables does not exist in closed form. Moreover, such an equation would require as many phase-space variables as there are states in the RODE system. To alleviate this curse of dimensionality, we instead derive exact, albeit unclosed, reduced-order PDF (RoPDF) equations for low-dimensional observables/quantities of interest. The unclosed terms take the form of state-dependent conditional expectations, which are directly estimated from data at sparse observation times. However, for systems exhibiting stiff, multiscale dynamics, data sparsity introduces regression discrepancies that compound during RoPDF evolution. This is overcome by introducing a kinetic-like defect term to the RoPDF equation, which is learned by assimilating in sparse, low-fidelity RoPDF estimates. Two assimilation methods are considered, namely nudging and deep neural networks, which are successfully tested against Monte Carlo simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Radionuclide Screening Analysis and Transport Parameters for Pahute Mesa Detonations, Nevada National Security Site

A “screening analysis” is implemented at Pahute Mesa (PM), Nevada National Security Site (NNSS), to determine subsets of the 43 radionuclides listed in the radionuclide inventory of Finnegan et al. (2016) that are of relevance or potential relevance to the hydrologic source term (HST) for assessment of radionuclide transport in groundwater. Consideration is added to how levels of contamination are defined and whether any other radionuclides not in the inventory are relevant to the HST. A model is developed to estimate a range of possible source concentrations in groundwater that account for uncertainty in partitioning into melt glass and sorption into surrounding rock of the exchange zone for the 82 PM underground nuclear tests detonated in vertical shafts. The model also accounts for the varied hydrogeochemical settings. Transport parameters needed for the screening model calculations are developed from databases for hydrogeologic units, chemistry, mineralogy, fracture spacing, fracture aperture, fracture openness, matrix porosity, bulk density, saturation, and alteration. Of key importance is consideration of diagenetic zonation of the mineralogy related to hydrothermal alteration. The screening analysis compares model results with available groundwater radiochemistry data. Consistency between screening model results and available data helps validate the model for application to all 43 radionuclides, most of which have no measurements or only non-detect measurements of concentration in groundwater. The screening analysis determines that ten radionuclides – tritium, Sr-90, I-129, Cs-137, U-232, U-233, U-234, U-238, and the total of Pu-239 and Pu-240 are relevant to the HST. Determination of relevance is primarily based on data and/or model results indicating source concentration exceeding a maximum contaminant level (MCL) out to 100 years since the time of detonation. U-238 is relevant to assessment of contamination relative to MCLs for uranium and gross alpha particle activity. Five other radionuclides – C-14, Cl-36, Tc-99, Np-237, and Pu-238 – are determined potentially relevant to the HST based on a 0.1 MCL threshold. The screening analysis recommends additional attention to the natural daughter radionuclides of the uranium series (U-238) and thorium series (Th-232) decay chains, particularly Ra-226 and Ra-228, that were not included the inventory.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Improvement of Quantification in Non-Destructive Characterization of Radioactive Waste Packages -Introduction to the Ideas of a Currently Started Project - 20197

For the characterization of radioactive waste packages (e.g. 200-L drums) with inhomogeneous activity and density distributions preferably non-destructive methods like segmented gamma scanning and transmission measurements are applied. While the first is used for identification of the radioactive inventory, both are necessary for a reliable activity quantification. Although at nuclear facilities, segmented gamma scanner are state-of-the-art and are mostly equipped with a vertically moveable transmission source the combining of this data and further a-priori information improving accuracy in data estimation giving realistic, reliable and traceable uncertainty estimations are improvable. This challenge is tackled by a R and D project funded by the German Federal Ministry of Education and Research (BMBF). The boundary conditions for this project are (i) taking into account the typical instrumentation available at nuclear facilities, (ii) minimize the overall measurement time and (iii) combining data from emission and transmission measurements by Bayesian methods. For the R and D project, being still in its initial stage, the basic ideas, the data used for evaluation and the result of a preliminary investigation demonstrating the proof of principle are presented. Additionally, some first aspects on spatial resolution in emission and transmission measurements are addressed. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Radiative transition probabilities of neutral and singly ionized Europium estimated by laser-induced breakdown spectroscopy (LIBS)

Laser-induced breakdown spectroscopy (LIBS) is a versatile technique for compositional analysis for solids, liquids, or gasses. LIBS is an asset for the quantitative or qualitative analysis of resource limited materials like actinides and rare earths because it is quasi-nondestructive. Two Eu 2 O 3 pellets were synthesized to be a test and validation set, respectively. Spectral lines identified from the National Institute of Standards and Technology database with fundamental data reported were used to form Saha-Boltzmann plots. The Saha-Boltzmann plots were used to determine the plasma temperatures and electron densities of the laser-induced plasmas for both samples. These Saha-Boltzmann plots were then used to calculate previously unreported transition probabilities associated with identified peaks. Additionally, the transition probabilities presented in this paper provide the capability for calibration free LIBS to be performed more readily on europium samples and the spectroscopic analysis of stellar bodies. Eight previously unreported transition probabilities are presented in this paper; five for Eu I and three for Eu II. The transition probabilities for Eu I ranged from 0.172 to 7.38 × 10 7 s -1 and those for Eu II ranged from 1.56 to 6.75 × 10 7 s -1 .

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Science–based maize stover removal can be sustainable

Maize (Zea mays L.) stover can be harvested for multiple uses or left in the field to sustain soil organic carbon (SOC), cycle essential plant nutrients, and protect soil health. This 13-yr field study quantified effects of no (0 Mg ha –1 yr –1 ), low (1.0–1.4 Mg ha –1 yr –1 ), moderate (3.5–4.0 Mg ha –1 yr –1 ), or high rates (4.7–5.4 Mg ha –1 yr –1 ) of stover harvest from either continuous maize or maize–soybean [Glycine max (L.) Merr.] rotation on grain yield, plant nutrient concentrations, and multiple soil properties at two sites in Iowa. Stover harvest increased plant macro- and micro-nutrient removal, but did not affect average grain yields of either crops. Soil inorganic carbon (IC), SOC, bulk density, pH, and cation exchange capacity (CEC) showed no significant differences due to stover harvest. Plant tissue and soil-test nutrient concentration effects were also minor and site-specific. Stover harvest significantly (p < .05) decreased exchangeable K and Ca concentrations by 8.3–23.8% and 0.3–22.5% but overall soil health indicator effects were minimal. Overall, based on crop yields, plant nutrient and soil-test concentrations, soil health indicators, and carbon sequestration estimates, maize stover harvest can be sustainable provided: (a) grain yields consistently exceed 11 Mg ha –1 , (b) stover removal does not exceed 40% of the aboveground biomass (i.e., 3.5–4.0 Mg ha –1 yr –1 ), and (c) plant nutrients (especially K) are closely monitored.

54 ENVIRONMENTAL SCIENCES↗