Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “peak finding”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Probing quarkyonic matter in neutron stars with the Bayesian nuclear-physics multimessenger astrophysics framework

The interiors of neutron stars contain matter at the highest densities realized in our Universe. Interestingly, theoretical studies of dense matter, in combination with the existence of two-solar-mass neutron stars, indicate that the speed of sound $c_s$ has to increase to values well above the conformal limit ($c_s^2$ = 1/3) before decreasing again at higher densities. Further, the decrease could be explained by either a strong first-order phase transition or a crossover transition from hadronic to quark matter. The latter scenario leads to a pronounced peak in the speed of sound, reaching values above the conformal limit, naturally explaining the inferred behavior. In this work, we use the nuclear-physics multimessenger astrophysics (NMMA) framework to compare predictions of the quarkyonic matter model with astrophysical observations of neutron stars, with the goal of constraining model parameters. Assuming quarkyonic matter to be realized within neutron stars, we find that there can be a significant amount of quarks inside the cores of neutron stars with masses in the two-solar-mass range, amounting to up to ≈0.13$M$ ⊙ , contributing ≈ 5.9% of the total mass. Furthermore, for the quarkyonic matter model investigated here, the radius of a 1.4$M$ ⊙ neutron star would be $13.44_{–1.54}^{+1.69}(13. 54_{–1.04}^{+1.02})$ km, at 95% credibility, without (with) the inclusion of AT2017gfo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Fast Charging Infrastructure for Electrifying Road Trips to and from National Parks in the Western United States

National parks in the Western United States draw over 80 million visitors every year, and most visitors rely on personal cars for their road trips (or long-distance travels). Travel to national parks represents distinct travel demand, as they are typically located in remote areas necessitating long-distance trips. This study investigates the quantity and locations of on-route fast charging infrastructure needed by 2030 to enable seamless travel to/from national parks using electric vehicles in seven target states in the region, employing unprecedented high-resolution spatial and temporal analysis. We find that the required number of fast charging ports for on-route charging infrastructure ranges from 1,200 to 22,000, depending on different assumptions of key input parameters - vehicle electrification rate, charging behavior, average gap between charging stations, port utilization rate, and towing trailers. Our analysis also indicates that electrical load for on-route fast charging infrastructure would peak in the afternoon, in the range of 70-400 MW, varying with the key input parameters. This study illustrates how different input parameters result in different degrees of impact on various aspects of charging infrastructure. We also examine the characteristics of projected charging infrastructure in terms of land use type, relationship with traffic volume, size of stations, and other variables.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Causal relationships of vegetation productivity with root zone water availability and atmospheric dryness at the catchment scale

Abstract. This study explores the causal relationships between catchment water availability, vapor pressure deficit, and gross primary productivity (GPP) across 341 catchments in the contiguous US. Seasonal climatic, hydrological, and vegetation characteristics were represented using the Horton index, ecological aridity index, evaporative fraction index, and carbon uptake efficiency. Statistical methods, including circularity statistics, correlation analysis, and causality tests, were employed to determine the complex interactions between catchment wetness, atmospheric dryness, and vegetation carbon uptake. The results revealed a maximum lag of 2 months in the intra-annual variability of catchment water supply–productivity and atmospheric water demand–productivity relationships, with hysteresis patterns varying with the catchment's hydrological characteristics. In catchments not permanently under water-limited or energy-limited conditions, vegetation experiences hydrological stress during the peak growing period, coinciding with the highest gross primary productivity and carbon uptake efficiency being out of phase with the Horton index and in phase with the evaporative fraction index. Causality analysis highlights strong temporal continuity in GPP seasonal characteristics, with a cause–effect relationship between catchment water supply, atmospheric demand, and vegetation productivity spanning a maximum of 2 months. These findings underscore the need for a comprehensive functional framework that integrates catchment water supply, atmospheric demand, and vegetation productivity to enhance our understanding and predictive capabilities with regard to ecosystem responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Measurement of the Triple-Differential Cross Section of Muon Neutrino Charged Current Interactions with Low Hadronic Activity in the NO$\nu$A Near Detector

