Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “geoscience”

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,333 records · Page 74

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗

Towards the next generation of Geospatial Artificial Intelligence

Geospatial Artificial Intelligence (GeoAI), as the integration of geospatial studies and AI, has become one of the fastest-developing research directions in spatial data science and geography. This rapid change in the field calls for a deeper understanding of the recent developments and envision where the field is going in the near future. In this work, we provide a quantitative analysis of the GeoAI literature from the spatial, temporal, and semantic aspects. We briefly discuss the history of AI and GeoAI by highlighting some pioneering work. Then we discuss the current landscape of GeoAI by selecting five representative subdomains including remote sensing, urban computing, Earth system science, cartography, and geospatial semantics. Finally, we highlight several unique future research directions of GeoAI which are classified into two groups: GeoAI method development challenges and GeoAI Ethics challenges. Topics include heterogeneity-aware GeoAI, knowledge-guided GeoAI, spatial representation learning, geo-foundation models, fairness-aware GeoAI, privacy-aware GeoAI, as well as interpretable and explainable GeoAI. We hope our review of GeoAI’s past, present, and future is comprehensive and can enlighten the next generation of GeoAI research.

58 GEOSCIENCES↗

Detecting an underground tunnel by applying joint traveltime and waveform inversion

Underground tunnel and void detection is a challenging geophysical problem, and many methods have been proposed. Seismic techniques are promising because of the large seismic velocity contrasts between the air-filled void and the surrounding sediment and concrete. We apply a joint seismic traveltime and waveform inversion method to image a buried tunnel with concrete walls and a void space inside. The joint inversion images the top of the concrete tunnel as a high-velocity anomaly and the void space as a low-velocity anomaly. The location of the velocity anomalies predicted by the method agrees with the known location of the tunnel. As a comparison, the stand-alone full waveform inversion is also applied to the data. The first-arrival traveltime tomography shows weak nonlinearity but fails to image the hidden low-velocity layer. Full waveform inversion is able to image complex near-surface structures, but fitting waveforms may not honor the traveltime fit, especially when the data contain noise. Synthetic and real data tests show that the joint inversion method retains the advantages of both traveltime inversion and full waveform inversion and overcomes their respective drawbacks at the same time. The field example shows the joint inversion provides a better reconstruction of the high-velocity feature representing the top of the concrete wall in terms of its magnitude and location.

58 GEOSCIENCES↗

The effect of magnetic susceptibility and magnetic field strength on porosity estimates determined from low-field nuclear magnetic resonance

The presence of materials with high magnetic susceptibility are known to have an impact on 1H nuclear magnetic resonance (NMR) measurements and in laboratory data this often results in poor estimates of porosity from NMR data. To quantify and understand the cause of the poor NMR porosity estimations, in this laboratory study, we examine the effect of magnetic susceptibility, NMR measurement parameters, and NMR instrument design on NMR measurements. Data were collected with two instruments, each with a different Larmor frequency (2 MHz and 485 kHz), on water-saturated unconsolidated sediments with magnetic susceptibility values ranging from 3.6 × 10 -6 to 7020 × 10 -6 SI. The results show that for materials with low magnetic susceptibility (<378 × 10 -6 SI), the gravimetric porosity is accurately predicted from the NMR measurements. For the samples with high magnetic susceptibility (>987 × 10 -6 SI) the gravimetric porosity is poorly predicted from the 2 MHz NMR measurements made at all echo times (from 0.2 to 3.0 ms). In contrast, the gravimetric porosity is more accurately predicted at an echo time of 0.2 ms for measurements made using the 485 kHz instrument, although at larger echo times (>1.0 ms), the porosity estimate becomes poor. The 485 kHz NMR instrument has non-zero internal magnetic field gradients, similar to those found in borehole instruments, in contrast to the 2 MHz NMR instrument, which has a homogeneous applied magnetic field. Finally, we conclude that differences is the magnetic field strength and higher magnetic field inhomogeneities in the 485 kHz NMR instrument contribute to a reduction of the impact of inhomogeneities in the magnetic field caused by materials with high magnetic susceptibility, allowing for improved porosity estimation. These results indicate that NMR measurements collected at short echo times in low, inhomogeneous static fields, e.g., borehole instruments, may provide accurate estimates of porosity in water saturated sediments, even in the presence of magnetic minerals.

