Search NASA⌕ Search

SEARCH · Search NASA

Results for “Temperature variability”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Developing and Evaluating Metrics for Demand Flexibility in Buildings: Comparing Simulations and Field Data

Building demand flexibility (DF) has attracted significant attention in recent years among researchers, technology developers and control companies, aggregators, utilities, and many others. There are numerous challenges with today’s electricity systems such as managing peak demand capacity and integrating variable renewables into the grid. Flexible building loads can provide various grid services to help reduce electricity costs, smooth out renewables intermittency and balance supply and demand. Recognizing this, the US DOE is leading the Grid-interactive Efficient Buildings (GEB) initiative which includes research to evaluate the potential, availability and timing of flexible loads. In this paper we present load shed metrics for three building types – medium office, large office and retail store – and compare prototype simulation results with measured data from 12 actual buildings that participated in hot summertime utility demand response (DR) events. The DR strategies included varying zone temperature and reducing light levels. The magnitude of a key DF metric, “demand decrease intensity” (or “shed intensity”) (W/ft2), between the simulation results and field data are similar (14-32% differences) for both mean and median values, though the field data show much larger variation among DR events. The coefficient p-values from linear regression model tests showed that outside air temperature is a significant variable for the whole building shed intensity when the resetting zone temperature strategy is deployed. These findings support the concept of using prototype building simulation to estimate building DF and expanding future simulation research to additional building types and climate zones.

Liu, Jingjing↗
DOE OSTI2023

Response of Extreme Rainfall to Atmospheric Warming and Wetting: Implications for Hydrologic Designs Under a Changing Climate

Understanding the processes of rainfall extremes and their response to anthropogenic climate change is pivotal for improved adaptation of unprecedented flood hazards around the world. Here we take the record-breaking 20 July 2021 storm over central China as an example. We investigate the response of this particular storm to atmospheric warming (i.e., increase in air temperature) and wetting (i.e., increase in atmospheric moisture content) based on a series of convection-permitting model simulations. Our results show non-monotonic changes of the space-time rainfall variability to either increased temperature or atmospheric moisture content. The most extreme rain rate is produced when relative humidity is increased by 20~40 % or temperature is increased by less than 2°C. The non-monotonic rainfall response is more clearly revealed at fine spatial (100-1000 km 2 ) and temporal scales (less than 6 hours) rather than over the entire domain (~10 4 km 2 ) and aggregated over the storm duration (around two days). This is mainly attributable to the distinct feedbacks from atmospheric dynamics (i.e., moisture convergence and interaction with regional topography) rather than regulated by thermodynamic changes alone. Furthermore, atmospheric warming poses notable changes in the vertical structure of storm cells, contributing to reduced areal reduction factors at small spatial scales and short durations, while atmospheric wetting additionally modifies storm evolution properties. Our modeling analyses challenge the existing practices for hydrologic designs under a changing climate, highlighting particular vulnerability for cities or small basins to short-duration rainfall extremes and the resultant flash flood hazards.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparison of two pressure–temperature equilibration methods

We compare and contrast the traditionally used method of solving the pressure–temperature equilibration problem in hydrodynamics, where specific internal energy and density are considered independent variables, with a different method where pressure and temperature are independent variables. With the goal of examining the robustness of the two methods as the number of components increases, we examine 2-, 4-, 6-, and 8-component systems. After equilibrating more than 10 4 initial conditions for each system using both methods, we demonstrate that the latter method constrains the search space by lowering its dimensionality and forces a better initial guess, resulting in a higher probability of convergence to solution with fewer, cheaper iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Thermodynamically consistent semi-compressible fluids: a variational perspective

