Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Index”

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 127 records · Page 7

Fuel Effects on Advanced Compression Ignition Load Limits

In order to increase the efficiency of light-duty gasoline engines, the Co-Optimization of Fuels and Engines (Co-Optima) initiative from the U.S. Department of Energy is investigating multi-mode combustion strategies. Multi-mode combustion can be describe as using conventional spark-ignited combustion at high loads, and at the part-load operating conditions, various advanced compression ignition (ACI) strategies are being investigated to increase efficiency. Of particular interest to the Co-Optima initiative is the extent to which optimal fuel properties and compositions can enable higher efficiency ACI combustion over larger portions of the operating map. Extending the speed-load range of these ACI modes can enable greater part-load efficiency improvements for multi-mode combustion strategies. In this manuscript, we investigate fuel effects for six different fuels, including four with a research octane number (RON) of 98 and differing fuel chemistries, iso-octane, and a market representative E10 fuel, on the load limits for two different ACI strategies: spark-assisted compression ignition (SACI) and partial fuel stratification-gasoline compression ignition (PFS-GCI) operation. Experimental results show that limits to intake boosting limit high load operation for most fuels, but high smoke emissions for high particulate matter index (PMI) fuels under SACI conditions could also be a limitation. Contrastingly, low load is limited by combustion efficiency, but these effects have more pronounced variation with fuel chemistry for PFS-GCI than with SACI. Additional, distinct effects affecting autoignition timing and peak heat release at higher speeds were identified for fuels having different low temperature heat release (LTHR) propensities for both ACI modes.

33 ADVANCED PROPULSION SYSTEMS↗

Method to produce colorless, high porosity, transparent polymer aerogels

A method to produce a polymer gel includes dissolving precursors in a solvent to form a precursor solution, the precursors including polymer precursors, a stable free radical, one or more initiating radicals, and one or more stable free radical control agents, and heating the precursor solution to a temperature of polymerization to produce a cross-linked gel. A dried polymer aerogel has a Brunauer-Emmett Teller (BET) surface area over 100 m2/g, porosity of greater than 10%, visible transparency greater than 20%, color rendering index of over 20%, and average pore size of less than 100 nm.

Chintapalli, Mahati↗

Analysis of MHD stability and active mode control on KSTAR for high confinement, disruption-free plasma

Long-pulse plasma operation at high normalized beta, $β_N$, above the $\textit{n}$= 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong $\textit{m/n}$= 2/1 tearing instabilities at $β_N$ lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. We report the experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different $β_N$ is computed by using the resistive DCON code and the M3D-C 1 code. Equilibrium reconstructions that include constraints from internal profile diagnostics, and computed fast particle pressure are used as input for reliable computation of stability. The classical tearing stability index, Δ', from resistive DCON is computed to be unstable when the island is fully saturated with large amplitude, while the unstable mode is not computed by M3D-C 1 . The modified Rutherford equation (MRE) describing the evolution of the island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. The MRE model estimates a saturated island width corresponding to ~10% of the plasma minor radius for equilibrium at high $β_N$ 3 having a stable Δ' from the resistive DCON. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm has been completed and enabled on KSTAR. To accurately determine the $\textit{n}$= 1component produced by RWMs, an algorithm has been developed that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors. Use of multiple toroidal sensor arrays is enabled by modifying the sensor toroidal angles assumed in mode decomposition to include the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on stability, transport, and control provides the required foundation for disruption prediction and avoidance research on KSTAR.

resistive wall mode (RWM)↗

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗

Dopant Segregation in Single-Crystal Optical Fiber Grown via the Laser-Heated Pedestal Growth Technique

Single crystal (SC) optical fiber has promising potential to be used for optical fiber sensing applications in harsh conditions due to its robustness to high temperature, high radioactivity, and resistance to chemical corrosion as compared to optical sensors using silica fiber. However, SC fiber grown via the laser-heated pedestal growth (LHPG) technique innately does not have a core-cladding structure found in standard optical fiber, resulting in optical losses. In this work we investigate optimization of the growth parameters of a two LHPG process used to grow SC fiber with a graded index via introduction of dopants to the feedstock material. Feedstock material is fabricated with the first LHPG device, then sol-gel dip-coated to create outer films of dopant material. The dip-coated feedstock is used to grow SC fiber in which segregation of the dopant constituents occurs, resulting in a graded index of refraction across the fiber, and an effective core-cladding structure. Hardware and software improvements to both LHPG systems are presented and the growth parameters for short pieces of ~320-330 μm diameter YAG fiber has been established. Characterization techniques/procedures have also been established for future grown SC fiber. These improvements and preparations are anticipated to result in a significant increase in grown fiber quality with a similar growth rate to that previously established.