58 GEOSCIENCES↗

Bedrock architecture, soil texture, and hyporheic zone characterization combining electrical resistivity and induced polarization imaging

The shallow (<15 m depth) subsurface environment of a short reach along a bedrock stream was investigated with electrical resistivity and induced polarization (IP) to map details of bedrock and soil sediments. The bedrock is mostly comprised of limestone, and is generally resistive compared to the overlying soil. The soil-bedrock interface was determined through a trial and error approach using a sharp boundary feature in the inversion model. The inferred bedrock surface determined from the inversions exhibited undulating patterns with troughs and ridges. A near continuous trough ran alongside the stream within the floodplain, and is suggestive of a paleochannel among other interpretations of this feature. The structure of the electrical resistivity above the bedrock showed small-scale elongated features. The chargeability from the IP method showed larger scale features. High values of chargeability were associated with the sediments in the floodplain, and low values were associated with bedrock, stream, and soil on the elevated banks above the stream. If the chargeability is associated with membrane polarization characteristic of clayey soils, then IP seems to highlight the mere existence of clay, while the resistivity may be more discernable of the relative proportion of clay. Vegetation differences may also explain the chargeability distribution, where parts of the survey with high chargeability had dense pine with no understory making the soils more organically rich.

54 ENVIRONMENTAL SCIENCES↗

Stratigraphic Identification with Airborne Electromagnetic Methods at the Hanford Site, Washington

Stratigraphic units can influence the fate and transport of subsurface contaminants within groundwater. Units having coarse-grained sediments act as preferential flow pathways, and therefore can accelerate the transport of contaminants to reach human and ecological receptors. At legacy waste sites, detailed knowledge of subsurface stratigraphy can be used for effective monitoring and remediation planning to help minimize risk to human health and the environment. Airborne electromagnetic (AEM) methods can non-invasively provide information on kilometer-scale or larger subsurface stratigraphic features and fill informational gaps in directly sampled data from sparsely located boreholes. In this paper, we present inversion results of a 412 line-km frequency-domain AEM survey to delineate subsurface stratigraphic features at the Hanford Site, located in southeastern Washington State. The inversion was performed using a massively parallel 3D electromagnetic modeling and inversion code, where the modeling is based on solving frequency-domain Maxwell’s equations using an unstructured-mesh finite-element method and the inversion employs a Gauss-Newton optimization scheme. The results are compared to an underlying geologic framework model (GFM), built by interpolating contact depths of stratigraphic units interpreted from site borehole datasets. In areas with good borehole coverage, the inversion results show a good match with the GFM to a depth of about 60 m. Outside of these areas, the inversion results exhibit inconsistencies from the assumptions made to create the GFM, demonstrating that the AEM survey results can be used to improve the understanding of the geological conceptual model.

47 OTHER INSTRUMENTATION↗

Interferometric SAR modelling of near surface data to improve geological model in the Surat Basin, Australia