This dissertation describes the analysis and measurement of the triple-differential cross-section of the $\nu_\mu CC$ interaction with Low Hadronic Activity (maximum kinetic energy of 250 MeV for protons and 175 MeV for pions) in the NOvA Near Detector. This cross-section is reported in bins of Muon Kinematics and available hadronic energy. The data corresponds to an accumulated $14.2\times10^{20}$ protons-on-target (POT) in the neutrino mode of the NuMI beam, with a narrow band of neutrino energies peaked at 1.8 GeV. The analysis provides a sample of neutrino–nucleus interactions with an enhanced fraction of quasi-elastic (QE) and two-particle-two-hole (2p2h) interactions and a suppression of resonant pion production (RES) and deep inelastic scattering (DIS) interactions. This enhancement of QE and 2p2h allows for quantitative comparisons with various models of these processes. We find strong disagreement between data and theory-based models in the forward muon direction, especially for the 0.2-0.4 GeV available hadronic energy region.

Lesmeister, James Ryan [Houston U.]↗
DOE OSTI2020

Characterizing NWP Model Errors Using Doppler-Lidar Measurements of Recurrent Regional Diurnal Flows: Marine-Air Intrusions into the Columbia-River Basin

Ground-based Doppler-lidar instrumentation provides atmospheric wind data at dramatically improved accuracies and spatial/temporal resolutions. These capabilities have provided new insights into atmospheric flow phenomena, but they also should have a strong role in NWP model improvement. Insight into the nature of model errors can be gained by studying recurrent atmospheric flows, here a regional summertime diurnal sea breeze and subsequent marine-air intrusion into the arid interior of Oregon-Washington, where these winds are an important windenergy resource. These marine intrusions were sampled by three scanning Doppler lidars in the Columbia-River Basin as part of the Second Wind Forecast Improvement Project (WFIP2), using data from summer 2016. Lidar time-height cross sections of wind speed identified eight days when the diurnal-flow cycle, (peak wind speeds at midnight, afternoon minima) was obvious and strong. Eight-day composite time-height cross sections of lidar wind speeds are used to validate those generated by the operational NCEP-HRRR model. HRRR simulated the diurnal wind cycle, but produced errors in the timing of onset and significant errors due to a premature nighttime demise of the intrusion flow, producing low-bias errors of 6 m s-1 . Day-to-day and in the composite, whenever a marine intrusion occurred, HRRR made these same errors. The errors occurred under a range of gradient wind conditions indicating that they resulted from the misrepresentation of physical processes within a limited region around the measurement locations. Because of their generation within a limited geographical area, field-measurement programs can be designed to find and address the sources of these NWP errors.

Banta, Robert M.↗
DOE OSTI2022

Temporal complementarity of marine renewables with wind and solar generation: Implications for GB system benefits

Wave and tidal energy have the potential to provide benefits to power systems with high proportions of stochastic renewable generation. This is particularly applicable with wind and solar photovoltaics, as the offsetting of these renewable resources results in more reliable renewable generation. This study utilises ten metrics to quantify the temporal complementarity and supply-demand balancing requirements of the energy mix in Great Britain, to investigate the potential magnitude of these system benefits. Wave and tidal generation profiles are created using historical resource data and hydrodynamic models. The results show that the inclusion of wave and tidal generation creates a renewable energy mix which is more available under multiple conditions: throughout a year of operation; at times of peak demand; for multiple consecutive hourly time periods; and at times when wind and solar generation are not available. Three regional case studies also show that the inclusion of marine energy allows for improved regional supply-demand matching, reducing instances of energy shortage and excess and potentially relieving transmission congestion at particularly constrained locations within GB. Finally, the implications of these findings are discussed in terms of GB wholesale market operation, system balancing and system security.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Understanding Solution State Conformation and Aggregate Structure of Conjugated Polymers via Small Angle X-ray Scattering