Lander, Gary↗

Dopant Segregation in Single-Crystal Optical Fiber Grown via the Laser-Heated Pedestal Growth Technique

Single crystal (SC) optical fiber has promising potential to be used for optical fiber sensing applications in harsh conditions due to its robustness to high temperature, high radioactivity, and resistance to chemical corrosion as compared to optical sensors using silica fiber. However, SC fiber grown via the laser-heated pedestal growth (LHPG) technique innately does not have a core-cladding structure found in standard optical fiber, resulting in optical losses. In this work we investigate optimization of the growth parameters of a two LHPG process used to grow SC fiber with a graded index via introduction of dopants to the feedstock material. Feedstock material is fabricated with the first LHPG device, then sol-gel dip-coated to create outer films of dopant material. The dip-coated feedstock is used to grow SC fiber in which segregation of the dopant constituents occurs, resulting in a graded index of refraction across the fiber, and an effective core-cladding structure. Hardware and software improvements to both LHPG systems are presented and the growth parameters for short pieces of ~320-330 μm diameter YAG fiber has been established. Characterization techniques/procedures have also been established for future grown SC fiber. These improvements and preparations are anticipated to result in a significant increase in grown fiber quality with a similar growth rate to that previously established.

Lander, Gary↗

Earth System Model Overestimation of Cropland Temperatures Scales With Agricultural Intensity

Intensive crop growth can modify regional climate by partitioning energy to latent heating through transpiration, cooling growing season temperatures. Recent work shows that cooling associated with agriculture can dampen anthropogenic warming over breadbasket regions. However, it is unknown whether climate models reproduce crop influences on regional climate, and thus the future risk of extreme climate events over global breadbasket regions. Here, we show that models overestimate growing season temperatures and underestimate evapotranspiration (ET) over global croplands, and that these differences increase with cropped area. We trace this warm and dry difference through each model's representation of the surface energy budget, showing that model differences in transpiration, leaf area index, and the ratio of transpiration to total ET drive the overall effect. While the implications of these model deficiencies for future projections are uncertain, they point to the importance of improving representations of crop-climate processes to better assess breadbasket vulnerability to climate change.

54 ENVIRONMENTAL SCIENCES↗

Seasonal stabilization effects slowed the greening of the Northern Hemisphere over the last two decades

Rising atmospheric CO₂ and warming spring temperatures increase vegetation growth and the terrestrial carbon sink. However, drought, heat stress, phenology, and resource limitations may stabilize or limit theses projected increases. We investigate the balance between these amplifying and stabilizing ecological factors by asking whether enhanced early-season growth leads to continued late-season growth. Using the Moderate Resolution Imaging Spectroradiometer (MODIS) leaf area index (LAI) dataset, we identify three seasonal growth patterns based on early- and peak-season positive LAI anomalies: (1) amplification, where late-season LAI anomalies exceed earlier ones; (2) weak stabilization, where late-season anomalies remain similar or slightly lower; and (3) strong stabilization, where late-season anomalies become negative. Weak and strong stabilization events dominate across 67% and 26% of Northern Hemisphere ecosystems above 30°N, respectively. The absence of any trend in amplifying or stabilizing events suggests stabilizing factors seasonally offset CO₂ and temperature-induced spring greening. Terrestrial biosphere models underestimate strong stabilization and overestimate amplification events. This inconsistency arises from the models’ underestimation late-season LAI sensitivity to precipitation in water-limited regions; overlook negative legacy effects of early enhanced LAI on late-season soil moisture via evapotranspiration losses in energy-limited regions. Our findings suggest water/heat stress and resource limitations limit greening and the land carbon sink.

Zhang, Wen↗

Changes in high-latitude surface energy balance driven by snowpack and vegetation dynamics under warmer climate