This paper presents (Lagrangian) variational formulations for single and multicomponent semi-compressible fluids with both reversible (entropy-conserving) and irreversible (entropy-generating) processes. Semi-compressible fluids are useful in describing low-Mach dynamics, since they are soundproof. These models find wide use in many areas of fluid dynamics, including both geophysical and astrophysical fluid dynamics. Specifically, the Boussinesq, anelastic and pseudoincompressible equations are developed through a unified treatment valid for arbitrary Riemannian manifolds, thermodynamic potentials and geopotentials. By design, these formulations obey the 1st and 2nd laws of thermodynamics, ensuring their thermodynamic consistency. This general approach extends and unifies existing work, and helps clarify the thermodynamics of semi-compressible fluids. To further this goal, evolution equations are presented for a wide range of thermodynamic variables: entropy density s, specific entropy η, buoyancy b, temperature T, potential temperature θ and a generic entropic variable χ; along with a general definition of buoyancy valid for all three semicompressible models and arbitrary geopotentials. Finally, the elliptic equation for the pressure perturbation (the Lagrange multiplier that enforces semi-compressibility) is developed for all three equation sets in the case of reversible dynamics, and for the Boussinesq/anelastic equations in the case of irreversible dynamics; and some discussion is given of the difficulty in formulating it for the pseudoincompressible equations with irreversible dynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Northern Hemisphere Extratropical Cyclone Activity in the Twentieth Century Reanalysis Version 3 (20CRv3) and Its Relationship with Continental Extreme Temperatures

In this study, we detect and track extratropical cyclones using 6-hourly mean sea level pressure data taken from the Twentieth Century Reanalysis version 3 (20CRv3) over the period 1951–2015 and compare them with those in the Interim and fifth generation of ECMWF reanalyses over the period 1979–2018. Three indices were employed to characterize cyclone activity, including cyclone count, cyclone intensity, and a cyclone activity index (CAI) that combines the count and intensity. The results show that the cyclone indices in the three datasets have comparable annual climatologies and seasonal evolution over the northern extratropical land and ocean in recent decades. Based on the cyclone indices over the period 1951–2010 in 80 ensemble members of 20CRv3, cyclone count and intensity are negatively correlated in winter and tend to be positively and weakly correlated in summer. The interannual CAI variability is dominated by the cyclone count variability. Regional mean cyclone activity can be well represented using the ensemble average cyclone index. We then examined the linkage of the cyclone activity in 20CRv3 and observed cold and warm extremes over Eurasia and North America over the period 1951–2010. In winter, the principal components of interannual cold and warm extreme anomalies are more correlated with the regional mean cyclone count index over Eurasia, while they are more correlated with the cyclone intensity index over North America. The temperature anomalies associated with the regional and ensemble mean cyclone count index explain about 10% (20%) of interannual cold (warm) extreme variances averaged over Eurasia. The temperature anomalies associated with the mean cyclone intensity explain about 10% of interannual cold and warm extreme variances over North America. Large-scale atmospheric circulation anomalies in association with cyclone activity and the induced temperature advection drive temperature anomalies over Eurasia and North America. In summer, circulation and thermal advection anomalies associated with cyclone activity are weak over the two continents. Hence, that season’s relationship between cyclone activity and extreme temperature variability is weak.

20CRv3↗
DOE OSTI2024

Node Distortions in UiO-66 Inform Negative Thermal Expansion Mechanisms: Kinetic Effects, Frustration, and Lattice Hysteresis

In metal–organic frameworks (MOFs) the interplay between the dynamics of individual components and how these are constrained by the extended lattice can yield unusual emergent phenomena. For the archetypal Zr-MOF, UiO-66, we explore the cooperative dynamics of a Zr-node transformation that gives rise to negative thermal expansion (NTE). Here, using in situ synchrotron X-ray scattering, with powder diffraction and pair distribution function (PDF) analyses, we identify lattice hysteresis and a thermal ramp-rate-dependence of the thermal expansion. Specifically, kinetic trapping of distorted node states formed at high temperature, leads to broad variability in the apparent thermal expansion which ranges from large positive to large negative thermal expansion with coefficients of thermal expansion (CTE) from +45 to –80 × 10 –6 K –1 . Time-resolved relaxation studies at selected temperatures suggest that when equilibrated UiO-66 is intrinsically NTE, with a CTE of –35 × 10 –6 K –1 . Kinetic trapping of the node-distorted state following high temperature activation has broad implications for characterization and applications of these Zr-MOFs; the nonequilibrium node state depends on the thermal history of the sample with quench vs slow cooling likely to impact gas binding, pore volume, and accessible catalytic sites.

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

