Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Geophysics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Isotopic gamma lines for identification of shielding materials

Identifying the constituting materials of concealed objects is crucial in a wide range of sectors, such as medical imaging, geophysics, nonproliferation, national security investigations, and so on. Existing methods face limitations, particularly when multiple materials are involved or when there are challenges posed by scattered radiation and large areal mass. Here we introduce a novel brute-force statistical approach for material identification using high spectral resolution detectors, such as HPGe. The method relies upon updated semianalytic formulae for computing uncollided flux from source of gamma radiation, shielded by a sequence of nested spherical or cylindrical materials. These semianalytical formulae make possible rapid flux estimation for material characterization via combinatorial search through all possible combinations of materials, using a high-resolution HPGe counting detector. An important prerequisite for the method is that the geometry of the objects is known (for example, from X-ray radiography). We demonstrate the viability of this material characterization technique in several use cases with both simulated and experimental data in spherical geometry.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Integrated geological, economic, and risk assessment of underground hydrogen storage

Underground hydrogen storage (UHS) is a promising option to buffer variable renewable power and support the hydrogen economy. Yet this technology is early-stage, and key uncertainties remain about how coupled processes such as geochemical and geomechanical impacts affect long-term storage security. In this paper, we seek to integrate the results from a large, multi-scale research program to quantify core feasibility metrics and situate these studies within ongoing efforts. Through this work, we evaluate hydrogen recoverability during storage by assessing diffusive losses, losses to geochemical interactions and their subsequent impact on geomechanical properties, and losses to residual saturation during injection and withdrawal. Molecular and laboratory studies show that lithologic heterogeneity, pore geometry, and mineral surface chemistry govern hydrogen–rock interactions, controlling loss pathways. High-pressure coreflooding is used to estimate relative permeabilities and measure parameters needed for reservoir-scale sensitivity analyses. These core-scale experiments found early hydrogen breakthrough at low saturations driven by capillary and viscous fingering, which limits pore-space utilization at early times, while repeated injection and withdrawal cycles improve deliverability over time. Technoeconomic factors and purity requirements for end uses are also considered as part of our comprehensive feasibility analysis. Finally, we underscore the importance of geophysical monitoring and tailored injection strategies to maintain integrity and efficiency. Together, these results provide a quantitative foundation for safe, scalable deployment of UHS and highlight priorities for future integrated studies.

08 HYDROGEN↗
DOE OSTI2024

Heuristic algorithms for design of integrated monitoring of geologic carbon storage sites

Designs for Risk Evaluation and Management (DREAM) is a tool developed under the National Risk Assessment Partnership (NRAP) to enhance geologic carbon storage safety and efficiency. Using potential leakage scenarios generated externally by the users preferred history-matching approach, DREAM constructs ideal combinations of sensor locations in the right place at the right time to detect as many leaks as possible, detect them as early as possible, and minimize cost. This user-friendly tool, developed in Java, features a window-based GUI for input and a 3D visualization tool for viewing the domain space and optimized monitoring plans. DREAM's latest version accommodates real-world usage by allowing for joint optimization of wellbore point sensor placements and surface geophysics survey geometries, and by using more efficient multi-objective optimization algorithms. We show an example where, these two improvements combined allow us to support containment assurance and go from detecting 80–90 % of the potential CO 2 leakage to +99.7 %, a step-change improvement that can make the deciding difference in whether a site is suitable for geologic carbon storage. Though developed for geologic carbon storage, this tool would be equally applicable in many surface or offshore environmental monitoring projects.

58 GEOSCIENCES↗
DOE OSTI2024

An assessment of controlled source EM for monitoring subsurface CO 2 injection at the wyoming carbonSAFE geologic carbon storage site

Here we evaluate if electromagnetic (EM) geophysical methods for monitoring geologic carbon storage (GCS) efforts at the Wyoming CarbonSAFE project adjacent to the Dry Fork Station power plant near Gillette, Wyoming. This first involved acquiring both electric and magnetic fields at eleven different locations ranging in distance from immediately adjacent to 4 km from the plant. Passive EM measurements were made to provide spectral EM noise measurements generated by electricity production at the plant and to determine if useful magnetotelluric (MT) data can be successfully collected in the region. The processed data indicate that useful MT data can be collected as long as the site is located more than 2km away from the power plant as well as active roads and rail lines. Controlled source EM data were collected using three different source configurations, two of which connected to steel casings used to complete the injection wells. Comparing the EM noise measurements to the CSEM data show measurable electric and magnetic field signals at all sites. Next a series of three-dimensional (3D) numerical models were built that simulate resistivity changes caused by the proposed CO2 injection at depths ranging from 2.4 to 3.0km. These models were used to simulate various EM measurement configurations. The modeling shows that casing-source CSEM monitoring can provide sensitivity to the injected CO 2 if source electrodes are connected to the bottom of one or both of the injection wells.

