Search NASA⌕ Search

SEARCH · Search NASA

Results for “data and proxy modeling”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Entropy-Assisted Quality Pattern Identification in Finance

Short-term patterns in financial time series form the cornerstone of many algorithmic trading strategies, yet extracting these patterns reliably from noisy market data remains a formidable challenge. In this paper, we propose an entropy-assisted framework for identifying high-quality, non-overlapping patterns that exhibit consistent behavior over time. We ground our approach in the premise that historical patterns, when accurately clustered and pruned, can yield substantial predictive power for short-term price movements. To achieve this, we incorporate an entropy-based measure as a proxy for information gain: patterns that lead to high one-sided movements in historical data yet retain low local entropy are more “informative” in signaling future market direction. Compared to conventional clustering techniques such as K-means and Gaussian Mixture Models (GMMs), which often yield biased or unbalanced groupings, our approach emphasizes balance over a forced visual boundary, ensuring that quality patterns are not lost due to over-segmentation. By emphasizing both predictive purity (low local entropy) and historical profitability, our method achieves a balanced representation of Buy and Sell patterns, making it better suited for short-term algorithmic trading strategies. This paper offers an in-depth illustration of our entropy-assisted framework through two case studies on Gold vs. USD and GBPUSD. While these examples demonstrate the method’s potential for extracting high-quality patterns, they do not constitute an exhaustive survey of all possible asset classes.

Physics↗
DOE OSTI2024

CCUS 2024, Interpreting the strain tensor Larry Murdoch Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring during CCUS. We have demonstrated this method by deploying strainmeters at shallow depths (30 to 40m) and then conducting injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. We have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS. Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and our objective was to evaluate opportunities for strain monitoring during characterization and monitoring for CCUS. Our approach was to deploy strainmeters at shallow depths (30 to 40m) and then conduct injection well tests in an underlying reservoir at 530m depth. The results indicate that the horizontal strain at shallow strainmeters was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 ne/d over a few days (1 nanostrain = 1 part per billion strain). We then used the strain data to estimate reservoir properties, geometry and pressure through inversion of poroelastic forward models using both numerical and novel analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed fast, closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. This approach markedly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. Bayesian inversion is important because it provides predictions with uncertainties, which makes the results useful for decision making. In conclusion, we have shown with field tests and simulations that the strain tensor in the caprock is sensitive to pressure in the reservoir, reservoir properties and boundaries, and pressure in the caprock caused by leaks. These results indicate that measuring and interpreting the shallow strain tensor could be a valuable tool for both initial reservoir characterization efforts and long-term monitoring during CCUS.

Murdoch, Larry↗
DOE OSTI2021

Nonmonotonic Energy Dependence of Net-Proton Number Fluctuations

Nonmonotonic variation with collision energy (√ s NN) of the moments of the net-baryon number distribution in heavy-ion collisions, related to the correlation length and the susceptibilities of the system, is suggested as a signature for the quantum chromodynamics critical point. Here, we report the first evidence of a nonmonotonic variation in the kurtosis times variance of the net-proton number (proxy for net-baryon number) distribution as a function of √ s NN with 3.1 σ significance for head-on (central) gold-on-gold (Au+Au) collisions measured solenoidal tracker at Relativistic Heavy Ion Collider. Data in noncentral Au+Au collisions and models of heavy-ion collisions without a critical point show a monotonic variation as a function of √ s NN.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Improving E3SM Land Model Photosynthesis Parameterization via Satellite SIF, Machine Learning, and Surrogate Modeling

