SEARCH · Engineering Papers
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.
This report describes the use of ComStock to model the implementation and impacts of upgrading the wall insulation, roof insulation, and windows in the U.S. commercial building stock to comply with the current state energy code in force.
Spin glasses, generally defined as disordered systems with randomized competing interactions, are a widely investigated complex system. Theoretical models describing spin glasses are broadly used in other complex systems, such as those describing brain function, error-correcting codes5 or stock-market dynamics. This wide interest in spin glasses provides strong motivation to generate an artificial spin glass within the framework of artificial spin ice systems. Here we present the experimental realization of an artificial spin glass consisting of dipolar coupled single-domain Ising-type nanomagnets arranged onto an interaction network that replicates the aspects of a Hopfield neural network. Using cryogenic X-ray photoemission electron microscopy (XPEEM), we performed temperature-dependent imaging of thermally driven moment fluctuations within these networks and observed characteristic features of a two-dimensional Ising spin glass. Specifically, the temperature dependence of the spin glass correlation function follows a power-law trend predicted from theoretical models on two-dimensional spin glasses. As a result, we observe clear signatures of the hard-to-observe rugged spin glass free energy in the form of sub-aging, out-of-equilibrium autocorrelations and a transition from stable to unstable dynamics.
A kinetic model for a solar-simulator pumped iodine laser system is developed and compared to an experiment in which the solar simulator output is dispersed over a large active volume (150 cu cm) with low simulator light intensity (approx. 200 solar constants). A trace foreign gas which quenches the upper level is introduced into the model. Furthermore, a constant representing optical absorption of the stimulated emission is introduced, in addition to a constant representing the scattering at each of the mirrors, via the optical cavity time constant. The non-uniform heating of the gas is treated as well as the pressure change as a function of time within the cavity. With these new phenomena introduced into the kinetic model, a best reasonable fit to the experimental data is found by adjusting the reaction rate coefficients within the range of known uncertainty by numerical methods giving a new bound within this range of uncertainty. The experimental parameters modeled are the lasing time, laser pulse energy, and time to laser threshold.
In recent years, record-breaking wildfire activities in the western US illustrate the need for fire mitigation efforts, such as forest fuels reduction treatments. Forests serve as crucial carbon sinks that combat the increasing effects of climate change while fuels reduction treatments may remove carbon from forested systems. As a result, forest managers need to find a balance between fire mitigation and carbon preservation. This project partnered with Boulder County Parks and Open Space (BCPOS) and the University of Colorado, Denver to investigate the 2020 Cal-Wood fire in Boulder County, Colorado. Using remote sensing data from Landsat 8 OLI, SRTM, Sentinel-2 MSI, and LiDAR, we mapped post-fire forest carbon pools and compared these values with values derived from measurements from plots on the ground. Results indicate that the correlations are R2=0.76 for aboveground live carbon, R2=0.44 for standing carbon, and R2=0.43 for aboveground dead carbon (R2=0.43), and R2=0.35 for total carbon (R2=0.35). Additionally, 38.5% of total carbon was stored in dead carbon and 14.4% was stored in live carbon. We next compared the post-fire pools between treated and untreated areas. Our analysis suggests fuels reduction treatments did not reduce carbon loss in the presence of wildfire enough to clearly distinguish the post-fire carbon in the treated and untreated areas. However, our final carbon maps still provide BCPOS and researchers with an opportunity to explore carbon estimation models based on remotely sensed data as well as a framework to evaluate fuels reduction treatment effectiveness and impact on forest carbon stocks for future wildfire events.
Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.
Climate warming is occurring fastest at high latitudes. Based on short-term field experiments, this warming is projected to stimulate soil organic matter decomposition, and promote a positive feedback to climate change. We show here that the tightly coupled, nonlinear nature of high-latitude ecosystems implies that short-term (<10 year) warming experiments produce emergent ecosystem carbon stock temperature sensitivities inconsistent with emergent multi-decadal responses. We first demonstrate that a well-tested mechanistic ecosystem model accurately represents observed carbon cycle and active layer depth responses to short-term summer warming in four diverse Alaskan sites. We then show that short-term warming manipulations do not capture the non-linear, long-term dynamics of vegetation, and thereby soil organic matter, that occur in response to thermal, hydrological, and nutrient transformations belowground. Our results demonstrate significant spatial heterogeneity in multi-decadal Arctic carbon cycle trajectories and argue for more mechanistic models to improve predictive capabilities.
This dataset contains quarterly soil core measurements from the Missouri Ozarks AmeriFlux (MOFLUX) site located at the University of Missouri’s Thomas H. Baskett Wildlife Research and Education Area near Ashland, Missouri. These data will be used to parameterize an ensemble of MOFLUX-optimized soil carbon-nitrogen models, used to simulate carbon (C) and nitrogen (N) cycling responses to future hydroclimatic scenarios and the trajectory of soil C stocks with concomitant forest decline. Beginning in 2017, eight soil cores were collected approximately quarterly near plot 1 of the southeast transect, near the automated soil respiration flux chambers, from 0–15 cm depth. Data are currently available through 2023 (2017-06-14 to 2023-11-13); additional observations will be appended to this dataset as they become available. Cores were analyzed for gravimetric moisture content, pH, total carbon and nitrogen, texture, microbial biomass carbon and nitrogen, and extractable dissolved organic carbon and nitrogen. This dataset contains one data file in comma separate (*.csv) format. Additional metadata are provided: one data dictionary and a file-level metadata file in comma separate (*.csv) format and a user guide in PDF (*.pdf) format.
Abstract Ecosystem connectivity tends to increase the resilience and function of ecosystems responding to stressors. Coastal ecosystems sequester disproportionately large amounts of carbon, but rapid exchange of water, nutrients, and sediment makes them vulnerable to sea level rise and coastal erosion. Individual components of the coastal landscape (i.e., marsh, forest, bay) have contrasting responses to sea level rise, making it difficult to forecast the response of the integrated coastal carbon sink. Here we couple a spatially-explicit geomorphic model with a point-based carbon accumulation model, and show that landscape connectivity, in-situ carbon accumulation rates, and the size of the landscape-scale coastal carbon stock all peak at intermediate sea level rise rates despite divergent responses of individual components. Progressive loss of forest biomass under increasing sea level rise leads to a shift from a system dominated by forest biomass carbon towards one dominated by marsh soil carbon that is maintained by substantial recycling of organic carbon between marshes and bays. These results suggest that climate change strengthens connectivity between adjacent coastal ecosystems, but with tradeoffs that include a shift towards more labile carbon, smaller marsh and forest extents, and the accumulation of carbon in portions of the landscape more vulnerable to sea level rise and erosion.
As part of the Communities Local Energy Action Program (CLEAP), the Lawrence Stakeholders Coalition (LSC) is interested in assessing and understanding the potential for and pathways to electrification for the City of Lawrence. The LSC's main questions are: What is the impact of various electrification packages on residential electricity bills and what types of buildings should the LSC target for electrification plus weatherization packages? This technical assistance, using ResStock tool modeling, aims to assist the Coalition's electrification and energy burden reduction planning by: 1. Providing cross-cutting data on housing stock characteristics, energy burden characteristics, fuel types, energy consumption, and system efficiency; and 2. Providing information on upgrade package costs, emissions reduction, and energy reductions by prioritized housing segment. The ResStock analysis presented here focuses on opportunities to reduce energy burden, energy consumption, and energy bills for single family homes, multifamily buildings, and mobile homes.
Despite sustained declines in Amazon deforestation, forest degradation from logging and firecontinues to threaten carbon stocks, habitat, and biodiversity in frontier forests along the Amazon arcof deforestation. Limited data on the magnitude of carbon losses and rates of carbon recoveryfollowing forest degradation have hindered carbon accounting efforts and contributed to incompletenational reporting to reduce emissions from deforestation and forest degradation (REDD+). Wecombined annual time series of Landsat imagery and high-density airborne lidar data to characterizethe variability, magnitude, and persistence of Amazon forest degradation impacts on abovegroundcarbon density (ACD) and canopy structure. On average, degraded forests contained 45.1% of thecarbon stocks in intact forests, and differences persisted even after 15 years of regrowth. Incomparison to logging, understory fires resulted in the largest and longest-lasting differences in ACD.Heterogeneity in burned forest structure varied by fire severity and frequency. Forests with a historyof one, two, and three or more fires retained only 54.4%, 25.2%, and 7.6% of intact ACD,respectively, when measured after a year of regrowth. Unlike the additive impact of successive fires,selective logging before burning did not explain additional variability in modeled ACD loss andrecovery of burned forests. Airborne lidar also provides quantitative measures of habitat structure thatcan aid the estimation of co-benefits of avoided degradation. Notably, forest carbon stocks recoveredfaster than attributes of canopy structure that are critical for biodiversity in tropical forests, includingthe abundance of tall trees. We provide the first comprehensive look-up table of emissions factors forspecific degradation pathways at standard reporting intervals in the Amazon. Estimated carbon lossand recovery trajectories provide an important foundation for assessing the long-term contributionsfrom forest degradation to regional carbon cycling and advance our understanding of the currentstate of frontier forests.
Large rainfall events over land can lead to a substantial flux of carbon and nutrients to estuaries and the coastal ocean. In the mid-Atlantic on the east coast of North America (35° - 42° N), these events often happen due to tropical storm activity as well as less predictable anomalously large midlatitude storms or abrupt spring snow melt runoff. Storms can directly impact the coastal carbon cycle via export of carbon from land to sea, while also stimulating phytoplankton production due to the influx of nutrients to the coastal ocean. To assess the impact of high precipitation and river flow events on the coastal carbon cycle, we have integrated multiple observing platforms in an analytical framework to dynamically observe carbon-related properties. The new observing system testbed (NOS-T) allows for an estimate of phytoplankton and organic carbon stocks in the surface ocean, with a goal of providing near real time analytical capability. A case study of high river discharge in the mid-Atlantic in the summer and fall of 2018 and 2021 were used to examine how riverine carbon manifests along the land-estuary-ocean continuum particularly in Chesapeake Bay. Chesapeake Bay was chosen as the study site because of long term monitoring by the Chesapeake Bay Program and a robust and well-developed regional biogeochemical modeling system, the Chesapeake Bay Environmental Forecast System (CBEFS), that is publicly available. There is also new capability in Chesapeake Bay to combine hyperspectral radiometric data from the Aerosol Robotic Network site that is online in 2022 with in-water observations of optical properties. Precipitation data from the Global Precipitation Measurement IMERG data set was used to establish triggering criteria for storm carbon flux observation. USGS and in-water carbon data were used to establish statistical models to estimate the mass flux of organic carbon into Chesapeake Bay using measured river discharge. The year 2018 was examined as there was record rainfall and near record river flow in the late summer. The CBEFS modeling was compared with in-water estimates of dissolved organic carbon (DOC) to assess the model’s ability to capture storm fluxes of carbon to the ecosystem. Multiple satellite platforms were also used to assess how remote sensing using passive Earth orbiting sensors can be used to observe carbon in this complex coastal region. A set of recommendations have been established to improve sensing capability to measure aquatic carbon during storm events. A unique challenge in these dynamic inland waters is how rapidly carbon fluxes can evolve in space and time, with many sources contributing to the water leaving reflectance that satellites can observe. To improve our ability to quantify carbon stocks and fluxes in near real time, a suite of satellite sensors and high-resolution modeling capability is needed, all supplemented by in situ monitoring. Future inclusion of in-water observations that would be deployed when the system reaches triggering criteria, as well as taskable orbital instruments, will improve estimates of ocean carbon properties and provide the ability to calculate major carbon stocks and fluxes in near real time following episodic storm events.
During normal assembly and disassembly of bolted-joint components, thread damage and/or deformation may occur. If threads are overloaded, thread damage/deformation can also be anticipated. Typical inspection techniques (e.g. using GO-NO GO gages) may not provide adequate visibility of the extent of thread damage. More detailed inspection techniques have provided actual pitch-diameter profiles of damaged-hardware holes. A method to predict the reduction in thread shear-out capacity of damaged threaded holes has been developed. This method was based on testing and analytical modeling. Test samples were machined to simulate damaged holes in the hardware of interest. Test samples containing pristine parent-holes were also manufactured from the same bar-stock material to provide baseline results for comparison purposes. After the particular parent-hole thread profile was machined into each sample a helical insert was installed into the threaded hole. These samples were tested in a specially designed fixture to determine the maximum load required to shear out the parent threads. It was determined from the pristine-hole samples that, for the specific material tested, each individual thread could resist an average load of 3980 pounds. The shear-out loads of the holes having modified pitch diameters were compared to the ultimate loads of the specimens with pristine holes. An equivalent number of missing helical coil threads was then determined based on the ratio of shear-out loads for each thread configuration. These data were compared with the results from a finite element model (FEM). The model gave insights into the ability of the thread loads to redistribute for both pristine and simulated damage configurations. In this case, it was determined that the overall potential reduction in thread load-carrying capability in the hardware of interest was equal to having up to three fewer threads in the hole that bolt threads could engage. One- half of this potential reduction was due to local pitch-diameter variations and the other half was due to overall pitch-diameter enlargement beyond Class 2 fit. This result was important in that the thread shear capacity for this particular hardware design was the limiting structural capability. The details of the method development, including the supporting testing, data reduction and analytical model results comparison will be discussed hereafter.
The outstanding mechanical and chemical properties of dental enamel emerge from its complex hierarchical architecture. An accurate, detailed multiscale model of the structure and composition of enamel is important for understanding lesion formation in tooth decay (dental caries), enamel development (amelogenesis) and associated pathologies (e.g., amelogenesis imperfecta or molar hypomineralization), and minimally invasive dentistry. Although features at length scales smaller than 100 nm (individual crystallites) and greater than 50 µm (multiple rods) are well understood, competing field of view and sampling considerations have hindered exploration of mesoscale features, i.e., at the level of single enamel rods and the interrod enamel (1 to 10 µm). Here, we combine synchrotron X-ray diffraction at submicrometer resolution, analysis of crystallite orientation distribution, and unsupervised machine learning to show that crystallographic parameters differ between rod head and rod tail/interrod enamel. This variation strongly suggests that crystallites in different microarchitectural domains also differ in their composition. Thus, we use a dilute linear model to predict the concentrations of minority ions in hydroxylapatite (Mg 2+ and CO 3 2- /Na + ) that plausibly explain the observed lattice parameter variations. While differences within samples are highly significant and of similar magnitude, absolute values and the sign of the effect for some crystallographic parameters show interindividual variation that warrants further investigation. In conclusion, by revealing additional complexity at the rod/interrod level of human enamel and leaving open the possibility of modulation across larger length scales, these results inform future investigations into mechanisms governing amelogenesis and introduce another feature to consider when modeling the mechanical and chemical performance of enamel.
A dataset of two global carbon fluxes: the loss of additional sink capacity (LASC) and present versus transient difference (PTD). The loss of additional sink capacity is the loss of indirect anthropogenic sink capacity in land ecosystems caused by the loss of ecosystems that can increase carbon in response to changing environments, primarily forests. The present transient difference refers to the direct anthropogenic flux in land ecosystems (from land cover and land use change) and the difference in that flux when calculated assuming "present day" (actually more like 1980) ecosystem carbon stocks or when calculated assuming transient carbon stocks (i.e. caused by environmental change). In this dataset, the LASC and PTD fluxes were calculated using the TRENDY ensemble of Dynamic Global Vegetation Models (DGVMs; Sitch et al., 2024) following the methods described therein and in Obermeier et al. (2021). For further details see Obermeier et al. (2021), Friedlingstein et al. (2019), Sitch et al. (2024), and Walker et al. (2025). The dataset includes ensemble means and standard deviations for annual fluxes, fluxes cumulated annually over the whole time period, and cumulated since 1959. These data were originally collected to compare different approaches (DGVMs versus bookkeeping models) for estimating land use and land cover change (LULCC) emissions (Obermeier et al. 2021). More recently, and similarly, these data were used to harmonize the estimates of the direct and indirect anthropogenic fluxes in the global carbon cycle and to recalculate the global carbon budget (Walker et al., 2025).
Our understanding of broad-scale forest disturbances under climatic extremes remains incomplete. Drought, as a typical extreme event, is a key driver of forest mortality but there have been no reports on continental-scale quantification of its impact on forest mortality or how it compares to other natural or anthropogenic drivers. Thus, our ability to understand and predict broad-scale carbon cycling in response to changing climate and extreme events is limited. In this study, we applied an attribution approach based on different sources of data to quantify the area and potential carbon loss/transfer in continental U.S. (CONUS) from four types of disturbance: (1) anthropogenic (especially timber harvest); (2) fire; (3) drought-associated; and (4) other from 2000 to 2014. Our results showed that anthropogenic disturbances, fire, drought-associated disturbances, and other disturbances accounted for 54.3, 10.7, 12.7, and 22.3% of total canopy area loss, respectively. Anthropogenic disturbance was the most important driver contributing to 58.1% potential carbon loss/transfer in CONUS for 2000–2014. The potential carbon loss/transfer from natural disturbances (fire, drought, and other) for the same study period accounted for approximately 41.9% of the total loss/transfer from all agents, suggesting that natural disturbances also played a very important role in forest carbon turnover. Potential carbon loss/transfer associated with drought accounted for approximately 11.6% of the total loss/transfer in CONUS, which was of similar magnitude to potential carbon loss/transfer from fire (~11.0%). The other natural disturbance accounted for 19.3% of potential carbon loss/transfer. Our results demonstrated the importance of the impacts of various disturbances on forest carbon stocks at the continental scale, and the drought-associated carbon loss/transfer data developed here could be used for evaluating the performance of predictive models of tree mortality under droughts.
Thanks to fast learning and sustained growth, solar photovoltaics (PV) is today a highly cost-competitive technology, ready to contribute substantially to CO 2 emissions mitigation. However, many scenarios assessing global decarbonization pathways, either based on integrated assessment models or partial-equilibrium models, fail to identify the key role that this technology could play, including far lower future PV capacity than that projected by the PV community. In this perspective, we review the factors that lie behind the historical cost reductions of solar PV and identify innovations in the pipeline that could contribute to maintaining a high learning rate. We also aim at opening a constructive discussion among PV experts, modelers, and policymakers regarding how to improve the representation of this technology in the models and how to ensure that manufacturing and installation of solar PV- can ramp up on time, which will be crucial to remain in a decarbonization path compatible with the Paris Agreement.
Increasing the amount of soil organic carbon (SOC) has agronomic benefits and the potential to mitigate climate change. Previous regional predictions of SOC trends under climate change often ignore or do not explicitly consider the effect of crop adaptation (i.e., changing planting dates and varieties). We used the DayCent biogeochemical model to examine the effect of adaptation on SOC for corn and soybean production in the U.S. Corn Belt using climate data from three models. Without adaptation, yields of both corn and soybean tended to decrease and the decomposition of SOC tended to increase leading to a loss of SOC with climate change compared to a baseline scenario with no climate change. With adaptation, the model predicted a substantially higher crop yield. The increase in yields and associated carbon input to the SOC pool counteracted the increased decomposition in the adaptation scenarios, leading to similar SOC stocks under different climate change scenarios. Consequently, we found that crop management adaptation to changing climatic conditions strengthen agroecosystem resistance to SOC loss. However, there are differences spatially in SOC trends. The northern part of the region is likely to gain SOC while the southern part of the region is predicted to lose SOC.
The characterstics of NITRONIC 40 were investigated in connection with its use in cryogenic wind tunnel models. In particular, the effects of carbide and sigma-phase precipitation resulting from heat treatment and the presence of delta ferrite were evaluated in relation to their effects on mechanical properties and the potential consequences of such degradation. Methods were examined for desensitizing the material and for possible removal of delta ferrite as a means of restoring the material to its advertised properties. It was found that heat treatment followed by cryogenic quenching is a technique capable of desensitizing NITRONIC 40. However, it was concluded that it is extremely difficult, if not impossible, to remove the delta ferrite from the existing stock of material. Furthermore, heat treatments for removing delta ferrite have to take place at temperatures that cause very large grain growth. The implications of using the degraded NITRONIC 40 material for cryogenic model testing were reviewed, and recommendations were submitted with regard to the acceptability of the material. The experience gained from the study of NITRONIC 40 clearly identifies the need to implement a policy for purchasing top-quality materials for cryogenic wind tunnel model applications.