Impacts of humic substances, elevated temperature, and UVB radiation on bacterial communities of the marine sponge Chondrilla sp

Abstract Sponges are abundant components of coral reefs known for their filtration capabilities and intricate interactions with microbes. They play a crucial role in maintaining the ecological balance of coral reefs. Humic substances (HS) affect bacterial communities across terrestrial, freshwater, and marine ecosystems. However, the specific effects of HS on sponge-associated microbial symbionts have largely been neglected. Here, we used a randomized-controlled microcosm setup to investigate the independent and interactive effects of HS, elevated temperature, and UVB radiation on bacterial communities associated with the sponge Chondrilla sp. Our results indicated the presence of a core bacterial community consisting of relatively abundant members, apparently resilient to the tested environmental perturbations, alongside a variable bacterial community. Elevated temperature positively affected the relative abundances of ASVs related to Planctomycetales and members of the families Pseudohongiellaceae and Hyphomonadaceae. HS increased the relative abundances of several ASVs potentially involved in recalcitrant organic matter degradation (e.g., the BD2-11 terrestrial group, Saccharimonadales, and SAR202 clade). There was no significant independent effect of UVB and there were no significant interactive effects of HS, heat, and UVB on bacterial diversity and composition. The significant, independent impact of HS on the composition of sponge bacterial communities suggests that alterations to HS inputs may have cascading effects on adjacent marine ecosystems.

Stuij, Tamara M. (ORCID:0000000152170798)↗
DOE OSTI2024

Measuring thermal diffusivity and gap conductance in uranium nitride and Zircaloy relevant for microreactor applications

Heat transfer across nuclear fuels and structural interfaces is an important factor for evaluating the performance of nuclear power systems. Specifically, heat generated as nuclear fuel fissions must be transported through the cladding material and through the reactor to reach the steam turbine for power generation. As new microreactor designs emerge, maximizing the efficiency of this heat transfer process becomes crucial to make them commercially viable. This article examines thermal diffusivity and gap conductance in uranium nitride (UN) fuel and Zircaloy-4 (Zry4) cladding using light flash analysis (LFA). Thermal diffusivity measurements were made on monolithic UN pellets and Zry4 exposed to carbon at peak operating temperatures of microreactors and show that carbon ingress has a minimal effect on thermal diffusivity when compared with identical materials not exposed to carbon. Evaluation of gap conductance at the UN-Zry4 interface was done using one-dimensional two-layer thermal transport models as a function of applied pressure. Here the results show that increasing pressure on the UN-Zry4 interface leads to gains in gap conductance per unit area in fuel-cladding assemblies at microreactor operating temperatures. While many other variables are expected to influence UN-Zry4 interfacial gap conductance (e.g. contact surface roughness, porosity, localized heating, environmental gas pressure), the work offers a demonstration of using a conventional LFA apparatus to determine this parameter at elevated temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Mars Upper Atmospheric Temperature and Atomic Oxygen Density Derived from the O i 130.4 nm Emission Observed by NASA’s MAVEN Mission

The abundance of atomic oxygen is a major unknown for the Mars upper atmosphere, and systematic measurements of this key species is a primary objective of NASA’s Mars Atmosphere and Volatile Evolution (MAVEN) mission. Moreover, the Martian dayside temperature structure and its variability has been a subject of considerable debate in the past. In this study, an inversion analysis of the periapsis limb scans of the O i 130.4 nm emission observed by MAVEN is performed to quantify the neutral temperature and atomic oxygen abundance in the Mars upper atmosphere. Our analysis yields unprecedented resolution of the spatial variation of Mars atomic oxygen, showing that its exobase density decreases with the solar zenith angle (SZA) from ∼0° to 60° by about an order of magnitude, from ∼10{sup 8} to ∼10{sup 7} cm{sup −3} near perihelion (2015 April) and from ∼10{sup 7} to ∼10{sup 6} cm{sup −3} near aphelion (2017 November). For a given SZA the exobase densities near perihelion are about an order of magnitude higher than those near aphelion. The [O]/[CO{sub 2}] mixing ratio at 130 km altitude can vary from ∼0.1% to 4.6%, which also decreases with increasing SZA. Moreover, it is shown that the Mars exobase temperature does not vary significantly with SZA from ∼0° to ∼60°, which is estimated to be ∼239.2 ± 27.6 K near perihelion and ∼162.4 ± 19.9 K near aphelion, with a difference of 76.8 K between these two seasonal extremes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Eulerian finite element implementations of a dislocation density-based continuum model