Donor-acceptor (D-A) conjugated polymers are high-performance organic electronic materials that exhibit complex aggregation behavior. Understanding the solution state conformation and aggregation of conjugated polymers is crucial for controlling morphology during thin-film deposition and the subsequent electronic performance. However, a precise multiscale structure of solution state aggregates is lacking. Here, we present an in-depth small-angle X-ray scattering (SAXS) analysis of the solution state structure of an isoindigo-bithiophene-based D-A polymer (PII-2T) in chlorobenzene and decane as our primary system. Modeling the system as a combination of hierarchical fibrillar aggregates mixed with dispersed polymers, we extract information about conformation and multiscale aggregation and also clarify the physical origin of features often observed but unaddressed or misinterpreted in small-angle scattering patterns of conjugated polymers. The persistence length of the D-A polymer extracted from SAXS agrees well with a theoretical model based on the dihedral potentials. Additionally, we show that the broad high q structure factor peak seen in scattering profiles can be attributed to lamellar stacking occurring within the fibril aggregates and that the low q aggregate scattering is strongly influenced by the polymer molecular weight. Overall, the SAXS profiles of D-A polymers in general exhibit a sensitive dependence on the co-existence of fibrillar aggregate and dispersed polymer chain populations. We corroborate our findings from SAXS with electron microscopy of freeze-dried samples for direct imaging of fibrillar aggregates. Finally, we demonstrate the generality of our approach by fitting the scattering profiles of a variety of D-A polymers based on thieno-isoindigo (PTII-2T), diketopyrrolopyrrole (DPP2T-TT, DPP-BTZ, PDPP2FT-C- 16 ), naphthalenediimide (P(NDI2OD-T2)), and a conjugated block copolymer P3HT-b-DPPT-T. The results presented here establish a picture of the D-A polymer solution state structure and provide a general method of interpreting and analyzing their scattering profiles.

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

A type-I van der Waals heterostructure formed by monolayer WS 2 and trilayer PdSe 2

Two-dimensional (2D) heterostructures, formed by stacking 2D semiconductors through the van der Waals force, have been extensively studied recently. However, the majority of the heterostructures discovered so far possess type-II interfaces that facilitate interlayer charge separation. Type-I interfaces, on the other hand, confine both electrons and holes in one layer, which is beneficial for optical applications that utilize electron–hole radiative recombination. So far, only a few type-I 2D heterostructures have been achieved, which has limited the construction of multilayer heterostructures with sophisticated band landscapes. Here, we report experimental evidence of a type-I interface between monolayer WS 2 and trilayer PdSe 2 . Two-dimensional PdSe 2 has emerged as a promising material for infrared optoelectronic and other applications. We fabricated the heterostructure by stacking an exfoliated monolayer WS 2 flake on top of a trilayer PdSe 2 film, synthesized by chemical vapor deposition. Photoluminescence spectroscopy measurements revealed that the WS 2 exciton peak is significantly quenched in the heterostructure, confirming efficient excitation transfer from WS 2 to PdSe 2 . Femtosecond transient absorption measurements with various pump/probe configurations showed that both electrons and holes photoexcited in the WS 2 layer of the heterostructure can efficiently transfer to PdSe 2 , while neither type of carriers excited in PdSe 2 can transfer to WS 2 . These experimental findings establish a type-I band alignment between monolayer WS 2 and trilayer PdSe 2 . Furthermore, our results further highlight PdSe 2 as an important 2D material for constructing van der Waals heterostructures with emergent electronic and optoelectronic properties.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

X-ray photoelectron spectroscopy of surfactants on sub-micron aqueous aerosols

Organic molecules at the aqueous interfaces of droplets and aerosols can influence multiphase chemistry and cloud condensation nuclei activity. While surface tension measurements on flat surfaces are commonly used to quantify organic partitioning and structure, extending these measurements to more realistic aerosols that might exist in metastable, supercooled and supersaturated states remains a challenge. Here, we use aerosol velocity map imaging X-ray photoelectron spectroscopy (A-VMI-XPS), a surface-sensitive and in situ technique, to study the partitioning and structure of a model surfactant (n-octyl β-D-thioglucopyranoside, OTG) in submicron aqueous aerosols with a mean radius of 116 nm. After accounting for finite size effects and evaporative cooling, we find the bulk concentrations of OTG in the aerosol are depleted 600-fold relative to macroscopic solutions. At the low temperatures of the experiment (ca. ∼188 K), we obtain a Langmuir equilibrium constant of 700 m3 mol-1, which corresponds to 4.9 m3 mol-1 at room temperature; a value that is consistent with prior literature reports. Analysis of photoelectron angular distributions, peak areas, and secondary electron escape barriers indicates a transition from a disordered surfactant layer at low [OTG] to a more structured layer at high [OTG]. These results link macroscopic surface tension measurements to nanoparticle surfactant behavior and demonstrate the utility of A-VMI-XPS in probing surfactant coating structure on free aerosols, with implications for understanding the impact of organic coatings on cloud formation and atmospheric chemistry.

