Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “depth”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2022

A Numerical Study of the Liquid Motion in Titan’s Subsurface Ocean

An ocean filled with liquid water lies beneath the icy surface of several Jovian and Saturnian moons. In such an ocean, the currents are driven by various phenomena such as the tidal forcing, the deformation of the ice shell lying at its top, the temperature gradient resulting from the surface and bottom heat fluxes…The flow induced by the first two forcings can be modelled by means of a 2D depth-averaged model, while the third one generates horizontal and vertical density variations whose effects can only be captured by a 3D baroclinic model. We study the tides of Titan’s subsurface ocean and the impact of the ice shell on the liquid motion by means of the Second-generation Louvain-la-Neuve Ice–ocean Model, SLIM (https://www.climate.be/slim). The impact of the ice shell lying at the top of the ocean is modelled by a surface friction term and surface pressure terms. The latter are a function of the difference between the ocean elevation and the vertical displacement of the shell and the time derivative of this difference. Because of Titan’s appreciable obliquity (0.306°), the tidal motion expected (and found) is similar to the Europa tidal scenario described by Tyler (2008): the surface elevation consists of two bulges rotating around Titan and the associated depth-averaged velocity field consists of two gyres, separated by an area of high speed flow, whose centre follows a sinusoidal path centred on the equator. The ice shell damps the surface motion, thus slowing down the flow, without significantly modifying the spatial patterns of these fields. The depth of the ocean and the mechanical characteristics of the ice shell being poorly constrained, a sensitivity analysis is conducted. The depth-averaged flow slows down when the depth is increased and a lag appears in the tidal phase but the tidal range remains similar. The ice shell mechanical characteristics influences both the elevation and depth-averaged velocity fields in terms of magnitude but does not modify the spatial patterns of these fields. The influence of the surface heat flux is studied by means of the 3D baroclinic version of SLIM. The heat flux derived from Titan’s topography by Kvorka et al. (2018) is used as surface boundary condition for the temperature equation while a uniform bottom heat flux is implemented. Its value is computed assuming that the heat budget of the ocean is at equilibrium. These boundary conditions cause density variations, which impact the hydrodynamics of the ocean. While the flow velocity induced by these variations is two orders of magnitude smaller than the tidal flow, its orientation is time-independent, hence impacting the orientation of the velocity field. Although the variations of ocean surface elevation and speed with respect to the shell mechanical properties can be larger than those induced by the surface heat flux, taking into account the latter results in large variations of the velocity field global patterns, which was not observed when modifying the shell mechanical properties. Future studies should therefore focus on modelling the surface and bottom heat fluxes while uncertainties about the mechanical characteristics of the shell can be tolerated.

Titan↗
NASA NTRSDate not supplied

TRIDENT Drill Validation at Mars and Lunar Analog Field Sites

Drilling on Earth is typically a human-intensive activity. Drilling on other planets is further complicated by the lack of prior local field surveys of their target area, hence blindly drilling into uncertain target rocks. Field conditions on the Moon or Mars are also different than for shallow drilling on Earth: lower temperatures and pressures, less power available, low masses (hence less weight-on-bit). Given the cost of transport from Earth, no drilling muds or working fluids are likely to be available to carry away cuttings. And impact-gardened regolith and dust vary mechanically and texturally from most terrestrial soils. The Regolith and Ice Drill for Exploration of New Terrains (TRIDENT) is a rotary-percussive 1m-class drill from Honeybee Robotics. It is low-power (rotary and percussive actuators are 200 W each) and lightweight (<20 kg) with the maximum weight on bit limited to 200 N. TRIDENT has been manifested for the Volatiles Investigating Polar Exploration Rover (VIPER) and PRIME-1 lunar south pole missions in 2024, has previously been field tested at a hot, dry analog site in the Atacama Desert, and in lunar conditions in thermal vacuum chamber tests. TRIDENT was also part of the 2019 Icebreaker Mars Discovery proposal, as well as in the Mars Life Explorer concept. During ARADS tests, drill control and fault recovery automation software enabled hands-off operations of a rover-mounted TRIDENT drill. TRIDENT Drill Analog Site Validation: Past TRIDENT tests in thermal vacuum (TVAC) chambers targeted containers of manufactured lunar simulants with added volatiles. 2022 TRIDENT ambient testing at NASA Ames drilled into cemented lunar simulant materials. Low cuttings-permeability led to cuttings buildup, and drill choking and binding was observed. The Atacama analog site in ARADS had desiccated unconsolidated sediments that did not challenge the TRIDENT design. However, lunar polar regolith is expected to be diverse and heterogeneous with varying clast sizes, with abundant impactites and perhaps subsurface ice deposits. Neither the simulants nor Atacama testing had completely covered the TRIDENT-targeted field characteristics, motivating further analog tests prior to the planned lunar missions. To gain more insight into the behavior of the TRIDENT hardware in diverse impactites and subsurface ice, and to verify the software automation in that environment, in August 2023 TRIDENT was brought to Haughton Crater, a field analog site in the Canadian Arctic. In September 2023 the same drill was brought to the Bishop Tuff in southern California to verify whether drilling binding behaviors previously seen in lunar simulant testing would be observed in naturally occurring fine-grained massive layers. The ~22 Ma Haughton Crater impact structure is located at 75 ̊22’ N, 89 ̊41’ W, on northwestern Devon Island, Nunavut, Canada. Numerous deposits of pale-grey crater-fill polymictic impact-melt breccia are found within the crater with a typical thickness reaching ~125 m or greater and covering ~60 km -2 . An approximately 600m-thick permafrost layer is also present with ice typically found within 0.5-0.6m of the surface. The volcanic tableland north of Bishop, CA exposes densely welded tuff laid down during the eruption that created the Long Valley Caldera at approximately 0.76 Ma. Extensional faults and the Owens River gorge expose cross-sections across the plateau. The area is viewed as an analog site for Mars features believed to be of pyroclastic origin. Results: Haughton Crater.Drilling tests were conducted 8-13 August 2023 at a previously undisturbed area separated 5-10 m from past years’ Drill Hill test sites (75.4208, -089.7613). In six days, TRIDENT drilled 8 holes to nearly 1 m depth each, totaling 7.80 m. The active layer/ice boundary was at ~67 cm depth, with a total of approximately 2.4 m drilled into ice or ice-cemented impact breccia. During drilling, five drill fault states were observed and successfully recovered. Holes 23-1, 4 and 7 were drilled under manual control, using Honeybee’s Thorax user interface. Holes 23-5, 6, and 8 were drilled with the Ames IBexec automated drilling control software. TRIDENT was observed to have little difficulty in the thawed uncemented impact breccia above the active layer boundary, but required percussion to make slower headway in the ice-cemented breccia. In Hole 23-7 (Fig.1), drilling slowed down in a massive unit just above the active-layer boundary (perhaps a large rock extending into the ice-cementation?), with only 7cm progress made in 27 minutes of high auger torque and constant percussion, leading to a choking fault and then a binding fault. A similar pattern had been observed in TRIDENT Rio Tinto test data from 2017 [6] as well as in the 2022laboratory tests.Bishop Tuff.A team from NASA Ames and the US Geological Survey deployed the same TRIDENT drill to Bishop Sites 1B and 1C (37.4203, -118.4289; 37.4265, -118.4215) on 13-16 September 2023, on the Bishop Tuff plateau. A third drill site was used 17-18 September 2023(37.4598, -118.3667) in an abandoned pumice mine. Four holes (totaling 2.5m depth) were drilled into the fine-grained, meters-thick tuff units at Sites 1B and 1C, and a further two boreholes (totaling 1.98m depth) were made at the pumice site. Drill behavior in the tuff below 10 cm depth was similar to that seen at 65-74 cm depth in Haughton Hole 23-7 (Fig. 1) and that seen in the 2022 lab simulant drilling. Drill safety torque limits were exceeded multiple times resulting in drill stops downhole. These freezes then required external added torques (with a pipe wrench) to resume rotation, to unstick the drill for withdrawal. To prevent this choking and binding behavior we found that more-frequent cuttings removal was necessary, e.g., reducing the “drill bite” size from the nominal 10 cm to 2 cm per bite --bringing the auger up to the surface more frequently, as seen in Bishop Site 1C Hole 2 (Fig.2). This permitted slow progress without drill binding and without external interventions. Conversely, TRIDENT drilling in the more porous pumice target material showed no cuttings buildup issue, and single bites as large as 40 cm were demonstrated. Discussion: We observed that TRIDENT easily penetrated unconsolidated heterogeneous soils (both above the active layer boundary at Haughton and previously in the Atacama). Cemented or consolidated targets that were cuttings-permeable (icy impact breccia, pumice) required more energy applied and percussion. However, in non-cuttings-permeable targets (welded microporous tuff, cemented simulants, boulder) TRIDENT was observed to be prone to excessive cuttings accumulation leading to choking/binding faults and stalling. The wedge cutting bit, used by TRIDENT in field tests and in its flight versions, pulverizes the target rock and creates fine cuttings that ideally are transported up the auger spirals for removal. In porous, fractured or vesicular target materials (such as at the Bishop pumice site) a significant portion of the cuttings are pushed aside, but for non-fractured, microporous targets the cuttings remain in the borehole and accumulate. Rock powder is relatively incompressible as a working fluid at only 100-200N downward force (TRIDENT limits) and hence eventual drilling progress slows or stops. Our recommended strategy for improving TRIDENT cuttings removal in massive target units with low cuttings-permeability is to reduce TRIDENT bite sizes when encountering these units, from 10cm to as little as 1-2cm, to effectively bail the accumulating cuttings. This approach was demonstrated to reduce choking and allowed slow progress to continue in cuttings-impermeable microporous target units (viz. the Bishop tuff in our September 2023 tests or cemented simulants in 2022 ambient tests).

robotic drilling↗
DOE OSTI2020

Borehole Disposal of Radioactive Waste in Israel - 20388

Israel is assessing borehole disposal of radioactive waste. With limited geological options for disposal, intermediate-depth borehole disposal is being considered in the arid Yamin Plain region of the northeastern Negev desert at depths of several hundred meters below ground surface in the vadose zone. Unlike deep borehole disposal of several kilometers, which relies on emplacement below the depth of recirculating groundwater, the safety case for intermediate-depth borehole disposal relies more on the aridity of the vadose zone, the robustness of the waste package, and the other engineered barriers (e.g., seals, backfill materials) in the disposal borehole. As a tradeoff, the shallower depth may allow the use of a larger diameter borehole than would be possible for a deeper borehole. To study the suitability of the Yamin Plain region for borehole disposal, a small-diameter characterization borehole is being planned to retrieve core samples and to better understand the vadose zone geomechanical and hydrogeochemical properties and percolation flux. The information from the characterization borehole will inform the safety case, and together with performance assessment analyses will help to identify key areas of uncertainty and guide future research and development activities aimed at demonstrating the feasibility of the intermediate-depth borehole disposal concept in Israel. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Heating effects on pyrogenic organic matter mass loss and temperature in a sand matrix in a lab setting from 2022

This dataset contains data associated with the preprint “A laboratory method for simulating the effects of subsequent fires on pyrogenic organic matter at different exposure depths in a sand matrix” (https://doi.org/10.1101/2024.07.30.605814), which has been accepted for publication in the International Journal of Wildland Fire. The dataset was generated to investigate the effects of subsequent fires on preexisting pyrogenic organic matter (PyOM) at different soil depths. Specifically, we designed a laboratory method to simulate soil heating from above, applying realistic heat flux profiles (High, Low, Control) derived from jack pine (Pinus banksiana Lamb.) log burns. Using a cone calorimeter, these profiles were applied to buried jack pine PyOM (produced at 350°C) in a sand matrix pyrocosm to examine PyOM mass loss at different depths (surface, 1 cm, 5 cm). Each treatment (e.g., High heat flux at 1 cm depth) was applied once to five replicates in a pyrocosm.The dataset includes mass loss fraction of PyOM and time-series temperature profiles (measured with type K thermocouples). Additionally, calibration datasets for high and low heat flux profiles are provided, including programmed heat flux profile data and three sets of calibration results. All datasets are in .csv format and were analyzed in R using code found at https://github.com/MengmengLuo/Reburning-pyrogenic-organic-matter-A-laboratory-method-for-dosing-dynamic-heat-fluxes-from-above. This dataset supports research on PyOM transformations during subsequent fires, and methodological advancements in simulating soil heating. It can be used to explore fire-induced PyOM alterations across different fire intensities and soil depths.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

How deep should we go to understand roots at the top of the world?

Harsh environmental conditions and the short summers of northern, high-latitude biomes impose unique constraints on the plants that live in the arctic tundra and the boreal forest. To escape the harsh aboveground environment, plants in these habitats often allocate a large portion of their biomass belowground to facilitate nutrient acquisition. In turn, the proximity of living plant roots to vast stores of sequestered soil carbon in these biomes means that shifts in rooting depth distribution and the size of the root–soil interface could significantly contribute to ongoing climate change. Indeed, plant ‘priming’ of rhizosphere decomposition via root exudation, particularly from shallowly distributed roots, can lead to losses of carbon from tundra soils. While we have a hard-won understanding of the distribution of plant communities across the arctic tundra and the boreal forest from direct field observations scaled to the landscape level using climate-informed mapping techniques (i.e. the Circumpolar Arctic Vegetation Map (CAVM); Walker et al. , 2005), these vegetation maps are only the tip of the iceberg (Iversen et al. , 2015). Root form and function remain hidden beneath the land surface. In an article recently published in New Phytologist , Blume-Werry et al. (2023, 10.1111/nph.18998) asked whether rooting depth distribution, and ensuing carbon emissions, could be inferred from commonly used vegetation mapping classifications across the pan-Arctic. An important question to guide our understanding, mapping, and prediction of belowground characteristics and ecosystem feedbacks at the top of the world. Unfortunately, they found that the answer was ‘not quite’. While rooting depth distribution varied demonstrably, in turn causing substantial changes in modeled carbon emissions via rhizosphere ‘priming’, variation across rooting depth profiles did not correspond with vegetation mapping classes. If we are unable to predict belowground rooting depth distributions across large spatial scales by leveraging aboveground vegetation community distributions, how then should belowground researchers proceed?

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Measuring and modeling soil moisture and runoff at solar farms using a disconnected impervious surface approach

Abstract Ground‐mounted photovoltaic sites are often treated as impervious surfaces in stormwater permits. This ignores the pervious soils beneath and between solar arrays and leads to an overestimation of runoff. Our objective was to improve solar farm stormwater hydrology models by explicitly considering the disconnected impervious nature of solar design and site characteristics. Experimental sites established on utility scale solar farms in Colorado, Georgia, Minnesota, New York, and Oregon had perennial vegetative plantings with mean precipitation ranging from 40.6 to 124.5 cm, and soil texture ranging from loamy sand to clay. Soil moisture measurements were collected beneath arrays, under drip edges, and in the vegetated area between arrays at each site. Hydrus‐3D models for soil moisture and stormwater hydrology were developed that accounted for precipitation falling on solar panels, drip edge redistribution of rainfall, infiltration, and runoff in the pervious areas between solar arrays and beneath panels. Drip edge runoff averaged 3‐ to 10‐times incident precipitation at the New York and Minnesota sites, respectively. Root mean square error values between measured sub‐hourly soil moisture and predicted moisture for large measured single storm events averaged 0.029 across all five sites. Predicted runoff depths were strongly affected by precipitation depth, soil texture, soil profile depth, and soil bulk density. Runoff depths across the five experimental sites averaged 13%, 25%, and 45% of the 2‐, 10‐, and 100‐year design storm depths, clearly showing that these solar farms do not behave like impervious surfaces, but rather as disconnected impervious surfaces with substantial infiltration of runoff in the vegetated areas between and beneath solar arrays.

Agriculture↗
DOE OSTI2020

Crack nucleation and growth during dynamic indentation of chemically-strengthened glass

Dynamic point impact loading is a primary cause of fracture of glass screens on mobile devices. An improved understanding of crack initiation and evolution under high-speed indentation could contribute to the development of materials with better performance, but experimental observations are challenging due to the short timescales and limited depth-of-field of optical microscopy. To address this need, we have observed fracture of a chemically-strengthened glass during dynamic indentation using in situ x-ray phase-contrast imaging (XPCI). Median cracks initiate below the surface of the glass, at a depth approximately corresponding to the depth at which the surface residual compressive stress diminishes to zero. These cracks initially propagate at ~10ms –1 , rapidly accelerating to > 100ms –1 . We also observe some evidence for rate-dependent behavior, in that indentation at the lowest rates studied here (below about ~0.15ms –1 ) fails to initiate cracks regardless of the depth of the indentation (up to 18μm), while indentation at higher rates produces median cracks that either arrest or cause complete fracture, depending on the depth of indentation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Forecasting Multi-Step-Ahead Street-Scale Nuisance Flooding using a seq2seq LSTM Surrogate Model for Real-Time Application in a Coastal-Urban City

In coastal-urban cities facing an elevated risk of nuisance flooding (by rain and tide) due to increased heavy rainfall, sea level rise, urbanization, and aging drainage systems, real-time flood forecasting at the street-scale can provide useful information to transportation decision-makers. Physics-Based Models (PBMs) that offer high accuracy come with high computational runtimes and costs that limit their application for real-time flood forecasting. To address this challenge, Machine Learning (ML) surrogate models trained from PBMs have been proposed to provide street-scale flood forecasts. Previous related studies have focused on using Long Short-Term Memory (LSTM) architectures to model hourly flood depth on streets. While LSTM models can capture input sequences effectively, they fall short in accurately preserving output sequences, limiting their suitability for multi-step-ahead forecasts. The seq2seq LSTM architecture offers a key advantage here by capturing the full sequence of input–output, making it potentially more suitable for multi-step-ahead flood forecasts compared to traditional LSTM models. However, seq2seq LSTM has not been tested for street-scale flood forecasting, particularly for rapidly fluctuating nuisance flooding events which require special attention to its temporal sequences. Hence, in this study, we applied the seq2seq LSTM model to explore multi-step-ahead street-scale nuisance flooding and compared its results to the traditional LSTM model as a benchmark model. LSTM and seq2seq LSTM surrogate models were applied to 22 flood-prone streets in Norfolk, Virginia, as a case study with a 4-hr (short-term) and 8-hr (long-term) lead time. The models were trained with environmental (rainfall and tide) and topographic (elevation, Topographic Wetness Index, and Depth-To-Water) features along with PBM-derived water depths for different storm events. The results demonstrated satisfactory performance of both LSTM and seq2seq LSTM surrogate models throughout the forecast period compared to the PBM. However, the seq2seq LSTM showed lower Mean Absolute Error (MAE)/ Root Mean Square Error (RMSE) and higher Nash–Sutcliffe Efficiency (NSE)/ correlation than the LSTM across most lead times, particularly for long-term forecasting due to its supremacy in handling both input–output sequences together, which is missing in the traditional LSTM. For example, in the long-term, the average RMSE ranges were 0.0268–0.0373 m for LSTM and 0.0226–0.0319 m for seq2seq LSTM, while in the short-term, they were 0.0263–0.0293 m and 0.0261–0.0283 m, respectively. Additionally, while both models exhibited similar performance in distinguishing flooded and non-flooded streets for flood depth ≥ 0.1 m, the seq2seq LSTM model demonstrated superior performance for higher flood depths (such as ≥ 0.2 m and ≥ 0.3 m). Once trained, inference took only 0.09 to 0.11 s (short-term) and 0.30 to 0.35 s (long-term) per storm event for the 22 streets, making the application highly suitable for real-time decision-making during nuisance flood events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

The nanoscale distribution of copper and its influence on charge collection in CdTe solar cells

For decades, Cu has been the primary dopant in CdTe photovoltaic absorbers. Typically, Cu acceptor concentrations in these devices are on the order of 1 1014 cm-3, which has made it notoriously difficult to directly correlate nanoscale Cu distributions to the local charge transport properties of these devices. To measure and correlate these properties, measurement techniques require high sensitivity to elemental concentration, large penetration depth, and operando compatibility. Techniques such as secondary-ion mass spectroscopy and X-ray energy dispersive spectroscopy are widely adopted to measure Cu concentrations, but they are limited by penetration depth, sensitivity, or spatial resolution. Additionally, they lack the operando capabilities required to correlate one-to-one Cu concentrations to electrical performance. In this work, correlative X-ray microscopy is used to investigate the spatial distribution of Cu and its impact on charge collection through the depth and breadth of CdTe photovoltaic devices. Plan-view, nanoscale X-ray fluorescence maps clearly demonstrate the spatial segregation of copper around regions thought to be CdTe grain boundaries. Complementary cross-section imaging unveils the transition of the maximum charge-collection efficiency from the ZnTe-CdTe interface to the CdS-CdTe interface as a function of Cu incorporation. The copper concentration through the depth of the CdTe layer is characterized by slow and fast diffusion components, and cross-section charge-transport modeling shows that the experimentally obtained charge collection can be explained by the modeled acceptor distribution through the depth of the CdTe layer.

14 SOLAR ENERGY↗
DOE OSTI2023

Collimated phase measuring deflectometry

Phase measuring deflectometry has been applied for free-form specular surface metrology, but its measured slope results are sensitive to the depth of sample positioning, which is also called the height-slope ambiguity. The objective of this work is to tackle this height-slope ambiguity problem. The main idea is to introduce collimated camera rays using a telecentric imaging lens and collimated structured-light illumination with a Fourier lens. This setup makes the fringe phases become only sensitive to the surface slopes and insensitive to the depth of the sample positioning. In this way, the slope calculation is theoretically independent of the sample depth. We call this new deflectometry technique Collimated Phase Measuring Deflectometry (CPMD). With our developed CPMD experimental setup, the measurement is insensitive to the depth of sample positioning, e.g., the measured height dispersion is less than 30 nm RMS within a 10 mm depth range when measuring a 50-mm-diameter spherical mirror with a 200 mm radius of curvature. In conclusion, the merits and limitations of the proposed CPMD technique are discussed, revealing its prospects in practical metrology applications and potential future investigations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Organic colloid composition in variable-redox porewaters within a mountainous floodplain

Redox gradients, often driven by changes in sediment moisture levels in porous, heterogeneous groundwater systems, create dynamic conditions that may promote the production and transport of colloids within natural waters. While much research has focused on the inorganic composition of colloids, the organic composition remains less well understood. Organic matter (OM) in colloids may associate with minerals, complex metal ions, and serve as an electron donor for microbial respiration; therefore, its composition is of high interest. We examined the composition of porewater OM along a redox gradient in a riparian soil located along the Slate River in Crested Butte, Colorado, USA as a function of depth (90, 130, 200, and 350 cm below ground surface). All depths were oxic to suboxic, except 200 cm, where the products of iron and sulfate reduction were observed concomitant with an increase in dissolved and/or colloidal OM, pH, alkalinity, and conductivity. We investigated the composition of porewater using correlated scanning transmission X-ray microscopy and transmission electron microscopy. We observed a change in the OM chemistry from carboxylate-rich at the 200 cm depth (where it was also enmeshed with non-crystalline iron) to phenol- and substituted-aromatic-rich at other depths. Radiocarbon dating revealed carbon in porewater at 200 cm was younger than depths above and below. Soil porewater can flow down into the underlaying gravel bed during baseflow conditions, thus we speculate whether riparian porewater could transport iron- and carboxylate-rich organic colloids into surrounding surface waters through the gravel bed conduit.

Alpine floodplain↗
DOE OSTI2025

Coalescence-Induced Spontaneous Shedding of Microdroplets on Superhydrophobic Surfaces Featuring Enclosed Micropillars with Hierarchical Roughness

This study investigated water vapor condensation on superhydrophobic surfaces (SHSs) featuring micropillars enclosed by wall lattices and having three-tier hierarchical roughness. A total of five samples were created with three (NW-J, W200-J and W400-J samples) having large micropillar depth (~6 μm) and two (NW-S and W200-S samples) having small micropillar depth (~1 μm). Two distinct condensate removal modes were observed during condensation: coalescence-induced jumping on samples with large micropillar depth and coalescence-induced shedding on samples with small micropillar depth. The results showed that the diameter of the shedding droplet on the W200-S sample having small micropillar depth could be as small as 107 μm, as compared to the theoretical critical diameter of 267 μm for gravitational shedding on the same sample. The enhanced functionality of the three-tier nanotextures on the W200-S sample could effectively suppress localized pinning of the three-phase contact line and Wenzel neck formation during the growth of condensate droplets. Consequently, during multidroplet coalescence, the released surface energy easily overcomes the solid–liquid adhesion, leading to spontaneous shedding of merged droplets. The inclusion of the wall lattice aids condensate growth by the droplet self-alignment along the walls and promoting coalescence. As a result, the W200-S sample exhibited the highest condensate collection as well. In conclusion, the proposed surface design has great potential for scaling up and implementation in heating, ventilation, and air-conditioning equipment due to the simplicity of the surface morphology and the facile spray-coating method used to achieve hierarchical roughness.

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

Microtopography Matters: Belowground CH 4 Cycling Regulated by Differing Microbial Processes in Peatland Hummocks and Lawns

Water table depth and vegetation are key controls of methane (CH 4 ) emissions from peatlands. Microtopography integrates these factors into features called microforms. Microforms often differ in CH 4 emissions, but microform-dependent patterns of belowground CH 4 cycling remain less clearly resolved. Here, to investigate the impact of microtopography on belowground CH 4 cycling, we characterized depth profiles of the community composition and activity of CH 4 -cycling microbes using 16S rRNA amplicon sequencing, incubations, and measurements of porewater CH 4 concentration and isotopic composition from hummocks and lawns at Sallie's Fen in NH, USA. Geochemical proxies of methanogenesis and methanotrophy indicated that microforms differ in dominant microbial CH 4 cycling processes. Hummocks, where water table depth is lower, had higher porewater redox potential (Eh) and higher porewater δ 13 C-CH 4 values in the upper 30 cm than lawns, where water table depth is closer to the peat surface. Porewater δ 13 C-CH 4 and δD-CH 3 D values were highest at the surface of hummocks where the ratio of methanotrophs to methanogens was also greatest. These results suggest that belowground CH 4 cycling in hummocks is more strongly regulated by methanotrophy, while in lawns methanogenesis is more dominant. We also investigated controls of porewater CH 4 chemistry. The ratio of the relative abundance of methanotrophs to methanogens was the strongest predictor of porewater CH 4 concentration and δ 13 C-CH 4 , while vegetation composition had minimal influence. As microbial community composition was strongly influenced by redox conditions but not vegetation, we conclude that water table depth is a stronger control of belowground CH 4 cycling across microforms than vegetation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A simple soil mass correction for a more accurate determination of soil carbon stock changes

Abstract Agricultural soils can act as a sink for large quantities of soil organic carbon (SOC) but can also be sources of carbon to the atmosphere. The international standard for assessing SOC stock and measuring stock change stipulates fixed depth sampling to at least 30 cm. The tendency of bulk density (BD) to decrease with decreasing disturbance and increasing SOC concentration and the assumption of constant SOC and BD within this depth profile promotes error in the estimates of SOC stock. A hypothetical but realistic change in BD from 1.5 to 1.1 g cm −3 from successive fixed depth sampling to 30 cm underestimates SOC stock change by 17%. Significant effort has been made to evaluate and reduce this fixed depth error by using the equivalent soil mass (ESM) approach, but with limited adoption. We evaluate the error in SOC stock assessment and change generated from fixed depth measurements over time relative to the ESM approach and propose a correction that can be readily adopted under current sampling and analytical methods. Our approach provides a more accurate estimate of SOC stock accumulation or loss that will help incentivize management practice changes that reduce the environmental impacts of agriculture and further legitimize the accounting practices used by the emerging carbon market and organizations that have pledged to reduce their supply chain greenhouse gas (GHG) footprints.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Equation of state for He bubbles in W and model of He bubble growth and bursting near W{100} surfaces derived from molecular dynamics simulations

Abstract Molecular dynamics (MD) simulations are performed to derive an equation of state (EOS) for helium (He) bubbles in tungsten (W) and to study the growth of He bubbles under a W(100) surface until they burst. We study the growth as a function of the initial nucleation depth of the bubbles. During growth, successive loop-punching events are observed, accompanied by shifts in the depth of the bubble towards the surface. Subsequently, the MD data are used to derive models that describe the conditions that cause the loop punching and bursting events. Simulations have been performed at 500, 933, 1500, 2000, and 2500 K to fit the parameters in the models. To compute the pressure in the bubble at the loop punching and bursting events from the models, we derive an EOS for He bubbles in tungsten with an accompanying volume model to compute the bubble volume for a given number of vacancies ( $$N_\text {V}$$ N V ), He atoms ( $$N_\text {He}$$ N He ), and temperature ( T ). To derive the bubble EOS, we firstly derive the EOS for a free He gas. The derived free-gas EOS can accurately predict all MD data included in the analysis (which span up to 54 GPa at 2500 K). Subsequently, the bubble EOS is derived based on the free-gas EOS by correcting the gas density to account for the interaction between He and W atoms. The EOS for the bubbles is fitted to data from MD simulations of He bubbles in bulk W that span a wide range of gas density and sizes up to about 3 nm in diameter. The pressure of subsurface bubbles at the loop punching events as calculated using the bubble-EOS and the volume model agrees well with the pressure obtained directly from the MD simulations. In the loop punching model, for bubbles consisting of $$N_\text {V}$$ N V vacancies and $$N_\text {He}$$ N He helium atoms, the $$N_\text {He}/N_\text {V}$$ N He / N V ratio that causes the event, the resulting increase in $$N_\text {V}$$ N V , and the associated shift of the bubble depth are formulated as a function of $$N_\text {V}$$ N V and T . In the bursting model, a bubble must simultaneously reach a certain depth and $$N_\text {He}/N_\text {V}$$ N He / N V ratio in order to burst. The burst depth and $$N_\text {He}/N_\text {V}$$ N He / N V are also modeled as a function of $$N_\text {V}$$ N V and T . The majority of the loop punching events occur at bubble pressures between 20 and 60 GPa, depending on the bubble size and temperature. The larger the bubble and the higher the temperature, the lower the bubble pressure. Furthermore, our results indicate that at a higher temperature, a bubble can burst from a deeper region.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Local-scale heterogeneity of soil thermal dynamics and controlling factors in a discontinuous permafrost region

In permafrost regions, the strong spatial and temporal variability in soil temperature cannot be explained by the weather forcing only. Understanding the local heterogeneity of soil thermal dynamics and their controls is essential to understand how permafrost systems respond to climate change and to develop process-based models or remote sensing products for predicting soil temperature. In this study, we analyzed soil temperature dynamics and their controls in a discontinuous permafrost region on the Seward Peninsula, Alaska. We acquired one-year temperature time series at multiple depths (at 5 or 10 cm intervals up to 85 cm depth) at 45 discrete locations across a 2.3 km 2 watershed. We observed a larger spatial variability in winter temperatures than that in summer temperatures at all depths, with the former controlling most of the spatial variability in mean annual temperatures. We also observed a strong correlation between mean annual ground temperature at a depth of 85 cm and mean annual or winter season ground surface temperature across the 45 locations. We demonstrate that soils classified as cold, intermediate, or warm using hierarchical clustering of full-year temperature data closely match their co-located vegetation (graminoid tundra, dwarf shrub tundra, and tall shrub tundra, respectively). We show that the spatial heterogeneity in soil temperature is primarily driven by spatial heterogeneity in snow cover, which induces variable winter insulation and soil thermal diffusivity. These effects further extend to the subsequent summer by causing variable latent heat exchanges. Finally, we discuss the challenges of predicting soil temperatures from snow depth and vegetation height alone by considering the complexity observed in the field data and reproduced in a model sensitivity analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Feature learning and generalization in deep networks with orthogonal weights

Fully-connected deep neural networks with weights initialized from independent Gaussian distributions can be tuned to criticality, which prevents the exponential growth or decay of signals propagating through the network. However, such networks still exhibit fluctuations that grow linearly with the depth of the network, which may impair the training of networks with width comparable to depth. We show analytically that rectangular networks with tanh activations and weights initialized from the ensemble of orthogonal matrices have corresponding preactivation fluctuations which are independent of depth, to leading order in inverse width. Moreover, we demonstrate numerically that, at initialization, all correlators involving the neural tangent kernel (NTK) and its descendants at leading order in inverse width—which govern the evolution of observables during training—saturate at a depth of ~20, rather than growing without bound as in the case of Gaussian initializations. We speculate that this structure preserves finite-width feature learning while reducing overall noise, thus improving both generalization and training speed in deep networks with depth comparable to width. We provide some experimental justification by relating empirical measurements of the NTK to the superior performance of deep non-linear orthogonal networks trained under full-batch gradient descent on the MNIST and CIFAR-10 classification tasks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

PRISM3D: a 3-D reference seismic model for Iberia and adjacent areas

SUMMARY We present PRISM3D, a 3-D reference seismic model of P- and S-wave velocities for Iberia and adjacent areas. PRISM3D results from the combination of the most up-to-date earth models available for the region. It extends horizontally from 15°W to 5°E in longitude, 34°N to 46°N in latitude and vertically from 3.5 km above to 200 km below sea level, and is modelled on a regular grid with 10 and 0.5 km of grid node spacing in the horizontal and vertical directions, respectively. It was designed using models inferred from local and teleseismic body-wave tomography, earthquake and ambient noise surface wave tomography, receiver function analysis and active source experiments. It includes two interfaces, namely the topography/bathymetry and the Mohorovičić (Moho) discontinuity. The Moho was modelled from previously published receiver function analysis and deep seismic sounding results. To that end we used a probabilistic surface reconstruction algorithm that allowed to extract the mean of the Moho depth surface along with its associated standard deviation, which provides a depth uncertainty estimate. The Moho depth model is in good agreement with previously published models, although it presents slightly sharper gardients in orogenic areas such as the Pyrenees or the Betic-Rif system. Crustal and mantle P- and S-wave wave speed grids were built separately on each side of the Moho depth surface by weighted average of existing models, thus allowing to realistically render the speed gradients across that interface. The associated weighted standard deviation was also calculated, which provides an uncertainty estimation on the average wave speed values at any point of the grid. At shallow depths (<10 km), low P and S wave speeds and high VP/VS are observed in offshore basins, while the Iberian Massif, which covers a large part of western Iberia, appears characterized by a rather flat Moho, higher than average VP and VS and low VP/VS. Conversely, the Betic-Rif system seems to be associated with low VP and VS, combined with high VP/VS in comparison to the rest of the study area. The most prominent feature of the mantle is the well known high wave speed anomaly related to the Alboran slab imaged in various mantle tomography studies. The consistency of PRISM3D with previous work is verified by comparing it with two recent studies, with which it shows a good general agreement.The impact of the new 3-D model is illustrated through a simple synthetic experiment, which shows that the lateral variations of the wave speed can produce traveltime differences ranging from –1.5 and 1.5 s for P waves and from –2.5 and 2.5 s for S waves at local to regional distances. Such values are far larger than phase picking uncertainties and would likely affect earthquake hypocentral parameter estimations. The new 3-D model thus provides a basis for regional studies including earthquake source studies, Earth structure investigations and geodynamic modelling of Iberia and its surroundings.

Arroucau, P.↗