The parameterization of key photosynthesis parameters is one of the key uncertain sources in modeling ecosystem gross primary productivity (GPP). Solar-induced chlorophyll fluorescence (SIF) offers a good proxy for GPP since it marks the actual process of photosynthesis; while machine learning (ML) provides a robust approach to model the GPP-SIF relationship. Here, we trained the boosted regressing tree (BRT) and the Random Forest ML models with Greenhouse Gases Observing Satellite SIF data and in situ GPP observations from 49 eddy covariance towers. These trained ML GPP-SIF models were fed into the Energy Exascale Earth System Model (E3SM) Land Model (ELM) to generate ELM-simulated global SIF estimates, which were then benchmarked against satellite SIF observations with a surrogate modeling approach. Our results indicated good modeling performance of the ML-based GPP-SIF relationship. The ELM model when fed with the ML GPP-SIF models also can well predict the spatial-temporal variations in SIF. We also found high model accuracy for the surrogate modeling. Model parameter sensitivity analysis suggested that the fraction of leaf nitrogen in RuBisCO (flnr) is the most sensitive parameter to the SIF; other sensitive parameters include the Ball-Berry stomatal conductance slope (mbbopt) and the vcmax entropy (vcmaxse). The posterior uncertainty in simulated GPP was greatly reduced after benchmarking, and the model produced improved spatial patterns of mean GPP relative to FLUXCOM GPP. Our integrated approach provides a new avenue for improving land models and using remote-sensing SIF, which can be further improved in the future with more ground- and satellite-based observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Potential of hotspot solar‐induced chlorophyll fluorescence for better tracking terrestrial photosynthesis

Abstract Remote sensing of solar‐induced fluorescence (SIF) opens a new window for quantifying a key ecological variable, the terrestrial ecosystem gross primary production (GPP), because of the revealed strong SIF–GPP correlation. However, similar to many other remotely sensed metrics, SIF observations suffer from the sun‐sensor geometry effects, which may have important impacts on the SIF–GPP relationship but remain poorly understood. Here we used remotely sensed SIF, globally distributed tower GPP data, and a mechanistic model to provide a systematic analysis. Our results reveal that leaf physiology, canopy structure, and sun‐sensor geometries all affect the SIF–GPP relationship. In particular, we found that SIF observations in the sun‐tracking hotspot direction can be a better proxy of GPP due to the similar responses of light use efficiency and SIF escaping probability in the hotspot direction to the increasing incoming solar radiation. Such conclusions are supported by a variety of modeling simulations and satellite observations over various plant function types, at different time scales and with satellite observational modes. This study demonstrates the potential and advantage of normalizing SIF observations to the hotspot direction for better global GPP estimations. This study also demonstrates the great potentials of current and future spaceborne sun‐tracking satellite missions for a significant improvement in measuring and monitoring, at a wide range of spatial and temporal scales, the changes in terrestrial ecosystem GPP in response to anticipated changes in the Earth's environmental conditions.

Hao, Dalei↗
DOE OSTI2024

An Insight-Centric Paradigm for Data Reduction and Inference Speed Improvement at the Scurry Area Canyon Reef Operator’s Committee (SACROC) Unit

The poster presents the work conducted under SMART focusing on using insight-centric approach to design a meaningful proxy for machine learning. Domain insights are critical not just in understanding the prediction results but also in designing the model. This study demonstrated that a single meaningful scaler (as an extreme case) can effectively replace full-size 3D geologic properties. The model's accuracies are on par with other models, and it is the fastest model to predict all test cases, 5,000 times faster than traditional simulations.

Shih, Chung Yan↗
DOE OSTI2020

A control-oriented combustion model framework for compression ignition engines operating on low-reactivity fuel

This work focuses on zero-dimensional modeling of the heat release rate in a compression ignition engine operating on gasoline-like fuels. Due to the properties of gasoline, such as high volatility and longer ignition delay than diesel, the injection strategies can vary significantly from the operation with conventional diesel fuel. Different injection strategies are commonly used to achieve varying degrees of in-cylinder stratification in order to shape the combustion event and maximize efficiency. The proposed zero-dimensional combustion model was developed to account for the different stages in combustion caused by the fuel stratification. As the ignition delay model is an integral part of the entire combustion process and significantly affects the prediction accuracy, special attention has been paid to local phenomena influencing ignition delay. A one-dimensional spray model by Musculus and Kattke was employed in conjunction with a Lagrangian tracking approach in order to estimate the local air–fuel ratio within the spray tip, as a proxy for reactivity. The local air–fuel ratio, in-cylinder temperature and pressure were used in an integral fashion to estimate the ignition delay. Heat release rates were modeled using first-order non-linear differential equations. The proposed combustion model was validated against experimental data of a heavy-duty compression ignition engine with up to three injection events at mostly 1038 r/min and 14 bar brake mean effective pressure. Further validation of the model was carried out at other engine loads and speeds. Model prediction errors in CA50 of less than 1 °CA across all conditions were found. Modeling results of other combustion metrics such as combustion duration and indicated mean effective pressure are also highly satisfactory. In addition, the model has been shown to be capable of estimating the ringing intensity for most conditions.