Deal, Alexandra M↗
DOE OSTI2021

Designing high-performance superconductors with nanoparticle inclusions: Comparisons to strong pinning theory

One of the most promising routes for achieving high critical currents in superconductors is to incorporate dispersed, non-superconducting nanoparticles to control the dissipative motion of vortices. However, these inclusions reduce the overall superconducting volume and can strain the interlaying superconducting matrix, which can detrimentally reduce T c . Consequently, an optimal balance must be achieved between the nanoparticle density n p and size d. Determining this balance requires garnering a better understanding of vortex–nanoparticle interactions, described by strong pinning theory. Here, we map the dependence of the critical current on nanoparticle size and density in (Y 0.77 , Gd 0.23 )Ba 2 Cu 3 O 7-δ films in magnetic fields of up to 35 T and compare the trends to recent results from time-dependent Ginzburg–Landau simulations. We identify consistency between the field-dependent critical current J c (B) and expectations from strong pinning theory. Specifically, we find that J c ∝ B -α , where α decreases from 0.66 to 0.2 with increasing density of nanoparticles and increases roughly linearly with nanoparticle size d/ξ (normalized to the coherence length). At high fields, the critical current decays faster (~B -1 ), suggesting that each nanoparticle has captured a vortex. When nanoparticles capture more than one vortex, a small, high-field peak is expected in J c (B). Due to a spread in defect sizes, this novel peak effect remains unresolved here. Finally, we reveal that the dependence of the vortex creep rate S on nanoparticle size and density roughly mirrors that of α, and we compare our results to low-T nonlinearities in S(T) that are predicted by strong pinning theory.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The distribution of neutral hydrogen in the colour–magnitude plane of galaxies

ABSTRACT We present the conditional H i (neutral hydrogen) mass function (HIMF) conditioned on observed optical properties, Mr (r-band absolute magnitude), and Cur (u − r colour), for a sample of 7709 galaxies from Arecibo Legacy Fast ALFA (40 per cent data release – α.40) which overlaps with a common volume in SDSS DR7. Based on the conditional HIMF, we find that the luminous red, luminous blue, and faint blue populations dominate the total HIMF at the high-mass end, knee, and the low-mass end, respectively. We use the conditional HIMF to derive the underlying distribution function of ΩH i (H i density parameter), p(ΩH i), in the colour–magnitude plane of galaxies. The distribution, p(ΩH i), peaks in the blue cloud at $\mathit{ M_{{r}}}^{\text{max}}=-19.25, \mathit{ C_{{ur}}}^{\text{max}}=1.44$ but is skewed. It has a long tail towards faint blue galaxies, and luminous red galaxies. We argue that p(ΩH i) can be used to reveal the underlying relation between cold gas, stellar mass, and the star formation rate in an unbiased way, that is, the derived relation does not suffer from survey or sample selection.

Dutta, Saili (ORCID:0000000288585845)↗
DOE OSTI2021

GT Flex: A Coordinated Multi-Building Pilot Study