With rapid climate warming, expected changes in snowpack and vegetation will alter the seasonal surface albedo of high-latitude ecosystems. The extent to which these albedo changes may affect surface energy balances and thus soil temperatures is uncertain, but represents a potentially important component of ecosystem feedbacks to climate change. Here, we apply a well-tested process-rich ecosystem model, ecosys , to examine changes in seasonal surface albedo and soil temperature driven by climate-induced snowpack and vegetation changes across Alaska under a warmer twenty-first century climate. Under the Representative Concentration Pathway 8.5 climate change scenario, the modeled changes in surface albedo exhibited large seasonal and spatial variations. We found spring albedo decreases driven by increases in snow-free periods (>20 d) and an extended growing season length that resulted in greater gains in leaf area index (LAI) in most parts of Alaska. In contrast, we modeled increases in summer and winter albedo (despite modeled increases in LAI) across much of the boreal forest due to an increased proportion of aspen, which has a higher leaf albedo than the currently dominant black spruce. Modeled latent heat fluxes generally increase across the twenty-first century, particularly during the spring and summer. Overall, climate warming and changes in surface energy fluxes resulted in a 3.5 ± 0.50 °C increase in spatial- and annual-averaged top 10 cm surface soil temperatures across Alaskan ecosystems by the year 2100, with larger increases in tundra than boreal forest regions. We conclude that under warmer climates, seasonal variations in albedo and surface energy fluxes are particularly pronounced during the spring and summer, driven by changes in snowpack and vegetation dynamics.

54 ENVIRONMENTAL SCIENCES↗

DEEPEN: Final 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

Part of the DEEPEN (DE-risking Exploration of geothermal Plays in magmatic ENvironments) project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology, detailed in the journal article below, is based on the method proposed by Poux & O'brien (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., (2022). For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty.

15 GEOTHERMAL ENERGY↗

Experimental investigation of turbulent wake flows in a helically wrapped rod bundle in presence of localized blockages

In nuclear sodium fast reactors, bundles of rods are tightly packed into a triangular lattice, enclosed in a hexagonal duct, and each pin is spirally wrapped with a thin wire. Flow blockages can potentially impact the local flow characteristics and heat transfer mechanisms in the bundle due to its small subchannel size. The effects of the blockage on the flow structures and heat transfer mechanisms are important aspects that require an accurate investigation. In this study, the flow-field characteristics in the vicinity of a blockage located in the exterior subchannel of rod bundles with helically wrapped wires were experimentally investigated. The velocity fields in the exterior subchannel were acquired by applying matched-index-of-refraction and time-resolved particle image velocimetry (TR-PIV) techniques for Reynolds numbers of Re 1 = 4000 and Re 2 = 17 000, i.e., equivalent to Rew 1 = 19 600 and Rew 2 = 83 200, respectively, based on the blockage width. The results from the TR-PIV measurements revealed an arch-shaped vortex with a large flow recirculation and a pair of counter-rotating vortices in the wake region downstream of the blockage, which is commonly observed in the wake flow of bluff bodies. The relative lateral distance and angle between the two vortices decreased when the Reynolds numbers increased. Profiles of maximum turbulence intensity along the shear layers illustrated the transition process including the growth, peak, and decay along the flow direction. From the spectral analysis of the turbulent velocities extracted at points along the shear layer, the Strouhal numbers ( St ) representing the vortex shedding frequency were found to be St = 0.25 and St = 0.56 for the left and right shear layers, respectively. Characteristics of shear layers generated by the blockage in the exterior subchannel were investigated via the two-point cross correlation of fluctuating velocities. The patiotemporal cross correlations of turbulent velocities, computed at points in the region where the left shear layer exhibited rolling effects and vortex breakdowns, were considerably wider and longer. The convection velocity U c was estimated to be ~0.82 U m to 0.93 U m . Proper orthogonal decomposition (POD) analysis was applied to the instantaneous velocity fields to extract the statistically dominant flow structures. It was found that POD modes 2–3 and 4–5 formed the pair modes when the corresponding POD temporal coefficients depicted sinusoidal shapes and exhibited nearly circular orbits in the phase space. Spectral analysis of the POD temporal coefficients confirmed the vortex shedding frequencies detected in the analysis of turbulent velocities.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Boron coordination change in barium borate melts and glasses and its contribution to configurational heat capacity, entropy, and fragility