Pamminger, Michael↗
DOE OSTI2021

Hamilton: Flexible, Open Source $10 Wireless Sensor System for Energy Efficient Building Operation

Sensors for improving building performance are rapidly populating the market, driven in part by the drive to reduce greenhouse gas emissions resulting from energy production as well as improve the interior environment for healthy and more productive spaces. UC Berkeley has led wireless sensor development over the past 25 years (e.g., Telos mote), with the Hamilton (named after Alexander Hamilton on the US $10 bill) as the most recent. The Hamilton sensor was designed as a low-cost high-performance sensor that is modular and interoperable. The objective of the Hamilton project was to create, evaluate and establish the technological foundations for secure and easy to deploy building energy efficiency applications utilizing pervasive, low-cost wireless sensors integrated with traditional Building Management Systems (BMS), consumer-sector building components, and powerful data analytics. The project included iterative hardware design, incorporating a high-performance database (BTrDb, http://btrdb.io/), creating and iterating the development of secure data middleware (BOSSwave, WAVE/WAVEMQ), working with and pushing the development of an open-source tiny operating system RiotOS, and implementing and improving protocols such as Thread/OpenThread and TCP/IP. The hardware benefited from careful design to drive down the cost; the design included a System-on-a-Chip (SoC), chip antenna, single crystal and five passive components. Careful design of the operating system created a low-power design to enable a long life with small batteries. The hardware included several sensors: temperature, radiant temperature, relative humidity, magnetometer, accelerometer, and light, with an optional occupancy (Passive InfraRed) sensor. The project was the basis of several applications, both internal to the research team and other researchers and professionals at other institutions. Several applications used the sensor hardware as the basis for other complex devices. Other applications used the sensors to improve building performance through interoperating with the building Heating Ventilation and Air-Conditioning (HVAC) system, such as using occupancy and/or distributed temperature sensing to reduce HVAC zone energy while still providing thermal comfort and to reduce peak loads in small commercial buildings. We demonstrated cloud-based energy analytics, implemented a schedule and a Model Predictive Controller in a small commercial building to optimize HVAC energy, occupancy and electricity price. Initial integration of these technological innovations was performed through the creation of execution containers containing the WAVE agent and various driver, proxy, or building system function logic. The research added to the understanding of efficient sensor hardware, secure middleware, time-series data management (high performance database), efficient communication protocols, and interoperating with applications and building systems. The project showed the technical effectiveness and economic feasibility of creating a low-cost, modular, and easy-to-deploy sensor. Through conversations with multiple end users, the research team discovered that many customers wanted data management and services in addition to the sensors. HamiltonIOT developed packages of sensors, border router, and data services to provide a seamless “plug-and-play” sensor deployment. Some customers were willing to pay for higher quality sensors (such as light); some customers wanted a robust enclosure (waterproof).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Personalized whole-brain neural mass models reveal combined Aβ and tau hyperexcitable influences in Alzheimer's disease.

Neuronal dysfunction and cognitive deterioration in Alzheimer's disease (AD) are likely caused by multiple pathophysiological factors. However, mechanistic evidence in humans remains scarce, requiring improved non-invasive techniques and integrative models. We introduce personalized AD computational models built on whole-brain Wilson-Cowan oscillators and incorporating resting-state functional MRI, amyloid-β (Aβ) and tau-PET from 132 individuals in the AD spectrum to evaluate the direct impact of toxic protein deposition on neuronal activity. This subject-specific approach uncovers key patho-mechanistic interactions, including synergistic Aβ and tau effects on cognitive impairment and neuronal excitability increases with disease progression. The data-derived neuronal excitability values strongly predict clinically relevant AD plasma biomarker concentrations (p-tau217, p-tau231, p-tau181, GFAP) and grey matter atrophy obtained through voxel-based morphometry. Furthermore, reconstructed EEG proxy quantities show the hallmark AD electrophysiological alterations (theta band activity enhancement and alpha reductions) which occur with Aβ-positivity and after limbic tau involvement. Microglial activation influences on neuronal activity are less definitive, potentially due to neuroimaging limitations in mapping neuroprotective vs detrimental activation phenotypes. Mechanistic brain activity models can further clarify intricate neurodegenerative processes and accelerate preventive/treatment interventions.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Using Knowledge-Guided Machine Learning To Assess Patterns of Areal Change in Waterbodies across the Contiguous United States

Lake and reservoir surface areas are an important proxy for freshwater availability. Advancements in machine learning (ML) techniques and increased accessibility of remote sensing data products have enabled the analysis of waterbody surface area dynamics on broad spatial scales. However, interpreting the ML results remains a challenge. While ML provides important tools for identifying patterns, the resultant models do not include mechanisms. Thus, the “black-box” nature of ML techniques often lacks ecological meaning. Using ML, we characterized temporal patterns in lake and reservoir surface area change from 1984 to 2016 for 103,930 waterbodies in the contiguous United States. We then employed knowledge-guided machine learning (KGML) to classify all waterbodies into seven ecologically interpretable groups representing distinct patterns of surface area change over time. Many waterbodies were classified as having “no change” (43%), whereas the remaining 57% of waterbodies fell into other groups representing both linear and nonlinear patterns. This analysis demonstrates the potential of KGML not only for identifying ecologically relevant patterns of change across time but also for unraveling complex processes that underpin those changes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using Knowledge-Guided Machine Learning To Assess Patterns of Areal Change in Waterbodies across the Contiguous United States

Lake and reservoir surface areas are an important proxy for freshwater availability. Advancements in machine learning (ML) techniques and increased accessibility of remote sensing data products have enabled the analysis of waterbody surface area dynamics on broad spatial scales. However, interpreting the ML results remains a challenge. While ML provides important tools for identifying patterns, the resultant models do not include mechanisms. Thus, the “black-box” nature of ML techniques often lacks ecological meaning. Using ML, we characterized temporal patterns in lake and reservoir surface area change from 1984 to 2016 for 103,930 waterbodies in the contiguous United States. We then employed knowledge-guided machine learning (KGML) to classify all waterbodies into seven ecologically interpretable groups representing distinct patterns of surface area change over time. Many waterbodies were classified as having “no change” (43%), whereas the remaining 57% of waterbodies fell into other groups representing both linear and nonlinear patterns. This analysis demonstrates the potential of KGML not only for identifying ecologically relevant patterns of change across time but also for unraveling complex processes that underpin those changes.

58 GEOSCIENCES↗
DOE OSTI2022

SDSS-IV MaNGA: Cannibalism Caught in the Act—On the Frequency of Occurrence of Multiple Cores in Brightest Cluster Galaxies

Although it is generally accepted that massive galaxies form in a two-phased fashion, beginning with a rapid mass buildup through intense starburst activities followed by primarily dry mergers that mainly deposit stellar mass at outskirts, the late time stellar mass growth of brightest cluster galaxies (BCGs), the most massive galaxies in the universe, is still not well understood. Several independent measurements have indicated a slower mass growth rate than predictions from theoretical models. We attempt to resolve the discrepancy by measuring the frequency of BCGs with multiple cores, which serve as a proxy of the merger rates in the central region and facilitate a more direct comparison with theoretical predictions. Using 79 BCGs at z = 0.06–0.15 with integral field spectroscopic data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) project, we obtain a multiple-core fraction of 0.11 ± 0.04 at z ≈ 0.1 within an 18 kpc radius from the center, which is comparable to the value of 0.08 ± 0.04 derived from mock observations of 218 simulated BCGs from the cosmological hydrodynamical simulation IllustrisTNG. We find that most cores that appear close to the BCGs from imaging data turn out to be physically associated systems. Anchoring on the similarity in the multiple-core frequency between the MaNGA and IllustrisTNG, we discuss the mass growth rate of BCGs over the past 4.5 Gyr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) (Technical Reference Manual V.4)