58 GEOSCIENCES↗
DOE OSTI2025

Multi-frequency progressive refinement for learned inverse scattering

Interpreting scattered acoustic and electromagnetic wave patterns is a computational task that enables remote imaging in a number of important applications, including medical imaging, geophysical exploration, sonar and radar detection, and nondestructive testing of materials. However, accurately and stably recovering an inhomogeneous medium from far-field scattered wave measurements is a computationally difficult problem, due to the nonlinear and non-local nature of the forward scattering process. We design a neural network, called Multi-Frequency Inverse Scattering Network (MFISNet), and a training method to approximate the inverse map from far-field scattered wave measurements at multiple frequencies. We consider three variants of MFISNet, with the strongest performing variant inspired by the recursive linearization method — a commonly used technique for stably inverting scattered wavefield data — that progressively refines the estimate with higher frequency content. MFISNet outperforms past methods in regimes with high-contrast, heterogeneous large objects, and inhomogeneous unknown backgrounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Joint state-parameter estimation for the reduced fracture model via the united filter

Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.

Bayesian inference↗
DOE OSTI2024

Spectral induced polarization (SIP) measurements across a PFAS-contaminated source zone

There is a pressing need for the development of field-scale, in situ screening technologies for assessing variations in aqueous film forming foam (AFFF) concentrations in soils at former fire training and storage sites. Field-scale Spectral Induced Polarization (SIP) geophysical measurements were acquired on a transect crossing an AFFF source zone. Soil samples were acquired from ten locations and used to determine variations in poly- and per-fluoroalkyl substances (PFAS) concentrations in soils and soil texture. These samples were also used to create triplicate soil columns for laboratory-grade SIP measurements. Field and laboratory observations provide evidence that SIP measurements are sensitive to the concentration of AFFF constituents associated with the pore surface in soils. The phase of the SIP measurements on the laboratory samples was linearly correlated with total soil-sorbed PFAS concentration. The phase from the field SIP measurements was highest over the location of maximum PFAS concentration measured on the laboratory samples, although a significant correlation between field-measured phase and laboratory-measured total PFAS concentration was not established. The sensitivity of the SIP response to the removal of soil PFAS using a methanol wash procedure (total PFAS concentration drop of 366 ppb) adds evidence for the case for SIP characterization of AFFF source zones. The results of these studies suggest that SIP might be developed into a field-scale technology for rapid, indirect assessment of AFFF source zones. Such a technology could improve the effectiveness of AFFF source zone characterization at reduced costs.

58 GEOSCIENCES↗
DOE OSTI2024

The Importance of Accounting for Landscape Position When Investigating Grasslands: A Multidisciplinary Characterisation of a California Coastal Grassland

Abstract Grasslands are one of the most common land‐cover types, providing important ecosystem services globally, yet few studies have examined grassland critical‐zone functioning throughout hillslopes. This study characterised a coastal grassland over a small hillslope at Point Reyes National Seashore, California, using multidisciplinary techniques, combining remotely‐sensed, geophysical, plant, and soil measurements. Clustering techniques delineated the study area into four landscape zones, up‐, mid‐, and down‐slope, and a bordering riparian ecotone, which had distinct environmental properties that varied spatially across the site, with depth, and time. Soil moisture increased with depth and down slope towards a bordering riparian zone, and co‐varied with soil CO 2 flux rates both spatially and temporally. This highlighted three distinct controls of soil moisture on soil respiration: CO 2 fluxes were inhibited by high moisture content in the down‐slope during the wet winter months, and converged across landscape positions in the dry summer months, while also displaying post‐rain pulses. The normalised difference vegetation index (NDVI) ranged from 0.32 (September)–0.80 (April) and correlated positively with soil moisture and aboveground biomass, moving down slope. Yet, NDVI, aboveground biomass, and soil moisture were not correlated to soil organic carbon (SOC) content (0.4%–4.5%), which was highest in the mid‐slope. The SOC content may instead be linked to shifts in dominant grassland species and their rhizosphere properties with landscape position. This multidisciplinary characterisation highlighted significant heterogeneity in grassland properties with landscape position, and demonstrated an approach that could be used to characterise other critical‐zone environments on hillslopes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sound Velocities of Stishovite at Simultaneous High Pressure and High Temperature Suggest an Eclogite‐Rich Layer Beneath the Hawaii Hotspot