Here we present a study on geophysical inverse modelling for subsurface structural properties of an unconventional hydrocarbon site that was monitored previously by Interferometric Synthetic Aperture Radar (InSAR) technology for surface deformation. A static three-dimensional geomodel along with extracted property maps replicates the depth of each underlying stratigraphic unit and structural feature with the density of each geological layer. We examine the hypothesis that integration of elastic properties of each formation layer with InSAR observations in a stratified elastic medium will lead to a viscoelastic geophysical inverse problem that can be solved to estimate fractional volume change at the reservoir level. Moreover, we examine synthetic scenarios in which the elastic properties of the formations are perturbed before determining the resulting impact on the rate of surface deformation. The results show that although the slope of underlying formations, their density and depth can define the extent and pattern of a deformation signal, their properties have a marginal impact on volumetric change compared to the dense network of shallow depth Coal Seam Gas (CSG) mining wells. Besides, it is also demonstrated that the inversion of InSAR deformation maps can resolve the uncertainties associated with low-resolution seismic interpretation as well as filling the data gaps within seismic acquisitions. A significant contribution of this investigation to the geological basin modelling involves a) introducing a remote and non-invasive technology such as InSAR to improve geophysical mapping of subsurface structures such as faults in areas with sparse or no reflective seismic information, and b) applying a multi-layer viscoelastic geophysical source model for an unconventional hydrocarbon reservoir such as CSG.

58 GEOSCIENCES↗

Deep learning inversion of gravity data for detection of CO 2 plumes in overlying aquifers

In this work, we developed an effective U-Net based deep learning (DL) model for inversion of surface gravity data on a rectangular grid to predict 2-D high-resolution subsurface CO 2 distribution along a vertical cross-section due to CO 2 leakage through a wellbore within a deep CO 2 storage reservoir. We used synthetic data to model two types of CO 2 leakage scenarios: one CO 2 plume in a shallow aquifer (single plume case), and two plumes present at different depths (double plume case). The 3-D synthetic plume samples were created by sampling among predetermined CO 2 plume depths, saturations, and volumes. The corresponding surface gravity data on a rectangular grid were generated by a 3-D forward model. The U-Net model detected 72% of single-plume samples, and one or both plumes in 75% of double-plume samples. Most of the undetected single plumes have small gravity field strengths below the typical noise level of 5 μGal. This model generated reproducible, reliable predictions with acceptable errors and demonstrated improved spatial resolution over the conventional least-squares inversion. In contrast to the conventional least-squares inversion, which often overestimates the size of its target and underestimates its density, this U-Net model accurately delineated the boundary of a target. Furthermore, this DL inversion detected deep, small, or low saturation CO 2 plumes that are often more difficult to resolve with conventional gravity inversion methods. We note the limitations of this feasibility study, including the use of synthetic data with regular CO 2 plume shapes, and the prediction of a 2-D plume cross-section rather than the full 3-D plume, as well, we recognize the lower detection fraction for double-plume scenarios. Nevertheless, this study demonstrates that DL gravity inversion is a promising and potentially superior method to conventional least-squares inversion. Our U-Net based deep learning inversion approach may be adapted for inversion of other types of geophysical data. DL inversion can facilitate near real-time monitoring of geologic carbon sequestration to provide site operators with prompt information about subsurface CO 2 distribution for risk management and mitigation.

58 GEOSCIENCES↗

SREMI: Super-resolution electromagnetic imaging with single-channel ground-penetrating radar

High-resolution near-surface imaging has important applications in civil engineering, infrastructure inspection, military threat detection, geological characterization, and lunar and planetary exploration. Zero-offset, singlechannel ground penetrating radar (GPR) imaging is an established technique for near-surface target imaging and sensing but often suffers from low spatial resolution and imaging artifacts, especially of deep structures. In response, we formulate the GPR imaging as a dual-sparsity optimization problem, and develop a super-resolution electromagnetic imaging method based on a fast iterative shrinkage-thresholding algorithm. We develop our GPR imaging method in the framework of electromagnetic exploding-reflectors simulation theory, therefore the imaging method is computationally efficient. In this work, we demonstrate through synthetic and field data examples that our method can produce sharper, more reliable images with fewer artifacts compared with single-pass reverse-time migration GPR method, thus leading to improved near-surface interpretation and object identification.

58 GEOSCIENCES↗

Comparison of machine learning and electrical resistivity arrays to inverse modeling for locating and characterizing subsurface targets