The HyRAM+ software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen, methane, and propane infrastructure and transportation systems. HyRAM+ is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of safety, hazards, and risk. HyRAM+ includes generic probabilities for equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of release and flame physics. HyRAM+ integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hazards (thermal effects from jet fires, overpressure effects from delayed ignition), and assessing impact on people and structures. HyRAM+ is developed at Sandia National Laboratories to support the development and revision of national and international codes and standards. HyRAM+ is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM+ version 4.0. The most significant change for HyRAM+ version 4.0 from HyRAM version 3.1 is the incorporation of other alternative fuels, namely methane (as a proxy for natural gas) and propane into the toolkit. This change necessitated significant changes to the installable graphical user interface as well as changes to the back-end Python models. A second major change is the inclusion of physics models for the overpressure associated with the delayed ignition of an unconfined jet/plume of flammable gas.

02 PETROLEUM↗
DOE OSTI2022

Can the Distance–Redshift Relation be Determined from Correlations between Luminosities?

We explore whether an independent determination of the distance–redshift relation, and hence cosmological model parameters, can be obtained from the apparent correlations between two different wave-band luminosities or fluxes, as has been claimed in recent works using the X-ray and ultraviolet luminosities and fluxes of quasars. We show that such an independent determination is possible only if the correlation between luminosities is obtained independently of the cosmological model and measured fluxes and redshifts, for example, based on sound theoretical models or unrelated observations. In particular, we show that if the correlation is determined empirically for two luminosities obtained from fluxes and redshifts, then the method suffers from circularity. In the case where the observed correlation between fluxes in very narrow redshift bins is used as a proxy for the luminosity correlation, we show that one is dealing with a pure tautology with no information on distances and cosmological model. We argue that the problem arises because of the incomplete treatment of the correlation, and we use numerical methods with a joint X-ray and ultraviolet quasar data set to demonstrate this shortcoming.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evaluation and improvement of the parameterization of aerosol hygroscopicity in global climate models using in-situ surface measurements (Final Report)