High-energy x-ray diffraction from molten and glassy BaB 2 O 4 and BaB 4 O 7 has been performed using aerodynamic levitation and laser heating over a wide range of temperatures. Remarkably, even in the presence of a heavy metal modifier dominating x-ray scattering, it was possible to extract accurate values for the tetrahedral, sp 3 , boron fraction, N 4 , which declines with increasing temperature, using bond valence-based mapping from the measured mean B–O bond lengths while accounting for vibrational thermal expansion. These are used within a boron-coordination-change model to extract enthalpies, ΔH, and entropies, ΔS, of isomerization between sp 2 and sp 3 boron. The results for BaB 4 O 7 , ΔH = 22(3) kJ mol -1 boron, ΔS = 19(2) J mol -1 boron K -1 , agree quantitatively with those found previously for Na 2 B 4 O 7 . Analytical expressions for N 4 (J, T) and associated configurational heat capacity, $C_{P}^{conf}$(J, T), and entropy, S conf (J, T), contributions are extended to cover a wide composition range 0 ≤ J = BaO/B 2 O 3 ≤ 3 using a model for ΔH(J) and ΔS(J) derived empirically for lithium borates. Maxima in the $C_{P}^{conf}$(J, T g ) and fragility index contributions are thereby predicted for J ≲ 1, higher than the maximum observed and predicted in N 4 (J, T g ) at J ≃ 0.6. We discuss the utility of the boron-coordination-change isomerization model in the context of borate liquids containing other modifiers and the prospect of neutron diffraction to aid in empirical determinations of modifier-dependent effects, illustrated by new neutron diffraction data on Ba 11 B 4 O 7 glass, its well-known α-polymorph, and lesser-known δ-phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ion Acceleration and the Development of a Power-law Energy Spectrum in Magnetic Reconnection

How charged particles are accelerated efficiently and form a power-law energy spectrum in magnetic reconnection is a problem that is not well understood. In a previous paper, it was shown that the electron Kelvin–Helmholtz instability (EKHI) in force-free magnetic reconnection generates fast-expanding vortices that can accelerate electrons in a few tens of ion gyroperiods (less than 1 ms in the solar corona) to form a power-law energy distribution. In this paper, we present a particle-in-cell (PIC) simulation study of ion acceleration in force-free magnetic reconnection in the presence of the EKHI-induced turbulence. We find that ions are not significantly accelerated by the EKHI-induced stochastic electric field until the magnetic vortices expand to sizes comparable to the ion gyroradius. The Alfvén waves generated by the EKHI couple with the magnetic vortices, leading to resonance between the ions inside the magnetic vortices and Alfvén waves and enhanced ion heating. The induced Alfvén wave resonance results in a broken power-law energy spectrum with a breakpoint at $\sim {m}_{i}{v}_{A}^{2}$, where v A is the Alfvén velocity. Here, we show that the process that forms the nonthermal tail is a second-order Fermi mechanism and the mean spectral index is α = (1 + 4 a 2 D/R) /2, where D is the spatial scale of the inductive electric field, R is that of vortices, and a = B g /B 0 , with ratio of guide field B g and asymptotic B 0 .

79 ASTRONOMY AND ASTROPHYSICS↗

Superconductivity in MgTa 2 S 5

Superconductivity in the new compound, MgTa 2 S 5 , is reported. Its crystal structure was indexed to space group $C$2/$m$, which is monoclinic with lattice parameters $a$ = 5.784(5) $\mathring{A}$, $b$ = 12.680(3) $\mathring{A}$, and $c$ = 6.106(5) $\mathring{A}$, with $β$ = 112.64°. Electrical resistivity reveals metallic behavior with $ρ$(296 K) = 5.2 mΩ cm and bulk superconductivity below $T$ $c$ = 1 K. Measurements of the specific heat C reveal the electronic specific heat coefficient $γ$ = 11.1(5) mJ/mol K 2 , a jump in $C$ at $T$ $c$ illustrating bulk superconductivity, and an energy gap associated with the superconducting state of $E$ $g$ = 0.25(3) meV. MgTa 2 S 5 can be classified as a BCS superconductor. In conclusion, an unidentified minority phase present in the sample (~ 4%) exhibits a superconducting phase transition near 2.6 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Will cities keep getting hotter? The interplay of urban expansion and greening reshapes future urban heat trajectories

Urban heat islands (UHIs) pose growing risks to public health, infrastructure, and resilience. While often assumed to intensify with urban growth, dynamic changes in urban expansion and vegetation greenness complicate UHI trajectories, which remain poorly understood. This study investigated the interplay of urban expansion and greenness change on UHI spatial profile across 36 Chinese megacities during 2003–2018 using multiple satellite products. We introduce a framework that classifies urban areas into four dynamic development pathways based on impervious surface area (ISA) and enhanced vegetation index (EVI) trends: urbanized-greening, urbanized-browning, urbanizing-greening (UingG), and urbanizing-browning (UingB). While most urbanized centers exhibited greening driven by targeted initiatives and urbanizing suburbs showed browning due to vegetation loss, about 30% urban areas showed the reverse pattern, revealing overlooked complexity in urban development. Urban expansion and browning strengthened UHI in suburban areas, whereas greening initiatives mitigated UHI in urban center and mitigated UHI enhancement in suburban areas. Slowed warming in urban centers together with accelerated warming in suburban areas flattened the temperature gradient between urban centers and suburbs. This dynamic expanded the spatial extent of elevated temperatures and reshaped the classic urban-to-rural UHI profile into a flatter form. In UingB areas, UHI intensification was jointly driven by increase in ISA, vegetation loss, and their interaction, while in UingG areas, EVI increases and a negative interaction together offset over half of the warming driven by urban expansion. These findings reveal that UHI evolution is not unidirectional but depends on localized urbanization and greening dynamics, offering pathways for strategic heat mitigation.