Compressional and shear wave velocities of polycrystalline stishovite (SiO 2 ) have been measured at simultaneous high pressures and temperatures up to 14.5 GPa and 800°C. By fitting velocities to the finite strain equations, the elastic moduli and density were determined to be K S0 = 306.6(46) GPa, $K$$^{′}_{S}$ = 4.92(10), ∂K S /∂T = −0.024(1) GPa/K, G 0 = 229.0(34) GPa, G′ = 1.07(10), ∂G/∂T = −0.017(1) GPa/K, ρ 0 = 4.287(2) g/cm 3 . Our modeling suggested that, in the eclogite, coesite-stishovite transition can increase P and S wave velocities by 2.4% and 3.5%, respectively. A comparison between geophysical observations and our model shows that the coesite-stishovite phase transition in the eclogite can potentially be responsible for the occurrence of the X discontinuity beneath Hawaii. In addition, our current results suggest an eclogite-rich layer between 340 and 450 km depth beneath Hawaii. The eclogite concentration at the top and bottom of the layer is 41–55 vol% and >77 vol%, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Reimagining Earth in the Earth System

Abstract Terrestrial, aquatic, and marine ecosystems regulate climate at local to global scales through exchanges of energy and matter with the atmosphere and assist with climate change mitigation through nature‐based climate solutions. Climate science is no longer a study of the physics of the atmosphere and oceans, but also the ecology of the biosphere. This is the promise of Earth system science: to transcend academic disciplines to enable study of the interacting physics, chemistry, and biology of the planet. However, long‐standing tension in protecting, restoring, and managing forest ecosystems to purposely improve climate evidences the difficulties of interdisciplinary science. For four centuries, forest management for climate betterment was argued, legislated, and ultimately dismissed, when nineteenth century atmospheric scientists narrowly defined climate science to the exclusion of ecology. Today's Earth system science, with its roots in global models of climate, unfolds in similar ways to the past. With Earth system models, geoscientists are again defining the ecology of the Earth system. Here we reframe Earth system science so that the biosphere and its ecology are equally integrated with the fluid Earth to enable Earth system prediction for planetary stewardship. Central to this is the need to overcome an intellectual heritage to the models that elevates geoscience and marginalizes ecology and local land knowledge. The call for kilometer‐scale atmospheric and ocean models, without concomitant scientific and computational investment in the land and biosphere, perpetuates the geophysical view of Earth and will not fully provide the comprehensive actionable information needed for a changing climate.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

4D Electrical Resistivity Imaging of Stress Perturbations Induced During High–Pressure Shear Stimulation Tests

Fluid flow through fractured media is typically governed by the distribution of fracture apertures, which are in turn governed by stress. Consequently, understanding subsurface stress is critical for understanding and predicting subsurface fluid flow. Although laboratory–scale studies have established a sensitive relationship between effective stress and bulk electrical conductivity in crystalline rock, that relationship has not been extensively leveraged to monitor stress evolution at the field scale using electrical or electromagnetic geophysical monitoring approaches. In this paper we demonstrate the use time–lapse 3–dimensional (4D) electrical resistivity tomography to image perturbations in the stress field generated by pressurized borehole packers deployed during shear–stimulation attempts in a 1.25 km deep metamorphic crystalline rock formation.

58 GEOSCIENCES↗
DOE OSTI2025

A Unified Analytical Model for Pressure Solution With Fully Coupled Diffusion and Reaction

Abstract Geophysical models for pressure solution are typically developed for diffusion‐controlled or reaction‐controlled scenarios. We present a unified analytical model that considers fully coupled diffusion and reaction during pressure solution. The model recovers the diffusion‐controlled and reaction‐controlled models in the literature as specific limiting cases. When diffusion and reaction exhibit comparable influences, we validate the proposed model against independent numerical simulations. The proposed model is then employed in interpreting experimental measurements, demonstrating a better agreement compared to previous interpretations.