Aerosols are tiny particles suspended on the atmosphere that can interact with incoming solar radiation and affect the Earth radiative budget. They do so by scattering and absorbing solar radiation, and these properties vary depending on the aerosol size and chemical composition. Moreover, by taking up water from the surrounding air, hygroscopic aerosol particles will grow in size and change their chemical composition, thus modifying their optical properties (scattering and absorption) and their final impact on radiative forcing calculations. An accurate knowledge of aerosol hygroscopicity is crucial for estimating the net radiative impact of aerosols. We took a three-pronged approach to improve our understanding of aerosol hygroscopicity and how it is implemented in Earth system models. In the first part of our project, we developed a benchmark dataset from existing aerosol hygroscopic growth measurements made by tandem nephelometer humidogram systems. We analyzed, using a standardized methodology, observations from 26 in-situ stations around the globe. Measurement data was collected from multiple data providers, reviewed and harmonized to create a consistent dataset of the scattering enhancement factor due to aerosol water uptake. This dataset is archived in several publicly available databases for use by interested researchers. In the second part of the project, we used the benchmark hygroscopicity dataset to perform a global study on aerosol hygroscopicity and aerosol optical properties. Measurements show a global picture of scattering enhancement with larger values for Arctic and marine sites and lower for urban and desert sites. We assessed the RH dependence of aerosol radiative forcing and showed that the overall effect of aerosol hygroscopicity on DARF is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 compared to dry conditions (RH<40%). Finally, we explored using aerosol single scattering albedo (SSA) and scattering Angstrom exponent (SAE) as possible proxies for aerosol hygroscopicity. SSA showed more promise as a surrogate for the scattering enhancement factor than SAE, but neither was ideal. In the third part of the project we evaluated the output of ten Earth system models (ESMs) against the benchmark hygroscopicity dataset. ESMs utilize various schemes for aerosol hygroscopicity which had not been previously tested against observations on a global scale. We found that ESMs currently overestimate scattering enhancement due to hygroscopic growth. Model parameterizations of hygroscopicity and model chemistry are two main factors driving the observed diversity in hygroscopicity simulations among the models. In addition, our study makes several suggestions for modelers, including improving the parameterizations of organic and sea salt aerosol hygroscopicity. Future hygroscopicity model evaluation experiments should include the model data related to particles size which was not available for our study.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

An adjoint-based optimization method for jointly inverting heterogeneous material properties and fault slip from earthquake surface deformation data