In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling Spatial Distribution of Snow Water Equivalent by Combining Meteorological and Satellite Data with Lidar Maps

Abstract An accurate characterization of the water content of snowpack, or snow water equivalent (SWE), is necessary to quantify water availability and constrain hydrologic and land surface models. Recently, airborne observations (e.g., lidar) have emerged as a promising method to accurately quantify SWE at high resolutions (scales of ∼100 m and finer). However, the frequency of these observations is very low, typically once or twice per season in the Rocky Mountains of Colorado. Here, we present a machine learning framework that is based on random forests to model temporally sparse lidar-derived SWE, enabling estimation of SWE at unmapped time points. We approximated the physical processes governing snow accumulation and melt as well as snow characteristics by obtaining 15 different variables from gridded estimates of precipitation, temperature, surface reflectance, elevation, and canopy. Results showed that, in the Rocky Mountains of Colorado, our framework is capable of modeling SWE with a higher accuracy when compared with estimates generated by the Snow Data Assimilation System (SNODAS). The mean value of the coefficient of determination R 2 using our approach was 0.57, and the root-mean-square error (RMSE) was 13 cm, which was a significant improvement over SNODAS (mean R 2 = 0.13; RMSE = 20 cm). We explored the relative importance of the input variables and observed that, at the spatial resolution of 800 m, meteorological variables are more important drivers of predictive accuracy than surface variables that characterize the properties of snow on the ground. This research provides a framework to expand the applicability of lidar-derived SWE to unmapped time points. Significance Statement Snowpack is the main source of freshwater for close to 2 billion people globally and needs to be estimated accurately. Mountainous snowpack is highly variable and is challenging to quantify. Recently, lidar technology has been employed to observe snow in great detail, but it is costly and can only be used sparingly. To counter that, we use machine learning to estimate snowpack when lidar data are not available. We approximate the processes that govern snowpack by incorporating meteorological and satellite data. We found that variables associated with precipitation and temperature have more predictive power than variables that characterize snowpack properties. Our work helps to improve snowpack estimation, which is critical for sustainable management of water resources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Characterizing soil water content variability across spatial scales from optimized high-resolution distributed temperature sensing technique

Fiber-optic Distributed Temperature Sensing, when combined with the Single-probe Heat-pulse technique can measure soil moisture (θ) across spatial scales. The key limitation of this system is in obtaining the relationship between soil thermal conductivity (λ) and θ for a specific field. Using the Department of Energy Atmospheric Radiation Measurement (ARM) site, this study tested a new methodology to account for the spatial variability in the λ-θ relationship using a Gaussian processes model. The resulting accurate θ measurements (RMSE = 0.03 m 3 m –3 ) were used to characterize the spatial variability of θ across scales and to develop an empirical equation that can correct for the changes in the θ spatial variability observed at different spatial resolutions. In addition, the number of required samples to accurately characterize θ and its variability over scales ranging from 5 m and 350 m were estimated. Finally, these findings provide key information to scale soil moisture from centimeters to hundreds of meters for process understanding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Evaluation of the Near-Surface Variables in the HRRR Weather Model Using Observations from the ARM SGP Site

