Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “k mean”

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 811 records · Page 45

Biogeochemical consequences of microbial evolution under drought (Final technical report)

As the climate changes, droughts may become more severe. Scientists are uncertain about how drought will affect the natural world, particularly the bacteria, fungi, and other microbes that live in soils and control the Earth’s flows of carbon and nutrients. Researchers at the University of California, Irvine, and Lawrence Berkeley National Laboratory teamed up to study how microbes at the soil surface change with drought and what that means for the carbon cycle. Since 2007, the researchers have used shelters with retractable roofs to prevent nearly half of normal rainfall from reaching grass and shrub landscapes, and their soil microbiomes, in Southern California.

54 ENVIRONMENTAL SCIENCES↗

I Know I'm Right, But Does My Phone?

Transportation is the largest source of green-house gas emissions in the United States. Reducing transportation emissions depends on human travel behavior, which relies on local land use and planning. Travel diaries, consisting of sequences of trips between places for a particular individual, are typically used to instrument human travel behavior. However, these diaries are only as accurate as the underlying methods used to construct them. Travel diary algorithms have been a popular research topic since the advent of GPS tracking surveys. Mode inference algorithms in particular have been well represented in literature. However, these algorithms have typically been validated using prompted recall of pre-segmented trips, which doesn't account for segmentation error, thus disregarding the continuity of mode inference. Furthermore, phone operating systems and applications have adopted battery-conserving techniques, but we are not aware of prior work that has characterized the resulting data collection errors or evaluated procedures to mitigate them. We introduce a framework to evaluate accuracy of trip length computations and mode inference. We develop a temporal alignment procedure in analyzing continuous mode-segmented trajectories for groups of trips. We then apply our framework to evaluate an example set of travel diary algorithms from the open-source OpenPATH travel diary platform against MobilityNet, a public dataset containing information from three artificial timelines that cover 15 different travel modes. Our results show that inference based on an integration with map features results in weighted F_1 scores of 0.60 (iOS) and 0.74 (android). We also show that OpenPATH tends to under count trip length, with mean of signed relative error of -0.0438 on android and -0.0704 on iOS. We hope that other travel diary algorithms will be evaluated using this standardized process, and that the results used to understand and improve the state-of-the-art in this field.

ADVANCED PROPULSION SYSTEMS,ENERGY PLANNING, POLIC↗

Deformation behavior of annealed Cu64Zr36 metallic glass via molecular dynamics simulations

In the present work, we investigate the isothermal annealing process of Cu 64 Zr 36 metallic glass (MG) by means of molecular dynamics (MD) simulations, and characterize the compression properties under different temperatures and compressive rates. The pair distribution functions, the bond-orientational order, and the cluster-type index method are modeled and analyzed to characterize changes in the structure. Results of the modeling and analysis reveal that the MgZn 2 -type Laves, MgCu 2 -type Laves and five-fold local symmetry structures can be formed during the isothermal holding at 950 K. Study on the compression properties shows that at a fixed strain rate, the yield strength and yield drop all increase with decreasing temperature. At a fixed temperature, however, the increase of the strain rate leads to a noticeable increase in the yield strength at 600 K, but has little effect on yield strength at 10 K. Moreover, the modeling and analysis of the structure at a temperature of 10 K and strain-rate of 5 × 10 7 s -1 demonstrate that the order-disorder transformation initiates the shear band.

36 MATERIALS SCIENCE↗

An Ensemble of Neural Networks for Moist Physics Processes, Its Generalizability and Stable Integration