SUMMARY Analysis of tectonic and earthquake-cycle associated deformation of the crust can provide valuable insights into the underlying deformation processes including fault slip. How those processes are expressed at the surface depends on the lateral and depth variations of rock properties. The effect of such variations is often tested by forward models based on a priori geological or geophysical information. Here, we first develop a novel technique based on an open-source finite-element computational framework to invert geodetic constraints directly for heterogeneous media properties. We focus on the elastic, coseismic problem and seek to constrain variations in shear modulus and Poisson’s ratio, proxies for the effects of lithology and/or temperature and porous flow, respectively. The corresponding nonlinear inversion is implemented using adjoint-based optimization that efficiently reduces the cost function that includes the misfit between the calculated and observed displacements and a penalty term. We then extend our theoretical and numerical framework to simultaneously infer both heterogeneous Earth’s structure and fault slip from surface deformation. Based on a range of 2-D synthetic cases, we find that both model parameters can be satisfactorily estimated for the megathrust setting-inspired test problems considered. Within limits, this is the case even in the presence of noise and if the fault geometry is not perfectly known. Our method lays the foundation for a future reassessment of the information contained in increasingly data-rich settings, for example, geodetic GNSS constraints for large earthquakes such as the 2011 Tohoku-oki M9 event, or distributed deformation along plate boundaries as constrained from InSAR.

Geochemistry & Geophysics↗
DOE OSTI2024

Permeability and Induced Polarization of Mudstones

Electrical measurements can be used to estimate hydraulic properties such as permeability ( k ) in sedimentary rocks. Previous work has focused on sandstones, siltstones, and carbonates, while investigations on mudstones have rarely been reported. In this study, we report on electrical geophysical measurements for 23 mudstone samples using an experimental approach designed to reliably saturate these low permeability mudstones. The modified Hagen-Poiseuille model linking permeability to the formation factor ( F ) and an effective pore radius ( r ) provides an excellent fit to the data set with a near-constant pore radius, indicating that the effective porosity (1/ F ) is the controlling factor on k . In these samples, the surface area normalized to pore volume ( S por ), frequently used in permeability estimation models, varies by 1–2 orders of magnitude and is thus not a reliable proxy of the inverse effective hydraulic radius. The formation factor also exerts the primary control on induced polarization (IP) parameters, whereas Spor shows no relation to the IP parameters. A strong linear relationship is found between IP parameters (imaginary conductivity and normalized chargeability) and surface conductivity, although the proportionality factor is significantly lower than those observed in more permeable rocks and sediments. Apparent relationships between the polarization strength-derived and time constant-derived geophysical length scales and the effective hydraulic radius appear to be driven by variations in the electrochemical parameters (i.e., specific polarizability and diffusion coefficient). Overall, these findings emphasize that predicting hydraulic properties from electrical measurements in fine-grained rocks remains challenging and requires further investigation into the electrochemical properties involved.

58 GEOSCIENCES↗
DOE OSTI2020

Ultra-short-period Planets Are Stable against Tidal Inspiral

It has been unambiguously shown both in individual systems and at the population level that hot Jupiters experience tidal inspiral before the end of their host stars’ main-sequence lifetimes. Ultra-short-period (USP) planets have orbital periods P < 1 day, rocky compositions, and are expected to experience tidal decay on similar timescales to hot Jupiters if the efficiency of tidal dissipation inside their host stars parameterized as Q{sub ∗}{sup ′} is independent of P and/or secondary mass M {sub p}. Any difference between the two classes of systems would reveal that a model with constant Q{sub ∗}{sup ′} is insufficient. If USP planets experience tidal inspiral, then USP planet systems will be relatively young compared to similar stars without USP planets. Because it is a proxy for relative age, we calculate the Galactic velocity dispersions of USP planet candidate host and non-host stars using data from Gaia Data Release 2 supplemented with ground-based radial velocities. We find that main-sequence USP planet candidate host stars have kinematics consistent with similar stars in the Kepler field without observed USP planets. This indicates that USP planet hosts have similar ages to field stars and that USP planets do not experience tidal inspiral during the main-sequence lifetimes of their host stars. The survival of USP planets requires that Q{sub ∗}{sup ′} ≳ 10{sup 7} at P ≈ 0.7 day and M{sub p}≈2.6 M{sub ⊕}. This result demands that Q{sub ∗}{sup ′} depend on the orbital period and/or mass of the secondary in the range 0.5 day≲P≲5 days and 1 M{sub ⊕}≲M{sub p}≲1000 M{sub ⊕}.

79 ASTRONOMY AND ASTROPHYSICS↗