greening initiative↗

New Record Ocean Temperatures and Related Climate Indicators in 2023

The global physical and biogeochemical environment has been substantially altered in response to increased atmospheric greenhouse gases from human activities. In 2023, the sea surface temperature (SST) and upper 2000 m ocean heat content (OHC) reached record highs. The 0–2000 m OHC in 2023 exceeded that of 2022 by 15 ± 10 ZJ (1 Zetta Joules = 10 21 Joules) (updated IAP/CAS data); 9 ± 5 ZJ (NCEI/NOAA data). The Tropical Atlantic Ocean, the Mediterranean Sea, and southern oceans recorded their highest OHC observed since the 1950s. Associated with the onset of a strong El Niño, the global SST reached its record high in 2023 with an annual mean of ~0.23°C higher than 2022 and an astounding > 0.3°C above 2022 values for the second half of 2023. The density stratification and spatial temperature inhomogeneity indexes reached their highest values in 2023.

54 ENVIRONMENTAL SCIENCES↗

Spatial heterogeneity effects on land surface modeling of water and energy partitioning

Abstract. Understanding the influence of land surface heterogeneity on surface water and energy fluxes is crucial for modeling earth system variability and change. This study investigates the effects of four dominant heterogeneity sources on land surface modeling, including atmospheric forcing (ATM), soil properties (SOIL), land use and land cover (LULC), and topography (TOPO). Our analysis focused on their impacts on the partitioning of precipitation (P) into evapotranspiration (ET) and runoff (R), partitioning of net radiation into sensible heat and latent heat, and corresponding water and energy fluxes. An initial set of 16 experiments were performed over the continental US (CONUS) using the E3SM land model (ELMv1) with different combinations of heterogeneous and homogeneous datasets. The Sobol' total and first-order sensitivity indices were utilized to quantify the relative importance of the four heterogeneity sources. Sobol' total sensitivity index measures the total heterogeneity effects induced by a given heterogeneity source, consisting of the contribution from its own heterogeneity (i.e., the first-order index) and its interactions with other heterogeneity sources. ATM and LULC are the most dominant heterogeneity sources in determining spatial variability of water and energy partitioning, mainly contributed by their own heterogeneity and slightly contributed by their interactions with other heterogeneity sources. Their heterogeneity effects are complementary, both spatially and temporally. The overall impacts of SOIL and TOPO are negligible, except TOPO dominates the spatial variability of R/P across the transitional climate zone between the arid western and humid eastern CONUS. Accounting for more heterogeneity sources improves the simulated spatial variability of water and energy fluxes when compared with ERA5-Land reanalysis dataset. An additional set of 13 experiments identified the most critical components within each heterogeneity source, which are precipitation, temperature, and longwave radiation for ATM, soil texture, and soil color for SOIL and maximum fractional saturated area parameter for TOPO.

54 ENVIRONMENTAL SCIENCES↗

Nature and Scalings of Density Fluctuations of Compressible Magnetohydrodynamic Turbulence with Applications to the Solar Wind

Abstract The solar wind is a magnetized and turbulent plasma. Its turbulence is often dominated by Alfvénic fluctuations and often deemed as nearly incompressible far away from the Sun, as shown by in situ measurements near 1 au. However, for solar wind closer to the Sun, the plasma β decreases (often lower than unity) while the turbulent Mach number M t increases (can approach unity, e.g., transonic fluctuations). These conditions could produce significantly more compressible effects, characterized by enhanced density fluctuations, as seen by several space missions. In this paper, a series of 3D MHD simulations of turbulence are carried out to understand the properties of compressible turbulence, particularly the generation of density fluctuations. We find that, over a broad range of parameter space in plasma β , cross helicity, and polytropic index, the turbulent density fluctuations scale linearly as a function of M t , with the scaling coefficients showing weak dependence on parameters. Furthermore, through detailed spatiotemporal analysis, we show that the density fluctuations are dominated by low-frequency nonlinear structures, rather than compressible MHD eigenwaves. These results could be important for understanding how compressible turbulence contributes to solar wind heating near the Sun.

79 ASTRONOMY AND ASTROPHYSICS↗