Abstract With the recent advances in data science, machine learning has been increasingly applied to convection and cloud parameterizations in global climate models (GCMs). This study extends the work of Han et al. (2020, https://doi.org/10.1029/2020MS002076 ) and uses an ensemble of 32‐layer deep convolutional residual neural networks, referred to as ResCu‐en, to emulate convection and cloud processes simulated by a superparameterized GCM, SPCAM. ResCu‐en predicts GCM grid‐scale temperature and moisture tendencies, and cloud liquid and ice water contents from moist physics processes. The surface rainfall is derived from the column‐integrated moisture tendency. The prediction uncertainty inherent in deep learning algorithms in emulating the moist physics is reduced by ensemble averaging. Results in 1‐year independent offline validation show that ResCu‐en has high prediction accuracy for all output variables, both in the current climate and in a warmer climate with +4K sea surface temperature. The analysis of different neural net configurations shows that the success to generalize in a warmer climate is attributed to convective memory and the 1‐dimensional convolution layers incorporated into ResCu‐en. We further implement a member of ResCu‐en into CAM5 with real world geography and run the neural‐network‐enabled CAM5 (NCAM) for 5 years without encountering any numerical integration instability. The simulation generally captures the global distribution of the mean precipitation, with a better simulation of precipitation intensity and diurnal cycle. However, there are large biases in temperature and moisture in high latitudes. These results highlight the importance of convective memory and demonstrate the potential for machine learning to enhance climate modeling.

Meteorology & Atmospheric Sciences↗

Direct observation of coherence transfer and rotational-to-vibrational energy exchange in optically centrifuged CO2 super-rotors

Abstract Optical centrifuges are laser-based molecular traps that can rotationally accelerate molecules to energies rivalling or exceeding molecular bond energies. Here we report time and frequency-resolved ultrafast coherent Raman measurements of optically centrifuged CO 2 at 380 Torr spun to energies beyond its bond dissociation energy of 5.5 eV ( J max = 364, E rot = 6.14 eV, E rot / k B = 71, 200 K). The entire rotational ladder from J = 24 to J = 364 was resolved simultaneously which enabled a more accurate measurement of the centrifugal distortion constants for CO 2 . Remarkably, coherence transfer was directly observed, and time-resolved, during the field-free relaxation of the trap as rotational energy flowed into bending-mode vibrational excitation. Vibrationally excited CO 2 ( ν 2 > 3) was observed in the time-resolved spectra to populate after 3 mean collision times as a result of rotational-to-vibrational (R-V) energy transfer. Trajectory simulations show an optimal range of J for R-V energy transfer. Dephasing rates for molecules rotating up to 5.5 times during one collision were quantified. Very slow decays of the vibrational hot band rotational coherences suggest that they are sustained by coherence transfer and line mixing.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Genome-Based Model to Predict the Virulence of Pseudomonas aeruginosa Isolates

ABSTRACT: Variation in the genome of Pseudomonas aeruginosa , an important pathogen, can have dramatic impacts on the bacterium’s ability to cause disease. We therefore asked whether it was possible to predict the virulence of P. aeruginosa isolates based on their genomic content. We applied a machine learning approach to a genetically and phenotypically diverse collection of 115 clinical P. aeruginosa isolates using genomic information and corresponding virulence phenotypes in a mouse model of bacteremia. We defined the accessory genome of these isolates through the presence or absence of accessory genomic elements (AGEs), sequences present in some strains but not others. Machine learning models trained using AGEs were predictive of virulence, with a mean nested cross-validation accuracy of 75% using the random forest algorithm. However, individual AGEs did not have a large influence on the algorithm’s performance, suggesting instead that virulence predictions are derived from a diffuse genomic signature. These results were validated with an independent test set of 25 P. aeruginosa isolates whose virulence was predicted with 72% accuracy. Machine learning models trained using core genome single-nucleotide variants and whole-genome k-mers also predicted virulence. Our findings are a proof of concept for the use of bacterial genomes to predict pathogenicity in P. aeruginosa and highlight the potential of this approach for predicting patient outcomes. IMPORTANCE Pseudomonas aeruginosa is a clinically important Gram-negative opportunistic pathogen. P. aeruginosa shows a large degree of genomic heterogeneity both through variation in sequences found throughout the species (core genome) and through the presence or absence of sequences in different isolates (accessory genome). P. aeruginosa isolates also differ markedly in their ability to cause disease. In this study, we used machine learning to predict the virulence level of P. aeruginosa isolates in a mouse bacteremia model based on genomic content. We show that both the accessory and core genomes are predictive of virulence. This study provides a machine learning framework to investigate relationships between bacterial genomes and complex phenotypes such as virulence.

59 BASIC BIOLOGICAL SCIENCES↗

Search for exotic decays of the Higgs boson into long-lived particles in pp collisions at $\sqrt{s}$ = 13 TeV using displaced vertices in the ATLAS inner detector

A novel search for exotic decays of the Higgs boson into pairs of long-lived neutral particles, each decaying into a bottom quark pair, is performed using 139 fb -1 of $\sqrt{s}$ = 13 TeV proton-proton collision data collected with the ATLAS detector at the LHC. Events consistent with the production of a Higgs boson in association with a leptonically decaying Z boson are analysed. Long-lived particle (LLP) decays are reconstructed from inner-detector tracks as displaced vertices with high mass and track multiplicity relative to Standard Model processes. The analysis selection requires the presence of at least two displaced vertices, effectively suppressing Standard Model backgrounds. The residual background contribution is estimated using a data-driven technique. No excess over Standard Model predictions is observed, and upper limits are set on the branching ratio of the Higgs boson to LLPs. Branching ratios above 10% are excluded at 95% confidence level for LLP mean proper lifetimes cτ as small as 4 mm and as large as 100 mm. For LLP masses below 40 GeV, these results represent the most stringent constraint in this lifetime regime.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design

Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets—from soil and seawater microbiomes—to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon’s abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.

59 BASIC BIOLOGICAL SCIENCES↗

Aspects of Crystal Growth and Atomic-Scale Characterization of U/Th Age-Dating Particles

Trace analysis of particulate material collected on environmental swipe samples continues to be one of the cornerstones of the IAEA’s process of verifying member state compliance as part of the Non-Proliferation Treaty (NPT). While complementary to traditional bulk analysis, the ability to measure individual particles can be far more impactful from a radiochronometry perspective. The 234U/230Th radiochronometer is being evaluated as a means of establishing model ages for discrete particles with sufficient confidence. To support the extension of current capabilities, reference materials with a suitable particle form factor and well-documented purification date are actively sought for development with large geometry secondary ion mass spectroscopy.

Nuclear Safeguards, IAEA, Particle Standards↗

Electronic structure and exchange interactions in altermagnetic MnGeP 2 in the quasiparticle self-consistent G W approach

The quasiparticle self-consistent 𝐺⁡𝑊 method is used to study the electronic band structure, optical dielectric function, and exchange interactions in chalcopyrite, 𝐼⁢ ‾ 4 ⁢2⁢𝑑, structure MnGeP 2 . The material is found to be an antiferromagnetic semiconductor with a lowest direct gap of 2.44 eV at the Γ point and a lower indirect gap of 1.87 eV from Γ to 𝑀. The material is an altermagnet because the two magnetic atoms of opposite spin are related by a twofold rotation operation perpendicular to the main fourfold rotation inversion axis. The spin splittings along a low symmetry line like 𝑃⁢𝑁 is sizable while at k points on the diagonal mirror planes or on the twofold symmetry axes the spin splitting is zero. The exchange interactions are calculated using a linear response approach. The antiferromagnetic exchange interaction between nearest neighbors in the primitive unit cell is dominating and found to be slightly decreasing upon carrier doping but not sufficiently to change the interaction to become ferromagnetic. The bare (noninteracting) and interacting transverse spin susceptibilities, which provide interatomic site exchange interactions after averaging over the muffin-tin spheres, are calculated from the 𝐺⁡𝑊 band structure and wave functions. From these exchange interactions, the spin wave spectra are obtained along the high symmetry lines and the Néel temperature is calculated using the mean-field and Tyablikov estimations. The dielectric function and the optical absorption spectra are calculated including excitonic effects using the Bethe Salpeter equation. The exchange interactions around Mn Ge defect sites is also studied. While we find it can generate ferromagnetic interactions with neighboring spins, we did not find direct evidence of producing an overall ferromagnetic phase. First, if Mn antisites are introduced by exchanging Mn with a nearby Ge, the interactions stay largely antiferromagnetic. Second, when we add additional Mn, in other words in Mn-rich stoichiometry, the Mn Ge antisites produce a strong ferromagnetic interaction primarily with the Mn in the same basal plane but weaker ferromagnetic interaction with adjacent plane Mn. The interactions between regular lattice Mn stay antiferromagnetic as before and thus favor keeping the antiferromagnetic order along the [001] direction. Adding Mn antisites, however, does lead to a metallic band structure.

36 MATERIALS SCIENCE↗

Interface roughness governed negative magnetoresistances in two-dimensional electron gases in AlGaN/GaN heterostructures

Negative magnetoresistances (NMRs) have been widely observed in two-dimensional electron gases (2DEGs). However, their origins are under debate. Here, we report on NMRs in the 2DEG in AlGaN/GaN heterostructures, aiming to uncover their origins by utilizing electric field gating. We systematically measured the magnetoresistances in magnetic fields up to 12 T and at temperatures between 3 and 260 K and observed NMRs over a wide range of temperatures from 3 to 170 K, which become more pronounced with decreasing temperature. We conducted electric field gating experiments to correlate the occurrence of NMRs with the relationship between electron mobility and density. The latter is governed by defect scattering sources in the sample and can be theoretically modeled. Comparison of the measured electron mobility and electron density relationship with theory reveals that interface roughness scattering plays a crucial role in obtaining large NMRs. Our work demonstrates that electric field gating provides a means not only to tune the values of NMRs but also to uncover their mechanisms.

2-dimensional systems↗

Multimodal Approach Reveals the Symmetry-Breaking Pathway to the Broken Helix in EuIn 2 ⁢As 2

Understanding and manipulating emergent phases, which are themes at the forefront of quantum-materials research, rely on identifying their underlying symmetries. This general principle has been particularly prominent in materials with coupled electronic and magnetic degrees of freedom, in which magnetic order influences the electronic band structure and can lead to exotic topological effects. However, identifying symmetry of a magnetically ordered phase can pose a challenge, particularly in the presence of small domains. Here we introduce a multimodal approach for determining magnetic structures, which combines symmetry-sensitive optical probes, scattering, and group-theoretical analysis. We apply it to EuIn 2 ⁢As 2 , a material that has received attention as a candidate axion insulator. While first-principles calculations predict this state on the assumption of a simple collinear antiferromagnetic structure, subsequent neutron-scattering measurements reveal a much more intricate magnetic ground state characterized by two coexisting magnetic wave vectors reached by successive thermal phase transitions. The proposed high- and low-temperature phases are a spin helix and a state with interpenetrating helical and Néel antiferromagnetic order termed a “broken helix,” respectively. Employing a multimodal approach, we identify the magnetic structure associated with these two phases of EuIn 2 ⁢As 2 . We find that the higher-temperature phase is characterized by a variation of the magnetic moment amplitude from layer to layer, with the moment vanishing entirely in every third Eu layer. The lower-temperature structure is similar to the broken helix, with one important difference: Because of local strain, the relative orientation of the magnetic structure and the lattice is not fixed. Consequently, the symmetry required to protect the axion phase is not generically protected in EuIn 2 ⁢As 2 , but we show that it can be restored if the magnetic structure is tuned with uniaxial strain. Finally, we present a spin Hamiltonian that identifies the spin interactions that account for the complex magnetic order in EuIn 2 ⁢As 2 . Our work highlights the importance of a multimodal approach in determining the symmetry of complex order parameters.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Further Evidence of Modified Spin-down in Sun-like Stars: Pileups in the Temperature–Period Distribution

We combine stellar surface rotation periods determined from NASA's Kepler mission with spectroscopic temperatures to demonstrate the existence of pileups at the long-period and short-period edges of the temperature–period distribution for main-sequence stars with temperatures exceeding ~5500 K. The long-period pileup is well described by a curve of constant Rossby number, with a critical value of Ro crit ≲ Ro ⊙ . The long-period pileup was predicted by van Saders et al. as a consequence of weakened magnetic braking, in which wind-driven angular momentum losses cease once stars reach a critical Rossby number. Stars in the long-period pileup are found to have a wide range of ages (~2–6 Gyr), meaning that, along the pileup, rotation period is strongly predictive of a star's surface temperature but weakly predictive of its age. The short-period pileup, which is also well described by a curve of constant Rossby number, is not a prediction of the weakened magnetic braking hypothesis but may instead be related to a phase of slowed surface spin-down due to core-envelope coupling. The same mechanism was proposed by Curtis et al. to explain the overlapping rotation sequences of low-mass members of differently aged open clusters. The relative dearth of stars with intermediate rotation periods between the short- and long-period pileups is also well described by a curve of constant Rossby number, which aligns with the period gap initially discovered by McQuillan et al. in M-type stars. These observations provide further support for the hypothesis that the period gap is due to stellar astrophysics, rather than a nonuniform star formation history in the Kepler field.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Effects of Aliovalent Substitutions in Thermoelectric Zintl Pnictides. Varied Polyanionic Dimensionality and Complex Structural Transformations–The Case of Sr 3 ZnP 3 vs Sr 3 Al x Zn 1– x P 3

The structures and the transport properties of a novel family of Zintl phosphides and arsenides with the formula AE 3 ZnPn 3 and the solid solutions AE 3 Al x Zn 1–x Pn 3 , AE 3 ZnAs y P 1–y (AE = Sr, Eu; Pn = P, As) are reported. Crystals of nine new phases have been obtained via Pb-flux reactions and used for structural work by means of single-crystal X-ray diffraction methods. The derived orthorhombic structure is without a direct analog, and features unusual structural units, where the Zn atoms are in both distorted tetrahedral and trigonal-planar coordination of pnictogens. Electronic structure calculations reveal moderately wide bandgaps for Sr 3 ZnP 3 and Sr 3 ZnAs 3 , on the order of 0.70 and 0.63 eV, respectively. Electrical transport measurements above room temperature indicate relatively high resistivity values above 500 K (ρ ≈ 4.8 Ω cm and above), but some of the samples exhibit very high Seebeck coefficients, as large as 300 μV/K at 560 K for Sr 3 ZnAs 3 . Aliovalent substitutions in AE 3 ZnPn 3 , achieved by the partial replacement of Zn 2+ with Al 3+ cations promote occupational and positional disorder, which causes structural transformation towards the disordered variant of the Sr 5 Al 2 Sb 6 structure type. Such substitutions also change the dimensionality of the polyanionic sub-lattice in the resulting quaternary AE 3 Al x Zn 1–x Pn 3 phases. Furthermore, preliminary transport property data on the latter reveal nine times lower electrical resistivity (ρ 500 ≈ 0.5 Ω cm) together with a significantly enhanced Seebeck coefficient, αmax ≈ 430 μV/K at 560 K.

36 MATERIALS SCIENCE↗

Muted extratropical low cloud seasonal cycle is closely linked to underestimated climate sensitivity in models

Abstract A large spread in model estimates of the equilibrium climate sensitivity (ECS), defined as the global mean near-surface air-temperature increase following a doubling of atmospheric CO 2 concentration, leaves us greatly disadvantaged in guiding policy-making for climate change adaptation and mitigation. In this study, we show that the projected ECS in the latest generation of climate models is highly related to seasonal variations of extratropical low-cloud fraction (LCF) in historical simulations. Marked reduction of extratropical LCF from winter to summer is found in models with ECS > 4.75 K, in accordance with the significant reduction of extratropical LCF under a warming climate in these models. In contrast, a pronounced seasonal cycle of extratropical LCF, as supported by satellite observations, is largely absent in models with ECS < 3.3 K. The distinct seasonality in extratropical LCF in climate models is ascribed to their different prevailing cloud regimes governing the extratropical LCF variability.

54 ENVIRONMENTAL SCIENCES↗

A strong blend in the morning: studying the circumgalactic medium before cosmic noon with strong, blended Ly α forest systems

ABSTRACT We study of the properties of a new class of circumgalactic medium absorbers identified in the Ly α forest: ‘Strong, Blended Lyman-α’ (or SBLA) absorption systems. We study SBLAs at 2.4 < z < 3.1 in SDSS-IV/eBOSS spectra by their strong extended Ly α absorption complexes covering 138 $\, \, {\rm km}\, {\rm s}^{-1}$ with an integrated $\log (N_{\rm H\, {\small I}}/\mathrm{cm}^{-2}) =16.04$$\substack{+0.05 \\ -0.06}$ and Doppler parameter b = 18.1$\substack{+0.7 \\ -0.4}$$\, \, {\rm km}\, {\rm s}^{-1}$. Clustering with the Ly α forest provides a large-scale structure bias of b = 2.34 ± 0.06 and halo mass estimate of $M_h \approx 10^{12}\, h^{-1}\, {\rm M_{\odot }}$ for our SBLA sample. We measure the ensemble mean column densities of 22 metal features in the SBLA composite spectrum and find that no single-population multiphase model for them is viable. We therefore explore the underlying SBLA population by forward modelling the SBLA absorption distribution. Based on covariance measurements and favoured populations we find that ≈25 per cent of our SBLAs have stronger metals. Using silicon only we find that our strong metal SBLAs trace gas with a log (nH/cm−3) > −2.40 for T = 103.5 K and show gas clumping on <210 parsec scales. We fit multiphase models to this strong subpopulation and find a low ionization phase with nH = 1 cm−3, T = 103.5 K, and [X/H] = 0.8, an intermediate ionization phase with log (nH/cm−3) = −3.05, T = 103.5 K and [X/H] = −0.8, and a poorly constrained higher ionization phase. We find that the low ionization phase favours cold, dense super-solar metallicity gas with a clumping scale of just 0.009 parsecs.

Morrison, Sean (ORCID:0000000267702627)↗

Improving Cathode Testing with a High-Gradient Cryogenic Normal Conducting RF Photogun

Future electron accelerator applications such as X-ray free electron lasers and colliders are dependent on significantly increasing beam brightness. With the observation that linac beam manipulation’s best preservation of max brightness is at the cathode, we are incentivized to create an environment where we can study how to achieve the highest possible photogun brightness. In order to do so, we intend to extract beams from high-brightness photocathodes with the highest achievable accelerating gradients we can manage in a klystron-powered radiofrequency (RF) photogun. We utilize here cryogenic normal conducting cavities to achieve ultra-high gradients via limitation of breakdown rates (BDR). The low temperatures should also reduce cathode emittance by reducing the mean transverse energy (MTE) of electrons near the photoemission threshold. To this end, we have designed and produced a new CrYogenic Brightness-Optimized Radiofrequency Gun (CYBORG) for use in a new beamline at UCLA. We will introduce the enabling RF and photoemission physics as a primer for the new regime of high field low temperature cathodes we intend to enter. We further report the current status of the beamline commissioning, including the cooling of the photogun to 100 K, and producing 0.5 MW of RF feed power, which corresponds to cathode accelerating fields in the range of 80–90 MV/m. We further plan iterative improvements to both to 77 K and 1 MW corresponding to our ultimate goal >120 MV/m. Our discussion will include future beamline tests and the consideration of the initial realization of an ultra-high-gradient photoinjector concept.

Lawler, Gerard Emile (ORCID:0000000160113435)↗

Phosphorus budgets of intensively managed row crops at a long–term agroecosystem research site in the upper US Midwest

Phosphorus (P) budgets for cropping systems provide insights for keeping soil P at optimal levels for crops while avoiding excess inputs. We quantified 12 years of P inputs (fertilizer and atmospheric deposition) and outputs (harvest and leaching losses) for replicated maize (Zea mays L.)—soybean (Glycine max L.)—wheat (Triticum aestivum) crop rotations under conventional, no-till, reduced input, and biologically based (organic without compost or manure) management systems at the Kellogg Biological Station LTAR site in southwest Michigan. Conventional, no-till, and reduced input systems were fertilized between 13 and 50 kg P ha –1 depending on year. Soil test phosphorus (STP) was measured at 0- to 25-cm depth every autumn. Leached P was measured as dissolved P in the soil solution sampled beneath the rooting depth and combined with modeled percolation. Fertilization and harvest were the predominant P fluxes in the fertilized systems, whereas only harvest dominated P flux in the unfertilized organic system. Leaching losses were minor terms in the budgets, but leachate concentrations were nevertheless close to the range of concern for downstream eutrophication. Over the 12-year study period, the organic system exhibited a negative P balance (–82.0 kg P ha –1 ), coinciding with suboptimal STP levels, suggesting a need for P supplementation. In contrast, the fertilized systems showed positive P balances (mean: 70.1 kg P ha –1 ) with STP levels well above agronomic optima. Results underscore the importance of tailored P management strategies to sustain crop productivity while mitigating environmental impacts.

60 APPLIED LIFE SCIENCES↗