Here, this study evaluates the performance of multiple machine learning (ML) algorithms and electrical resistivity (ER) arrays for inversion with comparison to a conventional Gauss-Newton numerical inversion method. Four different ML models and four arrays were used for the estimation of only six variables for locating and characterizing hypothetical subsurface targets. The combination of dipole-dipole with Multilayer Perceptron Neural Network (MLP-NN) had the highest accuracy. Evaluation showed that both MLP-NN and Gauss-Newton methods performed well for estimating the matrix resistivity while target resistivity accuracy was lower, and MLP-NN produced sharper contrast at target boundaries for the field and hypothetical data. Both methods exhibited comparable target characterization performance, whereas MLP-NN had increased accuracy compared to Gauss-Newton in prediction of target width and height, which was attributed to numerical smoothing present in the Gauss-Newton approach. MLP-NN was also applied to a field dataset acquired at U.S. DOE Hanford site.

54 ENVIRONMENTAL SCIENCES↗

Key Elements of Auroral Substorm Development and their Relationship to Recent Observations of Detached Sub-Auroral Phenomena Including STEVE-like Emissions

Since Akasofu introduced the concept of the substorm in 1964, numerous refinements have been added to the phenomenological model of how auroral substorms develop. Here, we review several of the most important of these including; polar cap patches, intensifications of the Polar Activated Band (PAB) at the poleward edge of the bulge, the development of new arcs poleward, the generation of auroral streamers, the evolution of streamers into auroral torches and omega bands, the association of streamers with particle injection and BBF activity in the tail, and wedgelet-like magnetic perturbations on the ground. We also review the concept of “contact breakups” which are breakups that appear to be triggered by the arrival of streamer-like forms in the equatorward regions of the oval. Finally, we remind the community that the addition of detached sub-auroral emissions was made to the phenomenological picture of substorm develop in the late 1960s along with the concept of auroral streamers. Newer observations of detached sub-auroral (STEVE-like emissions) from the Viking/UVI and POLAR/VIS imagers are presented which confirm that they are east-west aligned bands of emission that separate away from the equatorward edge of the auroral oval in response to intensifications of the poleward edge of the bulge and subsequent streamer production. We propose a new model for the formation of detached sub-auroral (STEVE-like) emissions in which a non-linear growth of a SAPS-driven instability (e.g. like the shear-flow/ballooning instability) at the plasmapause results in the disruption of the boundary separating hot plasma sheet particle from the cold plasmaspheric plasma. It is proposed that the resulting intermixing of plasma populations leads to both the observed STEVE-like emissions and also provides a source of cold plasma on open drift paths that can feed the long-lived drainage plumes that have recently been discovered.

58 GEOSCIENCES↗

Influence of support on Rh-Co bimetallic catalysts for ethylene hydroformylation

Bimetallic Rh- and Co-based catalysts are promising materials for the heterogenous ethylene hydroformylation reaction. Here, in this study, the influence of a mesoporous silica (SBA-15) support on Rh and Co bimetallic interactions was investigated through a comparison with silica gel and alumina supports. The bimetallic catalyst supported on mesoporous silica (RhCo 3 /SBA-15) showed the best C 3 oxygenate yield among the three bimetallic catalysts. In-situ vibrational studies suggested moderate binding of the gem-dicarbonyl and Rh(CO)(C 2 H 4 ) intermediates due to this bimetallic interaction that led to improved hydroformylation performance. Kinetic studies revealed a lower hydroformylation barrier for the bimetallic compared to a Rh monometallic catalyst, and in-situ X-ray absorption spectroscopy investigations showed clear Rh-Co alloy formation on RhCo 3 /SBA-15. The bimetallic enhancement effect from the interaction with the mesoporous silica support shown here can be further optimized to design olefin hydroformylation catalysts.

58 GEOSCIENCES↗

Molecular-level understanding of gibbsite particle aggregation in water