Buildings are a significant and untapped resource for providing utility electric grid services. Recent studies have estimated that buildings could reduce the peak demand on the electric grid in the U.S. by almost 25% through effective combinations of energy efficiency (EE) measures and load flexibility strategies (Langevin et al. 2021). The U.S. Department of Energy (DOE) has established a goal to triple energy efficiency and demand flexibility in both residential and commercial sections by 2030 compared to 2020 levels (Satchwell et al. 2021). Such findings place buildings alongside electric vehicles, photovoltaics, electric batteries, and other distributed energy resources (DERs) as primary technologies needed for supporting high renewable energy generation grids. Coordinating and optimizing multiple buildings and other DERs is more beneficial and valuable when compared with individual buildings and DERs operating as siloed resources, uncoordinated with others (Olgyay et al. 2020). A pilot study at the Georgia Institute of Technology (GIT) was conducted to evaluate value propositions of a multi-building scale project seeking carbon reduction, energy efficiency and grid-interactive capabilities, by demonstrating the means by which stakeholders can determine the technical and financial merits of grid-interactivity and energy efficiency technologies coordinated across multiple assets. The study focused on analyzing technical feasibility of deploying thermal load flexibility strategies at the multi-building scale, coordinated to not exceed existing infrastructure constraints at the pilot site. Results show that campus can provide 3-3.5 MW of potential load shed over a 4-hour event window through coordinated dispatch of thermal cooling load flexibility without exceeding existing infrastructure capacities. Under future high renewable scenarios, this thermal flexibility resource is also valuable when coordinated to reduce curtailment of intermittent renewables. Economic analyses were performed to effectively communicate various value propositions of grid-interactive and efficient building (GEB) thermal flexibility strategies. Load flexibility presents a financial value proposition to the campus today. By conducting rationalized, coordinated dispatch in response to real time price (RTP) fluctuations, the campus can benefit materially from daily price arbitrage. The RTP signal acts as an aggregating mechanism between the utility and customer to call on demand flexibility resources, with a large portion of the benefit deriving from a relatively small number of days. Realizing and maximizing this benefit with thermal load flexibility requires careful attention to the timing of pricing signals and parameterization of dispatch to overcome efficiency penalties. Grid value and signals are expected to evolve over time, and thermal load flexibility shows potential to adapt dispatch logic to support intermittent renewable generation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

High-Efficiency, Medium-Voltage-Input, Solid-State-Transformer-Based 400-kW/1000V/400A Extreme Fast Charger for Electric Vehicles

Range anxiety and long battery charging time continue to be critical challenges to mass adaptation of EVs. A major identified gap to wider adoption of BEVs is the ability and availability to refuel quickly or to fast charge. Studies have shown that in areas where drivers have access to 50-kW or 120-kW fast charge stations, annual electric vehicle (EV) miles traveled (i.e., eVMT) increased by over 25%, even in cases where fast charging was used for 1% to 5% of total charging events. Charge stations of higher power not only alleviate the “range anxiety” and reduce the driver’s waiting time, but also requires less investment. Michigan Energy Office completed a study in early 2019 titled “Electric Vehicle Charger Placement Optimization in Michigan: Phase I – Highways”. This study finds a system with 150kW chargers, though more expensive individually, actually has lower total system cost when compared to a 50kW charging system when serving the same battery size EV. To be truly competitive to the ICEV refueling experience, even higher power stations are necessary. However, high power charge stations would create large power draws from the grid. If this occurs during peak demand periods, grid capacity could be overloaded. This problem needs to be addressed to reduce the impact on the electric utility infrastructure. The main goal of this project is to develop a 400-kW/400-A XFC system targeting total efficiency of 96.5 percent from the MVAC grid to a vehicle. The novel SST power cell topology, combined with a new silicon carbide (SiC) MOSFET device, enables a 3.5 percent improvement in system efficiency, a 50-percent smaller equipment footprint, and four times less weight than today’s DCFC systems. The SST technology would directly utilize MVAC at 4.8-kV or 13.2-kV. This would eliminate the line frequency transformer (LFT), which steps down medium-voltage AC to 3-Phase 480-V line-to-line voltage in current DCFC systems.

42 ENGINEERING↗
DOE OSTI2020

SN 2018fif: The Explosion of a Large Red Supergiant Discovered in Its Infancy by the Zwicky Transient Facility

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. Here, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of $R={744.0}_{-128.0}^{+183.0}\,{R}_{\odot }$ and an ejected mass of ${M}_{\mathrm{ej}}={9.3}_{-5.8}^{+0.4}\,{M}_{\odot }$. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Rest-frame Submillimeter Spectrum of High-redshift, Dusty, Star-forming Galaxies from the SPT-SZ Survey