Abstract The performance of version 4 of the NOAA High-Resolution Rapid Refresh (HRRR) numerical weather prediction model for near-surface variables, including wind, humidity, temperature, surface latent and sensible fluxes, and longwave and shortwave radiative fluxes, is examined over the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) region. The study evaluated the model’s bias and bias-corrected mean absolute error relative to the observations on different time scales. Forecasts of near-surface geophysical variables at five SGP sites (HRRR at 3-km scale) were found to agree well with observations, but some consistent observation–forecast differences also occurred. Sensible and latent heat fluxes are the most challenging variables to be reproduced. The diurnal cycle is the main temporal scale affecting observation–forecast differences of the near-surface variables, and almost all of the variables showed different biases throughout the diurnal cycle. Results show that the overestimation of downward shortwave and the underestimation of downward longwave radiative flux are the two major biases found in this study. The timing and magnitude of downward longwave flux, wind speed, and sensible and latent heat fluxes are also different with contributions from model representations, data assimilation limitations, and differences in scales between HRRR and SGP sites. The positive bias in downward shortwave and negative bias in longwave radiation suggests that the model is underestimating cloud fraction in the study domain. The study concludes by showing a brief comparison with version 3 of the HRRR and shows that version 4 has better performance in almost all near-surface variables. Significance Statement A correct representation of the near-surface variables is important for numerical weather prediction models. This study investigates the capability of the latest NOAA High-Resolution Rapid Refresh (HRRRv4) model in simulating the near-surface variables by comparing against the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) in situ observations. Among others, we find that the surface heat fluxes, such as sensible and latent heat fluxes, are the most difficult variables to be reproduced. This study also shows that the diurnal cycle has the dominant impact on the model’s performance, which means the majority of the outputted near-surface variables have the strong diurnal cycle in their bias errors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Thermal transport properties of IrSbSe

Here we report a thermal transport study of IrSbSe, which crystallizes in a noncentrosymmetric cubic structure with the P2 1 3 space group and shows a narrow-gap semiconducting behavior. The large discrepancy between the activation energy for conductivity [E ρ = 128(2) meV] and for thermopower [E S = 17.7 (9) meV] from 200 to 300 K indicates a polaronic transport mechanism. The electrical resistivity varies as exp (T 0 /T) 1/4 and thermopower varies as T 1/2 at low temperatures, indicating that it evolves into Mott variable-range hopping dominant conduction with decreasing temperature. IrSbSe shows a relatively low value of thermal conductivity (~1.65 W/Km) and thermopower of about 0.24 mV/K around 100 K, yet poor electrical conductivity. On the other hand, a high vacancy defect concentration on both Ir and Sb atomic sites of up to 15% suggests a high defect tolerance and points to the possibility of a future improvement of carrier density by chemical substitution or defect optimization.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

The Phase-Locking of Tropical North Atlantic and the Contribution of ENSO

The Tropical North Atlantic (TNA) is characterized by significant interannual variability in sea surface temperature (SST), which is phase-locked to the boreal spring. In this study, the phase-locking of TNA is investigated by adopting a linear stochastic-dynamical model (SDM) using seasonally modulated TNA feedbacks together with the seasonal modulation of ENSO forcing. In the observations, the role of local TNA feedbacks and ENSO forcing in TNA phase-locking are equivalently important with both preferring the peak of TNA variability to appear in the boreal spring. Besides, the seasonal modulation of TNA feedbacks and ENSO forcing strength are both mainly controlled by thermodynamic processes. In most climate models, the contribution of ENSO on TNA phase-locking is weaker than that in observations. The strength of ENSO-related TNA phase-locking is highly correlated with the relationship between ENSO and TNA, which is mainly determined by the amplitude of ENSO and its teleconnection patterns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Mountain Basin Controls on the Snow-to-Streamflow Signal: An AIC-Weighted Multiple Linear Regression Framework