Wang, Ziyan [Energy Geosciences Division Lawrence ↗
DOE OSTI2024

Measurement of the Static Nonlinear Third-Order Elastic Moduli of Rocks: Problems and Applicability

The third-order elastic (TOE) model has been used to describe the widely observed nonlinear mechanical behaviors of earth materials. In addition to linear elastic constants ( λ , μ ), three nonlinear elastic moduli ( A , B , C ) are required for isotropic rocks. Contrary to previous research on dynamic TOE moduli, this study followed the protocol to measure strain and stress under uniaxial and hydrostatic compressive tests statically, which were later used to invert for the full set of TOE moduli for four standard rock types with differing pore structures of a porous oolitic limestone, quartz-rich sandstones, and a dense crystalline basalt. The applicability of the TOE model to characterize nonlinearity depends on the fulfillment of path-independence and small-strain assumptions. Using the measured static TOE moduli, the finite element model demonstrates that the stress in the vicinity of the wellbore is more amplified than the stress in the linear elastic case, which leads to a wider zone of rock failure around the wellbore. Due to the long-standing discrepancy between static and dynamic moduli, the rarely reported full set of static TOE moduli is necessary and will benefit future research in understanding the effect of rock nonlinearity on geophysical and geomechanical applications, such as long-term safe storage of CO 2 and generating process of geohazards.

58 GEOSCIENCES↗
DOE OSTI2025

Revisiting the Relationship Between Induced Polarization and Surface Conductivity: Ratios From Laboratory to Field

Abstract Among the subsurface geophysical methods used in the critical zone investigations, induced polarization (IP) shows great vitality thanks to its unique ability to assess porosity via bulk conduction and estimate permeability through surface conduction and/or polarization. However, such an advantageous separation between bulk and surface is mostly implemented by multi‐salinity experiments in the laboratory, which is incredibly difficult to realize in the field. One promising approach to address such an obstinate issue is to gauge the surface conductivity ( σ s ) from the quadrature conductivity ( σ ″) or normalized chargeability ( M n ) with the ratios between ( l = σ ″/ σ s , l mn = M n / σ s ). While these ratios are known not to be universal, the underlying principles are not fully understood and relevant theoretical studies are rare, which makes quantitative IP applications difficult. Here we scrutinize the conduction and polarization mechanisms of geomaterials and pinpoint that the two ratios are inherently functions of salinities and frequencies rather than only determined by the properties of the electrical double layer (EDL), hence representative samples from the investigated field must be calibrated in the laboratory and a characteristic frequency should be chosen for their usage. Besides the macro‐scale ratios l and l mn , we define two micro‐scale ratios χ and χ mn directly from the EDL, such that the new ratios exclude the effect of salinity and frequency and offer the opportunity to characterize and monitor changes of the EDL. Our study demonstrates that the existing macro‐scale ratios converge toward the values of novel micro‐scale ratios at high water salinity.

Qi, Youzheng [Earth and Environmental Sciences Are↗
DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗
DOE OSTI2025

Unsupervised Clustering of Microseismic Events and Focal Mechanism Analysis at the CO 2 Injection Site in Decatur, Illinois

Characterization of induced microseismicity at a carbon dioxide (CO 2 ) storage site is critical for preserving reservoir integrity and mitigating seismic hazards. We apply a multilevel machine learning (ML) approach that combines the nonnegative matrix factorization and hidden Markov model to extract spectral representations of microseismic events and cluster them to identify seismic patterns at the Illinois Basin-Decatur Project. Unlike traditional waveform correlation methods, this approach leverages spectral characteristics of first arrivals to improve event classification and detect previously undetected planes of weakness. By integrating ML-based clustering with focal mechanism analysis, we resolve small-scale fault structures that are below the detection limits of conventional seismic imaging. Our findings reveal temporal bursts of microseismicity associated with brittle failure, providing insights into the spatio-temporal evolution of fault reactivation during CO 2 injection. This approach enhances seismic monitoring capabilities at CO 2 injection sites by improving fault characterization beyond the resolution of standard geophysical surveys.

Willis, Rachel Marie [Sandia National Laboratories↗
DOE OSTI2025

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗
DOE OSTI2026

Anisotropic Melt Inclusions as a Confounding Signal for Ice‐Penetrating Radar Observations

Ice-penetrating radar is a powerful geophysical tool for understanding the subsurfaces of Earth, Mars, and icy moons. Radar reflectivity, attenuation, and birefringence are used to infer subsurface hydrology, englacial temperature, water content, and crystal orientation fabric. However, conventional radar sounding analyses either ignore melt or use classical mixing models which assume spherical melt inclusions, obscuring anisotropic contributions of melt. Here, we use geometric mixing models to calculate the reflectivity, attenuation, and birefringence of temperate ice containing anisotropic melt. We find that anisotropic melt can introduce significant deviations in radar measurements. For instance, melt anisotropy may impact reflectivity-based estimates of water content by up to 30% volume fraction, while attenuation-based estimates may vary by up to 43%. Critically, a melt volume fraction of just $f$ ~ 10 -5 can reproduce birefringence signals previously attributed to ice fabric. We conclude that melt anisotropy may significantly impact analyses of englacial hydrology, subglacial hydrology, and ice fabric.

58 GEOSCIENCES↗