Using molecular dynamics simulations, we investigate the molecular scale origin of crystal face selectivity when one gibbsite particle attaches to another in water. A comparison of the free energy per unit surface area of particle–particle attachment indicates that particle attachment through edge surfaces, where the edge surfaces are either (1 0 0) or (1 1 0) crystal faces, is more energetically favorable compared to attachment between two basal surfaces (i.e., (0 0 1) crystal faces) or between the basal surface of one particle and the edge surface of another. This result suggests that gibbsite crystals with low basal/edge surface area ratio will preferentially attach through edge surfaces, potentially helping the crystals grow laterally. However, for larger gibbsite particles (high basal/edge surface area ratio) the total free energy, not normalized by surface area, of particle attachment through the basal surfaces is lower (more negative) than attachment through the edge surfaces, indicating that larger gibbsite particles will preferentially aggregate through basal surface attachments. The short-range electrostatic interactions including the interparticle hydrogen bonds from surface –OH groups drive particle attachment, and the dominant contribution to the free energy minimum is enthalpic rather than entropic. However, the enthalpy of basal-edge attachment is significantly offset by the entropy leading to a higher free energy (less negative) compared to that of basal-basal attachment. Study of the free energy for a few imperfect attachments of two particles indicates a higher free energy (i.e., less negative, less stable), compared to a perfect attachment

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

La 3+ and Y 3+ interactions with the carboxylic acid moiety at the liquid/vapor interface: Identification of binding complexes, charge reversal, and detection limits

Specific interactions of yttrium and lanthanum ions with a fatty acid Langmuir monolayer were investigated using vibrational sum frequency spectroscopy. The trivalent ions were shown to interact with the charged form of the carboxylic acid group from nanomolar concentrations (<300 nM). Analysis of the spectral features from both the symmetric and the asymmetric carboxylate modes reveals the presence of at least three distinct coordination structures linked to specific binding configurations. Although the same species were identified for both La 3+ and Y 3+ , they display a different concentration dependence, highlighting the ion-specificity of the interaction. From the analysis of the response of interfacial water molecules, the reversal of the surface charge, as well as the formation of yttrium hydroxide complexes, were detected upon increasing the amount of salt in solution. The binding interaction and kinetics of absorption are sensitive to the solution pH, showing a distinct ion speciation in the interfacial region when compared to the bulk. Changing the subphase pH or adding a monovalent background electrolyte that promotes deprotonation of the carboxylic acid headgroup could further improve the detection limit of La 3+ and Y 3+ to concentrations < 100 nM. These findings demonstrate that nM concentrations of trace metals contaminants, typically found on monovalent salts, can significantly influence the binding structure and kinetics in Langmuir monolayers.

01 COAL, LIGNITE, AND PEAT↗

Impact of organic solutes on capillary phenomena in water-CO2-quartz systems

Hypothesis: The migration of supercritical CO2 (scCO2) injected into underground reservoirs as part of carbon capture and storage is influenced by organic contamination affecting mineral wettability. Molecular dynamics (MD) simulations of relevant systems that incorporate representative organic solutes allow detailed investigation of changes in fundamental interfacial and capillary properties. Experiments: We use MD simulations to explore the effects of four organic solutes (quinoline, decanoic acid, coronene, sorgoleone) on the wettability of quartz by water in the presence of scCO2. We examine the impacts of polar, alkyl, and aromatic moieties as well as fluid flow velocity at elevated temperatures and pressures. Findings: Organic molecules accumulate at the water-CO2 interface, where they distribute according to their size and functional groups. Certain organics penetrate the adsorbed water film at the quartz-CO2 interface, revealing two modes of hydrogen bonding between polar organic functional group, water, and quartz surface –OH groups. Interfacial energies and contact angles are minimally impacted by organic adsorption at the water-CO2 interface, possibly due to simultaneous CO2 desorption. Non- equilibrium MD simulations reveal that flow-induced redistribution of organic compounds modulates the radii of curvature of the advancing and receding water-CO2 interfaces. Our results indicate that organic adsorption on water surfaces is likely ubiquitous during multi-phase flow in soils and sedimen- tary rocks, with implications for the mobilization and transport of organic compounds.

36 MATERIALS SCIENCE↗