We present the average rest-frame spectrum of the final catalog of dusty star-forming galaxies (DSFGs) selected from the South Pole Telescope's SPT-SZ survey and measured with Band 3 of the Atacama Large Millimeter/submillimeter Array. This work builds on the previous average rest-frame spectrum, given in Spilker et al. (2014) for the first 22 sources, and is comprised of a total of 78 sources, normalized by their respective apparent dust masses. The spectrum spans 1.9 < z < 6.9 and covers rest-frame frequencies of 240–800 GHz. Combining this data with low-J CO observations from the Australia Telescope Compact Array, we detect multiple bright line features from 12 CO, [Ci], and H 2 O, as well as fainter molecular transitions from 13 CO, HCN, HCO + , HNC, CN, H 2 O + , and CH. We use these detections, along with limits from other molecules, to characterize the typical properties of the interstellar medium (ISM) for these high-redshift DSFGs. We are able to divide the large sample into subsets in order to explore how the average spectrum changes with various galaxy properties, such as effective dust temperature. We find that systems with hotter dust temperatures exhibit differences in the bright 12 CO emission lines, and contain either warmer and more excited dense gas tracers or larger dense gas reservoirs. These observations will serve as a reference point to studies of the ISM in distant luminous DSFGs (L IR > 10 12 L ⊙ ), and will inform studies of chemical evolution before the peak epoch of star formation at z = 2–3.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

SQUID Readout of a High-$Q$ Superconducting $LC$ Resonator

We demonstrate the readout of a superconducting $LC$ resonator, with $Q \approx 2 \times 10^{6}$ at approximately 250 kHz, using a dc SQUID followed by a SQUID series-array amplifier readout chain. We find that $Q$ depends on the SQUID flux-bias point, increasing near the shallow-slope point of the SQUID modulation curve and decreasing near the steep-slope point, consistent with the previously observed SQUID damping effects. From this variation, we infer the SQUID effective input impedance. We further infer the effective temperature of the resonator circuit from the resonance peak in the noise spectrum and show that it also depends on the SQUID flux-bias point, suggesting a possible contribution from SQUID back-action noise. This system provides a working prototype for future experiments based on lumped-element resonators with SQUID readout, in particular low-mass axion searches within the DMRadio program.

FOS: Physical sciences↗
DOE OSTI2022

Dryness controls temperature-optimized gross primary productivity across vegetation types

Temperature response of gross primary productivity (GPP) is a well-known property of ecosystem, but GPP at the optimum temperature (GPP_T opt ) has not been fully discussed. Our understanding of how GPP_T opt responds to warming and water availability is highly limited. Here, in this study, we analyzed data at 326 globally distributed eddy covariance sites (79°N-37°S), to identify controlling factors of GPP_T opt . Although GPP_T opt was significantly influenced by soil moisture, global solar radiation, mean annual temperature, and vapor pressure deficit in a non-linear pattern (R 2 = 0.47), the direction and magnitude of these climate variables’ effects on GPP_T opt depend on the dryness index (DI), a ratio of potential evapotranspiration to precipitation. The spatial pattern showed that soil moisture did not affect GPP_T opt across energy-limited sites with DI < 1 while dominated GPP_T opt across water-limited sites with DI >1. The temporal pattern showed that GPP_T opt was lowered by warming or low precipitation in water-limited sites while energy-limited sites tended to maintain a stable GPP_T opt regardless of changes in air temperature. Vegetation types in humid climates tended to have higher GPP_T opt and were more likely to benefit from a warmer climate since it was not restricted by water conditions. This study highlights that the response of GPP_T opt to global warming depends on the dryness conditions, which explains the nonlinear control of water and temperature over GPP_T opt . Our finding is essential to realistic prediction of terrestrial carbon uptake under future climate and vegetation conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Kinetic stability of Chapman–Enskog plasmas

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$ , and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$ , the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma ( viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$ , the electron inertial scale $d_e$ and the plasma $\beta$ . This work is of general consequence in emphasising the fact that high- $\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$ . In this situation, the determination of transport coefficients via the standard CE approach is not valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