A regression-based analysis quantifies how basin characteristics modulate the snow-to-streamflow signal. First, we use the ERA5-Land reanalysis gridded product (European Centre for Medium Range Weather Forecasts reanalysis 5 -Land component) for 4,655 hydrologic unit code - 10 (HUC10) mountain basins across the western United States (US) for water years 1987–2024. Linear regressions are performed for peak snow water equivalent (SWE) and annual streamflow for each mountain basin. Models use ordinary least squares in Python’s statsmodels package. After which, an Akaike Information Criterion (AIC)–weighted ensemble multiple linear regression (MLR) framework with 47 watershed traits is used to predict the linear regression coefficient of determination (r-squared) defining the ability of peak SWE to predict annual streamflow across all mountain basin. Predictor sets are constrained to avoid multicollinearity by excluding models with variance inflation factors (VIF) greater than 5. Mountain basin traits included in the MLR include seasonal climate, topography, vegetation type and structure, and bedrock geology. Accepted models are considered if their AIC is within 2.0 of the model with the minimum AIC, or best model. To compare predictor influence across acceptable models, we computed standardized regression coefficients. To evaluate structural redundancy among models, we constructed binary inclusion vectors for each acceptable model, denoting whether a predictor was present (1) or absent (0). Core predictor variables are defined as occurring in at least 67% of the acceptable models. For this regional analysis, only one model was found acceptable, with higher snow-to-streamflow translation (higher r-squared) occurring in colder mountain basins with higher relative winter precipitation, more snow accumulation and a lower fraction of annual precipitation that falls in the spring and summer. The second component of the data package uses previously published, high-resolution output from an integrated hydrological model of the East River watershed using the U.S. Geological Survey Groundwater and Surface water Flow model (GSFLOW, doi:10.15485/1998576). East River MLR expands upon the approach described above to explore the response of five streamflow metrics—annual streamflow, runoff efficiency, 7-day minimum flow, low-flow duration, and non-perennial stream fraction to snow system indicators including peak SWE, snow-covered area, snow disappearance date, and the fraction of basin area characterized by low-to-no snow, as well as seasonal precipitation and temperature, and annual hydrologic variables representing soil moisture, evapotranspiration (ET), the partitioning of incoming precipitation to evapotranspiration (ET/P), groundwater storage, and groundwater inflow to streams. MLR was done on all water years (P0: 1987-2024) and for each period as determined in the split analysis using pooled regression techniques (P1: 1987-2011 and P2: 2012-2024) to evaluate shifting predictor variable emphasis on streamflow generation. Results indicate that since 2012, peak SWE has lost statistical strength in its prediction of annual streamflow and runoff efficiency, and the indirect influence of spring temperature has emerged as critically important. Low-flow metrics remain largely influenced by soil moisture, vegetation water use and groundwater inflows with summer precipitation becoming a direct influence on minimum summer flow. Together, these data and Python-based analysis tools provide a framework for identifying the key watershed characteristics that control how streamflow responds to snow from year to year. The package also helps quantify uncertainty in statistical models and assess how snow–streamflow relationships vary across regions and over time. This dataset contains comma-separated values files (.csv), text files (.txt), python code files (.py), figure files (.png), and shapefiles (.cpg, .dbf, .prj, .sbn, .sbx, .shp, .xml). Further details on file contents and MLR execution can be found in the readme file and the FLMD files. Work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Modeling the distribution of the endangered Jemez Mountains salamander (Plethodon neomexicanus) in relation to geology, topography, and climate

The Jemez Mountains salamander (Plethodon neomexicanus; hereafter JMS) is an endangered salamander restricted to the Jemez Mountains in north-central New Mexico, United States. This strictly terrestrial and lungless species requires moist surface conditions for activities such as mating and foraging. Threats to its current habitat include fire suppression and ensuing severe fires, changes in forest composition, habitat fragmentation, and climate change. Forest composition changes resulting from reduced fire frequency and increased tree density suggest that its current aboveground habitat does not mirror its historically successful habitat regime. However, because of its limited habitat area and underground behavior, we hypothesized that geology and topography might play a significant role in the current distribution of the salamander. We modeled the distribution of the JMS using a machine learning algorithm to assess how geology, topography, and climate variables influence its distribution. The best habitat suitability model indicates that geology type and maximum winter temperature (November to March) were most important in predicting the distribution of the salamander (23.5% and 50.3% permutation importance, respectively). Minimum winter temperature was also an important variable (21.4%), suggesting this also plays a role in salamander habitat. Our habitat suitability map reveals low uncertainty in model predictions, and we found slight discrepancies between the designated critical habitat and the most suitable areas for the JMS. Because geological features are important to its distribution, we recommend that geological and topographical data are considered, both during survey design and in the description of localities of JMS records once detected.

59 BASIC BIOLOGICAL SCIENCES↗